1
|
Wang J, Fan L, Li L, Du Q, Jiao K. Ice Nucleation Mechanisms on Platinum Surfaces in PEM Fuel Cells: Effects of Surface Morphology and Wettability. ADVANCED SCIENCE (WEINHEIM, BADEN-WURTTEMBERG, GERMANY) 2024:e2406861. [PMID: 39116315 DOI: 10.1002/advs.202406861] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/20/2024] [Revised: 07/30/2024] [Indexed: 08/10/2024]
Abstract
Understanding the ice nucleation mechanism in the catalyst layers (CLs) of proton exchange membrane (PEM) fuel cells and inhibiting icing by designing the CLs can optimize the cold start strategies, which can enhance the performance of PEM fuel cells. Herein, mitigating the structural matching and templating effects by adjusting the surface morphology and wettability can inhibit icing on the platinum (Pt) catalyst surface effectively. The Pt(211) surface can inhibit icing because the atomic spacing of (211) crystalline surface is much larger than the characteristic distance of ice crystal, thereby mitigating the structural matching effects. A water overlayer on the Pt surface induced by the strong attraction of Pt can act as a template for ice layers and plays an important role in the icing process. Buckling of water overlayer due to the larger atomic spacing of (211) crystalline surface mitigates the templating effect and inhibits icing. Moreover, the water overlayer on the hydrophobic Pt(211) surface with fewer water molecules also mitigates the templating effect, which makes ice nucleation more difficult than homogeneous nucleation. These findings reveal the ice nucleation mechanisms on the Pt catalyst surface from the molecular level and are valuable for catalyst designs to inhibit icing in CL.
Collapse
Affiliation(s)
- Jiaqi Wang
- State Key Laboratory of Engines, Tianjin University, 135 Yaguan Road, Tianjin, 300350, China
- National Industry-Education Platform for Energy Storage, Tianjin University, 135 Yaguan Road, Tianjin, 300350, China
| | - Linhao Fan
- State Key Laboratory of Engines, Tianjin University, 135 Yaguan Road, Tianjin, 300350, China
- National Industry-Education Platform for Energy Storage, Tianjin University, 135 Yaguan Road, Tianjin, 300350, China
| | - Lincai Li
- State Key Laboratory of Engines, Tianjin University, 135 Yaguan Road, Tianjin, 300350, China
- National Industry-Education Platform for Energy Storage, Tianjin University, 135 Yaguan Road, Tianjin, 300350, China
| | - Qing Du
- State Key Laboratory of Engines, Tianjin University, 135 Yaguan Road, Tianjin, 300350, China
- National Industry-Education Platform for Energy Storage, Tianjin University, 135 Yaguan Road, Tianjin, 300350, China
| | - Kui Jiao
- State Key Laboratory of Engines, Tianjin University, 135 Yaguan Road, Tianjin, 300350, China
- National Industry-Education Platform for Energy Storage, Tianjin University, 135 Yaguan Road, Tianjin, 300350, China
| |
Collapse
|
2
|
Camarillo M, Oller-Iscar J, M Conde M, Ramírez J, Sanz E. Effect of substrate mismatch, orientation, and flexibility on heterogeneous ice nucleation. J Chem Phys 2024; 160:134505. [PMID: 38557847 DOI: 10.1063/5.0188929] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2023] [Accepted: 03/14/2024] [Indexed: 04/04/2024] Open
Abstract
Heterogeneous nucleation is the main path to ice formation on Earth. The ice nucleating ability of a certain substrate is mainly determined by both molecular interactions and the structural mismatch between the ice and the substrate lattices. We focus on the latter factor using molecular simulations of the mW model. Quantifying the effect of structural mismatch alone is challenging due to its coupling with molecular interactions. To disentangle both the factors, we use a substrate composed of water molecules in such a way that any variation on the nucleation temperature can be exclusively ascribed to the structural mismatch. We find that a 1% increase in structural mismatch leads to a decrease of ∼4 K in the nucleation temperature. We also analyze the effect of orientation of the substrate with respect to the liquid. The three main ice orientations (basal, primary prism, and secondary prism) have a similar ice nucleating ability. We finally assess the effect of lattice flexibility by comparing substrates where molecules are immobile to others where a certain freedom to fluctuate around the lattice positions is allowed. Interestingly, we find that the latter type of substrate is more efficient in nucleating ice because it can adapt its structure to that of ice.
Collapse
Affiliation(s)
- M Camarillo
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - J Oller-Iscar
- Department of Chemical Engineering, Universidad Politécnica de Madrid, José Gutiérrez Abascal 2, 28006 Madrid, Spain
| | - M M Conde
- Department of Chemical Engineering, Universidad Politécnica de Madrid, José Gutiérrez Abascal 2, 28006 Madrid, Spain
| | - J Ramírez
- Department of Chemical Engineering, Universidad Politécnica de Madrid, José Gutiérrez Abascal 2, 28006 Madrid, Spain
| | - E Sanz
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
3
|
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
|
4
|
Yuan T, DeFever RS, Zhou J, Cortes-Morales EC, Sarupria S. RSeeds: Rigid Seeding Method for Studying Heterogeneous Crystal Nucleation. J Phys Chem B 2023; 127:4112-4125. [PMID: 37130351 DOI: 10.1021/acs.jpcb.3c00910] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/04/2023]
Abstract
Heterogeneous nucleation is the dominant form of liquid-to-solid transition in nature. Although molecular simulations are most uniquely suited to studying nucleation, the waiting time to observe even a single nucleation event can easily exceed the current computational capabilities. Therefore, there exists an imminent need for methods that enable computationally fast and feasible studies of heterogeneous nucleation. Seeding is a technique that has proven to be successful at dramatically expanding the range of computationally accessible nucleation rates in simulation studies of homogeneous crystal nucleation. In this article, we introduce a new seeding method for heterogeneous nucleation called Rigid Seeding (RSeeds). Crystalline seeds are treated as pseudorigid bodies and simulated on a surface with metastable liquid above its melting temperature. This allows the seeds to adapt to the surface and identify favorable seed-surface configurations, which is necessary for reliable predictions of crystal polymorphs that form and the corresponding heterogeneous nucleation rates. We demonstrate and validate RSeeds for heterogeneous ice nucleation on a flexible self-assembled monolayer surface, a mineral surface based on kaolinite, and two model surfaces. RSeeds predicts the correct ice polymorph, exposed crystal plane, and rotation on the surface. RSeeds is semiquantitative and can be used to estimate the critical nucleus size and nucleation rate when combined with classical nucleation theory. We demonstrate that RSeeds can be used to evaluate nucleation rates spanning many orders of magnitude.
Collapse
Affiliation(s)
- Tianmu Yuan
- Department of Chemical and Biomolecular Engineering, Clemson University, Clemson, South Carolina 29634, United States
- Department of Chemical Engineering, The University of Manchester, Manchester, U.K. M13 9PL
| | - Ryan S DeFever
- Department of Chemical and Biomolecular Engineering, Clemson University, Clemson, South Carolina 29634, United States
| | - Jiarun Zhou
- Department of Chemical and Biomolecular Engineering, Clemson University, Clemson, South Carolina 29634, United States
| | | | - Sapna Sarupria
- Department of Chemistry, University of Minnesota Twin Cities, Minneapolis, Minnesota 55455, United States
| |
Collapse
|
5
|
Souda R, Aizawa T, Nagao T. Influences of 1-Propanol and Methanol Additives on Crystallization of Thin Amorphous Solid Water Films. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2022; 38:14422-14429. [PMID: 36383433 DOI: 10.1021/acs.langmuir.2c01978] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
In general, randomly oriented ice crystallites are formed by heating amorphous solid water (ASW) films at ∼160 K via homogeneous nucleation. Here, we demonstrate that monolayers of methanol and 1-propanol additives incorporated in the multilayer ASW film lead to heterogeneous nucleation at the substrate interface of Pt(111), as evidenced by the occurrence of epitaxial ice growth. The mobility of water in direct contact with the Pt(111) substrate is decreased relative to that in the bulk, but it can be increased via interactions with hydrophobic moieties of alcohols that are segregated to the interfacial region. As a result, heterogeneous nucleation occurs at ca. 160 K along with homogeneous nucleation in the film interior. However, the template effect is quenched when the alcohols are in direct contact with the substrate. The methanol adspecies deposited onto the ASW film surface induces heterogeneous nucleation at a temperature as low as 145 K, but the 1-propanol adspecies has no such an effect. Their different ability of heterogeneous nucleation at the free ASW film surface, as well as their uptake behaviors in the near surface region, is associated with the hydrophobic hydration of the alcohols resulting from different lengths of the aliphatic moiety.
Collapse
Affiliation(s)
- Ryutaro Souda
- Electron Microscopy Analysis Station, National Institute for Materials Science, 1-2-1 Sengen, Tsukuba, Ibaraki 305-0047, Japan
| | - Takashi Aizawa
- International Center for Materials Nanoarchitectonics, National Institute for Materials Science, 1-1 Namiki, Tsukuba, Ibaraki 305-0044, Japan
| | - Tadaaki Nagao
- International Center for Materials Nanoarchitectonics, National Institute for Materials Science, 1-1 Namiki, Tsukuba, Ibaraki 305-0044, Japan
| |
Collapse
|
6
|
Abstract
Crystal nucleation is one of the most fundamental processes in the physical sciences and almost always occurs heterogeneously with the aid of a nucleating substrate. No example of nucleation is more ubiquitous and impactful than the formation of ice, vital to fields as diverse as geology, biology, aeronautics, and climate science. However, despite considerable effort, we still cannot predict a priori the efficacy of a nucleating agent. Here we utilize deep learning methods to accurately predict nucleation ability from images of room temperature liquid water-generated from molecular dynamics simulations-on a broad range of substrates. The resulting model, named IcePic, can rapidly and accurately infer nucleation ability, eliminating the requirement for either notoriously expensive simulations or direct experimental measurement. In an online poll, IcePic was found to significantly outperform humans in predicting the ice nucleating efficacy of materials. By analyzing the typical errors made by humans, as well as the application of reverse interpretation methods, physical insights into the role the water contact layer plays in ice nucleation have been obtained. Moving forward, we suggest that IcePic can be used as an easy, cheap, and rapid way to discern the nucleation ability of substrates, also with potential for learning other properties related to interfacial water.
Collapse
|
7
|
Odendahl NL, Geissler PL. Local Ice-like Structure at the Liquid Water Surface. J Am Chem Soc 2022; 144:11178-11188. [PMID: 35696525 DOI: 10.1021/jacs.2c01827] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Experiments and computer simulations have established that liquid water's surfaces can deviate in important ways from familiar bulk behavior. Even in the simplest case of an air-water interface, distinctive layering, orientational biases, and hydrogen bond arrangements have been reported, but an overarching picture of their origins and relationships has been incomplete. Here we show that a broad set of such observations can be understood through an analogy with the basal face of crystalline ice. Using simulations, we demonstrate a number of structural similarities between water and ice surfaces, suggesting the presence of domains at the air-water interface with ice-like features that persist over 2-3 molecular diameters. Most prominent is a shared characteristic layering of molecular density and orientation perpendicular to the interface. Lateral correlations of hydrogen bond network geometry point to structural similarities in the parallel direction as well. Our results bolster and significantly extend previous conceptions of ice-like structure at the liquid's boundary and suggest that the much-discussed quasi-liquid layer on ice evolves subtly above the melting point into a quasi-ice layer at the surface of liquid water.
Collapse
Affiliation(s)
- Nathan L Odendahl
- Department of Chemistry, University of California, Berkeley, California 94720, United States.,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Phillip L Geissler
- Department of Chemistry, University of California, Berkeley, California 94720, United States.,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
8
|
Hussain S, Haji-Akbari A. How to quantify and avoid finite size effects in computational studies of crystal nucleation: The case of homogeneous crystal nucleation. J Chem Phys 2022; 156:054503. [DOI: 10.1063/5.0079702] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Affiliation(s)
- Sarwar Hussain
- Department of Chemical and Environmental Engineering, Yale University, New Haven, Connecticut 06520, USA
| | - Amir Haji-Akbari
- Department of Chemical and Environmental Engineering, Yale University, New Haven, Connecticut 06520, USA
| |
Collapse
|
9
|
Patterning Configuration of Surface Hydrophilicity by Graphene Nanosheet towards the Inhibition of Ice Nucleation and Growth. COATINGS 2022. [DOI: 10.3390/coatings12010052] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Abstract
Freezing of liquid water occurs in many natural phenomena and affects countless human activities. The freezing process mainly involves ice nucleation and continuous growth, which are determined by the energy and structure fluctuation in supercooled water. Herein, considering the surface hydrophilicity and crystal structure differences between metal and graphene, we proposed a kind of surface configuration design, which was realized by graphene nanosheets being alternately anchored on a metal substrate. Ice nucleation and growth were investigated by molecular dynamics simulations. The surface configuration could induce ice nucleation to occur preferentially on the metal substrate where the surface hydrophilicity was higher than the lateral graphene nanosheet. However, ice nucleation could be delayed to a certain extent under the hindering effect of the interfacial water layer formed by the high surface hydrophilicity of the metal substrate. Furthermore, the graphene nanosheets restricted lateral expansion of the ice nucleus at the clearance, leading to the formation of a curved surface of the ice nucleus as it grew. As a result, ice growth was suppressed effectively due to the Gibbs–Thomson effect, and the growth rate decreased by 71.08% compared to the pure metal surface. Meanwhile, boundary misorientation between ice crystals was an important issue, which also prejudiced the growth of the ice crystal. The present results reveal the microscopic details of ice nucleation and growth inhibition of the special surface configuration and provide guidelines for the rational design of an anti-icing surface.
Collapse
|
10
|
Nikiforidis VM, Datta S, Borg MK, Pillai R. Impact of surface nanostructure and wettability on interfacial ice physics. J Chem Phys 2021; 155:234307. [PMID: 34937379 DOI: 10.1063/5.0069896] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Ice accumulation on solid surfaces is a severe problem for safety and functioning of a large variety of engineering systems, and its control is an enormous challenge that influences the safety and reliability of many technological applications. The use of molecular dynamics (MD) simulations is popular, but as ice nucleation is a rare event when compared to simulation timescales, the simulations need to be accelerated to force ice to form on a surface, which affects the accuracy and/or applicability of the results obtained. Here, we present an alternative seeded MD simulation approach, which reduces the computational cost while still ensuring accurate simulations of ice growth on surfaces. In addition, this approach enables, for the first time, brute-force all-atom water simulations of ice growth on surfaces unfavorable for nucleation within MD timescales. Using this approach, we investigate the effect of surface wettability and structure on ice growth in the crucial surface-ice interfacial region. Our main findings are that the surface structure can induce a flat or buckled overlayer to form within the liquid, and this transition is mediated by surface wettability. The first overlayer and the bulk ice compete to structure the intermediate water layers between them, the relative influence of which is traced using density heat maps and diffusivity measurements. This work provides new understanding on the role of the surface properties on the structure and dynamics of ice growth, and we also present a useful framework for future research on surface icing simulations.
Collapse
Affiliation(s)
- Vasileios-Martin Nikiforidis
- School of Engineering, Institute for Multiscale Thermofluids, The University of Edinburgh, Edinburgh EH9 3FB, United Kingdom
| | - Saikat Datta
- School of Engineering, Institute for Multiscale Thermofluids, The University of Edinburgh, Edinburgh EH9 3FB, United Kingdom
| | - Matthew K Borg
- School of Engineering, Institute for Multiscale Thermofluids, The University of Edinburgh, Edinburgh EH9 3FB, United Kingdom
| | - Rohit Pillai
- School of Engineering, Institute for Multiscale Thermofluids, The University of Edinburgh, Edinburgh EH9 3FB, United Kingdom
| |
Collapse
|
11
|
Lu H, Xu Q, Wu J, Hong R, Zhang Z. Effect of interfacial dipole on heterogeneous ice nucleation. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2021; 33:375001. [PMID: 34181589 DOI: 10.1088/1361-648x/ac0f2c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/02/2021] [Accepted: 06/28/2021] [Indexed: 06/13/2023]
Abstract
In this work, we performed molecular dynamics simulations of ice nucleation on a rigid surface model of cubic zinc blende structure with different surface dipole strength and orientation. Our results show that, although substrates are excellently lattice-matched to cubic ice, ice nucleation merely occurred as the interfacial water molecules (IWs) show identical or similar orientations to that of water molecules in cubic ice. Free energy landscapes revealed that, as substrates have non-suitable dipole strength/orientation, there exist large free energy barriers for rotating dipole IWs to the right orientation to trigger ice formation. This study stresses that, beyond the traditional view of lattice match and the similarity of lattice length between the substrate and new-formed crystal, the similarity between molecular orientations of interfacial component and component in the specific new-formed crystalline face is also critical for promoting ice nucleation.
Collapse
Affiliation(s)
- Hao Lu
- Department of Physics, Research Institute for Biomimetics and Soft Matter, Fujian Provincial Key Laboratory for Soft Functional Materials Research, Jiujiang Research Institute, Xiamen University, Xiamen, Fujian 361005, People's Republic of China
| | - Quanming Xu
- Department of Physics, Research Institute for Biomimetics and Soft Matter, Fujian Provincial Key Laboratory for Soft Functional Materials Research, Jiujiang Research Institute, Xiamen University, Xiamen, Fujian 361005, People's Republic of China
| | - Jianyang Wu
- Department of Physics, Research Institute for Biomimetics and Soft Matter, Fujian Provincial Key Laboratory for Soft Functional Materials Research, Jiujiang Research Institute, Xiamen University, Xiamen, Fujian 361005, People's Republic of China
| | - Rongdun Hong
- Department of Physics, Research Institute for Biomimetics and Soft Matter, Fujian Provincial Key Laboratory for Soft Functional Materials Research, Jiujiang Research Institute, Xiamen University, Xiamen, Fujian 361005, People's Republic of China
| | - Zhisen Zhang
- Department of Physics, Research Institute for Biomimetics and Soft Matter, Fujian Provincial Key Laboratory for Soft Functional Materials Research, Jiujiang Research Institute, Xiamen University, Xiamen, Fujian 361005, People's Republic of China
| |
Collapse
|
12
|
Metya AK, Molinero V. Is Ice Nucleation by Organic Crystals Nonclassical? An Assessment of the Monolayer Hypothesis of Ice Nucleation. J Am Chem Soc 2021; 143:4607-4624. [PMID: 33729789 DOI: 10.1021/jacs.0c12012] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Potent ice nucleating organic crystals display an increase in nucleation efficiency with pressure and memory effect after pressurization that set them apart from inorganic nucleants. These characteristics were proposed to arise from an ordered water monolayer at the organic-water interface. It was interpreted that ordering of the monolayer is the limiting step for ice nucleation on organic crystals, rendering their mechanism of nucleation nonclassical. Despite the importance of organics in atmospheric ice nucleation, that explanation has never been investigated. Here we elucidate the structure of interfacial water and its role in ice nucleation at ambient pressure on phloroglucinol dihydrate, the paradigmatic example of outstanding ice nucleating organic crystal, using molecular simulations. The simulations confirm the existence of an interfacial monolayer that orders on cooling and becomes fully ordered upon ice formation. The monolayer does not resemble any ice face but seamlessly connects the distinct hydrogen-bonding orders of ice and the organic surface. Although large ordered patches develop in the monolayer before ice nucleates, we find that the critical step is the formation of the ice crystallite, indicating that the mechanism is classical. We predict that the fully ordered, crystalline monolayer nucleates ice above -2 °C and could be responsible for the exceptional ice nucleation by the organic crystal at high pressures. The lifetime of the fully ordered monolayer around 0 °C, however, is too short to account for the memory effect reported in the experiments. The latter could arise from an increase in the melting temperature of ice confined by strongly ice-binding surfaces.
Collapse
Affiliation(s)
- Atanu K Metya
- Department of Chemistry, The University of Utah, Salt Lake City, Utah 84112-0850, United States
| | - Valeria Molinero
- Department of Chemistry, The University of Utah, Salt Lake City, Utah 84112-0850, United States
| |
Collapse
|
13
|
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
|
14
|
Hussain S, Haji-Akbari A. How to quantify and avoid finite size effects in computational studies of crystal nucleation: The case of heterogeneous ice nucleation. J Chem Phys 2021; 154:014108. [DOI: 10.1063/5.0026355] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Affiliation(s)
- Sarwar Hussain
- Department of Chemical and Environmental Engineering, Yale University, New Haven, Connecticut 06520, USA
| | - Amir Haji-Akbari
- Department of Chemical and Environmental Engineering, Yale University, New Haven, Connecticut 06520, USA
| |
Collapse
|
15
|
Souda R, Aizawa T, Sugiyama N, Takeguchi M. Nucleation and growth of water ice on oxide surfaces: the influence of a precursor to water dissociation. Phys Chem Chem Phys 2020; 22:20515-20523. [PMID: 32966413 DOI: 10.1039/d0cp01897j] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
We have investigated how nucleation and growth processes of ice are influenced by interfacial molecular interactions on some oxide surfaces, such as rutile TiO2(110), TiO2(100), MgO(100), and Al2O3(0001), based on the diffraction patterns of electrons transmitted through ice crystallites under the experimental configuration of reflection high energy electron diffraction (RHEED). The cubic ice Ic grows on the TiO2(110) surface with the epitaxial relationship of (110)Ic//(110)TiO2 and [001]Ic//[11[combining macron]0]TiO2. The epitaxial ice growth tends to be disturbed on the TiO2(110) surface under the presence of oxygen vacancies and adatoms. The result is not simply ascribable to small misfit values between TiO2 and ice Ic lattices (∼2%) because ice grains are formed randomly on TiO2(100). No template effects are identified during ice nucleation on the pristine MgO(100) and Al2O3(0001) surfaces either. The water molecules are chemisorbed weakly on these surfaces as a precursor to dissociation via the acid-base interaction. Such anchored water species act as an inhibitor of epitaxial ice growth because the orientation flexibility of physisorbed water during nucleation is hampered at the interface by the preferential formation of hydrogen bonds.
Collapse
Affiliation(s)
- Ryutaro Souda
- Transmission Electron Microscopy Station, National Institute for Materials Science, 1-2-1 Sengen, Tsukuba, Ibaraki 305-0047, Japan.
| | - Takashi Aizawa
- Center for Functional Sensor & Actuator, National Institute for Materials Science, 1-1 Namiki, Tsukuba, Ibaraki 305-0044, Japan
| | - Naoyuki Sugiyama
- Transmission Electron Microscopy Station, National Institute for Materials Science, 1-2-1 Sengen, Tsukuba, Ibaraki 305-0047, Japan. and Toray Research Center Inc., 3-3-7 Sonoyama, Otsu, Shiga 520-8567, Japan
| | - Masaki Takeguchi
- Transmission Electron Microscopy Station, National Institute for Materials Science, 1-2-1 Sengen, Tsukuba, Ibaraki 305-0047, Japan.
| |
Collapse
|
16
|
Predicting heterogeneous ice nucleation with a data-driven approach. Nat Commun 2020; 11:4777. [PMID: 32963232 PMCID: PMC7509812 DOI: 10.1038/s41467-020-18605-3] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/06/2020] [Accepted: 08/28/2020] [Indexed: 01/05/2023] Open
Abstract
Water in nature predominantly freezes with the help of foreign materials through a process known as heterogeneous ice nucleation. Although this effect was exploited more than seven decades ago in Vonnegut's pioneering cloud seeding experiments, it remains unclear what makes a material a good ice former. Here, we show through a machine learning analysis of nucleation simulations on a database of diverse model substrates that a set of physical descriptors for heterogeneous ice nucleation can be identified. Our results reveal that, beyond Vonnegut's connection with the lattice match to ice, three new microscopic factors help to predict the ice nucleating ability. These are: local ordering induced in liquid water, density reduction of liquid water near the surface and corrugation of the adsorption energy landscape felt by water. With this we take a step towards quantitative understanding of heterogeneous ice nucleation and the in silico design of materials to control ice formation.
Collapse
|
17
|
Shao M, Zhang C, Qi C, Wang C, Wang J, Ye F, Zhou X. Hydrogen polarity of interfacial water regulates heterogeneous ice nucleation. Phys Chem Chem Phys 2019; 22:258-264. [PMID: 31808477 DOI: 10.1039/c9cp04867g] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Using all-atomic molecular dynamics (MD) simulations, we show that the structure of interfacial water (IW) induced by substrates characterizes the ability of a substrate to nucleate ice. We probe the shape and structure of ice nuclei and the corresponding supercooling temperatures to measure the ability of IW with various hydrogen polarities for ice nucleation, and find that the hydrogen polarization of IW even with the ice-like oxygen lattice increases the contact angle of the ice nucleus on IW, thus lifting the free energy barrier of heterogeneous ice nucleation. The results show that not only the oxygen lattice order but the hydrogen disorder of IW on substrates are required to effectively facilitate the freezing of top water.
Collapse
Affiliation(s)
- Mingzhe Shao
- College of Light Industry Science and Engineering, Tianjin University of Science and Technology, Tianjin, China
| | | | | | | | | | | | | |
Collapse
|
18
|
Soni A, Patey GN. Simulations of water structure and the possibility of ice nucleation on selected crystal planes of K-feldspar. J Chem Phys 2019; 150:214501. [DOI: 10.1063/1.5094645] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
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
|
19
|
Cao B, Xu E, Li T. Anomalous Stability of Two-Dimensional Ice Confined in Hydrophobic Nanopores. ACS NANO 2019; 13:4712-4719. [PMID: 30892864 DOI: 10.1021/acsnano.9b01014] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
The freezing of water mostly proceeds via heterogeneous ice nucleation, a process in which an effective nucleation medium not only expedites ice crystallization but also may effectively direct the polymorph selection of ice. Here, we show that water confined within a hydrophobic slit nanopore exhibits a freezing behavior strongly distinguished from its bulk counterpart. Such a difference is reflected by a strong, non-monotonic pore-size dependence of freezing temperature but, more surprisingly, by an unexpected stacking ordering of crystallized two-dimensional ice containing just a few ice layers. In particular, confined trilayer ice is found to exclusively crystallize into a well-ordered, hexagonal stacking sequence despite the fact that nanopore exerts no explicit constraint on stacking order. The absence of cubic stacking sequence is found to be originated from the intrinsically lower thermodynamic stability of cubic ice over hexagonal ice at the interface, which contrasts sharply the nearly degenerated stability of bulk hexagonal and cubic ices. Detailed examination clearly reveals that the divergence is attributed to the inherent difference between the two ice polymorphs in their surface phonon modes, which is further found to generically occur at both hydrophobic and hydrophilic surfaces.
Collapse
Affiliation(s)
- Boxiao Cao
- Department of Civil and Environmental Engineering , George Washington University , Washington , D.C. 20052 , United States
| | - Enshi Xu
- Department of Civil and Environmental Engineering , George Washington University , Washington , D.C. 20052 , United States
| | - Tianshu Li
- Department of Civil and Environmental Engineering , George Washington University , Washington , D.C. 20052 , United States
| |
Collapse
|
20
|
Wu S, He Z, Zang J, Jin S, Wang Z, Wang J, Yao Y, Wang J. Heterogeneous ice nucleation correlates with bulk-like interfacial water. SCIENCE ADVANCES 2019; 5:eaat9825. [PMID: 30993196 PMCID: PMC6461451 DOI: 10.1126/sciadv.aat9825] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/24/2018] [Accepted: 02/14/2019] [Indexed: 05/31/2023]
Abstract
Establishing a direct correlation between interfacial water and heterogeneous ice nucleation (HIN) is essential for understanding the mechanism of ice nucleation. Here, we study the HIN efficiency on polyvinyl alcohol (PVA) surfaces with different densities of hydroxyl groups. We find that the HIN efficiency increases with the decreasing hydroxyl group density. By explicitly considering that interfacial water molecules of PVA films consist of "tightly bound water," "bound water," and "bulk-like water," we reveal that bulk-like water can be correlated directly to the HIN efficiency of surfaces. As the density of hydroxyl groups decreases, bulk-like water molecules can rearrange themselves with a reduced energy barrier into ice due to the diminishing constraint by the hydroxyl groups on the PVA surface. Our study not only provides a new strategy for experimentally controlling the HIN efficiency but also gives another perspective in understanding the mechanism of ice nucleation.
Collapse
Affiliation(s)
- Shuwang Wu
- Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
| | - Zhiyuan He
- Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
| | - Jinger Zang
- Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
| | - Shenglin Jin
- Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
| | - Zuowei Wang
- School of Mathematical, Physical and Computational Sciences, University of Reading, Whiteknights, Reading RG6 6AX, UK
| | - Jianping Wang
- Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
| | - Yefeng Yao
- Physics Department and Shanghai Key Laboratory of Magnetic Resonance, East China Normal University, North Zhongshan Road 3663, Shanghai 200062, China
| | - Jianjun Wang
- Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
| |
Collapse
|
21
|
|
22
|
Sayer T, Cox SJ. Stabilization of AgI's polar surfaces by the aqueous environment, and its implications for ice formation. Phys Chem Chem Phys 2019; 21:14546-14555. [DOI: 10.1039/c9cp02193k] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
AgI is a potent inorganic ice nucleating particle, a feature often attributed to the lattice match between its {0001} surfaces and ice. Dissolved ions are found to be essential to the stability of these polar surfaces, and crucial to ice formation.
Collapse
Affiliation(s)
- Thomas Sayer
- Department of Chemistry
- University of Cambridge
- Cambridge CB2 1EW
- UK
| | - Stephen J. Cox
- Department of Chemistry
- University of Cambridge
- Cambridge CB2 1EW
- UK
| |
Collapse
|
23
|
Pedevilla P, Fitzner M, Sosso GC, Michaelides A. Heterogeneous seeded molecular dynamics as a tool to probe the ice nucleating ability of crystalline surfaces. J Chem Phys 2018; 149:072327. [PMID: 30134662 DOI: 10.1063/1.5029336] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Ice nucleation plays a significant role in a large number of natural and technological processes, but it is challenging to investigate experimentally because of the small time scales (ns) and short length scales (nm) involved. On the other hand, conventional molecular simulations struggle to cope with the relatively long time scale required for critical ice nuclei to form. One way to tackle this issue is to take advantage of free energy or path sampling techniques. Unfortunately, these are computationally costly. Seeded molecular dynamics is a much less demanding alternative that has been successfully applied already to study the homogeneous freezing of water. However, in the case of heterogeneous ice nucleation, nature's favourite route to form ice, an array of suitable interfaces between the ice seeds and the substrate of interest has to be built, and this is no trivial task. In this paper, we present a Heterogeneous SEEDing (HSEED) approach which harnesses a random structure search framework to tackle the ice-substrate challenge, thus enabling seeded molecular dynamics simulations of heterogeneous ice nucleation on crystalline surfaces. We validate the HSEED framework by investigating the nucleation of ice on (i) model crystalline surfaces, using the coarse-grained mW model, and (ii) cholesterol crystals, employing the fully atomistic TIP4P/ice water model. We show that the HSEED technique yields results in excellent agreement with both metadynamics and forward flux sampling simulations. Because of its computational efficiency, the HSEED method allows one to rapidly assess the ice nucleation ability of whole libraries of crystalline substrates-a long-awaited computational development in, e.g., atmospheric science.
Collapse
Affiliation(s)
- Philipp Pedevilla
- Thomas Young Centre, London Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - Martin Fitzner
- Thomas Young Centre, London Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - Gabriele C Sosso
- Department of Chemistry and Centre for Scientific Computing, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - Angelos Michaelides
- Thomas Young Centre, London Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| |
Collapse
|
24
|
|
25
|
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
|
26
|
Jin S, Liu J, Lv J, Wu S, Wangs J. Interfacial Materials for Anti-Icing: Beyond Superhydrophobic Surfaces. Chem Asian J 2018. [DOI: 10.1002/asia.201800241] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Shenglin Jin
- Key Laboratory of Green Printing; Institute of Chemistry; Chinese Academy of Sciences; Zhongguancun North First Street 2 100190 Beijing P. R. China
| | - Jie Liu
- Key Laboratory of Green Printing; Institute of Chemistry; Chinese Academy of Sciences; Zhongguancun North First Street 2 100190 Beijing P. R. China
- Max-Planck Institute for Polymer Research.; Ackermannweg 10 55128 Mainz Germany
| | - Jianyong Lv
- Key Laboratory of Green Printing; Institute of Chemistry; Chinese Academy of Sciences; Zhongguancun North First Street 2 100190 Beijing P. R. China
| | - Shuwang Wu
- Key Laboratory of Green Printing; Institute of Chemistry; Chinese Academy of Sciences; Zhongguancun North First Street 2 100190 Beijing P. R. China
| | - Jianjun Wangs
- Key Laboratory of Green Printing; Institute of Chemistry; Chinese Academy of Sciences; Zhongguancun North First Street 2 100190 Beijing P. R. China
| |
Collapse
|
27
|
Shevkunov SV. Water Structure in the Contact Layer on the Surface of Crystalline Silver Iodine. J STRUCT CHEM+ 2018. [DOI: 10.1134/s0022476618030137] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
28
|
Shevkunov SV. The Effect of Temperature on Nucleation of Condensed Water Phase on the Surface of a β-AgI Crystal. 2. Formation Work. COLLOID JOURNAL 2018. [DOI: 10.1134/s1061933x18020102] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
29
|
Shevkunov SV. The Effect of Temperature on Nucleation of Condensed Water Phase on the Surface of a β-AgI Crystal. 1. Structure. COLLOID JOURNAL 2018. [DOI: 10.1134/s1061933x18020096] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
30
|
Cox SJ, Taylor DJF, Youngs TGA, Soper AK, Totton TS, Chapman RG, Arjmandi M, Hodges MG, Skipper NT, Michaelides A. Formation of Methane Hydrate in the Presence of Natural and Synthetic Nanoparticles. J Am Chem Soc 2018; 140:3277-3284. [PMID: 29401390 PMCID: PMC5860788 DOI: 10.1021/jacs.7b12050] [Citation(s) in RCA: 26] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
Natural gas hydrates occur widely
on the ocean-bed and in permafrost
regions, and have potential as an untapped energy resource. Their
formation and growth, however, poses major problems for the energy
sector due to their tendency to block oil and gas pipelines, whereas
their melting is viewed as a potential contributor to climate change.
Although recent advances have been made in understanding bulk methane
hydrate formation, the effect of impurity particles, which are always
present under conditions relevant to industry and the environment,
remains an open question. Here we present results from neutron scattering
experiments and molecular dynamics simulations that show that the
formation of methane hydrate is insensitive to the addition of a wide
range of impurity particles. Our analysis shows that this is due to
the different chemical natures of methane and water, with methane
generally excluded from the volume surrounding the nanoparticles.
This has important consequences for our understanding of the mechanism
of hydrate nucleation and the design of new inhibitor molecules.
Collapse
Affiliation(s)
- Stephen J Cox
- Department of Chemistry , University College London , 20 Gordon Street , London WC1H 0AJ , United Kingdom.,Thomas Young Centre and London Centre for Nanotechnology , 17-19 Gordon Street , London WC1H 0AH , United Kingdom
| | - Diana J F Taylor
- Thomas Young Centre and London Centre for Nanotechnology , 17-19 Gordon Street , London WC1H 0AH , United Kingdom.,Department of Physics and Astronomy , University College London , Gower Street , London WC1E 6BT , United Kingdom
| | - Tristan G A Youngs
- ISIS Facility , STFC Rutherford Appleton Laboratory , Harwell Oxford , Didcot OX11 0QX , United Kingdom
| | - Alan K Soper
- ISIS Facility , STFC Rutherford Appleton Laboratory , Harwell Oxford , Didcot OX11 0QX , United Kingdom
| | - Tim S Totton
- BP Exploration Operating Co. Ltd , Chertsey Road , Sunbury-on-Thames TW16 7LN , United Kingdom
| | - Richard G Chapman
- BP Exploration Operating Co. Ltd , Chertsey Road , Sunbury-on-Thames TW16 7LN , United Kingdom
| | - Mosayyeb Arjmandi
- BP Exploration Operating Co. Ltd , Chertsey Road , Sunbury-on-Thames TW16 7LN , United Kingdom
| | - Michael G Hodges
- BP Exploration Operating Co. Ltd , Chertsey Road , Sunbury-on-Thames TW16 7LN , United Kingdom
| | - Neal T Skipper
- Thomas Young Centre and London Centre for Nanotechnology , 17-19 Gordon Street , London WC1H 0AH , United Kingdom.,Department of Physics and Astronomy , University College London , Gower Street , London WC1E 6BT , United Kingdom
| | - Angelos Michaelides
- Thomas Young Centre and London Centre for Nanotechnology , 17-19 Gordon Street , London WC1H 0AH , United Kingdom.,Department of Physics and Astronomy , University College London , Gower Street , London WC1E 6BT , United Kingdom
| |
Collapse
|
31
|
Qiu Y, Lupi L, Molinero V. Is Water at the Graphite Interface Vapor-like or Ice-like? J Phys Chem B 2018; 122:3626-3634. [PMID: 29298058 DOI: 10.1021/acs.jpcb.7b11476] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Graphitic surfaces are the main component of soot, a major constituent of atmospheric aerosols. Experiments indicate that soots of different origins display a wide range of abilities to heterogeneously nucleate ice. The ability of pure graphite to nucleate ice in experiments, however, seems to be almost negligible. Nevertheless, molecular simulations with the monatomic water model mW with water-carbon interactions parameterized to reproduce the experimental contact angle of water on graphite predict that pure graphite nucleates ice. According to classical nucleation theory, the ability of a surface to nucleate ice is controlled by the binding free energy between ice immersed in liquid water and the surface. To establish whether the discrepancy in freezing efficiencies of graphite in mW simulations and experiments arises from the coarse resolution of the model or can be fixed by reparameterization, it is important to elucidate the contributions of the water-graphite, water-ice, and ice-water interfaces to the free energy, enthalpy, and entropy of binding for both water and the model. Here we use thermodynamic analysis and free energy calculations to determine these interfacial properties. We demonstrate that liquid water at the graphite interface is not ice-like or vapor-like: it has similar free energy, entropy, and enthalpy as water in the bulk. The thermodynamics of the water-graphite interface is well reproduced by the mW model. We find that the entropy of binding between graphite and ice is positive and dominated, in both experiments and simulations, by the favorable entropy of reducing the ice-water interface. Our analysis indicates that the discrepancy in freezing efficiencies of graphite in experiments and the simulations with mW arises from the inability of the model to simultaneously reproduce the contact angle of liquid water on graphite and the free energy of the ice-graphite interface. This transferability issue is intrinsic to the resolution of the model, and arises from its lack of rotational degrees of freedom.
Collapse
Affiliation(s)
- Yuqing Qiu
- Department of Chemistry , The University of Utah , 315 South 1400 East , Salt Lake City , Utah 84112-0850 , United States
| | - Laura Lupi
- Department of Chemistry , The University of Utah , 315 South 1400 East , Salt Lake City , Utah 84112-0850 , United States
| | - Valeria Molinero
- Department of Chemistry , The University of Utah , 315 South 1400 East , Salt Lake City , Utah 84112-0850 , United States
| |
Collapse
|
32
|
Glatz B, Sarupria S. Heterogeneous Ice Nucleation: Interplay of Surface Properties and Their Impact on Water Orientations. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2018; 34:1190-1198. [PMID: 29020452 DOI: 10.1021/acs.langmuir.7b02859] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Ice is ubiquitous in nature, and heterogeneous ice nucleation is the most common pathway of ice formation. How surface properties affect the propensity to observe ice nucleation on that surface remains an open question. We present results of molecular dynamics studies of heterogeneous ice nucleation on model surfaces. The models surfaces considered emulate the chemistry of kaolinite, an abundant component of mineral dust. We investigate the interplay of surface lattice and hydrogen bonding properties in affecting ice nucleation. We find that lattice matching and hydrogen bonding are necessary but not sufficient conditions for observing ice nucleation at these surfaces. We correlate this behavior to the orientations sampled by the metastable supercooled water in contact with the surfaces. We find that ice is observed in cases where water molecules not only sample orientations favorable for bilayer formation but also do not sample unfavorable orientations. This distribution depends on both surface-water and water-water interactions and can change with subtle modifications to the surface properties. Our results provide insights into the diverse behavior of ice nucleation observed at different surfaces and highlight the complexity in elucidating heterogeneous ice nucleation.
Collapse
Affiliation(s)
- Brittany Glatz
- Department of Chemical & Biomolecular Engineering, Clemson University , Clemson, South Carolina 29634, United States
| | - Sapna Sarupria
- Department of Chemical & Biomolecular Engineering, Clemson University , Clemson, South Carolina 29634, United States
| |
Collapse
|
33
|
Lanaro G, Patey GN. The influence of ion hydration on nucleation and growth of LiF crystals in aqueous solution. J Chem Phys 2018; 148:024507. [DOI: 10.1063/1.5001521] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Affiliation(s)
- G. Lanaro
- 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
|
34
|
Zubeltzu J, Artacho E. Simulations of water nano-confined between corrugated planes. J Chem Phys 2017; 147:194509. [PMID: 29166107 DOI: 10.1063/1.5011468] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023] Open
Abstract
Water confined to nanoscale widths in two dimensions between ideal planar walls has been the subject of ample study, aiming at understanding the intrinsic response of water to confinement, avoiding the consideration of the chemistry of actual confining materials. In this work, we study the response of such nanoconfined water to the imposition of a periodicity in the confinement by means of computer simulations, both using empirical potentials and from first-principles. For that we propose a periodic confining potential emulating the atomistic oscillation of the confining walls, which allows varying the lattice parameter and amplitude of the oscillation. We do it for a triangular lattice, with several values of the lattice parameter: one which is ideal for commensuration with layers of Ih ice and other values that would correspond to more realistic substrates. For the former, the phase diagram shows an overall rise of the melting temperature. The liquid maintains a bi-layer triangular structure, however, despite the fact that it is not favoured by the external periodicity. The first-principles liquid is significantly affected by the modulation in its layering and stacking even at relatively small amplitudes of the confinement modulation. Beyond some critical modulation amplitude, the hexatic phase present in flat confinement is replaced by a trilayer crystalline phase unlike any of the phases encountered for flat confinement. For more realistic lattice parameters, the liquid does not display higher tendency to freeze, but it clearly shows inhomogeneous behaviour as the strength of the rugosity increases. In spite of this expected inhomogeneity, the structural and dynamical response of the liquid is surprisingly insensitive to the external modulation. Although the first-principles calculations give a more triangular liquid than the one observed with empirical potentials (TIP4P/2005), both agree remarkably well for the main conclusions of the study.
Collapse
Affiliation(s)
- Jon Zubeltzu
- CIC nanoGUNE, 20018 Donostia-San Sebastián, Spain
| | | |
Collapse
|
35
|
Fitzner M, Joly L, Ma M, Sosso GC, Zen A, Michaelides A. Communication: Truncated non-bonded potentials can yield unphysical behavior in molecular dynamics simulations of interfaces. J Chem Phys 2017; 147:121102. [DOI: 10.1063/1.4997698] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Martin Fitzner
- Thomas Young Centre, London Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - Laurent Joly
- Univ. Lyon, Université Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 Villeurbanne,
France
| | - Ming Ma
- Department of Mechanical Engineering, State Key Laboratory of Tribology and Center for Nano and Micro Mechanics, Tsinghua University, Beijing 100084, China
| | - Gabriele C. Sosso
- Department of Chemistry and Centre for Scientific Computing, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL,
United Kingdom
| | - Andrea Zen
- Thomas Young Centre, London Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - Angelos Michaelides
- Thomas Young Centre, London Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| |
Collapse
|
36
|
Shevkunov SV. Water-vapor clustering on the surface of β-AgI crystal in the field of defects with a disordered structure. COLLOID JOURNAL 2017. [DOI: 10.1134/s1061933x1705012x] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
37
|
Balance between hydration enthalpy and entropy is important for ice binding surfaces in Antifreeze Proteins. Sci Rep 2017; 7:11901. [PMID: 28928396 PMCID: PMC5605524 DOI: 10.1038/s41598-017-11982-8] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2017] [Accepted: 08/29/2017] [Indexed: 11/21/2022] Open
Abstract
Antifreeze Proteins (AFPs) inhibit the growth of an ice crystal by binding to it. The detailed binding mechanism is, however, still not fully understood. We investigated three AFPs using Molecular Dynamics simulations in combination with Grid Inhomogeneous Solvation Theory, exploring their hydration thermodynamics. The observed enthalpic and entropic differences between the ice-binding sites and the inactive surface reveal key properties essential for proteins in order to bind ice: While entropic contributions are similar for all sites, the enthalpic gain for all ice-binding sites is lower than for the rest of the protein surface. In contrast to most of the recently published studies, our analyses show that enthalpic interactions are as important as an ice-like pre-ordering. Based on these observations, we propose a new, thermodynamically more refined mechanism of the ice recognition process showing that the appropriate balance between entropy and enthalpy facilitates ice-binding of proteins. Especially, high enthalpic interactions between the protein surface and water can hinder the ice-binding activity.
Collapse
|
38
|
Perspectives on the Future of Ice Nucleation Research: Research Needs and Unanswered Questions Identified from Two International Workshops. ATMOSPHERE 2017. [DOI: 10.3390/atmos8080138] [Citation(s) in RCA: 49] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
39
|
Enhanced heterogeneous ice nucleation by special surface geometry. Nat Commun 2017; 8:15372. [PMID: 28513603 PMCID: PMC5442314 DOI: 10.1038/ncomms15372] [Citation(s) in RCA: 74] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2016] [Accepted: 03/23/2017] [Indexed: 12/01/2022] Open
Abstract
The freezing of water typically proceeds through impurity-mediated heterogeneous nucleation. Although non-planar geometry generically exists on the surfaces of ice nucleation centres, its role in nucleation remains poorly understood. Here we show that an atomically sharp, concave wedge can further promote ice nucleation with special wedge geometries. Our molecular analysis shows that significant enhancements of ice nucleation can emerge both when the geometry of a wedge matches the ice lattice and when such lattice match does not exist. In particular, a 45° wedge is found to greatly enhance ice nucleation by facilitating the formation of special topological defects that consequently catalyse the growth of regular ice. Our study not only highlights the active role of defects in nucleation but also suggests that the traditional concept of lattice match between a nucleation centre and crystalline lattice should be extended to include a broader match with metastable, non-crystalline structural motifs. Understanding ice nucleation is important for the development of accurate cloud models. Here Bi et al. show that sharp wedges can enhance ice nucleation both when the wedge geometry matches the ice lattice and when such matching is absent, in which case nucleation is promoted by topological defects.
Collapse
|
40
|
Strength of Alkane–Fluid Attraction Determines the Interfacial Orientation of Liquid Alkanes and Their Crystallization through Heterogeneous or Homogeneous Mechanisms. CRYSTALS 2017. [DOI: 10.3390/cryst7030086] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
41
|
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
|
42
|
Qiu Y, Odendahl N, Hudait A, Mason R, Bertram AK, Paesani F, DeMott PJ, Molinero V. Ice Nucleation Efficiency of Hydroxylated Organic Surfaces Is Controlled by Their Structural Fluctuations and Mismatch to Ice. J Am Chem Soc 2017; 139:3052-3064. [PMID: 28135412 DOI: 10.1021/jacs.6b12210] [Citation(s) in RCA: 112] [Impact Index Per Article: 16.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Heterogeneous nucleation of ice induced by organic materials is of fundamental importance for climate, biology, and industry. Among organic ice-nucleating surfaces, monolayers of long chain alcohols are particularly effective, while monolayers of fatty acids are significantly less so. As these monolayers expose to water hydroxyl groups with an order that resembles the one in the basal plane of ice, it was proposed that lattice matching between ice and the surface controls their ice-nucleating efficiency. Organic monolayers are soft materials and display significant fluctuations. It has been conjectured that these fluctuations assist in the nucleation of ice. Here we use molecular dynamic simulations and laboratory experiments to investigate the relationship between the structure and fluctuations of hydroxylated organic surfaces and the temperature at which they nucleate ice. We find that these surfaces order interfacial water to form domains with ice-like order that are the birthplace of ice. Both mismatch and fluctuations decrease the size of the preordered domains and monotonously decrease the ice freezing temperature. The simulations indicate that fluctuations depress the freezing efficiency of monolayers of alcohols or acids to half the value predicted from lattice mismatch alone. The model captures the experimental trend in freezing efficiencies as a function of chain length and predicts that alcohols have higher freezing efficiency than acids of the same chain length. These trends are mostly controlled by the modulation of the structural mismatch to ice. We use classical nucleation theory to show that the freezing efficiencies of the monolayers are directly related to their free energy of binding to ice. This study provides a general framework to relate the equilibrium thermodynamics of ice binding to a surface and the nonequilibrium ice freezing temperature and suggests that these could be predicted from the structure of interfacial water.
Collapse
Affiliation(s)
- Yuqing Qiu
- Department of Chemistry, The University of Utah , 315 South 1400 East, Salt Lake City, Utah 84112-0850, United States
| | - Nathan Odendahl
- Department of Chemistry, The University of Utah , 315 South 1400 East, Salt Lake City, Utah 84112-0850, United States
| | - Arpa Hudait
- Department of Chemistry, The University of Utah , 315 South 1400 East, Salt Lake City, Utah 84112-0850, United States
| | - Ryan Mason
- Department of Chemistry, University of British Columbia , Vancouver, British Columbia V6T 1Z1, Canada
| | - Allan K Bertram
- Department of Chemistry, University of British Columbia , Vancouver, British Columbia V6T 1Z1, Canada
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California, San Diego , La Jolla, California 92093, United States
| | - Paul J DeMott
- Department of Atmospheric Science, Colorado State University , Fort Collins, Colorado 80523-1371, United States
| | - Valeria Molinero
- Department of Chemistry, The University of Utah , 315 South 1400 East, Salt Lake City, Utah 84112-0850, United States
| |
Collapse
|
43
|
Dreischmeier K, Budke C, Wiehemeier L, Kottke T, Koop T. Boreal pollen contain ice-nucleating as well as ice-binding 'antifreeze' polysaccharides. Sci Rep 2017; 7:41890. [PMID: 28157236 PMCID: PMC5291224 DOI: 10.1038/srep41890] [Citation(s) in RCA: 57] [Impact Index Per Article: 8.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2016] [Accepted: 12/28/2016] [Indexed: 11/09/2022] Open
Abstract
Ice nucleation and growth is an important and widespread environmental process. Accordingly, nature has developed means to either promote or inhibit ice crystal formation, for example ice-nucleating proteins in bacteria or ice-binding antifreeze proteins in polar fish. Recently, it was found that birch pollen release ice-nucleating macromolecules when suspended in water. Here we show that birch pollen washing water exhibits also ice-binding properties such as ice shaping and ice recrystallization inhibition, similar to antifreeze proteins. We present spectroscopic evidence that both the ice-nucleating as well as the ice-binding molecules are polysaccharides bearing carboxylate groups. The spectra suggest that both polysaccharides consist of very similar chemical moieties, but centrifugal filtration indicates differences in molecular size: ice nucleation occurs only in the supernatant of a 100 kDa filter, while ice shaping is strongly enhanced in the filtrate. This finding may suggest that the larger ice-nucleating polysaccharides consist of clusters of the smaller ice-binding polysaccharides, or that the latter are fragments of the ice-nucleating polysaccharides. Finally, similar polysaccharides released from pine and alder pollen also display both ice-nucleating as well as ice-binding ability, suggesting a common mechanism of interaction with ice among several boreal pollen with implications for atmospheric processes and antifreeze protection.
Collapse
Affiliation(s)
- Katharina Dreischmeier
- Bielefeld University, Faculty of Chemistry, Atmospheric and Physical Chemistry, D-33615 Bielefeld, Germany
| | - Carsten Budke
- Bielefeld University, Faculty of Chemistry, Atmospheric and Physical Chemistry, D-33615 Bielefeld, Germany
| | - Lars Wiehemeier
- Bielefeld University, Faculty of Chemistry, Atmospheric and Physical Chemistry, D-33615 Bielefeld, Germany
- Bielefeld University, Faculty of Chemistry, Physical and Biophysical Chemistry, D-33615 Bielefeld, Germany
| | - Tilman Kottke
- Bielefeld University, Faculty of Chemistry, Physical and Biophysical Chemistry, D-33615 Bielefeld, Germany
| | - Thomas Koop
- Bielefeld University, Faculty of Chemistry, Atmospheric and Physical Chemistry, D-33615 Bielefeld, Germany
- Bielefeld University, Center for Molecular Materials, D-33615 Bielefeld, Germany
| |
Collapse
|
44
|
Bourque AJ, Locker CR, Rutledge GC. Heterogeneous Nucleation of an n-Alkane on Tetrahedrally Coordinated Crystals. J Phys Chem B 2017; 121:904-911. [PMID: 28071905 DOI: 10.1021/acs.jpcb.6b12590] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Heterogeneous nucleation refers to the propensity for phase transformations to initiate preferentially on foreign surfaces, such as vessel walls, dust particles, or formulation additives. In crystallization, the form of the initial nucleus has ramifications for the crystallographic form, morphology, and properties of the resulting solid. Nevertheless, the discovery and design of nucleating agents remains a matter of trial and error because of the very small spatiotemporal scales over which the critical nucleus is formed and the extreme difficulty of examining such events empirically. Using molecular dynamics simulations, we demonstrate a method for the rapid screening of entire families of materials for activity as nucleating agents and for characterizing their mechanism of action. The method is applied to the crystallization of n-pentacontane, a model surrogate for polyethylene, on the family of tetrahedrally coordinated crystals, including diamond and silicon. A systematic variation of parameters in the interaction potential permits a comprehensive, physically based screening of nucleating agents in this class of materials, including both real and hypothetical candidates. The induction time for heterogeneous nucleation is shown to depend strongly on crystallographic registry between the nucleating agent and the critical nucleus, indicative of an epitaxial mechanism in this class of materials. Importantly, the severity of this registry requirement weakens with decreasing rigidity of the substrate and increasing strength of attraction to the surface of the nucleating agent. Employing this method, a high-throughput computational screening of nucleating agents becomes possible, facilitating the discovery of novel nucleating agents within a broad "materials genome" of possible additives.
Collapse
Affiliation(s)
- Alexander J Bourque
- Department of Chemical Engineering, Massachusetts Institute of Technology , Cambridge, Massachusetts 02139, United States
| | - C Rebecca Locker
- ExxonMobil Research and Engineering Company , Annandale, New Jersey 08801, United States
| | - Gregory C Rutledge
- Department of Chemical Engineering, Massachusetts Institute of Technology , Cambridge, Massachusetts 02139, United States
| |
Collapse
|
45
|
|
46
|
He X, Shen Y, Hung FR, Santiso EE. Heterogeneous nucleation from a supercooled ionic liquid on a carbon surface. J Chem Phys 2016; 145:211919. [DOI: 10.1063/1.4963336] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
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
- Center for Computation & Technology, 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
|
47
|
Glatz B, Sarupria S. The surface charge distribution affects the ice nucleating efficiency of silver iodide. J Chem Phys 2016; 145:211924. [DOI: 10.1063/1.4966018] [Citation(s) in RCA: 49] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Affiliation(s)
- Brittany Glatz
- Department of Chemical and Biomolecular Engineering, Clemson University, Clemson, South Carolina 29634, USA
| | - Sapna Sarupria
- Department of Chemical and Biomolecular Engineering, Clemson University, Clemson, South Carolina 29634, USA
| |
Collapse
|
48
|
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
|
49
|
Yang H, Ma C, Li K, Liu K, Loznik M, Teeuwen R, van Hest JCM, Zhou X, Herrmann A, Wang J. Tuning Ice Nucleation with Supercharged Polypeptides. ADVANCED MATERIALS (DEERFIELD BEACH, FLA.) 2016; 28:5008-5012. [PMID: 27119590 DOI: 10.1002/adma.201600496] [Citation(s) in RCA: 41] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/27/2016] [Indexed: 06/05/2023]
Abstract
Supercharged unfolded polypeptides (SUPs) are exploited for controlling ice nucleation via tuning the nature of charge and charge density of SUPs. The results show that positively charged SUPs facilitate ice nucleation, while negatively charged ones suppress it. Moreover, the charge density of the SUP backbone is another parameter to control it.
Collapse
Affiliation(s)
- Huige Yang
- Institute of Chemistry, Chinese Academy of Sciences, Beijing, 100190, P. R. China
- School of Physical Sciences, University of Chinese Academy of Sciences, Beijing, 100049, P. R. China
| | - Chao Ma
- Department of Polymer Chemistry and Bioengineering, Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747AG, Groningen, The Netherlands
| | - Kaiyong Li
- School of Materials Science and Engineering, Luoyang Institute of Science and Technology, Luoyang, 471023, P. R. China
| | - Kai Liu
- Department of Polymer Chemistry and Bioengineering, Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747AG, Groningen, The Netherlands
| | - Mark Loznik
- Department of Polymer Chemistry and Bioengineering, Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747AG, Groningen, The Netherlands
| | - Rosalie Teeuwen
- Institute for Molecules and Materials, Radboud University Nijmegen, 6525, AJ, Nijmegen, The Netherlands
| | - Jan C M van Hest
- Institute for Molecules and Materials, Radboud University Nijmegen, 6525, AJ, Nijmegen, The Netherlands
| | - Xin Zhou
- School of Physical Sciences, University of Chinese Academy of Sciences, Beijing, 100049, P. R. China
| | - Andreas Herrmann
- Department of Polymer Chemistry and Bioengineering, Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747AG, Groningen, The Netherlands
| | - Jianjun Wang
- Institute of Chemistry, Chinese Academy of Sciences, Beijing, 100190, P. R. China
- School of Physical Sciences, University of Chinese Academy of Sciences, Beijing, 100049, P. R. China
| |
Collapse
|
50
|
Sosso G, Chen J, Cox SJ, Fitzner M, Pedevilla P, Zen A, Michaelides A. Crystal Nucleation in Liquids: Open Questions and Future Challenges in Molecular Dynamics Simulations. Chem Rev 2016; 116:7078-116. [PMID: 27228560 PMCID: PMC4919765 DOI: 10.1021/acs.chemrev.5b00744] [Citation(s) in RCA: 390] [Impact Index Per Article: 48.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2015] [Indexed: 11/28/2022]
Abstract
The nucleation of crystals in liquids is one of nature's most ubiquitous phenomena, playing an important role in areas such as climate change and the production of drugs. As the early stages of nucleation involve exceedingly small time and length scales, atomistic computer simulations can provide unique insights into the microscopic aspects of crystallization. In this review, we take stock of the numerous molecular dynamics simulations that, in the past few decades, have unraveled crucial aspects of crystal nucleation in liquids. We put into context the theoretical framework of classical nucleation theory and the state-of-the-art computational methods by reviewing simulations of such processes as ice nucleation and the crystallization of molecules in solutions. We shall see that molecular dynamics simulations have provided key insights into diverse nucleation scenarios, ranging from colloidal particles to natural gas hydrates, and that, as a result, the general applicability of classical nucleation theory has been repeatedly called into question. We have attempted to identify the most pressing open questions in the field. We believe that, by improving (i) existing interatomic potentials and (ii) currently available enhanced sampling methods, the community can move toward accurate investigations of realistic systems of practical interest, thus bringing simulations a step closer to experiments.
Collapse
Affiliation(s)
- Gabriele
C. Sosso
- Thomas Young Centre, London
Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street WC1E
6BT London, U.K.
| | - Ji Chen
- Thomas Young Centre, London
Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street WC1E
6BT London, U.K.
| | | | - Martin Fitzner
- Thomas Young Centre, London
Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street WC1E
6BT London, U.K.
| | - Philipp Pedevilla
- Thomas Young Centre, London
Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street WC1E
6BT London, U.K.
| | - Andrea Zen
- Thomas Young Centre, London
Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street WC1E
6BT London, U.K.
| | - Angelos Michaelides
- Thomas Young Centre, London
Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street WC1E
6BT London, U.K.
| |
Collapse
|