Pith. sign in

REVIEW 3 major objections 6 minor 31 references

Benchmarking the Plane-Wave Born and Distorted Waves approximations for electron-impact collision strength computations: the sample case of Sr II

T0 review · 3 major / 6 minor · reviewed 2026-08-04 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read The paper benchmarks two low-cost approximations for electron-impact excitation in Sr II against close-coupling R-matrix data, finding they reproduce the correct order of magnitude of effective collision strengths, with distorted waves with

desk verdict Useful benchmark of cheap collisional approximations for Sr II, but the R-matrix reference's 5–8% level-energy error is never propagated into the central accuracy claim, so the factor-2.5 headline is not yet established. read the letter →

arxiv 2509.07684 v2 pith:5GUUEWNF submitted 2025-09-09 physics.atom-ph astro-ph.HEphysics.comp-ph

Jérôme Deprince , Lucas Maison This is my paper · ORCID
classification physics.atom-phastro-ph.HEphysics.comp-ph
keywords kilonovaelectron-impactexcitationcollisionstrengthsSrIIplane-waveBorndistortedwavesforbiddentransitionsnebularphase
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper tests whether two computationally cheap approximations, the plane-wave Born (PWB) method implemented in the Hartree-Fock code and the distorted-waves (DW) method implemented in AUTOSTRUCTURE, can replace much costlier close-coupling R-matrix calculations for electron-impact excitation in Sr II. Benchmarked against published R-matrix effective collision strengths for the five lowest energy levels, both approximations give the right order of magnitude at the kilonova-relevant temperature of 1000 K, and the DW results agree within a factor of about 2.5. This matters because nebular-phase kilonova modeling needs collisional atomic data for dozens of heavy elements, for which complete R-matrix data are not available. The paper argues that these fast methods can supply large-scale datasets that are far better than the crude empirical formulas currently used, especially for forbidden transitions.

What carries the argument

The central objects are the collision strength and its Maxwellian-averaged effective collision strength, computed with two continuum-electron models: the Plane-Wave Born approximation, which neglects the interaction between the continuum electron and the target potential (with an empirical correction in the Cowan HFR code), and the Distorted-Waves approximation, which includes a static scattering potential but neglects channel coupling and resonances. These approximate rates are benchmarked against CC/R-matrix results at 1000 K for the five lowest levels of Sr II, the comparison being the mechanism that supports the claim.

What would settle it

Compute HFR-PWB and AS-DW effective collision strengths for a lanthanide ion such as Ce II or Sm II for its five lowest levels at 1000 K and compare with an R-matrix close-coupling calculation; if a substantial fraction of transitions deviate by more than an order of magnitude from the R-matrix values, the paper's central claim of reliable large-scale applicability to heavy elements would be falsified for those elements.

Watch

Extended reading notes

Core claim

Using Sr II as a sample case, the authors compute effective electron-impact collision strengths for the ten transitions among the five lowest energy levels, at temperatures below 10,000 K, with two approximate treatments of the continuum electron: plane-wave Born in the pseudo-relativistic Hartree-Fock method, and distorted waves in AUTOSTRUCTURE. Comparing to the CC/R-matrix calculations of Mulholland et al. [18], they find that the DW results agree within a factor of 2.5 or less, and even the PWB results capture the correct order of magnitude in all but one isolated transition. They also compute radiative rates for the two forbidden M1/E2 transitions from the 4d metastable levels, obtainin

Load-bearing premise

The conclusion that these cheap approximations work well enough is drawn from one alkali-like ion with a single valence electron and five low-lying levels, and is assumed to extend to lanthanides and actinides with open f-shell structures, dense level mixing, and strong resonances.

Editorial extensions

