Generalized Replica Exchange Method Research Articles

Studying vapor-liquid transition using a generalized ensemble.

Homogeneous vapor-liquid nucleation is studied using the generalized Replica Exchange Method (gREM). The generalized ensemble allows the study of unstable states that cannot directly be studied in the canonical ensemble. Along with replica exchange, this allows for efficient sampling of the multiple states in a single simulation. Statistical Temperature Weighted Histogram Analysis Method is used for postprocessing to get a continuous free energy curve from bulk vapor to bulk liquid. gREM allows the study of planar, cylindrical, and spherical interfaces in a single simulation. The excess Gibbs free energy for the formation of a spherical liquid droplet in vapor for a Lennard-Jones system is calculated from the free energy curve and compared against the umbrella sampling results. The nucleation free energy barrier obtained from gREM is then used to calculate the nucleation rate without relying on any classification scheme for separating the vapor and liquid.

Enhanced Sampling of Phase Transitions in Coarse-Grained Lipid Bilayers.

Freezing and melting of dipalmitoylphosphatidylcholine (DPPC) bilayers are simulated in both the explicit (Wet) and implicit solvent (Dry) coarse-grained MARTINI force fields with enhanced sampling, via the isobaric, molecular dynamics version of the generalized replica exchange method (gREM). Phase transitions are described with the entropic viewpoint, based upon the statistical temperature as a function of enthalpy, TS(H) = 1/(dS(H)/dH), where S is the configurational entropy. Bilayer thickness, area per lipid, and the second-rank order parameter (P2) are calculated vs temperature in the transition range. In a 32-lipid Wet MARTINI system, transitions in the lipid and water subsystems are strongly coupled, giving rise to considerable structure in TS(H) and the need to specify the state of the water when reporting a lipid transition temperature. For gel lipid + liquid water → fluid lipid + liquid water, we find 292.4 K. The small system is influenced by finite-size effects, but it is argued that the entropic approach is well suited to revealing them, which will be particularly relevant for studies of finite nanosystems where there is no thermodynamic limit. In a 390-lipid Dry MARTINI system, two-dimensional analogues of the topographies of coexisting states ("subphases") seen in pure fluids are found. They are not seen in the 32-lipid Wet or Dry system, but the Dry lipids show a new type of state with gel in one leaflet and tilted gel in the other. Dry bilayer transition temperatures are 333.3 K (390 lipids) and 338 K (32 lipids), indicating that the 32-lipid system is not too small for a qualitative study of the transition. Physical arguments are given for Dry lipid system size dependence and for the difference between Wet and Dry systems.

Freezing Transitions of Nanoconfined Coarse-Grained Water Show Subtle Dependence on Confining Environment.

The solid-to-liquid phase transition in water nanofilms confined between plates, with varying separations and water-plate interactions ranging from strongly hydrophobic to strongly hydrophilic, was simulated using a coarse-grained monatomic water model (mW) and the generalized replica exchange method (gREM). Extensive gREM simulations combined with the statistical temperature weighted histogram analysis method (ST-WHAM) provide a detailed description of the thermodynamic properties intrinsic to the phase transition, including the transition temperature, isobaric heat capacity, phase change enthalpy, entropy, and their dependence on the interplate distance and the plate-water interaction. The ice structure of water nanofilms was characterized at various conditions using the transverse density profile and the distribution of angles formed by hydrogen-bonded neighboring molecules. Flat bilayer ice was observed to be the dominant solid phase at close interplate distance, while puckered bilayer ice, similar to a slab of ice Ih, is the predominant structure at larger interplates. Stable puckered bilayer ice, previously observed to have a low melting point, is observed to have enhanced stability with high melting temperature when confined between hydrophilic plates. These results demonstrate the strong dependence of phase stability and coexistence in nanoconfined systems on the geometry and physical properties of the confining environment.

Water Freezing and Ice Melting.

The generalized replica exchange method (gREM) is designed to sample states with coexisting phases and thereby to describe strong first order phase transitions. The isobaric MD version of the gREM is presented and applied to the freezing of liquid water and the melting of hexagonal and cubic ice. It is confirmed that coexisting states are well-sampled. The statistical temperature as a function of enthalpy, TS(H), is obtained. Hysteresis between freezing and melting is observed and discussed. The entropic analysis of phase transitions is applied and equilibrium transition temperatures, latent heats, and surface tensions are obtained for hexagonal ice ↔ liquid and cubic ice ↔ liquid with excellent agreement with published values. A new method is given to assign water molecules among various symmetry types. Pathways for water freezing, ultimately leading to hexagonal ice, are found to contain intermediate layered structures built from hexagonal and cubic ice.

Isobaric Molecular Dynamics Version of the Generalized Replica Exchange Method (gREM): Liquid–Vapor Equilibrium

A prescription for sampling isobaric generalized ensembles with molecular dynamics is presented and applied to the generalized replica exchange method (gREM), which was designed to simulate first-order phase transitions. The properties of the isobaric gREM ensemble are discussed, and a study is presented for the liquid-vapor equilibrium of the guest molecules given for gas hydrate formation with the mW water model. Phase diagrams, critical parameters, and a law of corresponding states are obtained.

Generalized ensemble method applied to study systems with strong first order transitions

