Heterogeneous condensation of the Lennard-Jones vapor onto a nanoscale seed particle

The heterogeneous condensation of a Lennard-Jones vapor onto a nanoscale seed particle is studied using molecular dynamics simulations. Measuring the nucleation rate and the height of the free energy barrier using the mean first passage time method shows that the presence of a weakly interacting seed has little effect on the work of forming very small cluster embryos but accelerates the rate by lowering the barrier for larger clusters. We suggest that this results from a competition between the energetic and entropic features of cluster formation in the bulk and at the heterogeneity. As the interaction is increased, the free energy of formation is reduced for all cluster sizes. We also develop a simple phenomenological model of film formation on a small seed that captures the general features of the nucleation process for small heterogeneities. A comparison of our simulation results with the model shows that heterogeneous classical nucleation theory provides a good estimate of the critical size of the film but significantly overestimates the size of the barrier.     

Kinetics of ion-induced nucleation in a vapor-gas mixture

A general solution for the steady-state ion-induced nucleation kinetics has been derived, considering the differences between ion-induced nucleation and homogeneous nucleation. This solution includes a new effect for nucleation kinetics, the interaction of charged clusters with vapor molecules. Analytical expressions for the ion-induced nucleation rate have been obtained for the limiting cases of high and low thermodynamic barriers. The physical explanation of the so-called sign effect is proposed based on multipole expansion of an electric field of the cluster ion. This theory gives good agreement with experiments and is used to elucidate experimentally observed phenomena.     

Hydration of the bisulfate ion: atmospheric implications

Using molecular dynamics configurational sampling combined with ab initio energy calculations, we determined the low energy isomers of the bisulfate hydrates. We calculated the CCSD(T) complete basis set (CBS) binding electronic and Gibbs free energies for 53 low energy isomers of HSO(4)(-)(H(2)O)(n=1-6) and derived the thermodynamics of adding waters sequentially to the bisulfate ion and its hydrates. Comparing the HSO(4)(-)/H(2)O system to the neutral H(2)SO(4)/H(2)O cluster, water binds more strongly to the anion than it does to the neutral molecules. The difference in the binding thermodynamics of HSO(4)(-)/H(2)O and H(2)SO(4)/H(2)O systems decreases with increasing number of waters. The thermodynamics for the formation of HSO(4)(-)(H(2)O)(n=1-5) is favorable at 298.15 K, and that of HSO(4)(-)(H(2)O)(n=1-6) is favorable for T < 273.15 K. The HSO(4)(-) ion is almost always hydrated at temperatures and relative humidity values encountered in the troposphere. Because the bisulfate ion binds more strongly to sulfuric acid than it does to water, it is expected to play a role in ion-induced nucleation by forming a strong complex with sulfuric acid and water, thus facilitating the formation of a critical nucleus.     

Measurement of the thermodynamics of the hydrated dimer and trimer of sulfuric acid

We report measurements of the equilibrium constants, K(p), for the water mediated clustering of two and three sulfuric acid molecules for atmospheric temperatures and relative humidities. Limits for K(p) for the formation of the sulfuric acid tetramer and higher clusters as well as the kinetics of small sulfuric acid cluster growth are also presented. These results put strong constraints on the rate of nucleation in the atmosphere via the H(2)SO(4)/H(2)O neutral mechanism. We show that the neutral nucleation of H(2)SO(4) and H(2)O is slower than ion-induced nucleation of H(2)SO(4) and H(2)O for most conditions found in the middle and upper troposphere. These laboratory-based upper limits to the nucleation rates are also much lower than the predictions of the liquid drop model/classical nucleation theory.     

A limit of stability in supercooled liquid clusters

We examine the metastable liquid phase of a supercooled gold nanocluster by studying the free energy landscape using the largest solidlike embryo as an order parameter. Just below freezing, the free energy exhibits a local minimum at small embryo sizes and a maximum at a larger critical embryo size. At T=660 K the free energy becomes a monotonically decreasing function of the order parameter as the liquid phase becomes unstable, indicating that we have reached a limit of stability. In contrast to the mean-field theory predictions for a spinodal, the size of the critical embryo remains finite as the limit of stability is approached. We also calculate the rate of nucleation, independently from our free energy calculations, and observe a rapid increase in its temperature dependence when the free energy barrier is on the order of kT. We suggest that this supports the idea that freezing becomes a barrierless process at low temperatures.     

