Monte Carlo simulations for the surface properties of the classical one-component plasma

Abstract
The authors present the results of Monte Carlo (MC) simulations for the surface properties of the classical one-component plasma. They have used a system of slab geometry for simulating a planar surface and made some improvements over previous work. They have devised a practical method of obtaining the surface energy in the weak coupling regime, where the particle configurations far outside the surface are quite difficult to simulate, and they have improved the sampling efficiency using a modified MC method, which is particularly useful and even indispensable for simulating strongly inhomogeneous systems. Their results for the surface energy show, on contrast to previous results, a proper behavior in the weak coupling regime and a systematic variation with increasing coupling constants, although the density profiles are hardly distinguishable from the previous results.