pith. sign in

arxiv: 2505.13148 · v4 · pith:QQLJ277Rnew · submitted 2025-05-19 · ⚛️ physics.class-ph

High-Discretization Method of Moments for Capacitance Calculation: A Cube and a Hollow Cylinder

Haiyong Gu , Liyuan Huang , Peide Yang , Tianshu Luo , Han Dong This is my paper

Pith reviewed 2026-05-22 14:32 UTC · model grok-4.3

classification ⚛️ physics.class-ph
keywords method of momentscapacitance calculationcubehollow cylinderdiscretizationelectrostaticsnumerical methodssymmetry
0
0 comments X

The pith

The capacitance of a unit cube reaches a maximum of 73.519014 pF at 90 by 90 sub-areas per face in method of moments calculations, after which further refinement lowers the result.

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

This paper applies the method of moments to capacitance calculations for a cube and a hollow cylinder using very fine surface discretizations. For the unit cube each face is divided into as many as 600 by 600 sub-areas, with geometric symmetry and parallel computing used to keep the computation feasible. The resulting capacitance first rises with added sub-areas and then falls, reaching its highest value at 90 by 90 divisions. The same procedure applied to a hollow cylinder produces numbers that line up with both Lekner's theoretical formula and Cavendish's experimental measurements. A reader would care because the work shows that simply making the grid finer does not guarantee better accuracy in electrostatic numerical methods.

Core claim

By dividing each face of the unit cube into up to 600 by 600 sub-areas while fully exploiting geometric symmetry between sub-areas and using parallel computing, the method of moments produces a capacitance that increases to a peak of 73.519014 pF at 90 by 90 sub-areas and then decreases. This pattern indicates that higher accuracy cannot be obtained merely by increasing the number of discretized sub-areas indefinitely. When the same method is applied to a hollow cylinder the computed capacitances show good agreement with independent numerical solutions based on Lekner's theoretical formula and with Cavendish's experimental values.

What carries the argument

Method of moments applied to discretized conductor surfaces, with geometric symmetry exploitation to reduce the number of independent integrals needed for capacitance.

If this is right

  • The unit cube yields a reference capacitance of 73.519014 pF at the observed optimum discretization.
  • Capacitance results for the hollow cylinder align with Lekner's theoretical formula and Cavendish's experimental data.
  • Symmetry exploitation together with parallel computing makes high-discretization runs computationally practical.
  • Beyond the optimum number of sub-areas the computed capacitance declines rather than continuing to improve.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Round-off or accumulation of floating-point errors at very high discretization counts may be responsible for the observed drop in capacitance.
  • The non-monotonic behavior with grid density could appear in capacitance calculations for other symmetric conductors.
  • The reported reference values can serve as benchmarks for testing new numerical schemes in electrostatics.

Load-bearing premise

That fully exploiting geometric symmetry between sub-areas combined with parallel computing preserves numerical accuracy without introducing systematic bias or round-off effects.

What would settle it

Running the unit-cube calculation at a discretization of 200 by 200 sub-areas per face or finer and checking whether the capacitance value remains below 73.519014 pF or begins to increase again.

Figures

Figures reproduced from arXiv: 2505.13148 by Haiyong Gu, Han Dong, Liyuan Huang, Peide Yang, Tianshu Luo.

