skip to main content

Title: Data-driven Radiative Magnetohydrodynamics Simulations with the MURaM Code

We present a method of conducting data-driven simulations of solar active regions and flux emergence with the MURaM radiative magnetohydrodynamics (MHD) code. The horizontal electric field that is derived from the full velocity and magnetic vectors is implemented at the photospheric (bottom) boundary to drive the induction equation. The energy equation accounts for thermal conduction along magnetic fields, optically thin radiative loss, and heating of coronal plasma by viscous and resistive dissipation, which allows for a realistic presentation of the thermodynamic properties of coronal plasma that are key to predicting the observational features of solar active regions and eruptions. To validate this method, the photospheric data from a comprehensive radiative MHD simulation of solar eruption (the ground truth) are used to drive a series of numerical experiments. The data-driven simulation reproduces the accumulation of free magnetic energy over the course of flux emergence in the ground truth with an error of 3%. The onset time is approximately 8 minutes delayed compared to the ground truth. However, a precursor-like signature can be identified at the correct onset time. The data-driven simulation captures key eruption-related emission features and plasma dynamics of the ground truth flare over a wide temperature span, fromlog10T=4.5tolog10T>8. The evolution of the flare and coronal mass ejection as seen in synthetic extreme ultraviolet images is also reproduced with high fidelity. This method helps to understand the evolution of magnetic field in a more realistic coronal environment and to link the magnetic structures to observable diagnostics.

