skip to main content


Title: Nonequilibrium Ionization Modeling of Petschek-type Shocks in Reconnecting Current Sheets in Solar Eruptions
Abstract

Nonequilibrium ionization (NEI) is essentially required for astrophysical plasma diagnostics once the plasma status departs from the assumption of ionization equilibrium. In this work, we perform fast NEI calculations combined with magnetohydrodynamic (MHD) simulations and analyze the ionization properties of a Petschek-type magnetic reconnection current sheet during solar eruptions. Our simulation reveals Petschek-type slow-mode shocks in the classical Spitzer thermal conduction models and conduction flux-limitation situations. The results show that under-ionized features can be commonly found in shocked reconnection outflows and thermal halo regions outside the shocks. The departure from equilibrium ionization strongly depends on plasma density. In addition, this departure is sensitive to the observable target temperature: the high-temperature iron ions are strongly affected by the effects of NEI. The under-ionization also affects the synthetic SDO/AIA intensities, which indicates that the reconstructed hot reconnection current sheet structure may be significantly underestimated either for temperature or apparent width. We also perform an MHD-NEI analysis on the reconnection current sheet in the classical solar flare geometry. Finally, we show the potential reversal between the under-ionized and over-ionized states at the lower tip of reconnection current sheets where the downward outflow collides with closed magnetic loops, which can strongly affect multiple SDO/AIA band ratios along the reconnection current sheet.

 
more » « less
Award ID(s):
2108438
NSF-PAR ID:
10394552
Author(s) / Creator(s):
; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
943
Issue:
2
ISSN:
0004-637X
Format(s):
Medium: X Size: Article No. 111
Size(s):
["Article No. 111"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Magnetic reconnection is the key mechanism for energy release in solar eruptions, where the high-temperature emission is the primary diagnostic for investigating the plasma properties during the reconnection process. Non-thermal broadening of high-temperature lines has been observed in both the reconnection current sheet (CS) and flare loop-top regions by UV spectrometers, but its origin remains unclear. In this work, we use a recently developed three-dimensional magnetohydrodynamic (MHD) simulation to model magnetic reconnection in the standard solar flare geometry and reveal highly dynamic plasma flows in the reconnection regions. We calculate the synthetic profiles of the Fe XXI 1354 Å line observed by the Interface Region Imaging Spectrograph (IRIS) spacecraft by using parameters of the MHD model, including plasma density, temperature, and velocity. Our model shows that the turbulent bulk plasma flows in the CS and flare loop-top regions are responsible for the non-thermal broadening of the Fe XXI emission line. The modeled non-thermal velocity ranges from tens of km s −1 to more than two hundred km s −1 , which is consistent with the IRIS observations. Simulated 2D spectral line maps around the reconnection region also reveal highly dynamic downwflow structures where the high non-thermal velocity is large, which is consistent with the observations as well. 
    more » « less
  2. Magnetic reconnection is an important process in astrophysical environments, as it reconfigures magnetic field topology and converts magnetic energy into thermal and kinetic energy. In extreme astrophysical systems, such as black hole coronae and pulsar magnetospheres, radiative cooling modifies the energy partition by radiating away internal energy, which can lead to the radiative collapse of the reconnection layer. In this paper, we perform two- and three-dimensional simulations to model the MARZ (Magnetic Reconnection on Z) experiments, which are designed to access cooling rates in the laboratory necessary to investigate reconnection in a previously unexplored radiatively cooled regime. These simulations are performed in GORGON, an Eulerian two-temperature resistive magnetohydrodynamic code, which models the experimental geometry comprising two exploding wire arrays driven by 20 MA of current on the Z machine (Sandia National Laboratories). Radiative losses are implemented using non-local thermodynamic equilibrium tables computed using the atomic code Spk, and we probe the effects of radiation transport by implementing both a local radiation loss model and$P_{1/3}$multi-group radiation transport. The load produces highly collisional, super-Alfvénic (Alfvén Mach number$M_A \approx 1.5$), supersonic (Sonic Mach number$M_S \approx 4-5$) strongly driven plasma flows which generate an elongated reconnection layer (Aspect Ratio$L/\delta \approx 100$, Lundquist number$S_L \approx 400$). The reconnection layer undergoes radiative collapse when the radiative losses exceed the rates of ohmic and compressional heating (cooling rate/hydrodynamic transit rate =$\tau _{\text {cool}}^{-1}/\tau _{H}^{-1}\approx 100$); this generates a cold strongly compressed current sheet, leading to an accelerated reconnection rate, consistent with theoretical predictions. Finally, the current sheet is also unstable to the plasmoid instability, but the magnetic islands are extinguished by strong radiative cooling before ejection from the layer.

     
    more » « less
  3. Context. Magnetic flux emergence from the solar interior has been shown to be a key mechanism for unleashing a wide variety of phenomena. However, there are still open questions concerning the rise of the magnetized plasma through the atmosphere, mainly in the chromosphere, where the plasma departs from local thermodynamic equilibrium (LTE) and is partially ionized. Aims. We aim to investigate the impact of the nonequilibrium (NEQ) ionization and recombination and molecule formation of hydrogen, as well as ambipolar diffusion, on the dynamics and thermodynamics of the flux emergence process. Methods. Using the radiation-magnetohydrodynamic Bifrost code, we performed 2.5D numerical experiments of magnetic flux emergence from the convection zone up to the corona. The experiments include the NEQ ionization and recombination of atomic hydrogen, the NEQ formation and dissociation of H 2 molecules, and the ambipolar diffusion term of the generalized Ohm’s law. Results. Our experiments show that the LTE assumption substantially underestimates the ionization fraction in most of the emerged region, leading to an artificial increase in the ambipolar diffusion and, therefore, in the heating and temperatures as compared to those found when taking the NEQ effects on the hydrogen ion population into account. We see that LTE also overestimates the number density of H 2 molecules within the emerged region, thus mistakenly magnifying the exothermic contribution of the H 2 molecule formation to the thermal energy during the flux emergence process. We find that the ambipolar diffusion does not significantly affect the amount of total unsigned emerged magnetic flux, but it is important in the shocks that cross the emerged region, heating the plasma on characteristic times ranging from 0.1 to 100 s. We also briefly discuss the importance of including elements heavier than hydrogen in the equation of state so as not to overestimate the role of ambipolar diffusion in the atmosphere. 
    more » « less
  4. Abstract

    We analyze the structure and evolution of ribbons from the M7.3 SOL2014-04-18T13 flare using ultraviolet images from the Interface Region Imaging Spectrograph and the Solar Dynamics Observatory (SDO)/Atmospheric Imaging Assembly (AIA), magnetic data from the SDO/Helioseismic and Magnetic Imager, hard X-ray (HXR) images from the Reuven Ramaty High Energy Solar Spectroscopic Imager, and light curves from the Fermi/Gamma-ray Burst Monitor, in order to infer properties of coronal magnetic reconnection. As the event progresses, two flare ribbons spread away from the magnetic polarity inversion line. The width of the newly brightened front along the extension of the ribbon is highly intermittent in both space and time, presumably reflecting nonuniformities in the structure and/or dynamics of the flare current sheet. Furthermore, the ribbon width grows most rapidly in regions exhibiting concentrated nonthermal HXR emission, with sharp increases slightly preceding the HXR bursts. The light curve of the ultraviolet emission matches the HXR light curve at photon energies above 25 keV. In other regions the ribbon-width evolution and light curves do not temporally correlate with the HXR emission. This indicates that the production of nonthermal electrons is highly nonuniform within the flare current sheet. Our results suggest a strong connection between the production of nonthermal electrons and the locally enhanced perpendicular extent of flare ribbon fronts, which in turn reflects the inhomogeneous structure and/or reconnection dynamics of the current sheet. Despite this variability, the ribbon fronts remain nearly continuous, quasi-one-dimensional features. Thus, although the reconnecting coronal current sheets are highly structured, they remain quasi-two-dimensional and the magnetic energy release occurs systematically, rather than stochastically, through the volume of the reconnecting magnetic flux.

     
    more » « less
  5. Abstract

    The formation, development, and impact of slow shocks in the upstream regions of reconnecting current layers are explored. Slow shocks have been documented in the upstream regions of magnetohydrodynamic (MHD) simulations of magnetic reconnection as well as in similar simulations with thekglobalkinetic macroscale simulation model. They are therefore a candidate mechanism for preheating the plasma that is injected into the current layers that facilitate magnetic energy release in solar flares. Of particular interest is their potential role in producing the hot thermal component of electrons in flares. During multi-island reconnection, the formation and merging of flux ropes in the reconnecting current layer drives plasma flows and pressure disturbances in the upstream region. These pressure disturbances steepen into slow shocks that propagate along the reconnecting component of the magnetic field and satisfy the expected Rankine–Hugoniot jump conditions. Plasma heating arises from both compression across the shock and the parallel electric field that develops to maintain charge neutrality in a kinetic system. Shocks are weaker at lower plasmaβ, where shock steepening is slow. While these upstream slow shocks are intrinsic to the dynamics of multi-island reconnection, their contribution to electron heating remains relatively minor compared with that from Fermi reflection and the parallel electric fields that bound the reconnection outflow.

     
    more » « less