{"id":"a5ce44c4-f005-4801-813e-c1a751d914c2","arxiv_id":"2607.07648","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":8,"one_line_summary":"Solving the GP93 inverse problem yields a semilocal exchange factor with hydrogenic exactness and a -1/r tail that, when switched by kinetic-energy and Laplacian indicators, recovers Rydberg-like bound states in s-shell atoms where PBE and SCAN bind none.","lead":"This paper solves a 1993 equation by Gill and Pople to build a semilocal exchange functional that gives the hydrogen atom its exact exchange potential, then switches it on only in one-electron-like regions to recover Rydberg bound states that PBE and SCAN miss entirely. A generalist might read it because it shows a purely semilocal route to the -1/r exchange tail that standard DFT functionals lack, with potential applications to surface physics and excited states.","discovery_kind":"unclear","skeptic_critique":{"model":"glm-5.2","headline":"Rydberg binding in He relies on the bare F_mix (alpha-switch only, no q-switch) in a regime where the paper's own analysis shows the alpha-switch is insufficient for molecular stability — but the atomic He case is never tested with the q-switch active, leaving open whether q-switch interference or a","rationale":"The reader correctly identified the load-bearing concern: the molecular/solid-state extension via the q-switch is demonstrated only non-self-consistently, and the self-consistent implementation is deferred. I agree this is the weakest link. However, I would sharpen the concern slightly: the issue is not only about molecules — it is that the central atomic result (He Rydberg states) is obtained with the alpha-only form, and the q-switch that is supposed to extend the framework without damaging atomic tails has never been tested self-consistently even on atoms. The paper argues analytically (Eq. 17, 21) that q→+∞ in atomic tails makes D(q)→1, so the q-switch is saturated and harmless there. This argument is sound in principle, but the non-self-consistent verification (Sec. VI.D, Table V) covers molecules, not the self-consistent atomic Rydberg states themselves. The verdict should remain CONDITIONAL. The mathematical core (GP93 solution, Theorems 1-2, Proposition 3) is sound and independently checkable. The hydrogenic cation predictions (He+, Li2+) are genuine independent checks. The qualitative Rydberg binding in s-shell atoms is well-supported with appropriate caveats. But the paper's own honesty about limitations (H2 dissociation failure, basis dependence, p_c=5 ad hoc damping, no code/data release) supports the CONDITIONAL verdict. The paper would move toward ACCEPT with: (1) a self-consistent real-space implementation of the triple switch on at least one atom (He) and one molecule (H2), showing the Rydberg count is preserved and the bond-center divergence is suppressed; (2) code/data release. The 8-parameter count (p, s_0, w_alpha, q_c, w_q, p_c, a, b) is reasonable for a functional with this scope, though p_c=5 is acknowledged as ad hoc. No internal inconsistency was found; the concern is about unverified self-consistency, not incorrectness.","tokens_in":18429,"tokens_out":1033,"duration_ms":843631,"concrete_test":"Implement the triple-switched form (Eq. 18) self-consistently on a real-space grid for the He atom and recompute the Rydberg bound-state count and -epsilon_HOMO. If the count drops below 8 or -epsilon_HOMO shifts by more than ~1 eV relative to the alpha-only result (25.12 eV, 8 states), the claim that the q-switch preserves atomic tails under self-consistency fails. This is the minimal test that bridges the non-self-consistent discrimination evidence to the central atomic result.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central empirical result — 8 bound Rydberg states in He — is obtained with the alpha-switched form of Eq. (16), not the triple-switched form of Eq. (18). The q-switch is introduced specifically because alpha alone is insufficient to suppress the divergent F_mix branch in low-density, large-s, alpha≈0 regions (Sec. IV.A). The paper's own Table V shows He2+ at the bond center has switch=0.114 (borderline). But more critically: in the He atom tail, q→+∞ (Eq. 17), so D(q)→1 and the q-switch is saturated — meaning the triple switch should not perturb the He result. This is argued but never numerically verified self-consistently. The He Rydberg states are computed with the alpha-only form (Eq. 16), and the claim that the q-switch 'does not damage the atomic -1/r tail at all' (Eq. 21) is verified only non-self-consistently on fixed densities (Sec. VI.D). If a self-consistent q-switched calculation on He were to show fewer bound states or a shifted HOMO, the headline result would weaken. The paper acknowledges this gap (Sec. VIII: 'the self-consistent implementation of the triple switch belongs to real-space methods') but the reader's verdict correctly flags it. The deeper concern: the 8 bound states in He depend on the bare F_mix divergent branch acting in the tail where alpha≈0, but the same divergent branch causes SCF divergence in H2 dissociation and molecular bond centers. The discrimination between 'tail where divergence is good' and 'bond center where divergence is catastrophic' rests entirely on the sign of q, demonstrated only on fixed densities. The load-bearing question is whether this discrimination survives self-consistency.","agreement_with_reader":"agree"},"referee_report":{"model":"glm-5.2","summary":"This manuscript solves the Gill-Pople (GP93) ordinary differential equation for a GGA exchange functional that reproduces the exact exchange potential of the hydrogen atom, obtaining the enhancement factor $w(z)$ in closed form via variation of constants (Eq. 7) with Frobenius matching at $z=0$ (Eq. 10). The resulting GP93 factor is parameter-free, exhibits the hydrogen-exact large-gradient growth $O[s(¥ln s)^{2/3}]$ (Proposition 3), and is exact on hydrogenic 1s ions by coordinate scaling (Theorem 2). The factor is embedded in a switched meta-GGA form (Eq. 16) using the indicator $¥alpha=(¥tau-¥tau_W)/¥tau_{¥rm unif}$ to confine the divergent branch to one-electron-like regions. The central empirical result is that this functional produces Rydberg-like series of bound virtual Kohn-Sham states in all five tested s-shell systems (H, He, Li, Na, K) where PBE through SCAN bind none in He, while p-shell atoms (Ne, Ar) are correctly left unmodified. A triple switch adding the reduced Laplacian $q$ (Eq. 18) is introduced to suppress divergence at covalent bond centers, and an orbital-free replacement for $¥alpha$ via a Perdew-Constantin-type kinetic-energy-density functional (Eq. 28) is proposed to eliminate the need for a generalized Kohn-Sham solver.","tokens_in":18623,"tokens_out":1342,"duration_ms":256773,"significance":"The paper ships a genuinely parameter-free, analytically derived exchange ingredient (the GP93 factor) with a rigorous hydrogen-exactness theorem (Theorem 1) and a scale-covariance eigenvalue theorem (Theorem 2) verified numerically on He+ and Li2+ not used in calibration (Sec. VIII: 54.44 vs 54.42 eV). The asymptotic growth $O[s(¥ln s)^{2/3}]$ (Proposition 3, Eq. 12) is a concrete, falsifiable prediction distinct from the AK13 $s¥ln s$ form. The recovery of 4-8 bound Rydberg-like states across five s-shell systems (Tables II-III), with zero bound states in p-shell atoms by design, is a clear qualitative advance over existing semilocal functionals. The Frobenius matching is stable to eight digits (Appendix B), and implementation verification includes derivative checks to $10^{-9}$ and energy-minimum confirmation (Appendix C). The scope is honestly delimited: no thermochemistry claims, no p-shell HOMO improvement, no response properties.","major_comments":[{"comment":"§VI.B, Table II, and the He Rydberg result: The headline result of 8 bound Rydberg states in He is obtained with the alpha-switched form of Eq. (16), not the triple-switched form of Eq. (18). The q-switch is introduced (Sec. IV.A) because alpha alone is insufficient to suppress the divergent F_mix branch in low-density, large-s, alpha≈0 regions such as covalent bond centers. The manuscript argues (Eq. 21) that the q-switch should not perturb the He atomic tail because q→+∞ there, making D(q)→1. However, this argument is verified only non-self-consistently on fixed densities (Sec. VI.D). A self-consistent q-switched calculation on He—confirming that the bound-state count and HOMO are preserved when the full triple switch is active—would close the logical gap between the functional form advocated for molecules/solids (Eq. 18) and the functional form used to produce the central atomic Ryd伯格","section":null}],"minor_comments":[{"comment":"§V, Computational Methods: The basis-set dependence of the Rydberg state count is acknowledged but not quantified. A convergence study showing how the count of bound states in He changes from aug-cc-pVDZ through aug-cc-pV5Z would strengthen the claim that the qualitative contrast (zero vs. eight) is basis-robust.","section":null},{"comment":"Table I: Only H, He+, and Li2+ are listed. The text mentions Z=1-5 in Sec. VIII, but Table I shows only three systems. Adding Z=4,5 would complete the demonstration of Theorem 2.","section":null},{"comment":"§VII, Eq. (28): The large-p damping factor with p_c=5 is acknowledged as 'not a parameter-free consequence of the GP93 construction.' This is an additional empirical parameter that should be listed alongside the switching parameters in the parameter inventory for full transparency.","section":null},{"comment":"Fig. 1 caption: The y-axis label 'Mean absolute error' should specify units (eV) for consistency with Eq. (14).","section":null},{"comment":"§III, Eq. (13): The outer switch parameters (p, s_0) are described as 'calibrated' rather than 'empirical.' While the distinction is explained, the phrasing may invite confusion; a brief note that 'calibrated' means 'determined by the functional's design purpose' would help.","section":null},{"comment":"Reference [3] (Gill and Pople 1993) is central to the entire construction but is cited only in the introduction and Sec. II. A brief note on how the present solution relates to GP93's original (incomplete) treatment would provide useful historical context.","section":null}],"recommendation":"major_revision","confidential_remarks":"The core mathematical derivation (Sec. II) is clean and the GP93 factor itself is genuinely parameter-free. The central concern is the gap between the functional form used for the headline result (alpha-switched, Eq. 16) and the form advocated for general use (triple-switched, Eq. 18). The author acknowledges this gap (Sec. VIII) but the He Rydberg states are the paper's strongest selling point and they are computed with the simpler form. A single self-consistent q-switched He calculation would either confirm robustness or reveal a problem—either outcome substantially strengthens the paper. The orbital-free alpha_PC (Eq. 28) with its ad hoc p_c=5 damping is a secondary concern; it is presented as a pragmatic regularization rather than a fundamental result, which is appropriate. The paper is within scope for a computational physics journal if the major comment is addressed."},"author_rebuttal":{"model":"glm-5.2","summary":"The referee identifies one major issue: the headline He Rydberg result (8 bound states) is obtained with the alpha-switched form (Eq. 16), while the triple-switched form (Eq. 18) including the q-switch is advocated for molecules/solids but verified only non-self-consistently. The referee requests a self-consistent q-switched calculation on He to confirm the bound-state count and HOMO are preserved under the full triple switch. We agree this is a logical gap and will address it by performing the requested calculation or, where self-consistent implementation in Gaussian bases is not feasible, by providing a detailed non-self-consistent verification on the self-consistent alpha-switched density and revising the manuscript to state the limitation explicitly.","responses":[{"response":"The referee is correct that the headline He Rydberg result (Table II, 8 bound states) is obtained with the alpha-switched form of Eq. (16), and that the q-switch of Eq. (18) is verified only non-self-consistently. We agree that a self-consistent q-switched calculation on He would close the logical gap between the functional form advocated for general use and the one used to produce the central atomic result. We will address this in revision. The key question is whether the q-switch, when active self-consistently, preserves the He bound-state count and HOMO. Our analytical argument (Eq. 21) shows that in the He atomic tail q→+∞, so D(q)→1 and the q-switch is saturated; the derivative ∂D/∂q_c vanishes there, meaning the q-switch is inert in the region that determines the Rydberg levels. The referee's concern is whether this holds self-consistently, not just on fixed densities. We will attempt the self-consistent q-switched calculation on He. However, we note that the v_lapl term associated with q is numerically unstable in Gaussian bases (as discussed in Sec. IV.A and Sec. VIII), which is why the q-switch was demonstrated non-self-consistently in the original submission. If a stable self-consistent Gaussian implementation proves infeasible, we will: (1) provide a non-self-consistent evaluation of the full triple-switched functional on the self-consistent alpha-switched He density, showing that the exchange potential tail and bound-state spectrum are unchanged when D(q) is applied; (2) revise the manuscript to state explicitly that the self-consistent q-switched verification on He is a necessary next step requiring real-space implementation, and that the current evidence supports—but does not definitively confirm—that the triple switch preserves the atomic Rydberg result.","revision_made":"partial","referee_comment":"§VI.B, Table II, and the He Rydberg result: The headline result of 8 bound Rydberg states in He is obtained with the alpha-switched form of Eq. (16), not the triple-switched form of Eq. (18). The q-switch is introduced (Sec. IV.A) because alpha alone is insufficient to suppress the divergent F_mix branch in low-density, large-s, alpha≈0 regions such as covalent bond centers. The manuscript argues (Eq. 21) that the q-switch should not perturb the He atomic tail because q→+∞ there, making D(q)→1. However, this argument is verified only non-self-consistently on fixed densities (Sec. VI.D). A self-consistent q-switched calculation on He—confirming that the bound-state count and HOMO are preserved when the full triple switch is active—would close the logical gap between the functional form advocated for molecules/solids (Eq. 18) and the functional form used to produce the central atomic Ryd伯格"}],"tokens_in":18140,"tokens_out":932,"duration_ms":177677,"standing_objections":["A fully self-consistent q-switched calculation on He may not be achievable within the Gaussian-basis framework used in this work, because the v_lapl term (functional derivative of the exchange energy with respect to ∇²n) is numerically unstable in Gaussian bases. This is the same difficulty faced by Laplacian-level meta-GGAs and is the reason the q-switch was demonstrated non-self-consistently in the original submission. A definitive self-consistent verification likely requires a real-space grid implementation, which is beyond the scope of the current manuscript. If the self-consistent calculation cannot be stabilized, the logical gap identified by the referee can be narrowed but not fully closed in this revision."]},"desk_editor":{"model":"glm-5.2","letter":"The headline: Imoto solves the Gill-Pople (1993) ODE as an inverse problem and uses the resulting enhancement factor — the GP93 factor — as a nonempirical ingredient in a switched meta-GGA that produces bound Rydberg states in H, He, Li, Na, and K, where PBE through SCAN bind none in He. That is a real result. The GP93 factor is parameter-free, derived from the hydrogen-exactness condition, and the large-gradient growth O[s(ln s)^{2/3}] is distinct from AK13's s*ln(s). Theorems 1 and 2 (hydrogen exactness and hydrogenic eigenvalue scaling) are correctly stated and numerically verified on He+ and Li2+ — systems not used in calibration. The Frobenius matching at z=0 is stable to eight digits. The paper is unusually honest about scope: no thermochemistry claims, p-shell atoms excluded by design, H2 dissociation failure disclosed. That honesty earns credit. The soft spots are real but mostly where the paper admits them. The q-switch (reduced Laplacian) is demonstrated only non-self-consistently on fixed densities for molecules; the self-consistent version is unstable in Gaussian bases and deferred to real-space implementations. The orbital-free alpha_PC requires an ad hoc large-p damping factor (p_c=5) that the author explicitly calls 'not a parameter-free consequence of the GP93 construction.' No code or data is shipped. The HOMO-IP MAE of 0.80 eV is calibrated on the same five systems it is reported on, though the Rydberg binding is shown robust over s_0 = 0.13–0.28 and the hydrogenic cation checks are genuinely independent. On the stress-test concern about He Rydberg states being computed with the alpha-only form (Eq. 16) rather than the triple switch (Eq. 18): the paper's argument that q→+∞ in the atomic tail saturates D(q)→1, so the triple switch reduces to the alpha-only form, is analytically sound (Eq. 17). The concern that this hasn't been verified self-consistently is valid but secondary — the analytic argument is clean and the fixed-density evidence supports it. The deeper question — whether the discrimination between 'tail where divergence is good' and 'bond center where divergence is catastrophic' survives self-consistency — is the real open problem, and the paper says so. This paper is for DFT functional designers and anyone working on asymptotic potential correction in semilocal DFT. The core mathematical contribution (GP93 solution + theorems) is solid and new. The atomic Rydberg result is well-supported. The molecular extension is preliminary and correctly labeled as such. It deserves a serious referee who can check the derivation and assess whether the switching framework is viable enough for the community to pursue the real-space implementation.","headline":"First semilocal exchange functional to recover the -1/r asymptotic tail and bound Rydberg states in s-shell atoms, via a complete solution of the 1993 Gill-Pople equation.","tokens_in":19541,"tokens_out":686,"would_cite":true,"duration_ms":80537,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["31.15.E-","71.15.Mb"],"model":"glm-5.2","headline":"Hydrogen-exact exchange factor binds Rydberg states in semilocal DFT","keywords":[],"falsifier":"A self-consistent real-space implementation of the triple-switched functional that fails to converge, or that converges but does not reproduce the bound Rydberg states and HOMO-ionization-energy improvements shown in the fixed-density and Gaussian-basis calculations, would falsify the central claim that the GP93 factor can be embedded in a stable, practical semilocal functional for general systems.","tokens_in":18493,"feed_emoji":"⚛️","tokens_out":1283,"duration_ms":233840,"temperature":0.7,"pith_summary":"Semilocal density functionals such as PBE and SCAN produce exchange potentials that decay exponentially outside atoms and molecules, lacking the correct Coulombic tail of -1/r. This structural defect means no existing semilocal functional supports bound Rydberg-like virtual states in atoms like helium, where PBE through SCAN bind zero. The paper revisits a 1993 ordinary differential equation derived by Gill and Pople, which specifies the condition a gradient-corrected exchange functional must satisfy to reproduce the exact exchange potential on the hydrogen atom. The author solves this equation completely as an inverse problem, obtaining a unique enhancement factor called the GP93 factor. This factor is exact on the hydrogen 1s density by construction and, through uniform coordinate scaling, yields the exact eigenvalue -Z^2/2 for every hydrogenic 1s ion. Its large-gradient growth follows the form O[s(ln s)^{2/3}], which is the hydrogen-exact structure and is distinct from the s*ln(s) form used in prior attempts at improved asymptotics. The GP93 factor contains no empirical parameters. To prevent the unbounded enhancement factor from destabilizing many-electron systems, the author embeds it in a switched meta-GGA form. A kinetic-energy-density indicator alpha detects one-electron-like regions (where alpha equals zero) and activates the GP93 term only there, while reverting to PBE in many-electron regions. A reduced-Laplacian indicator q further distinguishes atomic tails (where q is large and positive) from covalent bond centers (where q is small or negative), suppressing the divergent branch at bond centers without disturbing the atomic tail. The resulting functional produces Rydberg-like series of bound virtual Kohn-Sham states in all five tested systems whose outermost shell is a one-electron-like s shell: H, He, and Li, Na, K with large-core pseudopotentials. For p-shell atoms Ne and Ar, the switch closes by design and no bound states appear. The HOMO eigenvalue moves toward the experimental ionization energy with a mean absolute error of 0.79 eV across the five s-shell systems, compared to 3.3-4.2 eV for PBE, SCAN, and B3LYP. The author also shows that the orbital-dependent kinetic-energy-density indicator can be replaced by an orbital-free kinetic-energy-density functional with a von Weizsacker lower bound, rendering the entire switch a pure density and","feed_headline":"Hydrogen-exact exchange factor binds Rydberg states in semilocal DFT","feed_subtitle":"Solving a 1993 equation as an inverse problem yields a nonempirical exchange ingredient that recovers bound Rydberg-like states where PBE to","key_machinery":"GP93 factor (solution of the Gill-Pople equation as inverse problem); alpha indicator (kinetic-energy-density ratio detecting one-electron character); q indicator (reduced Laplacian distinguishing atomic tails from covalent bond centers); triple-switched enhancement factor combining PBE with the GP93 factor; orbital-free alpha_PC replacing the orbital-dependent kinetic-energy density with a Perdew-Constantin-type functional","core_discovery":"The central object is the GP93 factor, defined by solving the Gill-Pople equation as an inverse problem. It is an enhancement factor for semilocal exchange that is exact on the hydrogen 1s density by construction, exact on all hydrogenic 1s ions by coordinate scaling, and has the hydrogen-exact large-gradient growth O[s(ln s)^{2/3}]. When embedded in a switched functional that activates the factor only in one-electron-like regions detected by the alpha indicator and further refined by the reduced Laplacian q, the functional produces bound Rydberg-like virtual Kohn-Sham states in all five tested s-shell systems (H, He, Li, Na, K) where PBE through SCAN bind none in He. The domain of applicabl","pith_inferences":[],"forward_implications":["A self-consistent real-space pseudopotential implementation of the triple-switched functional could recover Rydberg-like bound states and improved ionization energies at semilocal computational cost, without nonlocal exact exchange.","The alpha-switch design principle predicts that only systems with one-electron-like s-shell HOMOs benefit from the asymptotic correction, while p-shell atoms and many-electron regions revert to PBE, fixing the domain of applicability by construction.","Large-core pseudopotentials play a constructive role by rendering alkali valence electrons one-electron-like, suggesting a natural affinity between this functional and real-space pseudopotential codes.","The orbital-free replacement of the kinetic-energy-density indicator removes the need for a generalized Kohn-Sham solver, making the entire switch a pure density functional with a standard local exchange potential.","The -1/r correction is switched off in smooth bulk interiors and acts in low-density one-electron-like regions such as the vacuum side of surfaces, potentially affecting adsorption energetics where image-potential-like asymptotes matter."],"fun_headline_variants":["GP93 exchange factor recovers bound Rydberg states in semilocal DFT","Hydrogen-exact exchange ingredient binds Rydberg-like states in s-shell atoms","Switched semilocal functional with GP93 factor yields Rydberg series in H to K","Inverse-problem solution to Gill-Pople equation restores Rydberg states in DFT","Hydrogenic-exact exchange factor binds virtual Rydberg states where SCAN cannot"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The alpha-switch and q-switch, with a few calibrated parameters, are assumed to cleanly separate one-electron-like atomic tails from many-electron regions and covalent bond centers. This separation is demonstrated only non-self-consistently on fixed densities for molecules, and the self-consistent q-switch is unstable in Gaussian bases. The entire molecular and solid-state extension rests on the unverified premise that a real-space self-consistent implementation will converge","fun_headline_variants_meta":{"raw":{"variants":["GP93 exchange factor recovers bound Rydberg states in semilocal DFT","Hydrogen-exact exchange ingredient binds Rydberg-like states in s-shell atoms","Switched semilocal functional with GP93 factor yields Rydberg series in H to K","Inverse-problem solution to Gill-Pople equation restores Rydberg states in DFT","Hydrogenic-exact exchange factor binds virtual Rydberg states where SCAN cannot"]},"model":"glm-5.2","effort":"high","cost_usd":0.0,"raw_usage":{"total_tokens":912,"prompt_tokens":801,"completion_tokens":111,"prompt_tokens_details":null},"tokens_in":801,"tokens_out":111,"duration_ms":109021,"temperature":1.0,"reasoning_tokens":null,"cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-09T03:43:02.904501+00:00","model_set":{"reader":"glm-5.2"},"falsifier":"A self-consistent real-space implementation of the triple-switched functional that fails to converge, or that converges but does not reproduce the bound Rydberg states and HOMO-ionization-energy improvements shown in the fixed-density and Gaussian-basis calculations, would falsify the central claim that the GP93 factor can be embedded in a stable, practical semilocal functional for general systems.","supporting_citations":[],"review_version":1}