more » « less
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Medium: X Size: Article No. 118
["Article No. 118"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We present a comprehensive radiative magnetohydrodynamic simulation of the quiet Sun and large solar active regions. The 197 Mm wide simulation domain spans from 18(10) Mm beneath the photosphere to 113 Mm in the solar corona. Radiative transfer assuming local thermal equilibrium, optically thin radiative losses, and anisotropic conduction transport provide the necessary realism for synthesizing observables to compare with remote-sensing observations of the photosphere and corona. This model self-consistently reproduces observed features of the quiet Sun, emerging and developed active regions, and solar flares up to M class. Here, we report an overview of the first results. The surface magneto-convection yields an upward Poynting flux that is dissipated in the corona and heats the plasma to over 1 MK. The quiescent corona also presents ubiquitous propagating waves, jets, and bright points with sizes down to 2 Mm. Magnetic flux bundles emerge into the photosphere and give rise to strong and complex active regions with over 1023Mx magnetic flux. The coronal free magnetic energy, which is approximately 18% of the total magnetic energy, accumulates to approximately 1033erg. The coronal magnetic field is clearly non-force-free, as the Lorentz force needs to balance the pressure force and viscous stress as well as drive magnetic field evolution. The emission measure fromlog10T=4.5tolog10T>7provides a comprehensive view of the active region corona, such as coronal loops of various lengths and temperatures, mass circulation by evaporation and condensation, and eruptions from jets to large-scale mass ejections.

    more » « less
  2. Abstract

    We present a theory based on the conventional two-term (i.e. Lorentzian) approximation to the exact solution of the Boltzmann equation in non-magnetized weakly ionized plasma to efficiently obtain the electron rate and transport coefficients in a magnetized plasma for an arbitrary magnitude and direction of applied electric fieldEand magnetic fieldB. The proposed transcendental method does not require the two-term solution of the Boltzmann equation in magnetized plasma, based on which the transport parameters vary as a function of the reduced electric fieldE/N, reduced electron cyclotron frequencyωce/N, and angleE,BbetweenEandBvectors, whereNis the density of neutrals. Comparisons between the coefficients derived from BOLSIG+’s solution (obtained via the two-term expansion whenB0) and coefficients of the presented method are illustrated for air, a mixture of molecular hydrogen (H2) and helium (He) representing the giant gas planets of the Solar System, and pure carbon dioxide (CO2). The new approach may be used in the modeling of magnetized plasma encountered in the context of transient luminous events, e.g. sprite streamers in the atmosphere of Earth and Jupiter, in modeling the propagation of lightning’s electromagnetic pulses in Earth’s ionosphere, and in various laboratory and industrial applications of nonthermal plasmas.

    more » « less
  3. Abstract

    The best upper limit for the electron electric dipole moment was recently set by the ACME collaboration. This experiment measures an electron spin-precession in a cold beam of ThO molecules in their metastableH(3Δ1)state. Improvement in the statistical and systematic uncertainties is possible with more efficient use of molecules from the source and better magnetometry in the experiment, respectively. Here, we report measurements of several relevant properties of the long-livedQ(3Δ2)state of ThO, and show that this state is a very useful resource for both these purposes. TheQstate lifetime is long enough that its decay during the time of flight in the ACME beam experiment is negligible. The large electric dipole moment measured for theQstate, giving rise to a large linear Stark shift, is ideal for an electrostatic lens that increases the fraction of molecules detected downstream. The measured magnetic moment of theQstate is also large enough to be used as a sensitive co-magnetometer in ACME. Finally, we show that theQstate has a large transition dipole moment to theC(1Π1)state, which allows for efficient population transfer between the ground stateX(1Σ+)and theQstate viaXCQStimulated Raman Adiabatic Passage (STIRAP). We demonstrate 90 % STIRAP transfer efficiency. In the course of these measurements, we also determine the magnetic moment ofCstate, theXCtransition dipole moment, and branching ratios of decays from theCstate.

    more » « less
  4. Abstract

    The recent discovery of astrophysical neutrinos from the Seyfert galaxy NGC 1068 suggests the presence of nonthermal protons within a compact “coronal” region close to the central black hole. The acceleration mechanism of these nonthermal protons remains elusive. We show that a large-scale magnetic reconnection layer, of the order of a few gravitational radii, may provide such a mechanism. In such a scenario, rough energy equipartition between magnetic fields, X-ray photons, and nonthermal protons is established in the reconnection region. Motivated by recent 3D particle-in-cell simulations of relativistic reconnection, we assume that the spectrum of accelerated protons is a broken power law, with the break energy being constrained by energy conservation (i.e., the energy density of accelerated protons is at most comparable to the magnetic energy density). The proton spectrum isdnp/dEpEp1below the break anddnp/dEpEpsabove the break, with IceCube neutrino observations suggestings≃ 3. Protons above the break lose most of their energy within the reconnection layer via photohadronic collisions with the coronal X-rays, producing a neutrino signal in good agreement with the recent observations. Gamma rays injected in photohadronic collisions are cascaded to lower energies, sustaining the population of electron–positron pairs that makes the corona moderately Compton thick.

    more » « less
  5. Abstract

    Magnetic reconnection is often invoked as a source of high-energy particles, and in relativistic astrophysical systems it is regarded as a prime candidate for powering fast and bright flares. We present a novel analytical model—supported and benchmarked with large-scale three-dimensional kinetic particle-in-cell simulations in electron–positron plasmas—that elucidates the physics governing the generation of power-law energy spectra in relativistic reconnection. Particles with Lorentz factorγ≳ 3σ(here,σis the magnetization) gain most of their energy in the inflow region, while meandering between the two sides of the reconnection layer. Their acceleration time istaccγηrec1ωc120γωc1, whereηrec≃ 0.06 is the inflow speed in units of the speed of light andωc=eB0/mcis the gyrofrequency in the upstream magnetic field. They leave the region of active energization aftertesc, when they get captured by one of the outflowing flux ropes of reconnected plasma. We directly measuretescin our simulations and find thattesctaccforσ≳ few. This leads to a universal (i.e.,σ-independent) power-law spectrumdNfree/dγγ1for the particles undergoing active acceleration, anddN/dγγ2for the overall particle population. Our results help to shed light on the ubiquitous presence of power-law particle and photon spectra in astrophysical nonthermal sources.

    more » « less