{"id":"d8d9559c-af07-4657-bd05-525b1be18a9a","arxiv_id":"2507.21410","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A systematic relativistic Hartree-Fock calculation produces tabulated field and mass isotope shift coefficients for total electron binding energy of atoms and ions up to Z=120, with a simple power-law interpolation accurate to about 1%.","lead":"This paper computes how the total energy that binds electrons in an atom changes between isotopes, for neutral atoms and singly charged ions up to element 120. It provides tables of field and mass isotope shift coefficients that are useful for precise nuclear mass measurements and for searches for new physics using isotope shifts.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Percent-level field-shift accuracy rests on unquantified neglect of QED and correlation corrections to deep-shell contributions; a concrete estimate is needed.","rationale":"The reader's weakest_assumption correctly identifies the unsupported assertion that correlation corrections to the total-energy field shift are negligible. My independent analysis finds the same point is the most load-bearing: the paper's percent-level accuracy claim for F and G rests entirely on deep-shell dominance, and no quantitative estimate is provided for the corrections to the field shift itself. I sharpen the concern by noting that QED corrections, which are largest for deep s shells, have finite-nuclear-size derivatives that have not been estimated at all; without such an estimate the 'expected percent-level' uncertainty is an act of faith. The interpolation formula, while only weakly benchmarked, is a separate and less fundamental issue because it is an internal consistency check that could be repaired with a residual table. My recommendation of UNCHANGED reflects that the reader's CONDITIONAL verdict already captures exactly this weakness; no verdict adjustment is needed, but the conditionality is reinforced by the concrete QED test proposed.","tokens_in":7165,"tokens_out":8628,"duration_ms":111015,"concrete_test":"For one heavy closed-shell atom, e.g., Z=102 (No) or Z=118 (Og), recompute the field-shift coefficient F using the same RHF+Breit Hamiltonian but with model QED potentials for self-energy and vacuum polarization (e.g., the Flambaum-Ginges radiative potentials) added, using the same nuclear charge distributions. Compare the resulting F with the corresponding value in Table III. If the difference exceeds about 1%, the deep-shell justification is quantitatively insufficient and the uncertainty statement must be revised. Repeating at Z=54 (Xe) would also bound the effect in the interpolation range.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that the tabulated F and G coefficients, together with the bZ^k interpolation, describe the field shift of the total electron binding energy to about 1% for heavy atoms. The interpolation part is an empirical fit and could be validated by the paper's open-shell benchmark cases, but no residual comparison is shown. The more fundamental part is the assertion in Section III: 'Isotope shift in the total electron energy is dominated by the contribution of deep shells, therefore correlation corrections are expected to be small.' This assertion is the only justification for the percent-level accuracy of F and G themselves, and it is not quantified. The deepest shells are also where QED effects are largest: for Z~100 the 1s self-energy and vacuum polarization shift binding energies by hundreds of eV, and their finite-nuclear-size derivatives enter the field shift directly. The paper does not estimate whether these derivatives, or the derivatives of the correlation energy, are below 1% of F (which is a few eV/fm^2 for heavy atoms). If they are not, the tabulated coefficients carry a systematic error of unknown sign that propagates into every application. This is a load-bearing correctness risk, not a style issue.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper computes the isotope shift of the total electron binding energy for neutral atoms and selected ions up to Z=120 using relativistic Hartree-Fock including the Breit interaction. Field-shift coefficients F and G are extracted from total-energy calculations at three nuclear radii and tabulated for closed-shell systems, together with a coefficient a for an A^(1/3) parametrization. The normal mass shift is estimated from the total kinetic energy via the virial theorem, the specific mass shift is evaluated at the RHF level but described as unimportant, and a power-law interpolation F(Z)=bZ^k is proposed with a claimed accuracy of about 1%. The paper concludes that the field shift dominates above Z≈38 and that differences between neutral atoms and singly charged ions are at the few-percent level.","tokens_in":7357,"tokens_out":7476,"duration_ms":95131,"significance":"The paper fills a genuine gap: no comprehensive tabulation of isotope shifts of total electron binding energies currently exists, and the results would be useful for nuclear mass evaluations, drip-line studies, and isotope-shift searches for new physics. The method is transparent, the three-point fit retains a quadratic term for superheavy elements, and the bZ^k interpolation is simple and potentially useful. However, the central percent-level accuracy claim is not supported by the evidence actually presented: no residual analysis for the interpolation is shown, no numerical estimate is given for omitted correlation and QED contributions to the field shift, and at least one tabulated configuration appears not to be the neutral ground state. The approach is plausible, but the validation depth is insufficient for the stated accuracy.","major_comments":[{"comment":"The percent-level uncertainty claim for F and G rests entirely on the sentence \"Isotope shift in the total electron energy is dominated by the contribution of deep shells, therefore correlation corrections are expected to be small.\" This is load-bearing because F and G are extracted from RHF total energies that omit correlation and QED. For heavy atoms the 1s self-energy and vacuum polarization are hundreds of eV, and their finite-nuclear-size dependence enters the field shift alongside the RHF contribution; no estimate is given of whether the derivatives of these corrections are small compared with F, which is a few eV/fm^2 near Z~100. Please provide a quantitative estimate, for example by evaluating the finite-size derivative of the QED corrections from Ref. [19] for representative nuclei such as Pb, No, and Og, and state explicitly whether the tabulated F and G include any QED or correlation contribution. Without this, the claimed percent-level accuracy is unsupported.","section":"Section III (Conclusion)"},{"comment":"The claim that bZ^k reproduces calculated field shifts \"to within about 1%\" is not evidenced: no residual table, plot, or numerical comparison is shown, despite the text saying that open-shell benchmark cases were used to test the formula. Please include (F_fit - F_calc)/F_calc for all tabulated points and for the open-shell benchmark cases, and define precisely which neighboring points are used for each interpolation interval. This is needed because the 1% figure is the paper's main quantitative deliverable and is not currently verifiable from the manuscript.","section":"Section III and Eq. (5)"},{"comment":"The treatment of the specific mass shift is internally inconsistent for light elements. The text states that the SMS is 10-20% of the NMS in Xe and that it is not important because the mass shift is only relevant for light atoms, but Table II shows that for Z<38 the mass shift dominates the field shift. For a light atom such as Ne or Mg, neglecting a 10-20% SMS correction means the total isotope shift is wrong by 10-20%, which is not a modest error if the quantity of interest is the total binding-energy difference. Either tabulate KSMS values or explicitly restrict the paper's final formula to the field-shift part and remove the mass-shift term from Eq. (1) for light atoms.","section":"Section II, Eqs. (1)-(2) and Table II"},{"comment":"The listed ground-state configuration for Pt, [Xe]4f14 5d10, is not the neutral Pt ground state, which is [Xe]4f14 5d9 6s1. The listed configuration has 78 electrons but an empty 6s shell, and it is an excited configuration, not the ground state. If the calculation used this d10 configuration, the resulting F and G do not correspond to neutral Pt, and the few-percent effect of removing an outer s electron makes this matter for the table and for the interpolation around Z=78. If the configuration is a typographical error, please correct it and verify the reported values; otherwise, clarify why this configuration was used.","section":"Table III, Z=78 (Pt)"}],"minor_comments":[{"comment":"The abstract claims calculations for \"singly charged ions up to element Z=120,\" but Table IV lists only ten ions (In+, Cs+, La+, Lu+, Au+, Tl+, Fr+, Lr+, Nh+, E119+), not the full set of singly charged ions. Please either extend the table or revise the abstract and introduction to say \"selected singly charged ions.\"","section":"Abstract and Table IV"},{"comment":"The notation such as 4.86[-6] is not defined; please add a note in the Table I caption that [n] denotes multiplication by 10^n.","section":"Table I caption"},{"comment":"Equation (1) is split across two lines without a single equation number or a clear alignment; please reformat so that the two lines are visually one equation or are numbered separately.","section":"Eq. (1)"},{"comment":"The paper states that the SMS is evaluated but does not report any SMS values or a table of KSMS; since the mass shift is included in Eq. (1), please either provide the SMS values used or state clearly that Eq. (1) is given for completeness and that the numerical results do not include the SMS.","section":"Section II, mass shift discussion"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of an atomic-physics data journal and draws appropriately on the authors' prior work. The main risk is the unsupported percent-level accuracy claim; I would send the manuscript back for major revision with the request for quantitative estimates of omitted QED/correlation contributions and a residual analysis of the interpolation. I do not see a novelty-disclosure problem."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a useful data paper, the first systematic tabulation of isotope shift coefficients for the total electron binding energy up to Z=120. If you work on nuclear mass accounting, superheavy elements, or isotope shift applications, the tables of F, G, and a are worth having. The interpolation formula F = bZ^k with a Z-dependent exponent is simple and seems to work at the 1% level for reproducing their RHF numbers. The crossover between mass and field shift near Z=38 is cleanly presented.\n\nThe paper computes at the relativistic Hartree-Fock level with the Breit interaction, which is fine for a first pass. The quantity is genuinely new; they cite no prior calculation of it. Keeping the quadratic G term for large radius changes is sensible, and the A^(1/3) parametrization for cases where radius changes are unknown is a nice addition.\n\nThe soft spots are about what is not shown. The 1% interpolation accuracy is asserted rather than demonstrated: there is no residual table for the open-shell benchmarks, and those open-shell cases are exactly what the formula is supposed to predict. That is easy to fix. More important, the claim that the tabulated F and G are accurate at the percent level for heavy atoms rests on the assertion that deep shells dominate and correlation/QED corrections are small. They do not quantify either. For Z~100 the 1s QED corrections are hundreds of eV; the finite-nuclear-size derivative of those corrections can feed directly into F, and the reader gets no estimate of whether that is below 1% of F (which is a few eV/fm^2 for heavy atoms). Same for correlation derivatives. This is a load-bearing gap in the error budget, not a style issue. It does not kill the RHF-level tables, but it means the percent-level physical accuracy claim is unsubstantiated. A concrete estimate, even a rough one, would change the paper's value.\n\nThere is also a small inconsistency in the mass shift treatment. They say SMS is 10-20% of NMS, then neglect it in Table II. For light atoms, where mass shift dominates, 10-20% of the dominant term is not negligible. The heavy-atom conclusions are unaffected, but the wording should be tightened.\n\nBottom line: worth a serious referee. The dataset is new and useful; the requested changes are quantitative error estimates, a residual table for the interpolation, and a cleaned-up mass shift discussion. I would engage with it.","headline":"Useful first tabulation of total-binding-energy isotope shift coefficients up to Z=120, but the percent-level accuracy claim for the field shift factors needs quantitative support.","tokens_in":7884,"tokens_out":4027,"would_cite":true,"duration_ms":45661,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper establishes a systematic tabulation of isotope shifts of the total electron binding energy for neutral atoms and singly charged ions up to Z=120, with field-shift coefficients accurate to about 1% for heavy atoms, and shows…","keywords":["isotope shift","total electron binding energy","field shift","mass shift","superheavy elements","relativistic Hartree-Fock","nuclear radius"],"falsifier":"A direct test would be to recalculate the isotope shift of the total electron binding energy for a heavy closed-shell atom such as Hg or No using a high-precision relativistic many-body method that includes electron correlation and QED corrections, and compare the resulting F and G to the tabulated values; a deviation clearly exceeding the claimed percent-level accuracy would invalidate the assumption that inner-shell dominance makes such corrections negligible.","tokens_in":6941,"feed_emoji":"⚛️","tokens_out":7144,"duration_ms":71222,"temperature":0.7,"pith_summary":"This paper provides the first systematic tabulation of isotope shifts of the total electron binding energy for all neutral atoms and singly charged ions up to Z=120, computed with relativistic Hartree-Fock including the Breit interaction. The central result is that the field-shift part of the isotope shift is accurately captured by two tabulated coefficients F and G in the formula $F\\delta\\langle r^2\\rangle + G\\delta\\langle r^2\\rangle^2$, together with a simple power-law interpolation $F=bZ^k$ that reproduces the calculated values within about 1%. This matters because isotope mass differences enter nuclear mass measurements, searches for new forces via isotope shift spectroscopy, and predictions of superheavy element stability, where subtracting the electron binding energy at percent accuracy is needed. The paper also shows that the mass shift dominates below $Z\\approx 38$ and the field shift above, and that the difference between neutral atoms and singly charged ions stays at the few percent level, so the tables extend to higher charge states.","feed_headline":"Atomic binding-energy isotope shifts tabulated to Z=120","feed_subtitle":"A simple power law gives about 1% accuracy in field shifts, aiding nuclear mass and superheavy searches.","key_machinery":"The central object is the parametrization of the isotope shift as a field shift plus mass shift, with field-shift coefficients F and G extracted by computing total binding energies at three nuclear radii and fitting a parabola. The power-law interpolation $F(Z)=bZ^k$, with $k = \\ln(F_1/F_2)/\\ln(Z_1/Z_2)$, is the device that extends the closed-shell results to open-shell atoms with about 1% accuracy. The normal mass shift uses the virial theorem $E_k=-E_{total}$ from earlier total-energy calculations, and the specific mass shift $\\langle \\sum_{i<j} \\mathbf{p}_i\\cdot \\mathbf{p}_j\\rangle$ is evaluated in the relativistic Hartree-Fock ground state, relying on the fact that inner shells dominate so correlation corrections are small.","core_discovery":"The paper establishes that the isotope shift $\\Delta E_{IS}$ of the total electron binding energy can be parameterized as $\\Delta E_{IS} = F\\delta\\langle r^2\\rangle + G\\delta\\langle r^2\\rangle^2 + (m_e/M_1 - m_e/M_2)(K_{NMS} + K_{SMS})$, with tabulated F and G for closed-shell systems from Ne to Og and benchmark open-shell cases, and that the dependence of F on Z follows $F=bZ^k$ to within about 1% between neighboring closed shells. The effective exponent k grows from about 5 near $Z\\approx 50$ to 12.6 at $Z=118\\text{--}120$, reflecting the growing contribution of deep s and p_{1/2} shells. The calculations are done with relativistic Hartree-Fock including the Breit interaction, varying the nuclear radius to extract the field-shift coefficients; the normal mass shift is obtained from the virial theorem, while the specific mass shift is evaluated as a Hartree-Fock expectation value. The result makes the electron-binding contribution to isotope mass differences available for essentially the whole periodic table.","pith_inferences":["Because the paper does not quantify correlation or QED corrections, a natural next step is to recompute F and G for a few heavy atoms with a many-body method that includes those corrections; the comparison would either confirm the percent-level estimate or bound the omitted terms.","The steep rise of the exponent k from about 5 to 12.6 could serve as a compact diagnostic for how strongly the field shift concentrates into the innermost s and p_{1/2} orbitals, and might be used to extrapolate field-shift coefficients for elements beyond Z=120.","The near-independence of the isotope shift on ionization state suggests that the same tables could anchor total-energy isotope shifts for trace elements in astrophysical plasmas, where ionization levels are not singly charged, though explicit high-charge calculations would be needed to verify this extension."],"forward_implications":["For any isotope pair, the electron-binding contribution to the atomic mass difference can now be estimated to about 1% accuracy in heavy atoms using the tabulated F and G and the power-law interpolation, with no additional atomic-structure calculation.","Nuclear mass evaluations for elements beyond $Z\\approx 38$ will need to treat the field-shift part of the electron-binding correction as a leading term, since the mass-shift contribution becomes small there.","The tabulated coefficients extend directly to singly charged ions and, with few-percent adjustments, to higher charge states, because the isotope shift is dominated by inner shells.","The alternative coefficient a parametrizes the field shift in terms of $A^{1/3}$, providing a way to estimate isotope shifts for superheavy isotopes when only mass numbers are known, not radii.","The crossover near $Z\\approx 38$ identifies where the dominant uncertainty in isotope mass differences shifts from atomic electron motion to nuclear charge radii."],"supporting_citations":[{"why":"Supplies the total electron binding energies used to obtain kinetic energies via the virial theorem and sets the baseline for total-energy isotope shifts.","marker":"[19]"},{"why":"Provides the standard decomposition into field shift and mass shift that the paper parameterizes and tabulates.","marker":"[20]"},{"why":"Gives the $\\delta\\langle r^{2\\gamma}\\rangle$ representation of field shifts and the deep-shell dependence used to justify the power-law interpolation.","marker":"[21]"},{"why":"Supplies the $A^{1/3}$ parametrization of the field shift used for superheavy elements when radius changes are not known.","marker":"[24]"},{"why":"Reports the only prior calculation of total-energy isotope shift for several heavy atoms, which the present work generalizes.","marker":"[6]"},{"why":"Documents the sub-eV mass-difference precision needed in Penning-trap measurements, motivating the need for accurate electron-binding corrections.","marker":"[2]"}],"fun_headline_variants":["Isotope shifts of electron binding energy tabulated to Z=120","Atomic field shifts follow simple power law from Ne to Og","Binding-energy isotope shifts: 1% accuracy from a power law","Field shift coefficients for all closed-shell atoms to Og"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that correlation and QED corrections to the isotope shift are negligible because deep inner-shell electrons dominate the total-energy shift; the paper asserts this but provides no quantitative estimate, so if these corrections are not small the tabulated coefficients carry an unquantified systematic error.","fun_headline_variants_meta":{"raw":{"variants":["Isotope shifts of electron binding energy tabulated to Z=120","Atomic field shifts follow simple power law from Ne to Og","Binding-energy isotope shifts: 1% accuracy from a power law","Field shift coefficients for all closed-shell atoms to Og"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000678,"raw_usage":{"total_tokens":3088,"prompt_tokens":958,"completion_tokens":2130,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":574,"completion_tokens_details":{"reasoning_tokens":2059}},"tokens_in":574,"tokens_out":2130,"duration_ms":19022,"temperature":1.0,"reasoning_tokens":2059,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T12:46:36.310964+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct test would be to recalculate the isotope shift of the total electron binding energy for a heavy closed-shell atom such as Hg or No using a high-precision relativistic many-body method that includes electron correlation and QED corrections, and compare the resulting F and G to the tabulated values; a deviation clearly exceeding the claimed percent-level accuracy would invalidate the assumption that inner-shell dominance makes such corrections negligible.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the total electron binding energies used to obtain kinetic energies via the virial theorem and sets the baseline for total-energy isotope shifts."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the standard decomposition into field shift and mass shift that the paper parameterizes and tabulates."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the $\\delta\\langle r^{2\\gamma}\\rangle$ representation of field shifts and the deep-shell dependence used to justify the power-law interpolation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the $A^{1/3}$ parametrization of the field shift used for superheavy elements when radius changes are not known."},{"cited_title":"Further steps towards next generation of covariant energy density functionals","cited_arxiv_id":"2507.17082","evidence_quote":"Reports the only prior calculation of total-energy isotope shift for several heavy atoms, which the present work generalizes."},{"cited_title":"Dilling, K","cited_arxiv_id":null,"evidence_quote":"Documents the sub-eV mass-difference precision needed in Penning-trap measurements, motivating the need for accurate electron-binding corrections."}],"review_version":1}