pith. sign in

arxiv: 2601.21776 · v2 · submitted 2026-01-29 · ❄️ cond-mat.mtrl-sci · cond-mat.soft· physics.chem-ph· physics.class-ph· physics.comp-ph

Model density approach to Ewald summations

Chiara Ribaldone , Jacques Kontak Desmarais This is my paper

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

classification ❄️ cond-mat.mtrl-sci cond-mat.softphysics.chem-phphysics.class-phphysics.comp-ph
keywords Ewald summationmodel charge densitymultipole momentselectrostatic potentialperiodic systemsgallium arsenideconvergence accelerationfundamental gap
0
0 comments X

The pith

A model charge density cancels multipole moments of the crystal to accelerate Ewald sum convergence.

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

The paper introduces a model charge density constructed to cancel the multipole moments of the crystalline charge distribution up to any chosen order. This cancellation is used to speed convergence of the Ewald-type lattice sums that compute electrostatic potentials in periodic systems. A sympathetic reader would care because these sums are required for accurate energies and potentials in condensed-matter studies, and slow convergence restricts the size or precision of calculations that can be performed. The technique applies to bulk systems with arbitrary unit cells, in either classical or quantum settings, and with any choice of basis functions for the density. The authors confirm the acceleration numerically on the fundamental gap of bulk gallium arsenide.

Core claim

By introducing a compensating model charge density that matches and cancels the multipole moments of the actual crystalline charge distribution up to a desired order, the Ewald sums converge substantially faster while leaving the physical electrostatic potential unchanged.

What carries the argument

The model charge density, built to cancel multipole moments of the crystalline charge distribution up to a chosen order, thereby shortening the range of the real-space Ewald sum.

If this is right

  • Electrostatic potentials in periodic bulk systems can be evaluated with fewer terms in the lattice sums.
  • Properties such as the fundamental gap in semiconductors become accessible with higher numerical precision at the same computational effort.
  • The method works for both classical and quantum calculations that employ arbitrary unit cells and any representation of the charge density.
  • Existing implementations of Ewald sums in crystal codes can be clarified and improved by this construction.

Where Pith is reading between the lines

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

  • The same cancellation idea could be tested on surfaces or defects once a suitable model density for non-periodic boundaries is defined.
  • If the extra cost of building the model density stays low, the approach may become standard for large-scale molecular-dynamics runs of ionic materials.
  • Analogous model densities might accelerate other conditionally convergent sums, such as those appearing in dispersion or magnetic dipole interactions.

Load-bearing premise

A suitable model charge density can be constructed for arbitrary unit cells and arbitrary basis functions without introducing uncontrolled errors or prohibitive extra cost.

What would settle it

Numerical calculations of the gallium arsenide fundamental gap that show no meaningful acceleration in Ewald-sum convergence when the model density is added would falsify the central claim.

read the original abstract

The evaluation of the electrostatic potential is fundamental to the study of condensed phase systems. We discuss the calculation of the relevant lattice summations by Ewald-type techniques. A model charge density is introduced, that cancels multipole moments of the crystalline charge distribution up to a desired order, for accelerating convergence of the Ewald sums. The method is applicable to calculations of bulk systems, employing arbitrary unit cells in a classical or quantum context, and with arbitrary basis functions to represent the charge density. The efficacy of the method is demonstrated on the calculation of the fundamental gap of the gallium arsenide bulk semiconductor, as a prototype example, where significantly accelerated convergence is numerically confirmed. The approach clarifies a decades-old implementation in the CRYSTAL code.

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 proposes introducing a model charge density that cancels multipole moments of the crystalline charge distribution up to a chosen order, thereby accelerating convergence of Ewald-type lattice summations for the electrostatic potential. The approach is presented as general for bulk periodic systems with arbitrary unit cells, applicable in both classical and quantum contexts, and compatible with arbitrary basis functions for representing the charge density. Efficacy is demonstrated via numerical results on the fundamental gap of bulk GaAs, where accelerated convergence is reported, and the method is said to clarify a decades-old implementation in the CRYSTAL code.

