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.

Investigations of planetary boundary layer processes and particle formation in the atmosphere of planet Mars

The planet Mars is the Earth's neighbour in the Solar System. Planetary research stems from a fundamental need to explore our surroundings, typical for mankind. Manned missions to Mars are already being planned, and understanding the environment to which the astronauts would be exposed is of utmost importance for a successful mission. Information of the Martian environment given by models is already now used in designing the landers and orbiters sent to the red planet. In particular, studies of the Martian atmosphere are crucial for instrument design, entry, descent and landing system design, landing site selection, and aerobraking calculations. Research of planetary atmospheres can also contribute to atmospheric studies of the Earth via model testing and development of parameterizations: even after decades of modeling the Earth's atmosphere, we are still far from perfect weather predictions. On a global level, Mars has also been experiencing climate change. The aerosol effect is one of the largest unknowns in the present terrestrial climate change studies, and the role of aerosol particles in any climate is fundamental: studies of climate variations on another planet can help us better understand our own global change. In this thesis I have used an atmospheric column model for Mars to study the behaviour of the lowest layer of the atmosphere, the planetary boundary layer (PBL), and I have developed nucleation (particle formation) models for Martian conditions. The models were also coupled to study, for example, fog formation in the PBL. The PBL is perhaps the most significant part of the atmosphere for landers and humans, since we live in it and experience its state, for example, as gusty winds, nightfrost, and fogs. However, PBL modelling in weather prediction models is still a difficult task. Mars hosts a variety of cloud types, mainly composed of water ice particles, but also CO2 ice clouds form in the very cold polar night and at high altitudes elsewhere. Nucleation is the first step in particle formation, and always includes a phase transition. Cloud crystals on Mars form from vapour to ice on ubiquitous, suspended dust particles. Clouds on Mars have a small radiative effect in the present climate, but it may have been more important in the past. This thesis represents an attempt to model the Martian atmosphere at the smallest scales with high resolution. The models used and developed during the course of the research are useful tools for developing and testing parameterizations for larger-scale models all the way up to global climate models, since the small-scale models can describe processes that in the large-scale models are reduced to subgrid (not explicitly resolved) scale.

On measurements of formation, growth, hygroscopicity, and volatility of atmospheric ultrafine particles

Atmospheric aerosol particles affect the global climate as well as human health. In this thesis, formation of nanometer sized atmospheric aerosol particles and their subsequent growth was observed to occur all around the world. Typical formation rate of 3 nm particles at varied from 0.01 to 10 cm-3s-1. One order of magnitude higher formation rates were detected in urban environment. Highest formation rates up to 105 cm-3s-1 were detected in coastal areas and in industrial pollution plumes. Subsequent growth rates varied from 0.01 to 20 nm h-1. Smallest growth rates were observed in polar areas and the largest in the polluted urban environment. This was probably due to competition between growth by condensation and loss by coagulation. Observed growth rates were used in the calculation of a proxy condensable vapour concentration and its source rate in vastly different environments from pristine Antarctica to polluted India. Estimated concentrations varied only 2 orders of magnitude, but the source rates for the vapours varied up to 4 orders of magnitude. Highest source rates were in New Delhi and lowest were in the Antarctica. Indirect methods were applied to study the growth of freshly formed particles in the atmosphere. Also a newly developed Water Condensation Particle Counter, TSI 3785, was found to be a potential candidate to detect water solubility and thus indirectly composition of atmospheric ultra-fine particles. Based on indirect methods, the relative roles of sulphuric acid, non-volatile material and coagulation were investigated in rural Melpitz, Germany. Condensation of non-volatile material explained 20-40% and sulphuric acid the most of the remaining growth up to a point, when nucleation mode reached 10 to 20 nm in diameter. Coagulation contributed typically less than 5%. Furthermore, hygroscopicity measurements were applied to detect the contribution of water soluble and insoluble components in Athens. During more polluted days, the water soluble components contributed more to the growth. During less anthropogenic influence, non-soluble compounds explained a larger fraction of the growth. In addition, long range transport to a measurement station in Finland in a relatively polluted air mass was found to affect the hygroscopicity of the particles. This aging could have implications to cloud formation far away from the pollution sources.

