1
|
Lin M, Xiong Z, Cao H. Bridging classical nucleation theory and molecular dynamics simulation for homogeneous ice nucleation. J Chem Phys 2024; 161:084504. [PMID: 39206829 DOI: 10.1063/5.0216645] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2024] [Accepted: 08/14/2024] [Indexed: 09/04/2024] Open
Abstract
Water freezing, initiated by ice nucleation, occurs widely in nature, ranging from cellular to global phenomena. Ice nucleation has been experimentally proven to require the formation of a critical ice nucleus, consistent with classical nucleation theory (CNT). However, the accuracy of CNT quantitative predictions of critical cluster sizes and nucleation rates has never been verified experimentally. In this study, we circumvent this difficulty by using molecular dynamics (MD) simulation. The physical properties of water/ice for CNT predictions, including density, chemical potential difference, and diffusion coefficient, are independently obtained using MD simulation, whereas the calculation of interfacial free energy is based on thermodynamic assumptions of CNT, including capillarity approximation among others. The CNT predictions are compared to the MD evaluations of brute-force simulations and forward flux sampling methods. We find that the CNT and MD predicted critical cluster sizes are consistent, and the CNT predicted nucleation rates are higher than the MD predicted values within three orders of magnitude. We also find that the ice crystallized from supercooled water is stacking-disordered ice with a stacking of cubic and hexagonal ices in four representative types of stacking. The prediction discrepancies in nucleation rate mainly arise from the stacking-disordered ice structure, the asphericity of ice cluster, the uncertainty of ice-water interfacial free energy, and the kinetic attachment rate. Our study establishes a relation between CNT and MD to predict homogeneous ice nucleation.
Collapse
Affiliation(s)
- Min Lin
- Key Laboratory for Thermal Science and Power Engineering of Ministry of Education, Department of Energy and Power Engineering, Tsinghua University, Beijing 100084, People's Republic of China
| | - Zhewen Xiong
- Key Laboratory for Thermal Science and Power Engineering of Ministry of Education, Department of Energy and Power Engineering, Tsinghua University, Beijing 100084, People's Republic of China
| | - Haishan Cao
- Key Laboratory for Thermal Science and Power Engineering of Ministry of Education, Department of Energy and Power Engineering, Tsinghua University, Beijing 100084, People's Republic of China
| |
Collapse
|
2
|
Chew PY, Joseph JA, Collepardo-Guevara R, Reinhardt A. Thermodynamic origins of two-component multiphase condensates of proteins. Chem Sci 2023; 14:1820-1836. [PMID: 36819870 PMCID: PMC9931050 DOI: 10.1039/d2sc05873a] [Citation(s) in RCA: 10] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2022] [Accepted: 01/06/2023] [Indexed: 01/26/2023] Open
Abstract
Intracellular condensates are highly multi-component systems in which complex phase behaviour can ensue, including the formation of architectures comprising multiple immiscible condensed phases. Relying solely on physical intuition to manipulate such condensates is difficult because of the complexity of their composition, and systematically learning the underlying rules experimentally would be extremely costly. We address this challenge by developing a computational approach to design pairs of protein sequences that result in well-separated multilayered condensates and elucidate the molecular origins of these compartments. Our method couples a genetic algorithm to a residue-resolution coarse-grained protein model. We demonstrate that we can design protein partners to form multiphase condensates containing naturally occurring proteins, such as the low-complexity domain of hnRNPA1 and its mutants, and show how homo- and heterotypic interactions must differ between proteins to result in multiphasicity. We also show that in some cases the specific pattern of amino-acid residues plays an important role. Our findings have wide-ranging implications for understanding and controlling the organisation, functions and material properties of biomolecular condensates.
Collapse
Affiliation(s)
- Pin Yu Chew
- Yusuf Hamied Department of Chemistry, University of Cambridge Cambridge CB2 1EW UK
| | - Jerelle A. Joseph
- Yusuf Hamied Department of Chemistry, University of CambridgeCambridgeCB2 1EWUK,Department of Physics, University of CambridgeCambridgeCB3 0HEUK,Department of Genetics, University of CambridgeCambridgeCB2 3EHUK
| | - Rosana Collepardo-Guevara
- Yusuf Hamied Department of Chemistry, University of Cambridge Cambridge CB2 1EW UK .,Department of Physics, University of Cambridge Cambridge CB3 0HE UK.,Department of Genetics, University of Cambridge Cambridge CB2 3EH UK
| | - Aleks Reinhardt
- Yusuf Hamied Department of Chemistry, University of Cambridge Cambridge CB2 1EW UK
| |
Collapse
|
3
|
Chew PY, Reinhardt A. Phase diagrams-Why they matter and how to predict them. J Chem Phys 2023; 158:030902. [PMID: 36681642 DOI: 10.1063/5.0131028] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
Understanding the thermodynamic stability and metastability of materials can help us to, for example, gauge whether crystalline polymorphs in pharmaceutical formulations are likely to be durable. It can also help us to design experimental routes to novel phases with potentially interesting properties. In this Perspective, we provide an overview of how thermodynamic phase behavior can be quantified both in computer simulations and machine-learning approaches to determine phase diagrams, as well as combinations of the two. We review the basic workflow of free-energy computations for condensed phases, including some practical implementation advice, ranging from the Frenkel-Ladd approach to thermodynamic integration and to direct-coexistence simulations. We illustrate the applications of such methods on a range of systems from materials chemistry to biological phase separation. Finally, we outline some challenges, questions, and practical applications of phase-diagram determination which we believe are likely to be possible to address in the near future using such state-of-the-art free-energy calculations, which may provide fundamental insight into separation processes using multicomponent solvents.
Collapse
Affiliation(s)
- Pin Yu Chew
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Aleks Reinhardt
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
4
|
Mukhtyar AJ, Escobedo FA. Computing free energy barriers for the nucleation of complex network mesophases. J Chem Phys 2022; 156:034502. [DOI: 10.1063/5.0079396] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Ankita J. Mukhtyar
- Robert Frederick Smith School of Chemical and Biomolecular Engineering, Cornell University, Ithaca, New York 14850, USA
| | - Fernando A. Escobedo
- Robert Frederick Smith School of Chemical and Biomolecular Engineering, Cornell University, Ithaca, New York 14850, USA
| |
Collapse
|
5
|
Kar A, Bhati A, Lokanathan M, Bahadur V. Faster Nucleation of Ice at the Three-Phase Contact Line: Influence of Interfacial Chemistry. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2021; 37:12673-12680. [PMID: 34694119 DOI: 10.1021/acs.langmuir.1c02044] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Controlling the nucleation of ice is important in many areas including atmospheric sciences, cryopreservation, food science, and infrastructure protection. Presently, we conduct controlled experiments and analysis to uncover the influence of surface chemistry at the three-phase line on ice nucleation. We show that ice nucleation is faster upon replacing the air at the water-air interface with oils like silicone oil and almond oil. We show via statistically meaningful and carefully designed experiments that ice nucleation occurs at a higher temperature at an aluminum-water-silicone oil (or almond oil) interface as compared to an aluminum-water-air interface. We show that the location of ice nucleation can be controlled (in situations with multiple locations for ice nucleation) by controlling the interfacial chemistry at the three-phase line. We develop a model (which utilizes classical nucleation theory) to study the combined influence of two interfaces on a seed crystal of ice originating at the three-phase contact line. This model can evaluate the thermodynamic competition between nucleation at the three -phase line and heterogeneous nucleation at an interface. The model shows that three-phase contact lines usually result in a higher driving force than heterogeneous nucleation, which speeds up nucleation kinetics. Overall, our experiments and modeling uncover several useful insights into the influence of three-phase lines on nucleation during contact freezing.
Collapse
Affiliation(s)
- Aritra Kar
- Walker Department of Mechanical Engineering, The University of Texas at Austin, 204 E. Dean Keeton Street, Austin, Texas 78712, United States
| | - Awan Bhati
- Walker Department of Mechanical Engineering, The University of Texas at Austin, 204 E. Dean Keeton Street, Austin, Texas 78712, United States
| | - Manojkumar Lokanathan
- Walker Department of Mechanical Engineering, The University of Texas at Austin, 204 E. Dean Keeton Street, Austin, Texas 78712, United States
| | - Vaibhav Bahadur
- Walker Department of Mechanical Engineering, The University of Texas at Austin, 204 E. Dean Keeton Street, Austin, Texas 78712, United States
| |
Collapse
|
6
|
Reinhardt A, Cheng B. Quantum-mechanical exploration of the phase diagram of water. Nat Commun 2021; 12:588. [PMID: 33500405 PMCID: PMC7838264 DOI: 10.1038/s41467-020-20821-w] [Citation(s) in RCA: 25] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2020] [Accepted: 12/21/2020] [Indexed: 11/10/2022] Open
Abstract
The set of known stable phases of water may not be complete, and some of the phase boundaries between them are fuzzy. Starting from liquid water and a comprehensive set of 50 ice structures, we compute the phase diagram at three hybrid density-functional-theory levels of approximation, accounting for thermal and nuclear fluctuations as well as proton disorder. Such calculations are only made tractable because we combine machine-learning methods and advanced free-energy techniques. The computed phase diagram is in qualitative agreement with experiment, particularly at pressures ≲ 8000 bar, and the discrepancy in chemical potential is comparable with the subtle uncertainties introduced by proton disorder and the spread between the three hybrid functionals. None of the hypothetical ice phases considered is thermodynamically stable in our calculations, suggesting the completeness of the experimental water phase diagram in the region considered. Our work demonstrates the feasibility of predicting the phase diagram of a polymorphic system from first principles and provides a thermodynamic way of testing the limits of quantum-mechanical calculations.
Collapse
Affiliation(s)
- Aleks Reinhardt
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge, CB2 1EW, UK.
| | - Bingqing Cheng
- Accelerate Programme for Scientific Discovery, Department of Computer Science and Technology, 15 J.J. Thomson Avenue, Cambridge, CB3 0FD, UK. .,Cavendish Laboratory, University of Cambridge, J.J. Thomson Avenue, Cambridge, CB3 0HE, UK.
| |
Collapse
|
7
|
Tang X, Tian F, Xu T, Li L, Reinhardt A. Numerical calculation of free-energy barriers for entangled polymer nucleation. J Chem Phys 2020; 152:224904. [PMID: 32534553 DOI: 10.1063/5.0009716] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The crystallization of entangled polymers from their melt is investigated using computer simulation with a coarse-grained model. Using hybrid Monte Carlo simulations enables us to probe the behavior of long polymer chains. We identify solid-like beads with a centrosymmetry local order parameter and compute the nucleation free-energy barrier at relatively high supercooling with adaptive-bias windowed umbrella sampling. Our results demonstrate that the critical nucleus sizes and the heights of free-energy barriers do not significantly depend on the molecular weight of the polymer; however, the nucleation rate decreases with the increase in molecular weight. Moreover, an analysis of the composition of the critical nucleus suggests that intra-molecular growth of the nucleated cluster does not contribute significantly to crystallization for this system.
Collapse
Affiliation(s)
- Xiaoliang Tang
- National Synchrotron Radiation Lab, CAS Key Laboratory of Soft Matter Chemistry, Anhui Provincial Engineering Laboratory of Advanced Functional Polymer Film, University of Science and Technology of China, Hefei 230026, China
| | - Fucheng Tian
- National Synchrotron Radiation Lab, CAS Key Laboratory of Soft Matter Chemistry, Anhui Provincial Engineering Laboratory of Advanced Functional Polymer Film, University of Science and Technology of China, Hefei 230026, China
| | - Tingyu Xu
- National Synchrotron Radiation Lab, CAS Key Laboratory of Soft Matter Chemistry, Anhui Provincial Engineering Laboratory of Advanced Functional Polymer Film, University of Science and Technology of China, Hefei 230026, China
| | - Liangbin Li
- National Synchrotron Radiation Lab, CAS Key Laboratory of Soft Matter Chemistry, Anhui Provincial Engineering Laboratory of Advanced Functional Polymer Film, University of Science and Technology of China, Hefei 230026, China
| | - Aleks Reinhardt
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
8
|
Leoni F, Shi R, Tanaka H, Russo J. Crystalline clusters in mW water: Stability, growth, and grain boundaries. J Chem Phys 2019; 151:044505. [DOI: 10.1063/1.5100812] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Fabio Leoni
- School of Mathematics, University of Bristol, Bristol BS8 1TW, United Kingdom
| | - Rui Shi
- Department of Fundamental Engineering, Institute of Industrial Science, University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| | - Hajime Tanaka
- Department of Fundamental Engineering, Institute of Industrial Science, University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| | - John Russo
- School of Mathematics, University of Bristol, Bristol BS8 1TW, United Kingdom
- Department of Fundamental Engineering, Institute of Industrial Science, University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| |
Collapse
|
9
|
Schiffer J, Mael LE, Prather KA, Amaro RE, Grassian VH. Sea Spray Aerosol: Where Marine Biology Meets Atmospheric Chemistry. ACS CENTRAL SCIENCE 2018; 4:1617-1623. [PMID: 30648145 PMCID: PMC6311946 DOI: 10.1021/acscentsci.8b00674] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/21/2018] [Indexed: 05/25/2023]
Abstract
Atmospheric aerosols have long been known to alter climate by scattering incoming solar radiation and acting as seeds for cloud formation. These processes have vast implications for controlling the chemistry of our environment and the Earth's climate. Sea spray aerosol (SSA) is emitted over nearly three-quarters of our planet, yet precisely how SSA impacts Earth's radiation budget remains highly uncertain. Over the past several decades, studies have shown that SSA particles are far more complex than just sea salt. Ocean biological and physical processes produce individual SSA particles containing a diverse array of biological species including proteins, enzymes, bacteria, and viruses and a diverse array of organic compounds including fatty acids and sugars. Thus, a new frontier of research is emerging at the nexus of chemistry, biology, and atmospheric science. In this Outlook article, we discuss how current and future aerosol chemistry research demands a tight coupling between experimental (observational and laboratory studies) and computational (simulation-based) methods. This integration of approaches will enable the systematic interrogation of the complexity within individual SSA particles at a level that will enable prediction of the physicochemical properties of real-world SSA, ultimately illuminating the detailed mechanisms of how the constituents within individual SSA impact climate.
Collapse
Affiliation(s)
- Jamie
M. Schiffer
- Department of Chemistry and Biochemistry and Department of Nanoengineering, University of California, San Diego, 9500 Gilman Drive, La Jolla, California 92093-0378, United States
| | - Liora E. Mael
- Department of Chemistry and Biochemistry and Department of Nanoengineering, University of California, San Diego, 9500 Gilman Drive, La Jolla, California 92093-0378, United States
| | - Kimberly A. Prather
- Department of Chemistry and Biochemistry and Department of Nanoengineering, University of California, San Diego, 9500 Gilman Drive, La Jolla, California 92093-0378, United States
- Scripps
Institution of Oceanography, University
of California, San Diego, La Jolla, California 92093, United States
| | - Rommie E. Amaro
- Department of Chemistry and Biochemistry and Department of Nanoengineering, University of California, San Diego, 9500 Gilman Drive, La Jolla, California 92093-0378, United States
| | - Vicki H. Grassian
- Department of Chemistry and Biochemistry and Department of Nanoengineering, University of California, San Diego, 9500 Gilman Drive, La Jolla, California 92093-0378, United States
- Scripps
Institution of Oceanography, University
of California, San Diego, La Jolla, California 92093, United States
| |
Collapse
|
10
|
Espinosa JR, Navarro C, Sanz E, Valeriani C, Vega C. On the time required to freeze water. J Chem Phys 2018; 145:211922. [PMID: 28799362 DOI: 10.1063/1.4965427] [Citation(s) in RCA: 45] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
Abstract
By using the seeding technique the nucleation rate for the formation of ice at room pressure will be estimated for the TIP4P/ICE model using longer runs and a smaller grid of temperatures than in the previous work. The growth rate of ice will be determined for TIP4P/ICE and for the mW model of water. Although TIP4P/ICE and mW have a similar melting point and melting enthalpy, they differ significantly in the dynamics of freezing. The nucleation rate of mW is lower than that of TIP4P/ICE due to its higher interfacial free energy. Experimental results for the nucleation rate of ice are between the predictions of these two models when obtained from the seeding technique, although closer to the predictions of TIP4P/ICE. The growth rate of ice for the mW model is four orders of magnitude larger than for TIP4P/ICE. Avrami's expression is used to estimate the crystallization time from the values of the nucleation and growth rates. For mW the minimum in the crystallization time is found at approximately 85 K below the melting point and its value is of about a few ns, in agreement with the results obtained from brute force simulations by Moore and Molinero. For the TIP4P/ICE the minimum is found at about 55 K below the melting point, but its value is about ten microseconds. This value is compatible with the minimum cooling rate required to avoid the formation of ice and obtaining a glass phase. The crossover from the nucleation controlled crystallization to the growth controlled crystallization will be discussed for systems of finite size. This crossover could explain the apparent discrepancy between the values of J obtained by different experimental groups for temperatures below 230 K and should be considered as an alternative hypothesis to the two previously suggested: internal pressure and/or surface freezing effects. A maximum in the compressibility was found for the TIP4P/ICE model in supercooled water. The relaxation time is much smaller than the crystallization time at the temperature at which this maximum occurs, so this maximum is a real thermodynamic feature of the model. At the temperature of minimum crystallization time, the crystallization time is larger than the relaxation time by just two orders of magnitude.
Collapse
Affiliation(s)
- J R Espinosa
- Departamento de Quimica Fisica, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain and Departamento de Fisica Aplicada I , Facultad de Ciencias Fisicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - C Navarro
- Departamento de Quimica Fisica, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain and Departamento de Fisica Aplicada I , Facultad de Ciencias Fisicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - E Sanz
- Departamento de Quimica Fisica, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain and Departamento de Fisica Aplicada I , Facultad de Ciencias Fisicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - C Valeriani
- Departamento de Quimica Fisica, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain and Departamento de Fisica Aplicada I , Facultad de Ciencias Fisicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - C Vega
- Departamento de Quimica Fisica, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain and Departamento de Fisica Aplicada I , Facultad de Ciencias Fisicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
11
|
Guo J, Haji-Akbari A, Palmer JC. Hybrid Monte Carlo with LAMMPS. JOURNAL OF THEORETICAL & COMPUTATIONAL CHEMISTRY 2018. [DOI: 10.1142/s0219633618400023] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
We describe a strategy for performing canonical and isothermal-isobaric ensemble hybrid Monte Carlo (HMC) simulations with the widely-used Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) molecular dynamics (MD) software package. The overall workflow for the HMC simulations is handled using an external Python driver script, which invokes LAMMPS’ library interface to perform numerically intensive tasks such as MD integration. We document several rigorous consistency checks that have been used to validate our HMC implementation. We also demonstrate that our approach can be readily extended to implement biased HMC sampling schemes for computing free energies. Codes and input files from the documented examples are available on the web.
Collapse
Affiliation(s)
- Jingxiang Guo
- Department of Chemical and Biomolecular Engineering, University of Houston, Houston, Texas 77204, USA
| | - Amir Haji-Akbari
- Department of Chemical and Environmental Engineering, Yale University, New Haven, CT 06520, USA
| | - Jeremy C. Palmer
- Department of Chemical and Biomolecular Engineering, University of Houston, Houston, Texas 77204, USA
| |
Collapse
|
12
|
Zhang Z, Liu XY. Control of ice nucleation: freezing and antifreeze strategies. Chem Soc Rev 2018; 47:7116-7139. [DOI: 10.1039/c8cs00626a] [Citation(s) in RCA: 139] [Impact Index Per Article: 23.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Water freezing remains a perennial topic of great relevance to many important aspects of our lives; from the climate to human society and from economics to medicine, frozen water profoundly influences our living environment and life activities.
Collapse
Affiliation(s)
- Zhisen Zhang
- Research Institute for Biomimetics and Soft Matter
- Fujian Provincial Key Laboratory for Soft Functional Materials Research
- Department of Physics
- Department of Biomaterials
- Xiamen University
| | - Xiang-Yang Liu
- Research Institute for Biomimetics and Soft Matter
- Fujian Provincial Key Laboratory for Soft Functional Materials Research
- Department of Physics
- Department of Biomaterials
- Xiamen University
| |
Collapse
|
13
|
Ambler M, Vorselaars B, Allen MP, Quigley D. Solid–liquid interfacial free energy of ice Ih, ice Ic, and ice 0 within a mono-atomic model of water via the capillary wave method. J Chem Phys 2017; 146:074701. [DOI: 10.1063/1.4975776] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Affiliation(s)
- Michael Ambler
- Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - Bart Vorselaars
- Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - Michael P. Allen
- Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom
- H. H. Wills Physics Laboratory,
Royal Fort, Tyndall Avenue, Bristol BS8 1TL, United Kingdom
| | - David Quigley
- Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom
| |
Collapse
|
14
|
Koop T, Murray BJ. A physically constrained classical description of the homogeneous nucleation of ice in water. J Chem Phys 2016; 145:211915. [DOI: 10.1063/1.4962355] [Citation(s) in RCA: 73] [Impact Index Per Article: 9.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022] Open
Affiliation(s)
- Thomas Koop
- Faculty of Chemistry, Bielefeld University, Bielefeld, Germany
| | - Benjamin J. Murray
- Institute for Climate and Atmospheric Science, School of Earth and Environment, University of Leeds, Leeds, United Kingdom
| |
Collapse
|
15
|
Russo J, Tanaka H. Crystal nucleation as the ordering of multiple order parameters. J Chem Phys 2016; 145:211801. [DOI: 10.1063/1.4962166] [Citation(s) in RCA: 79] [Impact Index Per Article: 9.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Affiliation(s)
- John Russo
- Institute of Industrial Science, University of Tokyo, Meguro-ku, Tokyo 153-8505, Japan
- School of Mathematics, University of Bristol, Bristol BS8 1TW, United Kingdom
| | - Hajime Tanaka
- Institute of Industrial Science, University of Tokyo, Meguro-ku, Tokyo 153-8505, Japan
| |
Collapse
|
16
|
Lupi L, Peters B, Molinero V. Pre-ordering of interfacial water in the pathway of heterogeneous ice nucleation does not lead to a two-step crystallization mechanism. J Chem Phys 2016; 145:211910. [DOI: 10.1063/1.4961652] [Citation(s) in RCA: 49] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
Affiliation(s)
- Laura Lupi
- Department of Chemistry, The University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112-0850, USA
| | - Baron Peters
- Department of Chemical Engineering and Department of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106, USA
| | - Valeria Molinero
- Department of Chemistry, The University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112-0850, USA
| |
Collapse
|
17
|
Espinosa JR, Zaragoza A, Rosales-Pelaez P, Navarro C, Valeriani C, Vega C, Sanz E. Interfacial Free Energy as the Key to the Pressure-Induced Deceleration of Ice Nucleation. PHYSICAL REVIEW LETTERS 2016; 117:135702. [PMID: 27715079 DOI: 10.1103/physrevlett.117.135702] [Citation(s) in RCA: 35] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/04/2016] [Indexed: 06/06/2023]
Abstract
The avoidance of water freezing is the holy grail in the cryopreservation of biological samples, food, and organs. Fast cooling rates are used to beat ice nucleation and avoid cell damage. This strategy can be enhanced by applying high pressures to decrease the nucleation rate, but the physics behind this procedure has not been fully understood yet. We perform computer experiments to investigate ice nucleation at high pressures consisting in embedding ice seeds in supercooled water. We find that the slowing down of the nucleation rate is mainly due to an increase of the ice I-water interfacial free energy with pressure. Our work also clarifies the molecular mechanism of ice nucleation for a wide pressure range. This study is not only relevant to cryopreservation, but also to water amorphization and climate change modeling.
Collapse
Affiliation(s)
- Jorge R Espinosa
- Departamento de Quimica Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Alberto Zaragoza
- Departamento de Quimica Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Pablo Rosales-Pelaez
- Departamento de Quimica Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Caridad Navarro
- Departamento de Quimica Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Chantal Valeriani
- Departamento de Quimica Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
- Departamento de Fisica Aplicada I, Facultad de Ciencias Fisicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Carlos Vega
- Departamento de Quimica Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Eduardo Sanz
- Departamento de Quimica Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
18
|
He X, Shen Y, Hung FR, Santiso EE. Molecular simulation of homogeneous nucleation of crystals of an ionic liquid from the melt. J Chem Phys 2015; 143:124506. [PMID: 26429023 DOI: 10.1063/1.4931654] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The homogeneous nucleation of crystals of the ionic liquid [dmim(+)][Cl(-)] from its supercooled liquid phase in the bulk (P = 1 bar, T = 340 K, representing a supercooling of 58 K) was studied using molecular simulations. The string method in collective variables [Maragliano et al., J. Chem. Phys. 125, 024106 (2006)] was used in combination with Markovian milestoning with Voronoi tessellations [Maragliano et al., J. Chem. Theory Comput. 5, 2589-2594 (2009)] and order parameters for molecular crystals [E. E. Santiso and B. L. Trout, J. Chem. Phys. 134, 064109 (2011)] to sketch a minimum free energy path connecting the supercooled liquid and the monoclinic crystal phases, and to determine the free energy and the rates involved in the homogeneous nucleation process. The physical significance of the configurations found along this minimum free energy path is discussed with the help of calculations based on classical nucleation theory and with additional simulation results obtained for a larger system. Our results indicate that, at a supercooling of 58 K, the liquid has to overcome a free energy barrier of the order of 60 kcal/mol and to form a critical nucleus with an average size of about 3.6 nm, before it reaches the thermodynamically stable crystal phase. A simulated homogeneous nucleation rate of 5.0 × 10(10) cm(-3) s(-1) was obtained for our system, which is in reasonable agreement with experimental and simulation rates for homogeneous nucleation of ice at similar degrees of supercooling. This study represents our first step in a series of studies aimed at understanding the nucleation and growth of crystals of organic salts near surfaces and inside nanopores.
Collapse
Affiliation(s)
- Xiaoxia He
- Cain Department of Chemical Engineering, Louisiana State University, Baton Rouge, Louisiana 70803, USA
| | - Yan Shen
- Cain Department of Chemical Engineering, Louisiana State University, Baton Rouge, Louisiana 70803, USA
| | - Francisco R Hung
- Cain Department of Chemical Engineering, Louisiana State University, Baton Rouge, Louisiana 70803, USA
| | - Erik E Santiso
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina 27695, USA
| |
Collapse
|
19
|
Fraux G, Doye JPK. Note: Heterogeneous ice nucleation on silver-iodide-like surfaces. J Chem Phys 2015; 141:216101. [PMID: 25481172 DOI: 10.1063/1.4902382] [Citation(s) in RCA: 52] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
Affiliation(s)
- Guillaume Fraux
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, Oxford OX1 3QZ, United Kingdom
| | - Jonathan P K Doye
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, Oxford OX1 3QZ, United Kingdom
| |
Collapse
|
20
|
Zaragoza A, Conde MM, Espinosa JR, Valeriani C, Vega C, Sanz E. Competition between ices Ih and Ic in homogeneous water freezing. J Chem Phys 2015; 143:134504. [DOI: 10.1063/1.4931987] [Citation(s) in RCA: 55] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Affiliation(s)
- Alberto Zaragoza
- Departamento de Quimica Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Maria M. Conde
- Departamento de Quimica Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
- Laboratoire des Fluides Complexes et leurs Réservoirs, UMR 5150, Université de Pau et des Pays de l’Adour, B. P. 1155, Pau-Cedex 64013, France
| | - Jorge R. Espinosa
- Departamento de Quimica Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Chantal Valeriani
- Departamento de Quimica Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
- Departamento de Fisica Aplicada I, Facultad de Ciencias Fisicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Carlos Vega
- Departamento de Quimica Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Eduardo Sanz
- Departamento de Quimica Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
21
|
Quigley D, Alfè D, Slater B. Communication: On the stability of ice 0, ice i, and I(h). J Chem Phys 2015; 141:161102. [PMID: 25362263 DOI: 10.1063/1.4900772] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
Using ab initio methods, we examine the stability of ice 0, a recently proposed tetragonal form of ice implicated in the homogeneous freezing of water [J. Russo, F. Romano, and H. Tanaka, Nat. Mater. 13, 670 (2014)]. Vibrational frequencies are computed across the complete Brillouin Zone using Density Functional Theory (DFT), to confirm mechanical stability and quantify the free energy of ice 0 relative to ice I(h). The robustness of this result is tested via dispersion corrected semi-local and hybrid DFT, and Quantum Monte-Carlo calculation of lattice energies. Results indicate that popular molecular models only slightly overestimate the stability of ice zero. In addition, we study all possible realisations of proton disorder within the ice zero unit cell, and identify the ground state as ferroelectric. Comparisons are made to other low density metastable forms of ice, suggesting that the ice i structure [C. J. Fennel and J. D. Gezelter, J. Chem. Theory Comput. 1, 662 (2005)] may be equally relevant to ice formation.
Collapse
Affiliation(s)
- D Quigley
- Department of Physics and Centre for Scientific Computing, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - D Alfè
- Department of Earth Sciences and Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - B Slater
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| |
Collapse
|
22
|
Espinosa JR, Vega C, Valeriani C, Sanz E. The crystal-fluid interfacial free energy and nucleation rate of NaCl from different simulation methods. J Chem Phys 2015; 142:194709. [DOI: 10.1063/1.4921185] [Citation(s) in RCA: 49] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
Affiliation(s)
- Jorge R. Espinosa
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Carlos Vega
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Chantal Valeriani
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Eduardo Sanz
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
23
|
Malkin TL, Murray BJ, Salzmann CG, Molinero V, Pickering SJ, Whale TF. Stacking disorder in ice I. Phys Chem Chem Phys 2015; 17:60-76. [PMID: 25380218 DOI: 10.1039/c4cp02893g] [Citation(s) in RCA: 141] [Impact Index Per Article: 15.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Traditionally, ice I was considered to exist in two well-defined crystalline forms at ambient pressure: stable hexagonal ice (ice Ih) and metastable cubic ice (ice Ic). However, it is becoming increasingly evident that what has been called cubic ice in the past does not have a structure consistent with the cubic crystal system. Instead, it is a stacking-disordered material containing cubic sequences interlaced with hexagonal sequences, which is termed stacking-disordered ice (ice Isd). In this article, we summarise previous work on ice with stacking disorder including ice that was called cubic ice in the past. We also present new experimental data which shows that ice which crystallises after heterogeneous nucleation in water droplets containing solid inclusions also contains stacking disorder even at freezing temperatures of around -15 °C. This supports the results from molecular simulations, that the structure of ice that crystallises initially from supercooled water is always stacking-disordered and that this metastable ice can transform to the stable hexagonal phase subject to the kinetics of recrystallization. We also show that stacking disorder in ice which forms from water droplets is quantitatively distinct from ice made via other routes. The emerging picture of ice I is that of a very complex material which frequently contains stacking disorder and this stacking disorder can vary in complexity depending on the route of formation and thermal history.
Collapse
Affiliation(s)
- Tamsin L Malkin
- Institute for Climate and Atmospheric Science, School of Earth and Environment, University of Leeds, Woodhouse Lane, Leeds, LS2 9JT, UK.
| | | | | | | | | | | |
Collapse
|
24
|
Wang X, Wang S, Xu Q, Mi J. Thermodynamics of Ice Nucleation in Liquid Water. J Phys Chem B 2015; 119:1660-8. [DOI: 10.1021/jp512280p] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Xin Wang
- State Key Laboratory of Organic-Inorganic
Composites, Beijing University of Chemical Technology, Beijing 100029, China
| | - Shui Wang
- State Key Laboratory of Organic-Inorganic
Composites, Beijing University of Chemical Technology, Beijing 100029, China
| | - Qinzhi Xu
- Institute of Microelectronics of Chinese Academy of Sciences, Beijing, 100029, China
| | - Jianguo Mi
- State Key Laboratory of Organic-Inorganic
Composites, Beijing University of Chemical Technology, Beijing 100029, China
| |
Collapse
|
25
|
Ickes L, Welti A, Hoose C, Lohmann U. Classical nucleation theory of homogeneous freezing of water: thermodynamic and kinetic parameters. Phys Chem Chem Phys 2015; 17:5514-37. [DOI: 10.1039/c4cp04184d] [Citation(s) in RCA: 114] [Impact Index Per Article: 12.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Different formulations of the kinetic and thermodynamic parameters of CNT are evaluated against measured nucleation rates.
Collapse
Affiliation(s)
- Luisa Ickes
- Institute of Atmospheric and Climate Science
- ETH Zurich
- Zurich
- Switzerland
| | - André Welti
- Institute of Atmospheric and Climate Science
- ETH Zurich
- Zurich
- Switzerland
| | - Corinna Hoose
- Institute for Meteorology and Climate Research
- Karlsruhe Institute of Technology
- Karlsruhe
- Germany
| | - Ulrike Lohmann
- Institute of Atmospheric and Climate Science
- ETH Zurich
- Zurich
- Switzerland
| |
Collapse
|
26
|
Espinosa JR, Sanz E, Valeriani C, Vega C. Homogeneous ice nucleation evaluated for several water models. J Chem Phys 2014; 141:18C529. [DOI: 10.1063/1.4897524] [Citation(s) in RCA: 108] [Impact Index Per Article: 10.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- J. R. Espinosa
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - E. Sanz
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - C. Valeriani
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - C. Vega
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
27
|
Reinhardt A, Doye JPK. Effects of surface interactions on heterogeneous ice nucleation for a monatomic water model. J Chem Phys 2014; 141:084501. [DOI: 10.1063/1.4892804] [Citation(s) in RCA: 49] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022] Open
Affiliation(s)
- Aleks Reinhardt
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, Oxford OX1 3QZ, United Kingdom
| | - Jonathan P. K. Doye
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, Oxford OX1 3QZ, United Kingdom
| |
Collapse
|
28
|
Sanz E, Vega C, Espinosa JR, Caballero-Bernal R, Abascal JLF, Valeriani C. Homogeneous Ice Nucleation at Moderate Supercooling from Molecular Simulation. J Am Chem Soc 2013; 135:15008-17. [DOI: 10.1021/ja4028814] [Citation(s) in RCA: 214] [Impact Index Per Article: 19.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Affiliation(s)
- E. Sanz
- Departamento
de Química
Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - C. Vega
- Departamento
de Química
Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - J. R. Espinosa
- Departamento
de Química
Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - R. Caballero-Bernal
- Departamento
de Química
Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - J. L. F. Abascal
- Departamento
de Química
Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - C. Valeriani
- Departamento
de Química
Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|