pith. sign in

arxiv: 2605.22916 · v1 · pith:ADNXLZDLnew · submitted 2026-05-21 · 🌌 astro-ph.SR

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

P. Ondratschek , D. Przybylski , H.N. Smitha , J. Leenaarts , R. Cameron , S.K. Solanki This is my paper

Pith reviewed 2026-05-25 02:12 UTC · model grok-4.3

classification 🌌 astro-ph.SR
keywords Mg II h&k lines3D radiative transfersolar chromosphereMURaM-ChE simulationenhanced networkIRIS observationspartial frequency redistribution
0
0 comments X

The pith

3D radiative transfer produces Mg II h&k core intensities and distributions that match solar observations more closely than 1.5D methods.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper forward-models the Mg II h and k lines using full 3D radiative transfer with partial frequency redistribution inside a 3D radiative MHD simulation of an enhanced network region that includes non-equilibrium hydrogen ionization. The resulting spatially averaged line profiles are compared to IRIS data, and the 3D results are shown to reproduce observed core intensities and their spatial distribution on the Sun better than standard 1.5D calculations performed on the same atmosphere. The larger 1.5D–3D difference relative to earlier models is partly traced to horizontal velocities that appear only in the full 3D synthesis. Correlations between line features and atmospheric quantities remain qualitatively similar to those found in prior Bifrost snapshots, but exhibit greater scatter because the present simulation is more dynamic.

Core claim

The Mg II h&k lines computed with 3D RT match the observations better in the core intensities and their distribution on the Sun compared to 1.5D computations. This underlines the importance of 3D RT in the forward modeling of Mg II h&k.

What carries the argument

Full 3D radiative transfer with partial frequency redistribution applied to a self-consistent MURaM-ChE rMHD snapshot that includes NLTE energy transport and non-equilibrium hydrogen ionization.

If this is right

  • Spatially averaged Mg II h&k profiles from 3D RT approximate a typical IRIS observation, though peak separation remains slightly too small.
  • Horizontal velocities naturally present in 3D RT contribute substantially to the difference from 1.5D results.
  • Line–atmosphere correlations remain similar to those in earlier Bifrost models but display more scatter due to greater atmospheric dynamics.
  • The qualitative gap between 1.5D and 3D RT is larger in this MURaM-ChE snapshot than reported for public Bifrost models.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Accurate forward modeling of chromospheric lines in dynamic regions may require full 3D RT whenever horizontal flows are significant.
  • Improved peak-separation matches could be tested by varying the treatment of velocity fields or by running the same RT code on multiple independent simulations.
  • The greater scatter in correlations suggests that statistical relations derived from 3D models may need larger sample sizes to be robust for diagnostic use.

Load-bearing premise

The MURaM-ChE simulation accurately represents real solar chromospheric conditions including non-equilibrium ionization and NLTE effects.

What would settle it

If new observations of Mg II h&k core intensities in enhanced network regions are shown to agree more closely with 1.5D RT than with 3D RT on the same atmospheric model, the claimed superiority of 3D RT would be contradicted.

Figures

Figures reproduced from arXiv: 2605.22916 by D. Przybylski, H.N. Smitha, J. Leenaarts, P. Ondratschek, R. Cameron, S.K. Solanki.

