Molecular Dynamics Simulations of Homogeneous Nucleation

Nucleation is the first step in a phase transition where small nuclei of the new phase start appearing in the metastable old phase, such as the appearance of small liquid clusters in a supersaturated vapor. Nucleation is important in various industrial and natural processes, including atmospheric new particle formation: between 20 % to 80 % of atmospheric particle concentration is due to nucleation. These atmospheric aerosol particles have a significant effect both on climate and human health. Different simulation methods are often applied when studying things that are difficult or even impossible to measure, or when trying to distinguish between the merits of various theoretical approaches. Such simulation methods include, among others, molecular dynamics and Monte Carlo simulations. In this work molecular dynamics simulations of the homogeneous nucleation of Lennard-Jones argon have been performed. Homogeneous means that the nucleation does not occur on a pre-existing surface. The simulations include runs where the starting configuration is a supersaturated vapor and the nucleation event is observed during the simulation (direct simulations), as well as simulations of a cluster in equilibrium with a surrounding vapor (indirect simulations). The latter type are a necessity when the conditions prevent the occurrence of a nucleation event in a reasonable timeframe in the direct simulations. The effect of various temperature control schemes on the nucleation rate (the rate of appearance of clusters that are equally able to grow to macroscopic sizes and to evaporate) was studied and found to be relatively small. The method to extract the nucleation rate was also found to be of minor importance. The cluster sizes from direct and indirect simulations were used in conjunction with the nucleation theorem to calculate formation free energies for the clusters in the indirect simulations. The results agreed with density functional theory, but were higher than values from Monte Carlo simulations. The formation energies were also used to calculate surface tension for the clusters. The sizes of the clusters in the direct and indirect simulations were compared, showing that the direct simulation clusters have more atoms between the liquid-like core of the cluster and the surrounding vapor. Finally, the performance of various nucleation theories in predicting simulated nucleation rates was investigated, and the results among other things highlighted once again the inadequacy of the classical nucleation theory that is commonly employed in nucleation studies.

Peptide hairpin nucleation with the obligatory Type I' beta-turn Aib-(D)Pro segment

The alpha-aminoisobutyric acid-D-proline (Aib-(D)Pro) dipeptide is an obligatory Type I' beta-turn forming segment that nucleates hairpin formation.

How the development of condensation particle counters is reforming our view on atmospheric nucleation

Aerosol particles deteriorate air quality, atmospheric visibility and our health. They affect the Earth s climate by absorbing and scattering sunlight, forming clouds, and also via several feed-back mechanisms. The net effect on the radiative balance is negative, i.e. cooling, which means that particles counteract the effect of greenhouse gases. However, particles are one of the poorly known pieces in the climate puzzle. Some of the airborne particles are natural, some anthropogenic; some enter the atmosphere in particle form, while others form by gas-to-particle conversion. Unless the sources and dynamical processes shaping the particle population are quantified, they cannot be incorporated into climate models. The molecular level understanding of new particle formation is still inadequate, mainly due to the lack of suitable measurement techniques to detect the smallest particles and their precursors. This thesis has contributed to our ability to measure newly formed particles. Three new condensation particle counter applications for measuring the concentration of nano-particles were developed. The suitability of the methods for detecting both charged and electrically neutral particles and molecular clusters as small as 1 nm in diameter was thoroughly tested both in laboratory and field conditions. It was shown that condensation particle counting has reached the size scale of individual molecules, and besides measuring the concentration they can be used for getting size information. In addition to atmospheric research, the particle counters could have various applications in other fields, especially in nanotechnology. Using the new instruments, the first continuous time series of neutral sub-3 nm particle concentrations were measured at two field sites, which represent two different kinds of environments: the boreal forest and the Atlantic coastline, both of which are known to be hot-spots for new particle formation. The contribution of ions to the total concentrations in this size range was estimated, and it could be concluded that the fraction of ions was usually minor, especially in boreal forest conditions. Since the ionization rate is connected to the amount of cosmic rays entering the atmosphere, the relative contribution of neutral to charged nucleation mechanisms extends beyond academic interest, and links the research directly to current climate debate.

On the charging state of atmospheric aerosols and ion-induced nucleation

Floating in the air that surrounds us is a number of small particles, invisible to the human eye. The mixture of air and particles, liquid or solid, is called an aerosol. Aerosols have significant effects on air quality, visibility and health, and on the Earth's climate. Their effect on the Earth's climate is the least understood of climatically relevant effects. They can scatter the incoming radiation from the Sun, or they can act as seeds onto which cloud droplets are formed. Aerosol particles are created directly, by human activity or natural reasons such as breaking ocean waves or sandstorms. They can also be created indirectly as vapors or very small particles are emitted into the atmosphere and they combine to form small particles that later grow to reach climatically or health relevant sizes. The mechanisms through which those particles are formed is still under scientific discussion, even though this knowledge is crucial to make air quality or climate predictions, or to understand how aerosols will influence and will be influenced by the climate's feedback loops. One of the proposed mechanisms responsible for new particle formation is ion-induced nucleation. This mechanism is based on the idea that newly formed particles were ultimately formed around an electric charge. The amount of available charges in the atmosphere varies depending on radon concentrations in the soil and in the air, as well as incoming ionizing radiation from outer space. In this thesis, ion-induced nucleation is investigated through long-term measurements in two different environments: in the background site of Hyytiälä and in the urban site that is Helsinki. The main conclusion of this thesis is that ion-induced nucleation generally plays a minor role in new particle formation. The fraction of particles formed varies from day to day and from place to place. The relative importance of ion-induced nucleation, i.e. the fraction of particles formed through ion-induced nucleation, is bigger in cleaner areas where the absolute number of particles formed is smaller. Moreover, ion-induced nucleation contributes to a bigger fraction of particles on warmer days, when the sulfuric acid and water vapor saturation ratios are lower. This analysis will help to understand the feedbacks associated with climate change.

