🔭 astrophysics

Mg II h&k spectral line properties computed using 3D radiative transfer in an enhanced network region simulated with the MURaM-ChE code

This study utilizes 3D radiative transfer simulations of an enhanced network region with the MURaM-ChE code to demonstrate that including full 3D effects, particularly horizontal velocities, significantly improves the agreement between modeled and observed Mg II h&k spectral line properties compared to 1.5D approximations.

Original authors: P. Ondratschek, D. Przybylski, H. N. Smitha, J. Leenaarts, R. Cameron, S. K. Solanki

Published 2026-05-25

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: P. Ondratschek, D. Przybylski, H. N. Smitha, J. Leenaarts, R. Cameron, S. K. Solanki

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the Sun's atmosphere as a bustling, multi-layered city. Deep down, near the surface, it's a hot, dense crowd. Higher up, in the "chromosphere" (the middle layer), things get chaotic, with gas moving up and down like a stormy sea, and magnetic fields acting like invisible highways guiding the flow.

Scientists want to understand this city, but they can't send a probe there. Instead, they look at the light coming from the Sun. Specifically, they look at two very bright "streetlights" in the Sun's spectrum called Mg ii h&k lines. These lights are formed in the middle of the chromosphere and hold clues about the temperature and speed of the gas around them.

For a long time, scientists tried to understand these lights by looking at the Sun as if it were a flat, calm pancake. They would take a vertical slice of the atmosphere and calculate how light travels through it, assuming that the gas to the left and right didn't matter. This is like trying to understand the wind in a city by only looking at one single street and ignoring the wind blowing in from the neighboring avenues.

The New Experiment: A 3D City Map
In this paper, the authors decided to stop pretending the Sun is a flat pancake. They used a super-computer simulation called MURaM-ChE to create a realistic, 3D model of a "network region" on the Sun (a place where magnetic fields are strong). This model is like a full 3D video game map of the solar atmosphere, complete with swirling winds, hot and cold pockets, and magnetic loops.

They then ran a new calculation to see how light travels through this 3D model. This is called 3D Radiative Transfer (RT). They compared this to the old "flat pancake" method (1.5D RT) and to actual photos taken by a satellite called IRIS.

What They Found: The "Side-Wind" Effect
The results were surprising and important:

The Flat Model Overestimates the Brightness: When the scientists used the old "flat pancake" method, the simulated streetlights (the Mg ii lines) were way too bright. It was like their model predicted the city was glowing with the intensity of a thousand suns, while the real satellite photos showed a more modest glow.
The 3D Model is More Accurate: When they used the full 3D calculation, the simulated lights dimmed down and matched the real satellite photos much better.
The Secret Ingredient: Horizontal Winds: Why did the 3D model work better? The authors discovered that horizontal winds (gas moving sideways, not just up and down) play a huge role. In the old flat model, these side-winds were ignored. In the 3D model, the sideways movement of the gas acts like a "smear" that spreads the light out, making the peaks of the spectral lines less intense and the overall picture more realistic. It's like how a strong crosswind on a highway blurs the headlights of cars, making them look less sharp and bright than they would in still air.

New Clues from the Light
The paper also dug deeper into what the different parts of the light tell us:

Blue vs. Red Peaks: The Mg ii lines have two peaks, one on the blue side and one on the red side. The authors found that the blue peak tends to form in gas that is rushing upward, while the red peak forms in gas that is sinking downward. It's as if the blue light is a flag waving from a rising balloon, and the red light is a flag from a falling one.
Temperature Check: They tried to use the brightness of these lights to guess the temperature of the gas. They found that while it works okay for the brightest parts of the light, it gets messy for the dimmer parts because the gas is so thin and chaotic that the light stops behaving like a simple thermometer.

The Bottom Line
This paper is a victory for "thinking in 3D." It shows that to understand the Sun's atmosphere, we can't just look at vertical slices; we have to account for the sideways chaos. By including the full 3D effects and the horizontal winds, the scientists created a simulation that looks much more like the real Sun. This helps them interpret the data from satellites like IRIS more accurately, allowing us to better understand the invisible storms and magnetic fields that drive our star.

Technical Summary: Mg ii h&k Spectral Line Properties Computed Using 3D Radiative Transfer in an Enhanced Network Region

Problem and Context
The solar chromosphere is a complex, inhomogeneous interface between the photosphere and corona, where spectral line formation occurs under non-local thermodynamic equilibrium (NLTE) conditions. Strong resonance lines, such as Mg ii h&k, exhibit partial frequency redistribution (PRD) effects. While 3D radiative magnetohydrodynamics (rMHD) simulations provide realistic atmospheric models, the synthesis of spectral lines from these models often relies on 1.5D radiative transfer (RT) approximations. In 1.5D RT, each vertical column is treated as an independent plane-parallel atmosphere, neglecting horizontal radiative transfer and horizontal velocity fields. Previous studies using the public Bifrost simulation snapshot indicated discrepancies between 1.5D and 3D RT results, but the magnitude of these effects in more dynamic, higher-resolution models remained to be fully quantified. Specifically, the formation details of Mg ii h&k lines in enhanced network (EN) regions and the robustness of established atmospheric diagnostics in such dynamic environments require further verification.

