{"id":"8b65161d-ae50-4b4b-b1ac-902e8301cc87","arxiv_id":"1908.02324","paper_version":2,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A comprehensive set of D=3 and D=3-2epsilon Coulomb expectation values, with regularized divergent cases and general Laguerre integral formulas, is derived and tabulated.","lead":"This paper derives and tabulates expectation values of Coulomb operators in three dimensions and in 3-2epsilon dimensions, including dimensionally regularized expressions for divergent ones. It also gives general Laguerre polynomial integral formulas meant to simplify QED bound-state energy calculations.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (III.17) omits a factor: the coefficient of <(V')^2> should contain 3(1-2ε)^2, not 3(1-2ε); as printed it invalidates the exact-recursion check at the ε^0 order of the central S-state results.","rationale":"The reader's identified weakest assumption, the short-distance truncation of Ln0 around Eqs. (III.4)-(III.5), is not where the main risk lies: the Ltail terms become genuinely finite after the Lhat subtraction, so evaluating them at ε=0 is reasonable. The load-bearing problem is instead that Eq. (III.17), presented as an exact D-dimensional recursion and used to cross-check the tables, has an incorrect coefficient at O(ε). Deriving it from the paper's own Eq. (III.13) yields a factor (1-2ε)^2 where the paper prints 3(1-2ε). Because divergent S-state expectation values carry 1/ε poles, this O(ε) error changes O(1) finite parts, precisely the quantities that the paper tabulates as regularized values. A spot check of (A.13c) and (A.13e) against the printed relation for ℓ=0 fails at O(1), while the corrected coefficient restores consistency, indicating the values may be right but the supporting identity is wrong as printed. This warrants a conditional verdict: the paper should be accepted only after the authors correct Eq. (III.17), re-verify all identities and Appendix A entries that depend on it, and confirm that no tabulated expectation value shifts. This is not an ad hominem or an attack on the central construction; it is a specific, checkable algebraic defect in the paper's own derivation.","tokens_in":43018,"tokens_out":25467,"duration_ms":269215,"concrete_test":"Independently re-derive Eq. (III.17) from Eq. (III.13) at s = -2+4ε and compare the coefficient of <(Vbar')^2>. If the correct coefficient is [3(1-2ε)^2 - 4ℓ(ℓ+1-2ε)]/(1-2ε), then substitute the S-state entries (A.13c) and (A.13e), including the O(ε) expansion of φbar^2 and μ^{2ε}, into the corrected ℓ=0 identity at O(ε^0) and verify that the finite parts cancel. Then search all identities derived from III.17, including (A.15), for propagated missing-square errors.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The paper's central regularized S-state values for <Vbar^3> and <(Vbar')^2> are supported in part by the claim that Eq. (III.17) is exact in D dimensions. That equation is not correct as printed. Starting from Eq. (III.13) with s = -2+4ε and using s+1 = -(1-4ε), 2s+1+2ε = -(3-10ε), <r^s> = <Vbar^2>/β^2, <r^{s-1+2ε}> = -<Vbar^3>/β^3, and <r^{s-2}> = <(Vbar')^2>/[(1-2ε)^2 β^2], then multiplying by -β^2/2 gives the coefficient of <(Vbar')^2> as -s[s^2-(1-2ε)^2-4ℓ(ℓ+1-2ε)]/[2(1-2ε)^2] = [3(1-2ε)^2 - 4ℓ(ℓ+1-2ε)]/(1-2ε). The printed Eq. (III.17) has 3(1-2ε) in the numerator, missing the square. For ℓ=0 the correct coefficient is 3(1-2ε) = 3 - 6ε, not 3. This is not a cosmetic difference: <(Vbar')^2> has a -2/ε pole, so the missing -6ε term shifts the O(ε^0) finite part by +12C, where C is the common prefactor of the S-state entries. A direct O(ε^0) substitution of (A.13c) and (A.13e) into the printed identity fails by exactly this amount, while the corrected identity cancels. Thus the tabulated values appear consistent with the corrected relation, but any derivation or consistency check using (III.17) as written is invalid at the order at which the S-state results are quoted. This is the main load-bearing weakness: not the short-distance truncation of Ln0, which is defensible, but an incorrect exact identity used as support for the central regularized expectation values.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper develops the D=3−2ε dimensional Coulomb problem as a tool for regularizing divergent expectation values in atomic bound-state calculations. It constructs the radial wave function as a double power series in ρ and ε, derives perturbative energy and wave-function corrections, and tabulates a large set of coordinate-space and momentum-space expectation values in Appendix A. The central technical claims are that the dimensionally regularized divergent S-state expectation values such as ⟨V̄³⟩ and ⟨(V̄′)²⟩ are correctly obtained, that several finite D-dimensional expectation values such as ⟨r^{-2+4ε}∂_r²⟩ have ε→0 limits different from their direct three-dimensional limits, and that the tabulated entries are consistent with recursion relations, the Feynman-Hellmann theorem, and numerical integration. The paper also presents general integral formulas for ordinary and subtracted associated Laguerre polynomials, including logarithmic and double-logarithmic cases, and applies the results to an NRQED energy-level correction.","tokens_in":43579,"tokens_out":8782,"duration_ms":93831,"significance":"If the tabulated results are correct, this is a very useful reference for dimensional-regularization calculations of Coulombic bound-state properties, particularly for mα⁶ NRQED corrections. The paper is systematic and unusually cross-checked: it verifies many entries against recursion relations, Feynman-Hellmann identities, and direct numerical integration in three dimensions, and it clearly warns that S-state divergent entries are quoted only through O(ε⁰). The new formulas for subtracted Laguerre integrals and for diharmonic sums are independently valuable. The main weakness identified below is a concrete coefficient error in an exact D-dimensional identity that is cited in support of the consistency claims; this is localized and correctable, but it must be fixed and re-verified before the paper can serve as a trusted reference.","major_comments":[{"comment":"The coefficient of ⟨(V̄′)²⟩ in Eq. (III.17) is misprinted. Starting from Eq. (III.13) with s=−2+4ε and using ⟨r^s⟩=⟨V̄²⟩/β², ⟨r^{s−1+2ε}⟩=−⟨V̄³⟩/β³, and ⟨r^{s−2}⟩=⟨(V̄′)²⟩/[(1−2ε)²β²], multiplication by −β²/2 gives the coefficient [3(1−2ε)²−4ℓ(ℓ+1−2ε)]/(1−2ε), not [3(1−2ε)−4ℓ(ℓ+1−2ε)]/(1−2ε). The printed form is therefore not exact in D dimensions as claimed immediately below the equation. For ℓ=0 the correct coefficient is 3−6ε rather than 3; since ⟨(V̄′)²⟩_{n0} has a −2/ε pole (A.13e), this changes the O(ε⁰) finite part of the identity by +12πφ̄²_n m_r(Zα)³μ̄^{2ε}. Because the identity is called exact and is cited again in the list of exact relations before Eq. (A.15), the error invalidates the stated recursion-relation consistency check at precisely the order at which the S-state results are quoted. The authors should correct the coefficient and confirm that all table entries satisfy the corrected relation, or amend the consistency claim.","section":"§III, Eq. (III.17)"}],"minor_comments":[{"comment":"The warning that S-state expectation values involving V̄ are shown only through O(ε⁰) is essential and should also be repeated at the first use of these entries, e.g., in Eq. (IV.5), so that readers do not mistake the displayed finite part for a complete all-orders D-dimensional result.","section":"Appendix A, after Eq. (A.13)"},{"comment":"The symbols H_n, H_n^{(2)}, and diH±(n,m) are defined only in Appendix D; a forward reference at the first occurrence would improve readability, since Eq. (II.30) is otherwise hard to parse.","section":"§II, Eq. (II.30)"},{"comment":"The caption uses diH⁺(7,5), diH⁺(10,−3), and diH⁻(7,12) without defining the symbols; either define diH± in the caption or explicitly refer the reader to Eq. (D.29).","section":"Fig. 2 caption"}],"recommendation":"major_revision","confidential_remarks":"The coefficient error in Eq. (III.17) is the kind of defect that does not negate the paper's central derivation but undermines a prominent cross-check claim. The correctable nature of the issue makes rejection inappropriate; however, because the paper's value lies in the reliability of its tables and identities, the authors should fix the typo and explicitly re-verify the affected consistency checks before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Thanks for passing this along. I read the paper and also the stress-test note. The note is right: Eq. (III.17) is wrong as printed. Starting from (III.13) with s = -2+4ε, the coefficient of ⟨(V̄')²⟩ comes out [3(1-2ε)² - 4ℓ(ℓ+1-2ε)]/(1-2ε), not the printed [3(1-2ε) - ...]/(1-2ε). For ℓ=0 that's 3-6ε, not 3. Since ⟨(V̄')²⟩ has a 1/ε pole, the missing -6ε shifts the finite O(ε⁰) part by 12 times the overall prefactor. The paper's claim that (III.17) is exact is central to the S-state regularization story, and any check using it as printed will fail. The tabulated values appear consistent with the corrected relation, so the numerical output may be fine, but the derivation as written doesn't support it.\n\nThat said, this is a genuinely useful reference. The D=3-2ε results are mostly new, the general formulas for logarithmic and subtracted Laguerre integrals are a real contribution, and the diharmonic sums make the expansions systematic. The authors are honest about limitations: leading-order-in-ε entries are labeled, and the short-distance truncation assumption is discussed. The cross-checks against recursion relations, Feynman-Hellmann, and numerical integration give confidence that the bulk of the tables are correct. I would not treat the typo as evidence of broader sloppiness; it's a single but load-bearing mistake.\n\nThe soft spots beyond (III.17): the Laguerre tables are wide but not machine-checked, so transcription errors are possible; and the S-state divergent values rely on the assumption that only the first few terms of the generalized Laguerre series matter at short distances. That's defensible physically, but it would be nice to see a more formal justification or an independent check.\n\nWho is this for? Anyone using dimensional regularization in NRQED or bound-state QED at mα⁶ and beyond. It's a reference work, not a breakthrough. It deserves a serious referee, but the referee needs to catch (III.17) and make the authors fix it before publication. If I were the editor, I'd send it out with a clear request to correct the identity and re-verify the affected S-state entries. I'd still cite it for the tables, but I'd stay away from the un-repaired recursion relation.","headline":"A valuable reference for D-dimensional Coulomb expectation values, but Eq. (III.17) has a real coefficient error that undercuts the printed consistency checks.","tokens_in":44019,"tokens_out":4751,"would_cite":true,"duration_ms":42816,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81Q05","81V45","33C45","81T15"],"pacs":["03.65.-w","31.15.-p","12.20.-m"],"model":"deepseek-v4-flash","headline":"The paper establishes that the dimensionally regularized Coulomb problem in D=3−2ε dimensions yields finite regularized values for divergent expectation values such as ⟨V̄³⟩ and ⟨(V̄′)²⟩, and that the tabulated finite values include ε→0…","keywords":["Coulomb expectation values","dimensional regularization","D=3−2ε dimensions","associated Laguerre polynomials","subtracted Laguerre integrals","momentum space brackets","NRQED bound states","diharmonic numbers"],"falsifier":"For a chosen state, say n=2, ℓ=0, evaluate ⟨(V̄′)²⟩ by numerically solving the D=3−2ε radial Schrödinger equation with a small-distance cutoff, extract the 1/ε and constant terms, and compare them with the Laurent expansion in (A.13e); a mismatch would show that the short-distance isolation of L_{n0} misses contributions. A second check: compute ⟨$r^{{−2+4ε}}$∂_r²⟩ for ℓ=0 using a different regulator such as a momentum cutoff and test whether the δℓ=0 term in (III.8) appears.","tokens_in":42841,"feed_emoji":"⚛️","tokens_out":6166,"duration_ms":61645,"temperature":0.7,"pith_summary":"The paper establishes that the Coulomb bound-state problem in D=3−2ε dimensions works as a dimensional-regularization scheme for expectation values that diverge in three dimensions. Divergent S-state quantities such as ⟨V̄³⟩ and ⟨(V̄′)²⟩ are shown to have well-defined Laurent expansions in ε, with the 1/ε pole encoding the short-distance divergence and the finite part serving as the regularized value. The paper also identifies a subtlety: for some finite operators such as ⟨$r^{{−2+4ε}}$∂_r²⟩, the ε→0 limit of the D-dimensional expectation value differs from the value computed directly in three dimensions. A systematic table of finite and regularized expectation values and momentum-space brackets is provided, along with general integration formulas for associated Laguerre polynomials and their subtracted versions. These results give practitioners a ready-made toolkit for dimensional regularization in QED corrections to Coulombic bound-state energies.","feed_headline":"Coulomb divergences tamed in D=3−2ε dimensions","feed_subtitle":"A full table of regularized expectation values puts QED bound-state corrections on solid ground.","key_machinery":"The key object is the D-dimensional radial wave function, written as φ̄_{nℓ} Ω_{D−1}^{1/2} $e^{{−ρ/2}}$ ρ^ℓ L_{nℓ}(ρ) with ρ=2γ̄_{nℓ}r. Here L_{nℓ}(ρ) is not a standard Laguerre polynomial but a generalized power series L_{nℓ}(ρ)=Σ_{j,k} a_{jk} n̄_{nℓ}^k $ρ^{{j+2εk}}$, with coefficients fixed by a two-index recursion relation; in the ε→0 limit it reduces to the usual associated Laguerre polynomial. For divergent expectation values, the paper isolates the short-distance part of L_{nℓ} (the first few terms) as \\widehat L_{nℓ}, treats its integrals analytically in ε, and evaluates the remainder at ε=0 using subtracted associated Laguerre polynomials—polynomials with low powers of the variable removed. General formulas for integrals of one or two (subtracted) associated Laguerre polynomials times powers of x and ln x, together with D-dimensional Fourier transforms of momentum-space brackets, carry the tabulation.","core_discovery":"The central claim is that every divergent Coulomb expectation value relevant to bound-state QED corrections can be assigned a finite regularized value by working in D=3−2ε dimensions. The paper shows, for example, that ⟨(V̄′)²⟩_{n0} = π m_r (Zα)³ φ̄_n² μ̄^{2ε} {−2/ε −8 ln(μ n/(2 m_r Zα)) + 8H_n + 4/(3n²) −4/n −16/3 + O(ε)}, and that ℓ>0 cases remain finite. Similarly, ⟨V̄³⟩ has a 1/ε pole for S states and a finite closed form for ℓ>0. The paper further demonstrates that taking D→3 does not always commute with taking the expectation value: entries like ⟨$r^{{−2+4ε}}$∂_r²⟩ and ⟨$r^{{4ε}}$p⁴⟩ differ from their three-dimensional counterparts in explicit δℓ=0 terms. It also gives exact D-dimensional relations, such as 2m_r⟨(V̄′)²⟩ = ⟨p²V̄p²⟩ − ⟨p⁴V̄⟩, and a large table of coordinate-space expectation values and momentum-space brackets in both three and D dimensions.","pith_inferences":["A testable extension is to apply the same subtracted-Laguerre technique to other D=3−2ε central potentials; the harmonic oscillator would show whether the δℓ=0 non-commutativity of limits is generic or Coulomb-specific.","A consequence the paper leaves implicit is that any NRQED mα⁶ calculation using these tables inherits the ε→0-limit subtleties; a different regulator such as point-splitting should reproduce the same finite parts, which would be a direct check.","One could extend the tables to relativistic corrections by inserting the D-dimensional wave function into the Dirac–Coulomb problem and taking the nonrelativistic expansion; the same short-distance split would regulate the newly divergent operators."],"forward_implications":["Any QED correction to Coulomb bound states that uses dimensional regularization can quote the tabulated expectation values directly, including the divergent S-state ⟨V̄³⟩ and ⟨(V̄′)²⟩.","Operators whose radial derivatives generate negative powers of ρ acquire δℓ=0 correction terms in the ε→0 limit, so dimensional regularization and the D→3 limit must be taken in the right order.","The exact D-dimensional identities from recursion relations, the Feynman–Hellmann theorem, and momentum-space brackets let a user check or rederive any tabulated entry without a separate integration.","A sample NRQED calculation shows the CX1 contribution to bound-state energies reduces to −4m_r(...)⟨(V̄′)²⟩, so the tables plug directly into order-mα⁶ energy calculations.","The subtracted-Laguerre integral formulas provide a general technique for evaluating short-distance divergent expectation values for any potential whose wave function has a generalized Laguerre expansion."],"supporting_citations":[{"why":"Supplies the momentum-space bracket method and the n=1 D-dimensional results that this paper generalizes to arbitrary n.","marker":"[1]"},{"why":"Introduced the generalized power-series solution for the D=3−2ε wave function that underlies the whole regularization scheme.","marker":"[2]"},{"why":"Sets out the NRQED Bethe–Salpeter formalism used to turn bound-state energy shifts into the tabulated expectation values.","marker":"[11]"},{"why":"Provides the D-dimensional Fourier-transform formula used to convert momentum-space brackets into coordinate-space expectation values.","marker":"[15]"},{"why":"Gives an exact D-dimensional Coulomb solution that serves as the lowest-order problem for the perturbative ε expansion.","marker":"[22]"},{"why":"Another exact D-dimensional Coulomb solution used as the perturbation basis for energies and wave functions near ε=0.","marker":"[23]"},{"why":"Earlier result for the ⟨ln q⟩ momentum-space bracket that the paper's general formula reproduces as a check.","marker":"[34]"},{"why":"Defines NRQED, the low-energy effective field theory used in the paper's sample energy-level calculation.","marker":"[37]"}],"fun_headline_variants":["Dimensional regularization tames Coulomb divergences","Coulomb expectation values regularized in D=3−2ε","D=3−2ε method for divergent Coulomb integrals","Noncommuting limits in Coulomb expectation values","General formulas for D-dimensional Coulomb expectation values"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that a divergent S-state expectation value can be computed by isolating the first few terms of the generalized Laguerre series as the only short-distance-sensitive part, then evaluating the remaining infinite series at ε=0.","fun_headline_variants_meta":{"raw":{"variants":["Dimensional regularization tames Coulomb divergences","Coulomb expectation values regularized in D=3−2ε","D=3−2ε method for divergent Coulomb integrals","Noncommuting limits in Coulomb expectation values","General formulas for D-dimensional Coulomb expectation values"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000992,"raw_usage":{"total_tokens":4316,"prompt_tokens":1168,"completion_tokens":3148,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":784,"completion_tokens_details":{"reasoning_tokens":3074}},"tokens_in":784,"tokens_out":3148,"duration_ms":23860,"temperature":1.0,"reasoning_tokens":3074,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:46:57.673762+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"For a chosen state, say n=2, ℓ=0, evaluate ⟨(V̄′)²⟩ by numerically solving the D=3−2ε radial Schrödinger equation with a small-distance cutoff, extract the 1/ε and constant terms, and compare them with the Laurent expansion in (A.13e); a mismatch would show that the short-distance isolation of L_{n0} misses contributions. A second check: compute ⟨$r^{{−2+4ε}}$∂_r²⟩ for ℓ=0 using a different regulator such as a momentum cutoff and test whether the δℓ=0 term in (III.8) appears.","supporting_citations":[{"cited_title":"Czarnecki, K","cited_arxiv_id":null,"evidence_quote":"Supplies the momentum-space bracket method and the n=1 D-dimensional results that this paper generalizes to arbitrary n."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduced the generalized power-series solution for the D=3−2ε wave function that underlies the whole regularization scheme."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Sets out the NRQED Bethe–Salpeter formalism used to turn bound-state energy shifts into the tabulated expectation values."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the D-dimensional Fourier-transform formula used to convert momentum-space brackets into coordinate-space expectation values."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives an exact D-dimensional Coulomb solution that serves as the lowest-order problem for the perturbative ε expansion."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Another exact D-dimensional Coulomb solution used as the perturbation basis for energies and wave functions near ε=0."},{"cited_title":"Titard and F","cited_arxiv_id":null,"evidence_quote":"Earlier result for the ⟨ln q⟩ momentum-space bracket that the paper's general formula reproduces as a check."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines NRQED, the low-energy effective field theory used in the paper's sample energy-level calculation."}],"review_version":1}