Figure 2
Figure 2. Figure 2: Detailed capacitance values are provided in Table V [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Capacitance of a Unit Cube (edge length: 1 m) with [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Surface Charge Density Distribution on a Cube (edge [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: The surface of a hollow cylinder is divided into [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: The Charge Density Distribution Curve of the [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7: Schematic diagram of the electric potential [PITH_FULL_IMAGE:figures/full_fig_p007_7.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: Schematic diagram of the electric potential [PITH_FULL_IMAGE:figures/full_fig_p007_6.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8: Proportion of Positive [PITH_FULL_IMAGE:figures/full_fig_p009_8.png] view at source ↗
read the original abstract

This paper employs the method of moments (MOM) to calculate the capacitances of a cube and a hollow cylinder. For the cube, each face was divided into a maximum of 600 x 600 sub-areas. By fully exploiting the geometric symmetry between sub-areas and incorporating parallel computing, computational resources were significantly conserved. Our results show that the calculated capacitance of the cube first increases and then decreases as the number of sub-areas increases. When each face was divided into 90 x 90 sub-areas, the capacitance of the unit cube (with an edge length of 1 m) reached a maximum reference value of 73.519014 pF. This indicates that higher accuracy cannot be achieved merely by indefinitely increasing the number of discretized sub-areas. Subsequently, the method was applied to compute the capacitance of a hollow cylinder. The results were compared with numerical solutions based on Lekner's theoretical formula and Cavendish's experimental values, showing good agreement among the three.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript applies the method of moments (MOM) with high discretization (up to 600×600 sub-areas per face) to compute the capacitance of a unit cube (edge length 1 m) and a hollow cylinder. Geometric symmetry between sub-areas is fully exploited together with parallel computing to reduce resource demands. For the cube the computed capacitance is reported to increase with discretization level, reach a peak of 73.519014 pF at 90×90 sub-areas per face, and then decrease; the authors interpret the peak as a maximum reference value and conclude that further refinement cannot improve accuracy. For the hollow cylinder the MOM results are compared with Lekner’s theoretical formula and Cavendish’s experimental data and stated to show good agreement.

Significance. The computational strategy of symmetry reduction plus parallelization enables unusually fine discretizations and is a practical contribution to numerical electrostatics. If the non-monotonic cube result is shown to be free of artifacts, the identification of an optimal patch density would be useful for similar integral-equation problems. The cylinder comparison supplies an independent check that lends credibility to the overall approach.

major comments (2)
  1. [Abstract] Abstract (cube results): the claim that the capacitance reaches a maximum reference value of 73.519014 pF at 90×90 sub-areas and that “higher accuracy cannot be achieved merely by indefinitely increasing the number of discretized sub-areas” is load-bearing for the central conclusion. In standard MOM for the cube capacitance integral equation the 1/r kernel yields a capacitance that approaches the continuum limit (literature value ≈73.51 pF) with monotonic or steadily decreasing error once singular self-terms are treated accurately; the reported reversal after 90×90 therefore requires explicit demonstration that symmetry exploitation and parallel summation preserve positive-definiteness and introduce no systematic bias or round-off accumulation.
  2. [Abstract] Abstract and cube section: no convergence plots, error bars, or description of singular-integral handling (self-term quadrature, symmetry-enforced matrix assembly) are supplied for the 90×90 and higher discretizations. Without these diagnostics it is impossible to distinguish a genuine discretization optimum from a numerical artifact, directly undermining the interpretation that the observed peak is the highest attainable accuracy.
minor comments (2)
  1. [Abstract] The term “maximum reference value” is ambiguous; clarify whether it is proposed as the converged physical capacitance or simply the highest numerical result obtained.
  2. A figure or table plotting capacitance versus number of sub-areas per face would make the increase-then-decrease behavior immediately visible and aid reproducibility.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive feedback on our manuscript. The points raised highlight important aspects of numerical validation for the method of moments implementation. We address each major comment below and outline the revisions we will make to strengthen the presentation of our results.

read point-by-point responses

  1. Referee: [Abstract] Abstract (cube results): the claim that the capacitance reaches a maximum reference value of 73.519014 pF at 90×90 sub-areas and that “higher accuracy cannot be achieved merely by indefinitely increasing the number of discretized sub-areas” is load-bearing for the central conclusion. In standard MOM for the cube capacitance integral equation the 1/r kernel yields a capacitance that approaches the continuum limit (literature value ≈73.51 pF) with monotonic or steadily decreasing error once singular self-terms are treated accurately; the reported reversal after 90×90 therefore requires explicit demonstration that symmetry exploitation and parallel summation preserve positive-definiteness and introduce no systematic bias or round-off accumulation.

    Authors: We acknowledge that the non-monotonic trend after the peak at 90×90 sub-areas per face deviates from the expected monotonic convergence in standard MOM applications and therefore merits explicit justification. In the revised manuscript we will expand the methods section to detail the symmetry exploitation algorithm (including how equivalent sub-area interactions are grouped and matrix entries are populated without duplication) and the parallel summation procedure. We will add a brief verification that the resulting system matrix remains positive definite (via eigenvalue checks on representative smaller matrices or condition-number monitoring) and discuss the floating-point precision used to mitigate round-off accumulation. These additions will support the claim that the observed peak is not an artifact of the computational strategy. revision: yes

  2. Referee: [Abstract] Abstract and cube section: no convergence plots, error bars, or description of singular-integral handling (self-term quadrature, symmetry-enforced matrix assembly) are supplied for the 90×90 and higher discretizations. Without these diagnostics it is impossible to distinguish a genuine discretization optimum from a numerical artifact, directly undermining the interpretation that the observed peak is the highest attainable accuracy.

    Authors: We agree that the current manuscript lacks sufficient diagnostic material to allow readers to evaluate the convergence behavior and numerical stability at the finer discretizations. In the revision we will insert a new figure showing capacitance versus number of sub-areas per face (covering the full range up to 600×600) together with any available uncertainty estimates derived from the solver tolerance. We will also add an explicit description of the singular self-term quadrature (including the specific integration technique used to handle the 1/r singularity) and how symmetry constraints are enforced during matrix construction. These changes will enable independent assessment of whether the 90×90 peak represents a genuine optimum. revision: yes

Circularity Check

0 steps flagged

No circularity: direct numerical MOM outputs compared to external references

full rationale

The paper applies the standard method of moments discretization to the capacitance integral equation for a cube and hollow cylinder. Capacitance values are obtained by direct numerical solution after symmetry reduction and parallel assembly; the reported peak of 73.519014 pF at 90×90 patches per face is an output of that solver, not a fitted parameter or quantity defined in terms of itself. The non-monotonic trend is presented as an empirical observation from the computation, and the conclusion that further refinement yields no gain is an interpretation of that trend rather than a deductive reduction. For the cylinder, results are cross-checked against the independent Lekner formula and Cavendish experiment. No equations, self-citations, or ansatzes are shown to be load-bearing in a way that collapses the central claims back to the inputs by construction. The derivation chain remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The paper rests on the standard integral-equation formulation of the method of moments for electrostatics and the assumption that square-patch discretization plus symmetry reduction is sufficient; no new free parameters or invented entities are introduced in the abstract.

axioms (1)
  • domain assumption The method of moments integral formulation accurately represents the electrostatic potential for the cube and hollow cylinder geometries when surfaces are discretized into square sub-areas.
    Invoked implicitly when applying MOM to compute capacitance from charge distributions.

pith-pipeline@v0.9.0 · 5716 in / 1359 out tokens · 47536 ms · 2026-05-22T14:32:19.382540+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.