Pith. sign in

REVIEW 2 major objections 5 minor 35 references

Simulation study of neutral tungsten emissions for fusion applications

T0 review · 2 major / 5 minor · reviewed 2026-08-05 · deepseek-v4-flash

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

Pith's one-line read This paper computes electron-impact excitation cross sections and rate coefficients for three neutral-tungsten lines (400.87, 429.46, 430.21 nm), reporting the latter two for the first time.

desk verdict Useful new rate data for two WI lines, but the 430.21 nm result rests on an untested wavefunction assumption after a large energy shift. read the letter →

arxiv 2509.00878 v1 pith:OAPGO7EY submitted 2025-08-31 physics.plasm-ph physics.atom-ph

Ritu Dey , Ayushi Agrawal , Reetesh Kumar Gangwar , Deepti Sharma , Rajesh Srivastava , Malay B. Chowdhuri , Joydeep Ghosh This is my paper · ORCID
classification physics.plasm-phphysics.atom-ph PACS 52.20.Fs52.55.Fa
keywords tungstenelectron-impactexcitationrelativisticdistortedwaveflexibleatomiccoderatecoefficientstokamakdiagnosticsWIemissionlinesfusionplasma
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

This paper calculates electron-impact excitation cross sections and Maxwellian-averaged rate coefficients for three visible emission lines of neutral tungsten — 400.87, 429.46, and 430.21 nm — using the relativistic distorted-wave (RDW) approximation inside the Flexible Atomic Code. Two of the three datasets are reported for the first time. The motivation is practical: tungsten is the wall and divertor material of current and next-generation tokamaks, and the brightness of neutral-tungsten lines is used to measure wall erosion, so accurate excitation data feed directly into spectroscopic diagnostics. The authors also demonstrate that correcting the code's energy levels to match measured NIST values is essential — uncorrected calculations move the rate coefficients by up to roughly 50 percent. If the results hold up, diagnostics gain a consistent set of excitation data across the full temperature range of fusion plasmas.

What carries the argument

The relativistic distorted-wave (RDW) approximation implemented in the Flexible Atomic Code (FAC): a perturbative treatment of electron–atom collisions in which the projectile is described by distorted waves and the target states come from diagonalizing a Dirac–Coulomb Hamiltonian over 1754 fine-structure levels. The collision strength is built from Slater integrals between initial and final target states, converted to a cross section, then patched near threshold with the factor [1 − (E_threshold/E)^3], and finally Maxwell-averaged over electron energy to give rate coefficients. The level-shift correction that moves the code's energies onto the measured NIST values is the step that makes the

What would settle it

Measure the absolute electron-impact excitation cross sections for the 429.46 and 430.21 nm lines with a crossed-beam experiment, or recalculate them with a non-perturbative close-coupling method such as DARC; if either disagrees with the RDW curves by more than the factor-of-two spread already seen at 400.87 nm between the two theoretical methods, the corrected-RDW claim is refuted. A cheaper check is to compare near-threshold rate coefficients against measured or R-matrix data below 30 eV, since the present 400.87 nm rates already run below a prior calculation in exactly that range.

Watch

Extended reading notes

Core claim

The paper's central claim is that the relativistic distorted-wave method, running inside the Flexible Atomic Code, yields usable electron-impact excitation cross sections and rate coefficients for three neutral-tungsten emission lines — 400.87, 429.46, and 430.21 nm — provided the code's energy levels are first corrected to the measured NIST values. For the 429.46 and 430.21 nm lines the data are reported for the first time; the 400.87 nm results are new but comparable against an existing calculation. The uncorrected FAC levels are badly placed — the common lower level sits near 2.66 eV instead of 0.365 eV, and the 430.21 nm upper level even falls below it — and the authors show that correct

Load-bearing premise

That the theoretical wavefunctions for the excited states remain trustworthy after the code's energy levels are forcibly shifted to the measured NIST values, even though the uncorrected levels are far off — the lower level of all three lines starts near 2.7 eV instead of 0.365 eV, and the 430.21 nm upper level sits below the lower one.

Editorial extensions

