More Like this

1. Data-driven Radiative Magnetohydrodynamics Simulations with the MURaM Code

   https://doi.org/10.3847/1538-4357/acc8c5  

   Chen, Feng ; Cheung, Mark C. M. ; Rempel, Matthias ; Chintzoglou, Georgios ( June 2023 , The Astrophysical Journal)

   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, from${\log }_{10}T=4.5$to${\log }_{10}T>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
2. Analytical Model of Disk Evaporation and State Transitions in Accreting Black Holes

   https://doi.org/10.3847/1538-4357/ac6d5c  

   Cho 조, Hyerin 혜린 ; Narayan, Ramesh ( June 2022 , The Astrophysical Journal)

   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${\dot{m}}_c$to 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 radiusr_b∼ 10³(α/0.3)⁻², whereαis the viscosity parameter. Gas evaporates from the disk to the corona forr>r_b, and condenses back forr<r_b. Atr_b,${\dot{m}}_c$reaches its maximum,${\dot{m}}_{c,\max }\approx 0.02{(\alpha /0.3)}^3$. If atr≫r_bthe thin disk accretes with${\dot{m}}_d<{\dot{m}}_{c,\max }$, then the disk evaporates fully before reachingr_b, 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, andr_bremains unchanged. This model predicts strong coronal winds forr>r_b, and aT∼ 5 × 10⁸K Compton-cooled corona forr<r_b. Two-temperature effects are ignored, but may be important at small radii.

    

   more » « less
3. Improving Coronal Hole Detections and Open Flux Estimates

   https://doi.org/10.3847/1538-4357/ad01b6  

   Caplan, Ronald M. ; Mason, Emily I. ; Downs, Cooper ; Linker, Jon A. ( November 2023 , The Astrophysical Journal)

   One systematic limitation of solar coronal hole (CH) detection at extreme ultraviolet (EUV) wavelengths is the obscuration of dark regions of the corona by brighter structures along the line of sight. Another problem arises when using CHs to compute the Sun’s open magnetic flux, where surface measurements of the radial magnetic field,$B_r^{\odot }$, are situated slightly below the effective height of coronal EUV emission. In this paper, we explore these two limitations utilizing a thermodynamic magnetohydrodynamic (MHD) model of the corona for Carrington rotation (CR) 2101, where we generate CH detections from EUV 193 Å images of the corona forward-modeled from the MHD solution, and where the modeled open field is known. We demonstrate a method to combine EUV images into a full Sun map that helps alleviate CH obscuration called theminimum intensity diskmerge(MIDM). We also show the variation in measured open flux and CH area that is due to the effective height differences between EUV and$B_r^{\odot }$measurements. We then apply the MIDM method to SDO/AIA 193 Å observations from CR 2101, and conduct an analogous analysis. In this case, the MIDM method uses time-varying images, the effects of which are discussed. We show that overall, the MIDM method and an appreciation of the effective height mismatch provide a useful new way to extract a broader view of CHs, especially near the poles. In turn, they enable improved estimates of the open magnetic flux, and help facilitate comparisons between models and observations.

    

   more » « less
4. Kinetic-scale Current Sheets in the Solar Wind at 1 au: Scale-dependent Properties and Critical Current Density

   https://doi.org/10.3847/2041-8213/ac4fc4  

   Vasko, I. Y. ; Alimov, K. ; Phan, T. ; Bale, S. D. ; Mozer, F. S. ; Artemyev, A. V. ( February 2022 , The Astrophysical Journal Letters)

   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⁻¹magnetic 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,J₀≈ 6 nAm⁻²· (λ/100 km)^−0.56, but does not statistically exceed a critical value,J_A,corresponding to the drift between ions and electrons of local Alvén speed. The observed trend holds in normalized units:$J_0/J_A\approx 0.17·{(\lambda /{\lambda }_p)}^{-0.51}$. The CSs are statistically force-free with magnetic shear angle correlated with CS spatial scale:$\mathrm{\Delta }\theta \approx 19°·{(\lambda /{\lambda }_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.

    

   more » « less
5. Limits on Astrophysical Antineutrinos with the KamLAND Experiment

   https://doi.org/10.3847/1538-4357/ac32c1  

   Abe, S. ; Asami, S. ; Gando, A. ; Gando, Y. ; Gima, T. ; Goto, A. ; Hachiya, T. ; Hata, K. ; Hayashida, S. ; Hosokawa, K. ; et al ( January 2022 , The Astrophysical Journal)

   We report on a search for electron antineutrinos (${\overline{\nu }}_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⁻²s⁻¹(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 of⁸B solar neutrinos converting into${\overline{\nu }}_e$,$P_{{\nu }_e\to {\overline{\nu }}_e}<3.5\times {10}^{-5}$(90% CL) assuming an undistorted${\overline{\nu }}_e$shape. This corresponds to a solar${\overline{\nu }}_e$flux of 60 cm⁻²s⁻¹(90% CL) in the analysis energy range.

    

   more » « less