Molecular dynamics simulations of ice nucleation by electric fields

Molecular dynamics simulations are used to investigate heterogeneous ice nucleation in model systems where an electric field acts on water molecules within 10-20 ? of a surface. Two different water models (the six-site and TIP4P/Ice models) are considered, and in both cases, it is shown that a surface field can serve as a very effective ice nucleation catalyst in supercooled water. Ice with a ferroelectric cubic structure nucleates near the surface, and dipole disordered cubic ice grows outward from the surface layer. We examine the influences of temperature and two important field parameters, the field strength and distance from the surface over which it acts, on the ice nucleation process. For the six-site model, the highest temperature where we observe field-induced ice nucleation is 280 K, and for TIP4P/Ice 270 K (note that the estimated normal freezing points of the six-site and TIP4P/Ice models are ～289 and ～270 K, respectively). The minimum electric field strength required to nucleate ice depends a little on how far the field extends from the surface. If it extends 20 ?, then a field strength of 1.5 × 10(9) V/m is effective for both models. If the field extent is 10 ?, then stronger fields are required (2.5 × 10(9) V/m for TIP4P/Ice and 3.5 × 10(9) V/m for the six-site model). Our results demonstrate that fields of realistic strength, that act only over a narrow surface region, can effectively nucleate ice at temperatures not far below the freezing point. This further supports the possibility that local electric fields can be a significant factor influencing heterogeneous ice nucleation in physical situations. We would expect this to be especially relevant for ice nuclei with very rough surfaces where one would expect local fields of varying strength and direction.     

Kinetic Monte Carlo simulations of nucleation and growth in electrodeposition

Nucleation and growth during bulk electrodeposition is studied using kinetic Monte Carlo (KMC) simulations. Ion transport in solution is modeled using Brownian dynamics, and the kinetics of nucleation and growth are dependent on the probabilities of metal-on-substrate and metal-on-metal deposition. Using this approach, we make no assumptions about the nucleation rate, island density, or island distribution. The influence of the attachment probabilities and concentration on the time-dependent island density and current transients is reported. Various models have been assessed by recovering the nucleation rate and island density from the current-time transients.     

Homogeneous ice nucleation from supercooled water

Homogeneous ice nucleation from supercooled water was studied in the temperature range of 220-240 K through combining the forward flux sampling method (Allen et al., J. Chem. Phys., 2006, 124, 024102) with molecular dynamics simulations (FFS/MD), based on a recently developed coarse-grained water model (mW) (Molinero et al., J. Phys. Chem. B, 2009, 113, 4008). The calculated ice nucleation rates display a strong temperature dependence, ranging from 2.148 ± 0.635 × 10(25) m(-3) s(-1) at 220 K to 1.672 ± 0.970 × 10(-7) m(-3) s(-1) at 240 K. These rates can be fitted according to the classical nucleation theory, yielding an estimate of the effective ice-water interface energy γ(ls) of 31.01 ± 0.21 mJ m(-2) for the mW water model. Compared to experiments, our calculation underestimates the homogeneous ice nucleation rate by a few orders of magnitude. Possible reasons for the discrepancy are discussed. The nucleating ice embryo contains both cubic ice Ic and hexagonal ice Ih, with the fraction of each structure being roughly 50% when the critical size is reached. In particular, a novel defect structure containing nearly five-fold twin boundaries is identified in the ice clusters formed during nucleation. The way such defect structure is formed is found to be different from mechanisms proposed for the formation of the same defect in metallic nanoparticles and thin film. The quasi five-fold twin boundary structure found here is expected to occur in the crystallization of a wide range of materials with the diamond cubic structure, including ice.     

Finite-size effects in simulations of nucleation

