{"id":"78fa56ae-aaf8-4e83-ac07-19a7aacd3d94","arxiv_id":"2607.19197","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Approximate analytical energies for the radial screened Coulomb potential, accurate to ~0.06–0.6% at c=0.1, extended to positronium.","lead":"This paper derives closed-form approximate formulas for the energy levels of atoms in a radially screened Coulomb potential, where screening acts near the nucleus. The best formulas agree with high-precision numerical results to about a tenth of a percent at weak screening.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Variational excited-state energies rely on an unproven one-parameter ansatz; an orthogonalized re-minimization would settle whether the claimed improvement is robust.","rationale":"The strongest claim concerns the accuracy of the analytical energies, and the derivations are internally consistent; I found no algebraic red flag in the Coulomb, Kratzer, or unified expectation-value formulas. The excited-state variational procedure is the only place where a false positive could enter: because the β-minimization is not constrained by orthogonality, the reported improvement for 2s–10s could in principle be an artifact of the ansatz rather than a genuine approach to the exact eigenstates. A concrete orthogonalized re-minimization would settle this. The bias-direction contradiction in §2.3 is real and should be corrected, but it does not change the numerical energy values, so the reader's CONDITIONAL verdict stands unchanged.","tokens_in":23156,"tokens_out":26531,"duration_ms":260634,"concrete_test":"Recompute the 2s–10s variational energies after Gram-Schmidt orthogonalizing each scaled-Kratzer trial against the optimized lower-state trial functions and re-minimizing β; if any energy shifts by more than ~0.1% (a significant fraction of the claimed 0.06–0.40% errors), the excited-state variational improvement is not robust. As a secondary check, collect E_Kr − E_GPS for all Tables 1–2 entries to confirm whether the Kratzer bias is indeed negative everywhere, contradicting §2.3.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—sub-0.63% Kratzer and sub-0.41% variational energies at c=0.1—is numerically credible: I re-checked Eqs. (14), (24), and (A.18)–(A.20) against the master integral (13) and the 2s reduction, and the algebra is consistent. The load-bearing weakness is the excited-state variational step. Eq. (28) introduces a single scaling parameter β, and β_opt is found by minimizing E_var(β) separately for each state. For n>1 the trial functions are not orthogonalized against lower states, so the variational upper-bound theorem does not apply; nothing in the derivation prevents β_opt from being pulled toward a lower-energy subspace, and the claim that the variational method 'provides the best accuracy' for the first ten s-states is therefore an empirical observation, not a theorem. This matters because the reported improvement over the Kratzer reference is small (0.06–0.40%), exactly the scale at which non-orthogonal contamination could masquerade as improvement. Separately, §2.3's assertion that the Kratzer reference overestimates binding is contradicted by Tables 1–2, where E_Kr > E_GPS (less bound) for every listed state; this undermines the 'complementary biases' error-estimation claim, though not the energy formulas themselves.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops three analytical approximation schemes for bound states of the radial screened Coulomb potential (RSCP), V(r) = -(1/r)e^{-c/r}, for hydrogen-like atoms: (i) expectation values evaluated with Coulomb and Kratzer reference eigenfunctions, (ii) a one-parameter variational optimization using a scaled Kratzer basis, and (iii) a Hellmann–Feynman route obtained by integrating the derivative of the expectation value. Closed-form expressions in terms of modified Bessel functions are derived for arbitrary n and l, with a unified master formula in Appendix A. The results are benchmarked against GPS data at c=0.1, reporting relative errors below 0.63% for the Kratzer expectation-value method and below 0.41% for the variational method for the first ten s-states; the extension to Positronium is presented in Section 7.","tokens_in":23546,"tokens_out":6334,"duration_ms":60286,"significance":"The RSCP is not exactly solvable analytically, and the paper provides a transparent, self-contained set of closed-form approximations that could be useful for plasma diagnostics and parametric spectral studies. The algebraic core appears sound: I re-derived the master integral (13), the general expressions (14), (15), (A.18)–(A.20), and the c→0 limits, and they are internally consistent. A strength is that no parameter is fitted to the GPS benchmark; β_opt is chosen by minimizing the expectation value of the exact RSCP Hamiltonian. However, the manuscript makes two load-bearing interpretive claims that are not supported by its own evidence: the excited-state variational calculation lacks a variational upper bound and the stated 'complementary biases' of the two reference methods are contradicted by Tables 1–2. These issues do not invalidate the energy formulas themselves but they do affect the significance of the claimed improvements and the error-estimation rationale.","major_comments":[{"comment":"The variational calculation for excited states is not protected by the variational theorem. The trial functions (28) for different n_r are not orthogonalized against lower states, so minimizing E_var(β) for n_r>0 can lower the energy artificially by admixing lower-state character. The manuscript acknowledges this in the text after Eq. (28), yet Section 5.2 and the abstract present the variational results as providing the 'best accuracy' for the first ten s-states. Since the reported improvement over the Kratzer reference is small (roughly 0.06–0.40% in Table 1), this is exactly the scale at which non-orthogonality contamination could masquerade as improvement. A concrete test would be to Gram–Schmidt orthogonalize each excited trial function against the optimized lower states and re-minimize β; if the improvement persists, the empirical claim is robust. As it stands, the improvement is a","section":"Section 3, Eqs. (28)–(34), Table 3"},{"comment":"The stated complementary bias is contradicted by the paper's own numerical results. For every listed state, the Kratzer reference energy is above E_GPS (less bound), not below it. For example, Table 1 gives E_Kr(1s) = -0.3769734 Ha versus E_GPS(1s) = -0.3793464 Ha, and E_Kr(2s) = -0.1079848 Ha versus E_GPS(2s) = -0.1083227 Ha; Table 2 shows the same sign at c=1.0 and c=10.0. Thus the Kratzer reference does not overestimate binding; it underestimates binding, and the statement in Section 3 that the true RSCP energy lies between the Coulomb and Kratzer predictions is not supported. Consequently, the 'complementary biases ... robust error estimation' claim in the Abstract and Conclusions is not supported. This does not affect the energy formulas themselves, but it invalidates the error-estimation interpretation and should be corrected.","section":"Section 2.3, bullet 'Kratzer reference', Tables 1–2"},{"comment":"Table 4's caption states that relative errors are computed with respect to the reference values of Ref. [34], but no [34] Positronium values are listed; the table's last column is labeled S-H (2021) [32]. Moreover, the parenthetical percentages are not consistent with the S-H column if that is the reference: for the 1s variational energy, |(-0.2131429) - (-0.2134)|/0.2134 ≈ 0.120%, not 0.106% as tabulated. Without the actual [34] reference energies, the claimed sub-0.1% accuracy for Positronium cannot be verified. Please either quote the reference values in the table or add a dedicated reference column.","section":"Section 7, Table 4"}],"minor_comments":[{"comment":"Several typographical and grammatical errors should be corrected: 'RSCP. using' in the abstract, 'prvide' in the Introduction, 'numeriacl', 'significally', and 'efficeint' later in the text.","section":"Abstract and Introduction"},{"comment":"The text claims 'β_opt > 1 for all states and screening parameters', but Table 3 only reports β_opt for c=0.1. Either provide data for additional screening parameters or rephrase the claim to refer to the computed range.","section":"Section 5.3, Table 3"},{"comment":"The observation that Eq. (42) coincides with the bare Kratzer eigenvalue in Eq. (24) is important and should be stated earlier, in Section 2.2, where the Kratzer reference is introduced, not only in the comparison with Ref. [34].","section":"Section 6, Eq. (42)"},{"comment":"The legend and symbol definitions in Figure 2 are crowded and ambiguous: '1s Kratzer■Var.' and '2s◦Coulomb▲HF-Kr' mix methods and symbols without clear visual separation. Please use distinct linestyles or separate panels.","section":"Figure 2"},{"comment":"The symbol I_RSCP is introduced in Eq. (25) but defined only in Appendix A. A pointer to Eq. (A.21) would help the reader.","section":"Eq. (25)"}],"recommendation":"major_revision","confidential_remarks":"The algebraic core of the paper is sound and the benchmark numbers for hydrogen at c=0.1 are credible. The main obstacles are the unsupported excited-state variational claim (Section 3), the contradicted 'complementary biases' statement (Section 2.3), and the unverifiable Positronium reference values (Section 7, Table 4). These are fixable within the manuscript's scope, but they are load-bearing for the paper's central claims and must be addressed before acceptance."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Core take: this is a carefully done analytical approximation paper for the radial screened Coulomb potential. The derivations are self-contained and the central claim—sub-0.63% errors for the first ten s-states at c=0.1 with the Kratzer reference, and sub-0.41% with the scaled variational method—is supported by the tables. I traced the master integral, the Bessel-function double sums, and the 1s reductions; everything is consistent.\n\nWhat's new: the bare Kratzer eigenvalue is exactly Xu et al.'s asymptotic formula, and the paper says so. The actual new content is the expectation-value correction from the full RSCP, the scaled-Kratzer variational method, and the Hellmann-Feynman route. Those are honest additions, not repackaged old results. The formalism for arbitrary n and ℓ is a plus, and the Positronium extension is a nice check.\n\nSoft spots, in order of seriousness:\n\n1. The bias claim in §2.3 is simply wrong. It says the Kratzer reference overestimates binding (energies too negative), but Tables 1 and 2 show E_Kr is above E_GPS—less bound—for every listed state at c=0.1. Both references underestimate binding; the Kratzer just does so by less. The 'complementary biases' story in the abstract and Section 3 is contradicted by their own data. This doesn't affect the energy formulas, but it does undermine the error-estimation argument and needs a correction.\n\n2. The variational excited-state energies are an empirical observation, not a theorem. The paper discloses this when it says orthogonality to lower states would be needed for a strict upper bound. That's fine, but the reported improvement over the Kratzer reference is small (0.06–0.40%), exactly the scale where non-orthogonal contamination could matter. A quick orthogonalized re-minimization or a comment that the GPS benchmark confirms the improvement would settle it. The stress-test note overstates this; I don't see it as load-bearing, just a caveat.\n\n3. Table 4's caption says errors are relative to reference [34], but the listed reference column is Stachura-Hancock (2021), and the stated percentages don't match that column. That's a mislabeled caption, not a numerical fraud, but it needs fixing.\n\nMinor: a few typos throughout, and no code/data. For a pure analytical paper with closed forms, the absence of code is acceptable.\n\nOverall, the math is solid, the benchmarking is honest, and the central numerical claims hold up. The bias sign error is the most consequential issue, but it's cosmetic relative to the physics. This deserves a serious referee—someone should review it and ask for the bias section to be rewritten and the table caption corrected. I'd cite it if I worked on screened potentials, and I'd read it in a journal club focused on approximation methods.","headline":"Solid analytical approximation paper for a niche potential; the math checks out and the sub-0.63% benchmark is credible, but the bias interpretation in §2.3 contradicts their own tables and the excited-state variational claim is empirical rather than proven.","tokens_in":23980,"tokens_out":2957,"would_cite":true,"duration_ms":32158,"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":"The paper establishes that approximate energies for atoms in the radial screened Coulomb potential can be obtained in closed form—using Coulomb, Kratzer, variational, and Hellmann-Feynman methods—with relative errors below 0.63% for the fir","keywords":["radial screened Coulomb potential","bound states","Kratzer potential","effective angular momentum","modified Bessel functions","variational method","Hellmann-Feynman theorem","positronium"],"falsifier":"A converged numerical solution of the radial Schrödinger equation for the 1s state at c=0.1 can settle the claim: if the Kratzer expectation value deviates from the exact energy by more than the stated 0.63%, or if the variational energy falls below the exact ground-state energy, the central accuracy claim fails.","tokens_in":23071,"feed_emoji":"⚛️","tokens_out":6018,"duration_ms":52757,"temperature":0.7,"pith_summary":"This paper aims to show that the bound-state spectrum of hydrogen-like atoms in the radial screened Coulomb potential—a potential that is nonsingular at the origin and Coulombic at large distances—can be captured by simple analytical formulas instead of requiring a full numerical solution of the Schrödinger equation. The authors derive closed-form energy expressions in terms of modified Bessel functions by evaluating the exact potential in Coulomb and Kratzer reference states, then refine the results with a one-parameter variational scaling and with the Hellmann-Feynman theorem. Benchmarked against high-precision numerical data, the Kratzer-based formulas achieve relative errors below 0.63% for the first ten s-states at screening parameter c=0.1, and the variational version improves the ground-state error to about 0.4%. The same construction extends to positronium, demonstrating that the approach works for arbitrary reduced-mass systems. A sympathetic reader would care because these formulas make rapid, reasonably accurate predictions for plasma-embedded atoms and for spectral diagnostics without heavy computation.","feed_headline":"Screened-atom energy formulas reach 0.6 percent accuracy","feed_subtitle":"Kratzer-based shortcuts match high-precision numeric spectra and extend to positronium.","key_machinery":"The central object is the Kratzer reference Hamiltonian -1/r + c/r^2, whose eigenfunctions are Coulomb-like functions with a shifted effective angular momentum ν_l = -1/2 + sqrt((l+1/2)^2 + 2c). This absorbs the leading short-distance screening correction exactly, making the residual potential small. The master integral ∫ x^{ν-1} e^{-β/x - γx} dx = 2(β/γ)^{ν/2} K_ν(2√(βγ)) converts all expectation values into finite double sums of modified Bessel functions of the second kind. A single unified expression with an interpolation parameter k and a scaling parameter β generates the Coulomb-reference, Kratzer-reference, and variational methods as special cases.","core_discovery":"The central claim is that the energy eigenvalues of the radial screened Coulomb potential V(r) = -e^{-c/r}/r can be approximated analytically with sub-percent accuracy. The key move is to use the Kratzer potential -1/r + c/r^2 as a reference: because the leading short-distance screening term c/r^2 has the same radial form as a centrifugal barrier, it is absorbed exactly into an effective angular momentum ν_l = -1/2 + sqrt((l+1/2)^2 + 2c). Evaluating the exact RSCP Hamiltonian in Kratzer eigenstates produces closed-form energies as double sums of modified Bessel functions of the second kind, valid for arbitrary n and l. For the first ten s-states at c=0.1, the relative error against high-prec","pith_inferences":["Applying the same Kratzer-reference trick to other screened potentials, such as the Yukawa or exponential-cosine screened Coulomb potential, would test whether the effective-angular-momentum absorption generalizes; the paper does not do this.","A Gram-Schmidt orthogonalized version of the scaled-Kratzer variational basis would restore a rigorous upper-bound property for excited states and likely improve the variational energies further.","Because the Coulomb reference underbinds and the Kratzer reference overbinds at moderate c, the arithmetic mean of the two reference energies is a cheap estimator whose error could be quantified systematically against exact data.","The predicted absence of critical screening at large c for s-states is a qualitative signature that could be checked in simulations or experiments of strongly screened plasma-embedded atoms."],"forward_implications":["Energies for any bound state (n, l) and any screening parameter c are available as closed-form Bessel-function expressions, so parametric scans over screening strength need no numerical Schrödinger solver.","The effective-angular-momentum construction works uniformly across l, so states with non-zero angular momentum are approximated with errors comparable to the s-states.","The extension to positronium shows the method transfers to any hydrogen-like reduced-mass system, giving mass-dependent screening shifts.","The Hellmann-Feynman result implies ∂E/∂c > 0, so energies rise monotonically with screening and s-states remain bound for all finite c—there is no critical screening in this potential, unlike the Yukawa case.","Combining the Coulomb reference (which underbinds) with the Kratzer reference (which overbinds) brackets the true energy, providing an error estimate without a separate high-precision calculation."],"fun_headline_variants":["Kratzer basis: 0.63% error for screened atoms","Screened Coulomb: 0.6% error with Kratzer basis","Positronium too: Kratzer-based screened atom energies","Screening mimics rotation: exact Kratzer trick","From hydrogen to positronium: analytic screened spectra"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The accuracy claims rest on the assumption that a Kratzer-shaped trial wave function with one adjustable scale β lies close to the true eigenfunction; for excited states the variational step is not a rigorous upper bound because the trial states are not orthogonalized against lower states.","fun_headline_variants_meta":{"raw":{"variants":["Kratzer basis: 0.63% error for screened atoms","Screened Coulomb: 0.6% error with Kratzer basis","Positronium too: Kratzer-based screened atom energies","Screening mimics rotation: exact Kratzer trick","From hydrogen to positronium: analytic screened spectra"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000883,"raw_usage":{"total_tokens":3628,"prompt_tokens":696,"completion_tokens":2932,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":440,"completion_tokens_details":{"reasoning_tokens":2843}},"tokens_in":440,"tokens_out":2932,"duration_ms":18102,"temperature":1.0,"reasoning_tokens":2843,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T13:09:03.104056+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A converged numerical solution of the radial Schrödinger equation for the 1s state at c=0.1 can settle the claim: if the Kratzer expectation value deviates from the exact energy by more than the stated 0.63%, or if the variational energy falls below the exact ground-state energy, the central accuracy claim fails.","supporting_citations":[],"review_version":1}