Journal
JOURNAL OF CHEMICAL PHYSICS
Volume 140, Issue 21, Pages -Publisher
AMER INST PHYSICS
DOI: 10.1063/1.4878665
Keywords
-
Funding
- (U.S.) Department of Energy (DOE) at the Lawrence Livermore National Laboratory [DE-AC52-07NA27344]
- National Science Foundation (NSF) [NSF CHE-1301314, NSF CHE-1012545, DMS07-08140]
- Office of Naval Research (ONR) [N00014-11-1-0345]
- National Institutes of Health (NIH) [R01GM100472]
- Direct For Mathematical & Physical Scien
- Division Of Chemistry [1301314] Funding Source: National Science Foundation
- Direct For Mathematical & Physical Scien
- Division Of Materials Research [1207432] Funding Source: National Science Foundation
Ask authors/readers for more resources
The problem of predicting polymorphism in atomic and molecular crystals constitutes a significant challenge both experimentally and theoretically. From the theoretical viewpoint, polymorphism prediction falls into the general class of problems characterized by an underlying rough energy landscape, and consequently, free energy based enhanced sampling approaches can be brought to bear on the problem. In this paper, we build on a scheme previously introduced by two of the authors in which the lengths and angles of the supercell are targeted for enhanced sampling via temperature accelerated adiabatic free energy dynamics [T. Q. Yu and M. E. Tuckerman, Phys. Rev. Lett. 107, 015701 (2011)]. Here, that framework is expanded to include general order parameters that distinguish different crystalline arrangements as target collective variables for enhanced sampling. The resulting free energy surface, being of quite high dimension, is nontrivial to reconstruct, and we discuss one particular strategy for performing the free energy analysis. The method is applied to the study of polymorphism in xenon crystals at high pressure and temperature using the Steinhardt order parameters without and with the supercell included in the set of collective variables. The expected fcc and bcc structures are obtained, and when the supercell parameters are included as collective variables, we also find several new structures, including fcc states with hcp stacking faults. We also apply the new method to the solid-liquid phase transition in copper at 1300 K using the same Steinhardt order parameters. Our method is able to melt and refreeze the system repeatedly, and the free energy profile can be obtained with high efficiency. (C) 2014 AIP Publishing LLC.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available