{"id":"7d709599-e3b5-4812-8263-e0452262a814","arxiv_id":"2507.21946","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Precision resonance ionization spectroscopy in a hypersonic gas jet yields the first ionization potential of thorium at 6.306879(14) eV and the second at 12.300(9) eV, using a new Rydberg-series field-ionization technique.","lead":"Laser spectroscopy in a cold, fast gas jet measured the ionization energies of thorium atoms and thorium ions far more precisely than before, and found an efficient laser ionization route relevant to the nuclear clock candidate isotope 229Th. The improved values and new technique bring the study of the 229Th isomer in its singly charged state a step closer.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"IP2 rests on an unvalidated empirical S-curve; the Th calibration shows a 21 cm^-1 offset and no systematic uncertainty is assigned to the 73 cm^-1 error.","rationale":"The reader's verdict of CONDITIONAL is appropriate, but their primary weakest assumption focused on the IP1 Rydberg series assignment. I find that concern less decisive: the constant-quantum-defect chain in Table II has 17 members, the residuals are small, and the resulting IP1 50868.41(11) cm^-1 agrees with the literature value 50867(2) cm^-1 within 1.4 cm^-1. The IP2 determination is more fragile because it is based solely on an empirical S-curve convention with no Rydberg-series cross-check. The quoted uncertainty of 73 cm^-1 covers the observed 21 cm^-1 offset in neutral Th, but there is no evidence that the offset transfers to Th+, where the spectral congestion, collisional broadening, and autoionizing structures are very different. The paper does not assign a systematic uncertainty for the threshold method or for the subjective choice of background points, so the 73 cm^-1 error bar may understate the true uncertainty. This is the most load-bearing concern because the factor-22 improvement in IP2 is a headline result and the threshold method is its only experimental basis. The conditional acceptance recommended by the reader should therefore explicitly require a robustness analysis of the S-curve definition and background-point selection, or an independent calibration of the threshold offset for a second species in the same apparatus.","tokens_in":29895,"tokens_out":8101,"duration_ms":92883,"concrete_test":"Reanalyze the raw Th+ threshold data (Sec. III C, Fig. 4) with the E_IP2 definition varied: (a) use the sigmoid inflection point x0 alone, (b) use the half-maximum energy, (c) use a linear extrapolation of the steepest rise to the upper plateau. Independently vary the set of background points by selecting different 0.5 cm^-1 windows in the same resonance-free regions. If the spread of extracted E_IP2 exceeds the quoted 73 cm^-1, the result is model-dependent and a systematic term must be added; if the spread is small, the empirical definition is robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central IP2 claim rests entirely on the empirical S-curve threshold defined in Eq. (1): E_IP2 = 2*beta + x0, fitted to manually selected background points in a congested Th+ spectrum (Sec. III C). This definition has no first-principles justification; it is a convention calibrated on other species. The only calibration presented in this paper is the neutral Th case (Sec. V B), where the threshold fit lies 21 cm^-1 below the Rydberg-Ritz IP1. That offset is species- and spectrum-dependent: it depends on the density, widths, and autoionizing character of states near threshold, which differ substantially between Th and Th+ (the Th+ spectrum contains more than 700 peaks and is heavily affected by collisional quenching). The quoted uncertainty 99 207(73) cm^-1 appears to be the statistical fit error; no systematic uncertainty is assigned for the method or for the point selection. If the systematic offset for Th+ is comparable to or larger than 73 cm^-1, the stated value 12.300(9) eV is biased beyond its error budget and the claimed 22-fold improvement is not established. Unlike IP1, there is no Rydberg-series measurement of IP2 to anchor the method, so the threshold fit is the sole evidence.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports laser ionization spectroscopy of 232Th and 232Th+ in an argon gas cell and in a hypersonic gas jet. The authors identify a 17-member Rydberg series in neutral Th, fit it with the Rydberg-Ritz formula, and obtain IP1(Th) = 50868.41(2)_stat(11)_sys cm^-1, a large improvement over the literature value. For Th+, they extract a photoionization threshold by fitting a sigmoid to selected background regions of a congested spectrum and quote IP2(Th+) = 99207(73) cm^-1 = 12.300(9) eV. The paper also reports autoionizing states with laser ionization efficiencies up to a few percent, MCDHF and HFR+CPOL atomic structure calculations, an in-gas-jet hyperfine spectrum of 229Th, and an unsuccessful search for laser ionization of 229Th+ from 233U recoil sources.","tokens_in":30194,"tokens_out":7688,"duration_ms":87852,"significance":"If the two ionization potentials are correct, the paper delivers the most precise first ionization potential of thorium and a substantially improved second ionization potential, together with the first in-gas-jet Rydberg series detected via field ionization. The strengths of the paper are concrete: the Rydberg-Ritz analysis uses 17 unperturbed members with a constant quantum defect and a clearly stated step-size systematic uncertainty; the field-ionization mechanism is supported by time-of-flight measurements; and the gas-jet performance is demonstrated by a 238(30) MHz hyperfine line. The IP2 result is, however, less secure because it relies on an empirical threshold convention whose only in-paper validation is a single species with a 21 cm^-1 offset. The paper would be an important experimental contribution after that systematic issue is addressed quantitatively.","major_comments":[{"comment":"The quoted IP2 uncertainty of 73 cm^-1 is presented as a statistical fit error, but no systematic uncertainty is assigned for the threshold method or for the manual selection of background points. The single-species validation in Sec. V B yields a 21 cm^-1 offset between the threshold result and the Rydberg-Ritz IP1; the authors state that this is within the 73 cm^-1 statistical uncertainty, but that does not bound the offset for Th+, where the spectrum is much more congested and no independent Rydberg-series anchor exists. Because the abstract and conclusion claim a 22-fold improvement on IP2, an explicit systematic error budget and a discussion of the transferability of the calibration to Th+ are required before this claim is established.","section":"Sec. III C and Eq. (1), with Sec. V B"},{"comment":"The threshold fit depends on choices that are not reflected in the stated statistical error: the 0.5 cm^-1 averaging regions, the selection of 'valleys' at high energy, and the empirical definition E_IP2 = 2*beta + x0. I request a sensitivity analysis showing how the extracted IP2 changes when these choices are varied, and ideally a Monte Carlo over the point-selection procedure. Without such an analysis, the reader cannot distinguish the quoted 73 cm^-1 from a fitting artifact.","section":"Sec. III C, Fig. 4, Eq. (1)"},{"comment":"The 17-member Rydberg series is identified by constancy of the quantum defect, with n = 56, 57, 58 excluded as perturbed. Please provide a quantitative membership criterion: residuals of the Rydberg-Ritz fit in units of the line-position uncertainties, and the variation of the fitted IP and delta when individual members are included or excluded. The current result is plausible, but the statistical uncertainty of 0.02 cm^-1 is likely dominated by the choice of series membership, which is not quantified.","section":"Sec. V B, Table II, Fig. 9"}],"minor_comments":[{"comment":"The sentence 'These investigations, led to the extraction' should read 'These investigations led to the extraction'.","section":"Sec. I"},{"comment":"The phrase 'The raise in count rate' should be 'The rise in count rate'.","section":"Sec. III C"},{"comment":"In the row beginning '51 50825.28(1)', the leading '51' appears to be a formatting artifact; the row should be labeled as n = 51.","section":"Table II"},{"comment":"The phrase 'an 22-fold improvement' should read 'a 22-fold improvement', and 'in forseseen cases' should be 'in foreseen cases'.","section":"Sec. VII"},{"comment":"The reference list contains several formatting problems, notably Ref. [29] ('A. Kramida, Yu. Ralchenko, J. Reader, and and NIST ASD Team') and Ref. [39], which has an embedded URL in the title; these should be cleaned.","section":"References"},{"comment":"The phrase 'The found IP2' should be 'The fitted IP2' or 'The determined IP2'.","section":"Fig. 4 caption"}],"recommendation":"major_revision","confidential_remarks":"The paper is a strong experimental contribution and is within the journal's scope. My main reservation is the IP2 systematic error budget; the IP1 Rydberg-Ritz result appears sound. I did not see a circularity problem in the IP1 determination: using the literature IP as an initial guess for the principal quantum number assignment is absorbed by the allowed integer offset in the quantum defect and does not bias the fitted limit. The in-gas-jet Rydberg spectroscopy and the demonstration of field ionization are notable achievements. The revision should focus on making the IP2 error budget credible, rather than on redoing the central IP1 analysis."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is a genuinely useful experimental contribution. The standout result is the first Rydberg series measured in a hypersonic gas jet with field ionization, and it delivers a clean, nine-fold improved IP1 for thorium: 50868.41(11) cm^-1. The Rydberg-Ritz fit over 17 members, with three excluded as perturbed, is standard and well documented; the step-size systematic dominates honestly. The MCDHF and HFR+CPOL calculations are competent and add value, and the new Th+ ionization scheme with measured efficiencies is practically useful.\n\nThe soft spot is IP2. The value 99207(73) cm^-1 comes from fitting an empirical sigmoid to manually selected background points in a congested spectrum. The only calibration shown is the neutral Th measurement, which lands 21 cm^-1 below the Rydberg-Ritz IP1. The paper notes this is within the 73 cm^-1 statistical uncertainty, which is true, but no systematic uncertainty is assigned for the method or the point selection. For Th+, where the spectrum is far more congested and autoionizing, the offset could plausibly be larger. The 22-fold improvement claim is therefore not as secure as the IP1 claim. I do not think this is fatal: even a 50 cm^-1 systematic would still leave the value far better than the previous 11.9–12.3 eV bracket. But the authors should either add a systematic term or quote a more conservative error.\n\nTwo minor issues: the factor-of-two inconsistency in the 229Th+ count rate is left unresolved, and the fit without the three excluded Rydberg members is not shown. Neither undermines the main results.\n\nOverall, this is a solid paper that deserves a serious referee. IP1 is publishable essentially as is; IP2 needs an honest systematic-error discussion before the precision claim is fully credible. I would send it to review with that condition.","headline":"Clean, significantly improved IP1 from a first in-gas-jet Rydberg series; IP2 is plausible but the threshold method's missing systematic error keeps it from being equally secure.","tokens_in":30816,"tokens_out":2901,"would_cite":true,"duration_ms":34141,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["32.10.-f","32.30.-r","32.80.Fb","32.80.Rm"],"model":"deepseek-v4-flash","headline":"Thorium's first two ionization potentials are pinned down to record precision","keywords":["thorium ionization potential","Rydberg series","hypersonic gas jet","laser ionization spectroscopy","autoionizing states","field ionization","nuclear clock isomer","MCDHF calculations"],"falsifier":"A future high-resolution measurement of a different Rydberg series of neutral Th, or of the Th+ Rydberg series in a gas jet, that fits a series limit incompatible with 50 868.41(11) $cm^{-1}$ beyond the combined uncertainties would falsify the IP1 claim; for IP2, a direct Rydberg-series limit for Th+ differing from 99 207(73) $cm^{-1}$ by more than the 21 $cm^{-1}$ calibration offset would falsify the threshold extrapolation.","tokens_in":29697,"feed_emoji":"⚛️","tokens_out":6881,"duration_ms":66804,"temperature":0.7,"pith_summary":"This paper reports the most precise values yet for the first and second ionization potentials of thorium: IP1 = 50 868.41(11) $cm^{-1}$ and IP2 = 99 207(73) $cm^{-1}$, improving on previous measurements by factors of 9 and 22. The key advance is performing resonance ionization spectroscopy in a hypersonic argon gas jet, where collision-induced quenching is suppressed and a clean Rydberg series converging to IP1 becomes visible. The same gas-cell data yield an efficient laser ionization scheme for Th+ based on an autoionizing state at 99 766.87 $cm^{-1}$, intended for future study of the 229mTh nuclear-clock isomer. The authors also validate their threshold-fitting method for IP2 by showing it reproduces the Rydberg-derived IP1 with a 21 $cm^{-1}$ offset.","feed_headline":"Thorium's ionization energies pinned with record precision","feed_subtitle":"First IP 50 868.41(11) cm^-1; second IP 99 207(73) cm^-1, from in-gas-jet Rydberg fits.","key_machinery":"The central objects are the Rydberg series of neutral thorium observed in the hypersonic jet, described by the Rydberg-Ritz formula $E_n = E_{\\mathrm{IP}} - R_\\mu/(n-\\delta)^2$ with a constant quantum defect $\\delta = 0.56(1)$, and the empirical S-curve threshold $E_{\\mathrm{IP2}} = 2\\beta + x_0$ applied to the ion. The hypersonic argon jet, with Mach number above 8, is the enabling mechanism: it reduces the collision rate by about four orders of magnitude relative to the gas cell, suppressing quenching so that the Rydberg series becomes resolvable, and the electric field of the RFQ ion guide subsequently field-ionizes the Rydberg atoms in a localized region, as confirmed by time-of-flight measurements.","core_discovery":"We determine the first ionization potential of thorium to be 50 868.41(2)_stat(11)_sys $cm^{-1}$, or 6.306 879(14) eV, from a Rydberg-Ritz fit to 17 unperturbed members of a single Rydberg series observed in a hypersonic gas jet, and the second ionization potential to be 99 207(73) $cm^{-1}$, or 12.300(9) eV, from an S-curve threshold fit to photoionization data in a gas cell. The gas jet suppresses collisional quenching and delivers a spectral resolution of 240(30) MHz, and time-of-flight measurements show that the Rydberg atoms are ionized by field ionization in the RFQ region. Multiconfigurational Dirac-Hartree-Fock calculations reproduce the experimental values within 0.06% (IP1) and 0.19% (IP2). A laser ionization scheme for Th+ based on the autoionizing state at 99 766.87(22) $cm^{-1}$ reaches an efficiency of at least 1.2%, but attempts to observe laser photoionization of 229Th+ from a 233U recoil source were unsuccessful.","pith_inferences":["If the in-gas-jet Rydberg method can be extended to Th+, the IP2 uncertainty, currently 73 cm^-1, could shrink by roughly two orders of magnitude to match the IP1 precision; the paper's inability to observe a Th+ Rydberg series in the gas jet makes that the direct next test.","The unsuccessful 229Th+ photoionization suggests that a substantial fraction of recoiling ions occupy long-lived metastable 'dark' states; if so, isomer studies may need either an ionization scheme targeting those states or a neutralization-reionization strategy.","The constant quantum defect $\\delta = 0.56(1)$ over the range n = 44 to 63 is consistent with a single 6d7s np Rydberg series; identifying its spectator configuration could tie the measured limit to a specific LS term and test the present level assignments.","If the threshold S-curve offset of 21 cm^-1 is found to be similar for other elements, the threshold method could serve as a semi-empirical predictor of ionization potentials for species where Rydberg series are experimentally inaccessible."],"forward_implications":["The first ionization potential of thorium is now known to 0.11 cm^-1, a 9-fold improvement, providing a precise anchor for Rydberg extrapolations and atomic-structure calculations across the actinides.","The second ionization potential at 99 207(73) cm^-1 constrains the energies of autoionizing states and the electronic-bridge decay estimates relevant to the 229mTh isomer in Th+.","The demonstrated combination of in-gas-jet Rydberg spectroscopy with subsequent field ionization establishes a route to high-resolution ionization-potential measurements for short-lived actinides at on-line laser-spectroscopy facilities.","The laser ionization scheme using the 99 766.87 cm^-1 autoionizing state offers a practical basis for efficient gas-cell laser ion sources producing thorium beams.","The observed 21 cm^-1 offset between the threshold method and the Rydberg-Ritz fit provides a species-specific calibration point for the threshold approach."],"supporting_citations":[{"why":"Supplies the previous best IP1 value of 50 867(2) cm^-1, the baseline that the Rydberg-Ritz result improves by a factor 9.","marker":"[35]"},{"why":"Provides the earlier energy-level survey of Th+ and the 11.9-12.3 eV range for IP2 that motivates the threshold determination.","marker":"[14]"},{"why":"Characterizes the hypersonic de Laval nozzle and the 17 K jet conditions used to obtain the high-resolution Rydberg spectrum.","marker":"[22]"},{"why":"Demonstrates the in-gas-jet spectroscopy method on 254No, proving the technique that is here applied to thorium.","marker":"[23]"},{"why":"Supplies the CODATA reduced Rydberg constant $R_\\mu$ used in the Rydberg-Ritz series fit.","marker":"[39]"},{"why":"Introduces the S-curve threshold analysis for the ionization potential of astatine, the method adapted for IP2.","marker":"[30]"},{"why":"Documents the systematic deviations of the threshold method for several species, used to justify the 21 cm^-1 offset seen for thorium.","marker":"[31]"},{"why":"Provides earlier MCDHF/pseudopotential IP1 and IP2 values (50 813 cm^-1 and 100 093 cm^-1) that the present calculations improve upon.","marker":"[50]"}],"fun_headline_variants":["Thorium gas-jet spectroscopy pins both ionization potentials","Hypersonic jet yields precise thorium ionization energies","Thorium IPs measured via hypersonic jet: 6.306879 eV and 12.300 eV"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The IP1 value assumes the 17 selected peaks form a single Rydberg series with a constant quantum defect $\\delta = 0.56(1)$ and that the three excluded members (n = 56, 57, 58) are genuinely perturbed; the IP2 value additionally assumes the empirical S-curve threshold offset, calibrated on neutral thorium with a 21 $cm^{-1}$ deviation, applies to Th+.","fun_headline_variants_meta":{"raw":{"variants":["Thorium gas-jet spectroscopy pins both ionization potentials","Hypersonic jet yields precise thorium ionization energies","Thorium IPs measured via hypersonic jet: 6.306879 eV and 12.300 eV"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001654,"raw_usage":{"total_tokens":6670,"prompt_tokens":1148,"completion_tokens":5522,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":764,"completion_tokens_details":{"reasoning_tokens":5459}},"tokens_in":764,"tokens_out":5522,"duration_ms":47602,"temperature":1.0,"reasoning_tokens":5459,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T12:12:50.975793+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A future high-resolution measurement of a different Rydberg series of neutral Th, or of the Th+ Rydberg series in a gas jet, that fits a series limit incompatible with 50 868.41(11) $cm^{-1}$ beyond the combined uncertainties would falsify the IP1 claim; for IP2, a direct Rydberg-series limit for Th+ differing from 99 207(73) $cm^{-1}$ by more than the 21 $cm^{-1}$ calibration offset would falsify the threshold extrapolation.","supporting_citations":[{"cited_title":"K¨ ohler, R","cited_arxiv_id":null,"evidence_quote":"Supplies the previous best IP1 value of 50 867(2) cm^-1, the baseline that the Rydberg-Ritz result improves by a factor 9."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the earlier energy-level survey of Th+ and the 11.9-12.3 eV range for IP2 that motivates the threshold determination."},{"cited_title":"Ferrer, M","cited_arxiv_id":null,"evidence_quote":"Characterizes the hypersonic de Laval nozzle and the 17 K jet conditions used to obtain the high-resolution Rydberg spectrum."},{"cited_title":"Lantis, A","cited_arxiv_id":null,"evidence_quote":"Demonstrates the in-gas-jet spectroscopy method on 254No, proving the technique that is here applied to thorium."},{"cited_title":"Rothe, A","cited_arxiv_id":null,"evidence_quote":"Introduces the S-curve threshold analysis for the ionization potential of astatine, the method adapted for IP2."},{"cited_title":"Rothe, An all-solid state laser system for the laser ion source RILIS and in-source laser spectroscopy of astatine at ISOLDE/CERN , Ph.D","cited_arxiv_id":null,"evidence_quote":"Documents the systematic deviations of the threshold method for several species, used to justify the 21 cm^-1 offset seen for thorium."},{"cited_title":"Froese Fischer, G","cited_arxiv_id":null,"evidence_quote":"Provides earlier MCDHF/pseudopotential IP1 and IP2 values (50 813 cm^-1 and 100 093 cm^-1) that the present calculations improve upon."}],"review_version":1}