REVIEW 3 major objections 4 minor 1 cited by
Kilonova ejecta opacity inferred from new large-scale HFR atomic calculations in all elements between Ca (Z = 20) and Lr (Z = 103)
T0 review · 3 major / 4 minor · reviewed 2026-08-11 · deepseek-v4-flash
Pith's one-line read Lanthanides do not dominate kilonova opacity, new atomic data show.
desk verdict A solid, citable full-periodic-table opacity resource whose main claims hold up; the lanthanide conclusion depends on one abundance model, and the composition-mixing procedure deserves a proper multi-element Saha test. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The machinery is the composition-weighted expansion opacity, built on the standard Sobolev expansion formula $\kappa_{\mathrm{exp}}(\lambda) = (1/ct\rho)\sum_l (\lambda_l/\Delta\lambda)(1-e^{-\tau_l})$ and its Planck mean. The defining move is to scale each element's lower-level number density inside the optical depth by its molar abundance, $n_{l,Z} = y_Z n_l$, so the mixture opacity is computed before summing over elements rather than by averaging completed single-element opacities. The atomic input comes from large multiconfiguration pseudo-relativistic Hartree-Fock calculations for ions I-IV of all elements from calcium to lawrencium; the paper argues that the level errors and inverted ground states these models leave (18-30% for most ions, 49% for Nd III) have little effect on expansion opacities, based on tests for Nd, U, and Er. The line-binned opacity $\kappa^{\mathrm{bin}}_\nu\propto \sum_l N_l f_l$ scales linearly with abundance and therefore does not suffer from the a posteriori averaging problem.
What would settle it
Recompute the full-ejecta Planck opacity at $t=3.5$ d, $T=6000$ K, $\rho=10^{-13}\ \mathrm{g\,cm^{-3}}$ with the paper's HFR tables but with r-process abundances from a different nuclear mass model or fission prescription; if the lanthanide and actinide contribution then exceeds the d-shell contribution, the average non-dominance claim fails. A second check is to compare the predicted $\sim$7-8 day bolometric peak with the observed AT2017gfo light curve; if the HFR-opacity model is conclusively excluded by the timing, the opacity shift is not the right description.
Extended reading notes
Core claim
The discovery is a change in who controls kilonova opacity when the ejecta composition is taken seriously. For the sym-n1-a6 merger model at 3.5 days ($T=6000\ \mathrm{K}$, $\rho=10^{-13}\ \mathrm{g\,cm^{-3}}$), the composition-weighted HFR Planck opacity is $1.13\ \mathrm{cm^2\,g^{-1}}$, about 25% below the $1.43\ \mathrm{cm^2\,g^{-1}}$ obtained from a parametric formula that links opacity to the lanthanide-plus-actinide molar fraction. The opacity budget is shared among lanthanides and actinides and the much more abundant elements near $Z\simeq 24$ and $Z\simeq 40$. The paper's second result is methodological: because the expansion opacity is a nonlinear Sobolev sum over lower-level populations, replacing those populations by abundance-scaled ones inside the formula is not the same as averaging single-element opacities; the a posteriori average gives $0.11\ \mathrm{cm^2\,g^{-1}}$, an order of magnitude too low. When the HFR opacities replace the parametric formula in light-curve simulations, the bolometric peak moves from roughly 1-2 days to about 7-8 days, with a slightly lower peak luminosity for the full ejecta.
Load-bearing premise
The load-bearing premise is that the single abundance pattern from the sym-n1-a6 merger simulation with the BSkG3 mass model is representative enough that the cancellation between low lanthanide abundance and high lanthanide opacity holds in real kilonova ejecta, because other r-process models vary that abundance by orders of magnitude.
Editorial extensions
If this is right
- Kilonova light curves computed with these opacities peak at roughly 7-8 days rather than 1-2 days, a shift that is directly testable against the observed evolution of AT2017gfo.
- Lanthanide and actinide molar fraction alone is not a reliable opacity proxy; the parametric formula built on it is 25% high for the full ejecta and a factor of two low for the dynamical ejecta.
- Single-element expansion opacity tables should not be combined with arbitrary abundance patterns a posteriori; the correct composition-weighted calculation changes the answer by an order of magnitude here.
- The published grid of HFR atomic data and expansion/line-binned opacity tables for $Z=20$-103, ions I-IV, over $t=1$-7 days, $\rho=10^{-17}$-$10^{-13}\ \mathrm{g\,cm^{-3}}$, and $T=1000$-10000 K provides a common baseline for kilonova modeling.
- When only the fast dynamical ejecta is considered, the HFR opacities give a brighter peak than the parametric opacity, so atomic data matter even in the lanthanide-dominated component.
Reading between the lines
- The paper's lanthanide-non-dominance conclusion should be read as abundance-pattern dependent: replacing the BSkG3-based abundances with another r-process mass model or fission prescription, which can change lanthanide fractions by orders of magnitude, could restore lanthanide dominance, and this test is not performed in the paper.
- The same composition-weighting failure likely affects any kilonova or transient model that stores Planck mean opacities per element and later mixes them with abundances; line-binned opacities, being linear in abundance, are the safer stored quantity for such post-processing.
- Recomputing the light curve with per-trajectory compositions instead of the uniform composition assumed here could sharpen or soften the 7-8 day peak shift; the paper identifies this as future work.
- If the non-dominance holds broadly, the atomic-data bottleneck shifts from all lanthanides to a short list of species (Cr, Fe, Zr, Mo, Sm, Nd, Dy, Er, U, Np, Pu), and improving those atoms experimentally and theoretically would be the highest-value next step.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper presents new pseudo-relativistic Hartree-Fock (HFR) atomic structure calculations for all elements from Ca (Z = 20) to Lr (Z = 103) in ionization stages I-IV, together with expansion and line-binned opacities for a kilonova-relevant grid (T = 1000-10000 K, rho = 1e-17 to 1e-13 g/cm^3, t = 1-7 days). The data are made publicly available on Zenodo. The authors compare their Planck mean opacities with Tanaka et al. (2020) and Fontes et al. (2020, 2023), and apply the opacities to the neutron star merger model sym-n1-a6 from Just et al. (2023). Their central result is that, for this model, the composition-weighted HFR expansion opacity is 1.13 cm^2/g, about 25% lower than the parametric prescription of Just et al. (2022), and that lanthanides are not the dominant contributors 'at least on average'; instead, 3d-shell (Z ~ 24) and 4d-shell (Z ~ 40) elements contribute substantially. They also find that abundance-weighted single-element Planck opacities (0.11 cm^2/g) underestimate the mixture opacity by an order of magnitude, and that HFR opacities shift the modeled bolometric peak from roughly 1-2 days to roughly 7-8 days.
Significance. The database is a major community resource: consistent HFR oscillator strengths and energy levels for 84 elements in four charge states are not currently available elsewhere with this coverage. The comparison with existing data sets and the public release are valuable. If the central claim survives scrutiny, it challenges the common practice of treating lanthanide opacity as the sole determinant of kilonova red/infrared emission and the practice of building mixture opacities by abundance-weighting single-element tables. The result is, however, conditional on the adopted nucleosynthesis composition and on the composition-mixing approximation discussed below. The manuscript is generally careful in hedging ('for a given model', 'at least on average'), but the headline result deserves stronger validation.
major comments (3)
- [§6, Eq. (9)] The composition-weighted opacity is constructed by scaling pure-element lower-level densities n_l by the elemental molar fraction y_Z, i.e., n_{l,Z} = y_Z n_l. This does not solve the Saha equation (Eq. 5) with a common electron density n_e and charge conservation for the mixture; it implicitly assumes that each element's ionization balance in the mixture is identical to that in a pure plasma of that element at the same bulk density. Because n_e in a mixture is dominated by the abundant lighter elements, while the trace lanthanides (molar fraction 8.4e-5 in the full ejecta) have their own pure-plasma n_e, the ion fractions of the important species are not established. The claimed totals (1.13 cm^2/g for the full ejecta, 5.42 cm^2/g for the dynamical component) and the factor-10 discrepancy with the abundance-weighted single-element value (0.11 cm^2/g) therefore rest on an untested approximation. Please either implement a self-consistent multi-element LTE calculation for the sym-n1-a6 composition, or provide a quantitative demonstration that the pure-element Saha fractions are insensitive to n_e over the relevant range.
- [§2] The paper states that HFR energy-level discrepancies of 18-30% (49% for Nd III) and inverted ground states in 23 of the 120 lanthanide/actinide ions leave expansion opacities essentially unaffected, but this claim is demonstrated only for Nd II/III, U II/III, and Er III. The extrapolation to all 429 ions is a load-bearing step because the main astrophysical conclusion depends on the relative opacities of d-shell, lanthanide, and actinide species. A sensitivity test covering representative ions from each group (e.g., Cr, Mo, Sm, U) with and without level-energy calibration, or a propagated uncertainty estimate on the total mixture opacity, would substantially strengthen the robustness of the result.
- [§6] The 'lanthanides are not dominant' conclusion is drawn from a single r-process abundance pattern, model sym-n1-a6 with the BSkG3 mass model. The lanthanide-plus-actinide molar fraction in this model is 8.4e-5 in the full ejecta and 8.7e-4 in the dynamical component. Because single-element lanthanide/actinide opacities exceed those of the abundant d-shell elements by several orders of magnitude (Section 3), the ranking is a cancellation between per-element opacity and abundance. Published r-process calculations with other mass models or fission treatments show lanthanide fractions varying by orders of magnitude. The paper should either test the sensitivity of the ranking and of the 1.13 cm^2/g total to variations in the lanthanide/actinide fraction within the range of current models, or explicitly narrow the claim to this specific nucleosynthesis model in the abstract.
minor comments (4)
- [Figure captions 2-5] The density unit in the figure captions is given as '10^-13 cm^-1'; it should be '10^-13 g cm^-3' as in the text.
- [Section 7] Typos: 'The purpose of this work is to is to build' and 'impoortance' should be corrected; 'Eriii' should be 'Er III'; in the reference list, 'Astronomy & Astroohysics' should be 'Astronomy & Astrophysics'.
- [Acknowledgements] The reference to 'KILONOV A' should be 'KILONOVA'.
- [Eq. (7)] The Planck mean opacity definition uses B_lambda(T) in both numerator and denominator; it may be worth noting explicitly that the integration is over the full wavelength range, as stated in the text, to avoid confusion with the line-binned version.
Circularity Check
No significant circularity: the opacities are forward-computed from HFR atomic structure calculations and are benchmarked against independent methods and experimental levels.
full rationale
The paper's central outputs—per-element expansion opacities, Planck-mean opacities, and the composition-weighted ejecta opacity—are obtained from ab initio atomic structure calculations (HFR) combined with an external r-process abundance pattern, with no parameter fitted to opacity or light-curve targets. The comparison between the Eq. (9) mixture treatment and abundance-averaged single-element opacities is a computed consequence of the nonlinear expansion-opacity formula, not an input assumed into the result. The self-citations to Flörs et al. (2023) and Deprince et al. (2023, 2024) serve as convergence and accuracy studies that include independent benchmarks (FAC and experimental energy levels), so they do not reduce the argument to an unverified self-citation. The lanthanide non-dominance conclusion depends on one nucleosynthesis model, which is a robustness limitation rather than a circularity. No load-bearing step in the derivation is equivalent by construction to its own output.
Assumptions & free parameters
free parameters (3)
- Slater integral scaling factor =
0.85
- log gf cut-off =
-5
- Configuration list restrictions =
per-ion lists in Table 1 (Appendix)
assumptions (6)
- domain assumption LTE holds in the KN photosphere from 1 day to about 1 week post-merger
- domain assumption Only ionization stages I-IV are relevant at T <= 10000 K
- standard math Sobolev expansion opacity formalism is valid for NS merger ejecta
- domain assumption The sym-n1-a6 abundance pattern (BSkG3 mass model) is representative for the opacity-dominance question
- domain assumption Expansion opacities are insensitive to the documented HFR level errors and ground-state inversions
- domain assumption CI models built by homologous sequence from studied ions are converged for every element
Cite this review
Pith. "Pith review of Kilonova ejecta opacity inferred from new large-scale HFR atomic calculations in all elements between Ca (Z = 20) and Lr (Z = 103)." pith.science (2026). https://pith.science/paper/5X3CIEKI
@misc{pith2026241216688,
author = {Pith},
title = {Pith review of: Kilonova ejecta opacity inferred from new large-scale HFR atomic calculations in all elements between Ca (Z = 20) and Lr (Z = 103)},
year = {2026},
howpublished = {\url{https://pith.science/paper/5X3CIEKI}},
note = {Machine review of arXiv:2412.16688}
}
read the original abstract
In the context of kilonova (KN) modeling, the present work focusses on large-scale atomic data and opacity computations for all heavy elements from Ca to Lr, with a special effort on lanthanides and actinides, for a grid of typical KN ejecta conditions between one day and one week after the merger (corresponding to the LTE photosphere phase of the KN ejecta). In order to do so, we used the pseudo-relativistic Hartree-Fock (HFR) method, in which the choice of the interaction configuration model is of crucial importance. In this paper, HFR atomic data and opacities for all elements between Ca (Z = 20) and Lr (Z = 103) are presented, with a special focus on lanthanides and actinides. Besides, we also discuss the contribution of every single element to the total KN ejecta opacity for a given neutron star merger model, depending on their Planck mean opacities and elemental abundances. An important result is that lanthanides are found to not be the dominant sources of opacity, at least on average. The impact on KN light curves of considering such atomic-physics based opacity data instead of typical crude approximation formulae is also evaluated. In addition, the importance of taking the ejecta composition into account directly in the expansion opacity determination (instead of estimating single-element opacities) is highlighted. A database containing all the relevant atomic data and opacity tables has also been created and published online along with this work.
Figures
Figures from the paper (8 more)
Forward citations
Cited by 1 Pith paper
-
Calibrated Lanthanide Atomic Data for Kilonova Radiative Transfer. I. Atomic Structure and Opacities
A new calibrated atomic dataset for all singly and doubly ionized lanthanides provides 66,591 experimentally anchored transition wavelengths for kilonova modeling.
Reference graph
Works this paper leans on
- [1]
- [2]
- [3]
-
[4]
Banerjee, S., Tanaka, M., Kato, D., & Gaigalas, G. 2024, ApJ, 968, 64
work page 2024
-
[5]
Bar-Shalom, A., Klapisch, M., & Oreg, J. 2001, Journal of Quantitative Spec- troscopy and Radiative Transfer, 71, 169 Ben Nasr, S., Carvajal Gallego, H., Deprince, J., Palmeri, P., & Quinet, P. 2023, Astronomy & Astroohysics, 678 Ben Nasr, S., Carvajal Gallego, H., Deprince, J., Palmeri, P., & Quinet, P. 2024, Astronomy & Astrophysics, 687
work page 2001
-
[6]
Blaise, J. & Wyart, J. F. 1992, Energy levels and atomic spectra of actinides. (Tables internationales de constantes), pages: 479
work page 1992
-
[7]
Blaise, J. & Wyart, J.-F. 1994 Carvajal Gallego, H., Berengut, J., Palmeri, P., & Quinet, P. 2022a, MNRAS, 509 Carvajal Gallego, H., Berengut, J., Palmeri, P., & Quinet, P. 2022b, MNRAS, 513 Carvajal Gallego, H., Deprince, J., Berengut, J., Palmeri, P., & Quinet, P. 2023a, MNRAS, 518 Carvajal Gallego, H., Deprince, J., Godefroid, M., & et al. 2023b, Eur. ...
work page 1994
- [8]
Show all 55 references
-
[9]
Cowan, R. D. 1981, The theory of atomic structure and spectra (University of California Press)
1981
-
[10]
2024, Eur
Deprince, J., Carvajal Gallego, H., Ben Nasr, S., et al. 2024, Eur. Phys. J. D, 78
2024
-
[11]
2023, Eur
Deprince, J., Carvajal Gallego, H., Godefroid, M., et al. 2023, Eur. Phys. J. D, 77, 93
2023
-
[12]
N., Kononov, E
Ding, M., Ryabtsev, A. N., Kononov, E. Y ., Ryabchikova, T., & Pickering, J. C. 2024, Astronomy & Astrophysics, 78
2024
-
[13]
Eastman, R. G. & Pinto, P. A. 1993, ApJ, 412, 731
1993
-
[14]
1989, Nature, 340, 126 Flörs, A., Silva, R
Eichler, D., Livio, M., Piran, T., & Schramm, D. 1989, Nature, 340, 126 Flörs, A., Silva, R. F., Deprince, J., et al. 2023, Monthly Notices of the Royal Astronomical Society, 524, 3083
1989
-
[15]
J., Fryer, C
Fontes, C. J., Fryer, C. L., Hungerford, A. L., Wollaeger, R. T., & Korobkin, O. 2020, MNRAS, 493, 4143
2020
-
[16]
J., Fryer, C
Fontes, C. J., Fryer, C. L., Wollaeger, R. T., Mumpower, M. R., & Sprouse, T. M. 2023, MNRAS, 519, 2862
2023
-
[17]
2019, Astrophys
Gaigalas, G., Kato, D., Rynkun, P., Radži¯ut˙e, L., & Tanaka, M. 2019, Astrophys. J. Suppl., 240, 29
2019
-
[18]
2022, MNRAS, 517, 281
Gaigalas, G., Rynkun, P., Banerjee, S., et al. 2022, MNRAS, 517, 281
2022
-
[19]
2020, Astrophys
Gaigalas, G., Rynkun, P., Radži¯ut˙e, L., & et al. 2020, Astrophys. J. Suppl., 248, 13
2020
-
[20]
H., Smartt, S
Gillanders, J. H., Smartt, S. J., Sim, S. A., Bauswein, A., & Goriely, S. 2022, Mon. Not. Roy. Astron. Soc., 515, 631
2022
-
[21]
2011, Astrophys
Goriely, S., Bauswein, A., & Janka, H.-T. 2011, Astrophys. J. Lett., 738, L32
2011
-
[22]
2010, Phys
Goriely, S., Chamel, N., & Pearson, J. 2010, Phys. Rev. C, 82, 035804
2010
-
[23]
2023, Eur
Grams, G., Ryssens, W., Scamps, G., Goriely, S., & Chamel, N. 2023, Eur. Phys. J. A, 59, 270
2023
-
[24]
Gu, M. F. 2008, Canadian Journal of Physics, 86, 675
2008
-
[25]
2017, Handbook of Supernovae, ed
Janka, H.-T. 2017, Handbook of Supernovae, ed. A. W. Alsabti & P. Murdin, 1095
2017
-
[26]
2015, MNRAS, 448, 541
Just, O., Bauswein, A., Ardevol Pulpillo, R., Goriely, S., & Janka, H.-T. 2015, MNRAS, 448, 541
2015
-
[27]
2022, MNRAS, 509, 1377
Just, O., Goriely, S., Janka, H.-T., Nagataki, S., & Bauswein, A. 2022, MNRAS, 509, 1377
2022
-
[28]
2022, MNRAS, 510, 2820
Just, O., Kullmann, I., Goriely, S., et al. 2022, MNRAS, 510, 2820
2022
-
[29]
2023, The Astrophysical Journal Letters, 951, L12
Just, O., Vijayan, V ., Xiong, Z., et al. 2023, The Astrophysical Journal Letters, 951, L12
2023
-
[30]
R., & Barnes, J
Kasen, D., Badnell, N. R., & Barnes, J. 2013, ApJ, 774, 25
2013
-
[31]
2017, Nature, 551, 80
Kasen, D., Metzger, B., Barnes, J., Quataert1, E., & Ramirez-Ruiz, E. 2017, Nature, 551, 80
2017
-
[32]
2023, Eur
Koning, A., Hilaire, S., & Goriely, S. 2023, Eur. Phys. J. A, 59, 131
2023
-
[33]
2012, MNRAS, 26, 1940
Korobkin, O., Rosswog, S., Arcones, A., & Winteler, C. 2012, MNRAS, 26, 1940
2012
-
[34]
2023, MNRAS, 523, 2551
Kullmann, I., Goriely, S., Just, O., Bauswein, A., & Janka, H.-T. 2023, MNRAS, 523, 2551
2023
-
[35]
& Schramm, D
Lattimer, J. & Schramm, D. 1974, ApJ, 192, L145 Lemaître, J.-F., Goriely, S., Bauswein, A., & Janka, H.-T. 2021, Phys. Rev. C, 103, 025806
1974
-
[36]
2022, Atoms, 10, 130
Maison, L., Carvajal Gallego, H., & Quinet, P. 2022, Atoms, 10, 130
2022
-
[37]
2016, Phys
Marketin, T., Huther, L., & Martinez-Pinedo, G. 2016, Phys. Rev. C, 93, 025805
2016
-
[38]
M., Sprouse, T
Miller, J. M., Sprouse, T. M., Fryer, C. L., et al. 2019, arXiv e-prints, arXiv:1912.03378
2019 arXiv
-
[39]
2024, Phys
Peng, Y ., Ristic, M., Kedia, A., & et al. 2024, Phys. Rev. Res., 6, 033078
2024
-
[40]
2022, MNRAS, 513, 5174
Pognan, Q., Jerkstrand, A., & Grumer, J. 2022, MNRAS, 513, 5174
2022
-
[41]
2018, Astrophys
Radice, D., Perego, A., Hotokezaka, K., et al. 2018, Astrophys. J, 869, 130 Radži¯ut˙e, L., Gaigalas, G., Kato, D., Rynkun, P., & Tanaka, M. 2020, Astrophys. J. Suppl., 248, 17 Radži¯ut˙e, L., Gaigalas, G., Kato, D., Rynkun, P., & Tanaka, M. 2021, Astrophys. J. Suppl., 257, 29
2018
-
[42]
K., et al
Rosswog, S., Liebendörfer, M., Thielemann, F. K., et al. 1999, A&A, 341, 499
1999
-
[43]
1996, Astron
Ruffert, M., Janka, H.-T., & Schaefer, G. 1996, Astron. Astrophys., 311, 532
1996
-
[44]
2022, Astron
Rynkun, P., Banerjee, S., Gaigalas, G., & et al. 2022, Astron. Astrophys., 658, A82
2022
-
[45]
2019, Nature, 569, 241
Siegel, D., Barnes, J., & Metzger, B. 2019, Nature, 569, 241
2019
-
[46]
Sobolev, V . V . 1960, Moving Envelopes of Stars (Harvard University Press)
1960
-
[47]
& Hotokezaka, K
Tanaka, M. & Hotokezaka, K. 2013, ApJ, 775, 113
2013
-
[48]
2020, MNRAS, 496, 1369
Tanaka, M., Kato, D., Gaigalas, G., & Kawaguchi, K. 2020, MNRAS, 496, 1369
2020
-
[49]
2018, Astrophys
Wanajo, S., Müller, B., Janka, H.-T., & Heger, A. 2018, Astrophys. J, 853, 40
2018
-
[50]
2014, Astrophys
Wanajo, S., Sekiguchi, Y ., Nishimura, N., et al. 2014, Astrophys. J. Lett, 789, L39
2014
-
[51]
& Burrows, A
Wang, T. & Burrows, A. 2023, Astrophys. J, 954, 114
2023
-
[52]
J., Selsing, J., et al
Watson, D., Hansen, C. J., Selsing, J., et al. 2019, Nature, 574, 497
2019
-
[53]
O., Kushnir, D., & Gal-Yam, A
Waxman, E., Ofek, E. O., Kushnir, D., & Gal-Yam, A. 2018, MNRAS, 481, 3423
2018
-
[54]
2012, Astrophys
Winteler, C., Käppeli, R., Perego, A., et al. 2012, Astrophys. J, 750, L22
2012
-
[55]
P., & Cowley, C
Wyart, J.-F., Blaise, J., Bidelman, W. P., & Cowley, C. R. 1997, Physica Scripta, 56, 446 Article number, page 11 of 31 A&A proofs: manuscript no. paper Appendix Article number, page 12 of 31 J. Deprince et al.: Kilonova ejecta opacity inferred from new large-scale HFR atomic ...
1997
Reviewed August 11, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.