Time dependent nucleation

Anomalously abrupt switching of wurtzite-structured ferroelectrics: simultaneous non-linear nucleation and growth model

Ferroelectric polarization switching is one common example of a process that occurs via nucleation and growth, and understanding switching kinetics is crucial for applications such as ferroelectric memory. Here we describe and interpret anomalous switching dynamics in the wurtzite-structured nitride thin film ferroelectrics Al0.7Sc0.3N and Al0.94B0.06N using a general model that can be directly applied to other abrupt transitions that proceed via nucleation and growth. When substantial growth and impingement occur while nucleation rate is increasing, such as in these wurtzite-structured ferroelectrics under high electric fields, abrupt polarization reversal leads to very large Avrami coefficients (e.g., n = 11), inspiring an extension of the KAI (Kolmogorov-Avrami-Ishibashi) model. We apply this extended model to two related but distinct scenarios that crossover between (typical) behavior described by sequential nucleation and growth and a more abrupt transition arising from significant growth prior to peak nucleation rate. This work therefore provides a more complete description of general nucleation and growth kinetics applicable to any system while specifically addressing the anomalously abrupt polarization reversal behavior in new wurtzite-structured ferroelectrics.

Fokker-Planck equation for Coulomb relaxation and wave-particle diffusion: Spectral solution and the stability of the Kappa distribution to Coulomb collisions

The present paper considers the time evolution of a charged test particle of mass m in a constant temperature heat bath of a second charged particle of mass M. The time dependence of the distribution function of the test particles is given by a Fokker-Planck equation with a diffusion coefficient for Coulomb collisions as well as a diffusion coefficient for wave-particle interactions. For the mass ratio m/M→0, the steady distribution is a Kappa distribution which has been employed in space physics to fit observed particle energy spectra. The time dependence of the distribution functions with some initial value is expressed in terms of the eigenvalues and eigenfunctions of the linear Fokker-Planck operator and also interpreted with the transformation to a Schrödinger equation. We also consider the explicit time dependence of the distribution function with a discretization of the Fokker-Planck equation. We study the stability of the Kappa distribution to Coulomb collisions.

Derivation of an exact, nonequilibrium framework for nucleation: Nucleation is a priori neither diffusive nor Markovian

We discuss the structure of the equation of motion that governs nucleation processes at first order phase transitions. From the underlying microscopic dynamics of a nucleating system, we derive by means of a nonequilibrium projection operator formalism the equation of motion for the size distribution of the nuclei. The equation is exact, i.e., the derivation does not contain approximations. To assess the impact of memory, we express the equation of motion in a form that allows for direct comparison to the Markovian limit. As a numerical test, we have simulated crystal nucleation from a supersaturated melt of particles interacting via a Lennard-Jones potential. The simulation data show effects of non-Markovian dynamics.

Diffusion in a bistable system: The eigenvalue spectrum of the Fokker-Planck operator and Kramers' reaction rate theory

The time-dependent solution of the Fokker-Planck equation with bistable potentials is considered in terms of the eigenfunctions and eigenvalues of the linear Fokker-Planck operator. The Fokker-Planck equation is the high friction limit of the corresponding Kramers' equation. Two different potentials are considered defined with a constant diffusion coefficient, ε, and position-dependent drift coefficients. The smallest nonzero eigenvalue of the Fokker-Planck operator, λ_{1}, provides the long-time rate coefficient for the transformation of the different species in the two stable states. A novel pseudospectral method with nonclassical polynomials is applied to this class of systems. The convergence of the eigenvalues and eigenfunctions of the Fokker-Planck operator versus the number of basis functions is studied and compared with previous results. The results are consistent with Kramers' theory, and a linear relationship between lnλ_{1} and 1/ε for sufficiently small ε values is verified. A comparison with analytic approximations to λ_{1} is provided.

Kappa and other nonequilibrium distributions from the Fokker-Planck equation and the relationship to Tsallis entropy

This paper considers two nonequilibrium model systems described by linear Fokker-Planck equations for the time-dependent velocity distribution functions that yield steady state Kappa distributions for specific system parameters. The first system describes the time evolution of a charged test particle in a constant temperature heat bath of a second charged particle. The time dependence of the distribution function of the test particle is given by a Fokker-Planck equation with drift and diffusion coefficients for Coulomb collisions as well as a diffusion coefficient for wave-particle interactions. A second system involves the Fokker-Planck equation for electrons dilutely dispersed in a constant temperature heat bath of atoms or ions and subject to an external time-independent uniform electric field. The momentum transfer cross section for collisions between the two components is assumed to be a power law in reduced speed. The time-dependent Fokker-Planck equations for both model systems are solved with a numerical finite difference method and the approach to equilibrium is rationalized with the Kullback-Leibler relative entropy. For particular choices of the system parameters for both models, the steady distribution is found to be a Kappa distribution. Kappa distributions were introduced as an empirical fitting function that well describe the nonequilibrium features of the distribution functions of electrons and ions in space science as measured by satellite instruments. The calculation of the Kappa distribution from the Fokker-Planck equations provides a direct physically based dynamical approach in contrast to the nonextensive entropy formalism by Tsallis [J. Stat. Phys. 53, 479 (1988)JSTPBS0022-471510.1007/BF01016429].

Transient effects in ice nucleation of a water drop impacting onto a cold substrate

The impact of water drops onto a solid surface at subfreezing temperatures has been experimentally studied. Drop nucleation has been observed using a high-speed video system. The statistics of nucleation allows the estimation of the average number of nucleation sites per unit area of the wetted part of the substrate. We have discovered that the nucleation rate in the impacting drop is not constant. The observed significant increase of the nucleation rate at small times after impact t<50 ms can be explained by the generation of nanobubbles at early times of drop impact. These bubbles serve as additional nucleation sites and enhance the nucleation rate.