We investigate the importance of finite-size effects in simulations of nucleation processes. Most molecular dynamics simulations of first order phase transitions, such as vapor-liquid nucleation, are performed in the canonical NVT ensemble where, owing to the fixed total number of molecules N, the growth of the new phase causes the depletion of the metastable phase. This effect may lead to significant errors in the simulation and even to the impossibility of observing nucleation in a small finite system. We present a theory to estimate the system size beyond which these finite-size effects are expected to be negligible. This optimization saves valuable calculation time and can extend the range of supersaturations and rates attainable by simulations by several orders of magnitude. Our results are applicable to diverse situations, such as crystallization, capillary condensation, or the melting of nanoclusters.     

Gas hydrate fast nucleation from melting ice and quiescent growth along vertical heat transfer tube

Gas hydrates, or clathrate hydrates, are ice-likecrystal, composed of host lattice (cavities) formed byhydrogen-bonded water molecules, and other guestmolecules called guest molecules. The guest mole-cules act with host lattice in weak van der Waals force…     

Kinetic aspects of the thermostatted growth of ice from supercooled water in simulations

In experiments, the growth rate of ice from supercooled water is seen to increase with the degree of supercooling, that is, the lower the temperature, the faster the crystallization takes place. In molecular dynamics simulations of the freezing process, however, the temperature is usually kept constant by means of a thermostat that artificially removes the heat released during the crystallization by scaling the velocities of the particles. This direct removal of energy from the system replaces a more realistic heat-conduction mechanism and is believed to be responsible for the curious observation that the thermostatted ice growth proceeds fastest near the melting point and more slowly at lower temperatures, which is exactly opposite to the experimental findings [M. A. Carignano, P. B. Shepson, and I. Szleifer, Mol. Phys. 103, 2957 (2005)]. This trend is explained by the diffusion and the reorientation of molecules in the liquid becoming the rate-determining steps for the crystal growth, both of which are slower at low temperatures. Yet, for a different set of simulations, a kinetic behavior analogous to the experimental finding has been reported [H. Nada and Y. Furukawa, J. Crystal Growth 283, 242 (2005)]. To clarify this apparent contradiction, we perform relatively long simulations of the TIP4P/Ice model in an extended range of temperatures. The temperature dependence of the thermostatted ice growth is seen to be more complex than was previously reported: The crystallization process is very slow close to the melting point at 270 K, where the thermodynamic driving force for the phase transition is weak. On lowering the temperature, the growth rate initially increases, but displays a maximum near 260 K. At even lower temperatures, the freezing process slows down again due to the reduced diffusivity in the liquid. The velocity of the thermostatted melting process, in contrast, shows a monotonic increase upon raising the temperature beyond the normal melting point. In this case, the effects of the increasing thermodynamic driving force and the faster diffusion at higher temperatures reinforce each other. In the context of this study, we also report data for the diffusion coefficient as a function of temperature for the water models TIP4P/Ice and TIP4P/2005.     

Kinetics of ice nucleation initiated by molecular crystals

Heterogeneous nucleation of ice in the presence of a number of organic molecular crystals is characterized by effective kinetic constants. We have found a correlation between the values of the kinetic constants and the amount of absorbed water within the volume of the crystal.Translated from Teoreticheskaya i Éksperimental'naya Khimiya, Vol. 30, No. 6, pp. 323–327, November–December, 1994.     

Dynamics of ice nucleation on water repellent surfaces

Prevention of ice accretion and adhesion on surfaces is relevant to many applications, leading to improved operation safety, increased energy efficiency, and cost reduction. Development of passive nonicing coatings is highly desirable, since current antiicing strategies are energy and cost intensive. Superhydrophobicity has been proposed as a lead passive nonicing strategy, yet the exact mechanism of delayed icing on these surfaces is not clearly understood. In this work, we present an in-depth analysis of ice formation dynamics upon water droplet impact on surfaces with different wettabilities. We experimentally demonstrate that ice nucleation under low-humidity conditions can be delayed through control of surface chemistry and texture. Combining infrared (IR) thermometry and high-speed photography, we observe that the reduction of water-surface contact area on superhydrophobic surfaces plays a dual role in delaying nucleation: first by reducing heat transfer and second by reducing the probability of heterogeneous nucleation at the water-substrate interface. This work also includes an analysis (based on classical nucleation theory) to estimate various homogeneous and heterogeneous nucleation rates in icing situations. The key finding is that ice nucleation delay on superhydrophobic surfaces is more prominent at moderate degrees of supercooling, while closer to the homogeneous nucleation temperature, bulk and air-water interface nucleation effects become equally important. The study presented here offers a comprehensive perspective on the efficacy of textured surfaces for nonicing applications.     

