Direct numerical simulation of homogeneous nucleation and growth in a phase-field model using cell dynamics method

The homogeneous nucleation and growth in a simplest two-dimensional phase field model is numerically studied using the cell dynamics method. The whole process from nucleation to growth is simulated and is shown to follow closely the Kolmogorov-Johnson-Mehl-Avrami (KJMA) scenario of phase transformation. Specifically the time evolution of the volume fraction of new stable phase is found to follow closely the KJMA formula. By fitting the KJMA formula directly to the simulation data, not only the Avrami exponent but the magnitude of nucleation rate and, in particular, of incubation time are quantitatively studied. The modified Avrami plot is also used to verify the derived KJMA parameters. It is found that the Avrami exponent is close to the ideal theoretical value m=3. The temperature dependence of nucleation rate follows the activation-type behavior expected from the classical nucleation theory. On the other hand, the temperature dependence of incubation time does not follow the exponential activation-type behavior. Rather the incubation time is inversely proportional to the temperature predicted from the theory of Shneidman and Weinberg [J. Non-Cryst. Solids 160, 89 (1993)]. A need to restrict thermal noise in simulation to deduce correct Avrami exponent is also discussed.     

Computational modelling of nematic phase ordering by film and droplet growth over heterogeneous substrates

This paper presents a computational study of defect nucleation associated with the kinetics of the isotropic-to-nematic phase ordering transition over heterogeneous substrates, as it occurs in new liquid crystal biosensor devices, based on the Landau-de Gennes model for rod-like thermotropic nematic liquid crystals. Two regimes are identified due to interfacial tension inequalities: (i) nematic surface film nucleation and growth normal to the heterogeneous substrate, and (ii) nematic surface droplet nucleation and growth. The former, known as wetting regime, leads to interfacial defect shedding at the moving nematic-isotropic interface. The latter droplet regime, involves a moving contact line, and exhibits two texturing mechanisms that also lead to interfacial defect shedding: (a) small and large contact angles of drops spreading over a heterogeneous substrate, and (b) small drops with large curvature growing over homogeneous patches of the substrate. The numerical results are consistent with qualitative defect nucleation models based on the kinematics of the isotropic-nematic interface and the substrate-nematic-isotropic contact line. The results extend current understanding of phase ordering over heterogeneous substrates by elucidating generic defect nucleation processes at moving interfaces and moving contact lines.     

Computational modelling of nematic phase ordering by film and droplet growth over heterogeneous substrates

This paper presents a computational study of defect nucleation associated with the kinetics of the isotropic‐to‐nematic phase ordering transition over heterogeneous substrates, as it occurs in new liquid crystal biosensor devices, based on the Landau–de Gennes model for rod‐like thermotropic nematic liquid crystals. Two regimes are identified due to interfacial tension inequalities: (i) nematic surface film nucleation and growth normal to the heterogeneous substrate, and (ii) nematic surface droplet nucleation and growth. The former, known as wetting regime, leads to interfacial defect shedding at the moving nematic‐isotropic interface. The latter droplet regime, involves a moving contact line, and exhibits two texturing mechanisms that also lead to interfacial defect shedding: (a) small and large contact angles of drops spreading over a heterogeneous substrate, and (b) small drops with large curvature growing over homogeneous patches of the substrate. The numerical results are consistent with qualitative defect nucleation models based on the kinematics of the isotropic–nematic interface and the substrate–nematic–isotropic contact line. The results extend current understanding of phase ordering over heterogeneous substrates by elucidating generic defect nucleation processes at moving interfaces and moving contact lines.     

Cell dynamics simulation of droplet and bridge formation within striped nanocapillaries