Significance. If the central construction can be shown to achieve exact multipole cancellation without residual errors, uncontrolled approximations, or prohibitive extra cost across arbitrary bases and cells, the technique would provide a practical enhancement to standard Ewald methods for condensed-matter electrostatics. The GaAs demonstration supplies initial numerical support, and clarifying the CRYSTAL implementation adds archival value to existing software practices.

major comments (2)
  1. [§2 (Model density construction)] The construction of the model charge density for arbitrary basis functions (localized, plane-wave, or mixed) must be specified explicitly, including any linear algebra steps or projections, to confirm that multipole cancellation remains exact and does not introduce basis-dependent truncation or conditioning issues that could affect the reported GaAs gap convergence at the level of the claimed acceleration.
  2. [§4 (Numerical results)] The numerical confirmation for the GaAs fundamental gap requires a direct comparison table or plot against a fully converged reference Ewald calculation (with error metrics such as ΔE in meV) and against standard Ewald at equivalent computational effort; without these, it is impossible to verify that the acceleration is free of hidden tuning or post-hoc adjustments.
minor comments (2)
  1. The abstract states 'significantly accelerated convergence' but does not specify the multipole order chosen for the GaAs example or the observed speedup factor; adding these quantitative details would improve clarity.
  2. A brief statement on the scaling of the extra work required to build the model density relative to the Ewald sum itself would help readers assess practicality for large cells.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their thorough review and valuable suggestions. We address the major comments point by point below and will incorporate the necessary clarifications and additional data in the revised manuscript.

read point-by-point responses

  1. Referee: [§2 (Model density construction)] The construction of the model charge density for arbitrary basis functions (localized, plane-wave, or mixed) must be specified explicitly, including any linear algebra steps or projections, to confirm that multipole cancellation remains exact and does not introduce basis-dependent truncation or conditioning issues that could affect the reported GaAs gap convergence at the level of the claimed acceleration.

    Authors: We agree that more explicit details are needed for generality. In the revised version, we will provide the explicit construction in §2: the model density is a linear combination of the same basis functions, with coefficients determined by solving a small linear system that enforces exact matching of the multipole moments (up to order N) of the total charge density. The multipole moments are computed via analytic integrals over the basis, ensuring exact cancellation independent of the basis type. No projections or truncations are involved beyond the chosen multipole order. For the GaAs case, we will add a note confirming the matrix is well-conditioned for the employed basis. revision: yes

  2. Referee: [§4 (Numerical results)] The numerical confirmation for the GaAs fundamental gap requires a direct comparison table or plot against a fully converged reference Ewald calculation (with error metrics such as ΔE in meV) and against standard Ewald at equivalent computational effort; without these, it is impossible to verify that the acceleration is free of hidden tuning or post-hoc adjustments.

    Authors: We acknowledge this point and will revise §4 to include the requested comparisons. A new table will present the GaAs gap as a function of the Ewald convergence parameter, with and without the model density, alongside a reference value from a highly converged standard Ewald sum (error < 0.1 meV). We will also report the number of lattice vectors required for convergence to a given tolerance, demonstrating the reduction in computational effort. No post-hoc adjustments were used; the parameters are as described in the text. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected; method extends standard Ewald without self-referential reduction

full rationale

The derivation introduces a compensating model density to cancel multipoles up to chosen order, then applies standard Ewald summation to the corrected neutral system. No equations or steps in the abstract or description reduce the claimed acceleration to a fitted parameter, self-citation load-bearing premise, or ansatz smuggled from prior work by the same authors. The GaAs gap demonstration is presented as numerical confirmation rather than a prediction forced by construction. The approach is framed as clarifying an existing CRYSTAL implementation without the central result depending on self-referential definitions or uniqueness theorems imported from the authors' own prior papers.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 1 invented entities

