A new method of semigrand canonical ensemble to calculate first-order phase transitions for binary mixtures

First-order phase transitions of binary mixtures at the given pressure (P) and temperature (T) are studied by taking into account the composition fluctuations. Isothermal-isobaric semigrand canonical ensemble is adopted to find the relations among the total number of molecules, the composition fluctuations and Gibbs free energy density. By combining two identical subsystems of mixtures successively, the free energy density is transformed until being stable and its linear segments represent phase transitions. A new method is developed to calculate the phase equilibriums of binary mixtures. The method handles multiple types and number of phase equilibriums at single time and its solutions are physically justified. One example is shown for calculating the phase diagram of binary Lennard-Jones mixture. It demonstrates that the fluctuations of the total number of molecules in mixtures are fundamental behind phase transitions and the van der Waals loops in Gibbs free energy are reasonable.