pith. sign in

arxiv: 2601.09009 · v2 · submitted 2026-01-13 · ❄️ cond-mat.stat-mech · cond-mat.mtrl-sci· physics.comp-ph

Divergent Fluctuations from a 2D Infrared Catastrophe

Richard G. Hennig , Clotilde S. Cucinotta This is my paper

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

classification ❄️ cond-mat.stat-mech cond-mat.mtrl-sciphysics.comp-ph
keywords periodic boundary conditionselectrostatic fluctuationsinterfaces2D periodicityBrownian bridgedipolar mediasimulation artifacts
0
0 comments X

The pith

Lateral periodic boundaries create unscreened modes that make electrostatic potential fluctuations diverge with depth.

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

The paper shows that periodic boundary conditions in the plane of an interface force every lateral replica to carry identical charge fluctuations. This produces an unscreened uniform mode that converts the plane-averaged potential into a stochastic integral along the perpendicular direction. In a semi-infinite slab the variance therefore grows linearly with depth. In a finite cell with fixed-potential boundaries the profile becomes a parabolic Brownian bridge pinned at both ends, with amplitude inversely proportional to the lateral area. These effects disappear as soon as lateral periodicity is removed, so the authors supply an analytic expression for dipolar media that lets researchers choose safe lateral cell sizes.

Core claim

The q_parallel=0 mode arising from lateral periodicity reduces the plane-averaged potential to a stochastic integral of the charge density along z. In a semi-infinite slab this produces linearly growing variance with depth; in a finite periodic cell pinned at fixed potential it produces a Brownian bridge whose variance is inversely proportional to the lateral cell area. These fluctuations are a direct consequence of identical charge fluctuations in every lateral replica and remain bounded only when lateral periodicity is removed.

What carries the argument

The unscreened q_parallel=0 mode from lateral periodic boundary conditions, which converts the plane-averaged potential into a Brownian bridge or linear random walk along z.

If this is right

  • The fluctuation amplitude in the center of a slab scales with the square of the slab thickness divided by the lateral area.
  • Electrostatic properties extracted from periodic interface simulations contain size-dependent artifacts that increase with perpendicular extent.
  • An analytic formula for dipolar media gives the minimum lateral area needed to keep the fluctuation magnitude below a chosen threshold.
  • The divergences vanish entirely in systems that are non-periodic or finite in the lateral directions.

Where Pith is reading between the lines

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

  • The same infrared mode may affect other long-range interactions under 2D periodicity, such as in gravitational or magnetic simulations of slabs.
  • A correction term could be added to Ewald or particle-mesh methods to subtract the contribution of this uniform mode without changing the lateral periodicity.
  • Direct numerical checks of the predicted parabolic profile versus lateral area would confirm the scaling without requiring changes to the simulation code.

Load-bearing premise

Every lateral replica carries exactly identical charge fluctuations at every instant, leaving the uniform mode completely unscreened.

What would settle it

Measure the variance of the plane-averaged potential in the same interfacial system simulated once with lateral periodic boundaries and once without; the variance should grow with depth only in the periodic case.

Figures

Figures reproduced from arXiv: 2601.09009 by Clotilde S. Cucinotta, Richard G. Hennig.

