A general method for calculating solid/liquid interfacial free energies from atomistic simulations: application to CaSO4.xH2O

We present a general method for computing interfacial free energies from atomistic simulations, which is particularly suitable for solid/liquid interfaces. Our method uses an Einstein crystal as a universal reference state and is more flexible than previous approaches. Surfaces with dipoles, complex reconstructions, and partially dissolved species are all easily accommodated within the framework. It may also be extended to calculating the relative free energies of different phases and other types of defect. We have applied our method to interfaces of bassanite and gypsum with water and obtained interfacial free energies of the order of 0.15 J/m2, of which approximately 50 % is due to entropic contributions. Our calculations of the interfacial free energy of NaCl with water obtained a value of 0.13 J/m2 of which only 19 % is from entropic contributions. We have also predicted equilibrium morphologies for bassanite and gypsum that compare well with experiments and previous calculations.