1
|
Walsh MR. Comparing brute force to transition path sampling for gas hydrate nucleation with a flat interface: comments on time reversal symmetry. Phys Chem Chem Phys 2024; 26:5762-5772. [PMID: 38214888 DOI: 10.1039/d3cp05059a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2024]
Abstract
Fluid to solid nucleation is often investigated with the rare event method transition path sampling (TPS). I claim that the inherent irreversibility of solid nucleation, even at stationary conditions, calls into question TPS's applicability for determining solid nucleation mechanisms, especially for pre-critical behavior. Even when applied to a phenomenon which displays time reversal asymmetry like solid nucleation, TPS is a good means of exploring phase space and giving trends in post-critical structure, and its ability to facilitate nucleation rate and free energy calculations remains outstanding. Forward-only splitting and ratcheting methods such as forward flux sampling are more attractive for understanding nucleation mechanisms as they do not require time reversal symmetry, but at low driving forces may suffer from the same limitations as brute force: they may never make it to the first ratchet. Here I briefly summarize the TPS method and gas hydrate nucleation simulation literature, focusing on topics within both to facilitate a comparison of brute force hydrate nucleation to transition path sampling of hydrate nucleation. Perhaps anecdotally, the brute force technique results in more crystalline trajectories despite having higher driving forces than TPS. I maintain this difference is because of the inherent irreversibility of hydrate nucleation, meaning its pre-critical behavior cannot accurately be determined by the melting trajectories that comprise approximately half of the configurations in TPS's path ensemble.
Collapse
|
2
|
Povey MJ, Ettelaie R, Lewtas K, Price A, Lai X, Sheng F. "Sounding" out crystal nuclei-A mathematical-physical and experimental investigation. J Chem Phys 2023; 158:2887908. [PMID: 37129304 DOI: 10.1063/5.0139811] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2022] [Accepted: 04/12/2023] [Indexed: 05/03/2023] Open
Abstract
We outline techniques for the control and measurement of the nucleation of crystalline materials. Small angle x-ray scattering/wide angle x-ray scattering x-ray diffraction measurements are presented that demonstrate the impact of low power, continuous, non-cavitational ultrasound on the nucleation and crystallization of a wax-n-eicosane dissolved in a heptane/toluene solvent. A mathematical-physical approach based on the rectification of heat and mass transport by such a low power oscillating pressure field is outlined, and it is suggested that this approach be combined with dissipative particle dynamics computational modeling to develop a predictive method capable of modeling the impact of low power oscillating pressure fields (acoustics and ultrasonics) on a wide range of nucleating systems. Combining the ultrasound pitch and catch speed of sound measurements with low power harmonically oscillating pressure fields to monitor and control nucleation presents the prospect of entirely new industrially significant methods of process control in crystallization. It also offers new insights into nucleation processes in general. However, for the acoustic control technique to be widely applied , further theoretical and modeling work will be necessary since, at present, we are unable to predict the precise effect of low power ultrasound in any given situation.
Collapse
Affiliation(s)
- Megan J Povey
- School of Food Science and Nutrition, The University of Leeds, Leeds LS2 9JT, United Kingdom
| | - Rammile Ettelaie
- School of Food Science and Nutrition, The University of Leeds, Leeds LS2 9JT, United Kingdom
| | - Ken Lewtas
- Lewtas Science and Technologies, Ltd., 264 Banbury Road, Oxford OX2 7DY, United Kingdom
| | - Andy Price
- Lewtas Science and Technologies, Ltd., 264 Banbury Road, Oxford OX2 7DY, United Kingdom
| | - Xiaojun Lai
- School of Chemical and Process Engineering, The University of Leeds, Leeds LS2 9JT, United Kingdom
| | - Fei Sheng
- School of Chemical and Process Engineering, The University of Leeds, Leeds LS2 9JT, United Kingdom
| |
Collapse
|
3
|
Zeng SY, Hsu CH, Wu TM. Bond Orientational Order Parameters for Classifying Solid-like Clusters in a Lennard-Jones System near Liquid-Solid Transition and at Solid States. J Phys Chem A 2022; 126:2018-2030. [PMID: 35297626 DOI: 10.1021/acs.jpca.1c09527] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In this paper, we introduced an order parameter, named the local structure similarity (LSS), to measure the resemblance of a cluster structure in a liquid with respect to a perfect crystal. The LSS is based on a dot product of two bond orientational order complex vectors, with one vector associated with a particle in a liquid and the other vector with a particle in a crystal. The calculation of the LSS should scan the entire space of the Euler angles determined by the two coordinate frames describing individually the liquid and the crystal. The effectiveness of the LSS was examined by solid-like clusters in a Lennard-Jones (LJ) system near its liquid-solid phase transition and at solid states below its melting point, where the thermodynamic states of the LJ system were obtained by simulation annealing. The LSS measure was utilized to scrutinize the fcc-like, hcp-like, and bcc-like clusters classified by criteria based on W4 and W6 order parameters. As indicated by our results, the two ways of classification are consistent for fcc-like and hcp-like clusters, which are in a close resemblance to their crystalline counterparts. However, the classification with positive W6 for bcc-like clusters is inconsistent with the results of the LSS measure, which was confirmed by clusters in a LJ system confined between two parallel slabs of particles in the bcc structure arrangement.
Collapse
Affiliation(s)
- Sheng-Yuan Zeng
- Institute of Physics, National Yang Ming Chiao Tung University, Hsinchu 30010, Taiwan, R.O.C
| | - Chih-Hao Hsu
- Institute of Physics, National Yang Ming Chiao Tung University, Hsinchu 30010, Taiwan, R.O.C
| | - Ten-Ming Wu
- Institute of Physics, National Yang Ming Chiao Tung University, Hsinchu 30010, Taiwan, R.O.C
| |
Collapse
|
4
|
Blow KE, Quigley D, Sosso GC. The seven deadly sins: When computing crystal nucleation rates, the devil is in the details. J Chem Phys 2021; 155:040901. [PMID: 34340373 DOI: 10.1063/5.0055248] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
The formation of crystals has proven to be one of the most challenging phase transformations to quantitatively model-let alone to actually understand-be it by means of the latest experimental technique or the full arsenal of enhanced sampling approaches at our disposal. One of the most crucial quantities involved with the crystallization process is the nucleation rate, a single elusive number that is supposed to quantify the average probability for a nucleus of critical size to occur within a certain volume and time span. A substantial amount of effort has been devoted to attempt a connection between the crystal nucleation rates computed by means of atomistic simulations and their experimentally measured counterparts. Sadly, this endeavor almost invariably fails to some extent, with the venerable classical nucleation theory typically blamed as the main culprit. Here, we review some of the recent advances in the field, focusing on a number of perhaps more subtle details that are sometimes overlooked when computing nucleation rates. We believe it is important for the community to be aware of the full impact of aspects, such as finite size effects and slow dynamics, that often introduce inconspicuous and yet non-negligible sources of uncertainty into our simulations. In fact, it is key to obtain robust and reproducible trends to be leveraged so as to shed new light on the kinetics of a process, that of crystal nucleation, which is involved into countless practical applications, from the formulation of pharmaceutical drugs to the manufacturing of nano-electronic devices.
Collapse
Affiliation(s)
- Katarina E Blow
- Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - David Quigley
- Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - Gabriele C Sosso
- Department of Chemistry, University of Warwick, Coventry CV4 7AL, United Kingdom
| |
Collapse
|
5
|
Menon S, Díaz Leines G, Drautz R, Rogal J. Role of pre-ordered liquid in the selection mechanism of crystal polymorphs during nucleation. J Chem Phys 2020; 153:104508. [DOI: 10.1063/5.0017575] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Sarath Menon
- Interdisciplinary Centre for Advanced Materials Simulation, Ruhr-Universität Bochum, 44801 Bochum, Germany
| | - Grisell Díaz Leines
- Interdisciplinary Centre for Advanced Materials Simulation, Ruhr-Universität Bochum, 44801 Bochum, Germany
| | - Ralf Drautz
- Interdisciplinary Centre for Advanced Materials Simulation, Ruhr-Universität Bochum, 44801 Bochum, Germany
| | - Jutta Rogal
- Interdisciplinary Centre for Advanced Materials Simulation, Ruhr-Universität Bochum, 44801 Bochum, Germany
| |
Collapse
|
6
|
Liang Y, Díaz Leines G, Drautz R, Rogal J. Identification of a multi-dimensional reaction coordinate for crystal nucleation in Ni3Al. J Chem Phys 2020; 152:224504. [DOI: 10.1063/5.0010074] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Affiliation(s)
- Yanyan Liang
- Interdisciplinary Centre for Advanced Materials Simulation, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Grisell Díaz Leines
- Interdisciplinary Centre for Advanced Materials Simulation, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Ralf Drautz
- Interdisciplinary Centre for Advanced Materials Simulation, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Jutta Rogal
- Interdisciplinary Centre for Advanced Materials Simulation, Ruhr-Universität Bochum, 44780 Bochum, Germany
| |
Collapse
|
7
|
Kuhnhold A, Meyer H, Amati G, Pelagejcev P, Schilling T. Derivation of an exact, nonequilibrium framework for nucleation: Nucleation is a priori neither diffusive nor Markovian. Phys Rev E 2019; 100:052140. [PMID: 31869953 DOI: 10.1103/physreve.100.052140] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2019] [Indexed: 06/10/2023]
Abstract
We discuss the structure of the equation of motion that governs nucleation processes at first order phase transitions. From the underlying microscopic dynamics of a nucleating system, we derive by means of a nonequilibrium projection operator formalism the equation of motion for the size distribution of the nuclei. The equation is exact, i.e., the derivation does not contain approximations. To assess the impact of memory, we express the equation of motion in a form that allows for direct comparison to the Markovian limit. As a numerical test, we have simulated crystal nucleation from a supersaturated melt of particles interacting via a Lennard-Jones potential. The simulation data show effects of non-Markovian dynamics.
Collapse
Affiliation(s)
- Anja Kuhnhold
- Physikalisches Institut, Albert-Ludwigs-Universität, 79104 Freiburg, Germany
| | - Hugues Meyer
- Physikalisches Institut, Albert-Ludwigs-Universität, 79104 Freiburg, Germany
- Research Unit in Engineering Science, Université du Luxembourg, L-4364 Esch-sur-Alzette, Luxembourg
| | - Graziano Amati
- Physikalisches Institut, Albert-Ludwigs-Universität, 79104 Freiburg, Germany
| | - Philipp Pelagejcev
- Physikalisches Institut, Albert-Ludwigs-Universität, 79104 Freiburg, Germany
| | - Tanja Schilling
- Physikalisches Institut, Albert-Ludwigs-Universität, 79104 Freiburg, Germany
| |
Collapse
|
8
|
DeFever RS, Targonski C, Hall SW, Smith MC, Sarupria S. A generalized deep learning approach for local structure identification in molecular simulations. Chem Sci 2019; 10:7503-7515. [PMID: 31768235 PMCID: PMC6839808 DOI: 10.1039/c9sc02097g] [Citation(s) in RCA: 34] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2019] [Accepted: 07/04/2019] [Indexed: 12/20/2022] Open
Abstract
Identifying local structure in molecular simulations is of utmost importance. The most common existing approach to identify local structure is to calculate some geometrical quantity referred to as an order parameter. In simple cases order parameters are physically intuitive and trivial to develop (e.g., ion-pair distance), however in most cases, order parameter development becomes a much more difficult endeavor (e.g., crystal structure identification). Using ideas from computer vision, we adapt a specific type of neural network called a PointNet to identify local structural environments in molecular simulations. A primary challenge in applying machine learning techniques to simulation is selecting the appropriate input features. This challenge is system-specific and requires significant human input and intuition. In contrast, our approach is a generic framework that requires no system-specific feature engineering and operates on the raw output of the simulations, i.e., atomic positions. We demonstrate the method on crystal structure identification in Lennard-Jones (four different phases), water (eight different phases), and mesophase (six different phases) systems. The method achieves as high as 99.5% accuracy in crystal structure identification. The method is applicable to heterogeneous nucleation and it can even predict the crystal phases of atoms near external interfaces. We demonstrate the versatility of our approach by using our method to identify surface hydrophobicity based solely upon positions and orientations of surrounding water molecules. Our results suggest the approach will be broadly applicable to many types of local structure in simulations.
Collapse
Affiliation(s)
- Ryan S DeFever
- Department of Chemical & Biomolecular Engineering , Clemson University , Clemson , SC 29634 , USA
| | - Colin Targonski
- Department of Electrical & Computer Engineering , Clemson University , Clemson , SC 29634 , USA .
| | - Steven W Hall
- Department of Chemical & Biomolecular Engineering , Clemson University , Clemson , SC 29634 , USA
| | - Melissa C Smith
- Department of Electrical & Computer Engineering , Clemson University , Clemson , SC 29634 , USA .
| | - Sapna Sarupria
- Department of Chemical & Biomolecular Engineering , Clemson University , Clemson , SC 29634 , USA
| |
Collapse
|
9
|
Paul S, Nair NN, Vashisth H. Phase space and collective variable based simulation methods for studies of rare events. MOLECULAR SIMULATION 2019. [DOI: 10.1080/08927022.2019.1634268] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- Sanjib Paul
- Department of Chemical Engineering, University of New Hampshire, Durham, NH, USA
| | - Nisanth N. Nair
- Department of Chemistry, Indian Institute of Technology, Kanpur, India
| | - Harish Vashisth
- Department of Chemical Engineering, University of New Hampshire, Durham, NH, USA
| |
Collapse
|
10
|
Díaz Leines G, Rogal J. Maximum Likelihood Analysis of Reaction Coordinates during Solidification in Ni. J Phys Chem B 2018; 122:10934-10942. [DOI: 10.1021/acs.jpcb.8b08718] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Grisell Díaz Leines
- Interdisciplinary Centre for Advanced Materials Simulation, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Jutta Rogal
- Interdisciplinary Centre for Advanced Materials Simulation, Ruhr-Universität Bochum, 44780 Bochum, Germany
| |
Collapse
|
11
|
Bi Y, Porras A, Li T. Free energy landscape and molecular pathways of gas hydrate nucleation. J Chem Phys 2018; 145:211909. [PMID: 28799352 DOI: 10.1063/1.4961241] [Citation(s) in RCA: 37] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Abstract
Despite the significance of gas hydrates in diverse areas, a quantitative knowledge of hydrate formation at a molecular level is missing. The impediment to acquiring this understanding is primarily attributed to the stochastic nature and ultra-fine scales of nucleation events, posing a great challenge for both experiment and simulation to explore hydrate nucleation. Here we employ advanced molecular simulation methods, including forward flux sampling (FFS), pB histogram analysis, and backward flux sampling, to overcome the limit of direct molecular simulation for exploring both the free energy landscape and molecular pathways of hydrate nucleation. First we test the half-cage order parameter (H-COP) which we developed for driving FFS, through conducting the pB histogram analysis. Our results indeed show that H-COP describes well the reaction coordinates of hydrate nucleation. Through the verified order parameter, we then directly compute the free energy landscape for hydrate nucleation by combining both forward and backward flux sampling. The calculated stationary distribution density, which is obtained independently of nucleation theory, is found to fit well against the classical nucleation theory (CNT). Subsequent analysis of the obtained large ensemble of hydrate nucleation trajectories show that although on average, hydrate formation is facilitated by a two-step like mechanism involving a gradual transition from an amorphous to a crystalline structure, there also exist nucleation pathways where hydrate crystallizes directly, without going through the amorphous stage. The CNT-like free energy profile and the structural diversity suggest the existence of multiple active transition pathways for hydrate nucleation, and possibly also imply the near degeneracy in their free energy profiles among different pathways. Our results thus bring a new perspective to the long standing question of how hydrates crystallize.
Collapse
Affiliation(s)
- Yuanfei Bi
- Department of Civil and Environmental Engineering, George Washington University, Washington DC 20052, USA
| | - Anna Porras
- Department of Civil and Environmental Engineering, George Washington University, Washington DC 20052, USA
| | - Tianshu Li
- Department of Civil and Environmental Engineering, George Washington University, Washington DC 20052, USA
| |
Collapse
|
12
|
Eslami H, Sedaghat P, Müller-Plathe F. Local bond order parameters for accurate determination of crystal structures in two and three dimensions. Phys Chem Chem Phys 2018; 20:27059-27068. [DOI: 10.1039/c8cp05248d] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
Local order parameters for the characterization of liquid and different two- and three-dimensional crystalline structures are presented.
Collapse
Affiliation(s)
- Hossein Eslami
- Department of Chemistry
- College of Sciences
- Persian Gulf University
- Boushehr 75168
- Iran
| | - Parvin Sedaghat
- Department of Chemistry
- College of Sciences
- Persian Gulf University
- Boushehr 75168
- Iran
| | - Florian Müller-Plathe
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie
- Technische Universität Darmstadt
- 64287 Darmstadt
- Germany
| |
Collapse
|
13
|
Cabriolu R, Skjelbred Refsnes KM, Bolhuis PG, van Erp TS. Foundations and latest advances in replica exchange transition interface sampling. J Chem Phys 2017; 147:152722. [DOI: 10.1063/1.4989844] [Citation(s) in RCA: 30] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Raffaela Cabriolu
- Department of Chemistry, Norwegian University of Science and Technology (NTNU), Høgskoleringen 5, 7491 Trondheim, Norway
| | - Kristin M. Skjelbred Refsnes
- Department of Chemistry, Norwegian University of Science and Technology (NTNU), Høgskoleringen 5, 7491 Trondheim, Norway
| | - Peter G. Bolhuis
- Van ’t Hoff Institute for Molecular Sciences, University of Amsterdam, P.O. Box 94157, 1090 GD Amsterdam, The Netherlands
| | - Titus S. van Erp
- Department of Chemistry, Norwegian University of Science and Technology (NTNU), Høgskoleringen 5, 7491 Trondheim, Norway
| |
Collapse
|
14
|
Díaz Leines G, Drautz R, Rogal J. Atomistic insight into the non-classical nucleation mechanism during solidification in Ni. J Chem Phys 2017; 146:154702. [DOI: 10.1063/1.4980082] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Grisell Díaz Leines
- Interdisciplinary Centre for Advanced Materials Simulation, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Ralf Drautz
- Interdisciplinary Centre for Advanced Materials Simulation, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Jutta Rogal
- Interdisciplinary Centre for Advanced Materials Simulation, Ruhr-Universität Bochum, 44780 Bochum, Germany
| |
Collapse
|
15
|
Eslami H, Khanjari N, Müller-Plathe F. A Local Order Parameter-Based Method for Simulation of Free Energy Barriers in Crystal Nucleation. J Chem Theory Comput 2017; 13:1307-1316. [DOI: 10.1021/acs.jctc.6b01034] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Hossein Eslami
- Eduard-Zintl-Institut
für Anorganische und Physikalische Chemie, Technische Universität Darmstadt, Alarich-Weiss-Str. 8, 64287 Darmstadt, Germany
- Department
of Chemistry, College of Sciences, Persian Gulf University, Boushehr 75168, Iran
| | - Neda Khanjari
- Department
of Chemistry, College of Sciences, Persian Gulf University, Boushehr 75168, Iran
| | - Florian Müller-Plathe
- Eduard-Zintl-Institut
für Anorganische und Physikalische Chemie, Technische Universität Darmstadt, Alarich-Weiss-Str. 8, 64287 Darmstadt, Germany
- Department
of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| |
Collapse
|
16
|
Dittmar HR, Kusalik PG. Ordering effects of conjugate thermal fields in simulations of molecular liquids: Carbon dioxide and water. J Chem Phys 2016; 145:134504. [DOI: 10.1063/1.4963767] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
Affiliation(s)
- Harro R. Dittmar
- Department of Chemistry, University of Calgary, Calgary, Alberta T2N 1N4, Canada
| | - Peter G. Kusalik
- Department of Chemistry, University of Calgary, Calgary, Alberta T2N 1N4, Canada
| |
Collapse
|
17
|
Jungblut S, Dellago C. Pathways to self-organization: Crystallization via nucleation and growth. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2016; 39:77. [PMID: 27498980 DOI: 10.1140/epje/i2016-16077-6] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/13/2016] [Accepted: 07/13/2016] [Indexed: 06/06/2023]
Abstract
Crystallization, a prototypical self-organization process during which a disordered state spontaneously transforms into a crystal characterized by a regular arrangement of its building blocks, usually proceeds by nucleation and growth. In the initial stages of the transformation, a localized nucleus of the new phase forms in the old one due to a random fluctuation. Most of these nuclei disappear after a short time, but rarely a crystalline embryo may reach a critical size after which further growth becomes thermodynamically favorable and the entire system is converted into the new phase. In this article, we will discuss several theoretical concepts and computational methods to study crystallization. More specifically, we will address the rare event problem arising in the simulation of nucleation processes and explain how to calculate nucleation rates accurately. Particular attention is directed towards discussing statistical tools to analyze crystallization trajectories and identify the transition mechanism.
Collapse
Affiliation(s)
- S Jungblut
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, 1090, Wien, Austria
| | - C Dellago
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, 1090, Wien, Austria.
| |
Collapse
|
18
|
Abstract
Pathogenic protein fibrils have been shown in vitro to have nucleation-dependent kinetics despite the fact that one-dimensional structures do not have the size-dependent surface energy responsible for the lag time in classical theory. We present a theory showing that the conformational entropy of the peptide chains creates a free-energy barrier that is analogous to the translational entropy barrier in higher dimensions. We find that the dynamics of polymer rearrangement make it very unlikely for nucleation to succeed along the lowest free-energy trajectory, meaning that most of the nucleation flux avoids the free-energy saddle point. We use these results to construct a three-dimensional model for amyloid nucleation that accounts for conformational entropy, backbone H bonds, and side-chain interactions to compute nucleation rates as a function of concentration.
Collapse
Affiliation(s)
- Lingyun Zhang
- Department of Physics, Kansas State University, Manhattan, KS 66506, USA
- Beijing National Laboratory for Condensed Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
| | - Jeremy D. Schmit
- Department of Physics, Kansas State University, Manhattan, KS 66506, USA
| |
Collapse
|
19
|
Paul S, Taraphder S. Determination of the Reaction Coordinate for a Key Conformational Fluctuation in Human Carbonic Anhydrase II. J Phys Chem B 2015; 119:11403-15. [DOI: 10.1021/acs.jpcb.5b03655] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Affiliation(s)
- Sanjib Paul
- Department of Chemistry, Indian Institute of Technology, Kharagpur 721302, India
| | - Srabani Taraphder
- Department of Chemistry, Indian Institute of Technology, Kharagpur 721302, India
| |
Collapse
|
20
|
Mithen JP, Callison AJ, Sear RP. Nucleation of crystals that are mixed composites of all three polymorphs in the Gaussian core model. J Chem Phys 2015; 142:224505. [DOI: 10.1063/1.4922321] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023] Open
Affiliation(s)
- J. P. Mithen
- Department of Physics, University of Surrey, Guildford GU2 7XH, United Kingdom
| | - A. J. Callison
- Department of Physics, University of Surrey, Guildford GU2 7XH, United Kingdom
| | - R. P. Sear
- Department of Physics, University of Surrey, Guildford GU2 7XH, United Kingdom
| |
Collapse
|
21
|
Leitold C, Lechner W, Dellago C. A string reaction coordinate for the folding of a polymer chain. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2015; 27:194126. [PMID: 25923377 DOI: 10.1088/0953-8984/27/19/194126] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
We investigate the crystallization mechanism of a single, flexible homopolymer chain with short range attractions. For a sufficiently narrow attractive well, the system undergoes a first-order like freezing transition from an expanded disordered coil to a compact crystalline state. Based on a maximum likelihood analysis of committor values computed for configurations obtained by Wang-Landau sampling, we construct a non-linear string reaction coordinate for the coil-to-crystal transition. In contrast to a linear reaction coordinate, the string reaction coordinate captures the effect of different degrees of freedom controlling different stages of the transition. Our analysis indicates that a combination of the energy and the global crystallinity parameter Q6 provide the most accurate measure for the progress of the transition. While the crystallinity paramter Q6 is most relevant in the initial stages of the crystallization, the later stages are dominated by a decrease in the potential energy.
Collapse
Affiliation(s)
- Christian Leitold
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, 1090 Vienna, Austria
| | | | | |
Collapse
|
22
|
Barnes BC, Knott BC, Beckham GT, Wu DT, Sum AK. Reaction coordinate of incipient methane clathrate hydrate nucleation. J Phys Chem B 2014; 118:13236-43. [PMID: 25347748 DOI: 10.1021/jp507959q] [Citation(s) in RCA: 64] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
Nucleation from solution is a ubiquitous phenomenon with relevance to myriad scientific disciplines, including pharmaceuticals, biomineralization, and disease. One prominent example is the nucleation of clathrate hydrates, multicomponent crystalline inclusion compounds relevant to the energy industry where they block pipelines and also constitute a potential vast energy resource. Despite their importance, the molecular mechanism of incipient hydrate formation remains unknown. Herein, we employ advanced molecular simulation tools (pB histogram, equilibrium path sampling) to provide a statistical-mechanical basis for extracting physical insight into the molecular steps by which clathrates form. Through testing the Mutually Coordinated Guest (MCG) order parameter, we demonstrate that both guest (methane) and host (water) structuring are crucial to accurately describe the nucleation of hydrates and determine a critical nucleus size of MCG-1 = 16 at 255 K and 500 bar. Equipped with a validated (and novel) reaction coordinate, subsequent equilibrium path sampling simulations yield the free energy barrier and nucleation rate. The resulting quantitative nucleation process is described by the MCG clustering mechanism. This constitutes a significant advance in the field of hydrates research, as the fitness of a molecular descriptor has never been statistically verified. More broadly, this work has significance to a wide range of multicomponent nucleation contexts wherein the formation mechanism depends on contributions from both solute and solvent.
Collapse
Affiliation(s)
- Brian C Barnes
- Center for Hydrate Research, Chemical & Biological Engineering Department, Colorado School of Mines , Golden, Colorado 80401, United States
| | | | | | | | | |
Collapse
|
23
|
Wang GC, Wang Q, Li SL, Ai XG, Fan CG. Evidence of multi-step nucleation leading to various crystallization pathways from an Fe-O-Al melt. Sci Rep 2014; 4:5082. [PMID: 24866413 PMCID: PMC4035573 DOI: 10.1038/srep05082] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/16/2014] [Accepted: 05/07/2014] [Indexed: 01/26/2023] Open
Abstract
The crystallization process from a solution begins with nucleation, which determines the structure and size of the resulting crystals. Further understanding of multi-pathway crystallizations from solution through two-step nucleation mechanisms is needed. This study uses density functional theory to probe the thermodynamic properties of alumina clusters at high temperature and reveals the thermodynamic relationship between these clusters and the saturation levels of dissolved oxygen and aluminum in an Fe–O–Al melt. Based on the thermodynamics of cluster formation and the experimental evidence for both excess oxygen in the Fe-O-Al melt and for alumina with a polycrystalline structure in solidified iron, we demonstrate that the appearance of various types of clusters that depends on the saturation ratio determines the nucleation steps that lead to the various crystallization pathways. Such mechanisms may also be important in nucleation and crystallization from solution.
Collapse
Affiliation(s)
- G C Wang
- 1] Key Laboratory of Chemical Metallurgy Engineering, Liaoning Province, University of Science and Technology Liaoning, Anshan, Liaoning, 114051, China [2] School of Materials and Metallurgy, University of Science and Technology Liaoning, Anshan, 114051, China [3] Jiangxi University of Science and Technology, Ganzhou, 341000, China
| | - Q Wang
- 1] Key Laboratory of Chemical Metallurgy Engineering, Liaoning Province, University of Science and Technology Liaoning, Anshan, Liaoning, 114051, China [2] School of Materials and Metallurgy, University of Science and Technology Liaoning, Anshan, 114051, China
| | - S L Li
- 1] Key Laboratory of Chemical Metallurgy Engineering, Liaoning Province, University of Science and Technology Liaoning, Anshan, Liaoning, 114051, China [2] School of Materials and Metallurgy, University of Science and Technology Liaoning, Anshan, 114051, China
| | - X G Ai
- 1] Key Laboratory of Chemical Metallurgy Engineering, Liaoning Province, University of Science and Technology Liaoning, Anshan, Liaoning, 114051, China [2] School of Materials and Metallurgy, University of Science and Technology Liaoning, Anshan, 114051, China
| | - C G Fan
- Jiangxi University of Science and Technology, Ganzhou, 341000, China
| |
Collapse
|
24
|
Dittmar H, Kusalik PG. Driving ordering processes in molecular-dynamics simulations. PHYSICAL REVIEW LETTERS 2014; 112:195701. [PMID: 24877946 DOI: 10.1103/physrevlett.112.195701] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/06/2014] [Indexed: 05/28/2023]
Abstract
Self-organized criticality describes the emergence of complexity in dynamical nonequilibrium systems. In this paper we introduce a unique approach whereby a driven energy conversion is utilized as a sampling bias for ordered arrangements in molecular dynamics simulations of atomic and molecular fluids. This approach gives rise to dramatically accelerated nucleation rates, by as much as 30 orders of magnitude, without the need of predefined order parameters, which commonly employed rare-event sampling methods rely on. The measured heat fluxes suggest how the approach can be generalized.
Collapse
Affiliation(s)
- Harro Dittmar
- Department of Chemistry, University of Calgary, 2500 University drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Peter G Kusalik
- Department of Chemistry, University of Calgary, 2500 University drive NW, Calgary, Alberta T2N 1N4, Canada
| |
Collapse
|
25
|
Agarwal V, Peters B. Solute Precipitate Nucleation: A Review of Theory and Simulation Advances. ADVANCES IN CHEMICAL PHYSICS 2014. [DOI: 10.1002/9781118755815.ch03] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
26
|
Li W, Ma A. Recent developments in methods for identifying reaction coordinates. MOLECULAR SIMULATION 2014; 40:784-793. [PMID: 25197161 PMCID: PMC4152980 DOI: 10.1080/08927022.2014.907898] [Citation(s) in RCA: 44] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Abstract
In the study of rare events in complex systems with many degrees of freedom, a key element is to identify the reaction coordinates of a given process. Over recent years, a number of methods and protocols have been developed to extract the reaction coordinates based on limited information from molecular dynamics simulations. In this review, we provide a brief survey over a number of major methods developed in the past decade, some of which are discussed in greater detail, to provide an overview of the problems that are partially solved and challenges that still remain. A particular emphasis has been placed on methods for identifying reaction coordinates that are related to the committor.
Collapse
Affiliation(s)
- Wenjin Li
- Department of Bioengineering, The University of Illinois at Chicago, 851 South Morgan Street, Chicago, IL 60607, USA
| | - Ao Ma
- Department of Bioengineering, The University of Illinois at Chicago, 851 South Morgan Street, Chicago, IL 60607, USA
| |
Collapse
|
27
|
Cardenas AE, Elber R. Computational study of peptide permeation through membrane: Searching for hidden slow variables. Mol Phys 2013; 111:3565-3578. [PMID: 26203198 PMCID: PMC4507298 DOI: 10.1080/00268976.2013.842010] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Atomically detailed molecular dynamics trajectories in conjunction with Milestoning are used to analyze the different contributions of coarse variables to the permeation process of a small peptide (N-acetyl-L-tryptophanamide, NATA) through a 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) membrane. The peptide reverses its overall orientation as it permeates through the biological bilayer. The large change in orientation is investigated explicitly but is shown to impact the free energy landscape and permeation time only moderately. Nevertheless, a significant difference in permeation properties of the two halves of the membrane suggests the presence of other hidden slow variables. We speculate, based on calculation of the potential of mean force, that a conformational transition of NATA makes significant contribution to these differences. Other candidates for hidden slow variables may include water permeation and collective motions of phospholipids.
Collapse
Affiliation(s)
- Alfredo E. Cardenas
- Institute for Computational Engineering and Sciences, University of Texas at Austin, Austin TX 78712, USA
| | - Ron Elber
- Institute for Computational Engineering and Sciences, University of Texas at Austin, Austin TX 78712, USA
- Department of Chemistry and Biochemistry, University of Texas at Austin, Austin TX 78712, USA
| |
Collapse
|
28
|
Affiliation(s)
- Allan S. Myerson
- Department of Chemical Engineering, Massachusetts Institute of Technology, 77 Massachusetts Avenue, E19-502, Cambridge MA 02130, USA
| | - Bernhardt L. Trout
- Department of Chemical Engineering, Massachusetts Institute of Technology, 77 Massachusetts Avenue, E19-502, Cambridge MA 02130, USA
| |
Collapse
|
29
|
Du WN, Bolhuis PG. Adaptive single replica multiple state transition interface sampling. J Chem Phys 2013; 139:044105. [DOI: 10.1063/1.4813777] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
30
|
Reinhardt A, Doye JPK, Noya EG, Vega C. Local order parameters for use in driving homogeneous ice nucleation with all-atom models of water. J Chem Phys 2012. [DOI: 10.1063/1.4766362] [Citation(s) in RCA: 81] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023] Open
|
31
|
Jas GS, Hegefeld W, Májek P, Kuczera K, Elber R. Experiments and comprehensive simulations of the formation of a helical turn. J Phys Chem B 2012; 116:6598-610. [PMID: 22335541 PMCID: PMC3361543 DOI: 10.1021/jp211645s] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
We investigate the kinetics and thermodynamics of a helical turn formation in the peptide Ac-WAAAH-NH(2). NMR measurements indicate that this peptide has significant tendency to form a structure of a helical turn, while temperature dependent CD establishes the helix fraction at different temperatures. Molecular dynamics and milestoning simulations agree with experimental observables and suggest an atomically detailed picture for the turn formation. Using a network representation, two alternative mechanisms of folding are identified: (i) a direct co-operative mechanism from the unfolded to the folded state without intermediate formation of hydrogen bonds and (ii) an indirect mechanism with structural intermediates with two residues in a helical conformation. This picture is consistent with kinetic measurements that reveal two experimental time scales of sub-nanosecond and several nanoseconds.
Collapse
Affiliation(s)
- Gouri S. Jas
- Department of Chemistry, Biochemistry, and Institute of Biomedical Studies, Baylor University, Waco, TX 76706
| | - Wendy Hegefeld
- Department of Chemistry, Biochemistry, and Institute of Biomedical Studies, Baylor University, Waco, TX 76706
| | - Peter Májek
- Institute of Computational Engineering and Sciences (ICES), University of Texas at Austin, Austin, TX 78712
| | - Krzysztof Kuczera
- Departments of Chemistry and Molecular Biosciences, The University of Kansas, Lawrence, KS 66045
| | - Ron Elber
- Institute of Computational Engineering and Sciences (ICES), University of Texas at Austin, Austin, TX 78712
- Department of Chemistry and Biochemistry, University of Texas at Austin, Austin, TX 78712
| |
Collapse
|