Direct calculation of solid-vapor coexistence points by thermodynamic integration: application to single component and binary systems

We present a new thermodynamic integration method that directly connects the vapor and solid phases by a reversible path. The thermodynamic integration in the isothermal-isobaric ensemble yields the Gibbs free energy difference between the two phases, from which the sublimation temperature can be easily calculated. The method extends to the binary mixture without any modification to the integration path simply by employing the isothermal-isobaric semigrand ensemble. The thermodynamic integration, in this case, yields the chemical potential difference between the solid and vapor phases for one of the components, from which the binary sublimation temperature can be calculated. The coexistence temperatures predicted by our method agree well with those in the literature for single component and binary Lennard-Jones systems.