Search

In this work, we investigate trajectories that require thrust to maintain periodic structure in the circular restricted three-body problem (CR3BP). We produce bounds in position and velocity space for the energy-constrained reachable set of initial conditions. Our trajectories are energy-optimal and analyzed via linear analysis. We provide validation for our technique and analyze the cost of deviating in various directions to the reference. For our given reference, we find that it is relatively expensive to decrease perilune distance for orbits in the Earth-Moon system., Comment: Submitted to 2024 AAS/AIAA Astrodynamics Specialist Conference Proceedings

In magnetic reconnection, the ion bulk outflow speed and ion heating have been shown to be set by the available reconnecting magnetic energy, i.e., the energy stored in the reconnecting magnetic field ($B_r$). However, recent simulations, observations, and theoretical works have shown that the released magnetic energy is inhibited by upstream ion plasma beta $\beta_{i}$ -- the relative ion thermal pressure normalized to magnetic pressure based on the reconnecting field -- for antiparallel magnetic field configurations. Using kinetic theory and hybrid particle-in-cell simulations, we investigate the effects of $\beta_{i}$ on guide field reconnection. While previous works have suggested that guide field reconnection is uninfluenced by $\beta_{i}$, we demonstrate that the reconnection process is modified and the outflow is reduced for sufficiently large $\beta_{i} > B_g^2/(B_r^2 + B_g^2)$. We develop a theoretical framework that shows that this reduction is consistent with an enhanced exhaust pressure gradient, which reduces the outflow speed as $v_0 \propto 1/\sqrt{\beta_{i}}$. These results apply to systems in which guide field reconnection is embedded in hot plasmas, such as reconnection at the boundary of eddies in fully developed turbulence like the solar wind or the magnetosheath as well as downstream of shocks such as the heliosheath or the mergers of galaxy clusters., Comment: 11 pages, 5 figures, submitted to ApJ

Collisionless shocks are frequently analyzed using the magnetohydrodynamics (MHD) formalism, even though MHD assumes a small mean free path. Yet, isotropy of pressure, fruit of binary collisions and assumed in MHD, may not apply in collisionless shocks. This is especially true within a magnetized plasma, where the field can stabilize an anisotropy. In a previous article \citep{BretJPP2022b}, a model was presented capable of dealing with the anisotropies that may arise at the front crossing. It was solved for any orientation of the field with respect to the shock front. Yet, for some values of the upstream parameters, several downstream solutions were found. Here, we complete the work started in \cite{BretJPP2022b} by showing how to pick the physical solution out of the ones offered by the algebra. This is achieved by 2 means: 1) selecting the solution that has the downstream field obliquity closest to the upstream one. This criterion is exemplified on the parallel case and backed up by Particle-in-Cell simulations. 2) Filtering out solutions which do not satisfy a criteria already invoked to trim multiple solutions in MHD: the evolutionarity criterion, that we assume valid in the collisionless case. The end result is a model in which a given upstream configuration results in a unique, or none (like in MHD), downstream configuration. The largest departure from MHD is found for the case of a parallel shock., Comment: 22 pages, 9 figures, to appear in Journal of Plasma Physics

With the successful impact of the NASA DART spacecraft in the Didymos-Dimorphos binary asteroid system, we provide an initial analysis of the post-impact perturbed binary asteroid dynamics. To compare our simulation results with observations, we introduce a set of "observable elements" calculated using only the physical separation of the binary asteroid, rather than traditional Keplerian elements. Using numerical methods that treat the fully spin-orbit-coupled dynamics, we estimate the system's mass and the impact-induced changes in orbital velocity, semimajor axis, and eccentricity. We find that the changes to the mutual orbit depend strongly on the separation distance between Didymos and Dimorphos at the time of impact. If Dimorphos enters a tumbling state after the impact, this may be observable through changes in the system's eccentricity and orbit period. We also find that any DART-induced reshaping of Dimorphos would generally reduce the required change in orbital velocity to achieve the measured post-impact orbit period and will be assessed by the ESA Hera mission in 2027., Comment: Accepted for publication in PSJ

