F. Calvo,J.P. Neirotti,David L. Freeman,J.D. Doll
Abstract:We study the 38-atom Lennard-Jones cluster with parallel tempering Monte Carlo methods in the microcanonical and molecular dynamics ensembles. A new Monte Carlo algorithm is presented that samples rigorously the molecular dynamics ensemble for a system at constant total energy, linear and angular momenta. By combining the parallel tempering technique with molecular dynamics methods, we develop a hybrid method to overcome quasi-ergodicity and to extract both equilibrium and dynamical properties from Monte Carlo and molecular dynamics simulations. Several thermodynamic, structural and dynamical properties are investigated for LJ$_{38}$, including the caloric curve, the diffusion constant and the largest Lyapunov exponent. The importance of insuring ergodicity in molecular dynamics simulations is illustrated by comparing the results of ergodic simulations with earlier molecular dynamics simulations.
What problem does this paper attempt to address?
The paper attempts to address the problem of how to effectively simulate the phase transition process and extract the equilibrium and dynamic properties in a Lennard-Jones (LJ) cluster composed of 38 atoms. Specifically, the paper focuses on how to overcome the quasi-ergodicity problem, which occurs when the system easily gets trapped in local energy minima during simulations on a complex potential energy surface, leading to the failure of the ergodic hypothesis on the simulation time scale.
To achieve this goal, the authors propose the following methods:
1. **Parallel Tempering**:
- Apply the parallel tempering method in the microcanonical ensemble and molecular dynamics ensemble, promoting the system to cross high energy barriers by exchanging configurations at different temperatures.
- This method ensures ergodicity in the microcanonical ensemble, which is a more challenging goal compared to achieving it in the canonical ensemble.
2. **Hybrid MC/MD Method**:
- Combine the Monte Carlo (MC) method and the Molecular Dynamics (MD) method to develop a new algorithm that can strictly sample the molecular dynamics ensemble under fixed total energy, linear momentum, and angular momentum.
- This method allows the extraction of the system's equilibrium and dynamic properties, such as heat capacity curves, diffusion constants, and the maximum Lyapunov exponent.
3. **Specific Application**:
- The paper applies these methods to a Lennard-Jones cluster of 38 atoms, studying its solid-solid and solid-liquid phase transition processes.
- By comparing the results of ergodic and non-ergodic simulations, the importance of ensuring ergodicity is demonstrated.
In summary, the main goal of the paper is to develop and validate a simulation method that can effectively achieve ergodicity in complex systems, thereby accurately extracting the thermodynamic and dynamic properties of the system.