On a direct method of calculating pressure in the canonical ensemble

Determining fluid pressure in Monte Carlo simulations can be a challenging task due to the reduced arsenal of computational tools available to perform such measurements. In addition, none of these tools is general enough to calculate the equilibrium pressure for a wide variety of models. The Gibbs–D...

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Veröffentlicht in:The Journal of chemical physics 2024-02, Vol.160 (6)
Hauptverfasser: Tanouye, Fernando Takeshi, Alves, Jozismar Rodrigues
Format: Artikel
Sprache:eng
Schlagworte:
Online-Zugang:Volltext
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Beschreibung
Zusammenfassung:Determining fluid pressure in Monte Carlo simulations can be a challenging task due to the reduced arsenal of computational tools available to perform such measurements. In addition, none of these tools is general enough to calculate the equilibrium pressure for a wide variety of models. The Gibbs–Duhem method, for example, is a very useful option, but only for pure (one component) systems. To enrich this arsenal, we propose here a direct method to calculate pressure in the canonical ensemble, which could easily be extended to mixtures at low densities. In analogy with the Widom method, our approach is based on the free energy variation with volume, described in terms of the removal of an empty or particle-occupied lattice column. We tested our approach for the lattice gas model and compared the results with exact Onsager solutions. Furthermore, the appearance of thermodynamic instabilities (loops) in the pressure isotherms during the phase transition and their relationship with interface effects are discussed. Finally, a phase diagram is obtained from these isotherms using the Hill construction.
ISSN:0021-9606
1089-7690
DOI:10.1063/5.0187585