Theoretical and computational approaches on heterogeneous nucleation

Nucleation is the first step of a first order phase transition. A new phase is always sprung up in nucleation phenomena. The two main categories of nucleation are homogeneous nucleation, where the new phase is formed in a uniform substance, and heterogeneous nucleation, when nucleation occurs on a pre-existing surface. In this thesis the main attention is paid on heterogeneous nucleation. This thesis wields the nucleation phenomena from two theoretical perspectives: the classical nucleation theory and the statistical mechanical approach. The formulation of the classical nucleation theory relies on equilibrium thermodynamics and use of macroscopically determined quantities to describe the properties of small nuclei, sometimes consisting of just a few molecules. The statistical mechanical approach is based on interactions between single molecules, and does not bear the same assumptions as the classical theory. This work gathers up the present theoretical knowledge of heterogeneous nucleation and utilizes it in computational model studies. A new exact molecular approach on heterogeneous nucleation was introduced and tested by Monte Carlo simulations. The results obtained from the molecular simulations were interpreted by means of the concepts of the classical nucleation theory. Numerical calculations were carried out for a variety of substances nucleating on different substances. The classical theory of heterogeneous nucleation was employed in calculations of one-component nucleation of water on newsprint paper, Teflon and cellulose film, and binary nucleation of water-n-propanol and water-sulphuric acid mixtures on silver nanoparticles. The results were compared with experimental results. The molecular simulation studies involved homogeneous nucleation of argon and heterogeneous nucleation of argon on a planar platinum surface. It was found out that the use of a microscopical contact angle as a fitting parameter in calculations based on the classical theory of heterogeneous nucleation leads to a fair agreement between the theoretical predictions and experimental results. In the presented cases the microscopical angle was found to be always smaller than the contact angle obtained from macroscopical measurements. Furthermore, molecular Monte Carlo simulations revealed that the concept of the geometrical contact parameter in heterogeneous nucleation calculations can work surprisingly well even for very small clusters.

Monte Carlo simulations of molecular clusters in nucleation

A better understanding of the limiting step in a first order phase transition, the nucleation process, is of major importance to a variety of scientific fields ranging from atmospheric sciences to nanotechnology and even to cosmology. This is due to the fact that in most phase transitions the new phase is separated from the mother phase by a free energy barrier. This barrier is crossed in a process called nucleation. Nowadays it is considered that a significant fraction of all atmospheric particles is produced by vapor-to liquid nucleation. In atmospheric sciences, as well as in other scientific fields, the theoretical treatment of nucleation is mostly based on a theory known as the Classical Nucleation Theory. However, the Classical Nucleation Theory is known to have only a limited success in predicting the rate at which vapor-to-liquid nucleation takes place at given conditions. This thesis studies the unary homogeneous vapor-to-liquid nucleation from a statistical mechanics viewpoint. We apply Monte Carlo simulations of molecular clusters to calculate the free energy barrier separating the vapor and liquid phases and compare our results against the laboratory measurements and Classical Nucleation Theory predictions. According to our results, the work of adding a monomer to a cluster in equilibrium vapour is accurately described by the liquid drop model applied by the Classical Nucleation Theory, once the clusters are larger than some threshold size. The threshold cluster sizes contain only a few or some tens of molecules depending on the interaction potential and temperature. However, the error made in modeling the smallest of clusters as liquid drops results in an erroneous absolute value for the cluster work of formation throughout the size range, as predicted by the McGraw-Laaksonen scaling law. By calculating correction factors to Classical Nucleation Theory predictions for the nucleation barriers of argon and water, we show that the corrected predictions produce nucleation rates that are in good comparison with experiments. For the smallest clusters, the deviation between the simulation results and the liquid drop values are accurately modelled by the low order virial coefficients at modest temperatures and vapour densities, or in other words, in the validity range of the non-interacting cluster theory by Frenkel, Band and Bilj. Our results do not indicate a need for a size dependent replacement free energy correction. The results also indicate that Classical Nucleation Theory predicts the size of the critical cluster correctly. We also presents a new method for the calculation of the equilibrium vapour density, surface tension size dependence and planar surface tension directly from cluster simulations. We also show how the size dependence of the cluster surface tension in equimolar surface is a function of virial coefficients, a result confirmed by our cluster simulations.

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.