The central claim rests on the existence of a constructible model density that can match the real charge distribution's multipoles up to arbitrary order for any periodic cell and basis; this is an unproven domain assumption whose generality is asserted but not derived from first principles in the visible text.

axioms (2)
  • standard math Ewald-type lattice summation converges for neutral periodic charge distributions
    Standard background result invoked for all Ewald techniques
  • domain assumption A model density can be chosen to cancel multipole moments up to any finite order without altering the physical potential
    Core unproven premise of the acceleration method
invented entities (1)
  • model charge density no independent evidence
    purpose: Cancels multipole moments of the crystalline charge to accelerate Ewald convergence
    Newly postulated auxiliary density introduced by the paper

pith-pipeline@v0.9.0 · 5430 in / 1361 out tokens · 27340 ms · 2026-05-16T09:42:58.302567+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.

  • IndisputableMonolith/Foundation/AlexanderDuality.lean alexander_duality_circle_linking unclear
    ?
    unclear

    Relation between the paper passage and the cited Recognition theorem.

    A model charge density is introduced, that cancels multipole moments of the crystalline charge distribution up to a desired order, for accelerating convergence of the Ewald sums... η^m_ℓ[¯n](R) = η^m_ℓ[n](R) ... κ^m_ℓ(r-R) = R^m_ℓ(∥r-R∥) X^m_ℓ(r-R) with R^m_ℓ(r) = C^m_ℓ δ(r) r^{-2(ℓ+1)}

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

