skip to main content


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
NSF-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. Aims.Recurring jets are observed in the solar atmosphere. They can erupt intermittently over a long period of time. By the observation of intermittent jets, we wish to understand what causes the characteristics of the periodic eruptions.

    Methods.We report intermittent jets observed by the Goode Solar Telescope (GST) with the TiO Broadband Filter Imager (BFI), the Visible Imaging Spectrometer (VIS) in Hα, and the Near-InfraRed Imaging Spectropolarimeter (NIRIS). The analysis was aided and complemented by 1400 Å and 2796 Å data from the Interface Region Imaging Spectrograph (IRIS). These observational instruments allowed us to analyze the temporal characteristics of the jet events. By constructing the Hαdopplergrams, we found that the plasma first moves upward, but during the second phase of the jet, the plasma flows back. Working with time slice diagrams, we investigated the characteristics of the jet dynamics.

    Results.The jet continued for up to 4 h. The time-distance diagram shows that the peak of the jet has clear periodic-eruption characteristics (5 min) during 18:00 UT–18:50 UT. We also found a periodic brightening phenomenon (5 min) during the jet bursts in the observed bands in the transition region (1400 Å and 2796 Å), which may be a response to intermittent jets in the upper solar atmosphere. The time lag is 3 min. Evolutionary images in the TiO band revealed a horizontal movement of the granulation at the location of the jet. By comparison to the quiet region of the Sun, we found that the footpoint of the jet is enhanced at the center of the Hαspectral line profile, without significant changes in the line wings. This suggests prolonged heating at the footpoint of the jet. In the mixed-polarity magnetic field region of the jet, we observed the emergence of magnetic flux, its cancellation, and shear, indicating possible intermittent magnetic reconnection. This is confirmed by the nonlinear force-free field model, which was reconstructed using the magneto-friction method.

    Conclusions.The multiwavelength analysis indicates that the events we studied were triggered by magnetic reconnection that was caused by mixed-polarity magnetic fields. We suggest that the horizontal motion of the granulation in the photosphere drives the magnetic reconnection, which is modulated byp-mode oscillations.

     
    more » « less
  2. 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
  3. Abstract

    We investigate the thermal, kinematic, and magnetic structure of small-scale heating events in an emerging flux region (EFR). We use high-resolution multiline observations (including Caii8542 Å, CaiiK, and the Fei6301 Å line pair) of an EFR located close to the disk center from the CRISP and CHROMIS instruments at the Swedish 1 m Solar Telescope. We perform non-LTE inversions of multiple spectral lines to infer the temperature, velocity, and magnetic field structure of the heating events. Additionally, we use the data-driven Coronal Global Evolutionary Model to simulate the evolution of the 3D magnetic field configuration above the events and understand their dynamics. Furthermore, we analyze the differential emission measure to gain insights into the heating of the coronal plasma in the EFR. Our analysis reveals the presence of numerous small-scale heating events in the EFR, primarily located at polarity inversion lines of bipolar structures. These events not only heat the lower atmosphere but also significantly heat the corona. The data-driven simulations, along with the observed enhancement of currents and Poynting flux, suggest that magnetic reconnection in the lower atmosphere is likely responsible for the observed heating at these sites.

     
    more » « less
  4. ABSTRACT Discs are essential to the formation of both stars and planets, but how they form in magnetized molecular cloud cores remains debated. This work focuses on how the disc formation is affected by turbulence and ambipolar diffusion (AD), both separately and in combination, with an emphasis on the protostellar mass accretion phase of star formation. We find that a relatively strong, sonic turbulence on the core scale strongly warps but does not completely disrupt the well-known magnetically induced flattened pseudo-disc that dominates the inner protostellar accretion flow in the laminar case, in agreement with previous work. The turbulence enables the formation of a relatively large disc at early times with or without AD, but such a disc remains strongly magnetized and does not persist to the end of our simulation unless a relatively strong AD is also present. The AD-enabled discs in laminar simulations tend to fragment gravitationally. The disc fragmentation is suppressed by initial turbulence. The AD facilitates the disc formation and survival by reducing the field strength in the circumstellar region through magnetic flux redistribution and by making the field lines there less pinched azimuthally, especially at late times. We conclude that turbulence and AD complement each other in promoting disc formation. The discs formed in our simulations inherit a rather strong magnetic field from its parental core, with a typical plasma-β of order a few tens or smaller, which is 2–3 orders of magnitude lower than the values commonly adopted in magnetohydrodynamic simulations of protoplanetary discs. To resolve this potential tension, longer term simulations of disc formation and evolution with increasingly more realistic physics are needed. 
    more » « less
  5. null (Ed.)
    Context. Physical processes that govern the star and planet formation sequence influence the chemical composition and evolution of protoplanetary disks. Recent studies allude to an early start to planet formation already during the formation of a disk. To understand the chemical composition of protoplanets, we need to constrain the composition and structure of the disks from whence they are formed. Aims. We aim to determine the molecular abundance structure of the young disk around the TMC1A protostar on au scales in order to understand its chemical structure and any possible implications for disk formation. Methods. We present spatially resolved Atacama Large Millimeter/submillimeter Array observations of CO, HCO + , HCN, DCN, and SO line emission, as well as dust continuum emission, in the vicinity of TMC1A. Molecular column densities are estimated both under the assumption of optically thin emission from molecules in local thermodynamical equilibrium (LTE) as well as through more detailed non-LTE radiative transfer calculations. Results. Resolved dust continuum emission from the disk is detected between 220 and 260 GHz. Rotational transitions from HCO + , HCN, and SO are also detected from the inner 100 au region. We further report on upper limits to vibrational HCN υ 2 = 1, DCN, and N 2 D + lines. The HCO + emission appears to trace both the Keplerian disk and the surrounding infalling rotating envelope. HCN emission peaks toward the outflow cavity region connected with the CO disk wind and toward the red-shifted part of the Keplerian disk. From the derived HCO + abundance, we estimate the ionization fraction of the disk surface, and find values that imply that the accretion process is not driven by the magneto-rotational instability. The molecular abundances averaged over the TMC1A disk are similar to its protostellar envelope and other, older Class II disks. We meanwhile find a discrepancy between the young disk’s molecular abundances relative to Solar System objects. Conclusions. Abundance comparisons between the disk and its surrounding envelope for several molecular species reveal that the bulk of planet-forming material enters the disk unaltered. Differences in HCN and H 2 O molecular abundances between the disk around TMC1A, Class II disks, and Solar System objects trace the chemical evolution during disk and planet formation. 
    more » « less