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 studies of nanoparticle impacts

This thesis concerns the dynamics of nanoparticle impacts on solid surfaces. These impacts occur, for instance, in space, where micro- and nanometeoroids hit surfaces of planets, moons, and spacecraft. On Earth, materials are bombarded with nanoparticles in cluster ion beam devices, in order to clean or smooth their surfaces, or to analyse their elemental composition. In both cases, the result depends on the combined effects of countless single impacts. However, the dynamics of single impacts must be understood before the overall effects of nanoparticle radiation can be modelled. In addition to applications, nanoparticle impacts are also important to basic research in the nanoscience field, because the impacts provide an excellent case to test the applicability of atomic-level interaction models to very dynamic conditions. In this thesis, the stopping of nanoparticles in matter is explored using classical molecular dynamics computer simulations. The materials investigated are gold, silicon, and silica. Impacts on silicon through a native oxide layer and formation of complex craters are also simulated. Nanoparticles up to a diameter of 20 nm (315000 atoms) were used as projectiles. The molecular dynamics method and interatomic potentials for silicon and gold are examined in this thesis. It is shown that the displacement cascade expansionmechanism and crater crown formation are very sensitive to the choice of atomic interaction model. However, the best of the current interatomic models can be utilized in nanoparticle impact simulation, if caution is exercised. The stopping of monatomic ions in matter is understood very well nowadays. However, interactions become very complex when several atoms impact on a surface simultaneously and within a short distance, as happens in a nanoparticle impact. A high energy density is deposited in a relatively small volume, which induces ejection of material and formation of a crater. Very high yields of excavated material are observed experimentally. In addition, the yields scale nonlinearly with the cluster size and impact energy at small cluster sizes, whereas in macroscopic hypervelocity impacts, the scaling 2 is linear. The aim of this thesis is to explore the atomistic mechanisms behind the nonlinear scaling at small cluster sizes. It is shown here that the nonlinear scaling of ejected material yield disappears at large impactor sizes because the stopping mechanism of nanoparticles gradually changes to the same mechanism as in macroscopic hypervelocity impacts. The high yields at small impactor size are due to the early escape of energetic atoms from the hot region. In addition, the sputtering yield is shown to depend very much on the spatial initial energy and momentum distributions that the nanoparticle induces in the material in the first phase of the impact. At the later phases, the ejection of material occurs by several mechanisms. The most important mechanism at high energies or at large cluster sizes is atomic cluster ejection from the transient liquid crown that surrounds the crater. The cluster impact dynamics detected in the simulations are in agreement with several recent experimental results. In addition, it is shown that relatively weak impacts can induce modifications on the surface of an amorphous target over a larger area than was previously expected. This is a probable explanation for the formation of the complex crater shapes observed on these surfaces with atomic force microscopy. Clusters that consist of hundreds of thousands of atoms induce long-range modifications in crystalline gold.

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 Caches: A caching structure for dynamic creation of application-specific Heterogeneous cache regions

CMPs enable simultaneous execution of multiple applications on the same platforms that share cache resources. Diversity in the cache access patterns of these simultaneously executing applications can potentially trigger inter-application interference, leading to cache pollution. Whereas a large cache can ameliorate this problem, the issues of larger power consumption with increasing cache size, amplified at sub-100nm technologies, makes this solution prohibitive. In this paper in order to address the issues relating to power-aware performance of caches, we propose a caching structure that addresses the following: 1. Definition of application-specific cache partitions as an aggregation of caching units (molecules). The parameters of each molecule namely size, associativity and line size are chosen so that the power consumed by it and access time are optimal for the given technology. 2. Application-Specific resizing of cache partitions with variable and adaptive associativity per cache line, way size and variable line size. 3. A replacement policy that is transparent to the partition in terms of size, heterogeneity in associativity and line size. Through simulation studies we establish the superiority of molecular cache (caches built as aggregations of molecules) that offers a 29% power advantage over that of an equivalently performing traditional cache.

Myotonic dystrophy type 2 (DM2) : Diagnostic methods and molecular pathology