Towards understanding ice nucleation by long chain alcohols

We demonstrate that infrared spectra of water covered by a film of heptadecanol show a continuous spectral shift, from a band characteristic of liquid water to one characteristic of ice, as the temperature is ramped from -10 to -17 degrees C. Experiments with pure water and water covered by films of long chain alkanes showed no such spectral shift. Analysis of the CH2 stretching features in the alcohols' absorbance bands reveals simultaneous structural changes within the alcohol film. We hypothesize that the spectral shift in the water band is due to an increasing fraction of water molecules participating in icelike clusters and that these clusters are stabilized by the influence of the flexible alcohol film.     

Investigating the effects of solid surfaces on ice nucleation

Understanding the role played by solid surfaces in ice nucleation is a significant step toward designing anti-icing surfaces. However, the uncontrollable impurities in water and surface heterogeneities remain a great challenge for elucidating the effects of surfaces on ice nucleation. Via a designed process of evaporation, condensation, and subsequent ice formation in a closed cell, we investigate the ice nucleation of ensembles of condensed water microdroplets on flat, solid surfaces with completely different wettabilities. The water microdroplets formed on flat, solid surfaces by an evaporation and condensation process exclude the uncontrollable impurities in water, and the effects of surface heterogeneities can be minimized through studying the freezing of ensembles of separate and independent water microdroplets. It is found that the normalized surface ice nucleation rate on a hydrophilic surface is about 1 order of magnitude lower than that on a hydrophobic surface. This is ascribed to the difference in the viscosity of interfacial water and the surface roughness.     

Influence of thermostats and carrier gas on simulations of nucleation

We investigate the influence of carrier gas and thermostat on molecular dynamics (MD) simulations of nucleation. The task of keeping the temperature constant in MD simulations is not trivial and an inefficient thermalization may have a strong influence on the results. Different thermostating mechanisms have been proposed and used in the past. In particular, we analyze the efficiency of velocity rescaling, Nose-Hoover, and a carrier gas (mimicking the experimental situation) by extensive MD simulations. Since nucleation is highly sensitive to temperature, one would expect that small variations in temperature might lead to differences in nucleation rates of up to several orders of magnitude. Surprisingly, the results indicate that the choice of the thermostating method in a simulation does not have--at least in the case of Lennard-Jones argon--a very significant influence on the nucleation rate. These findings are interpreted in the context of the classical theory of Feder et al. [Adv. Phys. 15, 111 (1966)] by analyzing the temperature distribution of the nucleating clusters. We find that the distribution of cluster temperatures is non-Gaussian and that subcritically sized clusters are colder while postcritically sized clusters are warmer than the bath temperature. However, the average temperature of all clusters is found to be always higher than the bath temperature.     

Particle nucleation and growth models

This review concentrates on the progress of modeling the nucleation process of particles by the balanced nucleation-growth (BNG) process. The BNG model will be compared with other models that try to predict material nucleation. Compared to other models, the BNG model allows quantifying the nucleation rate, maximum growth rate, and supersaturation during the nucleation period as a function of nucleation efficiency and maximum growth rate of the crystals. From this model, equations are derived that correlate the number of stable crystals formed with molar addition rate of reactants, solubility of the crystals, and temperature. The BNG model predicts the experimental result that many crystallization processes result in a limited number of crystals followed by growth. The model also predicts that factors like diffusion and kinetically controlled growth process, Ostwald ripening agents and growth restrainers control the crystal number. Equations are given for each of the variables that agree with experiments. The BNG model predicts the conditions for renucleation (formation of new crystals during precipitation). It leads to new equations for the prediction of crystal number and crystal size during controlled continuous precipitation in the continuous stirred tank reactor (CSTR) as a function of precipitation conditions.     