Molecular Dynamics Studies of the Kinetics of Phase Changes in Clusters IV： Crystal Nucleation from Molten （NaCl）256 and （NaCl）500 Clusters

Molecular dynamics computer simulation based on the Born-Mayer-Huggins potential function has been carried out to study the effects of duster size and temperature on the nucleation rate of sodium chloride dusters in the temperature range of 580 K to 630 K. Clusters with 256 and 500 NaCl molecules have been studied and the results have been compared with those obtained from 108 molecule dusters. The melting point (MP) of the clusters were observed to increase with the size of the clusters and can be well described by a linear equation MP =1107(37)-1229(23)N^-1/3(N is the number of molecules in the duster).The nucleation rate was found to decrease with increasing the duster size or temperature. Various nucleation theories have been used to interpret the nucleation rates obtained from this molecular dynamics simulation. It is possible to use a constant diffuse interface thickness to interpret the nucleation rate from the diffuse interface theory in the temperature range of this study. However, the interfacinl free energy estimated from classical nucleation theory and diffuse interface theory increases too fast with increasing the temperature while that from Gran-Gunton theory does not change with changing temperatures.The sizes of critical nuclei estimated from all the theories are smaller than those estimated from our simulations.     

Effect of chlorine ions on the stability of nucleation cores in condensing water vapors

The Gibbs energy and the equilibrium work of nucleation of condensed-phase nuclei in water vapors in the presence of chlorine anions are calculated at the molecular level using the Monte Carlo method at a temperature of 273 K for sizes of up to 200 molecules. It is found that the growth of a nucleus leads to the displacement of a chloride ion onto its surface, which is accompanied by a loss of the thermodynamic stability of the microdroplet. It is found that the size dependence of the work of nucleation exhibits a minimum and a maximum with an inflection point that separates the region of stable and unstable equilibrium sizes.     

Surface tension and scaling of critical nuclei in diatomic and triatomic fluids

Density functional theory has been used to investigate surface tension and scaling of critical clusters in fluids consisting of diatomic and rigid triatomic molecules. The atomic sites are hard spheres with attractive interactions obtained from the tail part of the Lennard-Jones potential. Asymmetry in attractive interactions between the atomic sites has been introduced to cause molecular orientation and oscillatory density profiles at liquid-vapor interfaces. The radial dependence of cluster surface tension in fluids showing modest orientation in unimolecular layer at the interface or no orientation at all resembles the surface tension behavior of clusters in simple monoatomic fluids, although the surface tension maximum becomes more pronounced with increasing chain length of the molecule. Surface tension of clusters having multiple oscillatory layers at the interface shows a prominent maximum at small cluster sizes; however, the surface tension of large clusters is lower than the planar value. The scaling relation for the number of molecules in the critical cluster and the nucleation barrier height developed by McGraw and Laaksonen [Phys. Rev. Lett. 76, 2754 (1996)] are well obeyed for fluids with little structure at liquid-vapor interface. However, fluids having enhanced interfacial structure show some deviation from the particle number scaling, and the barrier height scaling breaks up seriously.     

Recent developments in the kinetic theory of nucleation

