Heat capacities of supercritical fluids via Grand Canonical ensemble simulations

In this work, suitable mathematical relationships to compute isobaric heat capacities from molecular simulations in the Grand Canonical (GC) ensemble are derived and tested via Monte Carlo methods. Using atomistic classical force fields, the residual isobaric heat capacities of pure carbon dioxide (CO₂) and pure methanol (MeOH) were obtained at supercritical conditions (with critical properties estimated from a finite-size scaling analysis). The total isobaric heat capacity was determined by combining the residual isobaric heat capacity obtained from molecular simulations with the ideal gas contributions obtained from experimental correlations. Isobaric heat capacities generated from both GC and Isothermal–Isobaric ensemble simulations were compared to predictions from accurate equations of state (EOS)s for CO₂ and MeOH at corresponding reduced temperatures and pressures. Isobaric heat capacities calculated from both ensembles were in good agreement with those obtained from the Span and Wagner EOS for CO₂ and the IUPAC EOS for MeOH. For comparable computation times, simulations run in the GC ensemble generate results with significantly lower statistical uncertainty than those run in the Isothermal–Isobaric ensemble.