{"id":"342751fe-8cbe-4e30-918e-95a0d0326d72","arxiv_id":"1908.03799","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"For the cubic radial anharmonic oscillator, a three-parameter interpolating wavefunction gives energies accurate to roughly seven to eight digits and wavefunctions within about 1e-4 across all distances.","lead":"This paper constructs a compact three-parameter wavefunction formula that approximates the ground and low excited states of a particle in a radial anharmonic potential, with very accurate energies. A generalist may care because it offers a single uniform approximation across all coupling strengths and dimensions, useful in many areas of quantum mechanics.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 1e-4 uniform wavefunction accuracy claim rests on an unverified proxy: pointwise smallness of y1 is not shown to bound |Psi-Psi_t|/Psi_t, and no independent wavefunction comparison is reported.","rationale":"The reader's weakest_assumption identifies exactly the gap I find: the energy accuracy is independently cross-checked, but the advertised 1e-4 uniform wavefunction accuracy is only inferred from the smallness of y1 computed by the same Non-Linearization Procedure that defines the Approximant. The manuscript explicitly notes that y1 is not bounded and substitutes the smallness of y1/y0, which is not a demonstrated bound on the cumulative wavefunction error. Since the central contribution includes a locally accurate wavefunction, this is a load-bearing omission. I would not harden the verdict, because the energy claims have independent support, the Approximant encodes the correct leading asymptotics through constraints (V.14)-(V.15), the variational parameters are smooth, and the missing wavefunction comparison is straightforward to perform. The verdict should remain CONDITIONAL, with the condition being an independent wavefunction check or a real bound connecting y1 to |Psi - Psi_t|/Psi_t.","tokens_in":31950,"tokens_out":5570,"duration_ms":66065,"concrete_test":"For D in {1,2,3,6} and g in {0.1,1,10}, solve the radial Schrodinger equation directly for the ground state with an independent high-order method (e.g., Lagrange mesh with N=100 plus Richardson extrapolation, or a high-order finite-element/spectral solve). Normalize both the exact wavefunction and the Approximant (V.16) in the same way (Psi(0)=1 or unit L2 norm), using the paper's optimized parameters or re-minimized values. Compute sup_{r >= 0} |Psi_exact(r) - Psi_t(r)| / Psi_t(r). If any sup exceeds about 1e-4 for the studied cases, Eq. (V.18) fails; if all are below 1e-4, the uniform-accuracy claim is supported by an independent check.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim splits into two parts: variational energies (independently checked against the Lagrange Mesh method) and the uniform wavefunction accuracy stated as Eq. (V.18). The latter is the load-bearing unsupported claim. Section V.C says the Non-Linearization Procedure 'allows us to estimate' the deviation of the Approximant (V.16) from the exact ground state and asserts |Psi - Psi_t|/Psi_t is less than or similar to 1e-4 for all r in [0, infinity), for the D and g values studied. No bound is proved. If y = y_t + y1 + ... is the logarithmic derivative, then Psi/Psi_t = exp(-integral_r^infinity y1 ds) at first order, so a pointwise-small y1 does not by itself control the cumulative integral that enters the wavefunction ratio. The text's replacement criterion, |y1/y0| bounded and small, is heuristic; it is not connected to a sup-norm wavefunction error, and the same Non-Linearization Procedure is used both to construct (V.16) and to estimate its error. This is precisely the claim that needs an external numerical reference, and none is provided for wavefunctions.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a unified perturbative and semiclassical framework for D-dimensional radial anharmonic oscillators, based on the Riccati-Bloch (RB) equation in r-space and a Generalized Bloch (GB) equation in (gr)-space. It constructs a compact three-parameter trial wavefunction, the Approximant, by interpolating the small- and large-distance expansions of the logarithmic derivative, and fixes the parameters variationally. For the cubic anharmonic oscillator V(r)=r^2+gr^3, the authors report variational energies accurate to 7-8 significant digits for low-lying states, verified by an independent Lagrange Mesh calculation, and claim that the ground-state Approximant deviates from the exact wavefunction pointwise by at most about 10^-4 for all r in [0,∞), all studied g≥0, and all integer D≥1.","tokens_in":32208,"tokens_out":3688,"duration_ms":45275,"significance":"If the wavefunction accuracy claim is secured, this would be a genuinely useful result: a compact analytic trial function providing both highly accurate energies and a uniform approximation to the ground state of a non-solvable anharmonic problem, across coupling and dimension. The energy part of the paper is credible and well supported: the variational energies in Tables I-IV are cross-checked against the independent Lagrange Mesh method to nine decimal digits, and the second-order perturbative corrections E2 are consistently of order 10^-7 to 10^-8. The formal development of the RB and GB equations, the generating-function interpretation of the coefficients c_k^(n), and the identification of the GB expansion with a loop/semiclassical expansion are interesting and appear carefully derived. The main weakness is that the headline wavefunction claim, Eq. (V.18), is not validated by any independent wavefunction comparison and rests on a heuristic error proxy; this claim needs additional numerical support or a proof before it can be accepted.","major_comments":[{"comment":"The paper cites the boundedness of the first correction, |Y1(v)|≤Const, as the convergence criterion of the Non-Linearization Procedure (Eq. (II.13)), but for the cubic case it later states that y1 is not bounded and only the ratio |y1/y0| is bounded and small. Thus the theoretical guarantee quoted in Section II.A does not apply to the present Approximant. This does not invalidate the independently checked energies, but it removes the only theoretical rationale given for the wavefunction error estimate and reinforces the need for an independent wavefunction check.","section":"§V.C, Eq. (V.18)"}],"minor_comments":[{"comment":"The phrase '1 /greaterorsimilarr ≥ 0' appears twice in the discussion of the domain dominating the energy integrals; this looks like a corrupted comparison symbol and should be corrected to something like '1 ≲ r' or 'r≳1'.","section":"§V.C, text after Table I"},{"comment":"The table layout is inconsistent: for D=1 only E_0^(1) is shown, while for D=2 the columns include E_0^(1), -E2, and E_0^(2). Please use the same column structure for all dimensions or add a note explaining why the D=1 corrections are omitted.","section":"Table III"},{"comment":"The statement that 'all printed digits are exact' is stronger than what the comparison shows; the Lagrange Mesh method confirms the displayed digits at the nine-decimal level. I suggest replacing 'exact' by 'confirmed by independent Lagrange Mesh calculation to nine decimal digits' to avoid overclaiming.","section":"Tables I-IV, §V.C"},{"comment":"The formulas for G3 and G4 contain log[(w-1)/(w+1)] with w=(1+gr)^{1/2}; for r>0 this argument is positive and less than 1, but it may be worth stating the branch choice explicitly to avoid ambiguity.","section":"Appendix A, Eqs. (A.3)"}],"recommendation":"major_revision","confidential_remarks":"The energy results are solid and independently verified, and the formal framework is interesting. The main obstacle is the unsupported uniform wavefunction accuracy claim in Eq. (V.18), which is also stated in the abstract. I would advise the editor that this claim should be either substantiated by independent wavefunction calculations (e.g., Lagrange Mesh or high-order shooting) or explicitly downgraded to a heuristic estimate before publication. This is fixable within the manuscript's scope, so I do not recommend rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth knowing: this paper does both more and less than the abstract advertises. More, because the analytic machinery is genuinely clever and the energy numbers are good. Less, because the headline wavefunction-accuracy claim is not actually demonstrated.\n\nThe new content that holds up: the D-dimensional radial version of the GB equation, the generating-function reading of Z_n, and the three-parameter cubic Approximant. The analytic expansions in Sections I–III are careful, and the derivation of the generating functions from the recurrence relations is elegant. The variational energies in Tables I–IV are impressive: 7–8 significant digits for the low-lying states, and the Lagrange Mesh cross-check with nine stable decimal digits makes the energy claim credible. This part deserves to be published.\n\nThe soft spot is exactly where the stress-test note points. The claim that |Psi - Psi_t|/Psi_t <~ 1e-4 uniformly in r, Eq. (V.18), is inferred from the smallness of the first-order Non-Linearization correction y1. That inference is not justified as stated. The ratio Psi/Psi_t is exp(-integral y1), and pointwise smallness of y1 does not control the cumulative integral without an additional argument. The paper offers none, and there is no independent wavefunction comparison. So the 1e-4 uniform accuracy statement is plausible but unsupported. The energy accuracy is not circular, because the variational parameters are fixed by minimization and the final energies are checked against Lagrange Mesh. The wavefunction claim is circular in the sense that the same procedure used to construct the Approximant is used to estimate its error.\n\nMinor issues: the Lagrange Mesh details are deferred to a later paper, no code or data are shipped, and the Conclusions overstate things a bit (\"solution of the problem\"). These are addressable in revision.\n\nThe citation pattern is fine. The early references to the Non-Linearization Procedure and to the authors' own quartic work are appropriate; this is a new application, not a repackaging.\n\nBottom line: the energy results and the analytic framework are solid and worth a serious referee. The wavefunction accuracy claim needs either an independent numerical check (e.g., a grid solution or Lagrange Mesh wavefunction comparison) or a careful bound relating y1 to the wavefunction error. Send it to peer review, but flag that specific claim for the referee.","headline":"Solid, useful variational energies with an independent cross-check; the uniform 1e-4 wavefunction claim is plausible but rests on an unverified proxy.","tokens_in":32726,"tokens_out":1426,"would_cite":true,"duration_ms":17211,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81Q15","81Q20"],"pacs":["03.65.-w","03.65.Sq"],"model":"deepseek-v4-flash","headline":"A three-parameter trial wavefunction approximates anharmonic ground states uniformly to 1e-4.","keywords":["radial anharmonic oscillator","Riccati-Bloch equation","generalized Bloch equation","semiclassical expansion","perturbation theory","variational trial function","cubic anharmonic oscillator","uniform wavefunction approximation"],"falsifier":"Compute the exact ground state of the radial Schrödinger equation for $V = r^2 + g r^3$ by an independent high-precision method, such as a Lagrange mesh or spectral basis with many points, for example at $D = 1, 3$ and $g = 0.1, 1, 10$, and evaluate the relative deviation $|\\Psi_{\\mathrm{exact}} - \\Psi_{\\mathrm{Approximant}}|/\\Psi_{\\mathrm{Approximant}}$ on a fine grid extending to large $r$; if any deviation exceeds about $10^{-4}$, the uniform wavefunction claim fails, while if $|y_1/y_0|$ is small but the wavefunction deviation is large, the error-proxy assumption is falsified.","tokens_in":31689,"feed_emoji":"⚛️","tokens_out":6359,"duration_ms":58666,"temperature":0.7,"pith_summary":"The paper aims to build a single compact approximate wavefunction, called the Approximant, that stays close to the exact ground state of the D-dimensional radial anharmonic oscillator for every radius and every coupling strength. It works by writing the wavefunction through its logarithmic derivative and combining perturbation expansions from two complementary first-order equations, one in physical radius r and one in scaled radius gr. For the cubic case $V = r^2 + g r^3$, the resulting three-parameter trial function is claimed to reproduce the ground state to relative deviation at most about $10^{-4}$ uniformly in $r$, and to yield variational energies of low-lying states with 7 to 8 significant digits for all $g \\geq 0$ and integer $D$. If true, this gives a practical, locally accurate solution for a standard model of quantum mechanics and field theory that normally requires heavy numerical treatment.","feed_headline":"A 3-parameter wavefunction matches anharmonic states to 1e-4","feed_subtitle":"For r^2 + g r^3, one trial wavefunction yields 7-8 digit energies for every coupling and dimension.","key_machinery":"The machinery is the logarithmic derivative $y = -\\partial_r \\log \\Psi$, which turns the radial Schrödinger equation into a first-order Riccati equation. The paper introduces two rescaled versions: the Riccati-Bloch equation in the variable $v = (2M/\\hbar^2)^{1/4} r$, whose weak-coupling perturbation theory generates the small-$r$ Taylor series of $y$, and the Generalized Bloch equation in the variable $u = gr$, whose weak-coupling perturbation theory generates a semiclassical expansion at large $u$ that coincides with the loop expansion of the Euclidean path integral. Strong-coupling perturbation theory in the two equations supplies the complementary large-$r$ and small-$u$ expansions. Interpolating all four expansions yields the Approximant phase; the key identity is that the corrections $Z_n(u)$ from the Generalized Bloch equation act as generating functions for the large-$r$ coefficients of the Riccati-Bloch expansion, which is what makes the interpolation systematic.","core_discovery":"On the paper's own terms, the central discovery is that four asymptotic expansions of the logarithmic derivative $y(r)$ can be knitted into a closed analytic interpolant whose three free parameters are fixed variationally. In the weak-coupling regime, perturbation theory in the Riccati-Bloch equation reproduces the small-$r$ Taylor expansion, while the same expansion in the Generalized Bloch equation produces a new semiclassical expansion valid at large $gr$, identified with the loop expansion of the Euclidean path integral. The strong-coupling expansions fill the complementary regions. The interpolating phase, Eq. (V.13), plus two constraints, gives the Approximant Eq. (V.16); with optimal parameters, the ground state wavefunction is claimed to deviate from the exact one by less than about $10^{-4}$ uniformly for $r \\in [0, \\infty)$, and low-lying variational energies reach 7 to 8 significant digits, with second-order perturbative corrections at the $10^{-7}$ level.","pith_inferences":["If the error-proxy assumption holds, the reported $10^{-4}$ uniform wavefunction accuracy is likely conservative: the observed fast decay of higher corrections suggests the Approximant could be iterated to substantially higher precision, producing essentially exact analytic-like wavefunctions.","The identification of the Generalized Bloch expansion with the Euclidean path-integral loop expansion suggests the Approximant might be re-derived by summing a particular subset of fluctuation loops; a testable extension is to check whether the optimal variational parameters correspond to a resummation of a specific diagram class.","The paper's conjecture on square-root branch points in the strong-coupling plane could be probed numerically with the Approximant: locate the nearest branch point in the $g^{-4/5}$ variable by analytic continuation of the variational energy and compare it with the apparent convergence radius of the strong-coupling expansion."],"forward_implications":["For any $g \\geq 0$ and integer $D$, the ground state of the cubic radial anharmonic oscillator can be approximated by a three-parameter elementary function with claimed uniform relative wavefunction error below about $10^{-4}$.","Variational energies of low-lying states reach 7 to 8 significant digits, with second-order perturbative corrections near $10^{-7}$, so energy levels can be obtained essentially exactly without large basis sets.","The new semiclassical expansion from the Generalized Bloch equation provides a compact route to higher-order WKB and loop corrections for radial anharmonic potentials.","The Non-Linearization Procedure around the Approximant converges rapidly, with successive corrections dropping by a factor near $10^{-2}$, so the Approximant can serve as a zeroth order for systematic perturbative improvement of other observables.","The same interpolation scheme is announced to extend to quartic and sextic radial anharmonic oscillators, indicating a family of accurate analytic-looking approximations for polynomial radial potentials."],"supporting_citations":[{"why":"Supplies the earlier successful construction of a uniform trial phase for the quartic anharmonic oscillator, the model this paper extends to the radial cubic case.","marker":"[3]"},{"why":"Continues the quartic-oscillator program and demonstrates the high-accuracy energy and wavefunction results that motivate the Approximant approach.","marker":"[4]"},{"why":"Defines the Non-Linearization Procedure used here to generate perturbative corrections for energy and logarithmic derivative around the Approximant.","marker":"[21]"},{"why":"Introduces the one-dimensional generalized Bloch equation whose D-dimensional extension is a central tool of the paper.","marker":"[22]"},{"why":"Connects the GB-equation expansion to the semiclassical WKB and Euclidean path-integral loop expansion, grounding the claim that the new expansion coincides with the loop expansion.","marker":"[23]"},{"why":"Further develops the GB-equation perturbation scheme and its interpretation as a true semiclassical expansion, which the present paper extends to radial potentials.","marker":"[24]"},{"why":"Establishes the link between variational energy and perturbation theory used to estimate and correct the accuracy of variational energies.","marker":"[38]"},{"why":"Provides the independent Lagrange-mesh energy values used to confirm the 7 to 8 significant digit claim.","marker":"[41]"}],"fun_headline_variants":["3-param trial wavefunction: 1e-4 accuracy","Cubic anharmonic: one wavefunction, 1e-4 accuracy","4 expansions, 3 parameters: 1e-4 anharmonic states","3-param wavefunction: 1e-4 accuracy, 7-8 digit energies"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The uniform wavefunction accuracy claim rests on the untested assumption that the smallness of the first perturbative correction $y_1$, computed around the Approximant, bounds the true relative error of the wavefunction; the paper does not prove such a bound or compare the Approximant against an independent wavefunction.","fun_headline_variants_meta":{"raw":{"variants":["3-param trial wavefunction: 1e-4 accuracy","Cubic anharmonic: one wavefunction, 1e-4 accuracy","4 expansions, 3 parameters: 1e-4 anharmonic states","3-param wavefunction: 1e-4 accuracy, 7-8 digit energies"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.002486,"raw_usage":{"total_tokens":9668,"prompt_tokens":1201,"completion_tokens":8467,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":817,"completion_tokens_details":{"reasoning_tokens":8382}},"tokens_in":817,"tokens_out":8467,"duration_ms":57055,"temperature":1.0,"reasoning_tokens":8382,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:02:28.402250+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the exact ground state of the radial Schrödinger equation for $V = r^2 + g r^3$ by an independent high-precision method, such as a Lagrange mesh or spectral basis with many points, for example at $D = 1, 3$ and $g = 0.1, 1, 10$, and evaluate the relative deviation $|\\Psi_{\\mathrm{exact}} - \\Psi_{\\mathrm{Approximant}}|/\\Psi_{\\mathrm{Approximant}}$ on a fine grid extending to large $r$; if any deviation exceeds about $10^{-4}$, the uniform wavefunction claim fails, while if $|y_1/y_0|$ is small but the wavefunction deviation is large, the error-proxy assumption is falsified.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the earlier successful construction of a uniform trial phase for the quartic anharmonic oscillator, the model this paper extends to the radial cubic case."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Continues the quartic-oscillator program and demonstrates the high-accuracy energy and wavefunction results that motivate the Approximant approach."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the Non-Linearization Procedure used here to generate perturbative corrections for energy and logarithmic derivative around the Approximant."},{"cited_title":"Ari and M","cited_arxiv_id":null,"evidence_quote":"Introduces the one-dimensional generalized Bloch equation whose D-dimensional extension is a central tool of the paper."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Connects the GB-equation expansion to the semiclassical WKB and Euclidean path-integral loop expansion, grounding the claim that the new expansion coincides with the loop expansion."},{"cited_title":"Shuryak and A","cited_arxiv_id":null,"evidence_quote":"Further develops the GB-equation perturbation scheme and its interpretation as a true semiclassical expansion, which the present paper extends to radial potentials."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the link between variational energy and perturbation theory used to estimate and correct the accuracy of variational energies."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the independent Lagrange-mesh energy values used to confirm the 7 to 8 significant digit claim."}],"review_version":1}