Myotonic dystrophies type 1 (DM1) and type 2 (DM2) are the most common forms of muscular dystrophy affecting adults. They are autosomal dominant diseases caused by microsatellite tri- or tetranucleotide repeat expansion mutations in transcribed but not translated gene regions. The mutant RNA accumulates in nuclei disturbing the expression of several genes. The more recently identified DM2 disease is less well known, yet more than 300 patients have been confirmed in Finland thus far, and the true number is believed to be much higher. DM1 and DM2 share some features in general clinical presentation and molecular pathology, yet they show distinctive differences, including disease severity and differential muscle and fiber type involvement. However, the molecular differences underlying DM1 and DM2 muscle pathology are not well understood. Although the primary tissue affected is muscle, both DMs show a multisystemic phenotype due to wide expression of the mutation-carrying genes. DM2 is particularly intriguing, as it shows an incredibly wide spectrum of clinical manifestations. For this reason, it constitutes a real diagnostic challenge. The core symptoms in DM2 include proximal muscle weakness, muscle pain, myotonia, cataracts, cardiac conduction defects and endocrinological disturbations; however, none of these is mandatory for the disease. Myalgic pains may be the most disabling symptom for decades, sometimes leading to incapacity for work. In addition, DM2 may cause major socio-economical consequences for the patient, if not diagnosed, due to misunderstanding and false stigmatization. In this thesis work, we have (I) improved DM2 differential diagnostics based on muscle biopsy, and (II) described abnormalities in mRNA and protein expression in DM1 and DM2 patient skeletal muscles, showing partial differences between the two diseases, which may contribute to muscle pathology in these diseases. This is the first description of histopathological differences between DM1 and DM2, which can be used in differential diagnostics. Two novel high-resolution applications of in situ -hybridization have been described, which can be used for direct visualization of the DM2 mutation in muscle biopsy sections, or mutation size determination on extended DNA-fibers. By measuring protein and mRNA expression in the samples, differential changes in expression patterns affecting contractile proteins, other structural proteins and calcium handling proteins in DM2 compared to DM1 were found. The dysregulation at mRNA level was caused by altered transciption and abnormal splicing. The findings reported here indicate that the extent of aberrant splicing is higher in DM2 compared to DM1. In addition, the described abnormalities to some extent correlate to the differences in fiber type involvement in the two disorders.

Droplet Epitaxy of InN Quantum Dots on Si(111) by RF Plasma-Assisted Molecular Beam Epitaxy

InN quantum dots (QDs) were fabricated on Si(111) substrate by droplet epitaxy using an RF plasma-assisted MBE system. Variation of the growth parameters, such as growth temperature and deposition time, allowed us to control the characteristic size and density of the QDs. As the growth temperature was increased from 100 C to 300 degrees C, an enlargement of QD size and a drop in dot density were observed, which was led by the limitation of surface diffusion of adatoms with the limited thermal energy. Atomic force microscopy (AFM) and scanning electron microscopy (SEM) were used to assess the QDs size and density. The chemical bonding configurations of InN QDs were examined by X-ray photo-electron spectroscopy (XPS). Fourier transform infrared (FTIR) spectrum of the deposited InN QDs shows the presence of In-N bond. Temperature-dependent photoluminescence (PL) measurements showed that the emission peak energies of the InN QDs are sensitive to temperature and show a strong peak emission at 0.79 eV.

Ions in water: Role of attractive interactions in size dependent diffusivity maximum

A molecular dynamics study of model ions in water is reported. The van der Waals diameter of both the cations and anions is varied. We have carried out two sets of simulations-with and without dispersion interaction-between the ion and water. Self-diffusivity of the ions exhibits an anomalous maximum as a function of the van der Waals diameter for both these sets. This existence of a maximum in self-diffusivity when there is no dispersion interaction between the ion and the water is attributed to the attractive term from electrostatic interactions. Detailed analysis of this effect shows that the solvent shell is more strongly defined in the presence of dispersion interactions. A smaller ion exhibits biexponential decay while a single exponential decay is seen for the ion with maximum diffusivity in the self-part of the intermediate scattering function. The solvent structure around the ion appears to determine much of the dynamics of the ion. Interesting trends are seen in the activation energies and these can be understood in terms of the levitation effect. (C) 2010 American Institute of Physics. doi:10.1063/1.3481656]