If this is right

  • Large-scale collisional data for a wide range of heavy ions can be produced at modest computational cost, replacing empirical formulas that can be off by several orders of magnitude for forbidden transitions.
  • The existing HFR atomic targets for elements from Z=20 to Z=103 can be directly reused to compute PWB collision strengths without constructing new models from scratch.
  • Nebular-phase non-LTE spectral models of kilonovae can incorporate these approximate collision rates, potentially improving line identifications and abundance estimates.
  • For applications where a factor of 2.5 in collision strength is acceptable, AS-DW offers a reliable compromise between accuracy and completeness.
  • The accurate HFR radiative rates for the Sr II forbidden lines support their use in kilonova spectral analysis and in atomic databases.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The benchmark result may not transfer directly to lanthanides and actinides, whose open f-shells produce dense level structures and strong resonances; the paper provides no test in such ions, so the extrapolation is an unproven but testable expectation.
  • The isolated PWB failure for the 2–3 forbidden transition suggests that for some low-energy transitions the plane-wave treatment may be unreliable even in simple ions, so a safety margin or a hybrid approach (PWB for allowed lines, DW for forbidden lines) might be prudent.
  • The same benchmarking protocol could be applied to a heavier ion with many low-lying levels, such as Ce II or Sm II, to map where the factor-of-2.5 agreement breaks down.
  • The effective collision strength comparisons at 1000 K are a single temperature point; extending the benchmark over a wider temperature range would clarify whether the approximations degrade at higher energies where resonances differ.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 6 minor

Summary. The paper benchmarks two computationally inexpensive approximations—HFR with Plane-Wave Born (HFR-PWB) and AUTOSTRUCTURE with Distorted Waves (AS-DW)—for computing electron-impact excitation effective collision strengths in Sr II. The benchmark case uses the five lowest levels and ten transitions, comparing against published CC/R-matrix calculations [18] at T=1000 K (Table 3). The paper reports that AS-DW agrees with the R-matrix reference within a factor of about 2.5 and that HFR-PWB generally reproduces the correct order of magnitude, with one isolated outlier. It also presents new forbidden M1/E2 radiative rates for the 4d-5s transitions. The stated motivation is to validate cheap methods for large-scale production of collisional data for heavy elements needed in nebular-phase kilonova modeling.

Significance. If the reported accuracy were robust, the paper would provide a practical pathway for filling a genuine gap in atomic data for non-LTE kilonova spectral modeling. The HFR+Fit target energies (0.02% agreement with NIST, Table 1) and the radiative rates (within 4% for HFR+Fit, Table 2) are useful results in their own right. The comparison with R-matrix data is honest in that no parameter is fitted to the reference collision strengths, and the paper explicitly flags the HFR-PWB outlier. However, the central quantitative benchmark is weakened by an energy-threshold inconsistency in the reference data and by the narrow scope of the extrapolation, as detailed below.

major comments (3)
  1. [Table 1, Eq. (2), Table 3] The R-matrix reference [18] and the present approximate calculations use different target energy thresholds, and this difference is large enough to affect the T=1000 K benchmark. Table 1 shows that the AS energies used in [18] for the 4d levels are 15686 and 16001 cm^-1 versus NIST values of 14556 and 14836 cm^-1, i.e. offsets of 1130 and 1165 cm^-1. At T=1000 K, kT is about 695 cm^-1, so the Boltzmann factor in Eq. (2) differs by roughly exp(1.6)~5 for the 1-2 and 1-3 excitations. The reference effective collision strengths are therefore suppressed relative to an energy-consistent target, and the ratios in Table 3 mix target-structure errors with the errors of the PWB/DW continuum approximations. For example, AS-DW is 2.16 times the reference for 1-4, but if the reference threshold is corrected upward by the exp(-deltaE/kT) factor the ratio changes substantially. The claim that AS-DW ag
  2. [Section 6, Conclusions] The extrapolation from Sr II to lanthanides and actinides is not supported by the evidence in the manuscript. The benchmark covers only five levels of one alkali-like ion with a single valence electron. Open f-shell lanthanides and actinides have much higher level densities, stronger channel coupling, and more prominent resonance effects, all of which are neglected in PWB and DW. The sentence in the Conclusions stating that the benchmark 'demonstrates that HFR-PWB and AS-DW approximations can provide reliable large-scale data sets ... for systematic applications to lanthanides, actinides' overstates what a single sample case can show. Either add a representative heavy open-f-shell test case or restrict the conclusion to Sr II and clearly state that transferability to complex ions remains to be validated.
  3. [Table 3, Section 5] The HFR-PWB outlier 2-3 (4d 2D3/2 - 4d 2D5/2) is a factor of about 143 below the R-matrix value (9.50E-02 versus 1.36E+01). It is not merely a numerical curiosity: this is a forbidden transition within the ground configuration and is likely to be relevant for nebular-phase level populations. The text repeatedly summarizes the HFR-PWB accuracy as 'correct order of magnitude' and calls the failure 'isolated,' but with only ten transitions tested, one failure is a 10% failure rate on the most physically important class of transitions. The conclusion that HFR-PWB can provide reliable large-scale data should be qualified, and the manuscript should either explain the origin of this discrepancy or soften the recommendation.
