{"id":"dde45527-4f63-4630-b96f-4d1b81cded6a","arxiv_id":"1908.07374","paper_version":2,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"The integral over two spherical Bessel functions and a Gaussian, when dressed with a 3j symbol, equals a finite sum of modified spherical Bessel functions.","lead":"This paper derives an analytic finite-sum formula for the integral of two spherical Bessel functions times a Gaussian, weighted by an angular momentum coupling symbol. The result may speed up nuclear scattering calculations that use harmonic oscillator wavefunctions.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Result (2.9) rests entirely on the cited three-Bessel identity (2.1), which is neither re-derived nor numerically checked for generic λ3; an undetected slip there would invalidate the finite-sum formula.","rationale":"The reader correctly identified the uncited (in this paper) three-Bessel identity as the weakest assumption. I agree with that assessment. However, this is a verification gap rather than an identified error: the derivation is explicit, the closure step is standard, and the λ3=0 limit recovers the known Tabakin–Davies result, which provides independent support for the prefactor and the i_l structure. Cited prior work is a normal foundation for a methods paper, so I would not change the ACCEPT verdict on this basis alone. The proposed numerical test is cheap and would convert the residual risk into confidence. Thus verdict remains UNCHANGED.","tokens_in":70,"tokens_out":29238,"duration_ms":330755,"concrete_test":"Choose a generic allowed triple, e.g. (λ1,λ2,λ3)=(2,3,1) (triangle |3−2|=1≤1≤5, even sum 6), with k1=1.0, k2=1.3, α=0.7. Evaluate Eq. (2.9) by computing the 3j/6j symbols and i_l(k1k2/2α) with a standard library; evaluate the left-hand side ∫0∞ r^3 e^{-0.7r^2} j2(r)j3(1.3r) dr by adaptive quadrature (or a specialized oscillatory method) to 12 digits. Repeat for (λ1,λ2,λ3)=(3,3,2) and (1,2,1) with different k1,k2,α. If all agree to ~1e-10 the concern is settled; if any deviate, re-check Eq. (2.4) before Eq. (2.1) to localize the failure.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The derivation of the main formula (2.9) has three steps: (i) insert the closure relation (2.3) to rewrite the target integral as a k3-integral of a product of a one-Bessel Gaussian integral (2.5) and a three-spherical-Bessel integral; (ii) substitute the author's earlier formula (2.1) for the three-Bessel integral; (iii) evaluate the k3 integral via the Gegenbauer-type expansion (2.8). Step (ii) is the only place where external, nontrivial input enters: Eq. (2.1) is imported from refs. [8,10,12] and is not re-derived or numerically spot-checked in this paper. The sole internal check, Eq. (2.10), sets λ3=0, which forces λ1=λ2 and collapses the L and l sums to a single term; it therefore does not exercise the 3j/6j recoupling structure, the (k2/k1)^L factors, or the finite-sum combinatorics for generic triples. Since an algebraic error in Eq. (2.1) (a phase, a binomial factor, a 6j argument, or an unstated restriction) would propagate unchanged into the claimed finite-sum result, and since the paper gives no evidence that Eq. (2.1) has been re-verified, this is the weakest load-bearing point. The formal closure step, while distributionally delicate for oscillatory integrals, is standard and is supported by the λ3=0 limit; it is not the main residual risk.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper evaluates the integral of two spherical Bessel functions multiplied by a Gaussian and a power factor, with a 3j symbol prefactor, obtaining the finite-sum expression in Eq. (2.9). The derivation uses the closure relation for spherical Bessel functions, the known one-Bessel Gaussian integral (Eq. 2.5), a three-spherical-Bessel integral taken from the author's prior work (Eq. 2.1), and the Legendre expansion of the exponential (Eq. 2.8). The λ3=0 limit reproduces a known result of Tabakin and Davies (Eq. 2.10).","tokens_in":7,"tokens_out":31744,"duration_ms":551781,"significance":"If correct, the finite-sum formula is a useful analytic tool for nuclear scattering and momentum-space potential calculations, avoiding the infinite series of the standard Gradshteyn–Ryzhik expression. The derivation is compact and logically sound, and the λ3=0 check provides an independent consistency test. The paper ships no numerical verification for generic triples, which is the main residual risk, but the central mathematical structure appears sound.","major_comments":[],"minor_comments":[{"comment":"The central result Eq. (2.9) relies on the three-spherical-Bessel identity Eq. (2.1), which is cited from refs. [8,10,12] but not re-derived or numerically checked in this manuscript. The λ3=0 limit only exercises the single-term L=0 contribution and does not test the 3j/6j recoupling structure, the (k2/k1)^L factors, or the finite-sum combinatorics. A numeric spot-check for a generic triple (e.g., λ1=2, λ2=3, λ3=1) would substantially increase confidence in the final formula.","section":"§2, Eq. (2.1)"},{"comment":"The typesetting of the square roots is ambiguous: Eq. (2.5) appears to read sqrt(π)/2 but must be sqrt(π/2) to be consistent with the known λ=0 integral and with Eq. (2.10). Similarly, the prefactor in Eq. (2.9) should read sqrt(π/2), not sqrt(π)/2. Please ensure the radical sign covers the full fraction.","section":"§2, Eqs. (2.5), (2.9)"},{"comment":"There are several typos: 'MUL TIPLIED' in the title, 'Refrences' at the start of Section 2, and 'applicabale' before Eq. (2.9). These should be corrected.","section":"Title and text"},{"comment":"The Legendre expansion of the exponential in modified spherical Bessel functions is standard, but a specific reference (e.g., Abramowitz and Stegun) would help readers unfamiliar with this identity.","section":"§2, Eq. (2.8)"}],"recommendation":"minor_revision","confidential_remarks":"The manuscript leans heavily on the author's own earlier results (refs. [8,10,12]) for the three-Bessel integral. Since the current paper does not re-derive or independently check that formula, a numerical spot-check of Eq. (2.9) for a nontrivial triple would be especially valuable and would address any reader skepticism about self-citation. The mathematical derivation otherwise appears sound."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The main formula (2.9) is genuinely new: a finite sum over modified spherical Bessel functions for an integral that Gradshteyn and Ryzhik only gives as an infinite series. The derivation is transparent — closure relation, one-Bessel Gaussian integral, the author's earlier three-Bessel formula (2.1), and a Legendre expansion. The λ3=0 limit reproduces Tabakin–Davies, which is a solid independent check.\n\nThe load-bearing piece is (2.1), imported from earlier papers and not re-derived here. That is a real gap in self-containedness. The check in (2.10) only exercises λ3=0, which collapses the sums and does not stress the 3j/6j structure, the powers of k1 and k2, or the finite-sum combinatorics for generic triples. A numerical spot-check for one or two nonzero λ3 would have settled this. The closure relation (2.3) is written with an asymmetric 1/r² factor, but that is harmless under the integrals. I checked the algebra of the Δ to k3 change of variable and the Legendre expansion; the factors are correct.\n\nNothing here is a fatal flaw. The citation of the author's own prior work is appropriate, since (2.1) is published in J. Phys. A and Phys. Rev. C. The paper is exactly what it appears to be: a useful addition to the toolbox, not a new mechanism or a resolution of a big question. The target audience is people doing nuclear scattering or momentum-space potentials with Gaussian forms.\n\nIt deserves a serious referee. My recommendation would be to accept after the author adds a short numerical verification for a generic triple and cleans up the typos. The absence of that verification does not make the result suspect, but it would make the paper stronger.","headline":"A clean, narrowly scoped formula paper: the new finite-sum result is likely right, but the author should have numerically checked a generic case rather than leaning entirely on an earlier three-Bessel identity.","tokens_in":4877,"tokens_out":11021,"would_cite":true,"duration_ms":96855,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["33C10"],"pacs":[],"model":"deepseek-v4-flash","headline":"The Gaussian-weighted product of two spherical Bessel functions is exactly a finite sum of modified spherical Bessel functions under the angular-momentum triangle and parity conditions.","keywords":["spherical Bessel functions","modified spherical Bessel functions","Gaussian integral","3j symbols","6j symbols","closure relation","nuclear scattering","Gaussian potential"],"falsifier":"Take an admissible nontrivial case, for example $\\alpha=1$, $k_1=1$, $k_2=2$, and $\\lambda_1=\\lambda_2=\\lambda_3=2$, evaluate the right-hand side of (2.9), and compare it with high-precision numerical quadrature of the original oscillatory integral over $r$. Agreement to machine precision would support the identity; any disagreement beyond the quadrature error would refute it. Applying the same numerical check to the quoted three-Bessel formula (2.1) would then show whether the trouble lies in the starting identity or in the later reduction.","tokens_in":3838,"feed_emoji":"⚛️","tokens_out":16652,"duration_ms":146534,"temperature":0.7,"pith_summary":"The paper establishes a closed-form, finite-sum evaluation of the 3j-weighted radial integral of two spherical Bessel functions times a Gaussian, an object that appears in nuclear scattering when harmonic-oscillator wavefunctions are used or when a Gaussian potential is evaluated in momentum space. The central claim is that the angular-momentum constraints carried by the 3j symbol, $|\\lambda_2-\\lambda_1| \\le \\lambda_3 \\le \\lambda_1+\\lambda_2$ with $\\lambda_1+\\lambda_2+\\lambda_3$ even, turn what would otherwise be an infinite hypergeometric series into a finite sum over modified spherical Bessel functions of the first kind. A sympathetic reader should take away that this class of oscillatory integrals has an exact, compact analytic representation, so slow and delicate numerical integration of rapidly oscillating Bessel integrands can be replaced by direct summation.","feed_headline":"Bessel-Gaussian integral shrinks to a finite sum","feed_subtitle":"A closed finite-sum formula replaces slow oscillatory integration in nuclear scattering and Gaussian-potential matrix elements.","key_machinery":"The argument is carried by the spherical Bessel closure relation, $$\\int_0^\\infty $k^{2}$ j_L(kr)j_L(kr')\\,dk=\\frac{\\pi}{$2r^{2}$}\\delta(r-r'),$$ which inserts a third spherical Bessel function $j_{\\lambda_3}(k_3r)$ into the integrand and converts the target integral into an integral over $k_3$ of a known three-spherical-Bessel integral (equation 2.1) times a one-dimensional Gaussian transform (equation 2.5). The change of variable $k_3\\to \\Delta=(k_1^2+k_2^2-k_3^2)/(2k_1k_2)$ restricts $\\Delta$ to $[-1,1]$ and brings in the generating expansion $$\\exp\\!\\left(\\frac{k_1k_2\\$\\Delta$}{2\\$\\alpha$}\\right)=\\sum_L (2L+1)P_L(\\$\\Delta$)\\,i_L\\!\\left(\\frac{k_1k_2}{2\\$\\alpha$}\\right),$$ where $i_L(x)=\\sqrt{\\pi/(2x)}\\,I_{L+1/2}(x)$ is the modified spherical Bessel function of the first kind. This expansion is what turns the leftover integral into the finite sum of modified spherical Bessel functions.","core_discovery":"The paper's central discovery is the closed identity (2.9): for nonnegative integers $\\lambda_1,\\lambda_2,\\lambda_3$ with $|\\lambda_2-\\lambda_1|\\le \\lambda_3\\le \\lambda_1+\\lambda_2$ and $\\lambda_1+\\lambda_2+\\lambda_3$ even, and for positive $k_1,k_2,\\alpha$, the 3j-weighted integral $$\\begin{pmatrix}\\lambda_1 & \\lambda_2 & \\lambda_3\\\\ 0 & 0 & 0\\end{pmatrix}\\int_0^\\infty $r^{{\\lambda_3+2}}$$e^{{-\\alpha r^2}}$ j_{\\lambda_1}(k_1r)j_{\\lambda_2}(k_2r)\\,dr$$ equals a finite sum over $L$ from $0$ to $\\lambda_3$ and over the finite range of $l$ allowed by the angular-momentum coefficients, with terms made of 3j symbols, 6j symbols, and modified spherical Bessel functions $i_l(k_1k_2/(2\\alpha))$, multiplied by the Gaussian factor $\\exp[-(k_1^2+k_2^2)/(4\\alpha)]$. The finiteness is the point: the 3j constraints truncate sums that would otherwise run forever as an infinite hypergeometric series. Setting $\\lambda_3=0$, which forces $\\lambda_1=\\lambda_2=\\lambda$, reduces the identity to the single-term formula with $i_\\lambda(k_1k_2/(2\\alpha))$, recovering the previously published special case.","pith_inferences":["The same closure-relation insertion should produce finite closed forms for Gaussian-weighted products of three or more spherical Bessel functions whenever angular-momentum constraints truncate the intermediate sums; the paper does not carry out such cases.","Because the 3j symbol vanishes outside the stated triangle and parity window, identity (2.9) does not directly give the unweighted integral; using the orthogonality of 3j symbols to unwind the weight would yield the bare integral for arbitrary integer $\\lambda_3$.","The right-hand side of (2.9) is analytic in $\\alpha$ on the half-plane $\\mathrm{Re}\\,\\alpha>0$, so the identity plausibly extends by analytic continuation from real positive $\\alpha$ to complex Gaussian parameters, covering damped or oscillating momentum-space potentials.","The finite-sum form is well suited to reuse in momentum-space codes that currently evaluate such integrals by adaptive quadrature, since it replaces the expensive oscillatory integration with a small number of spherical Bessel evaluations."],"forward_implications":["The 3j-weighted two-spherical-Bessel Gaussian integral can be evaluated exactly by summing finitely many 3j and 6j coefficients and modified spherical Bessel functions; no infinite series or oscillatory radial quadrature is needed.","The special case $\\lambda_3=0$ (hence $\\lambda_1=\\lambda_2$) reduces the general result to a known one-term formula, providing an internal consistency check.","Nuclear scattering calculations using harmonic-oscillator target wavefunctions and momentum-space Gaussian potentials can substitute this closed form for numerical integration of the radial integral.","The derivation shows a general route: inserting the spherical Bessel closure relation reduces a product of two Bessel functions to a known three-Bessel integral whose angular-momentum structure is already understood."],"supporting_citations":[{"why":"Supplies the three-spherical-Bessel integral (2.1), the starting identity into which the Gaussian Bessel factor is inserted via the closure relation.","marker":"[8-10, 12]"},{"why":"Provides the Gaussian transform (2.5) used to integrate out the inserted third Bessel function, and the infinite hypergeometric series (1.1) that motivates the search for a finite form.","marker":"[21]"},{"why":"Supplies the previously published lambda_3=0 special case against which the general result is checked, yielding agreement (2.10).","marker":"[19]"}],"fun_headline_variants":["Finite sum tames Bessel-Gaussian integral","Bessel-Gaussian integral now a closed sum","Closed form for two-Bessel Gaussian integral","Exact finite sum for Bessel-Gaussian overlap","Finite sum replaces infinite series for Bessel-Gaussian"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The derivation rests on the three-spherical-Bessel integral (2.1), quoted from the author's earlier papers without being re-derived or independently checked here; if that quoted identity has a hidden restriction or an error, the final closed form (2.9) would not stand.","fun_headline_variants_meta":{"raw":{"variants":["Finite sum tames Bessel-Gaussian integral","Bessel-Gaussian integral now a closed sum","Closed form for two-Bessel Gaussian integral","Exact finite sum for Bessel-Gaussian overlap","Finite sum replaces infinite series for Bessel-Gaussian"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00058,"raw_usage":{"total_tokens":2760,"prompt_tokens":1003,"completion_tokens":1757,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":619,"completion_tokens_details":{"reasoning_tokens":1683}},"tokens_in":619,"tokens_out":1757,"duration_ms":13603,"temperature":1.0,"reasoning_tokens":1683,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:04:12.207882+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take an admissible nontrivial case, for example $\\alpha=1$, $k_1=1$, $k_2=2$, and $\\lambda_1=\\lambda_2=\\lambda_3=2$, evaluate the right-hand side of (2.9), and compare it with high-precision numerical quadrature of the original oscillatory integral over $r$. Agreement to machine precision would support the identity; any disagreement beyond the quadrature error would refute it. Applying the same numerical check to the quoted three-Bessel formula (2.1) would then show whether the trouble lies in the starting identity or in the later reduction.","supporting_citations":[{"cited_title":"Gradshteyn and I.M","cited_arxiv_id":null,"evidence_quote":"Provides the Gaussian transform (2.5) used to integrate out the inserted third Bessel function, and the infinite hypergeometric series (1.1) that motivates the search for a finite form."},{"cited_title":"Tabakin and K.T.R","cited_arxiv_id":null,"evidence_quote":"Supplies the previously published lambda_3=0 special case against which the general result is checked, yielding agreement (2.10)."}],"review_version":1}