Observations on the first steps of atmospheric particle formation and growth

Atmospheric aerosol particles have a significant impact on air quality, human health and global climate. The climatic effects of secondary aerosol are currently among the largest uncertainties limiting the scientific understanding of future and past climate changes. To better estimate the climatic importance of secondary aerosol particles, detailed information on atmospheric particle formation mechanisms and the vapours forming the aerosol is required. In this thesis we studied these issues by applying novel instrumentation in a boreal forest to obtain direct information on the very first steps of atmospheric nucleation and particle growth. Additionally, we used detailed laboratory experiments and process modelling to determine condensational growth properties, such as saturation vapour pressures, of dicarboxylic acids, which are organic acids often found in atmospheric samples. Based on our studies, we came to four main conclusions: 1) In the boreal forest region, both sulphurous compounds and organics are needed for secondary particle formation, the previous contributing mainly to particle formation and latter to growth; 2) A persistent pool of molecular clusters, both neutral and charged, is present and participates in atmospheric nucleation processes in boreal forests; 3) Neutral particle formation seems to dominate over ion-mediated mechanisms, at least in the boreal forest boundary layer; 4) The subcooled liquid phase saturation vapour pressures of C3-C9 dicarboxylic acids are of the order of 1e-5 1e-3 Pa at atmospheric temperatures, indicating that a mixed pre-existing particulate phase is required for their condensation in atmospheric conditions. The work presented in this thesis gives tools to better quantify the aerosol source provided by secondary aerosol formation. The results are particularly useful when estimating, for instance, anthropogenic versus biogenic influences and the fractions of secondary aerosol formation explained by neutral or ion-mediated nucleation mechanisms, at least in environments where the average particle formation rates are of the order of some tens of particles per cubic centimeter or lower. However, as the factors driving secondary particle formation are likely to vary depending on the environment, measurements on atmospheric nucleation and particle growth are needed from around the world to be able to better describe the secondary particle formation, and assess its climatic effects on a global scale.

Analysis of Atmospheric Particle Formation Events

Atmospheric aerosol particle formation events can be a significant source for tropospheric aerosols and thus influence the radiative properties and cloud cover of the atmosphere. This thesis investigates the analysis of aerosol size distribution data containing particle formation events, describes the methodology of the analysis and presents time series data measured inside the Boreal forest. This thesis presents a methodology to identify regional-scale particle formation, and to derive the basic characteristics such as growth and formation rates. The methodology can also be used to estimate concentration and source rates of the vapour causing particle growth. Particle formation was found to occur frequently in the boreal forest area over areas covering up to hundreds of kilometers. Particle formation rates of boreal events were found to be of the order of 0.01-5 cm^-3 s^-1, while the nucleation rates of 1 nm particles can be a few orders of magnitude higher. The growth rates of over 3 nm sized particles were of the order of a few nanometers per hour. The vapor concentration needed to sustain such growth is of the order of 10^7--10^8 cm^-3, approximately one order of magnitude higher than sulphuric acid concentrations found in the atmosphere. Therefore, one has to assume that other vapours, such as organics, have a key role in growing newborn particles to sizes where they can become climatically active. Formation event occurrence shows a clear annual variation with peaks in summer and autumns. This variation is similar to the variation exhibited the obtained formation rates of particles. The growth rate, on the other hand, reaches its highest values during summer. This difference in the annual behavior, and the fact that no coupling between the growth and formation process could be identified, suggest that these processes might be different ones, and that both are needed for a particle formation burst to be observed.

Dynamical processes in the stratosphere and upper troposphere and their influence on the distribution of trace gases in the polar atmosphere

