pith. sign in

arxiv: 2511.04516 · v2 · submitted 2025-11-06 · ⚛️ physics.plasm-ph

Approaching the thermodynamic limit of a bounded one-component plasma

Pith reviewed 2026-05-17 23:55 UTC · model grok-4.3

classification ⚛️ physics.plasm-ph
keywords one-component plasmathermodynamic limitmolecular dynamicselectrostatic energyCoulomb couplingequation of statecompressibility factorphase transition
0
0 comments X

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.

The authors perform molecular dynamics simulations of the one-component plasma inside a spherical reflecting boundary for various system sizes. They analyze how the electrostatic energies depend on size and extrapolate those dependencies to the infinite size limit representing the thermodynamic limit. This method provides estimates of the total electrostatic energy per ion for coupling parameters Gamma from 0.03 to 1000 with roughly 0.1 percent relative error. The energies are about 0.5 percent lower than previous Monte Carlo data at lower Gamma values but agree closely at higher values. Introducing excess interatomic and ion-background energies allows calculation of the ionic compressibility factor and a wide-range equation of state useful for validating other plasma simulations and refining molecular dynamics parameters.

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

Figures reproduced from arXiv: 2511.04516 by D. I. Zhukhovitskii, E. E. Perevoshchikov (Joint Institute of High Temperatures, Russian Academy of Sciences).

