Heterogeneous ice nucleation is ubiquitous but its microscopic mechanisms can be extraordinarily complex even on a simple surface. Such complexity poses a challenge in modeling nucleation using advanced sampling methods. Here, we investigate heterogeneous ice nucleation on an FCC (211) surface by a forward flux sampling (FFS) method to understand how the complexity in nucleation pathways affects the accuracy of FFS. We first show the commonly adopted, size-based order parameter fails to describe heterogeneous ice nucleation on the FCC (211) surface. Inclusion of geometric anisotropy of ice nucleus as an additional descriptor is found to significantly improve the quality of the size-based order parameter for the current system. Subsequent application of this new order parameter in FFS identifies two competing ice nucleation pathways in the system: a primary-prism-planed (PPP) path and a secondary-prism-planed (SPP) path, both leading to the formation of hexagonal ice but with different crystalline orientations. Although the PPP pathway dominates ice nucleation on the FCC (211) surface, the occurrence of the less efficient SPP pathway, particularly its strong presence at the early stage of FFS, is found to yield a significant statistical uncertainty in the calculated FFS rate constant. We develop a two-path model that enables gaining a general, quantitative understanding of the impact of initial finite sampling on the reliability of FFS calculations in the presence of multiple nucleation pathways. Our study also suggests a few general strategies for improving the accuracy of FFS when exploring unknown but complex systems.
more »
« less
Practical guide to replica exchange transition interface sampling and forward flux sampling
Path sampling approaches have become invaluable tools to explore the mechanisms and dynamics of the so-called rare events that are characterized by transitions between metastable states separated by sizable free energy barriers. Their practical application, in particular to ever more complex molecular systems, is, however, not entirely trivial. Focusing on replica exchange transition interface sampling (RETIS) and forward flux sampling (FFS), we discuss a range of analysis tools that can be used to assess the quality and convergence of such simulations, which is crucial to obtain reliable results. The basic ideas of a step-wise evaluation are exemplified for the study of nucleation in several systems with different complexities, providing a general guide for the critical assessment of RETIS and FFS simulations.
more » « less- Award ID(s):
- 2224643
- NSF-PAR ID:
- 10367712
- Publisher / Repository:
- American Institute of Physics
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 156
- Issue:
- 20
- ISSN:
- 0021-9606
- Page Range / eLocation ID:
- Article No. 200901
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
more » « less
Summary The evaluation of functionality and its evolution in the aftermath of extreme events and during the restoration phase is a critical step in disaster resilience assessment. To this respect, this paper presents the ‘
Functionality‐Fragility Surface ’ (FFS ), which is a tool for probabilistic functionality and resilience evaluation of damaged structures, infrastructure systems, and communities.FFS integrates two well‐known tools, namely Fragility Curves and Restoration Functions, to present the probability of loss of functionality of a system as a function of the extreme‐event intensity, as well as the elapsed time from the initiation of the restoration process. Because of their versatility,FFS s can be applied to components and systems belonging to different infrastructure sectors (e.g., transportation, power distribution, and telecommunication), so they provide a common rigorous paradigm for integrated resilience analyses of multiple sectors, as well as for studies on interdependencies within and across sectors. While it is shown thatFFS s can be developed using available data and simple computations for different types of structures and infrastructure systems, this paper proposes also a sophisticated simulation‐based methodology to developFFS s for individual bridges, taking into account the uncertainties involved in the response, damage, and restoration scheduling of bridges. A Multi‐Span Simply Supported Steel Girder bridge is used to showcase the application of the proposed methodology. Copyright © 2017 John Wiley & Sons, Ltd. -
more » « less
Cation‐π interactions are noncovalent interactions between a π‐electron system and a positively charged ion that are regarded as a strong noncovalent interaction and are ubiquitous in biological systems. Similarly, though less studied, anion‐ring interactions are present in proteins along with in‐plane interactions of anions with aromatic rings. As these interactions are between a polarizing ion and a polarizable π system, the accuracy of the treatment of these interactions in molecular dynamics (MD) simulations using additive force fields (FFs) may be limited. In the present work, to allow for a better description of ion‐π interactions in proteins in the Drude‐2013 protein polarizable FF, we systematically optimized the parameters for these interactions targeting model compound quantum mechanical (QM) interaction energies with atom pair‐specific Lennard‐Jones parameters along with virtual particles as selected ring centroids introduced to target the QM interaction energies and geometries. Subsequently, MD simulations were performed on a series of protein structures where ion‐π pairs occur to evaluate the optimized parameters in the context of the Drude‐2013 FF. The resulting FF leads to a significant improvement in reproducing the ion‐π pair distances observed in experimental protein structures, as well as a smaller root‐mean‐square differences and fluctuations of the overall protein structures from experimental structures. Accordingly, the optimized Drude‐2013 protein polarizable FF is suggested for use in MD simulations of proteins where cation‐π and anion‐ring interactions are critical. © 2019 Wiley Periodicals, Inc.
-
more » « less
Ergodicity (or at least the tantalizing promise of it) is a core animating principle of molecular-dynamics (MD) simulations: Put simply, sample for long enough (in time), and you will make representative visits to states of a system all throughout phase space, consistent with the desired statistical ensemble. However, one is not guaranteed a priori that the chosen window of sampling in a production run is sufficiently long to avoid problematically non-ergodic observations; one is also not guaranteed that successive measurements of an observable are statistically independent of each other. In this paper, we investigate several particularly striking and troublesome examples of statistical correlations in MD simulations of nanoconfined fluids, which have profound implications on the quantification of uncertainty for transport phenomena in these systems. In particular, we show that these correlations can lead to confidence intervals on the fluid self-diffusion coefficient that are dramatically overconfident and estimates of this transport quantity that are simply inaccurate. We propose a simple approach—based on the thermally accelerated decorrelation of fluid positions and momenta—that ameliorates these issues and improves our confidence in MD measurements of nanoconfined fluid transport properties. We demonstrate that the formation of faithful confidence intervals for measurements of self-diffusion under nanoscale confinement typically requires at least 20 statistically independent samples, and potentially more depending on the sampling technique used.
-
more » « less
Abstract The high cost of ab initio molecular dynamics (AIMD) simulations to model complex physical and chemical systems limits its ability to address many key questions. However, new machine learning‐based representations of complex potential energy surfaces have been introduced in recent years to circumvent computationally demanding AIMD simulations while retaining the same level of accuracy. As these machine learning methods gain in popularity over the next decade, it is important to address the appropriate way to develop and integrate them with well‐established simulation methods. This paper details the parameterization and training, using accurate electronic structure calculations, of artificial neural network potentials (NNPs) to model the intermolecular interactions in a hydrogen clathrate hydrate, wherein hydrogen molecules are confined inside the cavities of the 3D crystalline framework of hydrogen‐bonded water molecules. This new NNP is used in conjunction with new path‐integral based enhanced sampling methods, for inclusion of nuclear quantum effects and promotion of free energy barrier crossing, in order to determine properties that are important for understanding the diffusion of hydrogen gas in the clathrate system. These simulations demonstrate the influence of cage occupancy on the free energy barriers that determine the diffusivity of hydrogen gas through the network of large cages.
