System-size dependence of the free energy of crystalline solids.

We investigate the system-size dependence of the Helmholtz free energy of crystalline solids from computer simulation. We employ a standard thermodynamic integration technique along a reversible path that links the crystalline solid with a noninteracting Einstein crystal with the same structure. The key contribution to the free energy is computed by using the so-called expanded-ensemble technique and the results are compared with those obtained from conventional integration of the derivative of the free energy along the path using Gaussian-Legendre quadrature.