simulation of chemical systems undergoing any combination of reaction and phase
equilibria. The method requires only a knowledge of the species intermolecular potentials
and their ideal‐gas properties, in addition to specification of the system stoichiometry and
thermodynamic constraints. It avoids the calculation of chemical potentials and fugacities, as
is similarly the case for the Gibbs ensemble method for phase equilibrium simulations. The …