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
Pseudo-grand canonical molecular dynamics via volumetrically controlled osmotic pressure
Molecular dynamics simulations are typically constrained to have a fixed number of particles, which limits our capability to simulate chemical and physical processes where the composition of the system changes during the simulation time. Typical examples are the calculation of nucleation and crystal growth rates in heterogeneous solutions where the driving force depends on the composition of the fluid. Constant chemical potential molecular dynamics simulations would instead be required to compute time-independent growth and nucleation rates. While this can, in principle, be achieved through the addition and deletion of particles using the grand canonical partition function, this is very inefficient in the condensed phase due to the low acceptance probability of these events. Adaptive resolution schemes, which use a reservoir of non-interacting particles that can be transformed into solute particles, circumvent this problem, but at the cost of relatively complicated code implementations. In this work, a simpler approach is proposed that uses harmonic volumetric restraints to control the solute osmotic pressure, which can be considered a proxy for the system’s chemical potential. The osmotic pressure regulator is demonstrated to reproduce the expected properties of ideal gases and ideal solutions. Using the mW water model, the osmotic pressure regulator is shown to provide a constant growth rate for ice in the presence of an electrolyte solution, unlike what standard molecular dynamics simulations would produce.
more »
« less
- Award ID(s):
- 2053235
- PAR ID:
- 10636212
- Publisher / Repository:
- American Institute of Physics
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 163
- Issue:
- 11
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We use linear stability analysis and direct numerical simulations to investigate the coupling between centrifugal instabilities, solute transport and osmotic pressure in a Taylor–Couette configuration that models rotating dynamic filtration devices. The geometry consists of a Taylor–Couette cell with a superimposed radial throughflow of solvent across two semi-permeable cylinders. Both cylinders totally reject the solute, inducing the build-up of a concentration boundary layer. The solute retroacts on the velocity field via the osmotic pressure associated with the concentration differences across the semi-permeable cylinders. Our results show that the presence of osmotic pressure strongly alters the dynamics of the centrifugal instabilities and substantially reduces the critical conditions above which Taylor vortices are observed. It is also found that this enhancement of the hydrodynamic instabilities eventually plateaus as the osmotic pressure is further increased. We propose a mechanism to explain how osmosis and instabilities cooperate and develop an analytical criterion to bound the parameter range for which osmosis fosters the hydrodynamic instabilities.more » « less
-
null (Ed.)The microscopic properties that determine hygroscopic behavior are complex. The importance of hygroscopicity to many areas, and particularly atmospheric chemistry, in terms of aerosol growth and cloud nucleation, mandate the need for robust models to understand this behavior. Toward this end, we have employed molecular dynamics simulations to calculate hygroscopicity from atomistic models using free energy perturbation. We find that currently available force fields may not be well-suited to modeling the extreme environments of aerosol particles. Nonetheless, the results illuminate some shortcomings in our current understanding of hygroscopic growth and cloud nucleation. The most widely used model of hygroscopicity, κ-Köhler Theory (κKT), breaks down in the case of deviations from ideal solution behavior and empirical adjustments within the simplified framework cannot account for non-ideal behavior. A revised model that incorporates non-ideal mixing rescues the general framework of κKT and allows us to understand our simulation results as well as the behavior of atmospheric aerosols over the full range of humidity. The revised model shows that non-ideal mixing dominates hygroscopic growth at subsaturation humidity. Thus, a model based on ideal mixing will fail to predict subsaturation growth from cloud condensation nucleus (CCN) activation or vice versa ; a single parameter model for hygroscopicity will generally be insufficient to extrapolate across wide ranges of humidity. We argue that in many cases, when data are limited to subsaturation humidity, an empirical model for non-ideal mixing may be more successful than one for ideal mixing.more » « less
-
Brownian Motion, with some persistence in the direction of motion, typically known as active Brownian Motion, has been observed in many significant chemical and biological transport processes. Here, we present a model of drifted Brownian Motion that considers a nonlinear stochastic drift with constant or fluctuating diffusivity. The interplay between nonlinearity and structural heterogeneity of the environment can explain three essential features of active transport. These features, which are commonly observed in experiments and molecular dynamics simulations, include transient superdiffusion, ephemeral non-Gaussian displacement distribution, and non-monotonic evolution of non-Gaussian parameter. Our results compare qualitatively well with experiments of self-propelled particles in simple hydrogen peroxide solutions and molecular dynamics simulations of self-propelled particles in more complex settings such as viscoelastic polymeric media.more » « less
-
Abstract The structure and composition of the crystal growth unit are of huge fundamental and practical consequence. We propose a method to identify the solute species that incorporates into the growth site on crystal surfaces, the kinks, which rests on the kinetics of the elementary reaction at the kinks. We use as model crystals olanzapine, an antipsychotic medication, and etioporphyrin I, a field‐effect transistor. We combine time‐resolvedin situatomic force microscopy with Raman and absorption spectroscopies, complemented by density functional theory and all‐atom molecular dynamics modeling of the solutions. We show that the structure of the growth unit cannot be deduced neither from the solute oligomers nor from the crystal structure. Chemical kinetics analyses reveal that if the dominant solute species is the one that incorporates into the crystal growth sites, then the kinetics of layer growth complies with a monomolecular rate law. By contrast, if the crystal growth unit assembles from two units of the dominant solute form, a bimolecular rate law ensues. Solutions of both olanzapine and etioporphyrin I are dominated by solute monomers, which exist in equilibrium with a minority of dimers. Whereas numerous olanzapine crystal structures incorporate dimer motifs, etioporphyrin I crystals organize as stacks of monomers. Olanzapine crystal grow by incorporation of dimers. One of the studied face of etioporphyrin I grows by incorporation of the majority monomers, whereas the other one selects the minority dimers as a growth unit. The results highlight the power of the crystallization kinetics analyses to identify the growth unit and illuminate one of the most challenging issues of crystal growth.more » « less
