{"id":"56e54983-d0aa-46da-8549-25c991a7e329","arxiv_id":"2411.12427","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Ultra-precise finite-element solutions of the two-center Dirac equation are reported for H2+ and Th2^179+, matching recent independent values but with disputed uncertainty estimates.","lead":"A researcher solved the Dirac equation for an electron around two nuclei with extreme numerical precision, using a minmax finite-element method on the hydrogen molecular ion and a highly charged thorium quasi-molecular ion. The results match a recent independent calculation and could improve future QED and g-factor predictions for molecular ions.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Claimed 1e-21 fractional uncertainty for Th2^179+ is unsupported by the paper's own convergence data and the 5e-18 discrepancy with [19].","rationale":"The most load-bearing claim is the benchmark precision for both systems. H2+ is independently corroborated, so the entire weight falls on Th2^179+. The paper's own text admits that the singularity error is not fully controlled for Z=90 and that error cancellation is less efficient, which undermines the extrapolation's assumed single-power-law behavior. The Table IV Dmax scatter and the 5e-18 difference from Nogueira and Karr provide concrete evidence that the true error is larger than the claimed 1e-21 fractional uncertainty. This concern directly affects the central claim: the quoted Th uncertainty is not defensible. The reader's weakest assumption identified the same point, so I agree. The appropriate verdict remains conditional: the H2+ result and the Th value in the digits that agree with [19] are useful benchmarks, but the uncertainty estimates for Th need to be revised to reflect the actual scatter or supported with a rigorous error analysis. This is consistent with the reader's conditional verdict, so no change is recommended.","tokens_in":13426,"tokens_out":5077,"duration_ms":46832,"concrete_test":"Refit the Th2^179+ energy sequence in Table III with a two-term error expansion, e.g., E(N) = E_inf + a N^{-q1} + b N^{-q2} (q1 approximately 9.3, q2 > q1), and also refit using only the first 10 and only the last 10 grid rows. If the extrapolated E_inf shifts by more than about 1e-18 atomic units, or the residuals reject the single-power-law model, then the claimed fractional uncertainty of 1e-21 must be revised upward to at least the observed scatter (~5e-18), matching the discrepancy with [19].","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim for H2+ is secure: it matches the independent Nogueira-Karr result to ~7e-30. The load-bearing weakness is the Th2^179+ uncertainty estimate. Section III and Fig. 2 state that for high Z 'the singularity error outweighs the convergence order' and that error cancellation is less efficient, yet the paper quotes a fractional uncertainty of ~1e-21 in the shift, i.e., ~5.7e-19 absolute. Table IV's Dmax scan scatters over ~3.6e-18 in Erel, and the independent result [19] differs by 5e-18, an order of magnitude larger than the implied absolute uncertainty. The extrapolation assumes a single power law q=9.3 fitted to the last few points; no rigorous bound or second-order term is given. The claimed precision therefore rests on the unsupported assumption that the singularity error and fit error are below the quoted level, contrary to the paper's own statements. The Th benchmark value itself may be correct in its agreed digits, but the uncertainty statement is overstated.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper computes the ground-state (1σg) solutions of the two-center Dirac equation for H2+ (R = 2) and Th2^179+ (R = 2/90) using the minmax functional and a p = 10 finite-element discretization in prolate spheroidal coordinates, with a singular coordinate transformation of order ν to regularize the Coulomb singularities. Systematic grid sequences (Tables I and III) yield Erel, Enrel, and the relativistic shift; power-law extrapolation with fitted orders q ≈ 9.7 and q ≈ 9.3 produces the final values in Table V, with claimed fractional uncertainties of ~10^-23 and ~10^-21 in the two shifts. The results agree with the independent calculation of Nogueira and Karr [19] to ~7×10^-30 (H2+) and ~5×10^-18 (Th2^179+). Additional cross-checks vary the domain size Dmax (Tables II, IV), the number of integration points per element (Fig. 2b), and the transformation order ν.","tokens_in":13585,"tokens_out":33694,"duration_ms":270807,"significance":"If the results hold at the stated precision, these are benchmark-quality reference values for two-center relativistic one-electron systems, of use for QED corrections, g-factor studies, and method benchmarking in molecular relativistic quantum chemistry. The manuscript's strengths are concrete: systematic grid-convergence tables (Tables I and III), cross-checks against domain size (Tables II and IV), integration-point count (Fig. 2b), and transformation order; verification against an independent calculation [19]; and a minmax formulation whose raw FEM values converge from above without spurious states. The H2+ claim is fully supported by the ~7×10^-30 agreement with [19]. The Th2^179+ uncertainty estimate is not supported by the manuscript's own data (Major Comment 1), and the α-uncertainty propagation is misstated (Major Comment 2); both issues are localized and fixable within the scope of the paper.","major_comments":[{"comment":"The claimed fractional uncertainty of ~10^-21 in the Th2^179+ relativistic shift (absolute ~5.7×10^-19 au in a shift of -573.42 au) is not supported by the manuscript's own data. In Table IV, the extrapolated Erel values scatter over Dmax ∈ [0.30, 0.40] by ~3.6×10^-18 au, and the independent calculation [19] differs from the Table V value by 5×10^-18 au in Erel, an order of magnitude above the claimed shift uncertainty. The two internally consistent routes to the final energy also disagree: the direct ν = 10 extrapolation (Table III, -9504.756648434009500748) and the ν = 2 nonrelativistic value plus the extrapolated shift (Table V, -9504.756648434009500737) differ by ~1.1×10^-17 au, which exceeds the stated ~10^-18 energy uncertainty by an order of magnitude. The extrapolation assumes a single power law with q ≈ 9.3 fitted to the last grid points (Fig. 2), and the text itself states that for high Z 'the singularity error outweighs the convergence order' and that error cancellation in the shift is less efficient (Sec. III). The sentence that the discrepancy scatters around ≲10^-18 also conflicts with the reported 5×10^-18 difference from [19]. The authors should provide an error budget that accounts for the Table IV scatter, the route dependence, and the comparison with [19], or reduce the claimed precision to a level consistent with these data.","section":"Sec. III; Tables III–V; Fig. 2"},{"comment":"The propagation of the CODATA 2018 α uncertainty is misstated. With α^-1 = 137.035999084(21), the relative uncertainty is δα/α ≈ 1.5×10^-10. Since the relativistic shift scales as α^2 at fixed Z at leading order, δ(ΔE)/ΔE ≈ 2δα/α ≈ 3×10^-10, giving δ(ΔE) ≈ 2×10^-15 au for H2+ and ≈ 1.7×10^-7 au for Th2^179+. This contradicts the text's claim of 'uncertainty of order 10^-20 in the obtained energies' and the statement that the α-uncertainty is smaller for Th2^179+ than for H2+. In fact the α-induced fractional uncertainty (~3×10^-10) exceeds the claimed computational fractional uncertainties (10^-23 and 10^-21) by many orders of magnitude for both systems. The authors should correct this analysis and state explicitly that the quoted fractional uncertainties are computational, valid at the fixed CODATA 2018 α value, rather than uncertainties of the physical values.","section":"Sec. III (α-uncertainty paragraph)"}],"minor_comments":[{"comment":"In Table III, the element counts '200/1020' and '1152/5808' appear to be missing a digit; they should read 200/10201 and 1152/58081 to match the grid sequence used for H2+ in Table I.","section":"Table III"},{"comment":"The provenance of each Table V entry should be specified: the tabulated Erel (-9504.756648434009500737) is not exactly the sum of the tabulated Enrel and the shift (-8931.3371374090663382226 - 573.4195110249431625138 = -9504.7566484340095007364), so the reader cannot reproduce the final values without knowing the parameter sets (ν, Dmax) used for each column.","section":"Table V"},{"comment":"There are several typos and grammatical slips: 'minamx' (Sec. II), 'the the subspace of positronic states' (Sec. II), an unbalanced parenthesis in 'Z < 1/α < 137.036..)' (Sec. II), 'whrer' (Fig. 2 caption), 'blow the precision' for 'below the precision' (Sec. III), and missing spaces in the Conclusion paragraph ('Weshowedsystematicallyaccuratevaluesbyinvestigationoftheconvergence').","section":"Secs. II–III; Fig. 2; Conclusion"},{"comment":"The paper reports the slopes of the linear fits in Figs. 1 and 2 but does not state how many grid points enter each fit or how the extrapolated value changes when the fit window is varied; specifying the fit range and the sensitivity to it would support reproducibility of the extrapolation step.","section":"Figs. 1–2"},{"comment":"The label 'Rel. eff.' in Table V is not defined; the text uses 'relativistic shift', 'relativistic effect', and 'Rel. eff.' interchangeably, and the notation should be unified.","section":"Table V"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a methods/benchmark contribution whose core novelty over the author's prior work [8] is the higher precision and the heavy-Z uncertainty claim; the method itself rests on refs [6,7,8], all authored or co-authored by the author, which makes the Th uncertainty estimate difficult to verify independently of the author's own error-analysis conventions. If the two major comments are addressed, the paper is a solid contribution for a molecular-physics or PRA-type journal. I would also encourage the author to add a data-availability statement and to resolve the inconsistency between the claimed Th precision and the 5×10^-18 agreement level with [19] before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nQuick take: this is a careful numerical paper that deserves a referee, but the lead claim for Th2^179+—a fractional uncertainty of 1e-21 in the relativistic shift—doesn't hold up against the paper's own data. The H2+ result, by contrast, is in excellent shape.\n\nWhat's new: the author applies his established minmax FEM machinery to push the two-center Dirac ground state to more digits, and to a heavy quasi-molecular ion (Th2^179+). The convergence study is serious: systematic grid sequences, domain-size scans (Tables II and IV), integration-point tests, and a comparison to the independent Nogueira-Karr calculation. The H2+ energy matches [19] to 7e-30, which is a strong independent confirmation. The non-relativistic limit and counterpoise-style error cancellation are handled sensibly. Quadruple precision is used. All of this is real, reproducible evidence, even though no code or data files are shipped.\n\nThe soft spot is the thorium uncertainty. The paper itself says (Section III, Fig. 2) that for high Z the singularity error outweighs the convergence order and that error cancellation is less efficient. Yet Table V quotes a fractional uncertainty of ~1e-21 in the shift—about 5.7e-19 absolute. The independent value [19] differs by 5e-18, an order of magnitude larger, and the Dmax scan in Table IV scatters over ~3.6e-18. The extrapolation uses a single power law q=9.3 fitted to a few points, with no estimate of higher-order terms or a rigorous bound. So the claimed uncertainty for Th2^179+ is not supported. The value in its middle digits is likely correct, but the error bar is too small.\n\nThe citation pattern is heavily self-referential (method papers [6,7,8]), which is acceptable since the method is the author's; it's not a hidden circularity. The H2+ result itself is already in [19] to more digits, so the novelty rests on the precision benchmark, which is only partly new.\n\nBottom line: this deserves peer review. A competent referee should ask the author to either lower the thorium uncertainty to match the actual spread and agreement with [19], or provide a real error analysis for the extrapolation and singularity treatment. The H2+ part, and the method demonstration for heavy ions, are solid enough for publication after revision.","headline":"H2+ benchmark is solid and agrees with independent work to 7e-30, but the Th2^179+ uncertainty claim is overstated by roughly an order of magnitude and needs revision.","tokens_in":14196,"tokens_out":2856,"would_cite":false,"duration_ms":27090,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65N30","81Q05","81V45"],"pacs":["31.15.-p","31.30.J-","02.70.Dh"],"model":"deepseek-v4-flash","headline":"The two-center Dirac equation is solved by the minmax finite-element method to benchmark accuracy, with estimated fractional uncertainties of about 1e-23 in the relativistic shift for H2+ and 1e-21 for Th2^179+.","keywords":["two-center Dirac equation","finite element method","minmax principle","relativistic shift","H2+ molecular ion","heavy quasi-molecular ion Th2^179+","Coulomb singularity","high-precision benchmark"],"falsifier":"Refine both systems beyond the densest grids used here (1800 elements, 90601 points), for ${\\rm Th}_2^{179+}$ also with a higher $\\nu$ value; if the extrapolated energies move by more than about $10^{-28}$ au for ${\\rm H}_2^+$ or $10^{-18}$ au for ${\\rm Th}_2^{179+}$, the claimed benchmark uncertainties are falsified.","tokens_in":13106,"feed_emoji":"⚛️","tokens_out":9741,"duration_ms":85864,"temperature":0.7,"pith_summary":"The paper reports a benchmark-quality solution of the two-center Dirac equation, the relativistic one-electron problem of two Coulomb centers that models molecular hydrogen ions and heavy quasi-molecular ions. The method combines the minmax variational principle, which projects out the negative-energy (positronic) continuum, with a high-order finite-element discretization and a singular coordinate transformation that regularizes the point-nucleus Coulomb singularity. For ${\\rm H}_2^+$ at $R=2$, the final relativistic ground-state energy is $-1.10264158103257716411812499995$ atomic units, and for ${\\rm Th}_2^{179+}$ at $R=2/90$ it is $-9504.756648434009500737$ atomic units; the corresponding relativistic shifts carry estimated fractional uncertainties of about $10^{-23}$ and $10^{-21}$. If these values are correct, they provide reference energies for QED corrections and bound-electron g-factor calculations in molecular ions.","feed_headline":"Minmax-FEM solves two-center Dirac equation to 1e-23","feed_subtitle":"Finite-element minmax method pins H2+ and heavy Th2^179+ energies with fractional uncertainties down to 1e-23.","key_machinery":"The load-bearing machinery is the minmax energy functional, where the small spinor component is eliminated and the large-component equation becomes a nonlinear eigenvalue problem solved iteratively. The two-dimensional finite-element calculation uses complete polynomial shape functions of order $p=10$ in prolate spheroidal coordinates $(\\xi,\\eta)$, together with the singular coordinate transformation of Eq. (A1), whose order $\\nu$ controls the point density near the nuclei and regularizes the Coulomb singularity of the form $r^{-1+\\gamma_{l,\\kappa}}$. The relativistic shift is computed as $E_{\\rm rel}-E_{\\rm nrel}$ on identical grids, which cancels most of the smooth discretization error, and a power-law extrapolation over the grid sequence with the fitted convergence order $q$ as the leading error exponent produces the final energies. This combination of projection, singularity regularization, and error cancellation is what carries the claim.","core_discovery":"The central claim is that the minmax finite-element method solves the two-center Dirac equation with far higher precision than earlier relativistic molecular calculations: absolute energy accuracy near $10^{-28}$ au for ${\\rm H}_2^+$ and $10^{-18}$ au for ${\\rm Th}_2^{179+}$, with the relativistic shift (the difference between Dirac and Schrödinger energies) even more accurate because the non-relativistic and relativistic computations on the same grids share the same systematic errors. The energies converge from above, no spurious states appear, and extrapolation over grid refinements with fitted convergence orders $q \\approx 9.7$ and $q \\approx 9.3$ yields the final values. The results agree with a recent independent high-precision calculation to about $7\\times10^{-30}$ au for ${\\rm H}_2^+$ and $5\\times10^{-18}$ au for ${\\rm Th}_2^{179+}$, and the paper argues that the minmax solution is cleaner in the strongly relativistic regime because it is free from positronic contamination.","pith_inferences":["One testable extension is to go beyond the densest grids used here (1800 elements, 90601 points) or to raise $\\nu$ for ${\\rm Th}_2^{179+}$; if the extrapolated energy moves by more than the claimed $10^{-18}$ au, the uncertainty estimate would need revision.","The paper's uncertainty model assumes a single leading power-law error; if higher-order or oscillatory terms contribute at the last digits, the extrapolation could be biased, which a deliberately varied grid-sequence study could reveal.","Applying the same method at several internuclear distances $R$ would turn the single-point benchmark into an interpolation table for molecular spectroscopy and QED calculations.","Because the point-nucleus model is used, the heavy-ion value is a pure mathematical benchmark; adding finite-nuclear-size corrections would shift it at levels likely larger than the claimed numerical uncertainty, so direct experimental comparison needs that separate step."],"forward_implications":["The ${\\rm H}_2^+$ ground-state energy at $R=2$ becomes a fixed reference with about 28 decimal digits, useful for benchmarking other relativistic methods and for QED one-loop self-energy calculations.","The ${\\rm Th}_2^{179+}$ energy at $R=2/90$ improves on earlier heavy quasi-molecular results by many orders and provides a point-nucleus benchmark for strongly relativistic two-center systems.","The combination of minmax projection and singular-coordinate FEM should transfer to other internuclear distances and angular momentum channels, giving controlled convergence from above.","The error-cancelled relativistic shift is the quantity entering g-factor and radiative-correction calculations in molecular hydrogen ions, so the improved precision directly tightens theory-experiment comparisons.","The computational cost (quadruple precision, grids up to 90601 points) is moderate enough that the paper expects extension to multi-electron systems in Hartree-Fock or DFT frameworks."],"supporting_citations":[{"why":"Establishes the minmax variational characterization of Dirac eigenvalues that excludes the negative continuum.","marker":"[1]"},{"why":"Presents the minmax FEM formulation, the nonlinear iteration expansion, and the extrapolation procedure used here.","marker":"[6]"},{"why":"Earlier minmax FEM two-center results and the error-cancellation idea for the relativistic shift.","marker":"[7]"},{"why":"Previous high-precision result for H2+ that the present work improves by many orders.","marker":"[8]"},{"why":"Introduces the singular coordinate transformation that regularizes the Coulomb singularity in the two-center FEM.","marker":"[11]"},{"why":"Supplies the CODATA 2018 value of alpha used in the calculations.","marker":"[15]"},{"why":"Details the singular transformation and grid-distribution strategy that controls the singularity error.","marker":"[16]"},{"why":"Independent recent high-precision result against which the final values are compared in Table V.","marker":"[19]"}],"fun_headline_variants":["Minmax finite-element solves two-center Dirac to 1e-23","Two-center Dirac solved to 1e-23 by minmax FEM","Minmax FEM reaches 1e-23 in two-center Dirac","High-precision two-center Dirac solved to 1e-23","Minmax FEM solves two-center Dirac equation at 1e-23"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The uncertainty estimates assume that the extrapolated finite-element energies converge to the exact eigenvalue as a single power law with the fitted order $q$, and that the singular coordinate transformation with $\\nu=10$ fully controls the Coulomb singularity for ${\\rm Th}_2^{179+}$; the paper itself notes that for high $Z$ the singularity error outweighs the convergence order.","fun_headline_variants_meta":{"raw":{"variants":["Minmax finite-element solves two-center Dirac to 1e-23","Two-center Dirac solved to 1e-23 by minmax FEM","Minmax FEM reaches 1e-23 in two-center Dirac","High-precision two-center Dirac solved to 1e-23","Minmax FEM solves two-center Dirac equation at 1e-23"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000814,"raw_usage":{"total_tokens":3541,"prompt_tokens":892,"completion_tokens":2649,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":508,"completion_tokens_details":{"reasoning_tokens":2555}},"tokens_in":508,"tokens_out":2649,"duration_ms":18923,"temperature":1.0,"reasoning_tokens":2555,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T17:32:09.738395+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Refine both systems beyond the densest grids used here (1800 elements, 90601 points), for ${\\rm Th}_2^{179+}$ also with a higher $\\nu$ value; if the extrapolated energies move by more than about $10^{-28}$ au for ${\\rm H}_2^+$ or $10^{-18}$ au for ${\\rm Th}_2^{179+}$, the claimed benchmark uncertainties are falsified.","supporting_citations":[{"cited_title":"Dolbeault, M","cited_arxiv_id":null,"evidence_quote":"Establishes the minmax variational characterization of Dirac eigenvalues that excludes the negative continuum."},{"cited_title":"Kullie, inDissertation (Thesis), Universität Kassel (http://nbn-resolving.de/urn:nbn:de:hebis:34-1835, 2004)","cited_arxiv_id":null,"evidence_quote":"Presents the minmax FEM formulation, the nonlinear iteration expansion, and the extrapolation procedure used here."},{"cited_title":"Kullie, D","cited_arxiv_id":null,"evidence_quote":"Earlier minmax FEM two-center results and the error-cancellation idea for the relativistic shift."},{"cited_title":"Kullie and S","cited_arxiv_id":null,"evidence_quote":"Previous high-precision result for H2+ that the present work improves by many orders."},{"cited_title":"Kullie and D","cited_arxiv_id":null,"evidence_quote":"Introduces the singular coordinate transformation that regularizes the Coulomb singularity in the two-center FEM."},{"cited_title":"Kullie and D","cited_arxiv_id":null,"evidence_quote":"Details the singular transformation and grid-distribution strategy that controls the singularity error."},{"cited_title":"Looking at table V on finds an excellent agreement with the result of Nogueira et","cited_arxiv_id":null,"evidence_quote":"Independent recent high-precision result against which the final values are compared in Table V."}],"review_version":1}