An official website of the United States government
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
On the challenge of sampling multiple nucleation pathways: A case study of heterogeneous ice nucleation on FCC (211) surface
Heterogeneous ice nucleation is ubiquitous but its microscopic mechanisms can be extraordinarily complex even on a simple surface. Such complexity poses a challenge in modeling nucleation using advanced sampling methods. Here, we investigate heterogeneous ice nucleation on an FCC (211) surface by a forward flux sampling (FFS) method to understand how the complexity in nucleation pathways affects the accuracy of FFS. We first show the commonly adopted, size-based order parameter fails to describe heterogeneous ice nucleation on the FCC (211) surface. Inclusion of geometric anisotropy of ice nucleus as an additional descriptor is found to significantly improve the quality of the size-based order parameter for the current system. Subsequent application of this new order parameter in FFS identifies two competing ice nucleation pathways in the system: a primary-prism-planed (PPP) path and a secondary-prism-planed (SPP) path, both leading to the formation of hexagonal ice but with different crystalline orientations. Although the PPP pathway dominates ice nucleation on the FCC (211) surface, the occurrence of the less efficient SPP pathway, particularly its strong presence at the early stage of FFS, is found to yield a significant statistical uncertainty in the calculated FFS rate constant. We develop a two-path model that enables gaining a general, quantitative understanding of the impact of initial finite sampling on the reliability of FFS calculations in the presence of multiple nucleation pathways. Our study also suggests a few general strategies for improving the accuracy of FFS when exploring unknown but complex systems.
more »
« less
- Award ID(s):
- 2053330
- PAR ID:
- 10430339
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 158
- Issue:
- 12
- ISSN:
- 0021-9606
- Page Range / eLocation ID:
- 124501
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)Heterogeneous ice nucleation is a common process in the atmosphere, but relatively little is known about the role of different surface characteristics on the promotion of ice nucleation. We have used a series of iron oxides as a model system to study the role of lattice mismatch and defects induced by milling on ice nucleation activity. The iron oxides include wüstite (FeO), hematite (Fe 2 O 3 ), magnetite (Fe 3 O 4 ), and goethite (FeOOH). The iron oxides were characterized by X-ray diffraction (XRD) and Brunauer–Emmett–Teller (BET) surface area measurements. The immersion freezing experiments were performed using an environmental chamber. Wüstite (FeO) had the highest ice nucleation activity, which we attribute to its low lattice mismatch with hexagonal ice and the exposure of Fe–OH after milling. A comparison study of MnO and wüstite (FeO) with milled and sieved samples for each suggests that physical defects alone result in only a slight increase in ice nucleation activity. Despite differences in the molecular formula and surface groups, hematite (Fe 2 O 3 ), magnetite (Fe 3 O 4 ), and goethite (FeOOH) had similar ice nucleation activities, which may be attributed to their high lattice mismatch to hexagonal ice. This study provides further insight into the characteristics of a good heterogeneous ice nucleus and, more generally, helps to elucidate the interactions between aerosol particles and ice particles in clouds.more » « less
-
Heterogeneous ice nucleation in the atmosphere impacts climate, but the magnitude of the effect of ice clouds on radiative forcing is uncertain. Surfaces that promote ice nucleation are varied. Because O, Si, and Al are the most abundant elements in the Earth's crust, understanding how the Si : Al ratio impacts the ice nucleation activity of aluminosilicates through exploration of synthetic ZSM-5 samples provides a good model system. This paper investigates the immersion freezing of ZSM-5 samples with varying Si : Al ratios. Ice nucleation temperature increases with increasing surface Al content. Additionally, when ammonium, a common cation in aerosol particles, is adsorbed to the zeolite surface, initial freezing temperatures are reduced by up to 6 °C in comparison to proton-terminated zeolite surfaces. This large decrease in ice nucleation activity in the presence of ammonium suggests that the cation can interact with the surface to block or modify active sites. Our results on synthetic samples in which the surface composition is tunable gives insight into the role of surfaces in heterogeneous ice nucleation processes in the atmosphere. We emphasize the importance of examining surface chemical heterogeneities in ice nucleating particles that could result from a variety of aging pathways for a deeper understanding of the freezing mechanism.more » « less
-
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
-
Path sampling approaches have become invaluable tools to explore the mechanisms and dynamics of the so-called rare events that are characterized by transitions between metastable states separated by sizable free energy barriers. Their practical application, in particular to ever more complex molecular systems, is, however, not entirely trivial. Focusing on replica exchange transition interface sampling (RETIS) and forward flux sampling (FFS), we discuss a range of analysis tools that can be used to assess the quality and convergence of such simulations, which is crucial to obtain reliable results. The basic ideas of a step-wise evaluation are exemplified for the study of nucleation in several systems with different complexities, providing a general guide for the critical assessment of RETIS and FFS simulations.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1063/5.0144712