A review of recent progress in the kinetics of nucleation is presented. In the conventional approach to the kinetic theory of nucleation, it is necessary to know the free energy of formation of a new-phase particle as a function of its independent variables at least for near-critical particles. Thus the conventional kinetic theory of nucleation is based on the thermodynamics of the process. The thermodynamics of nucleation can be examined by using various approaches, such as the capillarity approximation, density functional theory, and molecular simulation, each of which has its own advantages and drawbacks. Relatively recently a new approach to the kinetics of nucleation was proposed [Ruckenstein E, Nowakowski B. J Colloid Interface Sci 1990;137:583; Nowakowski B, Ruckenstein E. J Chem Phys 1991;94:8487], which is based on molecular interactions and does not employ the traditional thermodynamics, thus avoiding such a controversial notion as the surface tension of tiny clusters involved in nucleation. In the new kinetic theory the rate of emission of molecules by a new-phase particle is determined with the help of a mean first passage time analysis. This time is calculated by solving the single-molecule master equation for the probability distribution function of a surface layer molecule moving in a potential field created by the rest of the cluster. The new theory was developed for both liquid-to-solid and vapor-to-liquid phase transitions. In the former case the single-molecule master equation is the Fokker-Planck equation in the phase space which can be reduced to the Smoluchowski equation owing to the hierarchy of characteristic time scales. In the latter case, the starting master equation is a Fokker-Planck equation for the probability distribution function of a surface layer molecule with respect to both its energy and phase coordinates. Unlike the case of liquid-to-solid nucleation, this Fokker-Planck equation cannot be reduced to the Smoluchowski equation, but the hierarchy of time scales does allow one to reduce it to the Fokker-Plank equation in the energy space. The new theory provides an equation for the critical radius of a new-phase particle which in the limit of large clusters (low supersaturations) yields the Kelvin equation and hence an expression for the macroscopic surface tension. The theory was illustrated with numerical calculations for a molecular pair interaction potential combining the dispersive attraction with the hard-sphere repulsion. The results for the liquid-to-solid nucleation clearly show that at given supersaturation the nucleation rate depends on the cluster structure (for three cluster structures considered-amorphous, fcc, and icosahedral). For both the liquid-to-solid and vapor-to-liquid nucleation, the predictions of the theory are consistent with the results of classical nucleation theory (CNT) in the limit of large critical clusters (low supersaturations). For small critical clusters the new theory provides higher nucleation rates than CNT. This can be accounted for by the fact that CNT uses the macroscopic interfacial tension which presumably overpredicts the surface tension of small clusters, and hence underpredicts nucleation rates.     

EXAMINATION OF THE SURFACE FREE ENERGY AND ACID-BASE PROPERTIES OF CELLULOSE BY THE COLUMN WICKING TECHNIQUE AND THE CRITICAL PACKING HEIGHT/DENSITY

INTRODUCTIONThe wicking technique has been used to determine the surface energy components of solids[1-7]. Generally, thereare two wicking processes in practice, i.e. thin layer wicking and column wicking. A comparison of the twowicking techniques made by van Oss et al.[2], shows that column wicking has a problem of showing differentpacking density due to some kinds of powder materials, i.e. clay can not give good repacking reproducibility[2].But thin layer wicking has a problem in sample…     

Sensitivity of nucleation phenomena on range of interaction potential

Theoretical and computational investigations of nucleation have been plagued by the sensitivity of the phase diagram to the range of the interaction potential. As the surface tension depends strongly on the range of interaction potential and as the classical nucleation theory (CNT) predicts the free energy barrier to be directly proportional to the cube of the surface tension, one expects a strong sensitivity of nucleation barrier to the range of the potential; however, CNT leaves many aspects unexplored. We find for gas-liquid nucleation in Lennard-Jones system that on increasing the range of interaction the kinetic spinodal (KS) (where the mechanism of nucleation changes from activated to barrierless) shifts deeper into the metastable region. Therefore the system remains metastable for larger value of supersaturation and this allows one to explore the high metastable region without encountering the KS. On increasing the range of interaction, both the critical cluster size and pre-critical minima in the free energy surface of kth largest cluster, at respective kinetic spinodals, shift towards smaller cluster size. In order to separate surface tension contribution to the increase in the barrier from other non-trivial factors, we introduce a new scaling form for surface tension and use it to capture both the temperature and the interaction range dependence of surface tension. Surprisingly, we find only a weak non-trivial contribution from other factors to the free energy barrier of nucleation.     

氟化钾团簇的相变、成核和重结晶的分子动力学模拟

