Approaching the thermodynamic limit of a bounded one-component plasma
Pith reviewed 2026-05-17 23:55 UTC · model grok-4.3
The pith
Bounded one-component plasma simulations extrapolate electrostatic energies to the thermodynamic limit across wide coupling strengths.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
By conducting molecular dynamics simulations on a series of sufficiently large bounded one-component plasmas with spherical reflecting boundaries, the size dependencies of the electrostatic energies are established and extrapolated to the thermodynamic limit of infinite size. The total electrostatic energy per ion is thus estimated for Gamma from 0.03 to 1000 with 0.1% relative error, being 0.5% lower than modern Monte Carlo data at Gamma<30 and coinciding at Gamma>175. Two more converging characteristic energies, the excess interatomic electrostatic energy and the excess ion-background electrostatic energy, are introduced to calculate the ionic compressibility factor. This leads to a wide-r
What carries the argument
Extrapolation of size-dependent electrostatic energies obtained from molecular dynamics of the bounded one-component plasma in spherical reflecting-boundary geometry.
Load-bearing premise
The observed size dependencies of the electrostatic energies in the spherical reflecting-boundary geometry permit a reliable extrapolation to the thermodynamic limit without residual boundary artifacts that survive at infinite radius.
What would settle it
A calculation of the electrostatic energy using an independent method, such as very large-scale periodic boundary condition simulations at a specific Gamma like 10, that disagrees with the extrapolated value by more than the stated 0.1% error would challenge the reliability of the thermodynamic limit estimate.
Figures
read the original abstract
The classical one-component plasma (OCP) bounded by a spherical surface reflecting ions (BOCP) is studied using molecular dynamics (MD). Simulations performed for a series of sufficiently large BOCP's make it possible to establish the size dependencies for the investigated quantities and extrapolate them to the thermodynamic limit. In particular, the total electrostatic energy per ion is estimated in the limit of infinite BOCP in a wide range of the Coulomb coupling parameter $\Gamma$ from 0.03 to 1000 with the relative error of the order 0.1%. Calculated energies are by about 0.5% lower as compared to the modern Monte Carlo (MC) simulation data obtained by different authors at $\Gamma<30$ and almost coincide with the MC results at $\Gamma>175$. We introduce two more converging characteristic energies, the excess interatomic electrostatic energy and the excess ion-background electrostatic energy, which enable us to calculate the ionic compressibility factor inaccessible in conventional MC and MD simulation of the OCP with periodic boundary conditions. The derived wide-range ionic equation of state can be recommended for testing OCP simulations with various effective interaction potentials. Based on this equation, we propose an improved cutoff radius for the interionic forces implemented in LAMMPS and perform MD simulation of the OCP to demonstrate that location of the metastable region of the fluid-solid phase transition depends sensitively on this radius.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript uses molecular dynamics simulations of the bounded one-component plasma (BOCP) confined by spherical reflecting boundaries. Simulations of a sequence of large finite systems are used to determine size dependencies of electrostatic energies, which are then extrapolated to the thermodynamic limit. The central result is the total electrostatic energy per ion for Γ from 0.03 to 1000, reported with ~0.1% relative error; these values lie ~0.5% below modern Monte Carlo data at Γ<30 and agree at Γ>175. Two additional converging energies (excess interatomic and excess ion-background) are introduced to obtain the ionic compressibility factor. An improved LAMMPS cutoff radius is proposed and its effect on the location of the fluid-solid metastable region is demonstrated.
Significance. If the extrapolation is reliable, the work supplies high-precision thermodynamic-limit reference energies for the OCP over a wide Γ range that can benchmark other simulation techniques. The excess-energy definitions provide a route to the compressibility factor that is unavailable in conventional periodic-boundary simulations. Direct MD generation of the data and independent comparison to external MC results are strengths. The practical LAMMPS cutoff recommendation adds immediate utility for plasma simulations.
major comments (1)
- [Abstract and finite-size extrapolation procedure] The abstract and the section describing the finite-size extrapolation state that energies were obtained by extrapolation from a series of large systems with 0.1% relative error, yet no functional form for the size dependence (e.g., 1/R, 1/N^{1/3} plus higher-order terms), number of system sizes, fitting procedure, or statistical uncertainties on the fit coefficients are supplied. Because the claimed 0.5% offset relative to periodic MC data at Γ<30 rests entirely on the accuracy of this extrapolation, residual surface or curvature corrections specific to the reflecting spherical geometry could survive in the N→∞ limit and contribute to the reported difference.
minor comments (1)
- [Notation and equations] Notation for the various electrostatic energies (total, excess interatomic, excess ion-background) should be introduced once and used consistently in all equations and figures.
Simulated Author's Rebuttal
We thank the referee for the thorough review and constructive feedback on our manuscript. We address the major comment below and will revise the manuscript accordingly to strengthen the presentation of our finite-size extrapolation procedure.
read point-by-point responses
-
Referee: The abstract and the section describing the finite-size extrapolation state that energies were obtained by extrapolation from a series of large systems with 0.1% relative error, yet no functional form for the size dependence (e.g., 1/R, 1/N^{1/3} plus higher-order terms), number of system sizes, fitting procedure, or statistical uncertainties on the fit coefficients are supplied. Because the claimed 0.5% offset relative to periodic MC data at Γ<30 rests entirely on the accuracy of this extrapolation, residual surface or curvature corrections specific to the reflecting spherical geometry could survive in the N→∞ limit and contribute to the reported difference.
Authors: We agree that additional details on the extrapolation are necessary for full reproducibility and to address concerns about residual corrections. In the revised manuscript, we will explicitly state the functional form employed: a leading 1/R (or equivalently 1/N^{1/3}) term motivated by the surface-charge correction for a spherical reflecting boundary, supplemented by a constant term and, where statistically justified, a 1/N term. We simulated 6–8 system sizes per Γ value, with N ranging from approximately 5×10^4 to 2×10^6 ions. The fit is performed via weighted least-squares minimization, and we will report the resulting extrapolated values together with their standard errors obtained from the fit covariance matrix (typically yielding relative uncertainties of 0.05–0.15%). Regarding possible surviving surface or curvature corrections, the reflecting spherical geometry is chosen precisely because it permits a clean extrapolation in which leading finite-size effects vanish as 1/R; our data collapse across the simulated range and the smooth Γ dependence of the extrapolated energies indicate that higher-order corrections are negligible within the quoted 0.1% precision. The systematic 0.5% offset relative to periodic-boundary MC results at low Γ is therefore interpreted as a genuine physical distinction between the two ensembles rather than an artifact of incomplete extrapolation. revision: yes
Circularity Check
Finite-size MD extrapolation to thermodynamic limit is data-driven and independent of prior fits
full rationale
The paper performs direct molecular dynamics simulations on finite BOCP systems with spherical reflecting boundaries, measures size dependence of electrostatic energies across multiple N, and extrapolates the observed trends to N→∞. This numerical procedure does not reduce any claimed result to a fitted parameter by algebraic construction, nor does it invoke self-citations for uniqueness theorems or smuggle ansatzes. The reported energies (with ~0.1% relative error) and the 0.5% offset versus external MC data at low Γ are outputs of the simulation-plus-extrapolation pipeline rather than tautological redefinitions. Introduction of excess interatomic and ion-background energies is a post-processing step that enables compressibility calculations unavailable in periodic-boundary OCP, but does not create circularity in the primary energy estimates. The overall chain remains self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
free parameters (2)
- extrapolation function coefficients
- proposed LAMMPS cutoff radius
axioms (2)
- domain assumption Classical molecular dynamics with reflecting boundaries faithfully samples the canonical ensemble for the BOCP.
- domain assumption Finite-size corrections vanish as a smooth, extrapolatable function of radius.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
u(Γ,N)=u∞(Γ)+c1(Γ)/lnN+c2(Γ) (Eq. 19); Zi=1+Γ(upex−2ubex)/3 (Eq. 16); virial theorem application (Eqs. 15,20)
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
extrapolation of electrostatic energies to N→∞ thermodynamic limit; comparison with periodic MC data
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
-
[1]
are noticeably higher. AtΓ>174, our data reach asymptotically the bcc energyu0, and in this region, they are in a very good agreement with [43] (Fig. 5). At the same time, for0.1<Γ<174, our energies are appre- ciably lower than other data, which can be clearly seen in the inset in Fig. 3. Here, a typical difference is about 0.5%, which is much higher than...
-
[2]
and our data. Based on Eqs. (38) and 39), we can approximate the internal energy in a wide range ofΓin the form εiex(Γ) = ( ufit(Γ)−u 0,0.1≤Γ≤Γ m, 3 2Γ + β2 Γ2 ,Γ>Γ m, (40) 9 Figure 2. Size dependence of the BOCP total electrostatic energy from MD simulation (dots) and its fit by Eq. (19) (lines): (a) diamonds and solid lines correspond toΓ = 0.1 and circ...
-
[3]
[69], compiled with the ScaFaCoS library v1.0.4 [70]. The total number of particlesNwas fixed within each simulation, but varied between different simulations and the exact values will be specified later. Periodic bound- ary conditions were applied in all three directions of a cubic cell with the edge lengthL= (4πN/3) 1/3. The integration time step was ch...
-
[4]
Baus, Statistical mechanics of simple Coulomb sys- tems, Physics Reports59, 1 (1980)
M. Baus, Statistical mechanics of simple Coulomb sys- tems, Physics Reports59, 1 (1980)
work page 1980
-
[5]
S.Ichimaru,Stronglycoupledplasmas: high-densityclas- sical plasmas and degenerate electron liquids, Reviews of Modern Physics54, 1017 (1982)
work page 1982
-
[6]
E. V. Kholopov, Convergence problems of Coulomb and multipole sums in crystals, Physics-Uspekhi47, 965 (2004)
work page 2004
-
[7]
E. H. Lieb,The Stability of Matter: From Atoms to Stars: Selecta of Elliott H. Lieb, edited by W. Thirring (Springer Berlin Heidelberg, 2005)
work page 2005
-
[8]
S. Ichimaru, H. Iyetomi, and S. Ogata, Freezing tran- sition and phase diagram of dense carbon-oxygen mix- tures in white dwarfs, The Astrophysical Journal334, L17 (1988)
work page 1988
-
[9]
G. Chabrier, N. W. Ashcroft, and H. E. DeWitt, White dwarfs as quantum crystals, Nature360, 48 (1992)
work page 1992
-
[10]
G. Chabrier, Quantum effects in dense coulumbic matter - application to the cooling of white dwarfs, The Astro- physical Journal414, 695 (1993)
work page 1993
-
[11]
A. Kozhberov, Thermal evolution of old white dwarfs, Journal of Physics: Conference Series929, 012012 (2017)
work page 2017
-
[12]
S. Blouin and J. Daligault, Phase separation in ultra- massive white dwarfs, The Astrophysical Journal919, 87 (2021)
work page 2021
-
[13]
D. A. Baiko, Phase diagrams of binary ionic mixtures and white dwarf cooling, Monthly Notices of the Royal Astronomical Society517, 3962 (2022)
work page 2022
-
[14]
Z. Medin and A. Cumming, Crystallization of classical multicomponent plasmas, Physical Review E81, 036107 (2010)
work page 2010
-
[15]
M. Iwamatsu, Calculation of the free energy of liq- uid alkali metal using the non-local pseudopotential — The one-component-plasma and the hard-sphere refer- ence systems, Physica B+C138, 310 (1986)
work page 1986
-
[16]
D. J. Wineland, J. C. Bergquist, W. M. Itano, J. J. Bollinger, and C. H. Manney, Atomic-ion Coulomb clus- ters in an ion trap, Physical Review Letters59, 2935 (1987)
work page 1987
-
[17]
D. H. E. Dubin and T. M. O’Neil, Trapped nonneutral plasmas, liquids, and crystals (the thermal equilibrium states), Reviews of Modern Physics71, 87 (1999)
work page 1999
-
[18]
O. Arp, D. Block, A. Piel, and A. Melzer, Dust Coulomb balls: Three-dimensional plasma crystals, Physical Re- view Letters93, 165004 (2004)
work page 2004
-
[19]
O. Arp, D. Block, M. Klindworth, and A. Piel, Confine- ment of Coulomb balls, Physics of Plasmas12, 122102 (2005)
work page 2005
-
[20]
S. Käding and A. Melzer, Three-dimensional stereoscopy of Yukawa (Coulomb) balls in dusty plasmas, Physics of Plasmas13, 090701 (2006)
work page 2006
-
[21]
H. Totsuji, T. Ogawa, C. Totsuji, and K. Tsuruta, Struc- ture of spherical Yukawa clusters, Journal of Physics A: Mathematical and General39, 4545 (2006)
work page 2006
-
[22]
S. W. S. Apolinario, J. A. Aguiar, and F. M. Peeters, An- gular melting scenarios in binary dusty-plasma Coulomb balls: Magic versus normal clusters, Physical Review E 16 90, 063113 (2014)
work page 2014
-
[23]
D. I. Zhukhovitskii, V. N. Naumkin, A. I. Khusnulgatin, V. I. Molotkov, and A. M. Lipaev, Dust coupling param- eter of radio-frequency-discharge complex plasma under microgravity conditions, Physical Review E96, 043204 (2017)
work page 2017
- [24]
-
[25]
M. A. Pokrant, Thermodynamic properties of the nonzero-temperature, quantum-mechanical, one- component plasma, Physical Review A16, 413 (1977)
work page 1977
-
[26]
D. Ceperley, Ground state of the fermion one-component plasma: A Monte Carlo study in two and three dimen- sions, Physical Review B18, 3126 (1978)
work page 1978
-
[27]
D. A. Baiko, D. G. Yakovlev, H. E. De Witt, and W. L. Slattery, Coulomb crystals in the harmonic lattice ap- proximation, Physical Review E61, 1912 (2000)
work page 1912
-
[28]
A. Chugunov, D. Baiko, D. Yakovlev, H. De Witt, and W. Slattery, Pair distribution of ions in Coulomb crys- tals, Physica A: Statistical Mechanics and its Applica- tions323, 413 (2003)
work page 2003
-
[29]
A. Chugunov and D. Baiko, Anharmonic corrections to the electrostatic energy of a Coulomb crystal, Physica A: Statistical Mechanics and its Applications352, 397 (2005)
work page 2005
-
[30]
R. T. Farouki and S. Hamaguchi, Thermodynamics of strongly-coupled Yukawa systems near the one- component-plasma limit. II. molecular dynamics simula- tions, The Journal of Chemical Physics101, 9885 (1994)
work page 1994
-
[31]
S. Hamaguchi and R. T. Farouki, Thermodynamics of strongly-coupled Yukawa systems near the one- component-plasma limit. I. derivation of the excess en- ergy, The Journal of Chemical Physics101, 9876 (1994)
work page 1994
-
[32]
S. Hamaguchi, R. T. Farouki, and D. H. E. Dubin, Phase diagram of Yukawa systems near the one-component- plasma limit revisited, The Journal of Chemical Physics 105, 7641 (1996)
work page 1996
-
[33]
J. M. Caillol and D. Gilles, Numerical simulations of screened Coulomb systems. A comparison between hy- perspherical and periodic boundary conditions, Journal of Statistical Physics100, 905 (2000)
work page 2000
-
[34]
J. M. Caillol and D. Gilles, Monte Carlo simulations of theYukawaone-componentplasma,JournalofStatistical Physics100, 933 (2000)
work page 2000
-
[35]
S. A. Khrapak and A. G. Khrapak, Simple thermody- namics of strongly coupled one-component-plasma in two and three dimensions, Physics of Plasmas21, 104505 (2014)
work page 2014
-
[36]
D. A. Young, E. M. Corey, and H. E. DeWitt, Analytic fit to the one-component-plasma structure factor, Physical Review A44, 6508 (1991)
work page 1991
-
[37]
N. V. Brilliantov, Accurate first-principle equation of state for the one-component plasma, Contributions to Plasma Physics38, 489 (1998)
work page 1998
-
[38]
G. Chabrier and A. Y. Potekhin, Equation of state of fully ionized electron-ion plasmas, Physical Review E58, 4941 (1998)
work page 1998
-
[39]
J. M. Caillol, Numerical simulations of Coulomb sys- tems: A comparison between hyperspherical and periodic boundary conditions, The Journal of Chemical Physics 111, 6528 (1999)
work page 1999
-
[40]
J. M. Caillol, A new (and better) lower bound for the excess internal energy of the one-component plasma, The Journal of Chemical Physics111, 9695 (1999)
work page 1999
-
[41]
N. Brilliantov, V. Malinin, and R. Netz, Systematic field- theory for the hard-core one-component plasma, The Eu- ropean Physical Journal D - Atomic, Molecular and Op- tical Physics18, 339 (2002)
work page 2002
-
[42]
S. A. Khrapak and A. G. Khrapak, Internal energy of the classical two- and three-dimensional one-component- plasma, Contributions to Plasma Physics56, 270 (2016)
work page 2016
-
[43]
S. G. Brush, H. L. Sahlin, and E. Teller, Monte Carlo study of a one-component plasma. I, The Journal of Chemical Physics45, 2102 (1966)
work page 1966
-
[44]
J. P. Hansen, Statistical mechanics of dense ionized matter. I. Equilibrium properties of the classical one- component plasma, Physical Review A8, 3096 (1973)
work page 1973
-
[45]
W. L. Slattery, G. D. Doolen, and H. E. DeWitt, Im- proved equation of state for the classical one-component plasma, Physical Review A21, 2087 (1980)
work page 2087
-
[46]
W. L. Slattery, G. D. Doolen, and H. E. DeWitt,Nde- pendence in the classical one-component plasma monte carlo calculations, Physical Review A26, 2255 (1982)
work page 1982
-
[47]
S. Ogata and S. Ichimaru, Critical examination ofnde- pendence in the Monte Carlo calculations for a classi- cal one-component plasma, Physical Review A36, 5451 (1987)
work page 1987
-
[48]
J. M. Caillol, Thermodynamic limit of the excess internal energy of the fluid phase of a one-component plasma: A Monte Carlo study, The Journal of Chemical Physics 111, 6538 (1999)
work page 1999
-
[49]
J.-M. Caillol and D. Gilles, An accurate equation of state fortheone-componentplasmainthelowcouplingregime, Journal of Physics A: Mathematical and Theoretical43, 105501 (2010)
work page 2010
-
[50]
G. S. Demyanov and P. R. Levashov, One-component plasma of a million particles via angular-averaged Ewald potential: A Monte Carlo study, Physical Review E106, 015204 (2022)
work page 2022
-
[51]
R. T. Farouki and S. Hamaguchi, Thermal energy of the crystalline one-component plasma from dynamical simu- lations, Physical Review E47, 4330 (1993)
work page 1993
-
[52]
E. L. Pollock and J. P. Hansen, Statistical mechan- ics of dense ionized matter. II. Equilibrium properties and melting transition of the crystallized one-component plasma, Physical Review A8, 3110 (1973)
work page 1973
-
[53]
J. D. Weeks, Volume change on melting for systems with inverse-power-law interactions, Physical Review B24, 1530 (1981)
work page 1981
-
[54]
S. Nordholm, Simple analysis of the thermodynamic properties of the one-component plasma, Chemical Physics Letters105, 302 (1984)
work page 1984
-
[55]
J. W. Lee and F. H. Ree, Perturbation theory of a classi- cal one-component plasma, Physical Review A38, 5714 (1988)
work page 1988
-
[56]
C. N. Likos and N. W. Ashcroft, Self-consistent theory of freezing of the classical one-component plasma, Physical Review Letters69, 316 (1992)
work page 1992
-
[57]
N. Yu, D. Huang, S. Lu, S. Khrapak, and Y. Feng, Uni- versal scaling of transverse sound speed and its isomor- phic property in Yukawa fluids, Physical Review E109, 035202 (2024)
work page 2024
-
[58]
S. A. Khrapak, Y.-F. Wu, and C.-R. Du, A binary col- lision approach in one-component plasma: How close to the Frenkel line can we go?, Physics of Plasmas32, 17 092114 (2025)
work page 2025
-
[59]
P. P. Ewald, Die Berechnung optischer und elektro- statischer Gitterpotentiale, Annalen der Physik369, 253 (1921)
work page 1921
-
[60]
B. Nijboer and F. De Wette, On the calculation of lattice sums, Physica23, 309 (1957)
work page 1957
-
[61]
N. Karasawa and W. A. Goddard, Acceleration of conver- genceforlatticesums,TheJournalofPhysicalChemistry 93, 7320 (1989)
work page 1989
-
[62]
E. Yakub and C. Ronchi, An efficient method for compu- tation of long-ranged Coulomb forces in computer sim- ulation of ionic fluids, The Journal of Chemical Physics 119, 11556 (2003)
work page 2003
-
[63]
G. S. Demyanov, A. S. Onegin, and P. R. Levashov, N- convergence in one-component plasma: Comparison of Coulomb, Ewald, and angular-averaged Ewald poten- tials, Contributions to Plasma Physics64, e202300164 (2024). [61]The Stability of Matter: From Atoms to Stars: Selecta of Elliott H. Lieb(Springer Berlin Heidelberg, 2001)
work page 2024
-
[64]
D. I. Zhukhovitskii and E. E. Perevoshchikov, Structural transition in strongly coupled Coulomb clusters, High Temperature62, 421 (2024)
work page 2024
-
[65]
A. S. Onegin, G. S. Demyanov, and P. R. Levashov, Pres- sure of Coulomb systems with volume-dependent long- rangepotentials,JournalofPhysicsA:Mathematicaland Theoretical57, 205002 (2024)
work page 2024
-
[66]
E. S. Shpil’ko and D. I. Zhukhovitskii, Relevance of the Wigner–Seitz cell approximation for the Coulomb clus- ters, Plasma Physics Reports49, 1207 (2023)
work page 2023
-
[67]
L. D. Landau and E. M. Lifshitz,Statistical Physics, Vol. 5 (Elsevier, Oxford, 1980)
work page 1980
-
[68]
N. Itoh and S. Ichimaru, Harmonic-lattice model for the internal energy of the classical one-component plasma fluid near the crystallization point, Physical Review A 22, 1318 (1980)
work page 1980
-
[69]
R. Hockney and J. Eastwood,Computer Simulation Us- ing Particles(CRC Press, 2021)
work page 2021
-
[70]
D. I. Zhukhovitskii, Multiscale approach to the theory of nonisothermal homogeneous nucleation, The Journal of Chemical Physics160, 194505 (2024)
work page 2024
-
[71]
A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolin- tineanu, W. M. Brown, P. S. Crozier, P. J. in ’t Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton, LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales, Computer Physics Comm...
work page 2022
- [72]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.