Kirkwood-Buff integration (KBI) is implemented into the massively-parallel molecular simulation tool ms2 and assessed by molecular dynamics simulations of binary liquid mixtures. The formalism of Krüger et al. (P. Krüger et al., J. Phys. Chem. Lett. 4: 235–238, 2013) that adopts NVT ensemble data to the μVT ensemble is employed throughout. Taking advantage of its linear scaling with inverse system size, the extrapolation to the thermodynamic limit is analyzed. KBI are calculated with standard radial distribution functions (RDF) and two corrected RDF forms. Simulations in the NVT ensemble are carried out in the entire composition range for four Lennard-Jones mixtures, studying system size dependence by varying N = 4000, 8000 and 16000 molecules. Moreover, four mixtures of “real” components are considered with N = 4000. Thermodynamic factor, partial molar volumes and isothermal compressibility are calculated from KBI and compared with benchmark data from NpT ensemble simulations. The assessment shows that the formalism of Krüger et al. greatly improves KBI and that extrapolation is important, particularly for smaller systems.
Read full abstract