skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: RSeeds: Rigid Seeding Method for Studying Heterogeneous Crystal Nucleation
Heterogeneous nucleation is the dominant form of liquid-to-solid transition in na- ture. Although molecular simulations are most uniquely suited to study nucleation, the waiting time to observe even a single nucleation event can easily exceed 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 successful at dramatically expanding the range of computationally accessible nucleation rates in simulation studies of ho- mogeneous crystal nucleation. In this paper, we introduce a new seeding method for heterogeneous nucleation called Rigid Seeding (RSeeds). Crystalline seeds are treated as pseudo-rigid 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 nucle- ation 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.  more » « less
Award ID(s):
2224643
PAR ID:
10516624
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
ACS
Date Published:
Journal Name:
The Journal of Physical Chemistry B
Volume:
127
Issue:
18
ISSN:
1520-6106
Page Range / eLocation ID:
4112 to 4125
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract. Heterogeneous ice nucleation is thought to be the primary pathway for the formation of ice in mixed-phase clouds, with the number of active ice-nucleating particles (INPs) increasing rapidly with decreasing temperature. Here, molecular-dynamics simulations of heterogeneous ice nucleation demonstrate that the ice nucleation rate is also sensitive to pressure and that negative pressure within supercooled water shifts freezing temperatures to higher temperatures. Negative pressure, or tension, occurs naturally in water capillary bridges and pores and can also result from water agitation. Capillary bridge simulations presented in this study confirm that negative Laplace pressure within the water increases heterogeneous-freezing temperatures. The increase in freezing temperatures with negative pressure is approximately linear within the atmospherically relevant range of 1 to −1000 atm. An equation describing the slope depends on the latent heat of freezing and the molar volume difference between liquid water and ice. Results indicate that negative pressures of −500 atm, which correspond to nanometer-scale water surface curvatures, lead to a roughly 4 K increase in heterogeneous-freezing temperatures. In mixed-phase clouds, this would result in an increase of approximately 1 order of magnitude in active INP concentrations. The findings presented here indicate that any process leading to negative pressure in supercooled water may play a role in ice formation, consistent with experimental evidence of enhanced ice nucleation due to surface geometry or mechanical agitation of water droplets. This points towards the potential for dynamic processes such as contact nucleation and droplet collision or breakup to increase ice nucleation rates through pressure perturbations. 
    more » « less
  2. Ice nucleation in the atmosphere influences cloud properties, altering precipitation and the radiative balance, ultimately regulating Earth’s climate. An accepted ice nucleation pathway, known as deposition nucleation, assumes a direct transition of water from the vapor to the ice phase, without an intermediate liquid phase. However, studies have shown that nucleation occurs through a liquid phase in porous particles with narrow cracks or surface imperfections where the condensation of liquid below water saturation can occur, questioning the validity of deposition nucleation. We show that deposition nucleation cannot explain the strongly enhanced ice nucleation efficiency of porous compared with nonporous particles at temperatures below −40 °C and the absence of ice nucleation below water saturation at −35 °C. Using classical nucleation theory (CNT) and molecular dynamics simulations (MDS), we show that a network of closely spaced pores is necessary to overcome the barrier for macroscopic ice-crystal growth from narrow cylindrical pores. In the absence of pores, CNT predicts that the nucleation barrier is insurmountable, consistent with the absence of ice formation in MDS. Our results confirm that pore condensation and freezing (PCF), i.e., a mechanism of ice formation that proceeds via liquid water condensation in pores, is a dominant pathway for atmospheric ice nucleation below water saturation. We conclude that the ice nucleation activity of particles in the cirrus regime is determined by the porosity and wettability of pores. PCF represents a mechanism by which porous particles like dust could impact cloud radiative forcing and, thus, the climate via ice cloud formation. 
    more » « less
  3. Abstract. Ice growth from vapor deposition is an important process for the evolution of cirrus clouds, but the physics of depositional ice growth at the low temperatures (<235 K) characteristic of the upper troposphere and lower stratosphere is not well understood. Surface attachment kinetics, generally parameterized as a deposition coefficient αD, control ice crystal habit and also may limit growth rates in certain cases, but significant discrepancies between experimental measurements have not been satisfactorily explained. Experiments on single ice crystals have previously indicated the deposition coefficient is a function of temperature and supersaturation, consistent with growth mechanisms controlled by the crystal's surface characteristics. Here we use observations from cloud chamber experiments in the Aerosol Interactions and Dynamics in theAtmosphere (AIDA) aerosol and cloud chamber to evaluate surface kinetic models in realistic cirrus conditions. These experiments have rapidly changing temperature, pressure, and ice supersaturation such that depositional ice growth may evolve from diffusion limited to surface kinetics limited over the course of a single experiment. In Part 1, we describe the adaptation of a Lagrangian parcel model with the Diffusion Surface Kinetics Ice Crystal Evolution (DiSKICE) model (Zhang and Harrington, 2014) to the AIDA chamber experiments. We compare the observed ice water content and saturation ratios to that derived under varying assumptions for ice surface growth mechanisms for experiments simulating ice clouds between 180 and 235 K and pressures between 150 and 300 hPa. We found that both heterogeneous and homogeneous nucleation experiments at higher temperatures (>205 K) could generally be modeled consistently with either a constant deposition coefficient or the DiSKICE model assuming growth on isometric crystals via abundant surface dislocations. Lower-temperature experiments showed more significant deviations from any depositional growth model, with different ice growth rates for heterogeneous and homogeneous nucleation experiments. 
    more » « less
  4. The Classical Nucleation Theory (CNT) has played a key role in crystal nucleation studies since the 19th century and has significantly advanced the understanding of nucleation. However, certain key assumptions of CNT, such as a compact and spherical nucleating cluster and the concept of individual diffusive jumps are questionable. The results of molecular dynamics (MD) studies of crystal nucleation in a Al 20 Ni 60 Zr 20 metallic liquid demonstrate that the nucleating cluster is neither spherical nor compact. The seeding method was employed to determine the critical cluster size and nucleation parameters from CNT, which were then compared to those derived from the Mean First Passage Time (MFPT) method. While the CNT-based nucleation rate aligns well with experimental data from similar metallic liquids, the MFPT rate differs significantly. Further, contrary to the assumption of individual jumps for atoms to join the nucleating cluster, a cooperative mechanism of attachment or detachment is observed. This is accompanied by synchronized changes in the local potential energy. Similar cooperative motion also appeared in a non-classical nucleation process, particularly during the coalescence of nuclei. 
    more » « less
  5. Heterogeneous ice nucleation in the atmosphere regulates cloud properties, such as phase (ice versus liquid) and lifetime. Aerosol particles of marine origin are relevant ice nucleating particle sources when marine aerosol layers are lifted over mountainous terrain and in higher latitude ocean boundary layers, distant from terrestrial aerosol sources. Among many particle compositions associated with ice nucleation by sea spray aerosols are highly saturated fatty acids. Previous studies have not demonstrated their ability to freeze dilute water droplets. This study investigates ice nucleation by monolayers at the surface of supercooled droplets and as crystalline particles at temperatures exceeding the threshold for homogeneous freezing. Results show the poor efficiency of long chain fatty acid (C16, C18) monolayers in templating freezing of pure water droplets and seawater subphase to temperatures of at least −30 °C, consistent with theory. This contrasts with freezing of fatty alcohols (C22 used here) at nearly 20 °C warmer. Evaporation of μL-sized droplets to promote structural compression of a C19 acid monolayer did not favor warmer ice formation of drops. Heterogeneous ice nucleation occurred for nL-sized droplets condensed on 5 to 100 μm crystalline particles of fatty acid (C12 to C20) at a range of temperatures below −28 °C. These experiments suggest that fatty acids nucleate ice at warmer than −36 °C only when the crystalline phase is present. Rough estimates of ice active site densities are consistent with those of marine aerosols, but require knowledge of the proportion of surface area comprised of fatty acids for application. 
    more » « less