Engineering topics
Zhang, Qile
Publications and source records attributed to Zhang, Qile.
Plasma Dynamics and Nonthermal Particle Acceleration in 3D Nonrelativistic Magnetic Reconnection
Understanding plasma dynamics and nonthermal particle acceleration in 3D magnetic reconnection has been a long-standing challenge. In this paper, we explore these problems by performing large-scale fully kinetic simulations of multi-X-line plasmoid reconnection with various parameters in both the weak- and strong-guide-field regimes. In each regime, we have identified its unique 3D dynamics that lead to field-line chaos and efficient acceleration, and we have achieved nonthermal acceleration of both electrons and protons into power-law spectra. The spectral indices agree well with a simple Fermi acceleration theory that includes guide-field dependence. In the low-guide-field regime, the flux rope kink instability governs the 3D dynamics for efficient acceleration. The weak dependence of the spectra on the ion-to-electron mass ratio and β ($\ll$1) implies that the particles are sufficiently magnetized for Fermi acceleration in our simulations. While both electrons and protons are injected at reconnection exhausts, protons are primarily injected by perpendicular electric fields through Fermi reflections and electrons are injected by a combination of perpendicular and parallel electric fields. The magnetic power spectra agree with in situ magnetotail observations, and the spectral index may reflect a reconnection-driven size distribution of plasmoids instead of the Goldreich–Sridhar vortex cascade. As the guide field becomes stronger, the oblique flux ropes of large sizes capture the main 3D dynamics for efficient acceleration. Intriguingly, the oblique flux ropes can also experience flux rope kink instability, to drive extra 3D dynamics. This work has broad implications for 3D reconnection dynamics and particle acceleration in heliophysics and astrophysics.
Multispecies Ion Acceleration in 3D Magnetic Reconnection with Hybrid-Kinetic Simulations
Magnetic reconnection drives multispecies particle acceleration broadly in space and astrophysics. Here we perform the first 3D hybrid simulations (fluid electrons, kinetic ions) that contain sufficient scale separation to produce nonthermal heavy-ion acceleration, with fragmented flux ropes critical for accelerating all species. We demonstrate the acceleration of all ion species (up to Fe) into power-law spectra with similar indices, by a common Fermi acceleration mechanism. The upstream ion velocities influence the first Fermi reflection for injection. The subsequent onsets of Fermi acceleration are delayed for ions with lower charge-mass ratios (Q/M), until growing flux ropes magnetize them. This leads to a species-dependent maximum energy/nucleon ∝(Q/M) α . These findings are consistent with in situ observations in reconnection regions, suggesting Fermi acceleration as the dominant multispecies ion acceleration mechanism.
Comment on “Nonideal Fields Solve the Injection Problem in Relativistic Reconnection”
Sironi (hereafter S22) reported the correlation between high-energy particles and their crossings of regions with electric field larger than magnetic field (E > B regions) in kinetic simulations of relativistic magnetic reconnection. They claim that E > B regions (for vanishing guide fields) accelerate particles to the injection energy γ inj ~ σ (magnetization). S22 shows that if test particles are reset to low energies in E > B regions, injection is suppressed. Here, we reexamine these claims using a simulation resembling the reference case in S22 with no guide field. We show that E > B regions contribute very little to injection (~ 10% γ inj ) as they only host particles for a short duration. The energization before any E > B crossings has a comparable contribution, indicating E > B regions are not unique for pre-acceleration. A new test-particle simulation that zeroes out E during E > B does not strongly influence injection. We suggest that the procedure to exclude E > B acceleration in S22 partly removes acceleration outside E > B, leading to a false conclusion.
Particle Acceleration during Magnetic Reconnection in Solar Flares and Its Parameter Dependence [Slide]
We performed 3D hybrid simulations (particle ions and fluid electrons) of magnetic reconnection and we have achieved, for the first time, nonthermal acceleration of all ion species (up to Iron) into power-law energy spectra. The spectra have a scaling of maximum energy per nucleon as a function of species consistent with satellite observations by NASA’s Parker Solar Probe mission.
Properties and Acceleration Mechanisms of Electrons Up To 200 keV Associated With a Flux Rope Pair and Reconnection X-Lines Around It in Earth's Plasma Sheet
The properties and acceleration mechanisms of electrons (<200 keV) associated with a pair of tailward traveling flux ropes and accompanied reconnection X-lines in Earth's plasma sheet are investigated with MMS measurements. Energetic electrons are enhanced on both boundaries and core of the flux ropes. The power-law spectra of energetic electrons near the X-lines and in flux ropes are harder than those on flux rope boundaries. Theoretical calculations show that the highest energy of adiabatic electrons is a few keV around the X-lines, tens of keV immediately downstream of the X-lines, hundreds of keV on the flux rope boundaries, and a few MeV in the flux rope cores. The X-lines cause strong energy dissipation, which may generate the energetic electron beams around them. The enhanced electron parallel temperature can be caused by the curvature-driven Fermi acceleration and the parallel electric potential. Betatron acceleration due to the magnetic field compression is strong on flux rope boundaries, which enhances energetic electrons in the perpendicular direction. Electrons can be trapped between the flux rope pair due to mirror force and parallel electric potential. Electrostatic structures in the flux rope cores correspond to potential drops up to half of the electron temperature. The energetic electrons and the electron distribution functions in the flux rope cores are suggested to be transported from other dawn-dusk directions, which is a 3-dimensional effect. The acceleration and deceleration of the Betatron and Fermi processes appear alternately indicating that the magnetic field and plasma are turbulent around the flux ropes.
Particle Acceleration during Magnetic Reconnection in Solar Flares and Its Parameter Dependence [Slide]
Time evolution of the electron and proton energy spectra in our 3D reconnection simulation. For the first time, the spectra form power laws with stable slopes over time for both species.
Variable Ion Compositions of Solar Energetic Particle Events in the Inner Heliosphere: A Field Line Braiding Model with Compound Injections
We propose a model for interpreting highly variable ion composition ratios in solar energetic particle (SEP) events recently observed by the Parker Solar Probe (PSP) at 0.3–0.45 au. We use numerical simulations to calculate SEP propagation in a turbulent interplanetary magnetic field with a Kolmogorov power spectrum from large scales down to the gyration scale of energetic particles. We show that when the source regions of different species are offset by a distance comparable to the size of the source regions, the observed energetic particle composition He/H can be strongly variable over more than two orders of magnitude, even if the source ratio is at the nominal value. Assuming a 3 He/ 4 He source ratio of 10% in impulsive 3 He-rich events and the same spatial offset of the source regions, the 3 He/ 4 He ratio at observation sites also vary considerably. The variability of the ion composition ratios depends on the radial distance, which can be tested by observations made at different radial locations. We discuss the implications of these results on the variability of ion composition of impulsive events and on further PSP and Solar Orbiter observations close to the Sun.
Efficient Nonthermal Ion and Electron Acceleration Enabled by the Flux-Rope Kink Instability in 3D Nonrelativistic Magnetic Reconnection
The relaxation of field-line tension during magnetic reconnection gives rise to a universal Fermi acceleration process involving the curvature drift of particles. However, the efficiency of this mechanism is limited by the trapping of energetic particles within flux ropes. Using 3D fully kinetic simulations, we demonstrate that the flux-rope kink instability leads to strong field-line chaos in weak-guide-field regimes where the Fermi mechanism is most efficient, thus allowing particles to transport out of flux ropes and undergo further acceleration. As a consequence, both ions and electrons develop clear power-law energy spectra that contain a significant fraction of the released energy. The low-energy bounds are determined by the injection physics, while the high-energy cutoffs are limited only by the system size. Furthermore, these results have strong relevance to observations of nonthermal particle acceleration in space and astrophysics.
Magnetic Energy Release, Plasma Dynamics, and Particle Acceleration in Relativistic Turbulent Magnetic Reconnection
In strongly magnetized astrophysical plasma systems, magnetic reconnection is believed to be the primary process during which explosive energy release and particle acceleration occur, leading to significant high-energy emission. Past years have witnessed active development of kinetic modeling of relativistic magnetic reconnection, supporting this magnetically dominated scenario. A much less explored issue in studies of relativistic reconnection is the consequence of three-dimensional dynamics, where turbulent structures are naturally generated as various types of instabilities develop. This paper presents a series of three-dimensional, fully kinetic simulations of relativistic turbulent magnetic reconnection (RTMR) in positron–electron plasmas with system domains much larger than kinetic scales. Our simulations start from a force-free current sheet with several different modes of long-wavelength magnetic field perturbations, which drive additional turbulence in the reconnection region. Because of this, the current layer breaks up and the reconnection region quickly evolves into a turbulent layer filled with coherent structures such as flux ropes and current sheets. Additionally, we find that plasma dynamics in RTMR is vastly different from its 2D counterpart in many aspects. The flux ropes evolve rapidly after their generation, and can be completely disrupted by the secondary kink instability. This turbulent evolution leads to superdiffusive behavior of magnetic field lines as seen in MHD studies of turbulent reconnection. Meanwhile, nonthermal particle acceleration and the timescale for energy release can be very fast and do not depend strongly on the turbulence amplitude. Additionally, the main acceleration mechanism is a Fermi-like acceleration process supported by the motional electric field, whereas the nonideal electric field acceleration plays a subdominant role. We also discuss possible observational implications of three-dimensional RTMR in high-energy astrophysics.
Simultaneous Power-law Formation of Both Electrons and Protons in 3D Non-relativistic Magnetic Reconnection
Explore the source record for details and available documents.
Instabilities and turbulence in low-$β$ guide field reconnection exhausts with kinetic Riemann simulations
The role of turbulence in low-$β$, guide-field reconnection exhausts is explored in 2D reconnection and 2D and 3D Riemann simulations. The structure of the exhaust and associated turbulence is controlled by a pair of rotational discontinuities (RDs) at the exhaust boundary and a pair of slow shocks (SSs) that are generated by counterstreaming ions beams. In 2D, the exhaust develops large-amplitude striations at the ion Larmor radius scale that are produced by electron-beam-driven ion cyclotron waves. The electron beams driving the instability are injected into the exhaust from one of the RDs. However, in 3D Riemann simulations, the additional dimension (in the out-of-plane direction) results in strong Buneman and electron-electron streaming instabilities at the RD which suppress the electron beam formation and therefore the striations in the exhaust. The strength of the streaming instabilities at the RD is controlled by the ratio of the electron thermal speed to Alfvén speed, with the lower thermal speed being more unstable. In the 3D simulations, an ion-ion streaming instability acts to partially thermalize the counterstreaming ion beams at the SSs. This instability is controlled by the ratio of the sound speed to Alfvén speed and is expected to be stable in the low $β$ solar corona. The results suggest that in a guide field reconnection exhaust with 1 $\gg$ $β >$ $m_e/m_i$, the kinetic-scale turbulence that develops will be too weak to play a significant role in energy conversion and particle acceleration. Therefore, the energy conversion will be mostly controlled by laminar physics or multi-x-line reconnection.
Particle heating and energy partition in low- β guide field reconnection with kinetic Riemann simulations
Kinetic Riemann simulations have been completed to investigate particle heating during guide field reconnection in the low-β environment of the inner heliosphere and the solar corona. The reconnection exhaust is bounded by two rotational discontinuities (RDs), and two slow shocks (SSs) form within the exhaust as in magnetohydrodynamic (MHD) models. At the RDs, ions are accelerated by the magnetic field tension to drive the reconnection outflow as well as flows in the out-of-plane direction. The out-of-plane flows stream toward the midplane and meet to drive the SSs. The SSs differ greatly from those in the MHD model. The turbulence at the shock fronts and both upstream and downstream is weak, and so the shocks are laminar and produce little dissipation. Downstream of the SSs, the counterstreaming ion beams lead to higher density, which leads to a positive potential between the SSs which acts to confine the downstream electrons to maintain charge neutrality. The potential accelerates electrons from upstream of the SSs to the downstream region and traps a small fraction but only modestly increases the downstream electron temperature above the upstream value. In the low-β limit, the released magnetic energy is split between bulk flow and ion heating with little energy going to electrons. That the model produces neither strong electron heating nor an energetic electron component implicates that other mechanisms, such as multiple x-line reconnection, are required to explain energetic electron production in large flares. The model can be tested with the expected data from the Parker Solar Probe.