At strong first-order phase transitions, the entropy versus energy or, at constant pressure, enthalpy, exhibits convex behavior, and the statistical temperature curve correspondingly exhibits an S-loop or back-bending. In the canonical and isothermal-isobaric ensembles, with temperature as the control variable, the probability density functions become bimodal with peaks localized outside of the S-loop region. Inside, states are unstable, and as a result simulation of equilibrium phase coexistence becomes impossible. To overcome this problem, a method was proposed by Kim, Keyes and Straub [1], where optimally designed generalized ensemble sampling was combined with replica exchange, and denoted generalized replica exchange method (gREM). This new technique uses parametrized effective sampling weights that lead to a unimodal energy distribution, transforming unstable states into stable ones. In the present study, the gREM, originally developed as a Monte Carlo algorithm, was implemented to work with molecular dynamics in an isobaric ensemble and coded into LAMMPS, a highly optimized open source molecular simulation package. The method is illustrated in a study of the very strong solid/liquid transition in water.

Investigating the solid-liquid phase transition of water nanofilms using the generalized replica exchange method.

The generalized Replica Exchange Method (gREM) was applied to study a solid-liquid phase transition in a nanoconfined bilayer water system using the monatomic water (mW) model. Exploiting optimally designed non-Boltzmann sampling weights with replica exchanges, gREM enables an effective sampling of configurations that are metastable or unstable in the canonical ensemble via successive unimodal energy distributions across phase transition regions, often characterized by S-loop or backbending in the statistical temperature. Extensive gREM simulations combined with Statistical Temperature Weighted Histogram Analysis Method (ST-WHAM) for nanoconfined mW water at various densities provide a comprehensive characterization of diverse thermodynamic and structural properties intrinsic to phase transitions. Graph representation of minimized structures of bilayer water systems determined by the basin-hopping global optimization revealed heterogeneous ice structures composed of pentagons, hexagons, and heptagons, consistent with an increasingly ordered solid phase with decreasing density. Apparent crossover from a first-order solid-liquid transition to a continuous one in nanoconfined mW water with increasing density of the system was observed in terms of a diminishing S-loop in the statistical temperature, smooth variation of internal energies and heat capacities, and a characteristic variation of lateral radial distribution functions, and transverse density profiles across transition regions.

Order parameter free enhanced sampling of the vapor-liquid transition using the generalized replica exchange method

The generalized Replica Exchange Method (gREM) is extended into the isobaric-isothermal ensemble, and applied to simulate a vapor-liquid phase transition in Lennard-Jones fluids. Merging an optimally designed generalized ensemble sampling with replica exchange, gREM is particularly well suited for the effective simulation of first-order phase transitions characterized by "backbending" in the statistical temperature. While the metastable and unstable states in the vicinity of the first-order phase transition are masked by the enthalpy gap in temperature replica exchange method simulations, they are transformed into stable states through the parameterized effective sampling weights in gREM simulations, and join vapor and liquid phases with a succession of unimodal enthalpy distributions. The enhanced sampling across metastable and unstable states is achieved without the need to identify a "good" order parameter for biased sampling. We performed gREM simulations at various pressures below and near the critical pressure to examine the change in behavior of the vapor-liquid phase transition at different pressures. We observed a crossover from the first-order phase transition at low pressure, characterized by the backbending in the statistical temperature and the "kink" in the Gibbs free energy, to a continuous second-order phase transition near the critical pressure. The controlling mechanisms of nucleation and continuous phase transition are evident and the coexistence properties and phase diagram are found in agreement with literature results.

Exploring the Solid–Liquid Phase Change of an Adapted Dzugutov Model Using Generalized Replica Exchange Method

The generalized replica exchange method (gREM) is applied to simulate a solid-liquid phase change in an adapted Dzugutov model system. The idea of gREM is to incorporate the merit of sampling from a generalized ensemble into the replica exchange paradigm. The generalized ensemble weights are determined from effective temperatures tailored for this system, through an inverse mapping strategy. The ordered and disordered phases are smoothly joined together through a succession of unimodal energy distributions, providing a continuous transformation from an isotropic liquid into an amorphous solid phase under the given condition. The phase change temperature is determined by evaluation of various thermodynamic and structural order parameters. The increasing tendency for icosahedral local order under cooling is quantified using bond order parameters. The ergodic measure is used to demonstrate the advantage of gREM over the traditional temperature-based REM in sampling the phase change region.

Generalized replica exchange method.

We present a powerful replica exchange method, particularly suited to first-order phase transitions associated with the backbending in the statistical temperature, by merging an optimally designed generalized ensemble sampling with replica exchanges. The key ingredients of our method are parametrized effective sampling weights, smoothly joining ordered and disordered phases with a succession of unimodal energy distributions by transforming unstable or metastable energy states of canonical ensembles into stable ones. The inverse mapping between the sampling weight and the effective temperature provides a systematic way to design the effective sampling weights and determine a dynamic range of relevant parameters. Illustrative simulations on Potts spins with varying system size and simulation conditions demonstrate a comprehensive sampling for phase-coexistent states with a dramatic acceleration of tunneling transitions. A significant improvement over the power-law slowing down of mean tunneling times with increasing system size is obtained, and the underlying mechanism for accelerated tunneling is discussed.