skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

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

Abstract

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, 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. Lastly, the method is illustrated in a study of the very strong solid/liquid transition in water.

Authors:
 [1];  [2];  [1]
  1. Boston Univ., Boston, MA (United States)
  2. Broad Institute of Massachusetts Inst. of Technology (MIT), Cambridge, MA (United States); Harvard Univ., Cambridge, MA (United States)
Publication Date:
Research Org.:
Boston Univ., MA (United States)
Sponsoring Org.:
USDOE
Contributing Org.:
Boston University Center for Computational Science
OSTI Identifier:
1247484
Grant/Contract Number:
SC0008810
Resource Type:
Journal Article: Accepted Manuscript
Journal Name:
Journal of Physics. Conference Series
Additional Journal Information:
Journal Volume: 640; Journal Issue: 2; Conference: XXVI IUPAP Conference on Computational Physics (CCP2014) , Boston, MA (United States), 11-14 Aug 2014; Journal ID: ISSN 1742-6588
Publisher:
IOP Publishing
Country of Publication:
United States
Language:
English
Subject:
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CHEMISTRY; 36 MATERIALS SCIENCE; 97 MATHEMATICS AND COMPUTING; phase transitions; computer simulation; algorithm design; water; freezing; melting

Citation Formats

Malolepsza, E., Kim, J., and Keyes, T. Generalized ensemble method applied to study systems with strong first order transitions. United States: N. p., 2015. Web. doi:10.1088/1742-6596/640/1/012003.
Malolepsza, E., Kim, J., & Keyes, T. Generalized ensemble method applied to study systems with strong first order transitions. United States. doi:10.1088/1742-6596/640/1/012003.
Malolepsza, E., Kim, J., and Keyes, T. Mon . "Generalized ensemble method applied to study systems with strong first order transitions". United States. doi:10.1088/1742-6596/640/1/012003. https://www.osti.gov/servlets/purl/1247484.
@article{osti_1247484,
title = {Generalized ensemble method applied to study systems with strong first order transitions},
author = {Malolepsza, E. and Kim, J. and Keyes, T.},
abstractNote = {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, 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. Lastly, the method is illustrated in a study of the very strong solid/liquid transition in water.},
doi = {10.1088/1742-6596/640/1/012003},
journal = {Journal of Physics. Conference Series},
number = 2,
volume = 640,
place = {United States},
year = {Mon Sep 28 00:00:00 EDT 2015},
month = {Mon Sep 28 00:00:00 EDT 2015}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Save / Share:
  • We have generalized the Gutzwiller method to the cases of n{ge}2 correlated bands and report studies on a degenerate two-band model with Hund{close_quote}s rule type on-site interactions. At half band filling the metal-insulator transitions are usually of first order. {copyright} {ital 1997} {ital The American Physical Society}
  • Relying on the recent proposed multicanonical algorithm, we present a numerical simulation of the first-order phase transition in the 2D 10-state Potts model on lattices up to sizes 100{times}100. It is demonstrated that the new algorithm {ital lacks} an exponentially fast increase of the tunneling time between metastable states as a function of the linear size {ital L} of the system. Instead, the tunneling time diverges approximately proportional to {ital L}{sup 2.65}. On our largest lattice we gain more than 2 orders of magnitude as compared to a standard heat-bath algorithm. As a first physical application we report a high-precisionmore » computation of the interfacial free energy per unit area.« less
  • We present a generalized circle theorem which includes the Lee-Yang theorem for symmetric transitions as a special case. It is found that zeros of the partition function can be written in terms of discontinuities in the derivatives of the free energy. For asymmetric transitions, the locus of the zeros in tangent to the unit circle if the partition functions of the two phases are added up with unequal prefactors. This conclusion is substantiated by explicit calculation of zeroes of the partition function for the Blume-Capel model near and at the triple line at low temperatures.
  • We consider a spin system with nearest-neighbor antiferromagnetic pair interactions in a two-dimensional lattice. We prove that the free energy of this system is differentiable with respect to the uniform external field h, for all temperatures and all h. This implies the absence of a first-order phase transition in this system.