Molecular simulation of volume of mixing,Helmholtz free energy of mixing and entropy of mixing in bulk fluid mixtures

We describe a kinetic Monte Carlo molecular simulation procedure to calculate the Helmholtz free energy, the entropy and the chemical potentials of all components in a bulk fluid mixture. This allows us to derive the excess properties (volume, free energy and entropy) resulting from the mixing of homogeneous fluids of pure components at constant temperature and pressure. We have chosen neon–xenon mixtures to illustrate our method because of the large difference in collision diameter and well-depth of the interaction energy. When xenon is predominant in the mixture, the volume of mixing is larger. The excess entropy of mixing correlates with the volume of mixing, since a positive excess volume enables more configurations (more possible molecular distributions). The excess thermodynamic quantities as functions of the total density were found to be insensitive to temperature. To investigate the effects of the molecular parameters, we also studied argon–nitrogen and argon–krypton mixtures. The effect of the difference in molecular parameters is in the order: argon–nitrogen < argon–krypton < neon–xenon. A large difference in the well-depth of the interaction energies results in an increase in the excess thermodynamic variables, which is in agreement with the literature McDonald IR. NpT-ensemble Monte Carlo calculations for binary liquid mixtures. Mol Phys. 1972;23(1):41–58; Singer JVL, Singer K. Monte Carlo calculation of thermodynamic properties of binary mixtures of Lennard-Jones (12-6) liquids. Mol Phys. 1972;24(2):357–390.