If this is right

  • Spectroscopic diagnostics of neutral tungsten in tokamaks gain excitation data for the 429.46 and 430.21 nm lines, which had no previously reported cross sections or rate coefficients.
  • The roughly 50 percent shift in rate coefficients between corrected and uncorrected energy levels demonstrates that matching the code to measured levels is a required step in future tungsten atomic-data calculations.
  • With incident energies to 30 keV and temperatures to 300 eV, the data cover conditions from the divertor to the core, so the same tables can serve multiple diagnostic regimes.
  • The lower 400.87 nm rates relative to the earlier calculation give a quantitative benchmark showing how much the choice of theoretical method matters for neutral tungsten.

Reading between the lines

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

  • If the FAC wavefunctions for the excited states are as distorted as their uncorrected energies suggest — one upper level even falls below its lower level — the corrected cross sections could carry a systematic error of the same order as the ~50 percent shift the correction removes. The paper does not test this directly.
  • Applying the same level-correction recipe to the other visible WI lines already studied in the literature (488.69, 498.26, 522.47 nm) would produce a consistent RDW dataset across the full set of lines used in erosion diagnostics.
  • The rates that matter most for divertor conditions (1–50 eV) rely on the empirical [1 − (E_th/E)^3] near-threshold patch, which is unbenchmarked for these transitions; anchoring it with a close-coupling calculation or experiment would be a direct test.
  • Because the two new lines are also seen in laser-induced breakdown spectroscopy, a laboratory measurement of their relative line intensities could serve as a low-cost check of the computed rate-coefficient ratios.
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

2 major / 5 minor

Summary. The manuscript reports relativistic distorted wave (RDW) calculations, carried out with the Flexible Atomic Code (FAC), of electron-impact excitation cross-sections and Maxwellian-averaged rate coefficients for three neutral tungsten (W I) transitions: 400.87 nm, 429.46 nm, and 430.21 nm. The lower level of all three is 5d^5(^6S)6s ^7S_3. The authors state that FAC energy levels were corrected to match NIST values and that, to their knowledge, the cross-sections for 429.46 nm and 430.21 nm are reported for the first time. The cross-sections are given for incident electron energies up to 30 keV, and rate coefficients are presented for electron temperatures up to about 250 eV. A low-energy correction from prior literature is applied near threshold, and the effect of the energy-level correction on rate coefficients is shown graphically.

Significance. If the results are reliable, the paper fills a small but real gap: visible W I excitation data for two transitions not previously computed, relevant to tokamak and LIBS diagnostics. The use of FAC/RDW is standard, the energy-level anchoring to NIST is conceptually sensible, and the Maxwellian averaging is transparent. The main strengths are the explicit treatment of a complex open-shell target and the direct provision of rate coefficients in a usable range. However, the paper's central claim rests on the quality of the FAC target-state wavefunctions after the NIST energy correction, and this is not validated. Because Table 2 shows that the uncorrected FAC model is seriously wrong for one of the three transitions, the missing wavefunction validation is a load-bearing issue, not a presentation concern.

major comments (2)
  1. [§3, comparison with Ref. [18]] The only external comparison is for the 400.87 nm transition, and it is qualitative: the present rates are said to be 'lower' than the Dirac R-matrix/MCDF rates of Kwon et al. for energies up to 30 eV. A quantitative comparison is needed, e.g., a ratio or a plot overlaying the two rate coefficients as a function of Te. This is particularly important because RDW and R-matrix methods often differ near threshold, and the magnitude of the difference will tell the reader whether the discrepancy is within the expected accuracy or indicates a target-state problem.
  2. [§3, near-threshold region and Eq. (5)] The reported cross-section maxima occur at incident energies close to threshold (5.0, 4.5, and 4.0 eV), exactly where the RDW method is least reliable and where the ad hoc correction of Eq. (5), [1 - (E_th/E)^3], is applied. The paper does not quantify how much Eq. (5) changes the cross sections or rate coefficients, nor does it provide uncorrected and corrected cross sections for comparison. Since applications in fusion edge plasmas and LIBS sample low-energy electrons, the uncertainty introduced by this correction should be stated. Please show the effect of Eq. (5) on the reported maxima and, if possible, benchmark the corrected results against another method in the threshold region.