Figure 1
Figure 1. Figure 1: Spatially averaged spectra of the Mg ii h&k lines. Shown are the spectra from the MURaM-ChE simulation once computed in full 3D RT (red) and once using the plane-parallel 1.5D RT approximation (grey). For compari￾son, we show an observation from a qualita￾tively similar bipolar feature, which covers a similar area on the Sun (black). z = 0 Mm lies at the average τ500 = 1 height. The atmosphere extends up t… view at source ↗
Figure 3
Figure 3. Figure 3: Qualitative comparison with observations. We compare three different sets of intensity images. The first row shows the intensity map of the Mg ii k line from the observation taken at the line-center rest wavelength (panel a), the k2v feature (panel b), the k3 feature (panel c), and the k2r feature (panel d). The second row (panels e,f,g, and h) shows the same quantities but for the MURaM-ChE model using 3D… view at source ↗
Figure 4
Figure 4. Figure 4: Statistical comparison with observations. In panel (a) we show distributions of the k2 peak brightness temperature, in panel (b) the dis￾tribution of k2 peak separation, and in panel (c) the distribution of the k2 peak intensity ratio. We show data from the IRIS observation in black and data from the synthetic spectra computed by 3D RT in red, and 1.5D RT in grey after degradation to instrumental condition… view at source ↗
Figure 5
Figure 5. Figure 5: Atmospheric properties at the formation height of the Mg ii k spectral line features. Intensity maps (first row, panels a–c), formation heights (second row, panels d–f), vertical velocity maps (third row, panels g–i), temperature (fourth row, panels j–l), and density (fifth row, panels m–o). We show these quantities for the k2v (left column, panels a,d,g,j, and m), k3 (middle column, panels b,e,h,k, and n)… view at source ↗
Figure 6
Figure 6. Figure 6: Correlations between Doppler shift of spectral features of Mg ii k and the vertical velocity in the atmosphere. Panel (a): Correlation between Doppler shift of k2v and the vertical component of the velocity at the formation height. Panel (b) shows a similar relation, but for k3. Panel (c) is the same as panel (a) but for k2r. The Pearson correlation coefficient RP is given in each panel. The blue solid lin… view at source ↗
Figure 7
Figure 7. Figure 7: Correlation between peak intensity ratio and average vertical velocity. The average vertical velocity is measured from the minimum formation height of the two k2 features and the formation height of the k3 feature. The red, green, and blue contours enclose 25%, 50%, and 90% of the data. The blue solid lines indicate “x = 0” and “y = 0”. chromosphere where the source function is still sufficiently cou￾pled … view at source ↗
Figure 8
Figure 8. Figure 8: Correlation between peak brightness temperatures and temperature at the formation height in the atmosphere. Panel (a) shows the correlation for the blue peak (k2v) and panel (b) shows the correlation for the red peak (k2r). The red, green, and blue contours enclose regions of 25%, 50%, and 90% of the data. The orange lines indicate “x = y”. We normalized each column of brightness temperature to the maximum… view at source ↗
Figure 9
Figure 9. Figure 9: Time-dependent velocity structure in the atmosphere. Panel (a) shows the average velocity as a function of height. Panel (b) shows the time-dependent average height of the τ500 = 1 surface. Panel (c) shows the average vertical velocity between z = 1.7 Mm and z = 2.7 Mm, which are the average formation heights of the k2v feature and k3 feature in the here presented snapshot. the observation, our results dem… view at source ↗
read the original abstract

The Mg II h&k lines form in the middle to upper chromosphere and are well-suited to study the structure of the chromosphere. However, the details of their formation in the solar chromosphere are not fully understood. We aim to study the effects of 3D radiative transfer (RT) on the Mg II h&k line properties and to verify known correlations between the underlying atmosphere and spectral line features in a new model of the chromosphere. We forward model the Mg II h&k lines in 3D RT with partial frequency redistribution (PRD) in a self-consistent 3D radiative magnetohydrodynamics (rMHD) simulation with non-local-thermodynamic-equilibrium (NLTE) energy transport and non-equilibrium (NE) hydrogen ionization of an enhanced network (EN) region simulated with the chromospheric extension of MURaM (MURaM-ChE). The spatially averaged Mg II h&k spectral lines computed with 3D RT match approximately a typical IRIS observation. The peak separation is still slightly lower in the simulation. In the MURaM-ChE model, the qualitative difference between 1.5D and 3D RT results is even more pronounced than in the public Bifrost snapshot, as given in the literature. We found that this large discrepancy might partly be attributed to the horizontal velocities that are naturally included in the full 3D RT synthesis but not in typical 1.5D RT computations. We confirm that correlations between spectral line properties and the underlying atmosphere from the MURaM-ChE simulation are similar to those obtained from Bifrost, but show more scatter due to the more dynamic atmosphere. The Mg II h&k lines computed with 3D RT match the observations better in the core intensities and their distribution on the Sun compared to 1.5D computations. This underlines the importance of 3D RT in the forward modeling of Mg II h&k.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The paper forward-models the Mg II h&k lines using 3D radiative transfer with PRD in a MURaM-ChE enhanced-network rMHD snapshot that includes NLTE energy transport and non-equilibrium hydrogen ionization. It reports that the spatially averaged 3D profiles match a typical IRIS observation more closely than 1.5D results (especially in core intensities and spatial distribution), that the 1.5D–3D discrepancy is larger than previously found with Bifrost, and that part of the difference may arise from horizontal velocities naturally present in 3D but omitted in standard 1.5D calculations. Correlations between line properties and atmospheric quantities are shown to be qualitatively similar to Bifrost but with greater scatter.