The kinetics of droplet and bridge formation within striped nanocapillaries is studied when the wetting film grows via interface-limited growth. The phenomenological time-dependent Ginzburg-Landau (TDGL)-type model with thermal noise is used and numerically solved using the cell dynamics method. The model is two-dimensional and consists of undersaturated vapor confined within a nanocapillary made of two infinitely wide flat substrates. The surface of the substrate is chemically heterogeneous with a single stripe of lyophilic domain that exerts long-range attractive potential to the vapor molecule. The dynamics of nucleation and subsequent growth of droplet and bridge can be simulated and visualized. In particular, the evolution of the morphology from droplet or bump to bridge is clearly identified. The crucial role played by the substrate potential on the morphology of bridge of nanoscopic size is clarified. Nearly temperature-independent evolution of capillary condensation is predicted when the interface-limited growth dominates. In addition, it is shown that the dynamics of capillary condensation follows the scenario of capillary condensation proposed by Everett and Haynes three decades ago.     

Kolmogorov–Johnson–Mehl–Avrami kinetics for non-isothermal phase transformations ruled by diffusional growth

We report on the functional form of the rate of the transformed volume fraction in non-isothermal phase transitions occurring by nucleation and diffusional growth. The microscopic growth rate is computed by solving the diffusion problem for time-dependent diffusion coefficient. The growth law is further employed in the Kolmogorov–Johnson–Mehl–Avrami (KJMA) theory for describing the time dependence of the transformed volume at constant heating rate. It is demonstrated that the transformation rate separates in the product of volume fraction and actual temperature functions. In the framework of the KJMA approach this factorization is exact. It is also shown that for real systems (due to the high values of the reduced activation energies for nucleation and growth), the kinetics is in excellent agreement with the stretched exponential function appropriate for isothermal transformations.     

自由表面的Ni原子团簇的熔化

采用分子动力学模拟技术研究了不同尺寸的Ni原子团簇的熔化过程.团簇的最初构型为FCC结构.研究结果表明,原子团簇的熔化温度与原子团簇中原子的个数有关,团簇的熔化首先从表面开始,当外层原子成为液态后,整个团簇的熔化从液态层开始,直至核心区域.该熔化过程可以被称为非均质熔化,自由表面充当非均质形核位置.作为对比,对无自由表面的大块固态Ni的熔化过程也进行了模拟,其熔化温度高于实验温度约400 K.表明对无自由表面的大块固态的熔化过程,液相形成无非均质形核位置,熔化的本质过程受均质形核机理控制.     

Nucleation in binary polymer blends: effects of foreign mesoscopic spherical particles

We study nucleation in binary polymer blends in the presence of mesoscopic spherical particles using self-consistent field theory, considering both heterogeneous and homogeneous nucleation mechanisms. Heterogeneous nucleation is found to be highly sensitive to surface selectivity and particle size, with rather subtle dependence on the particle size. Particles that preferentially adsorb the nucleating species generally favor heterogeneous nucleation. For sufficiently strong adsorption, barrierless nucleation is possible. By comparing the free energy barrier for homogeneous and heterogeneous nucleation, we construct a kinetic phase diagram.     

Heterogeneous and homogeneous nucleation compared: rapid nucleation on microscopic impurities

We use computer simulation to calculate the rates of both homogeneous nucleation and heterogeneous nucleation on microscopic impurities. We do so in perhaps the simplest model of fluids and magnets: the two-dimensional Ising model. We expect our results to be qualitatively applicable to many simple and complex fluids. We find that heterogeneous nucleation on an impurity that is not only microscopic but also as small as possible, that is, a single fixed spin, is more than four orders of magnitude faster than homogeneous nucleation. The rate of heterogeneous nucleation then increases by a factor of approximately five for each additional fixed spin in the impurity. These results suggest that impurities as small as single molecules can result in homogeneous nucleation being irrelevant due to heterogeneous nucleation on these microscopic impurities being much faster.     

Two-component heterogeneous nucleation kinetics and an application to Mars