Adsorption-nucleation/growth based model for electrochemical phase formation and Avrami ansatz for a multicomponent competitive growth process

The existing models describing electrochemical phase formation involving both adsorption and a nucleation/growth process are modified. The limiting cases leading to the existing models are discussed. The characteristic features of the potentiostatic transients are presented. A generalization of the Avrami ansatz is given for two or more competitive irreversibly growing phases.

Heterogeneous nucleation of solidification of cadmium particles embedded in an aluminium matrix

A hypomonotectic alloy of Al-4.5wt%Cd has been manufactured by melt spinning and the resulting microstructure examined by transmission electron microscopy. As-melt spun hypomonotectic Al-4.5wt%Cd consists of a homogeneous distribution of faceted 5 to 120 nm diameter cadmium particles embedded in a matrix of aluminium, formed during the monotectic solidification reaction. The cadmium particles exhibit an orientation relationship with the aluminium matrix of {111}Al//{0001}Cd and lang110rangAlAl//lang11¯20> Cd, with four cadmium particle variants depending upon which of the four {111}Al planes is parallel to {0001}Cd. The cadmium particles exibit a distorted cuboctahedral shape, bounded by six curved {100}Al//{20¯23}Cd facets, six curved {111}Al/{40¯43}Cd facets and two flat {111}Al//{0001}Cd facets. The as-melt spun cadmium particle shape is metastable and the cadmium particles equilibrate during heat treatment below the cadmium melting point, becoming elongated to increase the surface area and decrease the separation of the {111}Al//{0001}Cd facets. The equilibrium cadmium particle shape and, therefore, the anisotropy of solid aluminium-solid cadmium and solid aluminium -liquid cadmium surface energies have been monitored by in situ heating in the transmission electron microscope over the temperature range between room temperature and 420 °C. The anisotropy of solid aluminium-solid cadmium surface energy is constant between room temperature and the cadmium melting point, with the {100}Al//{20¯23}Cd surface energy on average 40% greater than the {111}Al//{0001}Cd surface energy, and 10% greater than the {111}Al//{40¯43Cd surface energy. When the cadmium particles melt at temperatures above 321 °C, the {100}Al//{20¯23}Cd facets disappear and the {111}Al//{40¯43}Cd and {111}A1//{0001}Cd surface energies become equal. The {111}Al facets do not disappear when the cadmium particles melt, and the anisotropy of solid aluminium-liquid cadmium surface energy decreases gradually with increasing temperature above the cadmium melting point. The kinetics of cadmium solidification have been examined by heating and cooling experiments in a differential scanning calorimeter over a range of heating and cooling rates. Cadmium particle solidification is nucleated catalytically by the surrounding aluminium matrix on the {111}Al faceted surfaces, with an undercooling of 56 K and a contact angle of 42 °. The nucleation kinetics of cadmium particle solidification are in good agreement with the hemispherical cap model of heterogeneous nucleation.

Nucleation and growth of Al2O3/metal composites by oxidation of aluminum alloys

The nucleation and growth mechanisms during high temperature oxidation of liquid Al-3% Mg and Al-3% Mg-3% Si alloys were studied with the aim of enhancing our understanding of a new composite fabrication process. The typical oxidation sequence consists of an initial event of rapid but brief oxidation, followed by an incubation period of limited oxide growth after which bulk Al2O3/Al composite forms. A duplex oxide layer, MgO (upper) and MgAl2O4 (lower), forms on the alloy surface during initial oxidation and incubation. The spinel layer remains next to the liquid alloy during bulk oxide growth and is the eventual repository for most of the magnesium in the original alloy. Metal microchannels developed during incubation continuously supply alloy through the composite to the reaction interface. During the growth process, a layered structure exists at the upper extremity of the composite, consisting of MgO at the top surface, MgAl2O4 (probably discontinuous), Al alloy, and finally the bulk Al2O3 composite containing microchannels of the alloy. The bulk oxide growth mechanism appears to involve continuous formation and dissolution of the Mg-rich oxides at the surface, diffusion of oxygen through the underlying liquid metal, and epitaxial growth of Al2O3 on the existing composite body. The roles of Mg and Si in the composite growth process are discussed.

Bubble formation at closely spaced orifices in aqueous solutions

