{"id":"22e64a8e-16fe-4857-893c-30336e11aac9","arxiv_id":"2411.08213","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"QED corrections to ionization and transition energies are implemented in a two-component ZORA framework and shown to mostly match four-component reference calculations within a few percent.","lead":"This paper implements quantum electrodynamic corrections (vacuum polarization and electron self-energy) in a two-component relativistic quantum chemistry framework called ZORA, and benchmarks them against four-component calculations for atoms and molecules. The work matters because it provides an efficient way to estimate QED effects in heavy-element molecules like radium monofluoride, which are used in precision searches for new physics.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Unexplained 63.8% deviation for Au 5p1/2 in Table V contradicts the unqualified 'excellent agreement' claim and leaves the ZORA picture-change transformation unvalidated for p1/2 orbitals.","rationale":"The reader identified the same weakest assumption (ZORA picture-change transformation), and Table V is the right place to test it. I considered but did not elevate the cubic spline interpolation (reported mean errors <=6.7e-5) or the point-like nucleus choice for QED potentials, because neither explains a 63.8% orbital-specific deviation and neither affects the abstract's blanket claim as directly. The good agreement in Tables I, III, and IV for ionization and transition energies is genuine evidence that the implementation is mostly sound; this is why the appropriate verdict remains CONDITIONAL rather than REJECT. The request to soften the abstract and explain Table V is a reasonable condition, not a wholesale invalidation. My concern is essentially the reader's, so the verdict should be unchanged.","tokens_in":21480,"tokens_out":7039,"duration_ms":72867,"concrete_test":"Reproduce the Au B3LYP/dyall.3zp self-consistent FG+UE results in Table V, but also report the Flambaum-Ginges and Uehling contributions separately for each orbital. Then recompute the same orbital corrections with an exact two-component (X2C) transformation or with a four-component Dirac-Kohn-Sham code using identical QED potentials, basis, and nuclear model. If the 5p1/2 relative deviation collapses when the separate contributions are compared or when X2C is used, the 63.8% is a near-cancellation artifact or a ZORA picture-change defect; if it persists in the four-component comparison, the disagreement originates from the reference protocol rather than from the two-component approximation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim—that QED corrections can be obtained in a two-component ZORA framework in excellent agreement with four-component results—rests on the picture-change transformation of Ref. [60] (Eqs. 13-14) producing reliable matrix elements for all orbitals of interest. This condition is not met in the authors' own benchmark data. In Table V, self-consistent ZORA-B3LYP/dyall.3zp FG+UE orbital-energy corrections for Au deviate from four-component B3LYP/dyall.3zp reference by 29.6% (1s1/2), 15.3% (2p1/2), and 63.8% (5p1/2), with the text stating that the origin of the latter deviation is 'presently open'. These are not deep-core-only artifacts: the 6s1/2 HOMO correction deviates by 10.0-10.3%. The large relative error for 5p1/2 occurs on a small absolute correction, so cancellation between the Flambaum-Ginges and Uehling contributions may amplify the relative deviation, but the paper does not provide the separate contributions needed to distinguish a cancellation artifact from a genuine failure of the ZORA picture-change transformation for p1/2 orbitals. Without that analysis, the unqualified abstract claim overstates the evidence and the generality of the method for valence and core-level spectroscopy remains unestablished.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript implements vacuum-polarisation (Uehling) and electron-self-energy (Flambaum-Ginges, Pyykkö-Zhao) effective potentials as one-electron operators in a two-component ZORA framework, with picture-change corrections following Ref. [60]. It benchmarks QED contributions to atomic ionization energies, ionic transition energies, valence orbital energies, and selected molecular properties (group 2 monofluorides, BaF, RaF) against four-component reference data, and it compares perturbative and self-consistent treatments of the QED potentials for gold Kohn-Sham orbital energies. The central claim, stated in the abstract and conclusion, is that QED corrections can be obtained in this two-component framework efficiently and in excellent agreement with corresponding four-component results.","tokens_in":21830,"tokens_out":5522,"duration_ms":56097,"significance":"If the claimed agreement holds for the properties of interest, the approach would enable routine QED estimates for heavy-element molecular spectroscopy at two-component cost, which is valuable for precision studies such as those on RaF. The paper has concrete strengths: it benchmarks against four-component implementations of the same effective potentials (Refs. [29,38]), reports systematic data over a wide Z range, and does not fit parameters to the target data. However, the abstract's unqualified 'excellent agreement' claim is undermined by the authors' own Table V, which shows deviations up to 32.4% for the 1s1/2 orbital and 63.8% for the 5p1/2 orbital of gold, with the origin of the latter explicitly left open. The method's validity for core and some valence orbital energies is therefore not established by the present evidence, even though the valence and transition-energy benchmarks are largely good.","major_comments":[{"comment":"The abstract's blanket statement that QED corrections are obtained 'in excellent agreement with corresponding four-component results' is contradicted by Table V, which reports self-consistent FG+UE deviations of 29.6% (1s1/2), 15.3% (2p1/2), and 63.8% (5p1/2) for gold relative to four-component B3LYP/dyall.3zp calculations of Ref. [38] using the same functional and basis set. The text states that the origin of the 5p1/2 deviation is 'presently open.' Because these are direct benchmarks of the ZORA picture-change transformation, the central claim must be qualified. Please provide a decomposition of the 5p1/2 correction into the individual Flambaum-Ginges contributions (magnetic, high-frequency, low-frequency) and the Uehling term, with signs, to determine whether the large relative error is a cancellation artifact; also report the un-averaged Kramers-pair values for the entries affected by footnote a of Table V. The 10.0-10.3% deviation for the 6s1/2 valence orbital should also be addressed, as it is not a deep-core effect.","section":"Abstract; §IV.C, Table V"},{"comment":"The computational details state that the QED potentials are evaluated with a point-like nucleus after the SCF was performed with finite Gaussian nuclear charge distributions. This ad hoc combination is not justified in the manuscript. Since the largest deviations in Table V occur for core orbitals, where finite-nuclear-size effects are largest, please quantify the effect of using the finite nuclear charge distribution in the Uehling and Flambaum-Ginges potentials for gold, at least for the 1s1/2 and 5p1/2 orbital-energy corrections, or provide a reference demonstrating that this effect is negligible at the reported accuracy level.","section":"§III.B and §II.D"}],"minor_comments":[{"comment":"The statement that the s-contributions 'compare well' with Koziol and Aucar is not supported by the reported deviations of up to 14.9% (Zn 1s), 22.9% (Cd 1s), and 35.9% (Hg 1s); please rephrase to describe the actual level of agreement and note explicitly that this comparison involves a different self-energy model (Welton picture).","section":"§IV.C, Table VI"},{"comment":"The spline interpolation accuracy is reported only as a mean relative error; please also report the maximum relative error in the grid range, since failures near the nucleus could disproportionately affect core-orbital QED matrix elements.","section":"§III.A"},{"comment":"Footnote c, which discusses a correction to a GRASP routine in Ref. [29], is lengthy and detailed; consider moving this explanation to the text or supporting information, as it interrupts the table of results.","section":"§IV.A, Table I"},{"comment":"The conclusion states that the impact of different self-energy schemes is 'much larger' than residual deviations between the present two-component framework and four-component frameworks; this is only shown for select examples and should be qualified as applying to the systems and potentials studied here.","section":"§V, Conclusion"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of the journal and the implementation appears useful, but the abstract overstates the agreement relative to the authors' own Table V. The requested decomposition of the 5p1/2 deviation and the point-nucleus test are feasible and should be obtained before acceptance. The paper would also benefit from the authors making the un-averaged Kramers-pair data and any supporting grid-convergence tests available in the supporting information."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe short version: this is a real implementation of QED effective potentials in a two-component ZORA framework, with new molecular data and a mostly honest benchmark set. The abstract oversells the agreement, but the method itself looks sound.\n\nWhat's new: the ZORA picture-change transformation applied to the Uehling and self-energy potentials, tested on atoms, ions, and molecules, including group 2 monofluorides and BaF/RaF transitions. The perturbative vs. self-consistent comparison for Au orbital energies is a useful addition.\n\nWhat it does well: benchmarks against four-component DHF and AOC-HF are extensive and mostly excellent—ionization energies within 0.1–2.6%, transition energies within 0.0–2.4%, valence orbital energies within 0.0–1.7%. The spline interpolation for the Flambaum–Ginges integrals is documented with error estimates. The paper is also candid about the one big outlier: the 5p1/2 contribution in Table V deviates by 63.8% from the four-component reference, with the origin stated as 'presently open.'\n\nThat outlier is the soft spot. The abstract claims 'excellent agreement' without qualification, which is too strong. The 5p1/2 deviation is on a small absolute contribution, and it's a single case, so I don't think it sinks the method. But it does mean the picture-change transformation is not validated for all orbitals, and the abstract should say so. The paper itself flags the deviation; the problem is only the unqualified claim in the abstract.\n\nThe comparison with Skripnikov et al. for BaF/RaF transitions shows 16–24% differences in QED contributions, attributed to different wave functions and SE formulations. That's plausible but not fully resolved.\n\nBottom line: the central methodology holds up, the new data are useful, and the paper deserves serious peer review. I'd ask for a revised abstract that acknowledges the Table V deviations and a brief discussion of whether the 5p1/2 case reflects a systematic limitation. Who's this for? Anyone working on heavy-element or radioactive-molecule precision spectroscopy, especially RaF.\n\nYes, send it to review.","headline":"Solid two-component ZORA implementation of QED potentials with strong benchmarks; the abstract overstates agreement given the unexplained Au 5p1/2 deviation.","tokens_in":22339,"tokens_out":2191,"would_cite":true,"duration_ms":21255,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["31.30.J-","31.15.-p"],"model":"deepseek-v4-flash","headline":"This paper demonstrates that vacuum polarisation and electron self-energy can be added as effective potentials in a two-component ZORA framework, matching four-component accuracy for heavy atoms and molecules.","keywords":["quantum electrodynamics","ZORA","two-component relativistic methods","vacuum polarisation","electron self-energy","Uehling potential","Flambaum-Ginges potential","heavy-element molecules"],"falsifier":"Recompute the Flambaum-Ginges self-energy contribution to the 5p1/2 Kohn-Sham orbital energy of gold in the same basis and with the same effective potentials, but using a different two-component transformation or the full four-component Hamiltonian; if the calculated value moves from the four-component reference by much less than 63.8%, the ZORA picture-change transformation is the cause.","tokens_in":21281,"feed_emoji":"⚛️","tokens_out":5040,"duration_ms":51263,"temperature":0.7,"pith_summary":"The paper aims to show that the two dominant one-electron quantum electrodynamic corrections, vacuum polarisation and electron self-energy, can be included as effective potentials inside a two-component zeroth-order regular approximation (ZORA) framework. If true, QED corrections for heavy atoms and molecules no longer require a full four-component Dirac machinery; they can be estimated at two-component cost with deviations of a few percent from four-component benchmarks for most valence, ionisation, and transition properties. The paper reports such corrections for ionisation energies of group 1, 2, 11, 12, 13, and 18 atoms, for Li-, Na-, and Cu-like ionic transitions up to Z = 90, and for BaF and RaF molecules, and compares perturbative and self-consistent treatments for gold.","feed_headline":"Two-component ZORA reproduces four-component QED corrections","feed_subtitle":"Effective potentials add vacuum polarisation and self-energy to heavy-atom molecules without four-component Dirac machinery.","key_machinery":"The carrying mechanism is the two-component ZORA Hamiltonian supplemented by one-electron QED effective potentials. The Uehling potential handles vacuum polarisation; the Flambaum–Ginges potential supplies magnetic, high-frequency, and low-frequency self-energy terms, with the high-frequency part regularised by a fitted cut-off factor; the Pyykkö–Zhao potential is a fitted Gaussian describing s-level self-energy shifts. The integrals in the magnetic and high-frequency parts are evaluated by cubic spline interpolation, and the whole set of operators is mapped into two-component form by the ZORA picture-change transformation of Ref. [60], which reconstructs the small component from the ZORA wave function.","core_discovery":"The central claim is that the Uehling potential for vacuum polarisation and effective one-electron potentials for the electron self-energy, namely the Flambaum–Ginges potential in its magnetic, high-frequency, and low-frequency parts and the Pyykkö–Zhao Gaussian potential, can be transformed into two-component ZORA form through the picture-change transformation that generates the small component from the ZORA wave function. With these potentials, QED corrections to orbital energies, ionisation energies, and electronic transition energies of heavy atoms and molecules agree with four-component Dirac-Hartree-Fock and average-of-configuration Hartree-Fock results to within a few percent for most cases. The largest deviations appear for deep core orbitals of gold and especially for the 5p1/2 orbital, where the deviation reaches 63.8% and the paper states that the origin is presently open.","pith_inferences":["If the 63.8% deviation for the gold 5p1/2 orbital reflects a picture-change failure of ZORA rather than a peculiarity of that orbital, the general agreement claim is not universal; repeating the comparison for 5p1/2 orbitals of other heavy atoms would settle this.","The same two-component machinery could be used to estimate QED corrections to properties beyond energies, such as hyperfine fields or parity-violation matrix elements, where short-range behaviour may amplify the deviations seen in core orbitals.","Because the paper reports that basis-set choice affects the QED corrections more than the level of theory, extending the benchmark to larger basis sets may tighten or shift the few-percent agreement for valence properties.","The method offers a practical route to include QED contributions in molecular dynamics or property calculations where four-component treatments remain prohibitively expensive."],"forward_implications":["QED corrections can be included in routine molecular electronic-structure calculations at two-component cost, making heavy-element spectroscopic predictions more complete.","The perturbative expectation-value treatment is largely sufficient: once linear response is accounted for, differences from self-consistent inclusion become negligible for most orbitals.","The results provide benchmark QED contributions to BaF and RaF transition energies, relevant for precision spectroscopy of radioactive molecules.","The Z-scaling fits for group trends allow one to estimate where QED corrections become significant, with fastest growth for group 11 and 12 elements.","For gold, QED corrections lift the HOMO orbital energy by about 0.25%, an effect relevant for meV-accuracy predictions such as ionisation potentials and electron affinities."],"supporting_citations":[{"why":"Provides the four-component Dirac-Hartree-Fock reference values for QED contributions to ionisation energies and ionic transition energies that the paper compares against.","marker":"[29]"},{"why":"Provides the four-component B3LYP/dyall.3zp Kohn-Sham orbital QED data for gold and the AOC-HF valence orbital data used as benchmarks.","marker":"[38]"},{"why":"Source of the Flambaum-Ginges effective self-energy potential with its magnetic, high-frequency, and low-frequency contributions.","marker":"[49]"},{"why":"Source of the Pyykkö-Zhao fitted Gaussian self-energy potential.","marker":"[50]"},{"why":"Supplies the numerical method used to evaluate the Uehling vacuum-polarisation potential.","marker":"[46]"},{"why":"Provides the ZORA picture-change transformation and matrix-element implementation that maps the QED potentials into two-component form.","marker":"[60]"},{"why":"Supplies the model-potential formulation of ZORA used to generate the small component from the ZORA wave function.","marker":"[61]"},{"why":"Supplies the response formalism used for perturbative QED corrections to orbital energies.","marker":"[59]"}],"fun_headline_variants":["Two-component ZORA captures QED effects for heavy atoms","Efficient QED corrections via ZORA match four-component results","Uehling and self-energy potentials add QED to ZORA framework","QED corrections in ZORA match four-component benchmarks","Efficient two-component QED for heavy elements and molecules"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the ZORA picture-change transformation turns the QED potentials into matrix elements that match four-component results for every orbital of interest; the 63.8% deviation for gold's 5p1/2 orbital, whose origin the paper says is open, shows this premise does not hold in at least one case.","fun_headline_variants_meta":{"raw":{"variants":["Two-component ZORA captures QED effects for heavy atoms","Efficient QED corrections via ZORA match four-component results","Uehling and self-energy potentials add QED to ZORA framework","QED corrections in ZORA match four-component benchmarks","Efficient two-component QED for heavy elements and molecules"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000508,"raw_usage":{"total_tokens":2499,"prompt_tokens":993,"completion_tokens":1506,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":609,"completion_tokens_details":{"reasoning_tokens":1422}},"tokens_in":609,"tokens_out":1506,"duration_ms":13303,"temperature":1.0,"reasoning_tokens":1422,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T21:50:53.626558+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the Flambaum-Ginges self-energy contribution to the 5p1/2 Kohn-Sham orbital energy of gold in the same basis and with the same effective potentials, but using a different two-component transformation or the full four-component Hamiltonian; if the calculated value moves from the four-component reference by much less than 63.8%, the ZORA picture-change transformation is the cause.","supporting_citations":[{"cited_title":"Kozio l, C","cited_arxiv_id":null,"evidence_quote":"Provides the four-component B3LYP/dyall.3zp Kohn-Sham orbital QED data for gold and the AOC-HF valence orbital data used as benchmarks."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the numerical method used to evaluate the Uehling vacuum-polarisation potential."},{"cited_title":"van W¨ ullen, J","cited_arxiv_id":null,"evidence_quote":"Supplies the model-potential formulation of ZORA used to generate the small component from the ZORA wave function."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the response formalism used for perturbative QED corrections to orbital energies."}],"review_version":1}