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.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Nonequilibrium ionization and ambipolar diffusion in solar magnetic flux emergence processes
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
Award ID(s):
1714955
PAR ID:
10180269
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Astronomy & Astrophysics
Volume:
633
ISSN:
0004-6361
Page Range / eLocation ID:
A66
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; (, Astronomy & Astrophysics)
    Context. Ambipolar diffusion is a physical mechanism related to the drift between charged and neutral particles in a partially ionized plasma that is key to many different astrophysical systems. However, understanding its effects is challenging due to basic uncertainties concerning relevant microphysical aspects and the strong constraints it imposes on the numerical modeling. Aims. Our aim is to introduce a numerical tool that allows us to address complex problems involving ambipolar diffusion in which, additionally, departures from ionization equilibrium are important or high resolution is needed. The primary application of this tool is for solar atmosphere calculations, but the methods and results presented here may also have a potential impact on other astrophysical systems. Methods. We have developed a new module for the stellar atmosphere Bifrost code that improves its computational capabilities of the ambipolar diffusion term in the generalized Ohm’s law. This module includes, among other things, collision terms adequate to processes in the coolest regions in the solar chromosphere. As the main feature of the module, we have implemented the super time stepping (STS) technique, which allows an important acceleration of the calculations. We have also introduced hyperdiffusion terms to guarantee the stability of the code. Results. We show that to have an accurate value for the ambipolar diffusion coefficient in the solar atmosphere it is necessary to include as atomic elements in the equation of state not only hydrogen and helium, but also the main electron donors like sodium, silicon, and potassium. In addition, we establish a range of criteria to set up an automatic selection of the free parameters of the STS method that guarantees the best performance, optimizing the stability and speed for the ambipolar diffusion calculations. We validate the STS implementation by comparison with a self-similar analytical solution. 
    more » « less
  2. ; ; ; ; (, The Astrophysical journal)
    The presence of magnetic field is crucial in the transport of energy through the solar atmosphere. Recent ground- based and space-borne observations of the quiet Sun have revealed that magnetic field accumulates at photospheric heights, via a local dynamo or from small-scale flux emergence events. However, most of this small-scale magnetic field may not expand into the chromosphere due to the entropy drop with height at the photosphere. Here we present a study that uses a high-resolution 3D radiative MHD simulation of the solar atmosphere with non-gray and non-LTE radiative transfer and thermal conduction along the magnetic field to reveal that (1) the net magnetic flux from the simulated quiet photosphere is not sufficient to maintain a chromospheric magnetic field (on average), (2) processes in the lower chromosphere, in the region dominated by magnetoacoustic shocks, are able to convert kinetic energy into magnetic energy, (3) the magnetic energy in the chromosphere increases linearly in time until the rms of the magnetic field strength saturates at roughly 4–30 G (horizontal average) due to conversion from kinetic energy, (4) and that the magnetic features formed in the chromosphere are localized to this region. 
    more » « less
  3. ; ; ; ; ; (, The Astrophysical Journal)
    The ponderomotive force has been suggested to be the main mechanism to produce the so-called first ionization potential (FIP) effect—the enrichment of low-FIP elements observed in the outer solar atmosphere, in the solar wind, and in solar energetic events. It is well known that the ionization of these elements occurs within the chromosphere. Therefore, this phenomenon is intimately tied to the plasma state in the chromosphere and the corona. For this study, we combine IRIS observations, a single-fluid 2.5D radiative magnetohydrodynamics (MHD) model of the solar atmosphere, including ion–neutral interaction effects and nonequilibrium (NEQ) ionization effects, and a novel multifluid multispecies numerical model (based on the Ebysus code). Nonthermal velocities of Si IV measured from IRIS spectra can provide an upper limit for the strength of any high-frequency Alfvén waves. With the single-fluid model, we investigate the possible impact of NEQ ionization within the region where the FIP may occur, as well as the plasma properties in those regions. These models suggest that regions with strongly enhanced network and type II spicules are possible sites of large ponderomotive forces. We use the plasma properties of the single-fluid MHD model and the IRIS observations to initialize our multifluid models to investigate the multifluid effects on the ponderomotive force associated with Alfvén waves. Our multifluid analysis reveals that collisions and NEQ ionization effects dramatically impact the behavior of the ponderomotive force in the chromosphere, and existing theories may need to be revisited. 
    more » « less
  4. ; (, The Astrophysical Journal)
    Abstract Emerging dimming occurs in isolated solar active regions (ARs) during the early stages of magnetic flux emergence. Observed by the Atmospheric Imaging Assembly, it features a rapid decrease in extreme-ultraviolet (EUV) emission in the 171 Å channel images, and a simultaneous increase in the 211 Å images. Here, we analyze the coronal thermodynamic and magnetic properties to probe its physical origin. We calculate the time-dependent differential emission measures for a sample of 18 events between 2010 and 2012. The emission measure (EM) decrease in the temperature range is well correlated with the EM increase in over eight orders of magnitude. This suggests that the coronal plasma is being heated from the quiet-Sun, sub-MK temperature to 1–2 MK, more typical for ARs. Potential field extrapolation indicates significant change in the local magnetic connectivity: the dimming region is now linked to the newly emerged flux via longer loops. We conclude that emerging dimming is likely caused by coronal heating episodes, powered by reconnection between the emerging and the ambient magnetic fields. 
    more » « less
  5. ; ; ; ; ; ; ; ; (, The Planetary Science Journal)
    Abstract MESSENGER observations suggest a magma ocean formed on proto-Mercury, during which evaporation of metals and outgassing of C- and H-bearing volatiles produced an early atmosphere. Atmospheric escape subsequently occurred by plasma heating, photoevaporation, Jeans escape, and photoionization. To quantify atmospheric loss, we combine constraints on the lifetime of surficial melt, melt composition, and atmospheric composition. Consideration of two initial Mercury sizes and four magma ocean compositions determines the atmospheric speciation at a given surface temperature. A coupled interior–atmosphere model determines the cooling rate and therefore the lifetime of surficial melt. Combining the melt lifetime and escape flux calculations provides estimates for the total mass loss from early Mercury. Loss rates by Jeans escape are negligible. Plasma heating and photoionization are limited by homopause diffusion rates of ∼106kg s−1. Loss by photoevaporation depends on the timing of Mercury formation and assumed heating efficiency and ranges from ∼106.6to ∼109.6kg s−1. The material for photoevaporation is sourced from below the homopause and is therefore energy limited rather than diffusion limited. The timescale for efficient interior–atmosphere chemical exchange is less than 10,000 yr. Therefore, escape processes only account for an equivalent loss of less than 2.3 km of crust (0.3% of Mercury’s mass). Accordingly, ≤0.02% of the total mass of H2O and Na is lost. Therefore, cumulative loss cannot significantly modify Mercury’s bulk mantle composition during the magma ocean stage. Mercury’s high core:mantle ratio and volatile-rich surface may instead reflect chemical variations in its building blocks resulting from its solar-proximal accretion environment. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email