We develop a two-component heterogeneous nucleation model that includes exact calculation of the Stauffer-type [D. Stauffer, J. Aerosol Sci. 7, 319 (1976)] steady-state kinetic prefactor using the correct heterogeneous Zeldovich factor for a heterogeneous two-component system. The model, and a simplified version of it, is tested by comparing its predictions to experimental data for water-n-propanol nucleating on silver particles. The model is then applied to water-carbon dioxide system in Martian conditions, which has not been modeled before. Using the ideal mixture assumption, the model shows theoretical possibilities for two-component nucleation adjacent to the initial stages of one-component water nucleation, especially with small water vapor amounts. The numbers of carbon dioxide molecules in the critical cluster are small in the case of large water amounts (up to 300 ppm) in the gas phase, but larger when there is very little water vapor (1 ppm).     

Phase behavior and viscoelastic properties of trisilanolcyclohexyl-POSS at the air/water interface

A trisilanol polyhedral oligomeric silsesquioxane (POSS), trisilanolcyclohexyl-POSS (TCyP), has recently been reported to undergo a series of phase transitions from traditional Langmuir monolayers to unique rodlike hydrophobic aggregates in multilayer films that are different from "collapsed" morphologies seen in other systems at the air/water interface. This paper focuses on the phase transitions and morphology of films varying in average thickness from monolayers to trilayers and the corresponding viscoelastic properties of trisilanolcyclohexyl-POSS molecules at the air/water interface by means of surface pressure-area per molecule (Pi-A) isotherms, Brewster angle microscopy (BAM), and interfacial stress rheometry (ISR) measurements. The morphology studies by BAM reveal that the TCyP monolayer can collapse into different 3D structures by homogeneous or heterogeneous nucleation mechanisms. For homogeneous nucleation, analysis by Vollhardt et al.'s nucleation and growth model reveals that TCyP collapse is consistent with instantaneous nucleation with hemispherical edge growth at Pi = 3.7 mN.m(-1). Both surface storage (Gs') and loss (Gs") moduli obtained by ISR reveal three different non-Newtonian flow regimes that correlate with phase transitions in the Pi-A isotherms: (A) A viscous liquidlike "monolayer"; (B) a "biphasic regime"between a liquidlike viscous monolayer and a more rigid trilayer; and (C) an elastic solidlike "trilayer". These observations provide interesting insights into collapse mechanisms and structures in Langmuir films.     

Homogeneous nucleation of droplets from supersaturated vapor in a closed system

Kinetic equations describing homogeneous nucleation kinetics within standard model are solved numerically under the condition of a constant number of molecules in the considered system. It has consequences to decrease the supersaturation of the supersaturated vapor during the process of the formation of small droplets of a new phase. The decrease of supersaturation occurs in a short time and reaches some value which remains unchanged for a relatively long time (quasistationary regime), especially at lower initial supersaturations. This time interval decreases with increasing value of the initial supersaturation. In the quasistationary regime the nucleation rate reaches its stationary value. At higher initial supersaturation, the rate of formation of nuclei goes to some maximum value corresponding to the stationary nucleation rate and then decreases with time due to the decrease of supersaturation.     

Molecular dynamics analysis of multiple site growth and coalescence effects on homogeneous and heterogeneous nucleations

Homogeneous and heterogeneous nucleations were simulated by molecular dynamics (MD). The behavior of Lennard-Jones molecules was studied inside a liquid-gas system where all dimensions of the wall were periodic and a soft core carrier gas within the system controlled the temperature. In this study, the classical nucleation theory was found to underestimate the homogeneous nucleation rate by five orders of magnitude, which complies with other MD studies. The discrepancy in the nucleation rate between theory and simulation was mainly caused by the fundamental assumption that there are no volumetric interactions in the growth process. In this particular case, however, growth was observed at multiple sites due to Ostwald ripening and coalescence between nuclei by Brownian motion. Furthermore, even though the supersaturation ratio is inadequate for homogeneous nucleation, once a seed is introduced to the system, a cluster can be created. The addition of seeds not only enhances nucleation but also renders coalescence as an important nucleation mechanism in the earlier stages compared to homogeneous nucleation.     

Heterogeneous critical nucleation on a completely wettable substrate

Heterogeneous nucleation of a new bulk phase on a flat substrate can be associated with the surface phase transition called wetting transition. When this bulk heterogeneous nucleation occurs on a completely wettable flat substrate with a zero contact angle, the classical nucleation theory predicts that the free-energy barrier of nucleation vanishes. In fact, there always exists a critical nucleus and a free-energy barrier as the first-order prewetting transition will occur even when the contact angle is zero. Furthermore, the critical nucleus changes its character from the critical nucleus of surface phase transition below bulk coexistence (undersaturation) to the critical nucleus of bulk heterogeneous nucleation above the coexistence (oversaturation) when it crosses the coexistence. Recently, Sear [J. Chem. Phys. 129, 164510 (2008)] has shown, by a direct numerical calculation of nucleation rate, that the nucleus does not notice this change when it crosses the coexistence. In our work, the morphology and the work of formation of critical nucleus on a completely wettable substrate are re-examined across the coexistence using the interface-displacement model. Indeed, the morphology and the work of formation changes continuously at the coexistence. Our results support the prediction of Sear and will rekindle the interest on heterogeneous nucleation on a completely wettable substrate.     

Nonhydrolytic synthesis of high-quality anisotropically shaped brookite TiO2 nanocrystals

A surfactant-assisted nonaqueous strategy, relying on high-temperature aminolysis of titanium carboxylate complexes, has been developed to access anisotropically shaped TiO2 nanocrystals selectively trapped in the metastable brookite phase. Judicious temporal manipulation of precursor supply to the reaction mixture enables systematic tuning of the nanostructure geometric features over an exceptionally wide dimensional range (30-200 nm). Such degree of control is rationalized within the frame of a self-regulated phase-changing seed-catalyzed mechanism, in which homogeneous nucleation, on one side, and heterogeneous nucleation/growth processes, on the other side, are properly balanced while switching from the anatase to the brookite structures, respectively, in a continuous unidirectional crystal development regime. The time variation of the chemical potential for the monomer species in the solution, the size dependence of thermodynamic structural stability of the involved titania polymorphs, and the reduced activation barrier for brookite nucleation onto initially formed anatase seeds play decisive roles in the crystal-phase- and shape-tailored growth of titania nanostructures by the present approach.     

Inhomogeneous phase separation kinetics in liquid binary mixtures: Sensitivity to initial local composition 2

The kinetics of phase separation subsequent to a finite temperature quench is assumed to be driven by diffusion on the altered free energy surface and is generally assumed to be slow. The situation can be different in phase separating liquid binary mixtures, especially for systems characterized by the large difference in mutual interactions between solute and solvent molecules. In such cases, the phase separation kinetics could be fast and may get completed within a short time (ns) scale. As a result, in these systems, one may observe diverse dynamical features arising out of local heterogeneity leading to the onset of phase separation through pattern formation, spinodal decomposition, nucleation, and growth. By using a coarse-grained analysis, we examine phase separation kinetics in each spatial grid and indeed observe important effects of initial heterogeneity on the subsequent evolution. Interestingly, we observe slower separation kinetics for those regions that correspond to the composition at the minimum of the high-temperature surface. The heterogeneous dynamics has been captured here through the non-linear susceptibility function, which shows a pattern similar to what is observed in the supercooled liquid. Each grid shows somewhat different dynamics in the three-stage (exponential, power-law, and logarithmic regime) phase separation dynamics. The late stage of phase separation kinetics is usually attributed to the coarsening of the phase-separated domains. However, in a liquid binary mixture, the late-stage power-law decay undergoes a further change. A new dynamical regime arises characterized by a logarithmic time dependence, which is due to the “smoothening” of the rough interface of already well-separated phases. This can also be described as opposite to the roughening transition described by Chui and Weeks [Phys. Rev. Lett. 40, 733 (1978)]. This reverse roughening transition can explain the logarithmic time dependence observed in the simulation.     

On the influence of a patterned substrate on crystallization in suspensions of hard spheres

We present a computer simulation study on crystal nucleation and growth in supersaturated suspensions of mono-disperse hard spheres induced by a triangular lattice substrate. The main result is that compressed substrates are wet by the crystalline phase (the crystalline phase directly appears without any induction time), while for stretched substrates we observe heterogeneous nucleation. The shapes of the nucleated crystallites fluctuate strongly. In the case of homogeneous nucleation amorphous precursors have been observed [T. Schilling et al., Phys. Rev. Lett. 105(2), 025701 (2010)]. For heterogeneous nucleation we do not find such precursors. The fluid is directly transformed into highly ordered crystallites.     

Turbulent fronts in resonantly forced oscillatory systems

Phase fronts in the forced complex Ginzburg-Landau equation, a model of a resonantly forced oscillatory reaction-diffusion system, are studied in the 3:1 resonance regime. The focus is on the turbulent (Benjamin-Feir-unstable) regime of the corresponding unforced system; in the forced system, phase fronts between spatially uniform phase-locked states exhibit complex dynamics. In one dimension, for strong forcing, phase fronts move with constant velocity. As the forcing intensity is lowered there is a bifurcation to oscillatory motion, followed by a bifurcation to a regime in which fronts multiply via the nucleation of domains of the third homogeneous phase in the front. In two dimensional systems, rough fronts with turbulent, complex internal structure may arise. For a critical value of the forcing intensity there is a nonequilibrium phase transition in which the turbulent interface grows to occupy the entire system. The phenomena we explore can be probed by experiments on periodically forced light sensitive reaction-diffusion systems.     

Nucleation in the presence of slow microscopic dynamics

Nucleation of a new thermodynamic phase is often a slow process due to the need to overcome a high free-energy barrier. However, there are other sources of slow dynamics; for example, at high densities/low temperatures, the movement of individual molecules or spins may be slow. Here, we study nucleation in a simple phenomenological model that has this type of slow microscopic dynamics. We do this to better understand how the two sources of slow dynamics interact. We find that as nucleation is intrinsically slow, only very slow microscopic dynamics strongly affect how nucleation occurs. The composition of the nucleus at the top of the nucleation barrier is much less sensitive to slow microscopic dynamics than is the composition of the nucleus once it is postcritical. However, slow dynamics affects not only the rate but also the pathway, which no longer goes over the saddle point in the free energy. We also find that the slow microscopic dynamics can cause sampling problems in an algorithm developed to calculate nucleation rates, and so cause it to predict the rate incorrectly.     

Multiscale approach to CO2 hydrate formation in aqueous solution: phase field theory and molecular dynamics. Nucleation and growth

A phase field theory with model parameters evaluated from atomistic simulations/experiments is applied to predict the nucleation and growth rates of solid CO(2) hydrate in aqueous solutions under conditions typical to underwater natural gas hydrate reservoirs. It is shown that under practical conditions a homogeneous nucleation of the hydrate phase can be ruled out. The growth rate of CO(2) hydrate dendrites has been determined from phase field simulations as a function of composition while using a physical interface thickness (0.85+/-0.07 nm) evaluated from molecular dynamics simulations. The growth rate extrapolated to realistic supersaturations is about three orders of magnitude larger than the respective experimental observation. A possible origin of the discrepancy is discussed. It is suggested that a kinetic barrier reflecting the difficulties in building the complex crystal structure is the most probable source of the deviations.     

高分子体系相分离动力学及图样生成和选择

高分子共混物的混合熵很小导致其多为热力学不相容体系而发生相分离,形成特定的时空图样。本文根据多年来我们自己的研究工作并结合实例,基于时间分辨的Ginzberg-Landau 方法研究高分子复杂体系相分离动力学及图样选择,重点介绍剪切外场下高分子共混物及嵌段高分子的相分离,耦合化学反应的相分离,在弯曲曲面特别是球面上的相分离,以及TDGL与密度泛函理论的有机结合即动态自洽场理论在具有不同链拓朴结构的嵌段高分子体系中研究相分离动力学。