Monte Carlo simulations of glass-forming liquids beyond Metropolis

Ludovic Berthier,Federico Ghimenti Frédéric van Wijland
DOI: https://doi.org/10.1063/5.0225978
2024-09-20
Abstract:Monte Carlo simulations are widely employed to measure the physical properties of glass-forming liquids in thermal equilibrium. Combined with local Monte Carlo moves, the Metropolis algorithm can also be used to simulate the relaxation dynamics, thus offering an efficient alternative to molecular dynamics. Monte Carlo simulations are however more versatile, because carefully designed Monte Carlo algorithms can more efficiently sample the rugged free energy landscape characteristic of glassy systems. After a brief overview of Monte Carlo studies of glass-formers, we define and implement a series of Monte Carlo algorithms in a three-dimensional model of polydisperse hard spheres. We show that the standard local Metropolis algorithm is the slowest, and that implementing collective moves or breaking detailed balance enhances the efficiency of the Monte Carlo sampling. We use time correlation functions to provide a microscopic interpretation of these observations. Seventy years after its invention, the Monte Carlo method remains the most efficient and versatile tool to compute low-temperatures properties in supercooled liquids.
Statistical Mechanics,Disordered Systems and Neural Networks,Soft Condensed Matter,Chemical Physics,Computational Physics
What problem does this paper attempt to address?
The paper aims to address the issue of simulating the physical properties of glass-forming liquids in a supercooled state. Specifically, the researchers use the Monte Carlo (MC) simulation method to measure the physical properties of these systems in thermal equilibrium and explore how to accelerate the sampling of complex free energy landscapes by designing more efficient Monte Carlo algorithms. Although the traditional Metropolis algorithm can be used to simulate the dynamical relaxation process of the system, its efficiency is relatively low. Therefore, this paper proposes a series of improved Monte Carlo algorithms, such as collective moves or breaking detailed balance, and provides a microscopic explanation through time correlation functions, demonstrating that these methods can significantly enhance the efficiency of Monte Carlo sampling. Additionally, the study shows that introducing diameter swaps can greatly expand the density range at which the system reaches equilibrium, making it possible to perform simulations at higher densities. In summary, the goal of this paper is to achieve faster and more efficient Monte Carlo simulations in supercooled liquids to better understand and predict the physical properties of glass materials.