Efficient calculation of temperature dependence of solid-phase free energies by overlap sampling coupled with harmonically targeted perturbation

We examine a method for computing the change in free energy with temperature of a crystalline solid. In the method, the free-energy difference between nearby temperatures is calculated via overlap-sampling free-energy perturbation with Bennett's optimization. Coupled to this is a harmonically targeted perturbation that displaces the atoms in a manner consistent with the temperature change, such that for a harmonic system, the free-energy difference would be recovered with no error. A series of such perturbations can be assembled to bridge larger gaps in temperature. We test this harmonically targeted temperature perturbation (HTTP) method through the application to the inverse-power soft potential, u(r)=ε(σ/r)(n), over a range of temperatures up to the melting condition. Three exponent values (n=12, 9, and 6) for the potential are studied with different crystal structures, specifically face-centered cubic (fcc), body-centered cubic (bcc), and hexagonal close packing. Absolute free energies (classical only) for each system are obtained by implementing the series to near-zero temperature, where the harmonic model becomes very accurate. The HTTP method is shown to provide very precise results, with errors in the free energy smaller than two parts in 10(5). An analysis of the thermodynamic stability of the various structures in the infinite-system limit confirms previous findings. In particular, for n=12 and 9, the fcc structure is stable for all temperatures up to melting, and for n=6, the bcc crystal becomes stable relative to fcc for temperatures above kT/ε=0.802±0.001. The effects of vacancies and other defects are not considered in the analysis.