Pith. sign in

REVIEW 3 major objections 4 minor 41 references

Modeling the Electrostatic Potential of Disks with Arbitrary Radial Charge Profiles

T0 review · 3 major / 4 minor · reviewed 2026-08-07 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read The paper derives an exact single-integral representation, Eq.

desk verdict 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. read the letter →

arxiv 2506.08040 v1 pith:JCJUJTR5 submitted 2025-06-06 physics.class-ph

Jose Ricardo de Sousa , Marcio Gomes da Silva , Orion Ciftja This is my paper
classification physics.class-ph PACS 41.20.Cv
keywords electrostaticpotentialchargeddiskellipticintegralsBesselfunctionsradialchargedistributionaxialsymmetryequipotentialLegendre
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

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.

What carries the argument

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.

What would settle it

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.

Watch

Extended reading notes

Core claim

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.

Load-bearing premise

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.

Editorial extensions

If this is right

  • 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.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • 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.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 4 minor

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.

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 (3)
  1. [Eq. (11) and Eq. (14)] 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.
  2. [Eqs. (10)-(17), Section 3] 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.
  3. [Section 3, p_c approximately -0.34] 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.
minor comments (4)
  1. [Eqs. (7)-(9) and (18)] 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!.
  2. [Throughout] There are numerous typos, including 'propertied', 'mathemarical', 'baxkground', 'Rigdly', 'Boudaries', and 'different forms'. Please proofread the manuscript carefully.
  3. [References] Reference [39] (Sonine) is listed but not cited in the text.
  4. [Eq. (17) and Eq. (18)] The prefactor in Eq. (17) is typeset as 'eVo', while Eq. (18) defines tilde V_o; please unify the notation.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: Eq. (17) is derived from Coulomb's law and standard special-function identities; benchmarks and the critical p_c are outputs, not fitted inputs.

full rationale

The derivation chain is self-contained and non-circular. It starts from the fundamental Coulomb potential Eq. (1), inserts the Bessel-function addition theorem Eq. (2), integrates the azimuthal angle to obtain the general representation Eq. (4), then specializes to the charge profile sigma(u) = sigma_0 (1-u^2)^nu u^p, Eq. (5). The normalization constant sigma_0 is determined by the total charge condition, Eqs. (6)-(9), which is an independent physical constraint and not the potential itself. The main result, Eq. (17), is obtained by substituting the Gradshteyn-Ryzhik identity Eq. (11) into Eq. (10), then using the Byrd-Friedman relation Eq. (15) to convert the Legendre function Q_{-1/2} into the complete elliptic integral K. These identities are external, standard mathematical results; they do not assume the final potential formula. The two limiting cases, Eqs. (19) and (20), are presented as special cases of Eq. (17) and are used as consistency checks; even though Refs. [12] and [15] include the present author, these citations are not used to justify the derivation and no load-bearing argument reduces to them. The threshold p_c ≈ -0.34 is a numerical output of the derived formula, not a fitted parameter, so no fitted input is being renamed as a prediction. The paper's possible gap regarding interchange of integrations or the z -> 0 limit is a rigor/correctness concern, not a circularity concern; nothing in the paper defines its target in terms of itself.

Assumptions & free parameters 2 free parameters · 5 assumptions · 0 invented entities

The shape exponents nu and p are freely chosen inputs that define the charge family; they are not fitted to external measurements. The derivation relies on three standard special-function identities from handbooks plus a convergence and interchange assumption at the singular endpoints. No new physical entities are introduced.

free parameters (2)
  • nu
    Shape exponent controlling edge behavior in the charge model; an input parameter of the family, not fitted to data.
  • p
    Shape exponent controlling center behavior in the charge model; an input parameter, not fitted to data.