Generalised vortex motion of compressible fluid with and without magnetic field and suction

The laminar boundary layer over a stationary infinite disk induced by a rotating compressible fluid is considered. The free stream velocity has been taken as tangential and varies as a power of radius, i.e. v∞ ˜ r−n. The effect of the axial magnetic field and suction is also included in the analysis. An implicit finite difference scheme is employed to the governing similarity equations for numerical computations. Solutions are studied for various values of disk to fluid temperature ratio and for values of n between 1 and −1. In the absence of the magnetic field and suction, velocity profiles exhibit oscillations. It has been observed that for a hot disk in the presence of a magnetic field the boundary layer solutions decay algebraically instead of decaying exponentially. In the absence of the magnetic field and suction, the solution of the similarity equations exists only for a certain range of n.

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.

Effect of natural ventilation on the boundary layer separation and near-wake vortex shedding characteristics of a sphere

Experiments were conducted in water and wind tunnels on spheres in the Reynolds number range 6 x 10(3) to 6.5 x 10(5) to study the effect of natural ventilation on the boundary layer separation and near-wake Vortex shedding characteristics. In the subcritical range of Re (<2 x 10(5)), ventilation caused a marginal downstream shift in the location of laminar boundary layer separation; there was only a small change in the vortex shedding frequency. In the supercritical range (Re > 4 x 10(5)), ventilation caused a downstream shift in the mean locations of boundary layer separation and reattachment; these lines showed significant axisymmetry in the presence of venting. No distinct vortex shedding frequency was found. Instead, a dramatic reduction occurred in the wake unsteadiness at all frequencies. The reduction of wake unsteadiness is consistent with the reduction in total drag already reported. Based on the present results and those reported earlier, the effects of natural ventilation on the flow past a sphere can be categorized in two broad regimes, viz., weak and strong interaction regimes. In the weak interaction regime (subcritical Re), the broad features of the basic sphere are largely unaltered despite the large addition of mass in the near wake. Strong interaction is promoted by the closer proximity of the inner and outer shear layers at supercritical Re. This results in a modified and steady near-wake flow, characterized by reduced unsteadiness and small drag.

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.

Some observations of tip-vortex cavitation

Cavitation has been observed in the trailing vortex system of an elliptic planform hydrofoil. A complex dependence on Reynolds number and gas content is noted at inception. Some of the observations can be related to tension effects associated with the lack of sufficiently large-sized nuclei. Inception measurements are compared with estimates of pressure in the vortex obtained from LDV measurements of velocity within the vortex. It is concluded that a complete correlation is not possible without knowledge of the fluctuating levels of pressure in tip-vortex flows. When cavitation is fully developed, the observed tip-vortex trajectory shows a surprising lack of dependence on any of the physical parameters varied, such as angle of attack, Reynolds number, cavitation number, and dissolved gas content.

Dependence of the large-scale vortex instability on latitude, stratification and domain size

In an earlier study, we reported on the excitation of large-scale vortices in Cartesian hydrodynamical convection models subject to rapid enough rotation. In that study, the conditions for the onset of the instability were investigated in terms of the Reynolds (Re) and Coriolis (Co) numbers in models located at the stellar North pole. In this study, we extend our investigation to varying domain sizes, increasing stratification, and place the box at different latitudes. The effect of the increasing box size is to increase the sizes of the generated structures, so that the principal vortex always fills roughly half of the computational domain. The instability becomes stronger in the sense that the temperature anomaly and change in the radial velocity are observed to be enhanced. The model with the smallest box size is found to be stable against the instability, suggesting that a sufficient scale separation between the convective eddies and the scale of the domain is required for the instability to work. The instability can be seen upto the colatitude of 30 degrees, above which value the flow becomes dominated by other types of mean flows. The instability can also be seen in a model with larger stratification. Unlike the weakly stratified cases, the temperature anomaly caused by the vortex structures is seen to depend on depth.