1
|
Soni A, Patey GN. Using machine learning with atomistic surface and local water features to predict heterogeneous ice nucleation. J Chem Phys 2024; 160:124501. [PMID: 38530008 DOI: 10.1063/5.0177706] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2023] [Accepted: 03/04/2024] [Indexed: 03/27/2024] Open
Abstract
Heterogeneous ice nucleation (HIN) has applications in climate science, nanotechnology, and cryopreservation. Ice nucleation on the earth's surface or in the atmosphere usually occurs heterogeneously involving foreign substrates, known as ice nucleating particles (INPs). Experiments identify good INPs but lack sufficient microscopic resolution to answer the basic question: What makes a good INP? We employ molecular dynamics (MD) simulations in combination with machine learning (ML) to address this question. Often, the large amount of computational cost required to cross the nucleation barrier and observe HIN in MD simulations is a practical limitation. We use information obtained from short MD simulations of atomistic surface and water models to predict the likelihood of HIN. We consider 153 atomistic substrates with some surfaces differing in elemental composition and others only in terms of lattice parameters, surface morphology, or surface charges. A range of water features near the surface (local) are extracted from short MD simulations over a time interval (≤300 ns) where ice nucleation has not initiated. Three ML classification models, Random Forest (RF), support vector machine, and Gaussian process classification are considered, and the accuracies achieved by all three approaches lie within their statistical uncertainties. Including local water features is essential for accurate prediction. The accuracy of our best RF classification model obtained including both surface and local water features is 0.89 ± 0.05. A similar accuracy can be achieved including only local water features, suggesting that the important surface properties are largely captured by the local water features. Some important features identified by ML analysis are local icelike structures, water density and polarization profiles perpendicular to the surface, and the two-dimensional lattice match to ice. We expect that this work, with its strong focus on realistic surface models, will serve as a guide to the identification or design of substrates that can promote or discourage ice nucleation.
Collapse
Affiliation(s)
- Abhishek Soni
- Department of Chemistry, University of British Columbia, Vancouver, British Columbia V6T 1Z1, Canada
| | - G N Patey
- Department of Chemistry, University of British Columbia, Vancouver, British Columbia V6T 1Z1, Canada
| |
Collapse
|
2
|
Roudsari G, Lbadaoui-Darvas M, Welti A, Nenes A, Laaksonen A. The molecular scale mechanism of deposition ice nucleation on silver iodide. ENVIRONMENTAL SCIENCE: ATMOSPHERES 2024; 4:243-251. [PMID: 38371604 PMCID: PMC10867811 DOI: 10.1039/d3ea00140g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/21/2023] [Accepted: 12/02/2023] [Indexed: 02/20/2024]
Abstract
Heterogeneous ice nucleation is a ubiquitous process in the natural and built environment. Deposition ice nucleation, i.e. heterogeneous ice nucleation that - according to the traditional view - occurs in a subsaturated water vapor environment and in the absence of supercooled water on the solid, ice-forming surface, is among the most important ice formation processes in high-altitude cirrus and mixed-phase clouds. Despite its importance, very little is known about the mechanism of deposition ice nucleation at the microscopic level. This study puts forward an adsorption-based mechanism for deposition ice nucleation through results from a combination of atomistic simulations, experiments and theoretical modelling. One of the most potent laboratory surrogates of ice nucleating particles, silver iodide, is used as a substrate for the simulations. We find that water initially adsorbs in clusters which merge and grow over time to form layers of supercooled water. Ice nucleation on silver iodide requires at minimum the adsorption of 4 molecular layers of water. Guided by the simulations we propose the following fundamental freezing steps: (1) Water molecules adsorb on the surface, forming nanodroplets. (2) The supercooled water nanodroplets merge into a continuous multilayer when they grow to about 3 molecular layers thick. (3) The layer continues to grow until the critical thickness for freezing is reached. (4) The critical ice cluster continues to grow.
Collapse
Affiliation(s)
| | - Mária Lbadaoui-Darvas
- Laboratory of Atmospheric Processes and their Impacts, ENAC, Ecole Polytechnique Fédérale de Lausanne Lausanne Switzerland
- Institute of Chemical Engineering Sciences, Foundation for Research and Technology Hellas (FORTH/ICE-HT) 26504 Patras Greece
| | - André Welti
- Finnish Meteorological Institute FI-00101 Helsinki Finland
| | - Athanasios Nenes
- Laboratory of Atmospheric Processes and their Impacts, ENAC, Ecole Polytechnique Fédérale de Lausanne Lausanne Switzerland
- Institute of Chemical Engineering Sciences, Foundation for Research and Technology Hellas (FORTH/ICE-HT) 26504 Patras Greece
| | - Ari Laaksonen
- Finnish Meteorological Institute FI-00101 Helsinki Finland
- Department of Applied Physics, University of Eastern Finland Kuopio 70211 Finland
| |
Collapse
|
3
|
Dickbreder T, Sabath F, Reischl B, Nilsson RVE, Foster AS, Bechstein R, Kühnle A. Atomic structure and water arrangement on K-feldspar microcline (001). NANOSCALE 2024; 16:3462-3473. [PMID: 38214028 DOI: 10.1039/d3nr05585j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/13/2024]
Abstract
The properties of clouds, such as their reflectivity or their likelihood to precipitate, depend on whether the cloud droplets are liquid or frozen. Thus, understanding the ice nucleation mechanisms is essential for the development of reliable climate models. Most ice nucleation in the atmosphere is heterogeneous, i.e., caused by ice nucleating particles such as mineral dusts or organic aerosols. In this regard, K-feldspar minerals have attracted great interest recently as they have been identified as one of the most important ice nucleating particles under mixed-phase cloud conditions. The mechanism by which feldspar minerals facilitate ice nucleation remains, however, elusive. Here, we present atomic force microscopy (AFM) experiments on microcline (001) performed in an ultrahigh vacuum and at the solid-water interface together with density functional theory (DFT) and molecular dynamics (MD) calculations. Our ultrahigh vacuum data reveal features consistent with a hydroxyl-terminated surface. This finding suggests that water in the residual gas readily reacts with the surface. Indeed, the corresponding DFT calculations confirm a dissociative water adsorption. Three-dimensional AFM measurements performed at the mineral-water interface unravel a layered hydration structure with two features per surface unit cell. A comparison with MD calculations suggests that the structure observed in AFM corresponds to the second hydration layer rather than the first water layer. In agreement with previous computation results, no ice-like structure is seen, questioning an explanation of the ice nucleation ability by lattice match. Our results provide an atomic-scale benchmark for the clean and water-covered microcline (001) plane, which is mandatory for understanding the ice nucleation mechanism on feldspar minerals.
Collapse
Affiliation(s)
- Tobias Dickbreder
- Faculty of Chemistry, Physical Chemistry I, Bielefeld University, 33615 Bielefeld, Germany.
| | - Franziska Sabath
- Faculty of Chemistry, Physical Chemistry I, Bielefeld University, 33615 Bielefeld, Germany.
| | - Bernhard Reischl
- Institute for Atmospheric and Earth System Research (INAR)/Physics, Faculty of Science, University of Helsinki, 00014 Helsinki, Finland
| | - Rasmus V E Nilsson
- Institute for Atmospheric and Earth System Research (INAR)/Physics, Faculty of Science, University of Helsinki, 00014 Helsinki, Finland
| | - Adam S Foster
- Department of Applied Physics, Aalto University, Finland
- Nano Life Science Institute (WPI-NanoLSI), Kanazawa University, Kanazawa 920-1192, Japan
| | - Ralf Bechstein
- Faculty of Chemistry, Physical Chemistry I, Bielefeld University, 33615 Bielefeld, Germany.
| | - Angelika Kühnle
- Faculty of Chemistry, Physical Chemistry I, Bielefeld University, 33615 Bielefeld, Germany.
| |
Collapse
|
4
|
Piaggi PM, Selloni A, Panagiotopoulos AZ, Car R, Debenedetti PG. A first-principles machine-learning force field for heterogeneous ice nucleation on microcline feldspar. Faraday Discuss 2024; 249:98-113. [PMID: 37791889 DOI: 10.1039/d3fd00100h] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/05/2023]
Abstract
The formation of ice in the atmosphere affects precipitation and cloud properties, and plays a key role in the climate of our planet. Although ice can form directly from liquid water under deeply supercooled conditions, the presence of foreign particles can aid ice formation at much warmer temperatures. Over the past decade, experiments have highlighted the remarkable efficiency of feldspar minerals as ice nuclei compared to other particles present in the atmosphere. However, the exact mechanism of ice formation on feldspar surfaces has yet to be fully understood. Here, we develop a first-principles machine-learning model for the potential energy surface aimed at studying ice nucleation at microcline feldspar surfaces. The model is able to reproduce with high-fidelity the energies and forces derived from density-functional theory (DFT) based on the SCAN exchange and correlation functional. Our training set includes configurations of bulk supercooled water, hexagonal and cubic ice, microcline, and fully-hydroxylated feldspar surfaces exposed to a vacuum, liquid water, and ice. We apply the machine-learning force field to study different fully-hydroxylated terminations of the (100), (010), and (001) surfaces of microcline exposed to a vacuum. Our calculations suggest that terminations that do not minimize the number of broken bonds are preferred in a vacuum. We also study the structure of supercooled liquid water in contact with microcline surfaces, and find that water density correlations extend up to around 10 Å from the surfaces. Finally, we show that the force field maintains a high accuracy during the simulation of ice formation at microcline surfaces, even for large systems of around 30 000 atoms. Future work will be directed towards the calculation of nucleation free-energy barriers and rates using the force field developed herein, and understanding the role of different microcline surfaces in ice nucleation.
Collapse
Affiliation(s)
- Pablo M Piaggi
- Department of Chemistry, Princeton University, Princeton, NJ 08544, USA.
| | - Annabella Selloni
- Department of Chemistry, Princeton University, Princeton, NJ 08544, USA.
| | | | - Roberto Car
- Department of Chemistry, Princeton University, Princeton, NJ 08544, USA.
- Department of Physics, Princeton University, Princeton, NJ 08544, USA
| | - Pablo G Debenedetti
- Department of Chemical and Biological Engineering, Princeton University, Princeton, NJ 08544, USA
| |
Collapse
|
5
|
Franceschi G, Conti A, Lezuo L, Abart R, Mittendorfer F, Schmid M, Diebold U. How Water Binds to Microcline Feldspar (001). J Phys Chem Lett 2024; 15:15-22. [PMID: 38156776 PMCID: PMC10788961 DOI: 10.1021/acs.jpclett.3c03235] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2023] [Revised: 12/12/2023] [Accepted: 12/14/2023] [Indexed: 01/03/2024]
Abstract
Microcline feldspar (KAlSi3O8) is a common mineral with important roles in Earth's ecological balance. It participates in carbon, potassium, and water cycles, contributing to CO2 sequestration, soil formation, and atmospheric ice nucleation. To understand the fundamentals of these processes, it is essential to establish microcline's surface atomic structure and its interaction with the omnipresent water molecules. This work presents atomic-scale results on microcline's lowest-energy surface and its interaction with water, combining ultrahigh vacuum investigations by noncontact atomic force microscopy and X-ray photoelectron spectroscopy with density functional theory calculations. An ordered array of hydroxyls bonded to silicon or aluminum readily forms on the cleaved surface at room temperature. The distinct proton affinities of these hydroxyls influence the arrangement and orientation of the first water molecules binding to the surface, holding potential implications for the subsequent condensation of water.
Collapse
Affiliation(s)
| | - Andrea Conti
- Institute
of Applied Physics, TU Wien, 1040 Vienna, Austria
| | - Luca Lezuo
- Institute
of Applied Physics, TU Wien, 1040 Vienna, Austria
| | - Rainer Abart
- Department
of Lithospheric Research, Universität
Wien, 1090 Vienna, Austria
| | | | - Michael Schmid
- Institute
of Applied Physics, TU Wien, 1040 Vienna, Austria
| | - Ulrike Diebold
- Institute
of Applied Physics, TU Wien, 1040 Vienna, Austria
| |
Collapse
|
6
|
Zhu X, Vandamme M, Jiang Z, Brochard L. Molecular simulation of the confined crystallization of ice in cement nanopore. J Chem Phys 2023; 159:154704. [PMID: 37850696 DOI: 10.1063/5.0169783] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2023] [Accepted: 09/26/2023] [Indexed: 10/19/2023] Open
Abstract
Freezing of water under nanoconfinement exhibits physical peculiarities with respect to the bulk water. However, experimental observations are extremely challenging at this scale, which limits our understanding of the effect of confinement on water properties upon freezing. In this study, we use molecular dynamic simulations to investigate how confinement affects the kinetics of growth of ice and the thermodynamic equilibrium of ice-liquid coexistence. TIP4P/Ice water model and CSH-FF model were applied to simulate ice crystallization in a confined cement system at temperatures down to 220 K. We adapted an interface detection algorithm and reparameterized the CHILL/CHILL+ algorithm to capture ice growth. The confinement leads to a shift of the maximum growth rate of ice to a higher temperature than for bulk water. Both the confinement and surface impurities contribute to slowing down the ice growth. For the ice-liquid coexistence at equilibrium, we derive a formulation of Thomson's equation adapted to statistical physics quantities accessible by molecular simulation, and we show that this adapted equation predicts accurately the melting line of bulk and confined ice Ih as a function of pressure. The confinement decreases systematically the melting temperature of ice of about 5 K compared with bulk ice Ih. A premelted water film about 1 nm thick is observed between the solid wall and ice, and its thickness is found to decrease continuously as temperature is lowered. We note that the surface impurities are key to the formation of the premelted water nanofilm when the temperature is lower than 250 K.
Collapse
Affiliation(s)
- Xinping Zhu
- Key Laboratory of Advanced Civil Engineering Materials of Ministry of Education, School of Materials Science and Engineering, Tongji University, Shanghai 201804, China
- Navier, Ecole des Ponts, Univ. Gustave Eiffel, CNRS, Marne-la-Vallée, France
| | - Matthieu Vandamme
- Navier, Ecole des Ponts, Univ. Gustave Eiffel, CNRS, Marne-la-Vallée, France
| | - Zhengwu Jiang
- Key Laboratory of Advanced Civil Engineering Materials of Ministry of Education, School of Materials Science and Engineering, Tongji University, Shanghai 201804, China
| | - Laurent Brochard
- Navier, Ecole des Ponts, Univ. Gustave Eiffel, CNRS, Marne-la-Vallée, France
| |
Collapse
|
7
|
Marks SM, Vicars Z, Thosar AU, Patel AJ. Characterizing Surface Ice-Philicity Using Molecular Simulations and Enhanced Sampling. J Phys Chem B 2023. [PMID: 37378637 DOI: 10.1021/acs.jpcb.3c01627] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/29/2023]
Abstract
The formation of ice, which plays an important role in diverse contexts ranging from cryopreservation to atmospheric science, is often mediated by solid surfaces. Although surfaces that interact favorably with ice (relative to liquid water) can facilitate ice formation by lowering nucleation barriers, the molecular characteristics that confer icephilicity to a surface are complex and incompletely understood. To address this challenge, here we introduce a robust and computationally efficient method for characterizing surface ice-philicity that combines molecular simulations and enhanced sampling techniques to quantify the free energetic cost of increasing surface-ice contact at the expense of surface-water contact. Using this method to characterize the ice-philicity of a family of model surfaces that are lattice matched with ice but vary in their polarity, we find that the nonpolar surfaces are moderately ice-phobic, whereas the polar surfaces are highly ice-philic. In contrast, for surfaces that display no complementarity to the ice lattice, we find that ice-philicity is independent of surface polarity and that both nonpolar and polar surfaces are moderately ice-phobic. Our work thus provides a prescription for quantitatively characterizing surface ice-philicity and sheds light on how ice-philicity is influenced by lattice matching and polarity.
Collapse
Affiliation(s)
- Sean M Marks
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Zachariah Vicars
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Aniket U Thosar
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Amish J Patel
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| |
Collapse
|
8
|
Ren Y, Bertram AK, Patey GN. Influence of pH on Ice Nucleation by Kaolinite: Experiments and Molecular Simulations. J Phys Chem A 2022; 126:9227-9243. [DOI: 10.1021/acs.jpca.2c05323] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/02/2022]
Affiliation(s)
- Yi Ren
- Department of Chemistry, University of British Columbia, Vancouver, British Columbia, Canada V6T 1Z1
| | - Allan K. Bertram
- Department of Chemistry, University of British Columbia, Vancouver, British Columbia, Canada V6T 1Z1
| | - G. N. Patey
- Department of Chemistry, University of British Columbia, Vancouver, British Columbia, Canada V6T 1Z1
| |
Collapse
|
9
|
Keinert A, Deck K, Gädeke T, Leisner T, Kiselev AA. Mechanism of ice nucleation in liquid water on alkali feldspars. Faraday Discuss 2022; 235:148-161. [DOI: 10.1039/d1fd00115a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Crystallization of supercooled liquid water in most natural environments starts with heterogeneous nucleation of ice induced by a nucleation site. Mineral surfaces, which form the majority of aqueous interfaces in...
Collapse
|
10
|
Lata NN, Zhou J, Hamilton P, Larsen M, Sarupria S, Cantrell W. Multivalent Surface Cations Enhance Heterogeneous Freezing of Water on Muscovite Mica. J Phys Chem Lett 2020; 11:8682-8689. [PMID: 32955892 DOI: 10.1021/acs.jpclett.0c02121] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Heterogeneous ice nucleation is a crucial phenomenon in various fields of fundamental and applied science. We investigate the effect of surface cations on freezing of water on muscovite mica. Mica is unique in that the exposed ion on its surface can be readily and easily exchanged without affecting other properties such as surface roughness. We investigate freezing on natural (K+) mica and mica in which we have exchanged K+ for Al3+, Mg2+, Ca2+, and Sr2+. We find that liquid water freezes at higher temperatures when ions of higher valency are present on the surface, thus exposing more of the underlying silica layer. Our data also show that the size of the ion affects the characteristic freezing temperature. Using molecular dynamics simulations, we investigate the effects that the ion valency and exposed silica layer have on the behavior of water on the surface. The results indicate that multivalent cations enhance the probability of forming large clusters of hydrogen bonded water molecules that are anchored by the hydration shells of the cations. These clusters also have a large fraction of free water that can reorient to take ice-like configurations, which are promoted by the regions on mica devoid of the ions. Thus, these clusters could serve as seedbeds for ice nuclei. The combined experimental and simulation studies shed new light on the influence of surface ions on heterogeneous ice nucleation.
Collapse
Affiliation(s)
- Nurun Nahar Lata
- Atmospheric Sciences Program, Michigan Technological University, Houghton, Michigan 49931, United States
| | - Jiarun Zhou
- Department of Chemical and Biomolecular Engineering, Clemson University, Clemson, South Carolina 29634, United States
| | - Pearce Hamilton
- Department of Physics and Astronomy, College of Charleston, Charleston, South Carolina 29424, United States
| | - Michael Larsen
- Department of Physics and Astronomy, College of Charleston, Charleston, South Carolina 29424, United States
- Atmospheric Sciences Program and Department of Physics, Michigan Technological University, Houghton, Michigan 49931, United States
| | - Sapna Sarupria
- Department of Chemical and Biomolecular Engineering, Clemson University, Clemson, South Carolina 29634, United States
| | - Will Cantrell
- Atmospheric Sciences Program and Department of Physics, Michigan Technological University, Houghton, Michigan 49931, United States
| |
Collapse
|
11
|
Gallo P, Loerting T, Sciortino F. Supercooled water: A polymorphic liquid with a cornucopia of behaviors. J Chem Phys 2019; 151:210401. [DOI: 10.1063/1.5135706] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023] Open
Affiliation(s)
- Paola Gallo
- Department of Mathematics and Physics, University Roma Tre Via della Vasca Navale 84, I-00146 Rome, Italy
| | - Thomas Loerting
- Institute of Physical Chemistry, University of Innsbruck Innrain 52c, A-6020 Innsbruck, Austria
| | - Francesco Sciortino
- Department of Physics, Sapienza Universitá di Roma, Piazzale Aldo Moro 2, I-00185 Roma, Italy
| |
Collapse
|