minor comments (5)
  1. [Fig. 1 caption] Collision strength is dimensionless; the vertical axis label 'Collision Strength (a.u.)' is misleading.
  2. [Throughout] Typographical errors: 'collison' (Fig. 1), 'transtions' (§3), 'addtion' (§3), 'KeV' (§3), and 'Driac-R matrix' should be 'Dirac R-matrix'.
  3. [Fig. 4 caption] The caption text says 'blue-dashed curve' for the corrected curves, while the legend text says 'dash-dotted'; please make the line styles and captions consistent.
  4. [References] Reference [19] is a private communication about LIBS observations. If this is the only experimental motivation for the 429.46 nm and 430.21 nm lines, a citable source or more detail should be provided.
  5. [Data availability] The figures are difficult to use quantitatively. Please consider providing tabulated cross-sections and rate coefficients as supplementary material, which would increase the utility of the data for the fusion community.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the reported cross-sections and rate coefficients are independent outputs of an RDW/FAC calculation, with NIST energies used only as level/threshold corrections, not as fitted targets.

full rationale

The paper's central claim is the simulation of electron-impact excitation cross-sections and Maxwellian rate coefficients for three WI lines using FAC's RDW implementation (Eqs. 1-4, 6). The cross-sections are built from Dirac-Coulomb target wavefunctions and collision-strength matrix elements; no measured cross-section or rate coefficient is used as an input. NIST data enter only through the stated correction of energy levels (Table 1, Section 3) and through the threshold energy in the fixed literature correction Eq. 5, which is attributed to Ref. [20] and merely paralleled by Ref. [21]. The transition wavelengths are identities of the levels, not predictions, so matching NIST energies is not a claimed derived result. The external comparison with Kwon et al. is qualitative, but a discrepancy is an accuracy concern, not a circularity. Table 2 does expose a serious internal risk: the uncorrected FAC model places the 430.21 nm upper level 0.59 eV below the lower level, and an a posteriori energy shift does not rebuild the eigenvectors entering Eq. (1); however, this is an unvalidated-wavefunction/validation gap, not a reduction of the prediction to its inputs. Refs. [8] and [21] may be prior work by the same group, but they are non-load-bearing method citations, and FAC itself is an external code. No circular step can be exhibited from the paper's own equations or claims.

Assumptions & free parameters 0 free parameters · 5 assumptions · 0 invented entities

The central claim (EIE cross sections and rates for three WI lines) rests on the validity of the RDW approximation for neutral tungsten, the fidelity of FAC wavefunctions after NIST energy shifting, and the applicability of an empirical near-threshold correction. No free parameters are fitted to data, and no new entities are introduced.

assumptions (5)
  • domain assumption The relativistic distorted wave approximation is valid for computing electron-impact excitation of neutral W at the considered energies.
    Invoked in Section 2 (Eqs. 1-4) without justification specific to W I; the authors themselves note RDW is perturbative and less accurate near threshold.
  • domain assumption The Dirac-Coulomb Hamiltonian with the listed configurations, after shifting energies to NIST values, yields wavefunctions accurate enough for cross-section calculations.
    Section 3 states 1754 fine-structure levels are considered and energies are corrected to NIST; Table 2 shows large uncorrected deviations, so wavefunction reliability is assumed.
  • domain assumption The empirical low-energy correction (1 - (Ethreshold/E)^3) from Wunderlich et al. applies to neutral tungsten.
    Eq. (5) is taken from prior literature and applied without validation against measurement for W I.
  • domain assumption Maxwellian averaging (Eq. 6) is the appropriate way to derive rate coefficients for fusion plasma conditions.
    Standard practice in the field, stated in Section 2.
  • standard math NIST energy level values are accurate references.
    NIST data are used as the benchmark for level corrections in Section 3 and Table 2.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Simulation study of neutral tungsten emissions for fusion applications." pith.science (2026). https://pith.science/paper/OAPGO7EY

@misc{pith2026250900878,
  author       = {Pith},
  title        = {Pith review of: Simulation study of neutral tungsten emissions for fusion applications},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/OAPGO7EY}},
  note         = {Machine review of arXiv:2509.00878}
}
read the original abstract