Transport plays an important role in the distribution of long-lived gases such as ozone and water vapour in the atmosphere. Understanding of observed variability in these gases as well as prediction of the future changes depends therefore on our knowledge of the relevant atmospheric dynamics. This dissertation studies certain dynamical processes in the stratosphere and upper troposphere which influence the distribution of ozone and water vapour in the atmosphere. The planetary waves that originate in the troposphere drive the stratospheric circulation. They influence both the meridional transport of substances as well as parameters of the polar vortices. In turn, temperatures inside the polar vortices influence abundance of the Polar Stratospheric Clouds (PSC) and therefore the chemical ozone destruction. Wave forcing of the stratospheric circulation is not uniform during winter. The November-December averaged stratospheric eddy heat flux shows a significant anticorrelation with the January-February averaged eddy heat flux in the midlatitude stratosphere and troposphere. These intraseasonal variations are attributable to the internal stratospheric vacillations. In the period 1979-2002, the wave forcing exhibited a negative trend which was confined to the second half of winter only. In the period 1958-2002, area, strength and longevity of the Arctic polar vortices do not exhibit significant long-term changes while the area with temperatures lower than the threshold temperature for PSC formation shows statistically significant increase. However, the Arctic vortex parameters show significant decadal changes which are mirrored in the ozone variability. Monthly ozone tendencies in the Northern Hemisphere show significant correlations (|r|=0.7) with proxies of the stratospheric circulation. In the Antarctic, the springtime vortex in the lower stratosphere shows statistically significant trends in temperature, longevity and strength (but not in area) in the period 1979-2001. Analysis of the ozone and water vapour vertical distributions in the Arctic UTLS shows that layering below and above the tropopause is often associated with poleward Rossby wave-breaking. These observations together with calculations of cross-tropopause fluxes emphasize the importance of poleward Rossby wave breaking for the stratosphere-troposphere exchange in the Arctic.

Size-dependent magnetic properties of iron carbide nanoparticles embedded in a carbon matrix

Here we report on the magnetic properties of iron carbide nanoparticles embedded in a carbon matrix. Granular distributions of nanoparticles in an inert matrix, of potential use in various applications, were prepared by pyrolysis of organic precursors using the thermally assisted chemical vapour deposition method. By varying the precursor concentration and preparation temperature, compositions with varying iron concentration and nanoparticle sizes were made. Powder x-ray diffraction, transmission electron microscopy and Mossbauer spectroscopy studies revealed the nanocrystalline iron carbide (Fe3C) presence in the partially graphitized matrix. The dependence of the magnetic properties on the particle size and temperature (10 K < T < 300 K) were studied using superconducting quantum interference device magnetometry. Based on the affect of surrounding carbon spins, the observed magnetic behaviour of the nanoparticle compositions, such as the temperature dependence of magnetization and coercivity, can be explained.

Raman spectrum of ethyl chloroacetate

The Raman spectrum of ethyl chloroacetate has been studied at 13° C., 28° C. and 78° C. The carbonyl frequency was found to be split up into two due to the presence of rotational isomers. The higher frequency line due to thecis isomer was found to decrease in intensity with temperature. It appears that the gauche isomer will predominate in the vapour state. Altogether thirty-eight Raman lines have been recorded. Reasonable assignments for the observed Raman lines were made in comparison with ethyl acetate spectrum.

Vibrational analysis of the (B2P{cyrillic}→X2P{cyrillic}) system of NS

A few red degraded bands attributable to NS have been reported earlier by Fowler and Barker, Dressler and Barrow et al, and they occur in the same region (2300 to 2700 Å) as the bands of the known systems (C2∑+-X2P{cyrillic}) and (A2Δ-X2P{cyrillic}). Measurements made on the heads of some of these weak bands led Barrow et al. to believe that these bands may form a system analogous to the β-system of NO and be due to a2P{cyrillic}-2P{cyrillic} transition. The spectrum of NS has now been studied in a little more detail by means of an uncondensed discharge through dry nitrogen and sulphur vapour in the presence of argon and thirty three bands belonging to this system have been recorded in the region 2280 to 2760 Å. It has been found possible to represent the band heads by means of the equation {Mathematical expression}. Taking the lower state doublet interval as 223 cm-1, it is shown that the separation in the upper state is 94 cm-1. The ratio of the force constants in the upper and the ground states is found to be 0·39 and is nearly the same as that in the β-system of NO (0·30). The present vibrational analysis therefore supports the view that these new red degraded bands of NS arise from a (B2P{cyrillic}→X2P{cyrillic}) transition and the observed intensity distribution in the form of a wide parabola is also in qualitative agreement with what is expected from the moderately large Δ re (∼0·12Å) value.