Effect of Sorbate-Zeolite Interaction on Cluster Lifetime and Size of Sorbates: Xe in NaY

Molecular dynamics calculations are reported for Xe in sodium Y zeolite with varying strengths of sorbate-zeolite dispersion interaction. In the absence of any dispersion interaction between the sorbate and the zeolite, the presence of the zeolite has a purely geometrical role. Increase in the strength of the sorbate-zeolite interaction increases the monomer population and decreases the population of dimers and higher sized clusters. The lifetime of the monomers as well as dimers increases with the strength of the dispersion interaction. The observed variations in the lifetime and the population of the different sized clusters is explained in terms of the changes in the potential energy surface caused by the increase in the strength of the dispersion interaction.

Molecular expression for dielectric friction on a rotating dipole: Reduction to the continuum theory

Recently we presented a microscopic expression for dielectric friction on a rotating dipole. This expression has a rather curious structure, involving the contributions of the transverse polarization modes of the solvent and also of the molecular length scale processes. It is shown here that under proper limiting conditions, this expression reduces exactly to the classical continuum model expression of Nee and Zwanzig [J. Chem. Phys. 52, 6353 (1970)]. The derivation requires the use of the asymptotic form of the orientation‐dependent total pair correlation function, the neglect of the contributions of translational modes of the solvent, and also the use of the limit that the size of the solvent molecules goes to zero. Thus, the derivation can be important in understanding the validity of the continuum model and can also help in explaining the results of a recent computer simulation study of dielectric relaxation in a Brownian dipolar lattice.

Diffusion of sorbates in zeolites Y and A novel dependence on sorbate size and strength of sor ba t e-Zeoli t e interaction

A systematic investigation of monatomic spherical sorbates in the supercages of zeolites Y and A by molecular dynamics technique is presented. Rates of intercage diffusion, rates of cage visits, and the diffusion coefficients have been calculated as a function of the sorbate-zeolite interaction strength. These properties exhibit markedly different dependences on interaction strength for the two zeolites. The observed behavior is shown to be a consequence of the two principal mechanisms of intercage diffusion and the energetic barrier associated with them. The diffusion coefficient and other properties associated with intercage diffusion are found to be directly proportional to the reciprocal of the square of the sorbate diameter when the sorbate diameter is significantly smaller than the window diameter. As the sorbate diameter increases, a peak is observed in all the transport properties investigated including the diffusion coefficient. We call this surprising effect as the ring or levitation effect and it explains several anomalous results reported in the literature and suggests a breakdown of the geometrical criterion for diffusion of sorbates. It shows that under certain conditions nongeometrical factors play a major role and geometrical factors become secondary in the determination of the molecular sieve property. A generalized parameter has been proposed which suggests conditions under which one can expect the ring or levitation effect in any porous medium. Inverse size selectivity becomes operative under these conditions.

Synthesis and characterization of barium bis(citrato) oxozirconate(iv) tetrahydrate--A new molecular precursor for fine particle BaZrO3

Barium metazirconate (BaZrO3) fine powder has been produced by thermally decomposing a molecular precursor, barium bis(citrato)oxozirconate(IV) tetrahydrate at about 700-degrees-C. The precursor, Ba[ZrO(C6H6O7)2] . 4H2O (BZO) has been synthesized and characterized by employing a combination of spectroscopic and thermoanalytical techniques. The precursor undergoes thermal decomposition in three major stages: (i) dehydration to give an anhydrous barium zirconyl citrate, (ii) decomposition of the anhydrous citrate in a multistep process to form an ionic oxycarbonate intermediate, Ba2Zr2O5CO3, and (iii) decomposition of the oxycarbonate to produce BaZrO3 fine powder. The particle size of the resultant BaZrO3 is about 0.2 mum, and the surface area is found to be 4.0 m2 g-1.

Molecular organization of the gene encoding Xenopus laevis transforming growth factor-beta 5

