skip to main content

Title: A Comprehensive Radiative Magnetohydrodynamics Simulation of Active Region Scale Flux Emergence from the Convection Zone to the Corona

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 more » 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.

« less
; ;
Publication Date:
Journal Name:
The Astrophysical Journal
Page Range or eLocation-ID:
Article No. 91
DOI PREFIX: 10.3847
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    State transitions in black hole X-ray binaries are likely caused by gas evaporation from a thin accretion disk into a hot corona. We present a height-integrated version of this process, which is suitable for analytical and numerical studies. With radiusrscaled to Schwarzschild units and coronal mass accretion rateṁcto Eddington units, the results of the model are independent of black hole mass. State transitions should thus be similar in X-ray binaries and an active galactic nucleus. The corona solution consists of two power-law segments separated at a break radiusrb∼ 103(α/0.3)−2, whereαis the viscosity parameter. Gas evaporates from the disk to the corona forr>rb, and condenses back forr<rb. Atrb,ṁcreaches its maximum,ṁc,max0.02(α/0.3)3. If atrrbthe thin disk accretes withṁd<ṁc,max, then the disk evaporates fully before reachingrb, giving the hard state. Otherwise, the disk survives at all radii, giving the thermal state. While the basic model considers only bremsstrahlung cooling and viscous heating, we also discuss a more realistic model that includes Compton cooling and direct coronal heating by energy transport from the disk. Solutions are again independent of black hole mass, andrbremainsmore »unchanged. This model predicts strong coronal winds forr>rb, and aT∼ 5 × 108K Compton-cooled corona forr<rb. Two-temperature effects are ignored, but may be important at small radii.

    « less
  2. Abstract

    We present analysis of 17,043 proton kinetic-scale current sheets (CSs) collected over 124 days of Wind spacecraft measurements in the solar wind at 11 samples s−1magnetic field resolution. The CSs have thickness,λ,from a few tens to one thousand kilometers with typical values around 100 km, or within about 0.1–10λpin terms of local proton inertial length,λp. We found that the current density is larger for smaller-scale CSs,J0≈ 6 nAm−2· (λ/100 km)−0.56, but does not statistically exceed a critical value,JA,corresponding to the drift between ions and electrons of local Alvén speed. The observed trend holds in normalized units:J0/JA0.17·(λ/λp)0.51. The CSs are statistically force-free with magnetic shear angle correlated with CS spatial scale:Δθ19°·(λ/λp)0.5. The observed correlations are consistent with local turbulence being the source of proton kinetic-scale CSs in the solar wind, while the mechanisms limiting the current density remain to be understood.

  3. Abstract

    We report on a search for electron antineutrinos (ν¯e) from astrophysical sources in the neutrino energy range 8.3–30.8 MeV with the KamLAND detector. In an exposure of 6.72 kton-year of the liquid scintillator, we observe 18 candidate events via the inverse beta decay reaction. Although there is a large background uncertainty from neutral current atmospheric neutrino interactions, we find no significant excess over background model predictions. Assuming several supernova relic neutrino spectra, we give upper flux limits of 60–110 cm−2s−1(90% confidence level, CL) in the analysis range and present a model-independent flux. We also set limits on the annihilation rates for light dark matter pairs to neutrino pairs. These data improve on the upper probability limit of8B solar neutrinos converting intoν¯e,Pνeν¯e<3.5×105(90% CL) assuming an undistortedν¯eshape. This corresponds to a solarν¯eflux of 60 cm−2s−1(90% CL) in the analysis energy range.

  4. Abstract

    We present a toy model for the thermal optical/UV/X-ray emission from tidal disruption events (TDEs). Motivated by recent hydrodynamical simulations, we assume that the debris streams promptly and rapidly circularize (on the orbital period of the most tightly bound debris), generating a hot quasi-spherical pressure-supported envelope of radiusRv∼ 1014cm (photosphere radius ∼1015cm) surrounding the supermassive black hole (SMBH). As the envelope cools radiatively, it undergoes Kelvin–Helmholtz contractionRvt−1, its temperature risingTefft1/2while its total luminosity remains roughly constant; the optical luminosity decays asνLνRv2Tefft3/2. Despite this similarity to the mass fallback rateṀfbt5/3, envelope heating from fallback accretion is subdominant compared to the envelope cooling luminosity except near optical peak (where they are comparable). Envelope contraction can be delayed by energy injection from accretion from the inner envelope onto the SMBH in a regulated manner, leading to a late-time flattening of the optical/X-ray light curves, similar to those observed in some TDEs. Eventually, as the envelope contracts to near the circularization radius, the SMBH accretion rate rises to its maximum, in tandem with the decreasing optical luminosity. This cooling-induced (rather than circularization-induced) delay of up to several hundred days may account for themore »delayed onset of thermal X-rays, late-time radio flares, and high-energy neutrino generation, observed in some TDEs. We compare the model predictions to recent TDE light-curve correlation studies, finding both agreement and points of tension.

    « less
  5. Abstract

    The Parker Solar Probe (PSP) entered a region of sub-Alfvénic solar wind during encounter 8, and we present the first detailed analysis of low-frequency turbulence properties in this novel region. The magnetic field and flow velocity vectors were highly aligned during this interval. By constructing spectrograms of the normalized magnetic helicity, cross-helicity, and residual energy, we find that PSP observed primarily Alfvénic fluctuations, a consequence of the highly field-aligned flow that renders quasi-2D fluctuations unobservable to PSP. We extend Taylor’s hypothesis to sub- and super-Alfvénic flows. Spectra for the fluctuating forward and backward Elsässer variables (z±, respectively) are presented, showing thatz+modes dominatezby an order of magnitude or more, and thez+spectrum is a power law in frequency (parallel wavenumber)f−3/2(k3/2) compared to the convexzspectrum withf−3/2(k3/2) at low frequencies, flattening around a transition frequency (at which the nonlinear and Alfvén timescales are balanced) tof−1.25at higher frequencies. The observed spectra are well fitted using a spectral theory for nearly incompressible magnetohydrodynamics assuming a wavenumber anisotropykk3/4, that thez+fluctuations experience primarily nonlinear interactions, and that the minorityzfluctuations experience both nonlinear and Alfvénic interactions withz+fluctuations. The density spectrum is a powermore »law that resembles neither thez±spectra nor the compressible magnetic field spectrum, suggesting that these are advected entropic rather than magnetosonic modes and not due to the parametric decay instability. Spectra in the neighboring modestly super-Alfvénic intervals are similar.

    « less