minor comments (6)
  1. [Eq. (2)] The integration variable in the definition of the effective collision strength is written as epsilon_j. The standard definition integrates over the incident electron energy; please clarify the notation and identify which energy the subscript j refers to.
  2. [Section 5] Please state explicitly whether the collision calculations use the HFR or HFR+Fit target energies, and whether the AS calculations include the TEC. Table 1 gives both HFR and HFR+Fit, so the reader cannot tell which target is used in Table 3 without inferring from context.
  3. [Table 2 caption] There is a typo: 'metatable' should be 'metastable', and 'in comparison the those' should be 'in comparison with those'.
  4. [Section 3, Table 1] The statement that 'results from [18] show a difference of 5.29% with NIST values on average' is not directly supported by the five values in Table 1, which give an average absolute deviation of about 4.2%. Please specify the set of levels over which the 5.29% average is computed.
  5. [Eq. (3)] The empirical correction factor (X+3)/(1+X) is introduced without discussion of its range of validity or the data from which it was derived. A sentence with a reference or a qualitative explanation would help the reader assess whether its use is appropriate for low-temperature effective collision strengths.
  6. [General] No uncertainties or convergence checks are reported for the collision strengths in Table 3. Since the central claim is a factor-level accuracy statement, the paper should at least discuss numerical uncertainties of the present calculations and of the R-matrix reference.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the benchmark uses an independent R-matrix reference and no fitted collision-strength inputs are recycled as predictions.

full rationale

The paper's central claim is that HFR-PWB and AS-DW effective collision strengths reproduce the order of magnitude of CC/R-matrix results for Sr II, with AS-DW within ~2.5. This is tested against the independent R-matrix calculation of Mulholland et al. [18]. The quantities in Table 3 are not fitted to those reference values: the authors improve their target energy levels using experimental NIST energies (HFR+Fit, AS-TEC), which are independent of the collision-strength benchmark; the Cowan empirical scaling in Eq. (3) is a fixed formula from Cowan's book [15], not calibrated to the [18] collision strengths in this paper. The HFR and AS structure models are used to generate the collision data from first principles (within the approximations), and the R-matrix data are then used only as a comparison. The only self-citations (e.g., [11] for existing HFR targets) are practical, not load-bearing for the accuracy conclusion. The skeptic concern about the R-matrix reference having 5–8% high thresholds affects the reliability or interpretation of the benchmark, but it does not constitute circularity: the paper does not define its predictions in terms of the reference, nor does it fit any parameter to the benchmark data. No specific circular reduction can be exhibited from the text. Therefore the appropriate finding is no circularity, score 0.

Assumptions & free parameters 3 free parameters · 4 assumptions · 0 invented entities

The computations rely on standard atomic structure theory (Hartree-Fock with Slater-Condon, Thomas-Fermi-Dirac-Amaldi model potentials, distorted wave scattering theory). The free parameters tune the structure model to experimental energies, but the collision strength benchmark is against independent R-matrix data, so circularity burden is low.

free parameters (3)
  • HFR radial parameter scaling (semi-empirical fit) = slater parameters scaled to NIST energies (not tabulated)
    Used to bring HFR energies to within 0.02% of NIST; affects the target wavefunctions used in Table 1 and collision strength calculations.
  • AUTOSTRUCTURE TFDA scaling parameters = optimized values (not tabulated)
    Scaling parameters for each orbital optimized in a multiconfiguration variational procedure; energy levels then TEC-corrected. These tune the structure model.
  • Term Energy Correction (TEC) values = applied to AS energy levels (not tabulated)
    AS correction to match experimental levels; shifts the level structure and hence transition energies.