Symmetry and shape issues in nanostructure growth

Symmetry plays a key role in dictating the equilibrium morphology of crystals. However, several growth morphologies that deviate from the point group symmetry are routinely observed under several different growth conditions. In this article, we present a summary of symmetry-breaking mechanisms that are operative for crystals grown from the vapour phase as well as those formed as a result of wet chemical synthesis. This understanding is crucial for rationalizing the variety of morphologies observed during nanocrystal synthesis and also providesa rational framework for the synthesis of anisotropic nanostructures.

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.

On photooxidation of Ge in obliquely deposited Ge and Ge25X75(X = S, Se, Te) thin films

It has been established by photoemission studies that Ge in obliquely deposited pure Ge and Ge-chalcogenide thin films undergoes predominant photooxidation when irradiated with band gap photons. The role of Ge appears to be that of providing a highly porous low density microstructure and photooxidation seems to be a direct consequence of such large scale porosity in these films. The formation of low vapour pressure oxide fractions of Ge and Te and volatile high vapour pressure oxide fractions of S and Se is responsible for anomalous photoinduced transformations in these films.

Saturated liquid densities of cryogenic liquids and refrigerants

The saturated liquid density, varrholr, data along the liquid vapour coexistence curve published in the literature for several cryogenic liquids, hydrocarbons and halocarbon refrigerants are fitted to a generalized equation of the following form varrholr = 1 + A(1 − Tr + B(1 − Tr)β The values of β, the index in phase density differences power law, have been obtained by means of two approaches namely statistical treatment of saturated fluid phase density difference data and the existence of a maximum in T(varrho1 − varrhov) along the saturation curve. Values of the constants A and B are determined utilizing the fact that Tvarrho1 has a maximum at a characteristic temperature T. Values of A, B and β are tabulated for Ne, Ar, Kr, Xe, N2, O2, methane, ethane, propane, iso-butane, n-butane, propylene, ethylene, CO2, water, ammonia, refrigerants-11, 12, 12B1, 13, 13B1, 14, 21, 22, 23, 32, 40, 113, 114, 115, 142b, 152a, 216, 245 and azeotropes R-500, 502, 503, 504. The average error of prediction is less than 2%.

Growth and characterization of germanium nanowires by electron beam evaporation

Germanium nanowires were grown on Au coated Si substrates at 380 degrees C in a high vacuum (5 x 10(-5) Torr) by e-beam evaporation of Germanium (Ge). The morphology observation by a field emission scanning electron microscope (FESEM) shows that the grown nanowires are randomly oriented with an average length and diameter of 600 nm and 120 nm respectively for a deposition time of 60 min. The nanowire growth ratewas measured to be similar to 10 nm/min. Transmission electron microscope (TEM) studies revealed that the Ge nanowires were single crystalline in nature and further energy dispersive X-ray analysis(EDAX) has shown that the tip of the grown nanowires was capped with Au nanoparticles, this shows that the growth of the Ge nanowires occurs by the vapour liquid solid (VLS) mechanism. HRTEM studies on the grown Ge nanowire show that they are single crystalline in nature and the growth direction was identified to be along [110]. (C) 2010 Elsevier B.V. All rights reserved.

Adsorption of oxygen and reactivity with HCl on a barium-doped lead surface

X.p.s. studies on the adsorption of oxygen on a barium-covered Pb surface have shown the presence of two distinct types of oxygen species: oxidic, O2–, and the peroxo-like O2–2(ads), and the surface has been identified as a composite of PbO and BaPbO3. On a barium pre-covered surface, the sticking probability of oxygen on Pb is increased. The O2–(ads) species preferentially reacts with HCl forming PbCl2(ads)via proton abstraction, whereas O2–2(ads) is not reactive with HCl vapour. On the Pb surface, the PbCl2 overlayer reacts with excess HCl, forming a volatile compound believed to be Pb(ClHCl)2, while in the presence of coadsorbed barium, the stability of PbCl2 is increased and the activation energy for the reaction: PbCl2(ads)+ 2HCl(g) Pb(ClHCl)2(g) is increased. Stronger intermetallic interaction is suggested to be the reason for higher PbCl2 stability.