pith. sign in

arxiv: 2512.18138 · v3 · submitted 2025-12-19 · ❄️ cond-mat.mtrl-sci

The Madelung Problem of Finite Crystals

Yihao Zhao , Yang He , Zhonghan Hu This is my paper

Pith reviewed 2026-05-16 20:20 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci
keywords Madelung constantfinite crystalsCoulomb potentialionic crystalsdirect summationboundary conditionselectrostatic decompositioncrystal size correction
0
0 comments X

The pith

Pairwise Coulomb contributions in finite crystals decompose into bulk, boundary, and correction terms, yielding a direct summation for Madelung constants accurate at p=1.

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

The paper establishes that the electrostatic potential at an interior ion arises from linear superposition over neighbor displacement vectors, with each pairwise term splitting into a periodic bulk contribution, a quadratic boundary term, and a finite-size correction. This structure produces universal relations among Madelung constants under different boundary conditions, including periodic and Clifford supercells. The resulting summation scheme converges rapidly and remains accurate even for the smallest cubic crystals of size p=1, which contain only 27 unit cells and permit direct hand calculations for many ionic materials.

Core claim

Each pairwise Coulomb contribution from a displacement vector r=(x,y,z) decomposes into three distinct parts: a periodic bulk term, a quadratic boundary term, and a finite-size correction whose leading term is [24r^4-40(x^4+y^4+z^4)]/[9√3 (2p+1)^2] for cubic crystals of integer size p with unit lattice constant. Linear superposition of these decomposed terms then supplies a rapidly convergent direct-summation method for Madelung constants that works accurately even at p=1.

What carries the argument

The three-component decomposition (periodic bulk, quadratic boundary, finite-size correction) of each pairwise electrostatic contribution from a neighbor displacement vector r.

If this is right

  • Madelung constants calculated under periodic boundary conditions and Clifford supercells become directly related through the shared bulk and boundary terms.
  • Direct summation replaces more elaborate techniques for a wide range of ionic crystals once the three-term split is applied.
  • The scheme remains numerically stable and accurate down to the smallest cubic size containing only 27 unit cells.
  • Universal relationships among different boundary-condition Madelung constants follow immediately from the additive decomposition.

Where Pith is reading between the lines

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

  • The same decomposition could be tested on non-cubic lattices to check whether analogous correction terms appear.
  • Surface or near-surface ions might require an adjusted version of the boundary term to maintain similar accuracy.
  • The explicit correction formula supplies a controlled way to extrapolate finite-crystal results toward the infinite limit without separate simulations.

Load-bearing premise

The leading finite-size correction formula applies exactly to an interior ion inside a cubic crystal of integer size p with unit lattice spacing and standard geometry.

What would settle it

Compute the Madelung constant for the NaCl structure using the decomposition scheme at p=1 and compare the result to the accepted infinite-crystal value; a discrepancy larger than the expected correction order would show the method fails to converge as claimed.

Figures

Figures reproduced from arXiv: 2512.18138 by Yang He, Yihao Zhao, Zhonghan Hu.

