pith. sign in

arxiv: 2509.17003 · v3 · submitted 2025-09-21 · 🌌 astro-ph.HE

On Atomic Line Opacities for Modeling Astrophysical Radiative Transfer

Jonathan Morag This is my paper

Pith reviewed 2026-05-18 14:56 UTC · model grok-4.3

classification 🌌 astro-ph.HE
keywords atomic line opacityradiative transferexpansion opacityEP93 formulaphoton emissivityfrequency averagingastrophysical simulationsequation of state
0
0 comments X

The pith

The EP93 expansion opacity underestimates photon emissivity and reprocessing rates even when it matches photon mean free paths.

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

This paper examines how atomic line opacities are approximated in astrophysical radiative transfer simulations through the line-expansion formalism. It focuses on the Eastman and Pinto 1993 formula that is widely used in codes such as STELLA and shows that the formula underestimates how photons are emitted and reprocessed. The underestimation explains large differences in predicted emissions between simulation groups. The authors introduce a modified frequency-bin averaging approach for emissivity that incorporates expansion effects on line strength and stress the need for electron excitation cutoffs in equations of state. Accurate treatment of these opacities matters because they control radiation transport in high-energy events where direct frequency resolution of every line is impossible.

Core claim

The widely used EP93 expansion opacity substantially underestimates photon emissivity and reprocessing rates, even when it correctly captures photon mean-free-paths. The paper reproduces EP93 opacities from STELLA simulations that previously showed orders-of-magnitude emission discrepancies with other codes. A new emissivity calculation is proposed by modifying the simple frequency-bin averaged opacity method to include the effect of expansion on effective line strength. The importance of micro-plasma electron excitation level cutoffs in the equation of state for opacity calculations is also highlighted. No fully consistent coarse-frequency solution for line modeling currently exists.

What carries the argument

The EP93 line-expansion formalism for frequency-averaged opacities, together with a proposed modification of frequency-bin averaging that accounts for expansion effects on effective line strength.

If this is right

  • Codes that rely on the EP93 formula will underpredict photon emission and reprocessing in environments with dense line forests.
  • Adopting the modified bin-averaged emissivity method should bring emission predictions from different simulation codes into closer agreement.
  • Opacity tables must incorporate electron excitation level cutoffs from the equation of state to avoid systematic errors in line strengths.
  • High-resolution frequency-dependent opacity tables can serve as a benchmark for validating coarse-frequency approximations.

Where Pith is reading between the lines

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

  • Revised emissivity treatments could alter predicted light curves and spectra for explosive transients such as supernovae.
  • Other radiative transfer codes using similar expansion opacity approximations may share the same underestimation bias.
  • Full radiation-hydrodynamics runs with the new method would test whether improved emissivity changes overall energy transport and observed signatures.

Load-bearing premise

The orders-of-magnitude emission discrepancies between STELLA and other simulations arise mainly from differences in line opacity treatment rather than from differences in hydrodynamics or equation-of-state implementations.

What would settle it

A side-by-side computation of emissivity and reprocessing rates using the EP93 formula versus the proposed modified bin-averaged method in an isolated radiative transfer test, compared against a fully resolved high-frequency reference calculation.

Figures

Figures reproduced from arXiv: 2509.17003 by Jonathan Morag.

Figure 1
Figure 1. Figure 1: Our imitation of the bound-free opacity in B98 (black line) for the example case of = 10−13 cm−3 , = 15, 000 K for a solar mixture. In red dashed-line, we show the result when the Hummer & Mihalas (1988) factor is not included, and we limit the Hydrogen partition function to max = 400, finding reasonable agreement with B98. In blue solid lines we show the result of a converged H partition function, represe… view at source ↗
Figure 2
Figure 2. Figure 2: Bound-bound opacity example from B98 fig. 1, compared with our reproduction using M23 modified to employ the EP93 prescription. We finding good agreement to a factor of a few or better. The plasma parameters are = 10−13 cm−3 , = 15, 000 K, exp = 15 days. Similarly to B98, we use a coarse frequency grid with Δ/ ∼ 0.01. The lines that are in excellent agreement in the range 200 Å < < 400 Å are dominated by a… view at source ↗
read the original abstract