Collisionless plasma shocks are a common feature of many space and astrophysical systems and are sources of high-energy particles and non-thermal emission, channeling as much as 20\% of the shock's energy into non-thermal particles. The generation and acceleration of these non-thermal particles have been extensively studied, however, how these particles feed back on the shock hydrodynamics has not been fully treated. This work presents the results of self-consistent hybrid particle-in-cell simulations that show the effect of self-generated non-thermal particle populations on the nature of collisionless, quasi-parallel shocks. They contribute to a significant heat flux density upstream of the shock. Non-thermal particles downstream of the shock leak into the upstream region, taking energy away from the shock. This increases the compression ratio, slows the shock down, and flattens the non-thermal population's spectral index for lower Mach number shocks. We incorporate this into a revised theory for the Rankine-Hugoniot jump conditions that include this effect and it shows excellent agreement with simulations. The results have the potential to explain discrepancies between predictions and observations in a wide range of systems, such as inaccuracies of predictions of arrival times of coronal mass ejections and the conflicting radio and x-ray observations of intracluster shocks. These effects will likely need to be included in fluid modeling to accurately predict shock evolution., Comment: 7 pages, 3 figures, a lot of appendix

We analyze the generation of kinetic instabilities and their effect on the energization of ions in non-relativistic, oblique collisionless shocks using a 3D-3V simulation by $\texttt{dHybridR}$, a hybrid particle-in-cell code. At sufficiently high Mach number, quasi-perpendicular and oblique shocks can experience rippling of the shock surface caused by kinetic instabilities arising from free energy in the ion velocity distribution due to the combination of the incoming ion beam and the population of ions reflected at the shock front. To understand the role of the ripple on particle energization, we devise the new instability isolation method to identify the unstable modes underlying the ripple and interpret the results in terms of the governing kinetic instability. We generate velocity-space signatures using the field-particle correlation technique to look at energy transfer in phase space from the isolated instability driving the shock ripple, providing a viewpoint on the different dynamics of distinct populations of ions in phase space. We generate velocity-space signatures of the energy transfer in phase space of the isolated instability driving the shock ripple using the field-particle correlation technique. Together, the field-particle correlation technique and our new instability isolation method provide a unique viewpoint on the different dynamics of distinct populations of ions in phase space and allow us to completely characterize the energetics of the collisionless shock under investigation., Comment: 32 pages, 14 figures, accepted by the Journal of Plasma Physics

Examining energization of kinetic plasmas in phase space is a growing topic of interest, owing to the wealth of data in phase space compared to traditional bulk energization diagnostics. Via the field-particle correlation (FPC) technique and using multiple means of numerically integrating the plasma kinetic equation, we have studied the energization of ions in phase space within oblique collisionless shocks. The perspective afforded to us with this analysis in phase space allows us to characterize distinct populations of energized ions. In particular, we focus on ions which reflect multiple times off the shock front through shock-drift acceleration, and how to distinguish these different reflected populations in phase space using the FPC technique. We further extend our analysis to simulations of three-dimensional shocks undergoing more complicated dynamics, such as shock ripple, to demonstrate the ability to recover the phase space signatures of this energization process in a more general system. This work thus extends previous applications of the FPC technique to more realistic collisionless shock environments, providing stronger evidence of the technique's utility for simulation, laboratory, and spacecraft analysis., Comment: 9 pages, 5 figures

Heliophysics theory and modeling build understanding from fundamental principles to motivate, interpret, and predict observations. Together with observational analysis, they constitute a comprehensive scientific program in heliophysics. As observations and data analysis become increasingly detailed, it is critical that theory and modeling develop more quantitative predictions and iterate with observations. Advanced theory and modeling can inspire and greatly improve the design of new instruments and increase their chance of success. In addition, in order to build physics-based space weather forecast models, it is important to keep developing and testing new theories, and maintaining constant communications with theory and modeling. Maintaining a sustainable effort in theory and modeling is critically important to heliophysics. We recommend that all funding agencies join forces and consider expanding current and creating new theory and modeling programs--especially, 1. NASA should restore the HTMS program to its original support level to meet the critical needs of heliophysics science; 2. a Strategic Research Model program needs to be created to support model development for next-generation basic research codes; 3. new programs must be created for addressing mission-critical theory and modeling needs; and 4. enhanced programs are urgently required for training the next generation of theorists and modelers., Comment: White paper submitted to Heliophysics 2024 Decadal Survey