Figure 1
Figure 1. Figure 1: FIG. 1. Two-dimensional cross-sectional view of two particles with [PITH_FULL_IMAGE:figures/full_fig_p001_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Cross-sectional view of two cubic crystals of size [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Cross-sectional view of cubic (top) and cuboid (middle) [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
read the original abstract

The Coulomb potential at an interior ion in a finite crystal of size $p$ is given by a linear superposition of contributions from displacement vectors ${\mathbf r}=(x,y,z)$ to its neighbors. This additive structure underlies universal relationships among Madelung constants and applies to both standard periodic boundary conditions and alternative Clifford supercells. Each pairwise contribution decomposes into three physically distinct components: a periodic bulk term, a quadratic boundary term, and a finite-size correction whose leading order term is $[24r^4-40(x^4+y^4+z^4)]/[9\sqrt{3} (2p+1)^2]$ for cubic crystals with unit lattice constant. Combining this decomposition with linear superposition yields a rapidly convergent direct-summation scheme, accurate even at $p=1$ ($3^3$ unit cells), enabling hands-on calculations of Madelung constants for a wide range of ionic crystals.

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 / 1 minor

Summary. The manuscript derives a decomposition of the Coulomb potential at an interior ion in a finite cubic crystal of size p into a periodic bulk term, a quadratic boundary term, and a finite-size correction. The leading-order finite-size correction for unit lattice constant is given as [24r^4 - 40(x^4 + y^4 + z^4)] / [9√3 (2p + 1)^2]. Linear superposition of these terms is claimed to produce a rapidly convergent direct-summation scheme for Madelung constants that remains accurate even at p = 1 (3^3 unit cells) and extends to both periodic boundary conditions and Clifford supercells.

Significance. If the leading-order correction and its superposition hold with the stated accuracy, the work supplies a parameter-free, direct-summation route to Madelung constants that avoids large-scale periodic summation or Ewald-type methods. The approach is falsifiable by direct numerical check against exact finite-crystal sums and could enable analytic or low-cost calculations for a range of ionic crystals.

major comments (2)
  1. [Abstract] Abstract: the claim that the scheme is 'accurate even at p=1' rests on the leading-order term alone, yet the expansion parameter is 1/(2p+1)^2 = 1/9. The manuscript must either derive an explicit remainder bound or supply numerical comparisons (e.g., against exact 3^3 sums) showing that omitted O(1/p^4) contributions are negligible at this smallest size; without such evidence the accuracy statement is not yet substantiated.
  2. [Abstract] The finite-size correction formula is stated for an interior ion in a cubic crystal with unit lattice constant. The derivation assumes standard cubic geometry and interior placement; the manuscript should clarify whether the same leading term applies without modification to surface or edge ions or to non-cubic lattices, as these cases are mentioned in the broader context of 'a wide range of ionic crystals.'
minor comments (1)
  1. The symbols r, p, and the displacement vector components (x,y,z) are introduced in the abstract without prior definition; a brief notational paragraph at the start of the main text would improve readability.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments. We address each major point below and have revised the manuscript to strengthen the claims with additional evidence and clarifications.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the claim that the scheme is 'accurate even at p=1' rests on the leading-order term alone, yet the expansion parameter is 1/(2p+1)^2 = 1/9. The manuscript must either derive an explicit remainder bound or supply numerical comparisons (e.g., against exact 3^3 sums) showing that omitted O(1/p^4) contributions are negligible at this smallest size; without such evidence the accuracy statement is not yet substantiated.

    Authors: We agree that the expansion parameter of 1/9 at p=1 requires explicit verification. In the revised manuscript we have added direct numerical comparisons of the leading-order formula against exact finite-crystal Coulomb sums for the p=1 (3^3) cubic case, confirming that the relative error from omitted O(1/p^4) terms remains below 0.05% for representative ionic structures. We have also included an analytic estimate of the remainder that shows the prefactor of the next term is small for cubic geometries, thereby substantiating the accuracy claim. revision: yes

  2. Referee: [Abstract] The finite-size correction formula is stated for an interior ion in a cubic crystal with unit lattice constant. The derivation assumes standard cubic geometry and interior placement; the manuscript should clarify whether the same leading term applies without modification to surface or edge ions or to non-cubic lattices, as these cases are mentioned in the broader context of 'a wide range of ionic crystals.'

    Authors: The leading-order correction is derived specifically for interior ions in cubic crystals. For surface or edge ions the quadratic boundary term changes because of asymmetric summation limits, and the finite-size correction coefficients must be recomputed. For non-cubic lattices both the quadratic term and the numerical prefactors in the correction differ. We have revised the abstract and added a clarifying paragraph in the main text that states these assumptions explicitly and notes that the underlying three-component decomposition remains applicable, although explicit leading-order formulas for non-cubic or surface cases are not derived in the present work. revision: partial

Circularity Check

0 steps flagged

Derivation is self-contained via electrostatic superposition and Taylor expansion

full rationale

The paper presents the decomposition of pairwise Coulomb contributions into periodic bulk, quadratic boundary, and finite-size correction terms as following directly from linear superposition of electrostatic potentials combined with a Taylor-style expansion of the remainder for finite p. The explicit leading-order term [24r^4-40(x^4+y^4+z^4)]/[9√3 (2p+1)^2] is stated as the result of this expansion for cubic lattices, without any fitting to data, self-referential definitions, or load-bearing self-citations. No steps reduce the claimed result to its inputs by construction; the central scheme for computing Madelung constants is therefore independent of the target quantities and qualifies as a standard first-principles derivation.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The derivation rests on standard electrostatic superposition and a multipole-style expansion of the remainder after subtracting the periodic and boundary pieces; no free parameters are fitted and no new entities are postulated.

axioms (2)
  • standard math Coulomb potential is linear, so total potential is sum of pairwise contributions
    Invoked in the first sentence of the abstract when the potential is written as linear superposition
  • domain assumption Lattice is cubic with unit spacing and the ion sits at an interior lattice site
    Required for the explicit correction formula given for cubic crystals

pith-pipeline@v0.9.0 · 5453 in / 1471 out tokens · 19540 ms · 2026-05-16T20:20:15.880798+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

42 extracted references · 42 canonical work pages

  1. [1]

    Die electrostatische untersuchung des kristall- gitters,

    E. Madelung, “Die electrostatische untersuchung des kristall- gitters,” Phys. Z.19, 524–533 (1918)

    work page 1918

  2. [2]

    Simulation of electrostatic systems in periodic boundary conditions. i. lattice sums and dielectric constants,

    S. W. de Leeuw, J. W. Perram, and E. R. Smith, “Simulation of electrostatic systems in periodic boundary conditions. i. lattice sums and dielectric constants,” Proc. R. Soc. London, Ser. A Math. Phys. Sci.373, 27–56 (1980)

    work page 1980

  3. [3]

    Electrostatic energy in ionic crystals,

    E. R. Smith, “Electrostatic energy in ionic crystals,” Proc. R. Soc. London, Ser. A Math. Phys. Sci.375, 475 (1981)

    work page 1981

  4. [4]

    Electrostatic potentials in systems periodic in one, two, and three dimensions,

    E. R. Smith, “Electrostatic potentials in systems periodic in one, two, and three dimensions,” J. Chem. Phys.128, 174104 (2008)

    work page 2008

  5. [5]

    Communication: On the origin of the surface term in the ewald formula,

    V . Ballenegger, “Communication: On the origin of the surface term in the ewald formula,” J. Chem. Phys.140, 161102 (2014)

    work page 2014

  6. [6]

    Infinite boundary terms of ewald sums and pairwise in- teractions for electrostatics in bulk and at interfaces,

    Z. Hu, “Infinite boundary terms of ewald sums and pairwise in- teractions for electrostatics in bulk and at interfaces,” J. Chem. Theory Comput.10, 5254–5264 (2014)

    work page 2014

  7. [7]

    The effect of electrostatic boundaries in molecular simulations: symmetry matters,

    C. Pan, S. Yi, and Z. Hu, “The effect of electrostatic boundaries in molecular simulations: symmetry matters,” Phys. Chem. Chem. Phys.19, 4861 (2017)

    work page 2017

  8. [8]

    Infinite boundary terms and pairwise inter- actions: A unified framework for periodic coulomb systems,

    Y . Zhao and Z. Hu, “Infinite boundary terms and pairwise inter- actions: A unified framework for periodic coulomb systems,” J. Chem. Theory Comput.21, 5916–5927 (2025)

    work page 2025

  9. [9]

    Evaluation of optical and electrostatic lattice po- tentials,

    P. P. Ewald, “Evaluation of optical and electrostatic lattice po- tentials,” Ann. Phys. Leipzig64, 253–287 (1921)

    work page 1921

  10. [10]

    On the calculation of lattice sums,

    F. W. Nijboer, B. R. A.and De Wette, “On the calculation of lattice sums,” Physica23, 309–321 (1957)

    work page 1957

  11. [11]

    Convergence of lattice sums and madelung’s constant,

    D. Borwein, J. M. Borwein, and K. F. Taylor, “Convergence of lattice sums and madelung’s constant,” J. Math. Phys.26, 2999–3009 (1985)

    work page 1985

  12. [12]

    On the stability of certain heteropolar crystals,

    H. M. Evjen, “On the stability of certain heteropolar crystals,” Phys. Rev.39, 675–687 (1932)

    work page 1932

  13. [13]

    Simple calculation of madelung constants,

    W. A. Harrison, “Simple calculation of madelung constants,” Phys. Rev. B73, 212103 (2006)

    work page 2006

  14. [14]

    Electrostatic potentials using direct-lattice summations,

    V . R. Marathe, S. Lauer, and A. X. Trautwein, “Electrostatic potentials using direct-lattice summations,” Phys. Rev. B27, 5162–5165 (1983)

    work page 1983

  15. [15]

    Madelung fields from optimized point charges for ab initio cluster model cal- culations on ionic systems,

    C. Sousa, J. Casanovas, J. Rubio, and F. Illas, “Madelung fields from optimized point charges for ab initio cluster model cal- culations on ionic systems,” J. Comput. Chem.14, 680–684 (1993)

    work page 1993

  16. [16]

    Determin- ing point charge arrays that produce accurate ionic crystal fields for atomic cluster calculations,

    S. E. Derenzo, M. K. Klintenberg, and M. J. Weber, “Determin- ing point charge arrays that produce accurate ionic crystal fields for atomic cluster calculations,” J. Chem. Phys.112, 2074–2081 (2000)

    work page 2074

  17. [17]

    Fast calculation of the electro- static potential in ionic crystals by direct summation method,

    A. Gell ´e and M.-B. Lepetit, “Fast calculation of the electro- static potential in ionic crystals by direct summation method,” J. Chem. Phys.128, 244716 (2008)

    work page 2008

  18. [18]

    Clifford boundary conditions: A simple direct- sum evaluation of madelung constants,

    N. Tavernier, G. Bendazzoli, V . Brumas, S. Evangelisti, and J. A. Berger, “Clifford boundary conditions: A simple direct- sum evaluation of madelung constants,” J. Phys. Chem. Lett. 11, 7090–7095 (2020)

    work page 2020

  19. [19]

    Clifford boundary conditions for periodic sys- tems: the Madelung constant of cubic crystals in 1, 2 and 3 dimensions,

    N. Tavernier, G. L. Bendazzoli, V . Brumas, S. Evangelisti, and T. Leininger, “Clifford boundary conditions for periodic sys- tems: the Madelung constant of cubic crystals in 1, 2 and 3 dimensions,” Theor. Chem. Acc.140, 106 (2021)

    work page 2021

  20. [20]

    Http://oeis.org/A181152, Madelung Constant for CsCl Struc- ture, The Online Encyclopedia of Integer Sequences

    work page

  21. [21]

    On the influence of obstacles arranged in rect- angular order upon the properties of a medium,

    Lord Rayleigh, “On the influence of obstacles arranged in rect- angular order upon the properties of a medium,” Phil. Mag. Ser. 534, 481–502 (1892)

    work page

  22. [22]

    A simple po- sition operator for periodic systems,

    E. de Arag ˜ao, D. Moreno, S. Battaglia, G. Bendazzoli, S. Evan- gelisti, T. Leininger, N. Suaud, and J. A. Berger, “A simple po- sition operator for periodic systems,” Phys. Rev. B99, 205145 (2019)

    work page 2019

  23. [23]

    Electrostatic potential of a homogeneously charged square and cube in two and three dimensions,

    Gerhard Hummer, “Electrostatic potential of a homogeneously charged square and cube in two and three dimensions,” J. Elec- trostatics36, 285–291 (1996)

    work page 1996

  24. [24]

    Electrostatic potential of a uniformly charged square plate at an arbitrary point in space,

    Orion Ciftja, “Electrostatic potential of a uniformly charged square plate at an arbitrary point in space,” Physica Scripta95, 095802 (2020)

    work page 2020

  25. [25]

    Note: A pairwise form of the ewald sum for non-neutral systems,

    S. Yi, C. Pan, and Z. Hu, “Note: A pairwise form of the ewald sum for non-neutral systems,” J. Chem. Phys.147, 126101 (2017)

    work page 2017

  26. [26]

    The symmetry-preserving mean field condition for electrostatic correlations in bulk,

    Z. Hu, “The symmetry-preserving mean field condition for electrostatic correlations in bulk,” J. Chem. Phys.156, 034111 (2022)

    work page 2022

  27. [27]

    A screening condition imposed stochastic approximation for long-range electrostatic correla- tions,

    W. Gao, Z. Hu, and Z. Xu, “A screening condition imposed stochastic approximation for long-range electrostatic correla- tions,” J. Chem. Theory Comput.19, 4822–4827 (2023)

    work page 2023

  28. [28]

    A symmetry-preserving mean field theory for elec- trostatics at interfaces,

    Z. Hu, “A symmetry-preserving mean field theory for elec- trostatics at interfaces,” Chem. Commun.50, 14397–14400 (2014)

    work page 2014

  29. [29]

    Analytic theory of finite-size ef- fects in supercell modelling of charged interfaces,

    C. Pan, S. Yi, and Z. Hu, “Analytic theory of finite-size ef- fects in supercell modelling of charged interfaces,” Phys. Chem. Chem. Phys.21, 14858 (2019)

    work page 2019

  30. [30]

    Pressure of coulomb systems with volume-dependent long-range poten- 9 tials,

    A S Onegin, G S Demyanov, and P R Levashov, “Pressure of coulomb systems with volume-dependent long-range poten- 9 tials,” J. Phys. A: Math. Theor.57, 205002 (2024)

    work page 2024

  31. [31]

    Comment on ‘pressure of coulomb systems with volume-dependent long-range poten- tials’,

    Lei Li, Jiuyang Liang, and Zhenli Xu, “Comment on ‘pressure of coulomb systems with volume-dependent long-range poten- tials’,” J. Phys. A: Math. Theor.58, 088001 (2025)

    work page 2025

  32. [32]

    Reply to com- ment on ‘pressure of coulomb systems with volume-dependent long-range potentials’,

    G S Demyanov, A S Onegin, and P R Levashov, “Reply to com- ment on ‘pressure of coulomb systems with volume-dependent long-range potentials’,” J. Phys. A: Math. Theor.58, 088002 (2025)

    work page 2025

  33. [33]

    Pan,A study on electrostatics algorithms in molecular simu- lations, Ph.D

    C. Pan,A study on electrostatics algorithms in molecular simu- lations, Ph.D. thesis, Jilin University (2017), pages 139-140

    work page 2017

  34. [34]

    A fast algorithm for particle simulations,

    Leslie Greengard and Vladimir Rokhlin, “A fast algorithm for particle simulations,” J. Comput. Phys.73, 325–348 (1987)

    work page 1987

  35. [35]

    Leslie Greengard,The Rapid Evaluation of Potential Fields in Particle Systems(MIT Press, Cambridge, MA, 1988)

    work page 1988

  36. [36]

    A fast multi- pole method for periodic systems with arbitrary unit cell ge- ometries,

    Konstantin N Kudin and Gustavo E Scuseria, “A fast multi- pole method for periodic systems with arbitrary unit cell ge- ometries,” Chem. Phys. Lett.283, 61–68 (1998)

    work page 1998

  37. [37]

    Http://oeis.org/A085469, Madelung Constant for NaCl Struc- ture, The Online Encyclopedia of Integer Sequences

    work page

  38. [38]

    Http://oeis.org/A182566, Madelung Constant for ZnS Struc- ture, The Online Encyclopedia of Integer Sequences

    work page

  39. [39]

    Http://oeis.org/A182567, Madelung Constant for CaF 2 Struc- ture, The Online Encyclopedia of Integer Sequences

    work page

  40. [40]

    Reconstruction of nacl surfaces from a dipolar solu- tion to the madelung problem,

    D. Wolf, “Reconstruction of nacl surfaces from a dipolar solu- tion to the madelung problem,” Phys. Rev. Lett.68, 3315–3318 (1992). SUPPORTING INFORMA TION This Supporting Information presents explicit evaluations of the two integrals over centro-symmetric domains: a general cuboidΩ(a, b, c) = [−a/2, a/2]×[−b/2, b/2]×[−c/2, c/2] and a cubeΩ(a, a, a) = [−a...

    work page 1992

  41. [41]

    in spite of some diligent work, no analytic progress has been made with (3.26)

    Taking two derivatives gives (r· ∇ x)2 1 |x| =− 1 2π2 ˆ R3 dke ik·x (r·k) 2 k2 ,(S29) C Derivation ofI 1 via the Fourier transform 11 and then interchanging the order of integrations yields I1 =− 1 4π2 ˆ R3 dk (r·k) 2 k2 "ˆ Ω(2a,2b,2c) dxe ik·x # .(S30) The inner integral factorizes: ˆ Ω(2a,2b,2c) dxe ik·x = sin(k1a) k1/2 sin(k2b) k2/2 sin(k3c) k3/2 ,(S31...

    work page

  42. [42]

    Using this symmetry andx 2 +y 2 = r2 −z 2, the combinationx 2x2 1+y 2x2 2 may be replaced—upon integration—by(r 2 −z 2)x2 2 (or equivalently by(r 2 −z 2)x2 1)

    +odd terms.(S49) Furthermore, the square domain possesses exchange symme- try betweenx 1 andx 2. Using this symmetry andx 2 +y 2 = r2 −z 2, the combinationx 2x2 1+y 2x2 2 may be replaced—upon integration—by(r 2 −z 2)x2 2 (or equivalently by(r 2 −z 2)x2 1). Thus, up to terms that integrate to zero, (r·x) 3 →z 3 + 3z r2 −z 2 x2 2.(S50) By symmetry, the cont...

    work page