Significance. If the MURaM-ChE chromospheric stratification and dynamics are representative, the work supplies direct numerical evidence that 3D RT effects are essential for accurate Mg II forward modeling and can explain observed core properties better than 1.5D approximations. The consistency check against both IRIS data and an independent Bifrost snapshot, together with the use of a self-consistent NE-ionization simulation, adds concrete support for the importance of 3D RT in chromospheric diagnostics.

major comments (2)
  1. [Abstract] Abstract: the statement that the larger 1.5D–3D discrepancy “might partly be attributed to the horizontal velocities” is presented without any quantitative test, isolation procedure, or metric (e.g., no comparison run with horizontal velocities suppressed or any reported contribution fraction). Because this attribution is invoked to explain why the 3D improvement is more pronounced than in Bifrost, the lack of supporting analysis makes the causal claim load-bearing yet unsupported.
  2. [Abstract] Abstract: the residual mismatch in peak separation is noted but not accompanied by any diagnostic (temperature, velocity, or density stratification comparison) that would indicate whether the discrepancy originates in the MURaM-ChE model itself; without such diagnostics it is unclear whether the reported 3D advantage is general or specific to this simulation’s velocity field.
minor comments (1)
  1. [Abstract] The abstract would be strengthened by inclusion of at least one quantitative metric (mean absolute difference, Kolmogorov–Smirnov statistic, or similar) for the core-intensity distributions rather than the qualitative statement “match approximately.”

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments. We respond to each major comment below and indicate where revisions to the manuscript will be made.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the statement that the larger 1.5D–3D discrepancy “might partly be attributed to the horizontal velocities” is presented without any quantitative test, isolation procedure, or metric (e.g., no comparison run with horizontal velocities suppressed or any reported contribution fraction). Because this attribution is invoked to explain why the 3D improvement is more pronounced than in Bifrost, the lack of supporting analysis makes the causal claim load-bearing yet unsupported.

    Authors: We agree that the manuscript offers no dedicated quantitative test isolating the contribution of horizontal velocities (such as a suppressed-velocity 3D run). The statement rests on the fact that standard 1.5D calculations omit horizontal velocities by construction while full 3D RT includes them, together with the observation that the 1.5D–3D difference is larger in the MURaM-ChE snapshot than reported for Bifrost. We will revise the abstract to present the horizontal-velocity contribution as a plausible hypothesis rather than a firm attribution, and we will add a short clarifying sentence in the discussion section. No new numerical experiments are planned at this stage. revision: partial

  2. Referee: [Abstract] Abstract: the residual mismatch in peak separation is noted but not accompanied by any diagnostic (temperature, velocity, or density stratification comparison) that would indicate whether the discrepancy originates in the MURaM-ChE model itself; without such diagnostics it is unclear whether the reported 3D advantage is general or specific to this simulation’s velocity field.

    Authors: The abstract already notes that peak separation remains slightly lower than in the IRIS reference profile. The manuscript does not include direct comparisons of temperature, velocity or density stratifications against observations that would allow us to attribute the residual mismatch to specific model properties. We will revise the abstract and add a brief paragraph in the discussion to state that the 3D advantage demonstrated here is specific to the MURaM-ChE enhanced-network snapshot and its velocity field, and that further diagnostics would be needed to assess generality. revision: partial

Circularity Check

0 steps flagged

No significant circularity; forward modeling compared to external observations

full rationale

