On the computational calculation of surface free energies for the disordered semihexagonal reconstructed Au(100) surface

Previously we developed a general method for calculating the free energy of any surface constrained to a distinct surface excess number/density. In this paper we show how to combine a range of such surfaces, whose free energies have been calculated, to produce an ad hoc semigrand canonical ensemble of surfaces from which ensemble surface properties can be calculated, including the ensemble surface free energy. We construct such an ensemble for the disordered Au(100) semihexagonal reconstructed surface using a Glue model potential at 1000 K and calculate the ensemble surface free energy to be 0.088 18 eVA(2). The ensemble average surface lateral density was found to be 1.375 (with respect to the bulk), which is in agreement with previous grand canonical Monte Carlo studies.