assumptions (4)
  • standard math Validity of Slater-Condon theory and HF variational method for the Sr II target description
    Invoked in Sections 2.1 and 2.2 as the basis of both structure methods.
  • domain assumption Distorted-wave and plane-wave approximations are meaningful at nebular-phase electron temperatures (1000-10000 K)
    Section 2.3 defines the collision strength and assumes these scattering approximations for the continuum electron; the paper benchmarks rather than derives them.
  • domain assumption Only the first five energy levels (below 4 eV) can be populated in the KN nebular phase
    Section 3, citing Gillanders & Smartt [22]; justifies truncating the analysis to five levels.
  • domain assumption R-matrix close-coupling results of Mulholland et al. [18] are an accurate reference benchmark
    Section 5 and Table 3 use these as the accuracy standard; no independent validation is provided.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Benchmarking the Plane-Wave Born and Distorted Waves approximations for electron-impact collision strength computations: the sample case of Sr II." pith.science (2026). https://pith.science/paper/5GUUEWNF

@misc{pith2026250907684,
  author       = {Pith},
  title        = {Pith review of: Benchmarking the Plane-Wave Born and Distorted Waves approximations for electron-impact collision strength computations: the sample case of Sr II},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/5GUUEWNF}},
  note         = {Machine review of arXiv:2509.07684}
}
read the original abstract

The discovery of gravitational waves from a neutron star merger in 2017 (GW170817) and the associated kilonova (AT2017gfo) confirmed these events as key sites for heavy element production through the r-process. Subsequent observations, including late-time spectra with JWST, have highlighted the need for accurate modeling of kilonova ejecta. In the photospheric phase, atomic level populations can be estimated under LTE using Boltzmann and Saha relations, but about a week after the merger the ejecta enters the nebular phase where non-LTE effects dominate. Modeling nebular spectra therefore requires a detailed treatment of radiative and collisional processes that affect the population of atomic levels. This work focuses on electron-impact excitation in Sr II, a heavy ion relevant for kilonova spectra. Two computational approaches are employed: the Plane Wave Born approximation within the pseudo-relativistic Hartree-Fock method, and a Distorted Waves method using AUTOSTRUCTURE. The resulting collision strengths are compared against reference R-matrix data to evaluate the accuracy of these approximations and their suitability for large-scale applications to all heavy elements. In addition, radiative parameters for forbidden transitions are computed. These results provide an essential benchmark of approximations that could be used to compute atomic data for nebular-phase kilonova modeling.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

31 extracted references · 30 canonical work pages

  1. [18]

    We find that our AS transition rates underestimate the NIST values by about 10%, while the results from [18] overestimate them by 10%

    and to experimental levels from [23, 24] and available in the NIST database [19] is given in Table 2. We find that our AS transition rates underestimate the NIST values by about 10%, while the results from [18] overestimate them by 10%. Nevertheless, our HFR transition probabilities are found to be much closer to the observed values from NIST, since the a...

    work page

  2. [1]

    Abbott, R

    B. Abbott, R. Abbott, T. Abbott et al., Phys. Rev. Lett.119, 161101 (2017)

    work page 2017

  3. [2]

    Kasen, B

    D. Kasen, B. Metzger, J. Barnes et al., Nature551, 80 (2017)

    work page 2017

  4. [3]

    Given the acute lack of atomic data for the wide range of heavy elements expected in KN ejecta, the balance be- tween completeness and accuracy becomes crucial

    and to the van Regemorter formula. Given the acute lack of atomic data for the wide range of heavy elements expected in KN ejecta, the balance be- tween completeness and accuracy becomes crucial. The present benchmark demonstrates that HFR-PWB and AS- DW approximations can provide reliable large-scale data sets of collisional parameters at relatively low ...

    work page

  5. [4]

    A. J. Levan, B. P. Gompertz, O. S. Salafia et al., Nature626, 737 (2024)

    work page 2024

  6. [5]

    Metzger, G

    B.D. Metzger, G. Martinez-Pinedo, S. Darbha et al., MN- RAS406, 2650 (2010)

    work page 2010

  7. [6]

    Watson, C

    D. Watson, C. J. Hansen, J. Selsing et al., Nature574, 497 (2019)

    work page 2019

  8. [7]

    Tanaka, D

    M. Tanaka, D. Kato, G. Gaigalas, K. Kawaguchi, MNRAS 496, 1369 (2020)

    work page 2020