Figure 1
Figure 1. Figure 1: Ion density distribution fc(r) (8) and the radial distribution function fp(r) (7) of BOCP at different Γ: (a) and (b) Γ = 0.1, 30, and 200 (red, green, and blue lines, respectively), N = 5000; (c) fc(r) (red line) and fp(r) (blue line) at Γ = 1000 and N = 14000 [PITH_FULL_IMAGE:figures/full_fig_p007_1.png] view at source ↗
Figure 3
Figure 3. Figure 3: Thermodynamic limit of the BOCP electrostatic [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figure 2
Figure 2. Figure 2: Size dependence of the BOCP total electrostatic [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗
Figure 4
Figure 4. Figure 4: Thermal fraction of the BOCP electrostatic energy [PITH_FULL_IMAGE:figures/full_fig_p010_4.png] view at source ↗
Figure 6
Figure 6. Figure 6: Size dependences of the excess interatomic elec [PITH_FULL_IMAGE:figures/full_fig_p010_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Excess interatomic electrostatic energy (circles) and [PITH_FULL_IMAGE:figures/full_fig_p011_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Size dependence of the ionic compressibility factor [PITH_FULL_IMAGE:figures/full_fig_p011_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Ionic compressibility factor Zi determined from MD simulation using Eqs. (16) (circles) and (24) (squares). Dashed curve indicates the approximation (43) of Zi from energies; solid line, Zi for the entire system (21); and the dashed-dotted line, Zi in the Debye–Hückel apprximation. IV. APPLICATION TO LAMMPS A. Simulation procedure A common way to calculate properties of an OCP system is to use the Ewald po… view at source ↗
Figure 10
Figure 10. Figure 10: Ionic compressibility factor Zi for OCP as a func￾tion of the LAMMPS cutoff length rc for the cases of random initialization (fluid) at Γ = 0.1 (circles), Γ = 10 (squares), Γ = 160 (diamonds), and Γ = 200 (triangles up) and ini￾tialization as bcc (solid) at Γ = 160 (triangles down) and Γ = 200 (stars). Dashed and dotted lines illustrate Zi ob￾tained for BOCP at Γ = 0.1 and 200, respectively. of coupling p… view at source ↗
Figure 12
Figure 12. Figure 12: Thermal fraction of the potential energy as a func [PITH_FULL_IMAGE:figures/full_fig_p014_12.png] view at source ↗
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.

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

1 major / 1 minor

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)
  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)
  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

1 responses · 0 unresolved

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
  1. 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

0 steps flagged

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

2 free parameters · 2 axioms · 0 invented entities

The central claims rest on the assumption that molecular dynamics with reflecting spherical boundaries produces size-dependent energies that extrapolate cleanly to the bulk limit and that standard classical statistical mechanics applies without additional corrections.

free parameters (2)
  • extrapolation function coefficients
    Coefficients in the fit of energy versus system size that are determined from the simulated data points.
  • proposed LAMMPS cutoff radius
    A specific numerical value chosen to improve agreement with the derived equation of state.
axioms (2)
  • domain assumption Classical molecular dynamics with reflecting boundaries faithfully samples the canonical ensemble for the BOCP.
    Invoked when the authors treat the MD trajectories as direct estimators of thermodynamic quantities.
  • domain assumption Finite-size corrections vanish as a smooth, extrapolatable function of radius.
    Required for the thermodynamic-limit estimates to be valid.

pith-pipeline@v0.9.0 · 5567 in / 1334 out tokens · 47117 ms · 2026-05-17T23:55:12.356263+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

72 extracted references · 72 canonical work pages

  1. [1]

    AtΓ>174, our data reach asymptotically the bcc energyu0, and in this region, they are in a very good agreement with [43] (Fig

    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. [2]

    Based on Eqs

    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. [3]

    The total number of particlesNwas fixed within each simulation, but varied between different simulations and the exact values will be specified later

    [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. [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)

  5. [5]

    S.Ichimaru,Stronglycoupledplasmas: high-densityclas- sical plasmas and degenerate electron liquids, Reviews of Modern Physics54, 1017 (1982)

  6. [6]

    E. V. Kholopov, Convergence problems of Coulomb and multipole sums in crystals, Physics-Uspekhi47, 965 (2004)

  7. [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)

  8. [8]

    Ichimaru, H

    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)

  9. [9]

    Chabrier, N

    G. Chabrier, N. W. Ashcroft, and H. E. DeWitt, White dwarfs as quantum crystals, Nature360, 48 (1992)

  10. [10]

    Chabrier, Quantum effects in dense coulumbic matter - application to the cooling of white dwarfs, The Astro- physical Journal414, 695 (1993)

    G. Chabrier, Quantum effects in dense coulumbic matter - application to the cooling of white dwarfs, The Astro- physical Journal414, 695 (1993)

  11. [11]

    Kozhberov, Thermal evolution of old white dwarfs, Journal of Physics: Conference Series929, 012012 (2017)

    A. Kozhberov, Thermal evolution of old white dwarfs, Journal of Physics: Conference Series929, 012012 (2017)

  12. [12]

    Blouin and J

    S. Blouin and J. Daligault, Phase separation in ultra- massive white dwarfs, The Astrophysical Journal919, 87 (2021)

  13. [13]

    D. A. Baiko, Phase diagrams of binary ionic mixtures and white dwarf cooling, Monthly Notices of the Royal Astronomical Society517, 3962 (2022)

  14. [14]

    Medin and A

    Z. Medin and A. Cumming, Crystallization of classical multicomponent plasmas, Physical Review E81, 036107 (2010)

  15. [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)

  16. [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)

  17. [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)

  18. [18]

    O. Arp, D. Block, A. Piel, and A. Melzer, Dust Coulomb balls: Three-dimensional plasma crystals, Physical Re- view Letters93, 165004 (2004)

  19. [19]

    O. Arp, D. Block, M. Klindworth, and A. Piel, Confine- ment of Coulomb balls, Physics of Plasmas12, 122102 (2005)

  20. [20]

    Käding and A

    S. Käding and A. Melzer, Three-dimensional stereoscopy of Yukawa (Coulomb) balls in dusty plasmas, Physics of Plasmas13, 090701 (2006)

  21. [21]

    Totsuji, T

    H. Totsuji, T. Ogawa, C. Totsuji, and K. Tsuruta, Struc- ture of spherical Yukawa clusters, Journal of Physics A: Mathematical and General39, 4545 (2006)

  22. [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)

  23. [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)

  24. [24]

    Huang, M

    D. Huang, M. Baggioli, S. Lu, Z. Ma, and Y. Feng, Re- vealing the supercritical dynamics of dusty plasmas and their liquidlike to gaslike dynamical crossover, Physical Review Research5, 013149 (2023)

  25. [25]

    M. A. Pokrant, Thermodynamic properties of the nonzero-temperature, quantum-mechanical, one- component plasma, Physical Review A16, 413 (1977)

  26. [26]

    Ceperley, Ground state of the fermion one-component plasma: A Monte Carlo study in two and three dimen- sions, Physical Review B18, 3126 (1978)

    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)

  27. [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)

  28. [28]

    Chugunov, D

    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)

  29. [29]

    Chugunov and D

    A. Chugunov and D. Baiko, Anharmonic corrections to the electrostatic energy of a Coulomb crystal, Physica A: Statistical Mechanics and its Applications352, 397 (2005)

  30. [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)

  31. [31]

    Hamaguchi and R

    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)

  32. [32]

    Hamaguchi, R

    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)

  33. [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)

  34. [34]

    J. M. Caillol and D. Gilles, Monte Carlo simulations of theYukawaone-componentplasma,JournalofStatistical Physics100, 933 (2000)

  35. [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)

  36. [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)

  37. [37]

    N. V. Brilliantov, Accurate first-principle equation of state for the one-component plasma, Contributions to Plasma Physics38, 489 (1998)

  38. [38]

    Chabrier and A

    G. Chabrier and A. Y. Potekhin, Equation of state of fully ionized electron-ion plasmas, Physical Review E58, 4941 (1998)

  39. [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)

  40. [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)

  41. [41]

    Brilliantov, V

    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)

  42. [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)

  43. [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)

  44. [44]

    J. P. Hansen, Statistical mechanics of dense ionized matter. I. Equilibrium properties of the classical one- component plasma, Physical Review A8, 3096 (1973)

  45. [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)

  46. [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)

  47. [47]

    Ogata and S

    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)

  48. [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)

  49. [49]

    Caillol and D

    J.-M. Caillol and D. Gilles, An accurate equation of state fortheone-componentplasmainthelowcouplingregime, Journal of Physics A: Mathematical and Theoretical43, 105501 (2010)

  50. [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)

  51. [51]

    R. T. Farouki and S. Hamaguchi, Thermal energy of the crystalline one-component plasma from dynamical simu- lations, Physical Review E47, 4330 (1993)

  52. [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)

  53. [53]

    J. D. Weeks, Volume change on melting for systems with inverse-power-law interactions, Physical Review B24, 1530 (1981)

  54. [54]

    Nordholm, Simple analysis of the thermodynamic properties of the one-component plasma, Chemical Physics Letters105, 302 (1984)

    S. Nordholm, Simple analysis of the thermodynamic properties of the one-component plasma, Chemical Physics Letters105, 302 (1984)

  55. [55]

    J. W. Lee and F. H. Ree, Perturbation theory of a classi- cal one-component plasma, Physical Review A38, 5714 (1988)

  56. [56]

    C. N. Likos and N. W. Ashcroft, Self-consistent theory of freezing of the classical one-component plasma, Physical Review Letters69, 316 (1992)

  57. [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)

  58. [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)

  59. [59]

    P. P. Ewald, Die Berechnung optischer und elektro- statischer Gitterpotentiale, Annalen der Physik369, 253 (1921)

  60. [60]

    Nijboer and F

    B. Nijboer and F. De Wette, On the calculation of lattice sums, Physica23, 309 (1957)

  61. [61]

    Karasawa and W

    N. Karasawa and W. A. Goddard, Acceleration of conver- genceforlatticesums,TheJournalofPhysicalChemistry 93, 7320 (1989)

  62. [62]

    Yakub and C

    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)

  63. [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)

  64. [64]

    D. I. Zhukhovitskii and E. E. Perevoshchikov, Structural transition in strongly coupled Coulomb clusters, High Temperature62, 421 (2024)

  65. [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)

  66. [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)

  67. [67]

    L. D. Landau and E. M. Lifshitz,Statistical Physics, Vol. 5 (Elsevier, Oxford, 1980)

  68. [68]

    Itoh and S

    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)

  69. [69]

    Hockney and J

    R. Hockney and J. Eastwood,Computer Simulation Us- ing Particles(CRC Press, 2021)

  70. [70]

    D. I. Zhukhovitskii, Multiscale approach to the theory of nonisothermal homogeneous nucleation, The Journal of Chemical Physics160, 194505 (2024)

  71. [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...

  72. [72]

    Bolten, F

    M. Bolten, F. Fahrenberger, R. Halver, F. Heber, M. Hofmann, I. Kabadshow, O. Lenz, M. Pippig, and G. Sutmann, ScaFaCoS, C subroutine library, http://scafacos.github.com/