我们利用Born-Mayer-Huggins相互作用势函数对(KF)_N(N=108，256，500和864)团簇进行了分子动力学(MD)模拟。为了避免周期性边界条件对相变、成核和重结晶的干扰作用，对体系采用了自由边界。基于MD模拟结果，对团簇的熔化温度、熔化焓、扩散系数、成核速率、固液界面自由能、临界核大小等进行了计算和讨论。在对(KF)₈₆₄双晶团簇的热退火模拟中，观察到了固态的重结晶和晶粒的生长。经典的成核理论成功地解释了(KF)₈₆₄双晶团簇的重结晶MD模拟结果。     

Global potential energy minima of C60(H2O)n clusters

Likely candidates for the global potential energy minima of C60(H2O)n clusters with n < or = 21 are found using basin-hopping global optimization. The potential energy surfaces are constructed using the TIP4P intermolecular potential for the water molecules, a Lennard-Jones water-fullerene potential, and a water-fullerene polarization potential, which depends on the first few nonvanishing C60 multipole polarizabilities. This combination produces a rather hydrophobic water-fullerene interaction. As a consequence, the water component of the lowest C60(H2O)n minima is quite closely related to low-lying minima of the corresponding TIP4P (H2O)n clusters. In most cases, the geometrical substructure of the water molecules in the C60(H2O)n global minimum coincides with that of the corresponding free water cluster. Exceptions occur when the interaction with C60 induces a change in geometry. This qualitative picture does not change significantly if we use the TIP3P model for the water-water interaction. Structures such as C60@(H2O)60, in which the water molecules surround the C60 fullerene, correspond to local minima with much higher potential energies. For such a structure to become the global minimum, the magnitude of the water-fullerene interaction must be increased to an unphysical value.     

Box size effects are negligible for solvation free energies of neutral solutes

Hydration free energy calculations in explicit solvent have become an integral part of binding free energy calculations and a valuable test of force fields. Most of these simulations follow the conventional norm of keeping edge length of the periodic solvent box larger than twice the Lennard-Jones (LJ) cutoff distance, with the rationale that this should be sufficient to keep the interactions between copies of the solute to a minimum. However, for charged solutes, hydration free energies can exhibit substantial box size-dependence even at typical box sizes. Here, we examine whether similar size-dependence affects hydration of neutral molecules. Thus, we focused on two strongly polar molecules with large dipole moments, where any size-dependence should be most pronounced, and determined how their hydration free energies vary as a function of simulation box size. In addition to testing a variety of simulation box sizes, we also tested two LJ cut-off distances, 0.65 and 1.0 nm. We show from these simulations that the calculated hydration free energy is independent of the box-size as well as the LJ cut-off distance, suggesting that typical hydration free energy calculations of neutral compounds indeed need not be particularly concerned with finite-size effects as long as standard good practices are followed.     

Nucleation near the spinodal: limitations of mean field density functional theory

We investigate the diverging size of the critical nucleus near the spinodal using the gradient theory (GT) of van der Waals and Cahn and Hilliard and mean field density functional theory (MFDFT). As is well known, GT predicts that at the spinodal the free energy barrier to nucleation vanishes while the radius of the critical fluctuation diverges. We show numerically that the scaling behavior found by Cahn and Hilliard for these quantities holds quantitatively for both GT and MFDFT. We also show that the excess number of molecules Deltag satisfies Cahn-Hilliard scaling near the spinodal and is consistent with the nucleation theorem. From the latter result, it is clear that the divergence of Deltag is due to the divergence of the mean field isothermal compressibility of the fluid at the spinodal. Finally, we develop a Ginzburg criterion for the validity of the mean field scaling relations. For real fluids with short-range attractive interactions, the near-spinodal scaling behavior occurs in a fluctuation dominated regime for which the mean field theory is invalid. Based on the nucleation theorem and on Wang's treatment of fluctuations near the spinodal in polymer blends, we infer a finite size for the critical nucleus at the pseudospinodal identified by Wang.     

Reductions with SmI2: mechanistic probe for distinguishing between two operational modes of proton transfer