In astrophysics, atomic transition line opacity is a primary source of uncertainty in theoretical calculations of radiative transfer. Much of this uncertainty is dominated by the inability to resolve the lines in frequency, leading to the use of approximate frequency-averaged treatments, often employing the `line-expansion formalism'. In this short paper we assess the usage of this formalism in simulations, specifically the prominent Eastman \& Pinto 1993 formula (hereafter EP93). As a case study, we reproduce EP93 opacities from the commonly-used STELLA simulations. The latter previously yielded orders of magnitude discrepancy in observed emission relative to similar simulations from our group. The discrepancy is due to differences in line opacity treatment. We show that the widely used EP93 expansion opacity substantially underestimates photon emissivity and reprocessing rates, even when it correctly captures photon mean-free-paths. We also highlight the importance of introducing micro-plasma electron excitation level cutoffs in the equation of state (EOS) for calculating opacity. We propose a new method for calculating emissivity, based on a modification of the simple frequency-bin averaged opacity method, in a way that incorporates the effect of expansion on effective line strength. This formulation should reduce the overestimation of the opacity that may occur with the simple averaging method. To our knowledge, no fully-consistent coarse-frequency solution currently exists for line modeling in these systems. Finally, we describe new features in our updated publicly available high-resolution frequency-dependent opacity table.

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 / 2 minor

Summary. The manuscript evaluates the widely used Eastman & Pinto 1993 (EP93) expansion-opacity formalism for atomic line opacities in astrophysical radiative transfer, using supernova modeling as a case study. The authors reproduce EP93 opacities inside the STELLA code and compare the results to their own simulations, which previously showed orders-of-magnitude differences in observed emission. They conclude that EP93 substantially underestimates photon emissivity and reprocessing rates even while correctly capturing photon mean free paths. A new emissivity calculation method is proposed that modifies the simple frequency-bin averaged opacity approach to incorporate expansion effects on effective line strength. The paper also stresses the need for micro-plasma electron excitation level cutoffs in the EOS when computing opacities and describes new features in the authors' publicly available high-resolution frequency-dependent opacity tables.

Significance. If the central claims are substantiated, the work identifies a potentially consequential limitation in a standard opacity treatment that could affect radiative-transfer predictions in expanding media. The proposed emissivity modification and the release of updated public opacity tables are constructive steps toward more consistent coarse-frequency solutions. These elements, together with the emphasis on EOS cutoffs, address a recognized source of uncertainty in the field.

major comments (2)
  1. Abstract: The statement that 'the discrepancy is due to differences in line opacity treatment' and that EP93 'substantially underestimates photon emissivity and reprocessing rates' rests on a comparison between STELLA (with EP93) and the authors' code. No controlled experiment is described in which hydrodynamics, EOS implementation, grid, or time-stepping are held fixed while only the opacity formalism is swapped. Without such isolation the orders-of-magnitude emission difference cannot be unambiguously assigned to EP93.
  2. Abstract and the section describing the STELLA reproduction: The claim that EP93 'correctly captures photon mean-free-paths' while underestimating emissivity requires quantitative verification (e.g., explicit mean-free-path and emissivity ratios or spectra) that is not visible in the provided abstract; the strength of the evidence therefore depends on details that must be shown explicitly.
minor comments (2)
  1. The abstract refers to 'new features in our updated publicly available high-resolution frequency-dependent opacity table' but does not indicate the version number, repository location, or which specific improvements (e.g., frequency resolution, cutoff handling) are new.
  2. Clarify the precise mathematical modification to the frequency-bin averaged opacity that incorporates the expansion effect on effective line strength; an explicit formula or algorithm would aid reproducibility.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful and constructive review of our manuscript. We address each major comment below in a point-by-point manner and have revised the manuscript where appropriate to strengthen the presentation and evidence.

read point-by-point responses

  1. Referee: Abstract: The statement that 'the discrepancy is due to differences in line opacity treatment' and that EP93 'substantially underestimates photon emissivity and reprocessing rates' rests on a comparison between STELLA (with EP93) and the authors' code. No controlled experiment is described in which hydrodynamics, EOS implementation, grid, or time-stepping are held fixed while only the opacity formalism is swapped. Without such isolation the orders-of-magnitude emission difference cannot be unambiguously assigned to EP93.

    Authors: We acknowledge that a fully controlled experiment isolating only the opacity formalism within a single code would provide the strongest possible attribution. Our approach instead reproduces the EP93 opacities directly inside the STELLA code under the same conditions as the original STELLA runs that exhibited the emission discrepancy with our group's simulations. This reproduction allows direct comparison of the resulting opacities and derived emissivities. We have revised the abstract and the STELLA reproduction section to more precisely describe the scope of the comparison, to note the challenges of cross-code isolation, and to emphasize that the emissivity differences arise under matched opacity implementations. We believe this addresses the concern without overstating isolation. revision: partial

  2. Referee: Abstract and the section describing the STELLA reproduction: The claim that EP93 'correctly captures photon mean-free-paths' while underestimating emissivity requires quantitative verification (e.g., explicit mean-free-path and emissivity ratios or spectra) that is not visible in the provided abstract; the strength of the evidence therefore depends on details that must be shown explicitly.

    Authors: We agree that explicit quantitative support is necessary for the claim. The full manuscript already contains these comparisons in the section on the STELLA reproduction, including figures that display mean-free-path ratios (typically agreeing within a factor of ~2) alongside emissivity and reprocessing rate differences exceeding an order of magnitude, as well as sample spectra. To improve visibility, we have revised the abstract to reference these quantitative results directly and have added explicit statements of the ratios in the main text with clearer figure citations. revision: yes

Circularity Check

0 steps flagged

No circularity: derivation compares external STELLA implementation to group code without reducing claims to self-definition or fitted inputs

full rationale

The paper reproduces the external EP93 formalism inside the publicly documented STELLA code and attributes emission discrepancies to opacity differences by direct comparison with that external implementation. No equations are shown to be equivalent by construction to their own inputs, no parameters are fitted to a subset and then relabeled as predictions, and no load-bearing uniqueness theorem or ansatz is imported solely via self-citation. The proposed emissivity modification is presented as a new averaging adjustment rather than a renaming of a prior result. The manuscript therefore remains self-contained against the external STELLA benchmark.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claims rest on the assumption that line opacity differences dominate the simulation discrepancy and that micro-plasma excitation cutoffs are necessary for accurate EOS-based opacity. No free parameters or invented entities are introduced in the abstract; the work modifies an existing averaging method rather than postulating new physics.

axioms (2)
  • domain assumption Line opacity treatment is the dominant source of the orders-of-magnitude emission discrepancy between STELLA and the authors' simulations.
    Stated in the abstract as the explanation for the observed difference.
  • domain assumption Micro-plasma electron excitation level cutoffs must be included in the EOS to calculate opacity correctly.
    Highlighted as important for opacity calculations.

pith-pipeline@v0.9.0 · 5784 in / 1424 out tokens · 32178 ms · 2026-05-18T14:56:23.258443+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

32 extracted references · 32 canonical work pages

  1. [1]

    L., Lund K

    Barnes J., Zhu Y. L., Lund K. A., Sprouse T. M., Vassh N., McLaughlin G. C., Mumpower M. R., Surman R., 2021, @doi [The Astrophysical Journal] 10.3847/1538-4357/ac0aec , 918, 44

    work page doi:10.3847/1538-4357/ac0aec 2021

  2. [2]

    H., Nugent P., Branch D., 1996, @doi [Monthly Notices of the Royal Astronomical Society] 10.1093/mnras/283.1.297 , 283, 297

    Baron E., Hauschildt P. H., Nugent P., Branch D., 1996, @doi [Monthly Notices of the Royal Astronomical Society] 10.1093/mnras/283.1.297 , 283, 297

    work page doi:10.1093/mnras/283.1.297 1996

  3. [3]

    Ben Nasr S., Carvajal Gallego H., Deprince J., Palmeri P., Quinet P., 2023, @doi [Astronomy & Astrophysics] 10.1051/0004-6361/202346198 , 678, A67

    work page doi:10.1051/0004-6361/202346198 2023

  4. [4]

    I., 1997, in Ruiz-Lapuente P., Canal R., Isern J., eds, NATO ASI Series , Thermonuclear Supernovae

    Blinnikov S. I., 1997, in Ruiz-Lapuente P., Canal R., Isern J., eds, NATO ASI Series , Thermonuclear Supernovae . Springer Netherlands, Dordrecht, pp 589--605, @doi 10.1007/978-94-011-5710-0_36 , https://doi.org/10.1007/978-94-011-5710-0_36

    work page doi:10.1007/978-94-011-5710-0_36 1997

  5. [5]

    I., Eastman R., Bartunov O

    Blinnikov S. I., Eastman R., Bartunov O. S., Popolitov V. A., Woosley S. E., 1998, @doi [The Astrophysical Journal] 10.1086/305375 , 496, 454

    work page doi:10.1086/305375 1998

  6. [6]

    Bulla M., 2023, @doi [Monthly Notices of the Royal Astronomical Society] 10.1093/mnras/stad232 , 520, 2558

    work page doi:10.1093/mnras/stad232 2023

  7. [7]

    I., 2007, Radiation Hydrodynamics

    Castor J. I., 2007, Radiation Hydrodynamics . Cambridge, UK ; New York

    work page 2007

  8. [8]

    Colgan J., et al., 2016, @doi [The Astrophysical Journal] 10.3847/0004-637X/817/2/116 , 817, 116

    work page doi:10.3847/0004-637x/817/2/116 2016

  9. [9]

    J., 2015, @doi [Monthly Notices of the Royal Astronomical Society] 10.1093/mnras/stv609 , 449, 4304

    Dessart L., Audit E., Hillier D. J., 2015, @doi [Monthly Notices of the Royal Astronomical Society] 10.1093/mnras/stv609 , 449, 4304

    work page doi:10.1093/mnras/stv609 2015

  10. [10]

    Domoto N., Tanaka M., Kato D., Kawaguchi K., Hotokezaka K., Wanajo S., 2022, @doi [The Astrophysical Journal] 10.3847/1538-4357/ac8c36 , 939, 8

    work page doi:10.3847/1538-4357/ac8c36 2022

  11. [11]

    G., Pinto P

    Eastman R. G., Pinto P. A., 1993, @doi [The Astrophysical Journal] 10.1086/172957 , 412, 731

    work page doi:10.1086/172957 1993

  12. [12]

    B., Castor J

    Friend D. B., Castor J. I., 1983, @doi [The Astrophysical Journal] 10.1086/161289 , 272, 259

    work page doi:10.1086/161289 1983

  13. [13]

    Förster F., et al., 2018, @doi [Nature Astronomy] 10.1038/s41550-018-0563-4 , 2, 808

    work page doi:10.1038/s41550-018-0563-4 2018

  14. [14]

    C., Pain J.-C., Godefroid M., Palmeri P., Quinet P., 2024, @doi [Journal of Physics B: Atomic, Molecular and Optical Physics] 10.1088/1361-6455/ad2182 , 57, 035001

    Gallego H. C., Pain J.-C., Godefroid M., Palmeri P., Quinet P., 2024, @doi [Journal of Physics B: Atomic, Molecular and Optical Physics] 10.1088/1361-6455/ad2182 , 57, 035001

    work page doi:10.1088/1361-6455/ad2182 2024

  15. [15]

    2012, MNRAS, 423, 2533, doi:10.1111/j.1365-2966.2012.21057.x17

    Hillier D. J., Dessart L., 2012, @doi [Monthly Notices of the Royal Astronomical Society] 10.1111/j.1365-2966.2012.21192.x , 424, 252

    work page doi:10.1111/j.1365-2966.2012.21192.x 2012

  16. [16]

    G., Mihalas D., 1988, @doi [The Astrophysical Journal] 10.1086/166600 , 331, 794

    Hummer D. G., Mihalas D., 1988, @doi [The Astrophysical Journal] 10.1086/166600 , 331, 794

    work page doi:10.1086/166600 1988

  17. [17]

    H., Lasher G., Chan K

    Karp A. H., Lasher G., Chan K. L., Salpeter E. E., 1977, @doi [The Astrophysical Journal] 10.1086/155241 , 214, 161

    work page doi:10.1086/155241 1977

  18. [18]

    C., & Nugent, P

    Kasen D., Thomas R. C., Nugent P., 2006, @doi [The Astrophysical Journal] 10.1086/506190 , 651, 366

    work page doi:10.1086/506190 2006

  19. [19]

    Kawaguchi K., Shibata M., Tanaka M., 2020, @doi [The Astrophysical Journal] 10.3847/1538-4357/ab61f6 , 889, 171

    work page doi:10.3847/1538-4357/ab61f6 2020

  20. [20]

    Kozyreva A., Nakar E., Waldman R., Blinnikov S., Baklanov P., 2020a, @doi [Monthly Notices of the Royal Astronomical Society] 10.1093/mnras/staa924 , 494, 3927

    work page doi:10.1093/mnras/staa924

  21. [21]

    Kozyreva A., Shingles L., Mironov A., Baklanov P., Blinnikov S., 2020b, @doi [Monthly Notices of the Royal Astronomical Society] 10.1093/mnras/staa2704 , 499, 4312

    work page doi:10.1093/mnras/staa2704

  22. [22]

    B., & Bernardi, M

    Kromer M., Sim S. A., 2009, @doi [Monthly Notices of the Royal Astronomical Society] 10.1111/j.1365-2966.2009.15256.x , 398, 1809

    work page doi:10.1111/j.1365-2966.2009.15256.x 2009

  23. [23]

    A., 2021, Astrophysics Source Code Library, p

    Kromer M., Sim S. A., 2021, Astrophysics Source Code Library, p. ascl:2103.020

    work page 2021

  24. [24]

    L., 1995, in ASP Conference Series

    Kurucz R. L., 1995, in ASP Conference Series . p. 583, http://adsabs.harvard.edu/abs/1995ASPC...81..583K

    work page 1995

  25. [25]

    B., 2002, @doi [ ] 10.1051/0004-6361:20011756 , https://ui.adsabs.harvard.edu/abs/2002A&A...384..725L 384, 725

    Lucy L. B., 2002, @doi [Astronomy and Astrophysics] 10.1051/0004-6361:20011756 , 384, 725

    work page doi:10.1051/0004-6361:20011756 2002

  26. [26]

    B., 2003, @doi [ ] 10.1051/0004-6361:20030357 , https://ui.adsabs.harvard.edu/abs/2003A&A...403..261L 403, 261

    Lucy L. B., 2003, @doi [Astronomy and Astrophysics] 10.1051/0004-6361:20030357 , 403, 261

    work page doi:10.1051/0004-6361:20030357 2003

  27. [27]

    W., 1999, Foundations of Radiation Hydrodynamics

    Mihalas D., Mihalas B. W., 1999, Foundations of Radiation Hydrodynamics . Mineola, NY

    work page 1999

  28. [28]

    Morag J., 2023, Frequency Dependent Opacity Table , https://github.com/jon-morag/Freq_Dept_Opac_Table/

    work page 2023

  29. [29]

    Morag J., Irani I., Sapir N., Waxman E., 2024, @doi [Monthly Notices of the Royal Astronomical Society] 10.1093/mnras/stae374 , 528, 7137

    work page doi:10.1093/mnras/stae374 2024

  30. [30]

    A., Eastman R

    Pinto P. A., Eastman R. G., 2000, @doi [The Astrophysical Journal] 10.1086/308380 , 530, 757

    work page doi:10.1086/308380 2000

  31. [31]

    S., Blinnikov S

    Potashov M. S., Blinnikov S. I., Sorokina E. I., 2021, @doi [Astronomy Letters] 10.1134/S1063773721030051 , 47, 204

    work page doi:10.1134/s1063773721030051 2021

  32. [32]

    I., Baklanov P., Sorokina E

    Tominaga N., Morokuma T., Blinnikov S. I., Baklanov P., Sorokina E. I., Nomoto K., 2011, @doi [The Astrophysical Journal Supplement Series] 10.1088/0067-0049/193/1/20 , 193, 20

    work page doi:10.1088/0067-0049/193/1/20 2011