Transforming growth factor-beta s (TGF-beta 5) are multifunctional polypeptides, known to influence proliferation and differentiation of many cell types. TGF-beta 5 cDNA was cloned from Xenopus laevis and this isoform is unique to the amphibians. Here, we report the isolation and characterization of the TGF-beta 5 genomic clones to determine the structure of TGF-beta 5 gene. The gene consists of seven exons and all intron-exon boundaries follow the GT-AG consensus. The organization of TGF-beta 5 gene was identical to that of the mammalian TGF-beta isoforms, with the exception of position of the first splice junction. We determined the size of TGF-beta 5 gene to be approximately 20 kb.

Indium Nitride Nanometric-Objects on c-Sapphire Grown by Plasma-Assisted Molecular Beam Epitaxy

The indium nitride (InN)-based nanometric-objects were grown directly on a c-sapphire substrate by using plasma-assisted molecular beam epitaxy (PAMBE) at different substrate temperatures. High resolution X-ray diffraction (HRXRD) reveals the InN (0002) reflection and full width at half maximum (FWHM) found to be decreased with increasing the growth temperature. The size, height and density of the grown nanometric-objects studied by scanning electron microscopy (SEM) has remarkable differences, evidencing the decisive role of substrate temperature. Photoluminescence (PL) studies revealed that the emission energy is shifted towards the higher side from the bulk value, i.e., a blue shift in the PL spectra was observed. The temperature dependence of the PL peak position shows an ``S-shaped'' emission energy shift, which can be attributed to the localization of carriers in the nanometric-objects.

An ultra-stable Cu nanowire under thermo-mechanical loading: A molecular dynamic study

We report on the formation of a stable Body-Centered Heptahedral (BCH) crystalline nanobridge structure of diameter ~ 1nm under high strain rate tensile loading to a <100> Cu nanowire. Extensive Molecular Dynamics (MD) simulations are performed. Six different cross-sectional dimensions of Cu nanowires are analyzed, i.e. 0.3615 x 0.3615 nm2, 0.723 x 0.723 nm2, 1.0845 x 1.0845 nm2, 1.446 x 1.446 nm2, 1.8075 x 1.8075 nm2, and 2.169 x 2.169 nm2. The strain rates used in the present simulations are 1 x 109 s-1, 1 x 108 s-1, and 1 x 107 s-1. We have shown that the length of the nanobridge can be characterized by larger plastic strain. A large plastic deformation is an indication that the structure is highly stable. The BCH nanobridge structure also shows enhanced mechanical properties such as higher fracture toughness and higher failure strain. The effect of temperature, strain rate and size of the nanowire on the formation of BCH structure is also explained in details. We also show that the initial orientation of the nanowires play an important role on the formation of BCH crystalline structure. Results indicate that proper tailoring of temperature and strain rate during processing or in the device can lead to very long BCH nanobridge structure of Cu with enhanced mechanical properties, which may find potential application for nano-scale electronic circuits.

Temperature-pressure-induced solid-solid < 100 > to < 110 > reorientation in FCC metallic nanowire: a molecular dynamic study

Atomistic simulation of initial < 100 > oriented FCC Cu nanowires shows a novel coupled temperature-pressure dependent reorientation from < 100 > to < 110 > phase. A temperature-pressure-induced solid-solid < 100 > to < 110 > reorientation diagram is generated for Cu nanowire with varying cross-sectional sizes. A critical pressure is reported for Cu nanowires with varying cross-sectional sizes, above which an initial < 100 > oriented nanowire shows temperature independent reorientation into the < 110 > phase. The effect of surface stresses on the < 100 > to < 110 > reorientation is also studied. The results indicate that above a critical cross-sectional size for a given temperature-pressure, < 100 > to < 110 > reorientation is not possible. It is also reported here that for a given applied pressure, an increase in temperature is required for the < 100 > to < 110 > reorientation with increasing cross-sectional size of the nanowire. The temperature-pressure-induced solid-solid < 100 > to < 110 > reorientation diagram reported in the present paper could further be used as guidelines for controlling the reorientations/shape memory in nano-scale applications of FCC metallic nanowires.