The reduction of 1,1-diaryl-2,2-dicyanoethylenes with SmI(2) in THF was studied in the presence of four proton donors: H(2)O, MeOH, i-PrOH and trifluoroethanol (TFE). The kinetic order for the first two is nearly unity at low proton donor concentrations and approaches four at high concentrations, whereas, for i-PrOH and TFE, the log-log plot is linear with a slope smaller than one. Detailed analysis shows that a curved log-log plot such as for H(2)O and MeOH is indicative of a major contribution by protonation within an ion pair of the radical anion and Sm(+3) complexed to a variable number of proton donor molecules, whereas a linear plot is a result of protonation from the bulk solution.     

Heterogeneous nucleation in the crystallization of polyolefins. III. Theory and mechansim

The nucleating activity of a large number of nucleating agents for the crystallization of polyolefins is not based on a high surface free energy of the substrate or on epitaxy, or on a minimum particle size just compatible with the thickness of the crystallizing polymer lamellae, or on any persistance of crystalline material in holes of the nucleating particles. The relevant mechanism is thought to be nucleation at steps of limited length in the surface of the nucleating particles. The nucleation is enhanced by a high degree of accommodation of polymer molecules at the surface of the substrate through the presence of ditches in the latter which cause prealignment of polymer chains and thus facilitate crystallization.     

Thermodynamics of epitaxial calcite nucleation on self-assembled monolayers

Classical heterogeneous nucleation theory is used to describe the epitaxial nucleation of calcite on self-assembled monolayers (SAMs). Both spherical and faceted clusters are considered. The use of faceted clusters reveals a useful relation between the shape of very small crystals and the ratio of the heterogeneous and homogeneous nucleation barriers. The experimental approach of this paper concerns the measurement of the threshold driving forces for both homogeneous and heterogeneous nucleation of calcite. This is accomplished by preparing solutions with well-defined driving forces and by measuring the resulting types of nucleation that are observed after a fixed experimental time. The results of the experiments and the theoretical shape analysis are compared, and it is shown that in the experiments no homogeneous nucleation of calcite occurs for driving forces up to at least Deltamu/k(B)T approximately equal to 6.0. A calculation of the critical cluster size for heterogeneous nucleation results in a range of 2-28 growth units and faceted critical clusters from 3-28 growth units, depending on the value of the surface free energy of calcite. These sizes are 50-100 times smaller than the crystalline domain sizes of SAMs and therefore small enough to explain the promoting effect of the substrate.     

Polarization-corrected molecular electrostatic potential for the cation binding problem

The molecular electrostatic potential (MESP) and polarization-corrected MESP (PMESP) minima for some small molecules are calculated on the surface generated by rolling cations (Li⁺ and Na⁺) on their van der Waals surfaces. The cation binding energies of these molecules are obtained with HF/6-31G^** level ab initio calculations. A noteworthy outcome of the present study is that the plot of these binding energies and the corresponding PMESP surface minimum values turns out to be remarkably linear with a slope close to unity. The PMESP is thus found to work as a powerful tool for unearthing the patterns of cation binding sites and energetics for molecular systems.     

Effects of charge and size on condensation of supersaturated water vapor on nanoparticles of SiO2

The effects of size and charge on the condensation of a supersaturated water vapor on monodisperse nanoparticles of SiO(2) were investigated in a flow cloud chamber. The dependences of the critical supersaturation S(cr) on particle size at diameters of 10, 12, and 15 nm as well as on charge and charge polarity are determined experimentally. A novel electrospray aerosol generator was developed to generate a high concentration of SiO(2) nanoparticles of less than 10 nm by electrospraying silicon tetraethoxide (STE) ethanol solution followed by the thermal decomposition of STE. The effects of liquid flow rate, liquid concentration, flow rate of carrier gas, and liquid conductivity on the particle size distribution and concentration were examined. For charged particles, the nucleation occurs at a critical supersaturation S(cr) lower than that on neutral particles, and the charge effect fades away as particle size increases. The charge effect is stronger than the theoretical predictions. In addition, a sign preference is detected, i.e., water vapor condenses more readily on negatively charged particle, a trend consistent with those observed on ions. However, both effects of charge and charge polarity on S(cr) are stronger than that predicted by Volmer's theory for ion-induced nucleation.