The paper performs direct numerical 3D RT synthesis (with PRD) on the MURaM-ChE rMHD snapshot and compares the resulting Mg II h&k profiles and statistics to independent IRIS observations plus literature Bifrost results. No parameters are fitted to the target data, no predictions reduce to inputs by construction, and no self-citation chain is load-bearing for the central claim. The external observational benchmark keeps the derivation self-contained.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The work rests on standard radiative transfer assumptions and the domain validity of the rMHD simulation; no free parameters are fitted to the spectral data itself and no new entities are postulated.

axioms (2)
  • standard math Partial frequency redistribution (PRD) treatment is required and correctly implemented for Mg II line formation in the chromosphere.
    Invoked for the 3D RT synthesis step.
  • domain assumption The MURaM-ChE rMHD simulation with NLTE energy transport and NE hydrogen ionization produces a physically realistic enhanced network chromosphere.
    Central premise enabling the forward modeling and 3D vs 1.5D comparison.

pith-pipeline@v0.9.0 · 5926 in / 1683 out tokens · 31031 ms · 2026-05-25T02:12:07.841915+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

33 extracted references · 33 canonical work pages

  1. [1]

    2013, ApJ, 779, 155 Bjørgen, J

    Avrett, E., Landi, E., & McKillop, S. 2013, ApJ, 779, 155 Bjørgen, J. P., Leenaarts, J., Rempel, M., et al. 2019, A&A, 631, A33 Bjørgen, J. P., Sukhorukov, A. V ., Leenaarts, J., et al. 2018, A&A, 611, A62

    work page 2013

  2. [2]

    M., Decaudin, M., Bruner, Jr., E

    Bonnet, R. M., Decaudin, M., Bruner, Jr., E. C., Acton, L. W., & Brown, W. A. 1980, ApJ, 237, L47

    work page 1980

  3. [3]

    1963, Meth

    Carlson, B. 1963, Meth. in Comp. Phys

    work page 1963

  4. [4]

    Carlsson, M., De Pontieu, B., & Hansteen, V . H. 2019, ARA&A, 57, 189

    work page 2019

  5. [5]

    H., Gudiksen, B

    Carlsson, M., Hansteen, V . H., Gudiksen, B. V ., Leenaarts, J., & De Pontieu, B. 2016, A&A, 585, A4

    work page 2016

  6. [6]

    & Leenaarts, J

    Carlsson, M. & Leenaarts, J. 2012, A&A, 539, A39

    work page 2012

  7. [7]

    & Stein, R

    Carlsson, M. & Stein, R. F. 1997, ApJ, 481, 500 de la Cruz Rodríguez, J. & van Noort, M. 2017, Space Sci. Rev., 210, 109 De Pontieu, B., Title, A. M., Lemen, J. R., et al. 2014, Sol. Phys., 289, 2733 del Toro Iniesta, J. C. & Ruiz Cobo, B. 2016, Liv. Rev. in Sol. Phys., 13, 4 Goši´c, M., de la Cruz Rodríguez, J., De Pontieu, B., et al. 2018, ApJ, 857, 48

    work page 1997

  8. [8]

    V ., Carlsson, M., Hansteen, V

    Gudiksen, B. V ., Carlsson, M., Hansteen, V . H., et al. 2011, A&A, 531, A154

    work page 2011

  9. [9]

    G., Kleint, L., Krucker, S., et al

    Hannah, I. G., Kleint, L., Krucker, S., et al. 2019, ApJ, 881, 109

    work page 2019

  10. [10]

    H., Martinez-Sykora, J., Carlsson, M., et al

    Hansteen, V . H., Martinez-Sykora, J., Carlsson, M., et al. 2023, ApJ, 944, 131 Jaume Bestard, J., Trujillo Bueno, J., Št ˇepán, J., & del Pino Alemán, T. 2021, ApJ, 909, 183

    work page 2023

  11. [11]

    G., Kleint, L., Leenaarts, J., Sukhorukov, A

    Judge, P. G., Kleint, L., Leenaarts, J., Sukhorukov, A. V ., & Vial, J.-C. 2020, ApJ, 901, 32

    work page 2020

  12. [12]

    K., & Peter, H

    Kayshap, P., Tripathi, D., Solanki, S. K., & Peter, H. 2018, ApJ, 864, 21

    work page 2018

  13. [13]

    & Carlsson, M

    Leenaarts, J. & Carlsson, M. 2009, 415, 87

    work page 2009

  14. [14]

    2009, ApJ, 694, L128

    Leenaarts, J., Carlsson, M., Hansteen, V ., & Rouppe van der V oort, L. 2009, ApJ, 694, L128

    work page 2009

  15. [15]

    2012, ApJ, 749, 136

    Leenaarts, J., Carlsson, M., & Rouppe van der V oort, L. 2012, ApJ, 749, 136

    work page 2012

  16. [16]

    Leenaarts, J., Pereira, T. M. D., Carlsson, M., Uitenbroek, H., & De Pontieu, B. 2013b, ApJ, 772, 90 Martínez-Sykora, J., de la Cruz Rodríguez, J., Goši ´c, M., et al. 2023, ApJ, 943, L14

    work page 2023

  17. [17]

    Milkey, R. W. & Mihalas, D. 1974, ApJ, 192, 769

    work page 1974

  18. [18]

    N., et al

    Ondratschek, P., Przybylski, D., Smitha, H. N., et al. 2024, A&A, 692, A6

    work page 2024

  19. [19]

    K., Rüedi, I., Landi, E., & Schühle, U

    Pauluhn, A., Solanki, S. K., Rüedi, I., Landi, E., & Schühle, U. 2000, A&A, 362, 737

    work page 2000

  20. [20]

    W., Labrosse, N., Schmieder, B., & Barczynski, K

    Peat, A. W., Labrosse, N., Schmieder, B., & Barczynski, K. 2021, A&A, 653, A5

    work page 2021

  21. [21]

    Pereira, T. M. D., Leenaarts, J., De Pontieu, B., Carlsson, M., & Uitenbroek, H. 2013, ApJ, 778, 143

    work page 2013

  22. [22]

    Pereira, T. M. D. & Uitenbroek, H. 2015, A&A, 574, A3

    work page 2015

  23. [23]

    K., et al

    Przybylski, D., Cameron, R., Solanki, S. K., et al. 2025, A&A, 703, A148

    work page 2025

  24. [24]

    K., et al

    Przybylski, D., Cameron, R., Solanki, S. K., et al. 2022, A&A, 664, A91

    work page 2022

  25. [25]

    Rybicki, G. B. & Hummer, D. G. 1991, A&A, 245, 171

    work page 1991

  26. [26]

    Rybicki, G. B. & Hummer, D. G. 1992, A&A, 262, 209

    work page 1992

  27. [27]

    H., Schou, J., Bush, R

    Scherrer, P. H., Schou, J., Bush, R. I., et al. 2012, Sol. Phys., 275, 207

    work page 2012

  28. [28]

    2015, ApJ, 811, 127

    Schmit, D., Bryans, P., De Pontieu, B., et al. 2015, ApJ, 811, 127

    work page 2015

  29. [29]

    H., Bush, R

    Schou, J., Scherrer, P. H., Bush, R. I., et al. 2012, Sol. Phys., 275, 229 Štˇepán, J., del Pino Alemán, T., & Trujillo Bueno, J. 2022, A&A, 659, A137

    work page 2012

  30. [30]

    Sukhorukov, A. V . & Leenaarts, J. 2017, A&A, 597, A46

    work page 2017

  31. [31]

    2001, ApJ, 557, 389

    Uitenbroek, H. 2001, ApJ, 557, 389

    work page 2001

  32. [32]

    E., Avrett, E

    Vernazza, J. E., Avrett, E. H., & Loeser, R. 1981, ApJS, 45, 635

    work page 1981

  33. [33]

    Withbroe, G. L. & Noyes, R. W. 1977, ARA&A, 15, 363 Article number, page 12 of 13 Ondratschek et al.: Mgiih&k spectra, 3D RT results 2796 2798 2800 2802 2804 Wavelength (Å) 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4Intensity ( n J m 2 s 1 Hz 1 sr 1) 1.5D RT 3D RT no vh 3D RT 40 20 02040 v (km s 1) 0.0 0.5 1.0 MgII-k Fig. A.1.The effect of horizontal velocities on th...

    work page 1977