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.
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 \citepBretJPP2022b, 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 \citeBretJPP2022b 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.
The non-resonant (Bell) streaming instability driven by energetic particles is crucial for producing amplified magnetic fields that are key to the acceleration of cosmic rays (CRs) in supernova remnants, around Galactic and extra-galactic CR sources, and for the CR transport. We present a covariant theory for the saturation of the Bell instability, substantiated by self-consistent kinetic simulations, that can be applied to arbitrary CR distributions and discuss its implications in several heliospheric and astrophysical contexts.
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.
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.
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.
The nonresonant streaming instability (Bell instability) plays a pivotal role in the acceleration and confinement of cosmic rays (CRs); yet, the exact mechanism responsible for its saturation and the magnitude of the final amplified magnetic field have not been assessed from first-principles. Using a survey of hybrid simulations (with kinetic ions and fluid electrons), we study the evolution of the Bell instability as a function of the parameters of the CR population. We find that, at saturation, the magnetic pressure in the amplified field is comparable with the initial CR anisotropic pressure, rather than with the CR energy flux as previously argued. These results provide a predictive prescription for the total magnetic field amplification expected in the many astrophysical environments where the Bell instability is important.
Fan Guo, Spiro Antiochos, Paul Cassak, Bin Chen, Xiaohang Chen, Chuanfei Dong, Cooper Downs, Joe Giacalone, Colby C. Haggerty, Hantao Ji, Judith Karpen, James Klimchuk, Wen Li, Xiaocan Li, Mitsuo Oka, Katharine K. Reeves, Marc Swisdak, Weichao Tu 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.
The recent discoveries in the theory of diffusive shock acceleration (DSA) that stem from first-principle kinetic plasma simulations are discussed. When ion acceleration is efficient, the back-reaction of non-thermal particles and self-generated magnetic fields becomes prominent and leads to both enhanced shock compression and particle spectra significantly softer than those predicted by the standard test-particle DSA theory. These results are discussed in the context of the non-thermal phenomenology of astrophysical shocks, with a special focus on the remnant of SN1006.
We study the nature of pressure-strain interaction at reconnection sites, detected by NASA's Magnetospheric Multiscale (MMS) Mission. We employ data from a series of published case studies, including a large-scale reconnection event at the magnetopause, three small-scale reconnection events at the magnetosheath current sheets, and one example of the recently discovered electron-only reconnection. In all instances, we find that the pressure-strain shows signature of conversion into (or from) internal energy at the reconnection site. The electron heating rate is larger than the ion heating rate and the compressive heating is dominant over the incompressive heating rate in all cases considered. The magnitude of thermal energy conversion rate is close to the electromagnetic energy conversion rate in the reconnection region. Although in most cases the pressure-strain interaction indicates that the particle internal energy is increasing, in one case the internal energy is decreasing. These observations indicate that the pressure-strain interaction can be used as an independent measure of energy conversion and dynamics in reconnection regions, in particular independent of measures based on the electromagnetic work. Finally, we explore a selected reconnection site in a turbulent Particle-in-Cell (PIC) simulation which further supports the observational results.
Using fully-kinetic plasma simulations, we study the non-resonant (Bell) streaming instability driven by energetic leptons. We identify the necessary conditions to drive it and the differences from the standard proton-driven case in both linear and saturated stages. A simple analytic theory is presented to explain simulations. Our findings are crucial for understanding the phenomenology of astrophysical environments where only electrons may be accelerated (e.g., oblique shocks) or where relativistic pairs are produced (e.g., around pulsar wind nebulae).
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énic Mach numbers, with particular focus on systems with Alfénic 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.
We derive basic scaling equations for relativistic magnetic reconnection in the general case of asymmetric inflow conditions and obtain predictions for the outflow Lorentz factor and the reconnection rate. Kinetic Particle-in-Cell simulations show that the outflow speeds as well as the nonthermal spectral index are constrained by the inflowing plasma with the weaker magnetic energy per particle, in agreement with the scaling predictions. These results are significant for understanding non-thermal emission from reconnection in magnetically-dominated, astrophysical systems, many of which may be asymmetric in nature. The results provide a quantitative approach for including asymmetry on reconnection in the relativistic regime.
A strong super-Alfvé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.
Cosmic rays (CRs) leave their sources mainly along the local magnetic field; in doing so they excite both resonant and nonresonant modes through streaming instabilities. The excitation of these modes leads to enhanced scattering and in turn to a large pressure gradient that causes the formation of bubbles of gas, CRs, and self-generated magnetic fields expanding into the interstellar medium. By means of hybrid Particle-In-Cell simulations, we show that, by exciting the nonresonant instability, CRs excavate a cavity around their source where the diffusivity is strongly suppressed. This finding invalidates the so far largely adopted flux-tube assumption, under which particles move along magnetic lines even in the nonlinear regime. This phenomenon is general and is expected to occur around any sufficiently powerful CR source in the Galaxy. Our results might provide a physical explanation of the numerous claims of suppressed CR diffusion around Galactic sources such as supernova remnants, pulsar wind nebulae, and stellar clusters.
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.
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én 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].
The nonresonant cosmic ray instability, predicted by Bell (2004), is thought to play an important role in the acceleration and confinement of cosmic rays (CRs) close to supernova remnants. Despite its importance, the exact mechanism responsible for the saturation of the instability has not been determined, and there is no first-principle prediction for the amplitude of the saturated magnetic field. Using a survey of self-consistent kinetic hybrid simulations (with kinetic ions and fluid electrons), we study the saturation of the non-resonant streaming instability as a function of the parameters of both the thermal background plasma and the CR population. The strength of the saturated magnetic field has important implications for both CR acceleration in supernova remnants and CR diffusion in the Galaxy.
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.
Using hybrid simulations (kinetic ions--fluid electrons), we test the linear theory predictions of the cosmic ray (CR) streaming instability. We consider two types of CR distribution functions: a "hot" distribution where CRs are represented by a drifting power law in momentum and an anisotropic "beam" of monochromatic particles. Additionally, for each CR distribution we scan over different CR densities to transition from triggering the resonant to the non-resonant (Bell) streaming instability. We determine the growth rates of these instabilities in simulations by fitting an exponential curve during the linear stage, and we show that they agree well with the theoretical predictions as a function of wave number agree. We also examine the magnetic helicity as a function of time and wave number, finding a general good agreement with the predictions, as well as some unexpected non-linear features to the instability development.
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.
P. Sharma Pyakurel, M. A. Shay, T. D. Phan, W. H. Matthaeus, J. F. Drake, J. M. TenBarge, C. C. Haggerty, K. Klein, P. A. Cassak, T. N. Parashar, M. Swisdak, A. Chasapis Using kinetic particle-in-cell (PIC) simulations, we simulate reconnection conditions appropriate for the magnetosheath and solar wind, i.e., plasma beta (ratio of gas pressure to magnetic pressure) greater than 1 and low magnetic shear (strong guide field). Changing the simulation domain size, we find that the ion response varies greatly. For reconnecting regions with scales comparable to the ion Larmor radius, the ions do not respond to the reconnection dynamics leading to ''electron-only'' reconnection with very large quasi-steady reconnection rates. The transition to more traditional ''ion-coupled'' reconnection is gradual as the reconnection domain size increases, with the ions becoming frozen-in in the exhaust when the magnetic island width in the normal direction reaches many ion inertial lengths. During this transition, the quasi-steady reconnection rate decreases until the ions are fully coupled, ultimately reaching an asymptotic value. The scaling of the ion outflow velocity with exhaust width during this electron-only to ion-coupled transition is found to be consistent with a theoretical model of a newly reconnected field line. In order to have a fully frozen-in ion exhaust with ion flows comparable to the reconnection Alfvén speed, an exhaust width of at least several ion inertial lengths is needed. In turbulent systems with reconnection occurring between magnetic bubbles associated with fluctuations, using geometric arguments we estimate that fully ion-coupled reconnection requires magnetic bubble length scales of at least several tens of ion inertial lengths.
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én 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én 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.
Kinetic plasma turbulence cascade spans multiple scales ranging from macroscopic fluid flow to sub-electron scales. Mechanisms that dissipate large scale energy, terminate the inertial range cascade and convert kinetic energy into heat are hotly debated. Here we revisit these puzzles using fully kinetic simulation. By performing scale-dependent spatial filtering on the Vlasov equation, we extract information at prescribed scales and introduce several energy transfer functions. This approach allows highly inhomogeneous energy cascade to be quantified as it proceeds down to kinetic scales. The pressure work, $-\left( \boldsymbol{P} \cdot \nabla \right) \cdot \boldsymbol{u}$, can trigger a channel of the energy conversion between fluid flow and random motions, which is a collision-free generalization of the viscous dissipation in collisional fluid. Both the energy transfer and the pressure work are strongly correlated with velocity gradients.
Kinetic particle-in-cell simulations are used to identify signatures of the electron diffusion region (EDR) and its surroundings during asymmetric magnetic reconnection. A "shoulder" in the sunward pointing normal electric field (EN > 0) at the reconnection magnetic field reversal is a good indicator of the EDR, and is caused by magnetosheath electron meandering orbits in the vicinity of the x-line. Earthward of the X-line, electrons accelerated by EN form strong currents and crescent-shaped distribution functions in the plane perpendicular to B. Just downstream of the X-line, parallel electric fields create field-aligned crescent electron distribution functions. In the immediate upstream magnetosheath, magnetic field strength, plasma density, and perpendicular electron temperatures are lower than the asymptotic state. In the magnetosphere inflow region, magnetosheath ions intrude resulting in an Earthward pointing electric field and parallel heating of magnetospheric particles. Many of the above properties persist with a guide field of at least unity.
The physical processes that control the partition of released magnetic energy between electrons and ions during reconnection is explored through particle-in-cell simulations and analytical techniques. We demonstrate that the development of a large-scale parallel electric field and its associated potential controls the relative heating of electrons and ions. The potential develops to restrain heated exhaust electrons and enhances their heating by confining electrons in the region where magnetic energy is released. Simultaneously the potential slows ions entering the exhaust below the Alfvénic speed expected from the traditional counterstreaming picture of ion heating. Unexpectedly, the magnitude of the potential and therefore the relative partition of energy between electrons and ions is not a constant but rather depends on the upstream parameters and specifically the upstream electron normalized temperature (electron beta). These findings suggest that the fraction of magnetic energy converted into the total thermal energy may be independent of upstream parameters.
Electron bulk heating during magnetic reconnection with symmetric inflow conditions is examined using kinetic particle-in-cell (PIC) simulations. The degree of electron heating is well correlated with the inflowing Alfvén speed $c_{Ar}$ based on the reconnecting magnetic field through the relation $\Delta T_e = 0.033 \,m_i\,c_{Ar}^2$, where $\Delta T_{e}$ is the increase in electron temperature. For the range of simulations performed, the heating shows almost no correlation with inflow total temperature $T_{tot} = T_i + T_e$ or plasma $\beta$. An out-of-plane (guide) magnetic field of similar magnitude to the reconnecting field does not affect the total heating, but it does quench perpendicular heating, with almost all heating being in the parallel direction. These results are qualitatively consistent with a recent statistical survey of electron heating in the dayside magnetopause, which also found that $\Delta T_e$ was proportional to the inflowing Alfvén speed. The net electron heating varies very little with distance downstream of the x-line. The simulations show at most a very weak dependence of electron heating on the ion to electron mass ratio. In the antiparallel reconnection case, the largely parallel heating is eventually isotropized downstream due a scattering mechanism such as stochastic particle motion or instabilities. The study highlights key properties that must be satisfied by an electron heating mechanism: (1) Preferential heating in the parallel direction; (2) Heating proportional to $m_i\,c_{Ar}^2$; (3) At most a weak dependence on electron mass; and (4) An exhaust electron temperature that varies little with distance from the x-line.