Shocks waves are a ubiquitous feature of many astrophysical plasma systems, and an important process for energy dissipation and transfer. The physics of these shock waves are frequently treated/modeled as a collisional, fluid MHD discontinuity, despite the fact that many shocks occur in the collisionless regime. In light of this, using fully kinetic, 3D simulations of non-relativistic, parallel propagating collisionless shocks comprised of electron-positron plasma, we detail the deviation of collisionless shocks form MHD predictions for varying magnetization/Alfv\'enic Mach numbers, with particular focus on systems with Alf\'enic Mach numbers much smaller than sonic Mach numbers. We show that the shock compression ratio decreases for sufficiently large upstream magnetic fields, in agreement with the predictions of Bret & Narayan (2018). Additionally, we examine the role of magnetic field strength on the shock front width. This work reinforces a growing body of work that suggest that modeling many astrophysical systems with only a fluid plasma description omits potentially important physics., Comment: 7 pages, 5 figures + a diagram, submitted to MNRAS

A strong super-Alfv\'{e}nic drift of energetic particles (or cosmic rays, CRs) in a magnetized plasma can amplify the magnetic field significantly through non-resonant streaming instability (NRSI). While the traditional analysis is done for an ion current, here we use kinetic particle-in-cell simulations to study how the NRSI behaves when it is driven by electrons or by a mixture of electrons and positrons. In particular, we characterize growth rate, spectrum, and helicity of the unstable modes, as well the level of magnetic field at saturation. Our results are potentially relevant for several space/astrophysical environments (e.g, electron strahl in the solar wind, at oblique non-relativistic shocks, around pulsar wind nebulae) and also in laboratory experiments., Comment: 16 pages, 9 figures, 1 table; Submitted to ApJ; Comments welcome!

Diffusive shock acceleration is a prominent mechanism for producing energetic particles in space and in astrophysical systems. Such energetic particles have long been predicted to affect the hydrodynamic structure of the shock, in turn leading to CR spectra flatter than the test-particle prediction. However, in this work along with a companion paper, C. C. Haggerty and D. Caprioli, 2020, arXiv:2008.12308 [astro-ph.HE], we use self-consistent hybrid (kinetic ions-fluid electrons) simulations to show for the first time how CR-modified shocks actually produce steeper spectra. The steepening is driven by the enhanced advection of CRs embedded in magnetic turbulence downstream of the shock, in what we call the "postcursor". These results are consistent with multi-wavelength observations of supernovae and supernova remnants and have significant phenomenological implications for space/astrophysical shocks in general., Comment: 10 pages, 5 figures, submitted to ApJ, Paper II of a series of 2 (Paper I: Kinetic Simulations of Cosmic-Ray-Modified Shocks II: Hydrodynamics, arXiv:2008.12308 [astro-ph], arXiv:2008.12308)

Collisionless plasma shocks are efficient sources of non-thermal particle acceleration in space and astrophysical systems. We use hybrid (kinetic ions -- fluid electrons) simulations to examine the non-linear feedback of the self-generated energetic particles (cosmic rays, CRs) on the shock hydrodynamics. When CR acceleration is efficient, we find evidence of both an upstream precursor, where the inflowing plasma is compressed and heated, and a downstream postcursor, where the energy flux in CRs and amplified magnetic fields play a dynamical role. For the first time, we assess how non-linear magnetic fluctuations in the postcursor preferentially travel away from the shock at roughly the local Alfv\'en speed with respect to the downstream plasma. The drift of both magnetic and CR energy with respect to the thermal plasma substantially increases the shock compression ratio with respect to the standard prediction, in particular exceeding 4 for strong shocks. Such modifications also have implications for the spectrum of the particles accelerated via diffusive shock acceleration, a significant result detailed in a companion paper, Caprioli, Haggerty & Blasi 2020, arXiv:2009.00007 [astro-ph.HE]., Comment: Paper 1 of 2, 15 pages, 11 figures, companion paper to arXiv:2009.00007 [astro-ph.HE], arXiv:2009.00007

