1
|
Rogal J, Díaz Leines G. Controlling crystallization: what liquid structure and dynamics reveal about crystal nucleation mechanisms. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2023; 381:20220249. [PMID: 37211029 DOI: 10.1098/rsta.2022.0249] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/03/2022] [Accepted: 12/06/2022] [Indexed: 05/23/2023]
Abstract
Over recent years, molecular simulations have provided invaluable insights into the microscopic processes governing the initial stages of crystal nucleation and growth. A key aspect that has been observed in many different systems is the formation of precursors in the supercooled liquid that precedes the emergence of crystalline nuclei. The structural and dynamical properties of these precursors determine to a large extent the nucleation probability as well as the formation of specific polymorphs. This novel microscopic view on nucleation mechanisms has further implications for our understanding of the nucleating ability and polymorph selectivity of nucleating agents, as these appear to be strongly linked to their ability in modifying structural and dynamical characteristics of the supercooled liquid, namely liquid heterogeneity. In this perspective, we highlight recent progress in exploring the connection between liquid heterogeneity and crystallization, including the effects of templates, and the potential impact for controlling crystallization processes. This article is part of a discussion meeting issue 'Supercomputing simulations of advanced materials'.
Collapse
Affiliation(s)
- Jutta Rogal
- Department of Chemistry, New York University, New York, NY 10003, USA
- Fachbereich Physik, Freie Universität Berlin, 14195 Berlin, Germany
| | - Grisell Díaz Leines
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, UK
| |
Collapse
|
2
|
Duan T, Shen Y, Imhoff SD, Yi F, Voyles PM, Perepezko JH. Nucleation kinetics model for primary crystallization in Al-Y-Fe metallic glass. J Chem Phys 2023; 158:064504. [PMID: 36792527 DOI: 10.1063/5.0135730] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/11/2023] Open
Abstract
The high density of aluminum nanocrystals (>1021 m-3) that develop during the primary crystallization in Al-based metallic glasses indicates a high nucleation rate (∼1018 m-3 s-1). Several studies have been advanced to account for the primary crystallization behavior, but none have been developed to completely describe the reaction kinetics. Recently, structural analysis by fluctuation electron microscopy has demonstrated the presence of the Al-like medium range order (MRO) regions as a spatial heterogeneity in as-spun Al88Y7Fe5 metallic glass that is representative for the class of Al-based amorphous alloys that develop Al nanocrystals during primary crystallization. From the structural characterization, an MRO seeded nucleation configuration is established, whereby the Al nanocrystals are catalyzed by the MRO core to decrease the nucleation barrier. The MRO seeded nucleation model and the kinetic data from the delay time (τ) measurement provide a full accounting of the evolution of the Al nanocrystal density (Nv) during the primary crystallization under isothermal annealing treatments. Moreover, the calculated values of the steady state nucleation rates (Jss) predicted by the nucleation model agree with the experimental results. Moreover, the model satisfies constraints on the structural, thermodynamic, and kinetic parameters, such as the critical nucleus size, the interface energy, and the volume-free energy driving force that are essential for a fully self-consistent nucleation kinetics analysis. The nucleation kinetics model can be applied more broadly to materials that are characterized by the presence of spatial heterogeneities.
Collapse
Affiliation(s)
- Tianrui Duan
- Department of Materials Science and Engineering, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| | - Ye Shen
- Department of Materials Science and Engineering, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| | - Seth D Imhoff
- Materials Science and Technology Sigma Division, Los Alamos National Laboratory, P.O. Box 1663, M.S. G770, Los Alamos, New Mexico 87545, USA
| | - Feng Yi
- Nanolab Technologies, Inc., Milpitas, California 95035, USA
| | - Paul M Voyles
- Department of Materials Science and Engineering, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| | - John H Perepezko
- Department of Materials Science and Engineering, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| |
Collapse
|
3
|
Hall SW, Díaz Leines G, Sarupria S, Rogal J. Practical guide to replica exchange transition interface sampling and forward flux sampling. J Chem Phys 2022; 156:200901. [DOI: 10.1063/5.0080053] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Steven W. Hall
- Department of Chemical Engineering and Materials Science, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | - Grisell Díaz Leines
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridgeshire CB2 1EW, United Kingdom
| | - Sapna Sarupria
- Department of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455, USA
- Department of Chemical and Biomolecular Engineering, Clemson University, Clemson, South Carolina 29634, USA
| | - Jutta Rogal
- Department of Chemistry, New York University, New York, New York 10003, USA
- Fachbereich Physik, Freie Universität Berlin, 14195 Berlin, Germany
| |
Collapse
|
4
|
Díaz Leines G, Rogal J. Template-Induced Precursor Formation in Heterogeneous Nucleation: Controlling Polymorph Selection and Nucleation Efficiency. PHYSICAL REVIEW LETTERS 2022; 128:166001. [PMID: 35522521 DOI: 10.1103/physrevlett.128.166001] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/26/2021] [Accepted: 03/29/2022] [Indexed: 06/14/2023]
Abstract
We present an atomistic study of heterogeneous nucleation in Ni employing transition path sampling, which reveals a template precursor-mediated mechanism of crystallization. Most notably, we find that the ability of tiny templates to modify the structural features of the liquid and promote the formation of precursor regions with enhanced bond-orientational order is key to determining their nucleation efficiency and the polymorphs that crystallize. Our results reveal an intrinsic link between structural liquid heterogeneity and the nucleating ability of templates, which significantly advances our understanding toward the control of nucleation efficiency and polymorph selection.
Collapse
Affiliation(s)
- Grisell Díaz Leines
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridgeshire CB2 1EW, United Kingdom
| | - Jutta Rogal
- Department of Chemistry, New York University, New York, New York 10003, USA and Fachbereich Physik, Freie Universität Berlin, 14195 Berlin, Germany
| |
Collapse
|
5
|
Bal KM. Nucleation rates from small scale atomistic simulations and transition state theory. J Chem Phys 2021; 155:144111. [PMID: 34654300 DOI: 10.1063/5.0063398] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
Abstract
The evaluation of nucleation rates from molecular dynamics trajectories is hampered by the slow nucleation time scale and impact of finite size effects. Here, we show that accurate nucleation rates can be obtained in a very general fashion relying only on the free energy barrier, transition state theory, and a simple dynamical correction for diffusive recrossing. In this setup, the time scale problem is overcome by using enhanced sampling methods, in casu metadynamics, whereas the impact of finite size effects can be naturally circumvented by reconstructing the free energy surface from an appropriate ensemble. Approximations from classical nucleation theory are avoided. We demonstrate the accuracy of the approach by calculating macroscopic rates of droplet nucleation from argon vapor, spanning 16 orders of magnitude and in excellent agreement with literature results, all from simulations of very small (512 atom) systems.
Collapse
Affiliation(s)
- Kristof M Bal
- Department of Chemistry and NANOlab Center of Excellence, University of Antwerp, Universiteitsplein 1, 2610 Antwerp, Belgium
| |
Collapse
|
6
|
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
|
7
|
Arjun A, Bolhuis PG. Homogenous nucleation rate of CO 2 hydrates using transition interface sampling. J Chem Phys 2021; 154:164507. [PMID: 33940852 DOI: 10.1063/5.0044883] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Carbon dioxide and water can form solid clathrate structures in which water cages encapsulate the gas molecules. Such hydrates have sparked much interest due to their possible application in CO2 sequestration. How the solid structure forms exactly from the liquid phase via a homogenous nucleation process is still poorly understood. This nucleation event is rare on the molecular timescale even under moderate undercooling or supersaturation conditions because of the large free energy barrier toward crystallization, rendering a brute force simulation of hydrate nucleation unfeasible for moderate undercooling or supersaturation. Here, we perform transition interface sampling simulations to quantify the homogenous nucleation rate for CO2 hydrate formation using accurate atomistic force fields at 500 bars for three different temperatures between 260 and 273 K. Collecting more than 100 000 pathways comprising roughly two milliseconds of simulation time, we computed a nucleation rate in the amorphous phase of ∼1021 nuclei s-1 cm-3 for a temperature of 260 K and a rate of ∼1012 nuclei s-1 cm-3 for a temperature of 265 K. For a temperature of 273 K, we find that the hydrate forms an sI crystalline phase with a rate of order of ∼101 nuclei s-1 cm-3. We compare these rates to classical nucleation theory estimates as well as experiments, and to nucleation rate estimates for methane hydrates and discuss possible causes of the observed differences. Our findings shed light on the kinetics of this important clathrate and should assist in future hydrate formation investigation.
Collapse
Affiliation(s)
- A Arjun
- van 't Hoff Institute for Molecular Sciences, University of Amsterdam, P.O. Box 94157, 1090 GD Amsterdam, The Netherlands
| | - Peter G Bolhuis
- van 't Hoff Institute for Molecular Sciences, University of Amsterdam, P.O. Box 94157, 1090 GD Amsterdam, The Netherlands
| |
Collapse
|
8
|
Abstract
The freezing of water into ice is one of the most important processes in the physical sciences. However, it is still not understood at the molecular level. In particular, the crystallization of cubic ice ([Formula: see text])-rather than the traditional hexagonal polytype ([Formula: see text])-has become an increasingly debated topic. Although evidence for [Formula: see text] is thought to date back almost 400 y, it is only in the last year that pure [Formula: see text] has been made in the laboratory, and these processes involved high-pressure ice phases. Since this demonstrates that pure [Formula: see text] can form, the question naturally arises if [Formula: see text] can be made from liquid water. With this in mind, we have performed a high-throughput computational screening study involving molecular dynamics simulations of nucleation on over 1,100 model substrates. From these simulations, we find that 1) many different substrates can promote the formation of pristine [Formula: see text]; 2) [Formula: see text] can be selectively nucleated for even the mildest supercooling; 3) the water contact layer's resemblance to a face of ice is the key factor determining the polytype selectivity and nucleation temperature, independent of which polytype is promoted; and 4) substrate lattice match to ice is not indicative of the polytype obtained. Through this study, we have deepened understanding of the interplay of heterogeneous nucleation and ice I polytypism and suggest routes to [Formula: see text] More broadly, the substrate design methodology presented here combined with the insight gained can be used to understand and control polymorphism and stacking disorder in materials in general.
Collapse
|