A filter cloth with 182 holes per 10−4 m2 has been used to generate air bubbles both in pure water and in aqueous solutions of electrolytes and non-electrolytes at various air flow rates. Potassium bromide and ammonium perchlorate were the electrolytes used, while the non-electrolytes were isopropanol, urea and glycerol. Bubble diameters and their size distribution were measured from photographs. The role of solutes in affecting bubble sizes and their distribution compared to that of pure water is discussed in the light of a hypothesis. This hypothesis assumes that if the final bubble diameter is less than the inter-orifice distance, then bubbles do not coalesce; on the other hand, if it is greater, then coalescence occurs when tf greater-or-equal, slantedti+ts, but does not occur when t

Influence of various gases on single bubble sonoluminescence

Influence of various gases on the intensity of single bubble sonoluminescence has been studied. The gases used were air, oxygen, nitrogen, argon and helium. Among these oxygen gave the brightest intensity with nitrogen giving the least.

Observations on transition in plane bubble plumes

Flow visualization studies of plane laminar bubble plumes have been conducted to yield quantitative data on transition height, wavelength and wave velocity of the most unstable disturbance leading to transition. These are believed to be the first results of this kind. Most earlier studies are restricted to turbulent bubble plumes. In the present study, the bubble plumes were generated by electrolysis of water and hence very fine control over bubble size distribution and gas flow rate was possible to enable studies with laminar bubble plumes. Present observations show that (a) the dominant mode of instability in plane bubble plumes is the sinuous mode, (b) transition height and wavelength are related linearly with the proportionality constant being about 4, (c) wave velocity is about 40 % of the mean plume velocity, and (d) normalized transition height data correlate very well with a source Grashof number. Some agreement and some differences in transition characteristics of bubble plumes have been observed compared to those for similar single-phase flows.

Heterogeneous nucleation of pb particles embedded in a Zn matrix

Zinc-10 and 20 wt pct Pb alloys have been rapidly solidified by melt spinning to obtain a very fine scale dispersion of nanometer-sized Pb particles embedded in Zn matrix. The microstructure and crystallography of the Pb particles have been studied using transmission electron microscopy (TEM). Each embedded Pb particle is a single crystal, with a truncated hexagonal biprism shape with the 6/mmm Zn matrix point group symmetry surrounded by and { 0001 á },\text { \text10[`\text1] \text0 },\text and { \text10[`\text1] \text1 }0001 1010 and 1011 facets. The Pb particles solidify with a well-defined orientation relationship with the Zn matrix of ( 0001 )Zn ||(111)Pb\text and\text [ \text11[`\text2] \text0 ]Zn| ||[ 1[`1] 0 ]Pb 0001Zn(111)Pb and 1120Zn110Pb . The melting and solidification behavior of the Pb particle have been studied using differential scanning calorimetry (DSC). The Pb particles solidify with an undercooling of approximately 30 K, by heterogeneous nucleation on the {0001} facets of the surrounding Zn matrix, with an apparent contact angle of 23 deg.

Electrochemical phase formation involving multicomponents: hemispherical overlap models for varied nucleation and growth conditions

The phenomenological theory of hemispherical growth in the context of phase formation with more than one component is presented. The model discusses in a unified manner both instantaneous and progressive nucleation (at the substrate) as well as arbitrary growth rates (e.g. constant and diffusion controlled growth rates). A generalized version of Avrami ansatz (a mean field description) is used to tackle the ''overlap'' aspects arising from the growing multicentres of the many components involved, observing that the nucleation is confined to the substrate plane only. The time evolution of the total extent of macrogrowth as well as those of the individual components are discussed explicitly for the case of two phases. The asymptotic expressions for macrogrowth are derived. Such analysis depicts a saturation limit (i.e. the maximum extent of growth possible) for the slower growing component and its dependence on the kinetic parameters which, in the electrochemical context, can be controlled through potential. The significance of this model in the context of multicomponent alloy deposition and possible future directions for further development are pointed out.

Crossover dynamics at large metastability in gas-liquid nucleation

We have developed an alternate description of dynamics of nucleation in terms of an extended set of order parameters. The order parameters consist of an ordered set of kth largest clusters, ordered such that k = 1 is the largest cluster in the system, k = 2 is the second largest cluster, and so on. We have derived an analytic expression for the free energy for the kth largest cluster, which is in excellent agreement with the simulated results. At large supersaturation, the free energy barrier for the growth of the kth largest cluster disappears and the nucleation becomes barrierless. The major success of this extended theoretical formalism is that it can clearly explain the observed change in mechanism at large metastability P. Bhimalapuram et al., Phys. Rev. Lett. 98, 206104 (2007)] and the associated dynamical crossover. The classical nucleation theory cannot explain this crossover. The crossover from activated to barrierless nucleation is found to occur at a supersaturation where multiple clusters cross the critical size. We attribute the crossover as the onset of the kinetic spinodal. We have derived an expression for the rate of nucleation in the barrierless regime by modeling growth as diffusion on the free energy surface of the largest cluster. The model reproduces the slower increase in the rate of growth as a function of supersaturation, as observed in experiments.