Show all 31 references
  1. [8]

    C. J. Fontes, C. L. Fryer, A. L. Hungerford, A. L. et al., MN- RAS493, 4143 (2020)

    2020

  2. [9]

    Gaigalas, D

    G. Gaigalas, D. Kato, P. Rynkun, L. Radziute, M. Tanaka, Astrophys. J. Suppl.240, 29 (2019)

    2019

  3. [10]

    Fl ¨ors, R

    A. Fl ¨ors, R. F. Silva, J. Deprince et al. MNRAS524, 3083 (2023)

    2023

  4. [11]

    Deprince, G

    J. Deprince, G. Wagle, S. Ben Nasr et al., A&A696, A32 (2025)

    2025

  5. [12]

    C. J. Fontes, C. L. Fryer, C. L., R. T. Wollaeger et al., MN- RAS519, 2862 (2023)

    2023

  6. [13]

    M. M. Kasliwal, D. Kasan, R. M. Lau, MNRAS510, L7 (2022)

    2022

  7. [14]

    Gillanders, S.J

    J.H. Gillanders, S.J. Smartt, S.A. Sim et al., MNRAS515, 631 (2022)

    2022

  8. [15]

    Pognan, A

    Q. Pognan, A. Jerkstrand, J. Grumer, MNRAS510, 3806 (2022)

    2022

  9. [16]

    Badnell, CPC182, 1528 (2011)

    N. Badnell, CPC182, 1528 (2011)

    2011

  10. [17]

    P. G. Burke, R-Matrix Theory of Atomic Collisions: Ap- plication to Atomic, Molecular and Optical Processes, Springer, Berlin (2011)

    2011

  11. [19]

    For AS, a Term Energy Correction (TEC) [16] was applied to opti- mize the energy levels

    was also performed with HFR to obtain a better agree- ment and correct our radial parameter values. For AS, a Term Energy Correction (TEC) [16] was applied to opti- mize the energy levels. Since we are interested in the KN ejecta nebular phase that is characterized by a relati...

  12. [20]

    R. D. Cowan, The Theory of Atomic Structure and Spectra, University of California Press (1981)

    1981

  13. [21]

    L. P. Mulholland, N. E. McElroy, F. L. McNeill et al., MN- RAS532, 2289 (2024)

    2024

  14. [22]

    Kramida, Yu

    A. Kramida, Yu. Ralchenko, J. Reader and NIST ASD Team (2024). NIST Atomic Spectra Database (ver. 5.12), [Online]. Available: https://physics.nist.gov/asd [2025, August 19]. National Institute of Standards and Technology, Gaithers- burg, MD

    2024

  15. [23]

    Moore, Atomic Energy Levels As Derived From the Analyses optical Spectra, V ol

    C. Moore, Atomic Energy Levels As Derived From the Analyses optical Spectra, V ol. II, National Standard Refer- ence Data Series. National Bureau of Standards, U.S., p. 259 (1952)

    1952

  16. [24]

    J. E. Sansonetti, J. Phys. Chem. Ref. Data41, 013102 (2012)

    2012

  17. [25]

    J. H. Gillanders, S. J. Smartt, MNRAS538, 1663 (2025)

    2025

  18. [26]

    Letchumanan, M

    V . Letchumanan, M. A. Wilson, P. Gill, A. G. Sinclair, Phys. Rev. A72, 012509 (2005)

    2005

  19. [27]

    Bi ´emont, J

    E. Bi ´emont, J. Lidberg, S. Mannervik et al., Eur. Phys. J. D 11, 355 (2000)

    2000

  20. [28]

    Pognan, J

    Q. Pognan, J. Grumer, A. Jerkstrand, S. Wanajo, MNRAS 526, 5220 (2023)

    2023

  21. [29]

    van Regemorter, ApJ136, 906 (1962)

    H. van Regemorter, ApJ136, 906 (1962)

    1962

  22. [30]

    T. S. Axelrod, PhD thesis, California University, Santa Cruz (1980)

    1980

  23. [31]

    Pognan, M

    Q. Pognan, M. R. Wu, G. Martinez-Pinedo et al., MNRAS 536, 2973 (2025). I-20/A-11-04

    2025

Pith tools

Reviewed August 4, 2026 · model on record in the stance chip above.