{"id":"2fd98b7e-c2c1-4ebd-9adf-f47a9cfa3736","arxiv_id":"2504.18118","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A two-parameter variational wave function reproduces the ground state of a trapped electron-positive-particle system to within 0.03% energy error and provides exact benchmarks for five e-PCP correlation functionals.","lead":"This paper introduces a toy quantum model, Exotic Harmonium, in which one electron and one positively charged particle of variable mass are trapped in harmonic oscillator wells. The model is solved with a compact variational wave function and used to benchmark electron-positive-particle correlation functionals from multicomponent density functional theory.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"0.03% total-energy accuracy does not establish 'exact' correlation benchmarks: at m=1836, ω=1 the variational correlation energy is off by ~26% from the FEM reference, so density-sensitive benchmark quantities need direct validation.","rationale":"The reader correctly identifies the KS-inversion pipeline in §5.2.2 as load-bearing for the benchmark claim, but the specific v-representability concern is likely moot: with exactly one electron and one PCP, any positive, well-behaved density can be represented as the single orbital of a non-interacting KS system (φ=√ρ), and the KS potential is then uniquely determined up to a constant. The genuinely insecure step is upstream: the input densities and pair function come from a two-parameter variational wave function, not from the near-exact FEM solution. Although total energies agree to ~0.03%, the correlation energy is a small difference of two large numbers, and for heavy PCPs the variational/FEM disagreement in this difference reaches ~26%. Because W_c, T_e^c, and T_p^c are all functionals of these approximate densities, the 'exact' E_epc values used to judge the five functionals inherit an unquantified error. The proposed test directly measures this by re-running the inversion with FEM-origin densities; if the benchmark values are stable, the concern is resolved, but if they shift substantially, the central benchmark claim requires revision. The verdict remains conditional, pending this check.","tokens_in":74616,"tokens_out":10760,"duration_ms":109315,"concrete_test":"Recompute the KS inversion of §5.2.2 for m=1836, ω=1 (and m=1, ω=1 as a control) using densities and the pair function built from the FEM radial ground-state wave function (already available for Table 4-2's RMSD entries) instead of the analytic variational densities (2-168/2-169). Then re-evaluate T_e^c, T_p^c, W_c and E_epc in Eq. (3-29). If E_epc changes by more than ~1×10⁻⁴ a.u. (≈10% of the FEM-HF correlation energy at m=1836, ω=1) or if the relative ranking or spacings of the five tested correlation functionals change, then the variational wave function does not provide exact benchmarks and the assessment in §5.2.3 needs to be re-run with FEM-based densities.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Table 4-3 shows that the variational total-energy error relative to FEM is at most ~0.03%, but the quantity used to benchmark e-PCP functionals is the correlation energy, which is a small difference between total and MC-HF energies. For heavy PCPs this error is not small: at m=1836, ω=1, the variational energy exceeds FEM by 0.000236 a.u., while the FEM-HF correlation energy is only -0.000921 a.u., so the variational correlation energy (VAR-HF) is in error by ~26%. The abstract's 'exact benchmarks' claim rests on using the variational wave function's densities and pair functions in the KS inversion of §5.2.2 and Eq. (3-29). No comparison of variational densities/pair functions to FEM-based counterparts is provided in the visible text, so the propagation of this error into T_e^c, T_p^c, W_c, and the ranking of the five functionals is unquantified. The reader's v-representability worry is not the bottleneck: with one electron and one PCP, any positive density is non-interacting v-representable via φ=√ρ, making the KS orbital and T_s unique; the insecure step is the accuracy of the input densities themselves, not their representability.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This thesis-style manuscript develops the 'Exotic Harmonium Model' (EHM), a two-component system composed of one electron and one positively charged particle (PCP) of mass m (1 ≤ m ≤ 1836) confined by a common harmonic-oscillator trap of frequency ω (10^-4 ≤ ω ≤ 10^4 a.u.) with attractive Coulomb interaction. The central methodological claim is that a two-parameter variational wave function of the form Ψ = N exp(-αr - βr² - γR²), where r and R are the relative and center-of-mass coordinates, reproduces the finite-element (FEM) ground-state total energy to within about 0.03% for all 90 mass/frequency combinations considered. The manuscript further extends the harmonium correlation toolbox — correlation hill, two-particle distribution functions, Kato cusp condition — to the e-PCP case, and uses the variational densities and pair functions in a Kohn-Sham inversion (Section 5.2.2) to assess five existing e-PCP correlation functionals, with the stated aim of providing 'exact benchmarks' for e-PCP correlation. A regression analysis (Section 5.3) is claimed to yield a compact analytical form for α(ω,m) and β(ω,m). The numerical tables provide extensive data on variational parameters, total and component energies, and MC-HF comparisons.","tokens_in":74877,"tokens_out":7583,"duration_ms":70790,"significance":"If the benchmark claim were fully substantiated, the EHM would provide a useful test bed for two-component DFT correlation functionals, complementing existing harmonium models for electron-electron correlation. The strength of the manuscript is the clean separation of center-of-mass and relative motion, the explicit analytical integrals for energy components and single-particle densities, the systematic scan over 90 systems, and the numerical validation of total energies against FEM. The variational derivation itself is standard and the reported total-energy agreement is impressive. However, the added value for the DFT community hinges on the accuracy of correlation energies and densities, not total energies; the manuscript currently does not demonstrate that accuracy, so the significance of the functional assessment remains prospective.","major_comments":[{"comment":"The abstract and Section 1.8 claim that the EHM provides 'exact benchmarks' for e-PCP correlation functionals. The benchmark quantities, however, are correlation energies defined as small differences between total energies (Eq. 3-29 and the VAR-HF/FEM-HF columns of Table 4-3). The reported total-energy accuracy of ≤0.03% does not translate into accurate correlation energies: for m=1836 and ω=1, the variational and FEM total energies differ by 0.000236 a.u., whereas the FEM-HF correlation energy is only -0.000921 a.u.; the variational correlation energy therefore deviates from the FEM-HF value by roughly 26%. Because the Kohn-Sham inversion of Section 5.2.2 uses the variational single-particle densities and the pair function (via Eq. 2-115) as input, the extracted T_e^c, T_p^c, and W_c inherit errors of this size, and the ranking of the five tested functionals is not reliably established. The manuscript provides no direct comparison of the variational densities or pair functions with FEM-based counterparts; such a validation, plus a reporting of the correlation-energy components with error bars, is required before the 'exact benchmark' claim can be sustained. The v-representability of the densities is not the limiting issue here, because with one electron and one PCP any positive density yields a unique KS orbital φ=√ρ; the relevant uncertainty is the accuracy of the variational density itself.","section":"§5.2.2, Eq. (3-29), Table 4-3"},{"comment":"The regression of the variational parameters α and β against ω and m is presented as a key deliverable ('compact yet precise analytical form'), but the provided text contains no diagnostics for the fit: the functional form, the number of fitted coefficients, the R² values, the maximum relative errors in α and β, and the error that the fitted parameters introduce into the variational energy are all absent. Without these, the claim of precision cannot be checked, and the fitted wave function cannot be used reliably by others. The fit should be reported together with a table comparing optimized and fitted parameters for all 90 systems and the resulting energy errors, including the correlation-energy error metric from the preceding comment.","section":"§5.3"},{"comment":"The text states that the Kato cusp condition (∂Ψ/∂r)|_{r=0} = -μ Ψ(0) is 'largely valid' because α≈μ. Table 4-2 shows that at low frequencies the optimized α approaches μ (e.g., α=0.999456 for m=1836 at ω=10^-4), but at high frequencies α falls to 0.4046 for m=1 at ω=10^4 (μ=0.5), a 19% deviation from the cusp value; similar deviations occur for all masses. Since the cusp condition controls the small-r behavior of the wave function and hence the pair function at short distances, the effect of this deviation on the correlation energy components and on the functional assessment should be quantified. At minimum, the claim 'largely valid' should be replaced by a quantitative statement.","section":"§3.3.7, Eq. (3-61), Table 4-2"}],"minor_comments":[{"comment":"The paragraph describing the computational methods appears twice, once with 'three methods' and once with 'four methods'; this duplication is an editing oversight and should be resolved.","section":"§4.1"},{"comment":"The text refers to 'expansion (5-1)' when it means equation (1-5); several similar cross-reference inconsistencies appear throughout the manuscript and should be corrected in a careful proofread.","section":"§1.3.2"},{"comment":"The abbreviation 'e-PCP' is used before 'PCP' is defined; define the term at first occurrence in the abstract or in Section 1.5.","section":"Abstract and §1.8"},{"comment":"The long density expression appears to have unbalanced parentheses, making it difficult to verify; check the transcription against the derivation.","section":"Eq. (2-168)"},{"comment":"The virial ratio is defined as 2⟨T⟩/⟨r·∇V⟩ in Eq. (4-7), but the text states that for HF the ratio is reported as −⟨V⟩/⟨T⟩; clarify which definition is used in Table 4-4, since the two differ for a harmonic-oscillator potential.","section":"§4.3, Eq. (4-7)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a PhD thesis in translation; its structure and length are not those of a journal article. The refereed version should be substantially condensed, with the DFT benchmark section (5.2) and the regression section (5.3) reported in full detail. The title contains a typographical artifact ('M odel') that should be corrected. The authors should also reconsider the term 'exact benchmarks' given the correlated-energy errors identified in the major comments, and target the revised version to a journal focused on density functional theory or few-body quantum systems."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The EHM paper is worth a serious look. The extension of harmonium to an attractive, unequal-mass two-component system is a real new model, and the two-parameter variational wave function is a neat construction. Across 90 systems the total energies track FEM to 0.03%, and the analytic work on densities, pair functions, and the correlation hill is carefully done. The observation that the power-series ansatz lands on the first excited state rather than the ground state for the attractive case is a genuine piece of model physics. The MC-HF comparisons are also useful and honestly reported.\n\nThe soft spot is the claim of \"exact benchmarks.\" A 0.03% error in the total energy does not translate into a small error in the correlation energy, which is what the DFT functional assessment actually uses. Your Table 4-3 shows this directly: at m=1836, ω=1, the variational total energy is 0.000236 a.u. above FEM, while the FEM–HF correlation energy is only –0.000921 a.u. So the variational correlation energy is off by roughly a quarter. Feeding those densities and pair functions into a Kohn–Sham inversion and then ranking five e-PCP functionals without checking the densities against FEM means the ranking could be significantly distorted by the variational error. That is not a fatal flaw in the model, but it is a real gap between the claim and the evidence.\n\nThe reader's v-representability worry is not the bottleneck. With one electron and one PCP, any positive density is non-interacting v-representable via φ=√ρ, so the KS inversion is well-defined in principle. The stress-test note is right: the unquantified risk is the accuracy of the input densities, not their representability.\n\nMinor issues: the duplicated paragraph in Section 4.1 should be cleaned up, and the regression in Section 5.3 should be clearly labeled as a compact interpolation, not a new result.\n\nWho gets value from this? Researchers in NEO-DFT and multi-component DFT, especially those working on electron–proton and electron–positron correlation functionals, will find the model and the variational machinery useful. It deserves peer review, but the revision should tone down \"exact,\" validate the single-particle densities and pair functions against FEM, and report how the variational error propagates into the correlation energies and potentials used to judge the functionals. I would send it out, with revision expected.","headline":"Solid variational toy model for electron–PCP correlation, but the 'exact benchmark' claim overreaches: total energies are within 0.03% of FEM while correlation energies, the quantity actually used to judge functionals, can be off by ~26%.","tokens_in":75395,"tokens_out":2255,"would_cite":false,"duration_ms":29568,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A two-parameter variational wave function reproduces ground-state energies of an electron–positive-particle harmonic trap to within 0.03%, supplying exact benchmarks for two-component density functional theory.","keywords":["Exotic Harmonium","electron-positive particle correlation","two-component density functional theory","variational wave function","correlation hill","non-Born-Oppenheimer","harmonic oscillator trap","Kohn-Sham inversion"],"falsifier":"An independent high-accuracy numerical solution of the relative-motion Schrödinger equation for the hardest cases—for instance a positron mass with $\\omega = 10^4$ or a proton mass with $\\omega = 10^{-1}$—would settle the claim: if the variational energy falls below the finite-element energy for any of the 90 systems, the optimization is invalid; if a deviation above 0.03% appears at intermediate frequencies where neither asymptotic limit applies, the two-parameter ansatz is not sufficiently flexible.","tokens_in":30,"feed_emoji":"⚛️","tokens_out":8019,"duration_ms":140254,"temperature":0.7,"pith_summary":"The paper introduces a toy model, Exotic Harmonium, in which one electron and one positively charged particle interact by Coulomb attraction while both are confined in harmonic oscillator traps with a common frequency. It argues that standard power-series methods that solve the original harmonium model only reach excited states for this attractive, unequal-mass system, so the ground state must be approached variationally. The variational wave function, built from a Slater-type and Gaussian term in the inter-particle distance times the exact center-of-mass Gaussian, matches finite-element energies to within about 0.03% for all 90 systems studied, covering ten particle masses from positron to proton and nine oscillator frequencies. From this wave function the paper derives closed-form densities, pair functions, correlation hills, mean inter-particle distances, and a fitted analytical form of the wave function in terms of oscillator frequency and particle mass. If correct, the model provides a reference laboratory for testing and improving electron–positively charged particle correlation functionals in two-component density functional theory.","feed_headline":"Two-parameter wave function solves exotic atom trap to 0.03%","feed_subtitle":"Electron-plus-positive-particle benchmark densities and correlation energies for two-component DFT.","key_machinery":"The load-bearing object is a two-parameter variational wave function that interpolates between the two limits of the model: $\\exp(-\\mu r)$ behaviour at low oscillator frequency (hydrogen-atom limit) and $\\exp(-\\frac{1}{2}\\mu\\omega r^2)$ behaviour at high frequency (harmonic-oscillator limit). Written as $\\Psi = N \\exp(-\\alpha r - \\beta r^2 - \\gamma R^2)$, it separates after a coordinate transformation into a known center-of-mass Gaussian and a relative-motion component. This single ansatz carries the whole paper: all energy components, densities, pair distribution functions, and the fitted parameter formula are obtained from its closed-form integrals.","core_discovery":"The central claim is that the ground state of the Exotic Harmonium Hamiltonian is faithfully represented by the two-parameter wave function $\\Psi = N \\exp(-\\alpha r - \\beta r^2 - \\gamma R^2)$, where $r$ is the inter-particle distance, $R$ the center-of-mass coordinate, and $\\gamma = M\\omega/2$ fixed by the total mass and trap frequency. Only $\\alpha$ and $\\beta$ need variational optimization. The variational energies agree with finite-element reference values to about 0.03% in the worst case across the entire grid of ten masses and nine frequencies, while the wave function also satisfies the Kato cusp condition in the form $\\alpha \\approx \\mu$, the reduced mass. The paper uses this accurate wave function to compute single-particle densities, two-particle distribution functions, correlation hills, and mean inter-particle distances in closed form, and then applies these quantities to benchmark five existing electron–positively charged particle correlation functionals within two-component density functional theory. A regression of the optimized parameters produces a compact analytical expression for the wave function that depends explicitly on the oscillator frequency and the positively charged particle mass.","pith_inferences":["The closed-form correlation hills and pair functions could be used as training data for a machine-learned electron–positive-particle correlation functional, something the paper does not itself pursue.","The benchmark comparison hints that any e-PCP functional depending only on the product of single-particle densities will miss the mass-ratio dependence of the correlation hill; a direct test would hold densities fixed while changing the mass ratio.","The failure of the power-series construction for attractive unequal masses suggests that exact solvability of trapped two-particle systems is tied to a sign-symmetry that is absent here, so a symmetry classification of such traps could predict when analytical ground states exist.","Extending the ansatz to excited states by multiplying with a factor $r^l$ would give a comparable benchmark for vibrationally or rotationally excited states of the model, which the paper does not address."],"forward_implications":["Each of the 90 systems provides a ready benchmark: any two-component density functional or multi-component wavefunction method can be tested against the variational energies without repeating finite-element calculations.","The fitted analytical form of $\\alpha$ and $\\beta$ as functions of oscillator frequency and positive-particle mass lets future users reconstruct an accurate ground-state wave function for this model directly, without re-optimizing parameters.","The correlation hill, the positive analogue of the correlation hole, gives a diagnostic picture of electron–positive-particle correlation that differs qualitatively from electron–electron correlation, and can guide the design of new correlation functionals.","Because the model spans the adiabatic and strongly correlated regimes as the mass ratio and frequency are varied, it can map where non-adiabatic and correlation energy contributions become comparable for exotic atoms such as positronium and muonium."],"supporting_citations":[{"why":"Supplies the quasi-analytical harmonium ground state that the paper tries to generalize and finds gives excited states for the attractive unequal-mass case.","marker":"[79]"},{"why":"Provides the power-series truncation method whose coefficients show why the analytical route fails for attractive interactions.","marker":"[87]"},{"why":"Introduces the harmonium model that Exotic Harmonium extends to two components with arbitrary mass.","marker":"[80]"},{"why":"Defines the two-particle density and correlation-hole concepts whose attractive counterpart becomes the correlation hill.","marker":"[98]"},{"why":"One of the five electron–positive-particle correlation functionals benchmarked against the Exotic Harmonium reference data.","marker":"[69]"},{"why":"Another of the tested functionals, the parameter-improved variant of the 2017 functional, used in the accuracy comparison.","marker":"[70]"},{"why":"The first electron–positron correlation functional, included among the five functionals assessed on the model.","marker":"[73]"}],"fun_headline_variants":["Two-parameter fit nails exotic atom correlation to 0.03%","Exotic Harmonium wave function hits 0.03% energy accuracy","Benchmarking electron-positive correlation with a toy model","Compact wave function unlocks electron-positive correlation energies"],"cache_read_input_tokens":77568,"weakest_assumption_plain":"The benchmark of the five correlation functionals assumes that each of the model's single-particle densities can be reproduced by a fictitious system of two non-interacting particles, and that the inversion of the Kohn-Sham equations is unique for every mass and frequency in the grid.","fun_headline_variants_meta":{"raw":{"variants":["Two-parameter fit nails exotic atom correlation to 0.03%","Exotic Harmonium wave function hits 0.03% energy accuracy","Benchmarking electron-positive correlation with a toy model","Compact wave function unlocks electron-positive correlation energies"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000427,"raw_usage":{"total_tokens":2244,"prompt_tokens":1065,"completion_tokens":1179,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":681,"completion_tokens_details":{"reasoning_tokens":1110}},"tokens_in":681,"tokens_out":1179,"duration_ms":9620,"temperature":1.0,"reasoning_tokens":1110,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T10:24:04.181981+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"An independent high-accuracy numerical solution of the relative-motion Schrödinger equation for the hardest cases—for instance a positron mass with $\\omega = 10^4$ or a proton mass with $\\omega = 10^{-1}$—would settle the claim: if the variational energy falls below the finite-element energy for any of the 90 systems, the optimization is invalid; if a deviation above 0.03% appears at intermediate frequencies where neither asymptotic limit applies, the two-parameter ansatz is not sufficiently flexible.","supporting_citations":[{"cited_title":"An effective method for generating nonadiabatic many -body wave function using explicitly correlated Gaussian-type functions,","cited_arxiv_id":null,"evidence_quote":"One of the five electron–positive-particle correlation functionals benchmarked against the Exotic Harmonium reference data."},{"cited_title":"The idea of a potential energy surface,","cited_arxiv_id":null,"evidence_quote":"Supplies the quasi-analytical harmonium ground state that the paper tries to generalize and finds gives excited states for the attractive unequal-mass case."},{"cited_title":"Reduced explicitly correlated Hartree -Fock approach within the nuclear-electronic orbital framework: Theoretical formulation,","cited_arxiv_id":null,"evidence_quote":"Provides the power-series truncation method whose coefficients show why the analytical route fails for attractive interactions."},{"cited_title":"Elimination of translational and rotational motions in nuclear orbital plus molecular orbital theory: Application of Møller-Plesset perturbation theory,","cited_arxiv_id":null,"evidence_quote":"Introduces the harmonium model that Exotic Harmonium extends to two components with arbitrary mass."},{"cited_title":"Development of a practical multicomponent density functional for electron -proton correlation to produce accurate proton densities,","cited_arxiv_id":null,"evidence_quote":"Defines the two-particle density and correlation-hole concepts whose attractive counterpart becomes the correlation hill."}],"review_version":1}