Figure 1
Figure 1. Figure 1: FIG. 1 [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
read the original abstract

Molecular simulations of interfacial polar media routinely employ periodic boundary conditions parallel to the interface. We show that this lateral periodicity introduces a spatially uniform in-plane mode ($q_{\parallel}=0$) that is unscreened because every lateral replica carries identical charge fluctuations. This 2D mode reduces the plane-averaged potential to a stochastic integral of the plane-averaged charge density along $z$, so that in a semi-infinite slab the variance of the potential grows linearly with depth. In a finite or periodic cell along $z$, with boundaries held at fixed potential, it follows a parabolic profile--a Brownian bridge--pinned to zero at both ends, with amplitude inversely proportional to the lateral cell area. These diverging fluctuations are a pure artifact of the imposed 2D lateral periodicity: they remain bounded in systems that are non-periodic or of finite lateral extent. We provide an analytic expression for their magnitude in dipolar media, yielding a practical criterion for the choice of lateral cell dimensions.

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

Summary. The paper claims that lateral periodic boundary conditions in molecular simulations of interfacial polar media introduce an unscreened q_parallel=0 Fourier mode of the electrostatic potential, since identical charge fluctuations are enforced across replicas. This reduces the plane-averaged potential to a stochastic integral of the plane-averaged charge density along z. Consequently, the potential variance grows linearly with depth in a semi-infinite slab and follows a parabolic Brownian-bridge profile (pinned at fixed-potential boundaries) in a finite or periodic cell along z, with amplitude scaling as 1/A (A = lateral cell area). These diverging fluctuations are presented as a pure artifact of 2D periodicity; they remain bounded for non-periodic or finite-lateral-extent systems. An analytic expression for the magnitude in dipolar media is supplied to guide choice of lateral cell dimensions.

Significance. If the central electrostatic reduction holds, the result identifies a previously under-appreciated artifact that can produce unphysical, diverging potential fluctuations in standard PBC simulations of interfaces. The derivation is parameter-free and follows directly from Poisson electrostatics under lateral periodicity, supplying a concrete, falsifiable criterion (1/A scaling and Brownian-bridge shape) for simulation-cell design in polar media. This has immediate practical value for the large community performing interfacial MD.

major comments (2)
  1. [Theory / derivation of the q=0 mode] The reduction of the 3D Poisson equation to the claimed one-dimensional stochastic integral for the q_parallel=0 mode is asserted but the explicit steps (including the handling of the in-plane average and any boundary terms) are not shown; an equation-by-equation derivation in the main text is required to confirm that no screening term survives.
  2. [Setup of periodic boundary conditions] The statement that 'every lateral replica carries identical charge fluctuations' is presented as a direct consequence of periodicity, yet the manuscript must demonstrate that this holds for the fluctuating charge density (not just the average) when the simulation samples independent configurations; otherwise the unscreened-mode argument is incomplete.
minor comments (3)
  1. Notation for the plane-averaged quantities (e.g., how the lateral integral is normalized) should be defined once at first use and used consistently.
  2. The analytic expression for dipolar media should be written explicitly (with all constants) rather than merely referenced, so readers can apply the cell-size criterion without re-deriving it.
  3. A short comparison figure or table showing the predicted 1/A scaling versus a small set of test simulations would make the result more immediately usable.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive assessment and for highlighting the practical relevance of these electrostatic artifacts. We address the two major comments below and will revise the manuscript to incorporate the requested clarifications.

read point-by-point responses

  1. Referee: [Theory / derivation of the q=0 mode] The reduction of the 3D Poisson equation to the claimed one-dimensional stochastic integral for the q_parallel=0 mode is asserted but the explicit steps (including the handling of the in-plane average and any boundary terms) are not shown; an equation-by-equation derivation in the main text is required to confirm that no screening term survives.

    Authors: We agree that an explicit derivation will strengthen the presentation. In the revised manuscript we will insert a dedicated paragraph (or short appendix) that starts from the 3D Poisson equation, performs the lateral Fourier transform, isolates the q_parallel=0 component, applies the in-plane average, and shows step-by-step that the dielectric screening term vanishes identically for this uniform mode while boundary terms are handled by the fixed-potential or periodic conditions along z. This will confirm the reduction to the stochastic integral without any surviving screening. revision: yes

  2. Referee: [Setup of periodic boundary conditions] The statement that 'every lateral replica carries identical charge fluctuations' is presented as a direct consequence of periodicity, yet the manuscript must demonstrate that this holds for the fluctuating charge density (not just the average) when the simulation samples independent configurations; otherwise the unscreened-mode argument is incomplete.

    Authors: We accept the need for greater explicitness. Under lateral periodic boundary conditions the charge density of each sampled configuration is replicated identically in every lateral image by construction of the simulation cell; i.e., ρ(x+nLx,y+mLy,z)=ρ(x,y,z) for all integers n,m. Consequently the fluctuations themselves—not merely their average—are identical across replicas for every independent snapshot. We will add a concise clarifying sentence in the revised text that states this property for the fluctuating density and notes that it applies snapshot-by-snapshot in the MD ensemble. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The derivation proceeds directly from Poisson electrostatics under the imposed lateral periodic boundary conditions. The q_parallel=0 mode being unscreened follows immediately from the periodicity assumption that all lateral replicas carry identical charge fluctuations, reducing the plane-averaged potential to a one-dimensional stochastic integral without any fitted parameters or self-referential steps. The linear variance growth in semi-infinite slabs and the parabolic Brownian-bridge profile in finite cells are standard consequences of this integral, with the 1/A scaling arising from central-limit averaging over the lateral area. No self-citations, ansatzes, or uniqueness theorems are invoked to support the central claim, and the result remains self-contained against the electrostatic setup without reducing predictions to inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The claim rests on the standard treatment of electrostatics in laterally periodic systems; no new entities or fitted constants are introduced.

axioms (2)
  • domain assumption Lateral periodic boundary conditions force identical charge fluctuations in every replica.
    This is the setup that produces the unscreened q_parallel=0 mode.
  • standard math The electrostatic potential satisfies Poisson's equation in the periodic geometry.
    Used to reduce the plane-averaged potential to an integral over charge density.

pith-pipeline@v0.9.0 · 8801 in / 1221 out tokens · 98402 ms · 2026-05-16T14:10:33.280052+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

34 extracted references · 34 canonical work pages

  1. [1]

    A. A. Kornyshev, J. Phys. Chem. B111, 5545 (2007)

    work page 2007

  2. [2]

    M. V. Fedorov and A. A. Kornyshev, Chem. Rev.114, 2978 (2014)

    work page 2014

  3. [3]

    Merlet, B

    C. Merlet, B. Rotenberg, P. A. Madden, P.-L. Taberna, P. Simon, and M. Salanne, Nat. Mater.11, 306 (2012)

    work page 2012

  4. [4]

    D. T. Limmer, C. Merlet, M. Salanne, D. Chandler, P. A. Madden, R. V. Roij, and B. Rotenberg, Phys. Rev. Lett. 111, 106102 (2013)

    work page 2013

  5. [5]

    M. Z. Bazant, B. D. Storey, and A. A. Kornyshev, Phys. Rev. Lett.106, 046102 (2011)

    work page 2011

  6. [6]

    R. A. Marcus, Rev. Mod. Phys.65, 599 (1993)

    work page 1993

  7. [7]

    Blumberger, Chem

    J. Blumberger, Chem. Rev.115, 11191 (2015)

    work page 2015

  8. [8]

    Fragasso, S

    A. Fragasso, S. Schmid, and C. Dekker, ACS nano14, 1338 (2020)

    work page 2020

  9. [9]

    Dekker, Nat

    C. Dekker, Nat. Nanotechnol.2, 209 (2007)

    work page 2007

  10. [10]

    A. A. Gurtovenko and I. Vattulainen, J. Phys. Chem. B 111, 13554 (2007)

    work page 2007

  11. [11]

    Spohr, J

    E. Spohr, J. Chem. Phys.107, 6342 (1997)

    work page 1997

  12. [12]

    Stengel, D

    M. Stengel, D. Vanderbilt, and N. A. Spaldin, Nat. Phys. 5, 304 (2009)

    work page 2009

  13. [13]

    Khatib, A

    R. Khatib, A. Kumar, S. Sanvito, M. Sulpizi, and C. S. Cucinotta, Electrochim. Acta391, 138875 (2021)

    work page 2021

  14. [14]

    M. T. Darby and C. S. Cucinotta, Curr. Opin. Elec- trochem.36, 101118 (2022)

    work page 2022

  15. [15]

    Raffone, R

    F. Raffone, R. Khatib, M. Sulpizi, and C. Cucinotta, Commun. Chem.8, 58 (2025)

    work page 2025

  16. [16]

    Buraschi, A

    M. Buraschi, A. P. Horsfield, and C. S. Cucinotta, J. Phys. Chem. Lett.15, 4872 (2024)

    work page 2024

  17. [17]

    C. S. Ahart, S. K. Chulkov, and C. S. Cucinotta, Journal of Chemical Theory and Computation20, 6772 (2024)

    work page 2024

  18. [18]

    Deissenbeck, D

    F. Deissenbeck, D. Toroz, R. Wippermann, and J. Neuge- bauer, Phys. Rev. Lett.126, 136803 (2021)

    work page 2021

  19. [19]

    Todorova, R

    T. Todorova, R. Wippermann, and J. Neugebauer, arXiv preprint (2025), arXiv:2411.05925 [cond-mat.mtrl-sci]

    work page arXiv 2025

  20. [20]

    Yeh and M

    I.-C. Yeh and M. L. Berkowitz, J. Chem. Phys.111, 3155 (1999)

    work page 1999

  21. [21]

    A. Y. Toukmaji and J. A. B. Jr., Comput. Phys. Commun. 95, 73 (1996)

    work page 1996

  22. [22]

    S. W. de Leeuw, J. W. Perram, and E. R. Smith, Proc. R. Soc. A373, 27 (1980)

    work page 1980

  23. [23]

    S. W. de Leeuw, J. W. Perram, and E. R. Smith, Proc. R. Soc. A373, 57 (1980)

    work page 1980

  24. [24]

    S. Chen, C. S. Niu, and Z. Q. Li, Phys. Rev. B83, 165125 (2011)

    work page 2011

  25. [25]

    Otani and O

    M. Otani and O. Sugino, Phys. Rev. B73, 115407 (2006)

    work page 2006

  26. [26]

    Planck, Ann

    M. Planck, Ann. Phys.4, 553 (1901)

    work page 1901

  27. [27]

    G. G. Stokes, Trans. Cambridge Philos. Soc.9, 8 (1851)

    work page

  28. [28]

    Feller,An Introduction to Probability Theory and Its Applications, 2nd ed., Vol

    W. Feller,An Introduction to Probability Theory and Its Applications, 2nd ed., Vol. II (Wiley, New York, 1971)

    work page 1971

  29. [29]

    J. D. Jackson,Classical Electrodynamics, 3rd ed. (Wiley, New York, 1998)

    work page 1998

  30. [30]

    Hansen and I

    J.-P. Hansen and I. R. McDonald,Theory of Simple Liq- uids, 4th ed. (Academic Press, Oxford, 2013)

    work page 2013

  31. [31]

    Plimpton, J

    S. Plimpton, J. Comput. Phys.117, 1 (1995)

    work page 1995

  32. [32]

    D. J. Price and B. R. Brooks, J. Chem. Phys.121, 10096 (2004)

    work page 2004

  33. [33]

    W. L. Jorgensen, J. Chandrasekhar, J. D. Madura, R. W. Impey, and M. L. Klein, J. Chem. Phys.79, 926 (1983)

    work page 1983

  34. [34]

    Ballenegger, J

    V. Ballenegger, J. Chem. Phys.140, 161102 (2014)

    work page 2014