One--Component Plasma Equation of State Revisited via Angular--Averaged Ewald Potential

We present analytic fits of classical one–component plasma (OCP) internal energy over a wide range of coupling parameter $0.01\le \mathrm{\Gamma }\le 170$ using Monte–Carlo data in the thermodynamic limit. We extend the dataset obtained in [Demyanov and Levashov, Phys. Rev. E 106, 015204 (2022)] using the angular–averaged Ewald potential with additional points at strong coupling ($\mathrm{\Gamma }=120,150,170$). We then fit two frequently used functional forms for the OCP equation of state: (i) a five-parameter equation by Caillol [J. Chem. Phys. 111, 6538–6547 (1999)] and (ii) the equation by Potekhin and Chabrier [Phys. Rev. E 62, 8554 (2000)] that enforces the Debye–Hückel limit. The presented fits reproduce our MC data within statistical uncertainties, recovering the correct weak-coupling behavior. Coefficients, recommended validity ranges, and comparisons to prior analytical and simulation results are provided.