Coarse-grained molecular dynamics modeling of strongly associating fluids: thermodynamics, liquid structure, and dynamics of symmetric binary mixture fluids

Symmetric binary mixtures capable of strong association via a highly directional and saturable specific interaction between unlike molecules are investigated by canonical molecular dynamics simulations. The specific interaction of the molecules is defined in a new coarse-grained pair potential that can be applied in continuous molecular dynamics as well as in Monte Carlo simulations. The thermodynamic, structural, and dynamic properties of the associating mixture fluids are investigated as a function of density, temperature, and association strength of the specific interaction. Detailed analysis of the simulation data confirms a two-stage mechanism in the formation of specific bonds with increasing interaction strength, including a fast dimerization process and a subsequent stage of perfecting the bonds. A large heat capacity peak is found during the formation or breaking of the bonds, reflecting the large energy fluctuation introduced by the strong association. The fractions of nonbonded molecules obtained from the simulations as a function of density, temperature, and interaction strength are in excellent agreement with the predictions of Wertheim's thermodynamic perturbation theory. The translational and rotational dynamics of the Tmer mixture are effectively retarded with increasing association strength and are analyzed in terms of autocorrelation functions and a non-Gaussian parameter for the translational dynamics. The lifetimes of molecules in bonded and nonbonded states provide detailed information about the transformation of molecules between the bonded state and the nonbonded state. Finally, simulation sampling problems inherent to strongly interacting systems are easily overcome using the parallel tempering simulation technique. This latter result confirms that with the new continuous coarse-grained simulation potential we have a versatile and flexible interaction potential that can be used with many available molecular dynamics and Monte Carlo algorithms under various ensembles.