43 extracted references · 43 canonical work pages

  1. [1]

    H. M. Evjen, On the Stability of Certain Heteropolar Crystals, Phys. Rev.39, 675-687 (1932)

    work page 1932

  2. [2]

    L. Z. Stolarczyk and L. Piela, Direct calculation of lattice sums. A method to account for the crystal field effects, Int. J. Quantum Chem.22, 911-927 (1982)

    work page 1982

  3. [3]

    K. N. Kudin and G. E. Scuseria, A fast multipole method for pe- riodic systems with arbitrary unit cell geometries, Chem. Phys. Lett.283, 61-68 (1998)

    work page 1998

  4. [4]

    Kadek, M

    M. Kadek, M. Repisky, and K. Ruud, All-electron fully relativis- tic Kohn-Sham theory for solids based on the Dirac-Coulomb Hamiltonian and Gaussian-type functions, Phys. Rev. B99, 205103 (2019)

    work page 2019

  5. [5]

    D. M. Heyes and K. D. Hammonds, Calculating Coulomb inter- actions in molecular dynamics simulations: The Evjen method revisited, J. Chem. Phys.164, 014505 (2026)

    work page 2026

  6. [6]

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

    work page 1957

  7. [7]

    Challacombe, C

    M. Challacombe, C. White and M. Head-Gordon, Periodic boundary conditions and the fast multipole method, J. Chem. Phys.107, 10131-10140 (1997)

    work page 1997

  8. [8]

    Hautman, J

    J. Hautman, J. W. Halley, and Y .-J. Rhee, Molecular dynamics simulation of water between two ideal classical metal walls, J. Chem. Phys.91, 467-472 (1989)

    work page 1989

  9. [9]

    A. Y . Toukmaji and J. A. Board, Ewald summation techniques in perspective: a survey, Comput. Phys. Commun.95, 73-92 (1996)

    work page 1996

  10. [10]

    Hummer, N

    G. Hummer, N. Grønbech-Jensen, and M. Neumann, Pressure calculation in polar and charged systems using Ewald summation: Results for the extended simple point charge model of water, J. Chem. Phys.109, 2791-2797 (1998)

    work page 1998

  11. [11]

    Kawata, M

    M. Kawata, M. Mikami, Rapid calculation of two-dimensional Ewald summation, Chem. Phys. Lett.340, 157-164 (2001)

    work page 2001

  12. [12]

    D. J. Price and C. L. Brooks III, A modified TIP3P water poten- tial for simulation with Ewald summation, J. Chem. Phys.121, 10096-10103 (2004)

    work page 2004

  13. [13]

    Krieger, J

    E. Krieger, J. E. Nielsen, C. A. Spronk, G. Vriend, Fast empirical pKa prediction by Ewald summation, J. Mol. Graph. Model.25, 481-486 (2006)

    work page 2006

  14. [14]

    T. R. Gingrich, and M. Wilson, On the Ewald summation of Gaussian charges for the simulation of metallic surfaces, Chem. Phys. Lett.500, 178-183 (2010)

    work page 2010

  15. [15]

    B. A. Wells and A. L. Chaffee, Ewald summation for molecular simulations, J. Chem. Theory Comput.11, 3684-3695 (2015)

    work page 2015

  16. [16]

    Shukla, M

    A. Shukla, M. Dolg, H. Stoll, and P. Fulde, An ab initio embedded-cluster approach to electronic structure calculations on perfect solids: a Hartree-Fock study of lithium hydride, Chem. Phys. Lett.262, 213-218 (1996)

    work page 1996

  17. [17]

    Kwangho, G

    N. Kwangho, G. Jiali, and D. M. York, An efficient linear-scaling Ewald method for long-range electrostatic interactions in com- bined QM/MM calculations, J. Chem. Theor. Comput.1, 2-13 (2005). Page 7 of 8

    work page 2005

  18. [18]

    N. D. Hine, M. Robinson, P. D. Haynes, C. K. Skylaris, M. C. Payne, and A. A. Mostofi, Accurate ionic forces and geometry optimization in linear-scaling density-functional theory with local orbitals, Phys. Rev. B83, 195102 (2011)

    work page 2011

  19. [19]

    Zhang, D

    P. Zhang, D. G. Truhlar, and J. Gao, Fragment-based quantum mechanical methods for periodic systems with Ewald summation and mean image charge convention for long-range electrostatic interactions, Phys. Chem. Chem. Phys.14, 7821-7829 (2012)

    work page 2012

  20. [20]

    R. Zhao, Y . Zhang, Y . Xiao, and W. Liu, Exact two-component relativistic energy band theory and application, J. Chem. Phys. 144, 044105 (2016)

    work page 2016

  21. [21]

    C. H. Patterson, Density fitting in periodic systems: application to TDHF in diamond and oxides, J. Chem. Phys.153, 064107 (2020)

    work page 2020

  22. [22]

    J. P. Pederson and J. G. McDaniel, DFT-based QM/MM with particle-mesh Ewald for direct, long-range electrostatic embed- ding, J. Chem. Phys.156, 174105 (2022)

    work page 2022

  23. [23]

    M. A. G. Bl ´azquez and J. J. Palacios, First-principles excitons in periodic systems with Gaussian density fitting and Ewald poten- tial functions, Phys. Rev. Res.7, 013156 (2025)

    work page 2025

  24. [24]

    Alrakik, G

    A. Alrakik, G. L. Bendazzoli, S. Evangelisti, and J. A. Berger, Quantum Chemistry for Solids Made Simple on the Clifford Torus, Phys. Rev. Lett.136, 016402 (2026)

    work page 2026

  25. [25]

    Pisani, R

    C. Pisani, R. Dovesi, and C. Roetti, Hartree-Fock Ab Initio Treat- ment of Crystalline Systems, Springer Berlin, Heidelberg (1988)

    work page 1988

  26. [26]

    V . R. Saunders, C. Freyria-Fava, R. Dovesi, L. Salasco, and C. Roetti, On the electrostatic potential in crystalline systems where the charge density is expanded in Gaussian functions, Mol. Phys. 77, 629-665 (1992)

    work page 1992

  27. [27]

    Knopp, Infinite Sequences and Series, 1st Edition, Dover Publications, New York (1956)

    See Chapter 3, Section 3.3, Theorem 3 in the book by K. Knopp, Infinite Sequences and Series, 1st Edition, Dover Publications, New York (1956)

    work page 1956

  28. [28]

    See Chapter IV , Section 4.43, paragraph II in the book by E. T. Whittaker and G. N. Watson, A Course of Modern Analysis, 4th Edition, Cambridge University Press, Cambridge Mathematical Library (1927)

    work page 1927

  29. [29]

    P. P. Ewald, Die Berechnung optischer und elektrostatischer Git- terpotentiale, Ann. Phys. (Leipzig)369, 253-287 (1921)

    work page 1921

  30. [30]

    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. Lond. A373, 27-56 (1980)

    work page 1980

  31. [31]

    Makov and M

    G. Makov and M. C. Payne, Periodic boundary conditions in ab initio calculations, Phys. Rev. B51, 4014-4022 (1995)

    work page 1995

  32. [32]

    J. K. Desmarais, A. Erba, and J. P. Flament, Structural relaxation of materials with spin-orbit coupling: Analytical forces in spin- current DFT, Phys. Rev. B108, 134108 (2023)

    work page 2023

  33. [33]

    Maschio, B

    L. Maschio, B. Kirtman, M. R´erat, R. Orlando, and R. Dovesi, Ab initio analytical Raman intensities for periodic systems through a coupled perturbed Hartree-Fock/Kohn-Sham method in an atomic orbital basis. I. Theory, J. Chem. Phys.139, 164101 (2013)

    work page 2013

  34. [34]

    Maschio, B

    L. Maschio, B. Kirtman, R. Orlando, and M. R´erat, Ab initio ana- lytical infrared intensities for periodic systems through a coupled perturbed Hartree-Fock/Kohn-Sham method, J. Chem. Phys.137, 204113 (2012)

    work page 2012

  35. [35]

    K. Doll, V . R. Saunders, and N. M. Harrison, Analytical Hartree–Fock gradients for periodic systems, Int. J. Quantum Chem.82, 1-13 (2001)

    work page 2001

  36. [36]

    W. W. Bell, Special Functions for Scientists and Engineers, Butler and Tanner Ltd., Frome and London (1968)

    work page 1968

  37. [37]

    Barton, Elements of Green’s Functions and Propagation: Po- tentials, Diffusion, and Waves, Clarendon Press (1989)

    G. Barton, Elements of Green’s Functions and Propagation: Po- tentials, Diffusion, and Waves, Clarendon Press (1989)

    work page 1989

  38. [38]

    Arfken, Mathematical Methods for Physicists, 3rd edition, Academic Press, Inc., San Diego, California, (1985)

    G. Arfken, Mathematical Methods for Physicists, 3rd edition, Academic Press, Inc., San Diego, California, (1985)

    work page 1985

  39. [39]

    Ribaldone and J

    C. Ribaldone and J. K. Desmarais, Spherical to Cartesian coordi- nates transformation for solid harmonics revisited: Construction of the Hartree potential, J. Chem. Phys.163, 074102 (2025)

    work page 2025

  40. [40]

    Abramowitz and I

    M. Abramowitz and I. A. Stegun, Handbook of Mathemati- cal Functions with Formulas, Graphs, and Mathematical Tables, Dover, New York (1964)

    work page 1964

  41. [41]

    Hobson, The Theory of Spherical and Ellipsoidal Harmonics, Cambridge University Press (1931)

    E. Hobson, The Theory of Spherical and Ellipsoidal Harmonics, Cambridge University Press (1931)

    work page 1931

  42. [42]

    J. M. Ziman, Principles of the Theory of Solids, 2nd Edition, Cambridge University Press (1972)

    work page 1972

  43. [43]

    Gradshteyn and I

    I. Gradshteyn and I. Ryzhik, Table of Integrals, Series, and Prod- ucts, 7th Edition, Academic Press (2007). Page 8 of 8

    work page 2007