skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: A computational model for ion and electron energization during macroscale magnetic reconnection
A set of equations is developed that extends the macroscale magnetic reconnection simulation model kglobal to include particle ions. The extension from earlier versions of kglobal, which included only particle electrons, requires the inclusion of the inertia of particle ions in the fluid momentum equation. The new equations will facilitate the exploration of the simultaneous non-thermal energization of ions and electrons during magnetic reconnection in macroscale systems. Numerical tests of the propagation of Alfvén waves and the growth of firehose modes in a plasma with anisotropic electron and ion pressure are presented to benchmark the new model.  more » « less
Award ID(s):
2109083
PAR ID:
10591912
Author(s) / Creator(s):
; ;
Publisher / Repository:
American Institute of Physics
Date Published:
Journal Name:
Physics of Plasmas
Volume:
31
Issue:
6
ISSN:
1070-664X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Magnetic reconnection is a fundamental process in space and astrophysical plasmas that converts magnetic energy to particle energy. Recently, a novel kind of reconnection, called electron-only reconnection, has been observed in Earth's magnetosheath plasma. A defining characteristic of electron-only reconnection is that electron jets are observed but ion jets are absent. This is in contrast with traditional ion-coupled reconnection, where both ions and electrons exhibit outflowing velocity jets. Findings from the Magnetospheric Multiscale mission observations and particle-in-cell simulations show clear signatures of electron heating in electron-only reconnection events, while ions are not heated or cooled in these events. This result is unlike ion-coupled reconnection, where both ions and electrons are heated to varying degrees. The ratio of electron to ion dissipation increases with the local magnetic curvature, indicating that the partition of heat into ions and electrons is dependent on the current-sheet thickness. 
    more » « less
  2. Abstract During magnetospheric substorms, plasma from magnetic reconnection in the magnetotail is thought to reach the inner magnetosphere and form a partial ring current. We simulate this process using a fully kinetic 3D particle‐in‐cell (PIC) numerical code along with a global magnetohydrodynamics (MHD) model. The PIC simulation extends from the solar wind outside the bow shock to beyond the reconnection region in the tail, while the MHD code extends much further and is run for nominal solar wind parameters and a southward interplanetary magnetic field. By the end of the PIC calculation, ions and electrons from the tail reconnection reach the inner magnetosphere and form a partial ring current and diamagnetic current. The primary source of particles to the inner magnetosphere is bursty bulk flows (BBFs) that originate from a complex pattern of reconnection in the near‐Earth magnetotail at to . Most ion acceleration occurs in this region, gaining from 10 to 50 keV as they traverse the sites of active reconnection. Electrons jet away from the reconnection region much faster than the ions, setting up an ambipolar electric field allowing the ions to catch up after approximately 10 ion inertial lengths. The initial energy flux in the BBFs is mainly kinetic energy flux from the ions, but as they move earthward, the energy flux changes to enthalpy flux at the ring current. The power delivered from the tail reconnection in the simulation to the inner magnetosphere is  W, which is consistent with observations. 
    more » « less
  3. Abstract Magnetic reconnection regions in space and astrophysics are known as active particle acceleration sites. There is ample evidence showing that energetic particles can take a substantial amount of converted energy during magnetic reconnection. However, there has been a lack of studies understanding the backreaction of energetic particles at magnetohydrodynamical scales in magnetic reconnection. To address this, we have developed a new computational method to explore the feedback by nonthermal energetic particles. This approach considers the backreaction from these energetic particles by incorporating their pressure into magnetohydrodynamics (MHD) equations. The pressure of the energetic particles is evaluated from their distribution evolved through Parker’s transport equation, solved using stochastic differential equations (SDEs), so we coin the name MHD-SDE. Applying this method to low-βmagnetic reconnection simulations, we find that reconnection is capable of accelerating a large fraction of energetic particles that contain a substantial amount of energy. When the feedback from these particles is included, their pressure suppresses the compression structures generated by magnetic reconnection, thereby mediating particle energization. Consequently, the feedback from energetic particles results in a steeper power-law energy spectrum. These findings suggest that feedback from nonthermal energetic particles plays a crucial role in magnetic reconnection and particle acceleration. 
    more » « less
  4. ABSTRACT We perform 2D particle-in-cell simulations of magnetic reconnection in electron-ion plasmas subject to strong Compton cooling and calculate the X-ray spectra produced by this process. The simulations are performed for trans-relativistic reconnection with magnetization 1 ≤ σ ≤ 3 (defined as the ratio of magnetic tension to plasma rest-mass energy density), which is expected in the coronae of accretion discs around black holes. We find that magnetic dissipation proceeds with inefficient energy exchange between the heated ions and the Compton-cooled electrons. As a result, most electrons are kept at a low temperature in Compton equilibrium with radiation, and so thermal Comptonization cannot reach photon energies $$\sim 100\,$$ keV observed from accreting black holes. Nevertheless, magnetic reconnection efficiently generates $$\sim 100\,$$ keV photons because of mildly relativistic bulk motions of the plasmoid chain formed in the reconnection layer. Comptonization by the plasmoid motions dominates the radiative output and controls the peak of the radiation spectrum Epk. We find Epk ∼ 40 keV for σ = 1 and Epk ∼ 100 keV for σ = 3. In addition to the X-ray peak around 100 keV, the simulations show a non-thermal MeV tail emitted by a non-thermal electron population generated near X-points of the reconnection layer. The results are consistent with the typical hard state of accreting black holes. In particular, we find that the spectrum of Cygnus X-1 is well explained by electron-ion reconnection with σ ∼ 3. 
    more » « less
  5. Abstract Mercury possesses a miniature yet dynamic magnetosphere driven primarily by magnetic reconnection occurring regularly at the magnetopause and in the magnetotail. Using the newly developed Magnetohydrodynamics with Adaptively Embedded Particle‐in‐Cell (MHD‐AEPIC) model coupled with planetary interior, we have performed a series of global simulations with a range of upstream conditions to study in detail the kinetic signatures, asymmetries, and flux transfer events (FTEs) associated with Mercury's dayside magnetopause reconnection. By treating both ions and electrons kinetically, the embedded PIC model reveals crescent‐shaped phase‐space distributions near reconnection sites, counter‐streaming ion populations in the cusp region, and temperature anisotropies within FTEs. A novel metric and algorithm are developed to automatically identify reconnection X‐lines in our 3D simulations. The spatial distribution of reconnection sites as modeled by the PIC code exhibits notable dawn‐dusk asymmetries, likely due to such kinetic effects as X‐line spreading and Hall effects. Across all simulations, simulated FTEs occur quasi‐periodically every 4–9 s. The properties of simulated FTEs show clear dependencies on the upstream solar wind Alfvénic Mach number (MA) and the interplanetary magnetic field orientation, consistent with MESSENGER observations and previous Hall‐MHD simulations. FTEs formed in our MHD‐AEPIC model tend to carry a large amount of open flux, contributing ∼3%–36% of the total open flux generated at the dayside. Taken together, our MHD‐AEPIC simulations provide new insights into the kinetic processes associated with Mercury's magnetopause reconnection that should prove useful for interpreting spacecraft observations, such as those from MESSENGER and BepiColombo. 
    more » « less