{"id":"f99fb888-72e6-462d-8c2b-458521030ef1","arxiv_id":"2509.00878","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"New RDW/FAC excitation cross sections and Maxwellian-averaged rate coefficients for three neutral tungsten lines (400.87, 429.46, 430.21 nm), two of which were not previously computed.","lead":"This paper computes new electron-impact excitation cross sections and rate coefficients for three visible neutral tungsten lines using a standard relativistic atomic code. The data are intended to support tungsten erosion diagnostics in tokamak fusion devices.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Post-hoc NIST energy shifts do not repair the FAC target-state wavefunctions; for the 430.21 nm transition the uncorrected upper level sits 0.59 eV below the lower level, so the RDW cross sections inherit an untested wavefunction error.","rationale":"The reader's weakest assumption identifies exactly the load-bearing point: the FAC wavefunctions for the target states may be unreliable even after the NIST energy correction. Table 2 makes this concrete—the uncorrected FAC energies are not merely shifted by a few percent; the 430.21 nm upper level is inverted relative to the lower level. Because the RDW cross section depends on the electric-dipole-like transition matrix elements between these states, a wavefunction error directly propagates into the central numerical results. The paper does not supply the supporting evidence that would retire this concern: no eigenvector convergence study, no comparison of A-values with NIST or with the MCDF/DARC calculation by Kwon et al., and no quantitative comparison for the one overlapping line. The conditional verdict is therefore appropriate. I would not move to reject, because RDW/FAC is an established method and the energy-corrected results are plausible; the issue is one of missing validation, not demonstrated failure. The concrete test proposed—an expanded-CI rerun with eigenvector and A-value comparison—would determine whether the concern lands, and would either confirm the reported data or reveal the sensitivity that the current manuscript omits.","tokens_in":5892,"tokens_out":5062,"duration_ms":71538,"concrete_test":"Re-run the FAC calculation with an expanded active space that includes the dominant correlation configurations omitted from §3 (e.g., 5d4 6p2, 5d5 6d, 5d4 6s6d, 5d4 6s7p), applying the NIST energy correction only after diagonalization. Then compare the eigenvector composition of the 7D3 upper state and the A-values for all three transitions against the present model. If the leading CI weight of 7D3 shifts by more than 0.05 or any A-value shifts by more than 30%, the RDW cross sections are not robust to the target description, and the paper should present the spread instead of single curves. A secondary check: compute the 400.87 nm rate coefficient at 10–30 eV and compare quantitatively with Kwon et al.; a factor-of-two discrepancy would corroborate the wavefunction concern.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's deliverable is new RDW EIE cross sections and Maxwellian rate coefficients for three WI lines. For these data to be usable, the target-state wavefunctions in Eq. (1) must be physically reliable. Table 2 shows that the unshifted FAC Hamiltonian places the lower 5d5(6S)6s 7S3 level at 2.66 eV instead of 0.365 eV, and for the 430.21 nm transition it places the upper 5d5 6s(6D)6p 7D3 level at 2.07 eV, i.e. 0.59 eV below the lower level. The authors then correct the level energies to NIST values, but the RDW collision strength is built from eigenvectors, not only eigenvalues. If the correction is the standard a posteriori energy shift after diagonalization, the wavefunctions are not rebuilt; the poor model that produced the inversion can still feed wrong mixing coefficients into the matrix elements of Eq. (1). The paper reports A-values for the three transitions but does not validate them against NIST or MCDF/DARC values, and the only external comparison (with Kwon et al. for 400.87 nm) is qualitative: the present rates are said to be 'lower' without a quantitative factor. Thus the central claim depends on an unverified assumption that the energy-corrected target states are accurate. This is not a disagreement with consensus; it is an internal gap, because the uncorrected model is demonstrably inconsistent for one of the three target transitions.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":6268,"tokens_out":4459,"duration_ms":58317,"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":[{"comment":"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.","section":"§3, comparison with Ref. [18]"},{"comment":"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.","section":"§3, near-threshold region and Eq. (5)"}],"minor_comments":[{"comment":"Collision strength is dimensionless; the vertical axis label 'Collision Strength (a.u.)' is misleading.","section":"Fig. 1 caption"},{"comment":"Typographical errors: 'collison' (Fig. 1), 'transtions' (§3), 'addtion' (§3), 'KeV' (§3), and 'Driac-R matrix' should be 'Dirac R-matrix'.","section":"Throughout"},{"comment":"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.","section":"Fig. 4 caption"},{"comment":"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.","section":"References"},{"comment":"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.","section":"Data availability"}],"recommendation":"major_revision","confidential_remarks":"The paper is a straightforward FAC/RDW application with modest novelty, but the missing validation of the target-state wavefunctions is a legitimate reason to require major revision before publication. I do not recommend rejection because the energy-correction approach is a known practice and the required validation is feasible within the paper's scope. The authors should be asked to provide a quantitative comparison with the prior 400.87 nm calculation and to demonstrate that the corrected target states produce reliable radiative rates and stable collision strengths."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my take. The paper delivers what it says: new RDW/FAC electron impact excitation cross sections and Maxwellian rate coefficients for two WI lines (429.46, 430.21 nm) and a recalculation of 400.87 nm. That's genuinely useful for people modeling tungsten erosion in tokamaks, and the authors are transparent about the energy-level correction, even listing uncorrected FAC values in Table 2. The two new datasets are a legitimate extension of an established approach.\n\nThe soft spot is real and it's the biggest one: Table 2 shows the uncorrected FAC model puts the lower level at 2.66 eV instead of 0.365 eV, and for 430.21 nm it puts the upper level 0.59 eV below the lower. The authors then shift eigenvalue energies to NIST values, but that doesn't rebuild the eigenvectors. The RDW matrix elements in Eq. (1) depend on those eigenvectors, so if the model that produced the inversion has badly mixed states, the cross sections inherit that error. The paper does not test this. It reports A-values but doesn't compare them against NIST or the MCDF/DARC results, and the only comparison with Kwon et al. for 400.87 nm is qualitative ('lower') with no factor.\n\nThere are two smaller issues. No uncertainty estimates anywhere. And Eq. (5) is a generic low-energy correction imported from other work; its validity near threshold for these particular transitions isn't argued.\n\nNone of this makes the paper garbage. The method is standard, the calculation is what it says, and the two new rate coefficient sets could be useful even with a caveat. But the 430.21 nm data in particular rest on an untested wavefunction assumption, so the data should be labeled preliminary until validated against an independent method or experiment.\n\nMy verdict: send it to peer review, but with a referee who knows atomic structure and will demand a quantitative comparison with existing 400.87 nm data and some check on the eigenvectors (e.g., a CI expansion with more correlation, or A-value benchmarks). A published version that includes those checks would be worth citing.","headline":"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.","tokens_in":6749,"tokens_out":2358,"would_cite":false,"duration_ms":28463,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["52.20.Fs","52.55.Fa"],"model":"deepseek-v4-flash","headline":"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.","keywords":["tungsten","electron-impact excitation","relativistic distorted wave","flexible atomic code","rate coefficients","tokamak diagnostics","WI emission lines","fusion plasma"],"falsifier":"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.","tokens_in":5841,"feed_emoji":"⚛️","tokens_out":8799,"duration_ms":91670,"temperature":0.7,"pith_summary":"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.","feed_headline":"Two tungsten lines get their first excitation cross sections","feed_subtitle":"New rate coefficients for these tokamak-observed lines will sharpen erosion and impurity diagnostics in fusion devices.","key_machinery":"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","core_discovery":"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","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"The Flexible Atomic Code and its RDW implementation supply the method for all collision-strength and cross-section calculations in the paper.","marker":"[14]"},{"why":"Prior MCDF/Dirac-R-matrix rate coefficients for the 400.87 nm line that the paper extends and uses as its comparison benchmark.","marker":"[18]"},{"why":"Measured NIST energy levels that the FAC energies are shifted to match, grounding the level-correction step.","marker":"[22]"},{"why":"Provides the low-energy correction formula used to patch the perturbative RDW cross sections near threshold.","marker":"[20]"},{"why":"Earlier RDW application to tungsten that establishes the methodological lineage for the present calculation.","marker":"[8]"},{"why":"ASDEX-Upgrade erosion measurement using the 400.9 nm WI line, the diagnostic use case the new data serve.","marker":"[5]"},{"why":"Previous application of the corrected-RDW and rate-coefficient recipe to a heavy atom (cesium), the template for this calculation.","marker":"[21]"},{"why":"Spectroscopic-diagnostics framework for tungsten emissions that motivates the demand for accurate atomic data.","marker":"[7]"}],"fun_headline_variants":["First cross sections for two tungsten tokamak lines","Tungsten excitation data updated with NIST-corrected levels","New electron-impact rates for three tungsten emission lines","Accurate tungsten line data sharpens fusion diagnostics","Two neutral tungsten lines get first cross sections"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["First cross sections for two tungsten tokamak lines","Tungsten excitation data updated with NIST-corrected levels","New electron-impact rates for three tungsten emission lines","Accurate tungsten line data sharpens fusion diagnostics","Two neutral tungsten lines get first cross sections"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000704,"raw_usage":{"total_tokens":2947,"prompt_tokens":613,"completion_tokens":2334,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":357,"completion_tokens_details":{"reasoning_tokens":2259}},"tokens_in":357,"tokens_out":2334,"duration_ms":18885,"temperature":1.0,"reasoning_tokens":2259,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T13:06:05.069956+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The Flexible Atomic Code and its RDW implementation supply the method for all collision-strength and cross-section calculations in the paper."},{"cited_title":"Sang et al, Nuclear Fusion, 61, 066004 (2021)","cited_arxiv_id":null,"evidence_quote":"Prior MCDF/Dirac-R-matrix rate coefficients for the 400.87 nm line that the paper extends and uses as its comparison benchmark."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Measured NIST energy levels that the FAC energies are shifted to match, grounding the level-correction step."},{"cited_title":"Bucalossi et al, Nucl","cited_arxiv_id":null,"evidence_quote":"Provides the low-energy correction formula used to patch the perturbative RDW cross sections near threshold."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Earlier RDW application to tungsten that establishes the methodological lineage for the present calculation."},{"cited_title":"Thoma et al, Plasma Phys","cited_arxiv_id":null,"evidence_quote":"Previous application of the corrected-RDW and rate-coefficient recipe to a heavy atom (cesium), the template for this calculation."},{"cited_title":"Finally, Figure 4 illustrates the EIE rate coe fficients for the above mentioned neutral W wavelengths, simulated with and without corrected excited energy levels","cited_arxiv_id":null,"evidence_quote":"Spectroscopic-diagnostics framework for tungsten emissions that motivates the demand for accurate atomic data."}],"review_version":1}