We discuss the recent developments in the theory of diffusive shock acceleration (DSA) by using both first-principle kinetic plasma simulations and analytical theory based on the solution of the convection/diffusion equation. In particular, we show how simulations reveal that the spectra of accelerated particles are significantly steeper than the $E^{-2}$ predicted by the standard theory of DSA for strong shocks, in agreement with several observational pieces of evidence. We single out which standard assumptions of test-particle and non-linear DSA are violated in the presence of strong (self-generated) magnetic turbulence and put forward a novel theory in better agreement with the particle spectra inferred with multi-wavelength observations of young SN remnants, radio-supernovae, and Galactic cosmic rays., Comment: 8 pages, 2 figures, 36th International Cosmic Ray Conference (ICRC2019), held July 24th-August 1st, 2019 in Madison, WI, U.S.A

We present the first plasma simulations obtained with the code dHybridR, a hybrid particle-in-cell code with fluid electrons and both thermal and energetic ions that retain relativistic dynamics. dHybridR is constructed to study astrophysical and space-physics problems where a few energetic non-thermal particles (i.e., cosmic rays, CRs) affect the overall dynamics of a non-relativistic plasma, such as CR-driven instabilities, collisionless shocks, magnetic reconnection, turbulence, etc.In this method paper we provide some applications to linear (resonant/non-resonant CR streaming instability) and strongly non-linear (parallel shocks) problems that show the capabilities of the code. In particular, we provide the first self-consistent hybrid runs that show the acceleration of relativistic ions at non-relativistic shocks; CRs develop a power-law in momentum, which translates to a broken power law in energy that exhibits a steepening around the ion rest mass, as predicted by the theory of diffusive shock acceleration. We present examples of 2D dHybridR runs relevant for fast shocks in radio supernovae, whose evolution can be followed in real time, and 3D runs of low-Mach-number heliospheric shocks, which can be compared with in-situ spacecraft observations., Comment: 14 pages, 12 figures, 1 table

The outflow velocity of jets produced by collisionless magnetic reconnection is shown to be reduced by the ion exhaust temperature in simulations and observations. We derive a scaling relationship for the outflow velocity based on the upstream Alfv\'en speed and the parallel ion exhaust temperature, which is verified in kinetic simulations and observations. The outflow speed reduction is shown to be due to the firehose instability criterion, and so for large enough guide fields this effect is suppressed and the outflow speed reaches the upstream Alfv\'en speed based on the reconnecting component of the magnetic field.

Magnetic reconnection is a ubiquitous phenomenon in turbulent plasmas. It is an important part of the turbulent dynamics and heating of space and astrophysical plasmas. We examine the statistics of magnetic reconnection using a quantitative local analysis of the magnetic vector potential, previously used in magnetohydrodynamics simulations, and now generalized to fully kinetic PIC simulations. Different ways of reducing the particle noise for analysis purposes including multiple smoothing techniques are explored. We find that a Fourier filter applied at the Debye scale is an optimal choice for analyzing PIC data. Finlay, we find a broader distribution of normalized reconnection rates compared to the MHD limit with rates as large as 0.5 but with an average of approximately 0.1.

Heliophysics theory and modeling build understanding from fundamental principles to motivate, interpret, and predict observations. Together with observational analysis, they constitute a comprehensive scientific program in heliophysics. As observations and data analysis become increasingly detailed, it is critical that theory and modeling develop more quantitative predictions and iterate with observations. Advanced theory and modeling can inspire and greatly improve the design of new instruments and increase their chance of success. In addition, in order to build physics based space weather forecast models, it is important to keep developing and testing new theories, and maintaining constant communications with theory and modeling. Maintaining a sustainable effort in theory and modeling is critically important to heliophysics. We recommend that all funding agencies join forces and consider expanding current and creating new theory and modeling programs–especially, 1. NASA should restore the HTMS program to its original support level to meet the critical needs of heliophysics science; 2. a Strategic Research Model program needs to be created to support model development for next-generation basic research codes; 3. new programs must be created for addressing mission-critical theory and modeling needs; and 4. enhanced programs are urgently required for training the next generation of theorists and modelers.