Unification of classical nucleation theories via a unified Itô-Stratonovich stochastic equation

Classical nucleation theory (CNT) is the most widely used framework to describe the early stage of first-order phase transitions. Unfortunately, the different points of view adopted to derive it yield different kinetic equations for the probability density function, e.g., Zeldovich-Frenkel or Becker-Döring-Tunitskii equations. Starting from a phenomenological stochastic differential equation, a unified equation is obtained in this work. In other words, CNT expressions are recovered by selecting one or another stochastic calculus. Moreover, it is shown that the unified CNT thus obtained produces the same Fokker-Planck equation as that from a recent update of CNT [J. F. Lutsko and M. A. Durán-Olivencia, J. Chem. Phys. 138, 244908 (2013)10.1063/1.4811490] when mass transport is governed by diffusion. Finally, we derive a general induction-time expression along with specific approximations of it to be used under different scenarios, in particular, when the mass-transport mechanism is governed by direct impingement, volume diffusion, surface diffusion, or interface transfer.

Mesoscopic nucleation theory for confined systems: a one-parameter model

Classical nucleation theory has been recently reformulated based on fluctuating hydrodynamics [J. F. Lutsko and M. A. Durán-Olivencia, Classical nucleation theory from a dynamical approach to nucleation, J. Chem. Phys. 138, 244908 (2013). The present work extends this effort to the case of nucleation in confined systems such as small pores and vesicles. The finite available mass imposes a maximal supercritical cluster size and prohibits nucleation altogether if the system is too small. We quantity the effect of system size on the nucleation rate. We also discuss the effect of relaxing the capillary-model assumption of zero interfacial width resulting in significant changes in the nucleation barrier and nucleation rate.

A novel approach to the theory of homogeneous and heterogeneous nucleation

A new approach to the theory of nucleation, formulated relatively recently by Ruckenstein, Narsimhan, and Nowakowski (see Refs. [7-16]) and developed further by Ruckenstein and other colleagues, is presented. In contrast to the classical nucleation theory, which is based on calculating the free energy of formation of a cluster of the new phase as a function of its size on the basis of macroscopic thermodynamics, the proposed theory uses the kinetic theory of fluids to calculate the condensation (W(+)) and dissociation (W(-)) rates on and from the surface of the cluster, respectively. The dissociation rate of a monomer from a cluster is evaluated from the average time spent by a surface monomer in the potential well as obtained from the solution of the Fokker-Planck equation in the phase space of position and momentum for liquid-to-solid transition and the phase space of energy for vapor-to-liquid transition. The condensation rates are calculated using traditional expressions. The knowledge of those two rates allows one to calculate the size of the critical cluster from the equality W(+)=W(-) as well as the rate of nucleation. The developed microscopic approach allows one to avoid the controversial application of classical thermodynamics to the description of nuclei which contain a few molecules. The new theory was applied to a number of cases, such as the liquid-to-solid and vapor-to-liquid phase transitions, binary nucleation, heterogeneous nucleation, nucleation on soluble particles and protein folding. The theory predicts higher nucleation rates at high saturation ratios (small critical clusters) than the classical nucleation theory for both solid-to-liquid as well as vapor-to-liquid transitions. As expected, at low saturation ratios for which the size of the critical cluster is large, the results of the new theory are consistent with those of the classical one. The present approach was combined with the density functional theory to account for the density profile in the cluster. This approach was also applied to protein folding, viewed as the evolution of a cluster of native residues of spherical shape within a protein molecule, which could explain protein folding/unfolding and their dependence on temperature.

Quantum-mechanical solution to fundamental problems of classical theory of water vapor nucleation

The inconsistent temperature dependence of nucleation rates, disagreement of theoretical critical or onset supersaturations with experimental data, and insufficiently accurate predictions of nucleation rates are fundamental problems of the classical nucleation theory (CNT) of water vapors, which is a foundation of various multicomponent nucleation models widely used in the aerosol microphysics, physical chemistry, and chemical technology. In the present study, a correction to the CNT obtained from "first principles" has been derived and significant progress has been made in solving the fundamental problem of predicting nucleation rates of water vapors. The modified model with the quantum-mechanical correction incorporated is in very good agreement with experiments over the full range of temperatures (T=210-290K) , saturation ratios (S=2-100) , and nucleation rates (J= approximately 10{1}-10{17} cm-3).

Nucleation and size distribution of nucleus during induction period of polyethylene crystallization

The crystallization process from supercooled melt results in the formation of nanosize nuclei in the earlier stage (induction period) through subsequent attachment or detachment of repeating unit to nuclei. The size distribution of nucleus f(N(j),t) in the induction period of nucleation process from the melts has not been experimentally confirmed yet by direct observation. The reason is that the number density of nuclei nu is too small to be detected experimentally. In our previous work, we showed the direct evidence of nucleation experimentally by means of small angle x-ray scattering (SAXS) technique. Further we have succeeded to observe the nucleation and f(N(j),t) of polymer crystallization from the melts by SAXS using synchrotron radiation. We increased nu by adding a nucleating agent to a polymer (polyethylene). The time evolution of f(N(j),t) was observed for the first time.

Collapse of transient nucleation fluxes in a cold Ising ferromagnet

We report the time-dependent nucleation fluxes and associated nucleation rates in a metastable Ising ferromagnet on square lattice with Metropolis (Glauber-type) dynamics. It is discovered that, with lowering of the temperature, fluxes collapse into several representative transient curves corresponding to magic cluster sizes. Those can be associated with physical droplets, i.e., long-lived configurations which provide a link with the classical Becker-Döring picture.