Search for: All records

Abstract 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β_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β_ion guide field reconnection. While previous works have suggested that guide field reconnection is uninfluenced byβ_i, we demonstrate that the reconnection process is modified and the outflow is reduced for sufficiently large ${\beta }_i>(B_r^2+B_g^2)/B_r^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_{\mathrm{out}}\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. 

Abstract 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. 

Abstract We survey 20 reconnection outflow events observed by Magnetospheric MultiScale in the low-βand high-Alfvén-speed regime of the Earth’s magnetotail to investigate the scaling of ion bulk heating produced by reconnection. The range of inflow Alfvén speeds (800–4000 km s⁻¹) and inflow ionβ(0.002–1) covered by this study is in a plasma regime that could be applicable to the solar corona and flare environments. We find that the observed ion heating increases with increasing inflow (upstream) Alfvén speed,V_A, based on the reconnecting magnetic field and the upstream plasma density. However, ion heating does not increase linearly as a function of available magnetic energy per particle, ${{\mathit{m}}_{\mathrm{i}}{\mathit{V}}_{\mathrm{A}}}^2$. Instead, the heating increases progressively less as ${{\mathit{m}}_{\mathrm{i}}{\mathit{V}}_{\mathrm{A}}}^2$rises. This is in contrast to a previous study using the same data set, which found that electron heating in this high-Alfvén-speed and low-βregime scales linearly with ${{\mathit{m}}_{\mathrm{i}}{\mathit{V}}_{\mathrm{A}}}^2$, with a scaling factor nearly identical to that found for the low-V_Aand high-βmagnetopause. Consequently, the ion-to-electron heating ratio in reconnection exhausts decreases with increasing upstreamV_A, suggesting that the energy partition between ions and electrons in reconnection exhausts could be a function of the available magnetic energy per particle. Finally, we find that the observed difference in ion and electron heating scaling may be consistent with the predicted effects of a trapping potential in the exhaust, which enhances electron heating, but reduces ion heating. 

ABSTRACT Cosmic rays (CRs) are thought to escape their sources streaming along the local magnetic field lines. We show that this phenomenon generally leads to the excitation of both resonant and non-resonant streaming instabilities. The self-generated magnetic fluctuations induce particle diffusion in extended regions around the source, so that CRs build up a large pressure gradient. By means of two-dimensional (2D) and three-dimensional (3D) hybrid particle-in-cell simulations, we show that such a pressure gradient excavates a cavity around the source and leads to the formation of a cosmic ray dominated bubble, inside which diffusivity is strongly suppressed. Based on the trends extracted from self-consistent simulations, we estimate that, in the absence of severe damping of the self-generated magnetic fields, the bubble should keep expanding until pressure balance with the surrounding medium is reached, corresponding to a radius of ∼10–50 pc. The implications of the formation of these regions of low diffusivity for sources of Galactic CRs are discussed. Special care is devoted to estimating the self-generated diffusion coefficient and the grammage that CRs might accumulate in the bubbles before moving into the interstellar medium. Based on the results of 3D simulations, general considerations on the morphology of the γ-ray and synchrotron emission from these extended regions also are outlined. 

Collisionless shocks are frequently analysed using the magnetohydrodynamics (MHD) formalism, even though MHD assumes a small mean free path. Yet, isotropy of pressure, the 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 (Bret & Narayan,J. Plasma Phys., vol. 88, no. 6, 2022b, p. 905880615), 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 Bret & Narayan (J. Plasma Phys., vol. 88, no. 6, 2022b, p. 905880615) by showing how to pick the physical solution out of the ones offered by the algebra. This is achieved by 2 means: (i) 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. (ii) 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 no downstream configuration (as in MHD). The largest departure from MHD is found for the case of a parallel shock. 

We analyse the generation of kinetic instabilities and their effect on the energization of ions in non-relativistic, oblique collisionless shocks using a 3D-3V (three spatial with three velocity components) simulation by 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 a 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. 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. 

Abstract A strong super-Alfvénic drift of energetic particles (or cosmic rays) in a magnetized plasma can amplify the magnetic field significantly through nonresonant 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 the growth rate, spectrum, and helicity of the unstable modes, as well the level of the magnetic field at saturation. Our results are potentially relevant for several space/astrophysical environments (e.g., electron strahl in the solar wind, at oblique nonrelativistic shocks, around pulsar wind nebulae), and also in laboratory experiments. 

ABSTRACT 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/modelled as a collisional, fluid magnetohydrodynamic (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 theoretical predictions from previous works. 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 modelling many astrophysical systems with only a fluid plasma description omits potentially important physics.