{"id":"213c1445-ff17-4574-8474-526c397b9e05","arxiv_id":"2607.22488","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Symmetry-broken Perdew-Zunger self-interaction correction reproduces the exact classical-limit energy of two-electron harmonium, with good accuracy across the full range of Planck's constant from 0 to the physical value.","lead":"This paper tests whether breaking the spherical symmetry of a density functional calculation and correcting self-interaction makes the energy exact in the classical limit of two-electron harmonium. The authors show numerically that symmetry-broken self-interaction correction approaches the exact classical energy as Planck's constant goes to zero, and that spherical averaging restores the symmetric density.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The limited Gaussian ansatz is not the critical weak point; the critical weak point is uncontrolled numerical cancellation of divergent PZ-SIC self-interaction terms in the ℏ→0 limit, which already fails for r2SCAN at ℏ=0.03.","rationale":"The reader's weakest assumption points to the limited flexibility of the two-Gaussian ansatz. I agree that this is a legitimate concern for finite ℏ accuracy, but it is not the most load-bearing issue for the exactness-in-the-limit claim: the two-Gaussian family contains the classical delta-function configuration as α→∞, R→r0, so in the ℏ→0 limit the variational ansatz is sufficient in principle. The more concrete and observable weakness is numerical: the PZ-SIC exactness relies on cancellation of Hartree and exchange-correlation self-interaction terms that individually diverge as α^{1/2}. Table I documents a spurious r2SCAN-PZ-SIC energy at ℏ=0.03 that does not converge with grid size, forcing the authors to omit that point. This is a direct indication that the numerical demonstration of the headline claim is not yet controlled for at least one of the functionals named in the conclusion. The reader's CONDITIONAL verdict already captures this risk, so no verdict change is needed; the requested code/data and higher-precision limiting analysis would settle it.","tokens_in":10844,"tokens_out":14297,"duration_ms":166056,"concrete_test":"Recompute the symmetry-broken r2SCAN-PZ-SIC energy at ℏ=0.03 and ℏ=0.01 with the Gaussian ansatz of Eqs. (10)-(11) using adaptive quadrature and/or exact analytic integrals for U[n_i], exchange, and correlation terms, refining the radial grid and increasing α optimization. If the energy converges to U_cl=0.616908 Ha to within ~1e-3, the claim is supported; if it remains roughly 0.4 Ha below the smooth trend or fails to converge, the 'any DFA' exactness demonstration fails for r2SCAN. As a secondary check, repeat with a 4-Gaussian per orbital basis; a lower energy would signal that the two-parameter ansatz also matters at the smallest calculated ℏ.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central exactness claim is that broken-symmetry PZ-SIC energy tends to U_cl as ℏ→0. For Eqs. (10)-(11), the delta-function limit is contained in the ansatz (α→∞, R→r0), so limited variational flexibility is not fatal for the limiting claim; PZ-SIC exactly cancels the divergent U[n_i] and E_xc[n_i,0] for any semilocal functional. The real load-bearing step is the numerical evaluation of this cancellation. In the localized limit, α∼1/ℏ, and U[n_i] and E_xc[n_i,0] both diverge as α^{1/2} while their difference remains finite; this requires quadrature to resolve large cancellations. Table I shows exactly this breakdown for r2SCAN-PZ-SIC: at ℏ=0.03 the energy is spuriously near zero and does not converge with radial grid, so the point was omitted (Appendix 2). Because the conclusion claims 'any density functional approximation,' the absence of a stable r2SCAN estimate in the limit leaves the broadest form of the claim unsubstantiated for one of the tested functionals, and a more flexible orbital basis would not cure a quadrature failure.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies the two-electron harmonium model with a tunable Planck constant, building on the semiclassical exact densities of Li and Li (Ref. [10]). It evaluates LSDA, PBE, and r2SCAN on those densities, first without and then with Perdew-Zunger self-interaction correction. In a symmetry-preserving (restricted) implementation, PZ-SIC substantially improves the strong-correlation energetics but leaves residual deviations at intermediate ℏ. The central part of the paper is a symmetry-broken two-Gaussian ansatz, Eqs. (10)-(11), with two variational parameters, R and α; the PZ-SIC energy is minimized for each ℏ, and the resulting energy approaches the classical minimum U_cl = 0.616908 Ha as ℏ→0 for HF, LSDA, PBE, and r2SCAN. Spherical averaging of the broken-symmetry density at ℏ=0.1 reproduces the exact spherically symmetric density. The conclusion states that symmetry-broken PZ-SIC applied to any density functional approximation is exact in the strongly-correlated limit and usefully accurate for all 0<ℏ<1.","tokens_in":11163,"tokens_out":6490,"duration_ms":70857,"significance":"If fully supported, the central claim would provide a practical route to strong correlation in DFT: symmetry-broken PZ-SIC restores the classical limit of semilocal functionals and gives a physical interpretation of the symmetric density as an orientational average. The harmonium model is exactly solvable, and the paper uses the exact densities of Ref. [10] rather than self-consistent approximate densities, which isolates functional-driven error. The comparison across three rungs of the semilocal ladder is informative and the restricted-versus-broken-symmetry contrast is clearly drawn. However, the exactness claim is stated more categorically than the evidence warrants: the numerical demonstration uses a two-parameter ansatz, and for r2SCAN the PZ-SIC evaluation fails to converge in the very limit where the claim is made. The paper is therefore valuable as a benchmark and a proof-of-principle, but the broad 'any density functional approximation' conclusion needs either a rigorous argument (e.g., from Ref. [15]) made explicit, or a numerical demonstration that is stable for all tested functionals.","major_comments":[{"comment":"The broad claim that symmetry-broken PZ-SIC applied to 'any density functional approximation' yields the exact energy in the strongly-correlated limit is not supported by the numerical evidence for r2SCAN, one of the three tested functionals. At ℏ=0.03 the r2SCAN-PZSIC energy remains anomalously near zero and does not converge as R_POINTS is increased from 900 to 2000; the point is therefore omitted as a 'numerical singularity'. This is precisely the limit in which the cancellation of divergent Hartree and XC self-interaction terms must be resolved. Without a stable evaluation scheme for meta-GGAs, the exactness claim is demonstrated only for LSDA and PBE, not for r2SCAN. Please provide a convergent r2SCAN value or restrict the conclusion.","section":"Appendix 2 and Table I"},{"comment":"The exactness statement is stronger than the variational evidence. The two-Gaussian ansatz of Eqs. (10)-(11) has only two parameters, and the authors state that it 'cannot span the full exact Hilbert space even when the variational optimization is fully converged'. The numerics show only that the optimized energy of this restricted ansatz approaches U_cl; the actual PZ-SIC minimum over a more flexible broken-symmetry orbital space could, in principle, differ. If the exactness is meant to follow from the formal argument of Ref. [15], that argument should be stated explicitly and connected to the numerical limit. If it is meant as a numerical demonstration, a basis-convergence study (e.g., multiple Gaussians per center) is needed to show that the two-parameter result is not an artifact.","section":"Section III.D and Conclusion"},{"comment":"The reference line called 'exact energy' is a fit: E(ℏ)=U_cl+(1/2)ℏ(5.0)ω, where the coefficient 5.0 is 'fitted at the endpoints' and differs from the harmonic-oscillator-approximation value 4.732. Using this fitted curve to assess the claim that broken-symmetry PZ-SIC is 'usefully accurate for all values between 0 and the physical value' is partly circular. Please report the actual reference energies from the exact construction of Ref. [10] (or the endpoint constraints used in the fit) and compare against those. The limiting exactness at ℏ=0 is unaffected, but the intermediate-ℏ accuracy claim needs a non-fitted benchmark.","section":"Eq. (8) and figures"}],"minor_comments":[{"comment":"Typos and inconsistencies: 'Plank's constant' should be 'Planck's constant'; 'the the sum' appears in Section II.B; 'DF A' is written with an inconsistent space in several places; 'PZSIC' and 'PZ-SIC' are used interchangeably.","section":"Throughout"},{"comment":"Several symbols in the exact density formula are not defined in this manuscript (B, C, c, γ, G_k). Please define them or refer explicitly to the corresponding equations in the supplemental material of Ref. [10].","section":"Eq. (2)"},{"comment":"The captions describe the black dashed line as 'the exact energy of Eq. (8)', but Eq. (8) is a fitted interpolation, not an exact expression. Re-labeling it as 'fit to reference energies' or 'semiclassical fit' would avoid misleading readers.","section":"Figures 1-3"},{"comment":"The phrase 'truncated approach' is introduced before it is defined in Appendix A. Consider moving or expanding the definition to the first occurrence.","section":"Section III.B"}],"recommendation":"major_revision","confidential_remarks":"The paper is largely a numerical demonstration of an argument already made in Ref. [15], and its incremental novelty lies in the harmonium benchmark and the density-restoration observation. That is a legitimate contribution, but the conclusion overreaches. The r2SCAN instability at ℏ=0.03 is, in my view, the central technical weakness: it strikes at the 'any DFA' claim and needs to be fixed or explicitly bracketed. The authors should also be asked to separate the formal exactness statement from the two-parameter variational demonstration. I do not recommend rejection if these points are addressed in a revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"What you should know: this paper does exactly what its title says — it demonstrates, for two-electron harmonium, that broken-symmetry PZ-SIC energies approach the classical minimum U_cl as Planck's constant tends to zero, for LSDA, PBE, and r2SCAN. The demonstration is clean: they use the recent exact or near-exact densities from Li and Li, optimize a simple two-Gaussian broken-symmetry ansatz, and then show spherical averaging of the broken-symmetry density reproduces the exact symmetric density. That last part is a nice confirmation of the standard symmetry-breaking/restoration picture, and it is genuinely new as a numerical benchmark.\n\nCredit where due: the variational parameters (R, alpha) are optimized against the PZ-SIC energy, not against the exact energy, so the central result is not fitted. The paper is also transparent about its own limitation — the r2SCAN instability at hbar=0.03 is documented in Table I and Appendix 2, with grid-convergence data showing the meta-GGA SIC evaluation fails. That is honest, reproducible reporting.\n\nNow the soft spots. The stress-test note is right: the limited Gaussian ansatz is not the critical weak point. The delta-function limit is contained in the ansatz (alpha to infinity, R to r0), and PZ-SIC exactly cancels the divergent U[n_i] + Exc[n_i,0] terms for any semilocal functional. The real load-bearing step is the numerical evaluation of that cancellation. Table I shows exactly this breakdown: at hbar=0.03, r2SCAN-PZSIC gives an energy spuriously near zero and does not converge with radial grid, so the point was omitted. The paper's broad claim — \"any density functional approximation\" — is therefore not actually demonstrated for one of the three tested functionals at the smallest hbar. The authors call it a numerical singularity and exclude it; that's defensible, but the conclusion should be softened or the quadrature fixed.\n\nSecond, the benchmark itself is partly a fit: Eq. (8) uses a fitted coefficient 5.0 rather than the analytical HOA value 4.732. The classical limit U_cl is exact, so the main conclusion doesn't depend on the fit, but the intermediate-hbar \"exact\" guide is not exact. Minor.\n\nThird, this is a numerical demonstration, not a proof. The exactness claim rests on the variational energy approaching U_cl, plus the earlier formal arguments in Ref. [15]. That's fine for a demonstration, but the conclusion states exactness as a general theorem with numerical support. A careful referee should ask for either a more rigorous limiting analysis or a more flexible basis to rule out accidental agreement.\n\nWho is this for? People working on self-interaction correction, strong correlation, and semiclassical limits of DFT. It's a useful benchmark and a clear illustration of a known principle. It deserves a serious referee — the r2SCAN-PZSIC numerical instability is a real issue that needs airing, and the paper's own appendix already invites scrutiny. I'd send it to review with a request to either fix the r2SCAN point or qualify the \"any DFA\" claim, and to release code and data.","headline":"A real numerical demonstration that broken-symmetry PZ-SIC captures the classical limit in harmonium, but the r2SCAN quadrature failure, not the Gaussian ansatz, is the main caveat.","tokens_in":867,"tokens_out":1128,"would_cite":true,"duration_ms":32385,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Symmetry-broken self-interaction correction makes density functional theory exact in the classical limit of two-electron harmonium.","keywords":["harmonium","self-interaction correction","strong correlation","classical limit","symmetry breaking","density functional theory","Perdew-Zunger","semiclassical limit"],"falsifier":"A calculation that optimizes a more flexible broken-symmetry orbital (e.g., a multi-Gaussian or numerically represented localized orbital) for harmonium at small ℏ and finds a PZ-SIC energy below the classical minimum 0.616908 Ha, or a limiting energy different from U_cl, would falsify the exactness claim.","tokens_in":10718,"feed_emoji":"⚛️","tokens_out":2838,"duration_ms":29408,"temperature":0.7,"pith_summary":"This paper tries to establish that the Perdew–Zunger self-interaction correction, applied to broken-symmetry localized orbitals, gives the exact ground-state energy of harmonium in the classical limit as Planck's constant tends to zero, and remains accurate for all values between zero and the physical value. The significance is that standard approximate density functionals such as LDA, PBE, and r2SCAN diverge or fail in this semiclassical strong-correlation limit, while self-interaction correction restores the correct finite limit. The paper also shows that spherically averaging the broken-symmetry density reproduces the exact symmetric density, supporting a time-centered interpretation of symmetry breaking.","feed_headline":"Broken-symmetry correction hits exact classical energy for harmonium","feed_subtitle":"Two-electron trap: PZ-SIC restores the ℏ→0 limit where LDA, PBE, and r2SCAN diverge.","key_machinery":"The central machinery is the Perdew–Zunger self-interaction correction formula, which removes one-electron self-interaction orbital by orbital; the variational broken-symmetry ansatz of two Gaussian orbitals located at ±R with an optimized width parameter α; and the identification of the classical limit (ℏ→0) with the strongly correlated limit, where electrons localize at the classical equilibrium separation.","core_discovery":"The authors demonstrate numerically that for two-electron harmonium, variationally optimizing a broken-symmetry pair of Gaussian orbitals centered at ±R and applying PZ-SIC yields energies that approach the exact classical minimum U_cl = 0.616908 Ha as ℏ→0, for LSDA, PBE, and r2SCAN alike. They argue this follows because PZ-SIC is exact for any collection of non-overlapping, spin-polarized one-electron densities, and in the classical limit the two electrons become localized on opposite sides of the trap with vanishing orbital overlap. The same broken-symmetry PZ-SIC energy remains close to the exact energy for all ℏ between 0 and 1, with residual deviations at intermediate ℏ attributed to th","pith_inferences":["The two-Gaussian variational ansatz has only two free parameters (R and α), and the authors admit it cannot span the full Hilbert space; a more flexible localized orbital might produce a still lower PZ-SIC energy in the ℏ→0 limit, meaning the exactness could be partly coincidental.","The mechanism suggests a testable extension: for any two-electron system with a classical analog and a Hamiltonian bounded from below, symmetry-broken PZ-SIC should reach the corresponding classical minimum; harmonium is the first demonstration, and other traps could be checked.","The r2SCAN-PZSIC numerical instability reported around ℏ=0.03 in the appendix indicates that meta-GGA self-interaction corrections may need special handling in strongly localized regimes, which could affect the generality of the 'any density functional approximation' claim.","The spherical-averaging argument implies a constructive principle: one could compute broken-symmetry localized densities from SIC and orientationally average to obtain symmetry-preserving densities, potentially bypassing the need for symmetric self-consistent calculations."],"forward_implications":["If the central claim is correct, symmetry-broken PZ-SIC provides a practical route to the exact strong-correlation limit for any density functional approximation, not just for harmonium.","The spherical-averaging result implies that symmetry-broken localized densities can be averaged to recover symmetry-preserving exact densities, offering a way to reconcile broken-symmetry calculations with exact symmetric ground states.","The demonstrated accuracy from ℏ=0 to ℏ=1 suggests that PZ-SIC could serve as a reliable correction across correlation regimes, including weakly correlated systems.","The classical-limit exactness supplies a new exact constraint that density functional approximations should satisfy, and the paper motivates developing scaled-down SIC beyond the local spin density approximation.","For two-electron singlets, the Weizsäcker kinetic energy is exact, so the remaining DFA error is dominated by exchange–correlation; SIC's success here isolates self-interaction error as the key failure mechanism in the semiclassical regime."],"fun_headline_variants":["Broken-symmetry SIC exact in classical limit for harmonium","PZ-SIC hits exact classical energy for two-electron harmonium","Classical limit: symmetry-broken PZ-SIC exact on harmonium","As ℏ→0, broken-symmetry PZ-SIC becomes exact in harmonium","Harmonium in ℏ→0: PZ-SIC achieves exact ground-state energy"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The claim relies on the two-Gaussian variational ansatz (only two optimized parameters R and α) being flexible enough to reach the true minimum of the PZ-SIC energy in the ℏ→0 limit; the authors note this ansatz cannot span the full exact Hilbert space, so the approach to the exact classical value could be coincidental.","fun_headline_variants_meta":{"raw":{"variants":["Broken-symmetry SIC exact in classical limit for harmonium","PZ-SIC hits exact classical energy for two-electron harmonium","Classical limit: symmetry-broken PZ-SIC exact on harmonium","As ℏ→0, broken-symmetry PZ-SIC becomes exact in harmonium","Harmonium in ℏ→0: PZ-SIC achieves exact ground-state energy"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000711,"raw_usage":{"total_tokens":3049,"prompt_tokens":769,"completion_tokens":2280,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":513,"completion_tokens_details":{"reasoning_tokens":2175}},"tokens_in":513,"tokens_out":2280,"duration_ms":16288,"temperature":1.0,"reasoning_tokens":2175,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T04:33:05.709314+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A calculation that optimizes a more flexible broken-symmetry orbital (e.g., a multi-Gaussian or numerically represented localized orbital) for harmonium at small ℏ and finds a PZ-SIC energy below the classical minimum 0.616908 Ha, or a limiting energy different from U_cl, would falsify the exactness claim.","supporting_citations":[],"review_version":1}