Methodology
The authors utilized a self-consistent 3D rMHD simulation of an EN region generated with the MURaM-ChE code (the chromospheric extension of MURaM). This model features a bipolar magnetic structure with a resolution of 23.4 km horizontally and 20 km vertically, incorporating NLTE energy transport and non-equilibrium hydrogen ionization.

The study employed two primary synthesis approaches:

3D Radiative Transfer: Spectra were computed using the Multi3D code, solving the radiative transfer equation (RTE) and statistical equilibrium equations simultaneously in full 3D geometry with PRD treatment.
1.5D Radiative Transfer: For comparison, spectra were computed from the same atmospheric snapshot using the plane-parallel approximation (1.5D RT) within Multi3D, neglecting horizontal radiative coupling and horizontal velocities.

Synthetic spectra were degraded to match the instrumental conditions of the Interface Region Imaging Spectrograph (IRIS), including spatial and spectral convolution. These were compared against:

IRIS observations of a bipolar feature near disk center.
Previous results from the public Bifrost snapshot.
Controlled experiments isolating the effects of angular quadrature and horizontal velocities.

Key Contributions and Results

Impact of 3D RT on Line Profiles: The spatially averaged Mg ii h&k spectra computed with 3D RT match IRIS observations significantly better in the core intensities and their distribution than 1.5D RT computations. However, the 1.5D approach significantly overestimates the peak intensities (k2v and k2r) and the central minimum (k3). In the MURaM-ChE model, the discrepancy between 1.5D and 3D RT is more pronounced than in the public Bifrost snapshot; 1.5D intensities are up to 72% higher than 3D results in the line core, whereas the difference was approximately 12% in previous Bifrost studies.
Role of Horizontal Velocities: A controlled experiment revealed that horizontal velocities, naturally included in full 3D RT but absent in standard 1.5D RT, significantly broaden the line profiles and reduce peak intensities. The inclusion of horizontal velocities in a plane-parallel calculation brings the result closer to the full 3D solution, suggesting that the dynamic horizontal velocity field in the MURaM-ChE model is a primary driver of the large discrepancy between 1.5D and 3D results.
Formation Height and Velocity Correlations: The study confirmed that the k2v (blue) and k2r (red) peaks form at different atmospheric heights and trace distinct flow structures:
- k2v: Forms preferentially in upflows at lower atmospheric heights (higher densities, ~6.8 kK).
- k2r: Forms preferentially in downflows at higher atmospheric heights (lower densities, ~8 kK).
- k3: The Doppler shift of the k3 minimum shows a tight correlation (Pearson coefficient $R_P = 0.94$ ) with the vertical velocity at its formation height, confirming its utility as a velocity diagnostic, though the scatter is larger than in the Bifrost snapshot due to the more dynamic MURaM-ChE atmosphere.
Temperature Diagnostics: The correlation between peak brightness temperature and local atmospheric temperature was found to be weak ( $R_P \approx 0.12$ for k2v and $0.06$ for k2r) for the MURaM-ChE model. While a correlation exists for intensities above ~5.25 kK, the source function is often decoupled from the local gas temperature, limiting the reliability of Mg ii h&k peaks as direct temperature diagnostics in this highly dynamic model.
Statistical Distributions: The distribution of peak brightness temperatures and peak separations in the 3D RT computation aligns more closely with IRIS observations than the 1.5D results. However, the 3D RT simulation still underestimates the peak separation compared to observations, potentially indicating that the simulation's vertical motions are lower than those on the real Sun or that the simulation box size limits the interaction between network elements and the quiet Sun.

Significance and Conclusions
The paper concludes that 3D radiative transfer is essential for the accurate forward modeling of Mg ii h&k lines, particularly in dynamic chromospheric regions like enhanced networks. The study demonstrates that neglecting horizontal radiative transfer and horizontal velocities leads to a significant overestimation of line core intensities and a misrepresentation of the line profile shape.

The authors emphasize that the MURaM-ChE model, despite moderate resolution, produces a more dynamic atmosphere than the public Bifrost snapshot, resulting in more complex spectral profiles and larger discrepancies between 1.5D and 3D RT. While the model successfully reproduces observed line widths and general intensity distributions, remaining discrepancies in peak separation and intensity ratios suggest that further improvements in simulation resolution and horizontal domain size may be necessary to fully capture the mass loading and dynamics of the solar chromosphere. The work reinforces the diagnostic potential of Mg ii h&k lines for constraining velocity structures but cautions against using them as direct temperature diagnostics in highly dynamic, NLTE-dominated environments without accounting for the complex formation physics revealed by 3D RT.

Technical Summary: Mg ii h&k Spectral Line Properties Computed Using 3D Radiative Transfer in an Enhanced Network Region

More like this