Search for: All records

Abstract We demonstrate a novel setup for hybrid particle-in-cell simulations designed to isolate the physics of the shock precursor over long time periods for significantly lower computational cost than previous methods. This is achieved using a “faux-shock” or shock-like boundary condition on one edge of our simulation domain such that particles that interact with the boundary either pass through it or are reflected off of it with a change in momentum that mimics scattering in the downstream. We show that our faux-shock setup reproduces the same fluid quantities and phase spaces as traditional shock simulations, including those which could otherwise only be done in 3D, with higher particle resolution and for reduced computational cost. While the method involves an assumed boundary condition, it nonetheless captures the essential physics of interest, establishing it as a reliable and efficient tool for future self-consistent studies of instabilities driven by cosmic rays in a shock upstream medium. 

Abstract Collisionless plasma shocks are a common feature of many space and astrophysical systems. They are sources of high-energy particles and nonthermal emission, channeling as much as 20% of the shock’s energy into nonthermal particles. The generation and acceleration of these nonthermal particles have been previously studied and shown to affect shock hydrodynamics to the zeroth order. In this work, we use self-consistent hybrid particle-in-cell simulations to examine the effect of self-generated nonthermal ion populations on the nature of collisionless, quasi-parallel shocks. Accelerated nonthermal particles downstream of the shock diffuse into the upstream region, taking energy away from the shock, which increases the compression ratio, slows the shock down, and flattens the nonthermal population’s spectral index for lower-Mach-number shocks. We show that this enhances shock compressibility when the heat flux is included in the Rankine–Hugoniot jump conditions, results that are roughly consistent with previous theories of “cosmic-ray-modified shocks.” Additionally, the simulation data show that heat flux and enthalpy flux cancels out in the upstream region, yielding a relatively simple, alternative closure for the jump conditions which accurately predict for the shock speed and compression ratio. The results have the potential to explain discrepancies between predictions and observations in a wide range of systems, such as inaccuracies in predictions of the 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 predict shock evolution accurately. 

Abstract We investigate the kinetic effects of upstream, magnetic field-aligned flow shear on antiparallel magnetic reconnection using 2.5D particle-in-cell simulations. Our results demonstrate that flow shear significantly alters the reconnection process, leading to enhanced ion heating, reduced outflow speeds, and a modified reconnection geometry. In contrast to previous Hall magnetohydrodynamic studies, we find that reconnection becomes a more efficient plasma heating mechanism in the presence of sub-Alfvénic flow shear, with ion heating increasing by as much as 300%. This enhanced heating is achieved by efficiently converting the incoming flow shear energy into thermal energy through isotropization in the exhaust. The enhanced heating leads to a pressure gradient away form thex-line exerting a force that reduces the outflow jet speed and slows down the reconnection process. This conversion is due to beam selection effects, mixing, and scattering in the exhaust. A theoretical model is developed that predicts well the exhaust heating and outflow speed reduction. These results offer a potential explanation for recent Parker Solar Probe observations of suppressed reconnection in the presence of flow shear and carry significant implications for energy dissipation in turbulent plasmas. 

Abstract We investigate the process of diffusive shock acceleration of particles with mass number to charge number ratiosA/Q > 1, e.g., partially ionized heavy ions. To this end, we introduce helium- and carbon-like ions at solar abundances into two-dimensional hybrid (kinetic ions–fluid electrons) simulations of nonrelativistic collisionless shocks. This study yields three main results: (1) Heavy ions are preferentially accelerated compared to hydrogen. For typical solar abundances, the energy transferred to accelerated helium ions is comparable to, or even exceeds, that of hydrogen, thereby enhancing the overall shock acceleration efficiency. (2) Accelerated helium ions contribute to magnetic field amplification, which increases the maximum attainable particle energy and steepens the spectra of accelerated particles. (3) The efficient acceleration of helium significantly enhances the production of hadronicγ-rays and neutrinos, likely dominating the one due to hydrogen. These effects should be taken into account, especially when modeling strong space and astrophysical shocks. 

ABSTRACT The light curves of radioactive transients, such as supernovae and kilonovae, are powered by the decay of radioisotopes, which release high-energy leptons through $$\beta ^+$$ and $$\beta ^-$$ decays. These leptons deposit energy into the expanding ejecta. As the ejecta density decreases during expansion, the plasma becomes collisionless, with particle motion governed by electromagnetic forces. In such environments, strong or turbulent magnetic fields are thought to confine particles, though the origin of these fields and the confinement mechanism have remained unclear. Using fully kinetic particle-in-cell (PIC) simulations, we demonstrate that plasma instabilities can naturally confine high-energy leptons. These leptons generate magnetic fields through plasma streaming instabilities, even in the absence of pre-existing fields. The self-generated magnetic fields slow lepton diffusion, enabling confinement, and transferring energy to thermal electrons and ions. Our results naturally explain the positron trapping inferred from late-time observations of thermonuclear and core-collapse supernovae. Furthermore, they suggest potential implications for electron dynamics in the ejecta of kilonovae. We also estimate synchrotron radio luminosities from positrons for Type Ia supernovae and find that such emission could only be detectable with next-generation radio observatories from a Galactic or local-group supernova in an environment without any circumstellar material. 

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