assumptions (5)
  • standard math Bessel addition theorem expansion of the Coulomb kernel, Eq (2).
    Used without proof to convert the surface integral over the disk into a Bessel integral; it is a textbook identity for the free-space Green's function in cylindrical coordinates.
  • standard math Gradshteyn-Ryzhik identity 6.612.3, Eq (11).
    Replaces the product of two Bessel functions integrated against a decaying exponential by the Legendre function Q. The paper cites the handbook rather than proving the identity.
  • standard math Byrd and Friedman relation Q_{-1/2}(Z)=sqrt(m)K(m), Eq (15).
    Converts the Legendre function to a complete elliptic integral of the first kind. Cited as formula 560.01 of Ref [41].
  • domain assumption Total charge finite requires nu>-1 and p>-1, Eq (6).
    The normalization constant Eq (9) is derived from a beta-function integral that converges only in this range; the paper restricts the family accordingly.
  • domain assumption Interchange of the k-integration and u'-integration in Eqs. (10)-(14).
    The exact formula Eq (17) follows only if the iterated integral can be collapsed via Eq (11); the paper asserts validity for the whole parameter range without proving uniform convergence.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Modeling the Electrostatic Potential of Disks with Arbitrary Radial Charge Profiles." pith.science (2026). https://pith.science/paper/JCJUJTR5

@misc{pith2026250608040,
  author       = {Pith},
  title        = {Pith review of: Modeling the Electrostatic Potential of Disks with Arbitrary Radial Charge Profiles},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/JCJUJTR5}},
  note         = {Machine review of arXiv:2506.08040}
}
read the original abstract

In this work, we present an analytical study of the electrostatic potential generated by a charged disk with a surface charge distribution that possesses radial axial symmetry. We express the potential in cylindrical coordinates and apply a Bessel function representation for the Coulomb kernel to simplify its calculation. We consider a broad family of surface charge densities and derive an exact expression for the potential based on a classical integral identity involving Bessel functions. This general formulation includes several important special cases such as the uniformly charged disk and edge-concentrated charge density distributions. We explore how variations in various parameters affect the resulting potential, thereby illustrating the influence of radial charge modulation. It is found that the potential undergoes a qualitative change of behavior when the chosen set of parameters is varied. The results are analytically elegant and computationally efficient as all functions involved are implemented in standard scientific computing libraries.

Figures

Figures reproduced from arXiv: 2506.08040 by the authors.

Figure 1
Figure 1. Behavior of the dimensionless surface charge density, [PITH_FULL_IMAGE:figures/full_fig_p008_1.png] view at source ↗
Figure 2
Figure 2. Dimensionless electrostatic potential in the plane [PITH_FULL_IMAGE:figures/full_fig_p010_2.png] view at source ↗

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

41 extracted references · 41 canonical work pages

  1. [1]

    D. J. Griffiths,Introduction to Electrodynamics(Cambridge University Press, Cambridge, 2017)

    work page 2017

  2. [2]

    J. R. Reitz and F. J. Milford,Foundations of Electromagnetic Theory(Addison-Wesley, Reading, 1960)

    work page 1960

  3. [3]

    J. D. Jackson,Classical Electrodynamics(John Wiley & Sons, Nova Jersey, 1999)

    work page 1999

  4. [4]

    J. A. Stratton,Electromagnetic Theory(John Wiley & Sons, Hoboken, 2007)

    work page 2007

  5. [5]

    W. K. H. Panofsky and M. Phillips,Classical Electricity and Magnetism(Addison-Wesley, Reading, 1962). 11

    work page 1962

  6. [6]

    L. D. Landau, E. M. Lifshitz and L. P. Pitaevskii,Eletrodynamics of Continuous Media, Vol. 8 of Course of Theoretical Physics (Pergamon Press, Oxford, 1984)

    work page 1984

  7. [7]

    Eyges,The Classical Electrodynamic Field(Dover, New York, 2012)

    L. Eyges,The Classical Electrodynamic Field(Dover, New York, 2012)

    work page 2012

  8. [8]

    E. J. Konopinski,Eletromagnetic Fields and Relativistic Particles(McGraw-Hill, New York, 1981)

    work page 1981

