{"id":"4a22526d-c550-4042-b336-5e1497a5132d","arxiv_id":"2506.08040","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"An exact one-dimensional integral formula gives the electrostatic potential of a disk whose charge density is (1-u^2)^nu u^p, and numerical plots reveal a crossover in the in-plane potential shape around p=-0.34.","lead":"This paper derives an exact formula for the electric potential of a flat disk whose charge density varies with radius according to a two-parameter power-law family. It recovers known uniform and equipotential disk solutions as special cases and shows how the in-plane potential changes from flat to edge-peaked as the profile is varied.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (17)'s advertised range ν>-1, p>-1 is only as secure as the unproved interchange of the Bessel and u'-integrals at z=0; that step is the load-bearing gap.","rationale":"After working through the derivation, I find no obvious algebraic error in Eqs. (10)-(17). For z>0 the double integral in Eq. (10) is absolutely convergent, so applying identity (11) and Fubini is legitimate; the resulting integral (17) is well-defined for ν>-1, p>-1 and, by dominated convergence, is continuous down to z=0. This is why I do not reject the paper's main formula. However, the manuscript does not supply these arguments. The claim that Eq. (17) is exact for the entire range ν>-1, p>-1 therefore rests on an unproved interchange step at z=0. Since the paper's own qualitative discussion of Figure 2 contradicts the reported threshold p_c≈-0.34, the numerical side of the paper also needs substantiation. The reader's weakest assumption identifies the same convergence/interchange issue, and the CONDITIONAL verdict remains appropriate: the paper should either prove the interchange and limit, or restrict the claim and document the numerical threshold.","tokens_in":9696,"tokens_out":20709,"duration_ms":181999,"concrete_test":"For the most singular admissible pair, ν=-0.9, p=-0.9, evaluate Eq. (17) at z=0, u=0.5 and compare it with a direct high-precision quadrature of the original Coulomb integral, Eq. (1) (or Eq. (4) with a large k-cutoff and Richardson extrapolation). Agreement to 1e-6 relative would demonstrate that the interchange does not fail in the hardest parameter corner; disagreement would invalidate the claimed range. Additionally, recompute p_c from Eq. (17) by locating the sign change of dV_p/du at u=0 and verify whether p_c≈-0.34 is reproduced.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Equation (17) is reached by substituting identity (11) into the double integral (10) and reducing it to a single quadrature. For z>0 the double integral is absolutely convergent, so the substitution is formally justified. The unsupported step is the plane limit z→0: the inner Bessel integral is then only conditionally convergent (Weber–Schafheitlin type), and the paper gives no argument that the interchange of integration survives, nor that the limit z→0+ commutes with the u'-integral for the full range ν>-1, p>-1. The endpoint singularities at u'=0 and u'=1 are integrable but uncontrolled for p,ν near -1. If the interchange or the limit fails on any sub-range, the exactness of Eq. (17) and the derived threshold p_c≈-0.34 are not established there. A separate internal inconsistency reinforces this: the text states that for -1<p<0 the in-plane potential decreases monotonically, then reports p_c≈-0.34, which implies a non-monotonic regime for p in (-0.34,0); the numerical threshold therefore needs independent documentation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the electrostatic potential of a thin charged disk with surface density sigma(u)=sigma0 (1-u^2)^nu u^p. Starting from the Bessel-function expansion of the Coulomb kernel, the authors reduce the potential to a single integral over a complete elliptic integral of the first kind, Eq. (17), claimed for all real nu>-1 and p>-1. The uniform disk (p=nu=0) and the equipotential edge-concentrated disk (p=0, nu=-1/2) are recovered as special cases. The paper then analyzes the in-plane potential for nu=-1/2 as a function of p and reports a critical exponent p_c approximately -0.34 separating monotonic from non-monotonic radial behavior.","tokens_in":9854,"tokens_out":25567,"duration_ms":263028,"significance":"The final integral representation is elegant and computationally convenient, since only standard special functions are involved. The derivation is not fitted to benchmarks; the two limiting cases are genuine checks, and the algebra from Eq. (4) to Eq. (17) is internally consistent once the typo in Eq. (11) is corrected. If the z=0 limiting step is supplied and the p_c analysis is documented, the result would be a useful compact addition to the classical electrostatics literature. The claimed crossover in the in-plane potential is the main new physical observation, but it currently rests on an undocumented numerical estimate.","major_comments":[{"comment":"The identity (11) is misprinted. Gradshteyn and Ryzhik 6.612.3 requires the prefactor (u u')^(-1/2), not (u u')^(1/2). As printed, substituting (11) into (10) produces an extra factor u' and sqrt(u), so it does not lead to Eq. (14). Since Eqs. (14), (15), and (17) are mutually consistent and reproduce the uniform-disk and equipotential benchmarks, this is a transcribal error in a load-bearing identity; it must be corrected, with the validity conditions (u,u'>0 and Re z>0) stated explicitly.","section":"Eq. (11) and Eq. (14)"},{"comment":"The reduction of the double integral (10) to the single integral (17) is fully justified only for z>0, where the integrand is absolutely integrable. For z=0, the inner Bessel integral is only conditionally convergent and is logarithmically divergent at u'=u, so the interchange of the q- and u'-integrations and the passage to the plane cannot be taken for granted. Eq. (17) is used at z=0 throughout Section 3, and the determination of p_c depends on it. Please supply an explicit limiting argument (for example, dominated convergence for z to 0+ with the bound K(m) <= C(1+|ln(1-m)|) and 1-m >= c((u-u')^2+z^2)), and treat u=0 separately.","section":"Eqs. (10)-(17), Section 3"},{"comment":"The text states that for -1<p<0 the in-plane potential decreases monotonically from the center, but then reports a critical value p_c approximately -0.34 above which a maximum appears. For p in (-0.34,0) both statements cannot be true. The contradiction needs to be resolved, and the numerical estimation of p_c needs documentation (method, grid, tolerance), since this crossover is the paper's main new qualitative result.","section":"Section 3, p_c approximately -0.34"}],"minor_comments":[{"comment":"The factorial notation nu!, (p/2)!, and related expressions is used for real parameters; please either replace these with Gamma functions or explicitly state the convention Gamma(x+1)=x!.","section":"Eqs. (7)-(9) and (18)"},{"comment":"There are numerous typos, including 'propertied', 'mathemarical', 'baxkground', 'Rigdly', 'Boudaries', and 'different forms'. Please proofread the manuscript carefully.","section":"Throughout"},{"comment":"Reference [39] (Sonine) is listed but not cited in the text.","section":"References"},{"comment":"The prefactor in Eq. (17) is typeset as 'eVo', while Eq. (18) defines tilde V_o; please unify the notation.","section":"Eq. (17) and Eq. (18)"}],"recommendation":"major_revision","confidential_remarks":"The central formula appears correct and the benchmarks are reassuring, but the derivation has a typo in a key identity, the z=0 passage is not justified, and the p_c estimate is undocumented. These are fixable within the scope of the manuscript, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear [Colleague],\n\nThis paper derives an exact one-dimensional integral for the off-axis potential of a charged disk whose surface density is σ(u)=σ0 (1-u²)^ν u^p. The result, Eq. (17), is a single quadrature over a complete elliptic integral K; it reduces to the known uniform-disk and equipotential-disk cases at (ν=0,p=0) and (ν=-1/2,p=0). That is genuinely useful: instead of a different integration for each profile, you get one formula valid for all ν,p>-1. The derivation from the Bessel expansion is standard, and the benchmarks give me confidence the algebra is right.\n\nThe soft spots are real but manageable. First, the title promises 'arbitrary radial charge profiles' but the exact reduction covers only the two-parameter power-law family. Eq. (4) is general, but that is just the standard Bessel representation; the new content is the closed form for this family. Second, the paper does not discuss the z→0 limit of the Bessel identity. For z>0 the double integral is absolutely convergent and the interchange is fine; at z=0 one should appeal to continuity and dominated convergence, which is likely true for ν,p>-1 but is not stated. This is a rigor gap, not a fatal one. Third, the internal description of the in-plane potential is inconsistent: the text says Vp(u) decreases monotonically for -1<p<0, then reports a critical p_c≈-0.34 with a maximum appearing for p>p_c. Those two statements cannot both be true for p in (-0.34,0). The threshold p_c is also just 'estimated' with no method or error bar. A referee should ask for a proper definition and numerical procedure for p_c. Finally, the normalization uses factorial notation for real p/2 and ν; it should be gamma functions, a minor notation fix.\n\nOverall, the mathematics is sound in the main line, the formula is a real convenience, and the benchmarks are honest. It is not a landmark, but it is a competent contribution to a classical problem. I'd send it to a serious referee, with the request that they verify the p_c calculation and clean up the monotonicity statement.","headline":"A solid, useful but modest extension of known disk-potential results; the main formula checks out, but the reported crossover needs a proper numerical description and one description of the in-plane potential is internally inconsistent.","tokens_in":10443,"tokens_out":8206,"would_cite":true,"duration_ms":76031,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["41.20.Cv"],"model":"deepseek-v4-flash","headline":"The paper derives an exact single-integral representation, Eq.","keywords":["electrostatic potential","charged disk","elliptic integrals","Bessel functions","radial charge distribution","axial symmetry","equipotential disk","Legendre functions"],"falsifier":"Directly evaluate the double integral in Eq. (10) by numerical quadrature over q and u' for a parameter set near the reported crossover, say ν = −1/2, p = −0.34, at z = 0 and u = 0.5, and compare the result with the single-integral formula Eq. (17) to machine precision; any difference beyond roundoff that persists under adaptive refinement would falsify the exactness claim for that sub-range. A second check is to test endpoint convergence with p → −1⁺ or ν → −1⁺, where the integrand becomes singular at u' = 0 or u' = 1.","tokens_in":9445,"feed_emoji":"⚡","tokens_out":8098,"duration_ms":78213,"temperature":0.7,"pith_summary":"The paper aims to establish an exact, single-integral formula for the electrostatic potential, anywhere in space, of a thin charged disk whose radial charge profile is σ(u) = σ0(1−$u^{2}$)^ν u^p, for any real ν > −1 and p > −1. It matters because this family covers the uniformly charged disk, edge-concentrated distributions, and center-depleted profiles in one unified expression, and because the remaining integral involves only the complete elliptic integral of the first kind, which is available in standard scientific libraries. The derivation collapses the usual double Bessel integral to one quadrature using a classical product-of-Bessel-functions identity and a relation between a Legendre function and an elliptic integral. The paper also reports a qualitative transition in the in-plane potential at an estimated critical exponent p_c ≈ −0.34, where the radial maximum appears or disappears.","feed_headline":"One integral gives the potential of a radially charged disk","feed_subtitle":"Exact formula unifies uniform, edge-concentrated, and intermediate profiles; only the elliptic integral K remains.","key_machinery":"The load-bearing object is the product-of-Bessel-functions identity (Eq. 11, Gradshteyn–Ryzhik 6.612.3), ∫_0^∞ J_μ(xu)J_μ(xu')$e^{{−x|z|}}$ dx = (1/π)√(uu') Q_{μ−1/2}(Z), with Z=($u^{2}$+u'^2+$z^{2}$)/(2uu'). Applied at μ=0, it converts the iterated integral in Eq. (10) into a single quadrature. The second ingredient is the reduction (Byrd–Friedman 560.01) of the Legendre function Q_{−1/2}(Z) to √m K(m), with m=4u'u/((u+u')^2+$z^{2}$), which produces the elliptic-integral form of Eq. (17). The normalization of σ0 via the $\\beta$ function fixes the total charge Q and sets the prefactor Ṽ_o.","core_discovery":"Starting from Coulomb's law, the authors write the potential as a Bessel integral (Eq. 4) and, for the profile family σ(u)=σ0(1−$u^{2}$)^ν u^p, reduce it to the exact representation Φ_{p,ν}(ρ,z) = Ṽ_o ∫$_0^{1}$ (1−u'^2)^ν u'^{p+1} K(4u'u/((u+u')^2+$z^{2}$)) / $\\sqrt$((u+u')^2+$z^{2}$) du' (Eq. 17) for all real ν > −1 and p > −1. The uniform disk (ν = p = 0) and the edge-concentrated equipotential disk (ν = −1/2, p = 0) emerge as limiting cases matching earlier results. On the disk plane with ν = −1/2, the normalized potential V_p(u) is flat for p = 0, monotonically decreasing for −1 < p < 0, and non-monotonic with an edge-side maximum for p > 0, with an estimated crossover at p_c ≈ −0.34.","pith_inferences":["If the formula is exact for all ν,p > −1, then profiles with simultaneous center and edge singularities (e.g., ν = −1/2, p = −1/2) are also captured; the paper plots this case but does not emphasize that the integrand diverges at both endpoints while the integral remains finite — a useful stress test for numerical implementations.","The same identity chain should generalize to other axisymmetric geometries, such as annular disks or finite cylinders, where the radial integration limits change but the Bessel-product identity applies unchanged; the authors do not discuss this extension.","The estimated critical exponent p_c ≈ −0.34 likely corresponds to a condition on the derivative of the in-plane potential at the rim or center; deriving an analytic equation for p_c from Eq. (17) would remove the reliance on numerical estimation.","The method's reliance on the interchange of the q- and u'-integrals suggests the formula's domain might be narrower than the stated ν,p range if convergence is only conditional; testing at the boundary values (p → −1⁺ or ν → −1⁺) would map the true domain."],"forward_implications":["Off-axis potentials for the entire profile family (uniform, edge-concentrated, center-depleted) are computable with a single one-dimensional quadrature over standard functions.","The known benchmark results — the uniform disk of Ref. [12] and the equipotential edge-concentrated disk of Ref. [15] — are recovered as special cases, validating the formula.","The crossover exponent p_c ≈ −0.34 marks a qualitative change in the in-plane potential: for p > p_c a maximum exists near the rim, while for p < p_c the potential decreases monotonically from the center.","Because the formula is exact in u and z and depends only on ν and p, it provides a fast analytic handle for modeling engineered surface-charge profiles in electrostatics applications."],"supporting_citations":[{"why":"Supplies the classical identity 6.612.3 that collapses the product of Bessel functions times e^{−x|z|} to a Legendre function, the core reduction of the paper.","marker":"[40]"},{"why":"Supplies formula 560.01 expressing Q_{−1/2}(Z) as √m K(m), converting the Legendre-function form into the elliptic-integral form of Eq. (17).","marker":"[41]"},{"why":"Provides the uniform-disk potential benchmark that the authors recover by setting ν=p=0.","marker":"[12]"},{"why":"Provides the equipotential edge-concentrated disk benchmark recovered at ν=−1/2, p=0, and comparison for the crossover analysis.","marker":"[15]"},{"why":"Provides the cylindrical addition theorem (Bessel expansion of the Coulomb kernel) used to derive the starting integral representation Eq. (4).","marker":"[38]"}],"fun_headline_variants":["One integral unifies all disk charge profiles","Exact disk potential for any radial charge shape","A single formula for arbitrary disk charge patterns","Disk electrostatics solved with one integral"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper passes from the iterated integral in Eq. (10) to the single integral in Eq. (14) by interchanging the k-integration and the u'-integration, and assumes this interchange (and the convergence of the double integral) is valid for the entire advertised range ν > −1, p > −1 without stating or proving the uniformity conditions.","fun_headline_variants_meta":{"raw":{"variants":["One integral unifies all disk charge profiles","Exact disk potential for any radial charge shape","A single formula for arbitrary disk charge patterns","Disk electrostatics solved with one integral"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000196,"raw_usage":{"total_tokens":1351,"prompt_tokens":924,"completion_tokens":427,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":540,"completion_tokens_details":{"reasoning_tokens":371}},"tokens_in":540,"tokens_out":427,"duration_ms":5299,"temperature":1.0,"reasoning_tokens":371,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T10:14:35.679601+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Directly evaluate the double integral in Eq. (10) by numerical quadrature over q and u' for a parameter set near the reported crossover, say ν = −1/2, p = −0.34, at z = 0 and u = 0.5, and compare the result with the single-integral formula Eq. (17) to machine precision; any difference beyond roundoff that persists under adaptive refinement would falsify the exactness claim for that sub-range. A second check is to test endpoint convergence with p → −1⁺ or ν → −1⁺, where the integrand becomes singular at u' = 0 or u' = 1.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the classical identity 6.612.3 that collapses the product of Bessel functions times e^{−x|z|} to a Legendre function, the core reduction of the paper."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies formula 560.01 expressing Q_{−1/2}(Z) as √m K(m), converting the Legendre-function form into the elliptic-integral form of Eq. (17)."},{"cited_title":"Ciftja and I","cited_arxiv_id":null,"evidence_quote":"Provides the uniform-disk potential benchmark that the authors recover by setting ν=p=0."},{"cited_title":"Ciftja,Results for charged disks with differente forms of surface charge density,Results Phys.16, 102962 (2020)","cited_arxiv_id":null,"evidence_quote":"Provides the equipotential edge-concentrated disk benchmark recovered at ν=−1/2, p=0, and comparison for the crossover analysis."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the cylindrical addition theorem (Bessel expansion of the Coulomb kernel) used to derive the starting integral representation Eq. (4)."}],"review_version":1}