{"id":"2631daae-fb44-45a9-a7d3-c9ef512658c5","arxiv_id":"1910.02772","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":3.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":16,"one_line_summary":"Ground-state electronic energy is claimed to be nearly linear in the electron-electron coupling strength, so the real energy can be approximated from the a=0 non-interacting solution, but with accuracy below chemical standards.","lead":"This paper scales the electron-electron repulsion in the Schrödinger equation by a parameter a, from zero (no repulsion) to one (real world), and argues that electronic energy rises almost linearly across that range. The possible payoff is a cheap way to estimate real molecular energies from a simple non-interacting reference, although the paper's own numbers show the estimate is not yet accurate enough for chemistry.","discovery_kind":"incremental","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central TNRS transfer (Eq.20) is an uncontrolled first-order approximation; the paper's own observation that larger basis sets increase curvature means the quasi-linearity it depends on is not established in the physically relevant limit.","rationale":"The exact identities (Eqs.4, 13, 18) are mathematically sound; the Hellmann-Feynman decomposition is correct. The concern is not internal inconsistency but the unsupported step from the exact identity Eq.13 to the approximate Eq.20. The paper's own Fig. 1 and text admit that curvature increases with basis size, so the empirical basis for quasi-linearity is thinner than the central claim requires. Since the reader already identified this as the weakest assumption and assigned CONDITIONAL, my stress-test does not shift the verdict; it sharpens the test needed. I also note that the wavefunction-level 'HK extension' is not needed for Eq.20 and is complicated by a=0 degeneracies, but that is secondary. The Eq.33 coefficients are fitted in-sample, so they do not independently validate the transfer. If the proposed FCI/a scan shows small, basis-stable error, Eq.20 would be a useful parameter-free estimate; if not, the algebraic-transfer claim is not supported beyond minimal-basis examples.","tokens_in":18192,"tokens_out":4323,"duration_ms":43078,"concrete_test":"For H2 and H2O, compute FCI-quality energies with Hee scaled by a (0 ≤ a ≤ 1) in STO-3G, cc-pVDZ, and cc-pVQZ. From the a=0 HF wavefunction get Y0, eelectr,0, and <Y0|Hee|Y0>; compare E_exact(1) - eelectr,0 against Eq.20 and against the exact integral (N(N-1)/2)∫_0^1 <y0(a)|r12^-1|y0(a)> da, where y0(a) is the FCI ground state at coupling a. If the Eq.20 error exceeds ~1 kcal/mol or grows monotonically with basis size, the quasi-linear approximation fails in the limit that matters; if it stays below chemical accuracy across all three bases, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Eq.13 is exact: Eelectr,0 - eelectr,0 = (N(N-1)/2)∫_0^1 <y0(a)|r12^-1|y0(a)> da. Eq.20 replaces the integrand by its a=0 value, <Y0|Hee|Y0>, i.e. it assumes the Hellmann-Feynman derivative is constant over [0,1]. This is not proven; the paper's justification is only the empirical smallness of <y0|r12^-1|∂y0/∂a> (reasoning after Eq.17). The exact relation Eq.19 contains the overlap denominator <y0(a)|Y0> and an off-diagonal matrix element; Eq.20 silently sets both to their a=0 diagonal limits, and no error bound is given for either replacement. The paper itself states that the curve is 'less linear with increasing basis set' (after Eq.17), and Fig. 1 is STO-3G only. Since chemical predictions require basis-set limits, the load-bearing assumption is exactly where the evidence is weakest. The later Eq.33 fit is made to the same 149-molecule G3 set that is used for evaluation, so it cannot serve as out-of-sample validation of the transfer.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the family of non-relativistic electronic Hamiltonians H(a) = Hkin + Hne + aHee. It derives exact energy relations connecting the physical system (a=1) with the totally non-interacting reference system (TNRS, a=0), most notably Eqs. (4), (13), and (18). The central claim is that enrgelectr,0(a) is quasi-linear in a, so that the physical ground-state energy can be estimated from the a=0 wavefunction via Eq. (20), Eelectr,0 ≈ eelectr,0 + (N(N-1)/2)<Y0|r12^-1|Y0>. The paper presents STO-3G numerical evidence for quasi-linearity (Fig. 1), an empirical correction formula (Eq. 33) fitted to G3 energies, and discusses extensions of the first Hohenberg-Kohn theorem, Koopmans' theorem, Hund's rule, and the virial theorem to general coupling strength a.","tokens_in":18706,"tokens_out":4219,"duration_ms":41902,"significance":"The exact identities in Eqs. (4), (13), and (18) are correct and potentially useful; they quantify the interaction-energy contribution in an internally consistent way and connect the non-interacting and physical ground states through an exact adiabatic integral. If the quasi-linearity conjecture were established with controlled accuracy, the TNRS transfer in Eq. (20) would be a parameter-free first-order estimate of correlation-free but interaction-inclusive energies, which is an attractive idea. The manuscript also provides a concrete computational protocol (HF-SCF/basis/a=0) that is faster and more stable than standard a=1 SCF. However, the paper's central quantitative claims currently rest on a minimal-basis empirical observation and an in-sample fit; the main approximation is not proven and is admitted to worsen with larger basis sets. The correct exact identities are the strongest part of the contribution.","major_comments":[{"comment":"The central approximation, Eq. (20), replaces the exact integrand in Eq. (13) by its a=0 value, which is equivalent to assuming that <y0(a)|r12^-1|y0(a)> is constant on [0,1]. The manuscript's justification is the smallness of <y0|r12^-1|∂y0/∂a> in Eq. (17), but no error bound or quantitative measure of this smallness is provided. The paper itself states that the curve is less linear with increasing basis set, and Fig. 1 shows only STO-3G results for seven molecules. Since the basis-set limit is the relevant regime for chemical predictions, the load-bearing assumption is unsupported exactly where it matters most.","section":"§Calculating ground state with a=0 (Eq. 20 and the discussion following Eq. 17)"},{"comment":"The coefficients in Eq. (33) are obtained by least-squares fitting to the same 149 G3 ground-state molecular energies that are then used to report the average and maximum absolute deviations (1.615905 h and 7.015398 h for L=2). This is an in-sample fit, not a prediction; the reported 'improvement' is therefore not evidence for the validity of the TNRS transfer or for the correction formula. An out-of-sample test (e.g., cross-validation or a separate test set) is needed before Eq. (33) can be presented as a meaningful calibration.","section":"§Computation properties of TNRS (Eq. 33)"},{"comment":"The second empirical assumption, that 'the LCAO coefficients are close to each other between Y0 and Ψ0' (i.e., quasi-independence of the LCAO coefficients on a), is used to justify both Eq. (20) and the w-based representation in Eqs. (26)-(29), but no systematic measurement, error estimate, or basis-set study is given. This assumption is not derived and is at least as consequential as the quasi-linearity of the energy; the paper should either quantify its accuracy or explicitly mark it as a conjecture.","section":"§Calculating ground state with a=0 (paragraph beginning 'Importantly, ∂enrgelectr,0(a)/∂a ≈ const.')"},{"comment":"The proposed extension of the first Hohenberg-Kohn theorem, Y0(a=0) ⇔ Hne ⇔ Ψ0(a=1), relies on the statement that the a=0 ground-state density determines the nuclear framework. This requires conditions of non-interacting v-representability and non-degeneracy that are not discussed. If the theorem is intended as a rigorous extension, those conditions must be stated; if it is intended as a practical mapping, the limitations should be acknowledged.","section":"§Calculating ground state with a=0 (Eqs. 24-25)"}],"minor_comments":[{"comment":"The coupling strength parameter is denoted 'a' in the Hamiltonian, but Eq. (33) uses 'a' for expansion coefficients; this notational collision makes the text confusing and should be resolved.","section":"Introduction, Eq. (2)"},{"comment":"The caption admits that 'larger basis set yields slightly larger curvature (not shown)'; showing at least one larger-basis-set curve would allow the reader to assess the extent of curvature growth and is important for the quasi-linearity claim.","section":"Figure 1 caption"},{"comment":"The technical modification to the Gaussian SCF algorithm (changing 1/rij to a/rij with a single line) is described only in words; for reproducibility, the authors should provide the precise patch, the exact version of Gaussian, and the SCF convergence settings used.","section":"§Computational protocol"},{"comment":"Reference [1] is a bare link to an author's arXiv/chemrxiv page rather than a citable work; a specific published paper should be cited. Several other references are self-citations; the authors should ensure that the most relevant independent literature on adiabatic connection and coupling-strength integration is cited.","section":"References"},{"comment":"The table lists CI energies for C, N, and O in the third column, but the basis set and method for those CI values are not specified in the table or the text; this makes the comparison difficult to evaluate.","section":"Table 1"}],"recommendation":"major_revision","confidential_remarks":"The exact identities in the first half of the paper are sound and could be a useful reference for future work on coupling-strength integration. The main risk is that the paper's headline claims (quasi-linearity, algebraic transfer, and the correction formula) are supported only by STO-3G data and an in-sample fit. I would suggest the editor ask for (i) out-of-sample validation of Eq. (33), (ii) a quantitative study of curvature vs. basis set, and (iii) a clear statement of the status of the quasi-linearity assumption, before further consideration. The high density of self-citations and the non-standard formatting may also require editorial attention."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, what you should know: this paper is mostly a repackaging of first-order perturbation theory in the electron-electron interaction, dressed up as a coupling-strength transfer. The central equation, Eq.20, is simply the linear extrapolation from the non-interacting reference using the Hellmann–Feynman derivative at a=0. The paper's own data show it is not chemically accurate, and the later \"improvement\" (Eq.33) is a least-squares fit to the same 149 G3 molecular energies used to report average and maximum deviations. That makes the claimed improvement a fitting error, not a predictive result.\n\nThe paper does some things correctly. The exact identities—Eq.4, Eq.13, and Eq.18—are right and it is useful to have them collected. The computational trick of scaling 1/rij in the SCF is simple and could be handy. The empirical observation that vee(a) is roughly linear in a for a handful of small molecules in STO-3G is mildly interesting, though it is a thin basis for a general claim. The author honestly notes that larger basis sets increase curvature and that Eq.20 is not accurate enough for chemical purposes, and he calls the Ne back-transfer \"accidental.\" That honesty counts for something.\n\nThe soft spots are substantial. Eq.20 replaces the integrand in the exact Eq.13 by its a=0 value, an uncontrolled first-order approximation with no error bound. The quasi-linearity it relies on is only demonstrated for minimal-basis STO-3G, and the author's own admission undercuts it in the basis-set limit. The \"extension\" of the first Hohenberg-Kohn theorem (Eqs.24-25) is just a composition of two standard HK bijections; it is not new. The generalizations of virial, Koopmans, and Hund's rule are straightforward and do not add much. The citation pattern is heavily self-referential, with a half-dozen of the author's own papers cited; that is not disqualifying, but it does not help the novelty argument.\n\nWho is this for? Possibly someone teaching the adiabatic connection or testing first-order perturbation theory in a minimal basis. As a research contribution, the practical claim is unsupported and the evaluation circular. I would not send this to a serious referee; it deserves a bench rejection. If the author wants to make a case, the next version would need out-of-sample validation, larger-basis evidence, and a clear separation between fitted parameters and predicted energies.","headline":"Essentially first-order perturbation theory in disguise; the only empirical observation is limited to STO-3G, and the 'improved' formula is an in-sample fit.","tokens_in":19127,"tokens_out":6725,"would_cite":false,"duration_ms":62065,"reading_group":"no","serious_thinker":"yes","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper claims that switching on electron–electron repulsion changes electronic energy almost linearly, so a zero-repulsion starting state can estimate the interacting energy.","keywords":["coupling strength parameter","quasi-linear energy curve","totally non-interacting reference system","Hohenberg-Kohn theorem","electron-electron repulsion","Hellmann-Feynman theorem","Hund's rule","Koopmans theorem"],"falsifier":"Compute $E(a)$ for a small atom or molecule at $a=0,0.25,0.5,0.75,1$ with a large-basis correlated method, so that basis error no longer hides curvature; if the intermediate points deviate from the straight chord between $a=0$ and $a=1$ by more than chemical tolerance, the quasi-linearity behind Eq. (20) fails. The paper's own two-electron model $z=\\exp(r_{12}/2)$ offers a place where this comparison can be made against an exact analytic solution.","tokens_in":17924,"feed_emoji":"⚛️","tokens_out":12365,"duration_ms":110404,"temperature":0.7,"pith_summary":"This paper studies the family of electronic Hamiltonians $H(a)=H_{\\rm kin}+H_{\\rm ne}+aH_{\\rm ee}$ obtained by scaling the electron–electron repulsion by a coupling parameter $a$, with $a=1$ the physical case and $a=0$ a fully non-interacting system. It claims that the ground-state electronic energy $E(a)$ is quasi-linear in $a$: the derivative $dE/da$ is nearly constant on the interval $[0,1]$. If that is right, the exact identity $E_{\\rm electr,0}=e_{\\rm electr,0}+(N(N-1)/2)\\langle\\Psi_0|r_{12}^{-1}|Y_0\\rangle/\\langle\\Psi_0|Y_0\\rangle$ collapses to a parameter-free first-order estimate, $E_{\\rm electr,0}\\approx e_{\\rm electr,0}+(N(N-1)/2)\\langle Y_0|r_{12}^{-1}|Y_0\\rangle$, using only the non-interacting Slater determinant $Y_0$. This matters because the $a=0$ calculation is a single cheap SCF step, yet the estimate would give a useful starting point for the full correlated energy and a concrete bridge between the non-interacting and physical ground states.","feed_headline":"Interacting energy from a non-interacting wavefunction in one step","feed_subtitle":"Scaling the electron repulsion strength gives a nearly straight energy curve, so a zero-repulsion state can estimate the real one.","key_machinery":"The load-bearing object is the scaled Hamiltonian $H(a)=H_{\\rm kin}+H_{\\rm ne}+aH_{\\rm ee}$ and its ground-state energy as a function of the coupling strength $a$. The Hellmann–Feynman theorem fixes the slope of that curve as $dE/da=\\langle y_0(a)|H_{\\rm ee}|y_0(a)\\rangle$, so the quasi-linearity claim is precisely the statement that this slope barely changes over $[0,1]$. The second piece of machinery is the totally non-interacting reference system (TNRS), the $a=0$ endpoint where the wavefunction is an exact single Slater determinant $Y_0$ and the one-electron equations decouple; combined with the exact identity $E_{\\rm electr,0}=e_{\\rm electr,0}+(N(N-1)/2)\\langle\\Psi_0|r_{12}^{-1}|Y_0\\rangle/\\langle\\Psi_0|Y_0\\rangle$, it turns into the linear transfer Eq. (20). The paper also introduces a real symmetric correlation factor $w$ through $y_0(a)=wY_0$ to describe evolution between the endpoints, and an $L$th-order algebraic expansion in Eq. (33) for more accurate corrections on top of the same transfer.","core_discovery":"The central claim has two parts. First, the electron–electron repulsion contribution $v_{\\rm ee}(a)=a(N(N-1)/2)\\langle y_0(a)|r_{12}^{-1}|y_0(a)\\rangle$ is a quasi-linear function of $a$, so the ground-state energy curve is almost straight; the curvature is controlled by the small quantity $\\langle y_0(a)|r_{12}^{-1}|\\partial y_0(a)/\\partial a\\rangle$. Second, the first Hohenberg–Kohn theorem extends across the coupling strength: the non-interacting ground state $Y_0(a=0)$ and the physical ground state $\\Psi_0(a=1)$ carry equivalent information about the nuclear framework, summarized as $\\Psi_0(a=1)\\leftrightarrow H_{\\rm ne}\\leftrightarrow Y_0(a=0)$. The transfer formula, Eq. (20), writes the physical electronic energy as the non-interacting energy plus $(N(N-1)/2)\\langle Y_0|r_{12}^{-1}|Y_0\\rangle$, a Coulomb/exchange term evaluated entirely with the $a=0$ determinant. The paper reports numerical support from 149 G3-benchmark molecules and derives generalized versions of the virial theorem, Koopmans' theorem, and Hund's rule on the same $a$-ladder.","pith_inferences":["Editorial inference: because the paper reports that curvature grows with basis-set size, the most favorable domain for Eq. (20) is minimal or rigid basis sets; with flexible basis sets one would need an explicit second-order term proportional to $\\partial y_0/\\partial a$ to maintain accuracy.","Editorial inference: the same coupling-strength ladder suggests a practical warm-start recipe, using the converged $a=0$ orbitals as the initial guess for an $a=1$ SCF calculation; the paper mentions this in passing but does not develop it as a numerical method.","Editorial inference: the $a=0$ density is an attractive reference for building exchange-correlation approximations along the adiabatic connection, since at that endpoint both the wavefunction and density are known exactly within basis error; the paper's Hohenberg–Kohn extension points toward such a construction without formulating a specific functional."],"forward_implications":["A single $a=0$ SCF calculation, one diagonalization step with no iterative convergence, yields a parameter-free first-order estimate of the full electronic energy through the term $\\langle Y_0|H_{\\rm ee}|Y_0\\rangle$.","The extended Hohenberg–Kohn statement $\\Psi_0(a=1)\\leftrightarrow H_{\\rm ne}\\leftrightarrow Y_0(a=0)$ implies that the non-interacting density and LCAO coefficients carry the same nuclear-frame information as the physical density, allowing DFT functionals to be re-expressed around the exactly known $a=0$ endpoint.","From the generalized virial relation, the non-interacting kinetic energy $t$ exceeds the physical kinetic energy $T$, so any correction built on the TNRS must subtract a large kinetic excess when transferring to $a=1$.","Koopmans' theorem holds for arbitrary $a$ and is trivial at $a=0$, so orbital energies from the decoupled $a=0$ calculation can be transferred to $a=1$ after the Eq. (20) correction.","The fitted $L$th-order expansion of Eq. (33) reduces the average absolute deviation from G3 molecular energies from about 3.5 hartree for minimal-basis $a=1$ HF to about 1.6 hartree, showing that the TNRS transfer is a useful skeleton for empirical corrections, though not yet chemical accuracy."],"supporting_citations":[{"why":"Defines the non-interacting reference system and gives the Hohenberg–Kohn/KS formalism that the paper extends to a=0.","marker":"[3]"},{"why":"Supplies the HF-SCF formalism, Koopmans' theorem, and the proof of Hund's rule that the paper generalizes to arbitrary a.","marker":"[4]"},{"why":"Supplies the Hellmann–Feynman theorem and the standard DFT background used to derive the slope dE/da and the TNRS transfer.","marker":"[5]"},{"why":"Provide the G3 benchmark energies for 149 molecules against which the TNRS transfer and fitted expansions are compared.","marker":"[8-9]"},{"why":"Gives the exact two-electron solution z=exp(r12/2) used in Appendix 2 to illustrate the correlation factor w behind Eq. (26).","marker":"[11]"}],"fun_headline_variants":["Quasi-linear repulsion energy enables a one-step energy estimate","Non-interacting wavefunction predicts interacting energy quasi-linearly","Coulomb integral buildup is quasi-linear across coupling strength","Zero-coupling state gives physical energy via quasi-linear path","One-step energy from a zero-repulsion wavefunction: quasi-linear"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole transfer depends on the assumption that the non-interacting $a=0$ wavefunction is representative of the physical $a=1$ wavefunction, so the energy curve is nearly straight; the paper states this as the smallness of $\\langle y_0(a)|r_{12}^{-1}|\\partial y_0(a)/\\partial a\\rangle$, and notes that larger basis sets increase the curvature and weaken the approximation.","fun_headline_variants_meta":{"raw":{"variants":["Quasi-linear repulsion energy enables a one-step energy estimate","Non-interacting wavefunction predicts interacting energy quasi-linearly","Coulomb integral buildup is quasi-linear across coupling strength","Zero-coupling state gives physical energy via quasi-linear path","One-step energy from a zero-repulsion wavefunction: quasi-linear"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000441,"raw_usage":{"total_tokens":2314,"prompt_tokens":1100,"completion_tokens":1214,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":716,"completion_tokens_details":{"reasoning_tokens":1130}},"tokens_in":716,"tokens_out":1214,"duration_ms":9526,"temperature":1.0,"reasoning_tokens":1130,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:01:02.604654+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute $E(a)$ for a small atom or molecule at $a=0,0.25,0.5,0.75,1$ with a large-basis correlated method, so that basis error no longer hides curvature; if the intermediate points deviate from the straight chord between $a=0$ and $a=1$ by more than chemical tolerance, the quasi-linearity behind Eq. (20) fails. The paper's own two-electron model $z=\\exp(r_{12}/2)$ offers a place where this comparison can be made against an exact analytic solution.","supporting_citations":[],"review_version":1}