Show all 41 references
  1. [9]

    F. R. Zypman,Off-axis Electric Field of a Ring of Charge,Am. J. Phys.74, 295 (2006)

    2006

  2. [10]

    Datta,Electric and Magnetic Fields from a Circular Coil using Elliptic Integrals,Phys

    S. Datta,Electric and Magnetic Fields from a Circular Coil using Elliptic Integrals,Phys. Educ. India24, 203 (2007)

    2007

  3. [11]

    Ciftja, A

    O. Ciftja, A. Babineaux and N. Hafeez,The Electrostatic Potential of a Uniformly Charged Ring,Eur. J. Phys.30, 623 (2009)

    2009

  4. [12]

    Ciftja and I

    O. Ciftja and I. Hysi,The Electrostatic Potential of a Uniformly Charged Disk as the Source of Novel Mathematical Identities,Appl. Math. Lett.24, 1919 (2011)

    2011

  5. [13]

    Bochko and Z

    V. Bochko and Z. K. Silagadze,On the Electrostatic Potential and Electric Field of a Uniformly Charged Disk,Eur. J. Phys.41, 045201 (2020)

    2020

  6. [14]

    Escalante,Electrostatic Potential and Electric Field in the z Axis of a non Centred Circular Charged Ring,Eur

    F. Escalante,Electrostatic Potential and Electric Field in the z Axis of a non Centred Circular Charged Ring,Eur. J. Phys.42, 065703 (2021)

    2021

  7. [15]

    Ciftja,Results for charged disks with differente forms of surface charge density,Results Phys.16, 102962 (2020)

    O. Ciftja,Results for charged disks with differente forms of surface charge density,Results Phys.16, 102962 (2020)

    2020

  8. [16]

    A. E. Sagaydak and Z. K. Silagadze,Electrostatic potential of a uniformly charged disk through Green ’s theorem,Eur. J. Phys.46, 015203 (2025)

    2025

  9. [17]

    Morgan and L

    T. Morgan and L. Morgan,The Gravitational Field of a Disk,Phys. Rev.188, 2544 (1969); Errata:Phys. Rev. D1, 3522 (1970)

    1969

  10. [18]

    F. T. Krogh, E. W. Ng and W. V. Snyser,The Gravitational Field of a Disk,Celestial Mech.26, 395 (1982)

    1982

  11. [19]

    Lass and L

    H. Lass and L. Blitzer,The Gravitational Potential Due to Uniform Disks and Rings, Celestial Mech.30, 225 (1983)

    1983

  12. [20]

    Neugebauer and R

    G. Neugebauer and R. Meinel,General Relativistic Gravitational Field of a Rigdly Rotating Disk of Dust: Axis Potential, Disk Metric, and Surface Massa Density,Phys. Rev. Lett.73, 2166 (1994)

    1994

  13. [21]

    J. T. Conway,Analytical Solution for the Newtonian Gravitational Field Induced by Matter within Axisymmetric Boudaries,Mon. Not. R. Astron. Soc.316, 540 (2000)

    2000

  14. [22]

    Elipe, E

    A. Elipe, E. Tresaco and A. Riaguas,Gravitational Potential of a Massive Disk,Adv. Astr. Sci.143, 843 (2009). 12

    2009

  15. [23]

    J. M. Hur´ e e F. Hersant,The Newtonian Potential of Thin Disks,Astr. Astrophys.531, A6 (2011)

    2011

  16. [24]

    Schulz,The Gravitational Force and Potential of the Finite Mestel Disk,The Astrophys

    E. Schulz,The Gravitational Force and Potential of the Finite Mestel Disk,The Astrophys. J.747, 106 (2012)

    2012

  17. [25]

    J. D. Weiss,Certain Aspects of the Gravitational Field of a Disk,Appl. Math.9, 1360 (2018)

    2018

  18. [26]

    A. J. Silenko and Y. A. Tsalkou,Quasi-Uniform Gravitational Field of a Disk Revisited, Int. J. Mod. Phys. A34, 1950228 (2019)

    2019

  19. [27]

    Tiang, X

    B. Tiang, X. Zheng, T. J. Kempa, Y. Fang, N. Yu, G. Yu, J. Huang, and C. M. Lieber, Coaxial Nanowires as Solar Cells and Nanoelectronic Power Sources,Nature449, 885 (2007)

    2007

  20. [28]

    Nduwimana and X

    A. Nduwimana and X. -Q. Wange,Charge Carrier Separation in Modulation Doped Coaxial Semiconductor Nanowires,Nano. Lett.9, 283 (2009)

    2009

  21. [29]

    Levin,Electrostatic Correlations: from Plasma to Biology,Rep

    Y. Levin,Electrostatic Correlations: from Plasma to Biology,Rep. Prog. Phys.65, 1577 (2002)

    2002

  22. [30]

    A. L. Kholodenko and A. L. Beyerlein,Time-Dependent Density-Functional Theory for Multicomponent Systems,Phys. Rev. A43, 3309 (1986)

    1986

  23. [31]

    R. D. Coalson and A. Duncan,Systematic Ionic Screening Theory of Macroions,J. Chem. Phys.97, 5653 (1992)

    1992

  24. [32]

    A. Agra, E. Trizak and L. Bocquet,The Interplay Between Screening Properties and Col- loid Anisotropy: Towards a Reliable pair Potential for Disc-Like Charged Particles,Eur. Phys. E15, 345 (2004)

    2004

  25. [33]

    Morf and B

    R. Morf and B. I. Halperin,Monte Carlo Evaluation of Trial Wave Functions for the Fractional Quantized Hall Effect: Disk Geometry,Phys. Rev. B33, 2221 (1986)

    1986

  26. [34]

    Xia,An Estimate of the Ground State Energy of the Fractional Quantum Hall Effect,J

    J. Xia,An Estimate of the Ground State Energy of the Fractional Quantum Hall Effect,J. Math. Phys.40, 150 (1999)

    1999

  27. [35]

    Ciftja,Monte Carlo Study of Bose Laughlin Wave Function for Filling Factors1/2, 1/4and1/6,Europhys

    O. Ciftja,Monte Carlo Study of Bose Laughlin Wave Function for Filling Factors1/2, 1/4and1/6,Europhys. Lett.74, 486 (2006)

    2006

  28. [36]

    Ciftja,Exact Results for Systems of Electrons in the Fractional Quantum Hall Regime II,PhysicaB4042244 (2009); Errata:Physica B406, 2054 (2011)

    O. Ciftja,Exact Results for Systems of Electrons in the Fractional Quantum Hall Regime II,PhysicaB4042244 (2009); Errata:Physica B406, 2054 (2011)

    2009

  29. [37]

    Ciftja,Analytic Wave Functions for the Half-Filled Lowest Landau Level,Int

    O. Ciftja,Analytic Wave Functions for the Half-Filled Lowest Landau Level,Int. J. Mod. B24, 3489 (2010). 13

    2010

  30. [38]

    G. B. Arfken and H. J. Weber,Mathematical Methods for Physicists, 3 a edi¸ c˜ ao (Academic Press, California, 1985)

    1985

  31. [39]

    N. J. Sonine,Recherches sur les fonctions cylindriques ei le d´ eveloppement des fonctions continues en s´ eries,Math. Ann.16, 1 (1880)

  32. [40]

    I. S. Gradshteyn and I. M. Ryzhik,Table of Integrals, Series and Products, Academic Press, New York, 2007

    2007

  33. [41]

    P. F. Byrd and M. D. Friedman,Handbook of Elliptic Integrals for Engineers and Scientists, 2nd ed., Springer-Verlag, New York–Berlin, 1971. 14

    1971

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.