Computer simulations of nucleation of nanoparticle superclusters from solution

This paper presents simulation studies of nanoparticle supercluster (NPSC) nucleation from a temperature quenched system. The nanoparticles are represented as 5 nm, spherical gold nanoparticles ligated with alkane thiols. The pair potential accounts for the van der Waals interaction between the metallic cores and ligand-ligand and ligand-solvent interactions. Phenomena well-known for molecular systems are observed including a prenucleation induction period, fluctuating prenucleation clusters that predominately add monomers one at a time, a critical nucleus size, and growth of NPSCs from solution in the presence of an equilibrium supernatant, all consistent with classical nucleation theory. However, only the largest prenucleating clusters are dense, and the cluster size can occasionally range greater than the critical size in the prenucleation regime until a cluster with low enough energy occurs, then nucleation ensues. Late in the nucleation process, the clusters display a crystalline structure that is a random mix of face-centered cubic (fcc) and hexagonal close-packed (hcp) lattices and indistinguishable from a randomized icosahedra structure.     

Secondary nucleation and growth of ZnO

Recently we discovered that under certain conditions new crystal growth (branch) can be induced on specific crystalline planes of the same material. This is a new phenomenon and is in sharp contrast to typical nucleation and growth in which a crystal will simply grow larger in preferred directions depending on the surface energy of the specific crystalline planes. Based on our observation, we developed a sequential nucleation and growth technique offering the power to assemble complex hierarchical crystals step-by-step. However, the key questions of when and how the secondary nucleation takes place have not been answered. Here we systematically study secondary ZnO crystal growth using organic diamine additives with a range of chain lengths and concentration. We found that ZnO branches form for a narrow diamine concentration range with a critical lower and upper critical nucleation concentration limit, which increases by about a factor of 5 for each additional carbon in the diaminoalkane chain. Our results suggest that the narrow window for secondary growth is dictated by the solubility of the ZnO crystals, where the low critical nucleation concentration is determined by slight etching of the surface to produce new nucleation sites, and the upper critical concentration is determined by the supersaturation concentration. Kinetic measurements show that the induction time and growth rate increase with increasing diamine concentration and follow classical nucleation and growth theory. Observations of branch morphological evolution reveal the mechanisms guiding the tunable crystal size and morphology.     

Kinetics of electrochemical nucleation and growth

A united scheme for the kinetics of electrochemical nucleation and the growth of a new phase is presented. The peculiarities of ion-transfer kinetics during electrochemical phase formation are analysed. The influence of the exchange current density at the electrolyte/cluster of the new phase interface on the nucleation rate, the nucleation induction time and the growth rate is reported.     

The nucleation rate of crystalline ice in amorphous solid water

The kinetics of crystalline ice nucleation and growth in nonporous, molecular beam deposited amorphous solid water (ASW) films are investigated at temperatures near 140 K. We implement an experimental methodology and corresponding model of crystallization kinetics to decouple growth from nucleation and quantify the temperature dependence and absolute rates of both processes. Nucleation rates are found to increase from approximately 3x10(13) m(-3) s(-1) at 134 K to approximately 2x10(17) m(-3) s(-1) at 142 K, corresponding to an Arrhenius activation energy of 168 kJ/mol. Over the same temperature range, the growth velocity increases from approximately 0.4 to approximately 4 A s(-1), also exhibiting Arrhenius behavior with an activation energy of 47 kJ/mol. These nucleation rates are up to ten orders of magnitude larger than in liquid water near 235 K, while growth velocities are approximately 10(9) times smaller. Crystalline ice nucleation kinetics determined in this study differ significantly from those reported previously for porous, background vapor deposited ASW, suggesting the nucleation mechanism is dependent upon film morphology.