The article reports electron-impact excitation cross-sections and rate coefficients for neutral tungsten for three transitions (400.87 nm, 429.46 nm, and 430.21 nm) using the relativistic distorted wave approach within the flexible atomic code. Some of these lines are also observed in tokamak plasma. Cross-sections are computed for incident electron energy up to 30 keV. The energy levels in flexible atomic code were corrected to match the NIST database. The electron impact excitation rate coefficients are also provided.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

35 extracted references · 35 canonical work pages

  1. [1]

    The configurations that are taken for the present cal- culations are as follows: 5d 46s2, 5d 56[s,p], 5d 46s6[p,d], 5d46s7[s,p] and 5d 6. It is noticed that all three transi- tions are dipole allowed transitions and the transition prob- abilities are 1.63×107 s−1, 1.24×107 s−1, and 3.60×106 s−1, for wavelengths 400.87 nm, 429.46 nm, and 430.21 nm, respect...

    work page

  2. [2]

    0.365 3.252 430.21 5d5(6S)6s (7S 3) 5d56s(6D)6p (7D0

    work page

  3. [3]

    While for 430.21 nm wavelength, the exci- tation cross-section shows a maximum 2.50× 10−18 cm2 at ∼ 4.0 eV and falls∼ 103 times at incident electron energy of 30 keV

    0.365 3.247 eV for the transition wavelengths, 400.87 nm, 429.46 nm, respectively, and it falls∼ 103 times at incident electron en- ergy of 30 KeV . While for 430.21 nm wavelength, the exci- tation cross-section shows a maximum 2.50× 10−18 cm2 at ∼ 4.0 eV and falls∼ 103 times at incident electron energy of 30 keV . In addtion, Figure 3 depicts the EIE rat...

    work page

  4. [4]

    0.365 3.457 429.46 5d5(6S)6s (7S 3) 5d5(6S)6p (7P0

    work page

  5. [7]

    Finally, Figure 4 illustrates the EIE rate coe fficients for the above mentioned neutral W wavelengths, simulated with and without corrected excited energy levels

    (red dash-dotted curve), respectively. Finally, Figure 4 illustrates the EIE rate coe fficients for the above mentioned neutral W wavelengths, simulated with and without corrected excited energy levels. The en- ergy levels before correction for these three particular tran- E (eV) 10 10 2 10 3 10 4 Cross section (cm 2) 10 -22 10 -21 10 -20 10 -19 10 -18 10...

    work page

  6. [8]

    (blue-solid curve), 5d5(6S)6p ( 7P0

    work page

  7. [10]

    Electron temp

    (red dash-dotted curve), respectively. Electron temp. (eV) 25 50 75 100 150 200 250 Rate Coe ffi cient (cm 3/s) 10 -11 10 -10 10 -9 10 -8 7 P0 4 (400.87 nm) 7 P0 2 (429.46 nm) 7 D 0 3 (430.21 nm) Fig. 3 Notations are similar to the figure 1, simulated EIE rate- coefficients for the neutral W. sitions are listed in Table 2. It is noticed that, the maximum de...

    work page

  8. [11]

    It is hereby confirmed that energy level corrections are nec- essary to obtain the accurate collision strengths and cross sections for the neutral W

    and -5d56s(6D)6p (7D0 3)) considered here. It is hereby confirmed that energy level corrections are nec- essary to obtain the accurate collision strengths and cross sections for the neutral W. 4 Conclusion The electron impact excitation cross-sections for WI are simulated with RDW approximation within the FAC code. The Maxwellian averaged rate coe fficien...

    work page

Show all 35 references
  1. [12]

    and 5d 56s(6D)6p ( 7D0 3), respec- tively. The electron impact excitation rate coe fficients for transitions in the visible range are presented, which will contribute significantly to spectroscopic diagnostics in fu- sion devices. 01-02 Electrom temp.(eV) 25 50 75 100 150 200 ...

  2. [13]

    (blue-dashed curve), 5d5(6S)6p ( 7P0

  3. [14]

    (black dashed curve) and 5d 56s(6D)6p (7D0

  4. [15]

    The blue, black and red dash-dotted curves represent the FAC simulated EIE rate coefficients uncor- rected energy levels

    (red dashed curve), respectively for corrected en- ergy levels. The blue, black and red dash-dotted curves represent the FAC simulated EIE rate coefficients uncor- rected energy levels. Table 2 Energy levels of WI transitions from FAC and the NIST

  5. [16]

    Dyall et al, GRASP: a general-purpose relativistic atomic structure program

    K. Dyall et al, GRASP: a general-purpose relativistic atomic structure program. Comput. Phys. Commun. 55, 425 (1989)

    1989

  6. [17]

    R. A. Pitts et al, Nuclear Materials and Energy 20, 100696 (2019)

    2019

  7. [18]

    Sang et al, Nuclear Fusion, 61, 066004 (2021)

    C. Sang et al, Nuclear Fusion, 61, 066004 (2021)

    2021

  8. [19]

    Yihan Wu et al, Nuclear Materials and Energy 33, 101297 (2022)

    2022

  9. [20]

    Bucalossi et al, Nucl

    J. Bucalossi et al, Nucl. Fusion 64, 112022 (2024)

    2024

  10. [21]

    Thoma et al, Plasma Phys

    A. Thoma et al, Plasma Phys. Control. Fusion 39, 1487 (1997)

    1997

  11. [22]

    recommended values based on experimental data. Wavel. Ei (FAC) Ek (FAC) Ei (NIST) Ek (NIST) Lower Upper Lower Upper (nm) (eV) (eV) (eV) (eV) 400.87 2.66 4.01 0.365 3.457 429.46 2.66 3.70 0.365 3.252 430.21 2.66 2.07 0.365 3.247 5 Acknowledgements The authors, Ritu Dey, Reetesh...

    2023

  12. [23]

    G. F. Matthews et al, Phys. Scr. T128, 137-143 (2007)

    2007

  13. [24]

    Beigman et al, Plasma Phys

    I. Beigman et al, Plasma Phys. Control. Fusion 49, 1833 (2007)

    2007

  14. [25]

    Priti et al, Atoms 2015, 3, 53-75 (2015)

    2015

  15. [26]

    Quinet et al, J

    P. Quinet et al, J. Phys. B At. Mol. Opt. Phys. 44, 1883 (2011)

    2011

  16. [27]

    J. P. Desclaux, Computer Physics Communications 9, 31 (1975)

    1975

  17. [28]

    Eugene Rudd, Physical Rev

    Yong-Ki Kim and M. Eugene Rudd, Physical Rev. A 50, 3954 (1994)

    1994

  18. [29]

    R. D. Cowan, The Theory of Atomic Structure and Spec- tra; University of California Press: Berkeley, CA, USA, 1981

    1981

  19. [30]

    Kramida, Atoms 7, 64 (2019)

    A. Kramida, Atoms 7, 64 (2019)

    2019

  20. [31]

    The flexible atomic code

    Ming Feng Gu, “The flexible atomic code”, Can. J. Phys. 86, 675-689, 2008

    2008

  21. [32]

    Ballance, DARC, Dirac atomic R-matrix code (2005)

    C. Ballance, DARC, Dirac atomic R-matrix code (2005)

    2005

  22. [33]

    Indelicato, MDFGME, a multiconfiguration Dirac- Fock and general matrix elements program (2024)

    P. Indelicato, MDFGME, a multiconfiguration Dirac- Fock and general matrix elements program (2024). https://www.lkb.upmc.fr/metrologysimplesystems/mdfgme- a-general-purpose-muticonnfiguration-dirac-foc- program/

    2024

  23. [34]

    Duck-Hee Kwon and Paul Indelicato, Eur. Phys. J. D. 78 86 (2024)

    2024

  24. [35]

    Bharat Hegde et al, private communication, (2025)

    2025

  25. [36]

    W¨underlich et al, Journal of Quantitative Spectroscopy & Radiative Transfer 149, 360 (2014)

    D. W¨underlich et al, Journal of Quantitative Spectroscopy & Radiative Transfer 149, 360 (2014)

    2014

  26. [37]

    Priti et al, Journal of Quantitative Spectroscopy & Radia- tive Transfer 187, 426 (2017)

    2017

  27. [38]

    Kramida, Yu

    A. Kramida, Yu. Ralchenko, J. Reader, and NIST ASD Team (2024). NIST Database (ver. 5.12) . Available: https://physics.nist.gov/asd. 01-03

    2024

Pith tools

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