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.

Thermodynamic Properties and Phase Diagram for the System MoO2–TiO2

The activity of molybdenum dioxide (MoO2) in the MoO2–TiO2 solid solutions was measured at 1600 K using a solid-state cell incorporating yttria-doped thoria as the electrolyte. For two compositions, the emf was also measured as a function of temperature. The cell was designed such that the emf is directly related to the activity of MoO2 in the solid solution. The results show monotonic variation of activity with composition, suggesting a complete range of solid solutions between the end members and the occurrence of MoO2 with a tetragonal structure at 1600 K. A large positive deviation from Raoult's law was found. Excess Gibbs energy of mixing is an asymmetric function of composition and can be represented by the subregular solution model of Hardy as follows.The temperature dependence of the emf for two compositions is reasonably consistent with ideal entropy of mixing. A miscibility gap is indicated at a lower temperature with the critical point characterized by Tc (K)=1560 and . Recent studies indicate that MoO2 undergoes a transition from a monoclinic to tetragonal structure at 1533 K with a transition entropy of 9.91 J·(mol·K)−1. The solid solubility of TiO2 with rutile structure in MoO2 with a monoclinic structure is negligible. These features give rise to a eutectoid reaction at 1412 K. The topology of the computed phase diagram differs significantly from that suggested by Pejryd.

Thermodynamic activities in the Pb(Zr1-XTiX)O-3 solid solution at 1373 K

Although Pb(Zr1-XTiX)O-3 solid solution is the cornerstone of the piezoelectric ceramics, there is no information in the literature on thermodynamic activities of the component phases in the solid solution. Using inter-crystalline ion exchange equilibria between Pb(Zr1-XTiX)O-3 solid solution with cubic perovskite structure and (Zr1-YTiY)O-2 solid solutions with monoclinic and tetragonal structures, activities of PbTiO3 and PbZrO3 in the perovskite solid solution have been derived at 1373 K using the modified Gibbs-Duhem integration technique of Jacob and Jeffes. Tie-lines from the cubic solid solution are skewed towards the ZrO2 corner. Activities in the zirconia-rich (Zr1-YTiY)02 solid solutions are taken from a recent emf study. The results for the perovskite solid solution at 1373 K can be represented by a sub-regular solution model: Delta G(E.M) (J mol(-1)) = X-PbTiO3 X-PbZrO3(5280X(PbTiO3) - 1980X(PbZrO3)) where Delta G(E.M) is the excess Gibbs energy of mixing of the cubic solid solution and Xi represents the mole fraction of component i. There is a significant positive deviation from ideality for PbTiO3-rich compositions and mild negative deviation near the PbZrO3 corner. The cubic solid solution is intrinsically stable against composition fluctuations at temperatures down to 840 K. The results contrast sharply with the recent calorimetric data on enthalpy of mixing which signal instability of the cubic perovskite solid solution. (C) 2007 Elsevier B.V. All rights reserved.

A study of the nanostructure of the cell wall of the tracheids of conifer xylem by x-ray scattering

The purpose of this study was to develop practical and reliable x-ray scattering methods to study the nanostructure of the wood cell wall and to use these methods to systematically study the nanostructure of Norway spruce and Scots pine grown in Finland and Sweden. Methods to determine the microfibril angle (MFA) distribution, the crystallinity of wood, and the average size of cellulose crystallites using wide-angle x-ray scattering were developed and these parameters were determined as a function of the number of the year ring. The mean MFA in Norway spruce decreases rapidly as a function of the number of the year ring and after the 7th year ring it varies between 6° and 10°. The mean MFA of Scots pine behaves the same way as the mean MFA of Norway spruce. The thickness of cellulose crystallites for Norway spruce and Scots pine appears to be constant as a function of the number of the year ring. The obtained mean values are 32 Å for Norway spruce and 31 Å for Scots pine. The length of the cellulose crystallites was also quite constant as a function of the year ring. The mean length of the crystallites for Norway spruce was 364 Å, while the standard deviation was 27 Å. The mass fraction of crystalline cellulose in wood is the crystallinity of wood and the intrinsic crystallinity of cellulose is the crystallinity of cellulose. The crystallinity of wood increases from the 2nd year ring to the 10th year ring from the pith and is constant after the 10th year ring. The crystallinity of cellulose obtained using nuclear magnetic resonance spectroscopy was 52% for both species. The crystallinity of wood and the crystallinity of cellulose behave the same way in Norway spruce and Scots pine. The methods were also applied to studies on thermally modified Scots pine wood grown in Finland. Wood is modified thermally by heating and steaming in order to improve its properties such as biological resistance and dimensional stability. Modification temperatures varied from 100 °C to 240 °C. The thermal modification increases the crystallinity of wood and the thickness of cellulose crystallites but does not influence the MFA distribution. When the modification temperature was 230 °C and time 4 h, the thickness of the cellulose crystallites increased from 31 Å to 34 Å.

A Fair Contract for Managing Water Scarcity

In public utilities, under supply constraints, fairness considerations lead to a market failure. This paper characterizes a two-period principal-agent contract for demand management, that mitigates this market failure in urban water systems. The contract is designed as an extensive form mechanism using subgame perfect Nash equilibrium (SPNE) as the solution concept. The contract is fair; and is shown to be economically efficient if, in case of deviation by the agent, the gain to the agent and the loss to the principal are small. It is shown that the assumption can be avoided in an infinite horizon contract.

Metal-spinel-corundum three phase equilibria in the system Ni-Cr-Al-O at 1373 K

The tie-lines representing the inter-crystalline ion exchange equilibria between the NiCr2O4-NiAl2O4 spinet solid solution and Cr2O3-Al2O3 corundum solid solution are determined by electron microprobe andEDAX pointcountanalysis of the oxide phases equilibrated with metallic Ni at 1373 K. The component activities in the spinet solid solution are derived from the tie-lines and thermodynamic data for Cr2O3-Al2O3 solid solution available in the literature. The Gibbs energy of mixing of the spinet solid solution calculated from the experimental data is discussed in relation to the values derived from the cation distribution models which assume random mixing of cations on both tetrahedral and octahedral sites. Positive deviation from the models is observed indicating significant positive enthalpy contribution arising form the size mismatch between Al+3 and Ni+2 ions on the tetrahedral site and Al+3, Ni+2 and Cr+3 on the octahedral site. Variation of the oxygen potential for threephase equilibrium involving metallic nickel, spinet solid solution and corundum solid solution is computed as a function of composition of the solid solutions at 1373 K. The oxygen potential exhibits a minimum at aluminum cationic fraction eta(Al)/(eta(Al) + eta(Cr)) = 0.524 in the oxide solid solutions.

(2-Chlorobenzo[h]quinolin-3-yl)methanol

In the title molecule, C14H10ClNO, all non-H atoms are coplanar (r.m.s deviation = 0.0266 angstrom). In the crystal, symmetry-related molecules are hydrogen bonded via intermolecular O-H center dot center dot center dot O interactions, forming chains along the b axis.

1-[(2-Chloro-7-methyl-3-quinolyl)methyl]pyridin-2(1H)-one

In the title compound, C16H13ClN2O, the quinoline ring system is essentially planar, with a maximum deviation of 0.021 (2) angstrom. The pyridone ring is oriented at a dihedral angle of 85.93 (6)degrees with respect to the quinoline ring system. In the crystal structure, intermolecular C-H center dot center dot center dot O hydrogen bonds link the molecules along the b axis. Weak pi-pi stacking interactions [centroid-centroid distances = 3.7218 (9) and 3.6083 (9) angstrom] are also observed.

1-[(2-Chloro-7,8-dimethylquinolin-3-yl)methyl]pyridin-2(1H)-one

In the title compound, C17H15ClN2O, the quinoline ring system is nearly planar, with a maximum deviation from the mean plane of 0.074 (2) angstrom, and makes a dihedral angle of 81.03 (7)degrees with the pyridone ring. The crystal packing is stabilized by pi-pi stacking interactions between the pyridone and benzene rings of the quinoline ring system [centroid-centroid distance = 3.6754 (10) angstrom]. Furthermore, weak intermolecular C-H center dot center dot center dot O hydrogen bonding links molecules into supramolecular chains along [001].

Monoclinic polymorph of Boc-Trp-Ile-Ala-Aib-Ile-Val-Aib-Leu-Aib-Pro-OMe(anhydrous). Parallel packing of 3(10)-/alpha-helices and a transition of helix type

The structures of two crystal forms of Boc-Trp-Ile-Ala-Aib-Ile-Val-Aib-Leu-Aib-Pro-OMe have been determined. The triclinic form (P1, Z = 1) from DMSO/H2O crystallizes as a dihydrate (Karle, Sukumar & Balaram (1986) Proc, Natl, Acad. Sci. USA 83, 9284-9288). The monoclinic form (P2(1), Z = 2) crystallized from dioxane is anhydrous. The conformation of the peptide is essentially the same in both crystal system, but small changes in conformational angles are associated with a shift of the helix from a predominantly alpha-type to a predominantly 3(10)-type. The r.m.s. deviation of 33 atoms in the backbone and C beta positions of residues 2-8 is only 0.29 A between molecules in the two polymorphs. In both space groups, the helical molecules pack in a parallel fashion, rather than antiparallel. The only intermolecular hydrogen bonding is head-to-tail between helices. There are no lateral hydrogen bonds. In the P2(1) cell, a = 9.422(2) A, b = 36.392(11) A, c = 10.548(2) A, beta = 111.31(2) degrees and V = 3369.3 A for 2 molecules of C60H97N11O13 per cell.

Dehydrogenation of Butyl Alcohol in Fixed Catalyst Beds

The vapor-phase dehydrogenation of 1 -butanol to butyraldehyde was studied in a fixed bed of catalyst from 250° to 360° C. Of all the catalysts studied during preliminary investigation, the one containing 90% copper, 8% chromia, and 2% carbon supported on pumice was best, with high activity and selectivity. The data are expressed in the form of a first-order irreversible reaction rate equation. Single-site surface reaction (hydrogen adsorbed) is the rate-controlling mechanism at all the temperatures studied. The rate data obtained in the entire range of experimental conditions fit the rate equation based on this mechanism with a standard deviation of ± 22.8%.

Correlation of thermodynamic properties of chloromethanes

A correlation for entropy and enthalpy based on the two-parameter law of corresponding states and the reference substance concept has been obtained. For 199 points tested, the over-all average absolute deviation and the maximum deviations in the calculated values of entropy with the available data are 0.74 and 7.20%. The corresponding deviations of enthalpy are 1.86 and 15.0%, respectively. A compressibility chart for chloromethanes has been made and shown to be superior to existing charts. For 102 points tested, the average absolute and maximum deviations in the compressibilities were 1.80 and 19.5%, respectively.

Maan magneettikentän dipolimallin toimivuus inklinaatiojakauman perusteella

The magnetic field of the Earth is 99 % of the internal origin and generated in the outer liquid core by the dynamo principle. In the 19th century, Carl Friedrich Gauss proved that the field can be described by a sum of spherical harmonic terms. Presently, this theory is the basis of e.g. IGRF models (International Geomagnetic Reference Field), which are the most accurate description available for the geomagnetic field. In average, dipole forms 3/4 and non-dipolar terms 1/4 of the instantaneous field, but the temporal mean of the field is assumed to be a pure geocentric axial dipolar field. The validity of this GAD (Geocentric Axial Dipole) hypothesis has been estimated by using several methods. In this work, the testing rests on the frequency dependence of inclination with respect to latitude. Each combination of dipole (GAD), quadrupole (G2) and octupole (G3) produces a distinct inclination distribution. These theoretical distributions have been compared with those calculated from empirical observations from different continents, and last, from the entire globe. Only data from Precambrian rocks (over 542 million years old) has been used in this work. The basic assumption is that during the long-term course of drifting continents, the globe is sampled adequately. There were 2823 observations altogether in the paleomagnetic database of the University of Helsinki. The effect of the quality of observations, as well as the age and rocktype, has been tested. For comparison between theoretical and empirical distributions, chi-square testing has been applied. In addition, spatiotemporal binning has effectively been used to remove the errors caused by multiple observations. The modelling from igneous rock data tells that the average magnetic field of the Earth is best described by a combination of a geocentric dipole and a very weak octupole (less than 10 % of GAD). Filtering and binning gave distributions a more GAD-like appearance, but deviation from GAD increased as a function of the age of rocks. The distribution calculated from so called keypoles, the most reliable determinations, behaves almost like GAD, having a zero quadrupole and an octupole 1 % of GAD. In no earlier study, past-400-Ma rocks have given a result so close to GAD, but low inclinations have been prominent especially in the sedimentary data. Despite these results, a greater deal of high-quality data and a proof of the long-term randomness of the Earth's continental motions are needed to make sure the dipole model holds true.

Thermal Conductivities of Organic Liquids-A New Correlation.

A new equation for predicting the thermal conductivities of organic liquids using dimension-less analysis is given. The equation (Equation Presented) correlates 51 different liquids tested within 11% average error and 17% standard deviation. A comparison of the proposed equation with the available correlations and its application to some industrially important liquids show that this equation can be safely used to calculate the thermal conductivities at 20°C. and 1 atm. pressure for organic liquids of known molecular weight. Cp and ΔHv - the only two parameters for which experimental values must be known for making use of this equation - can be calculated using other well known correlations. The proposed equation is not applicable to inorganic liquids.

Fugacity coefficients. Enthalpy, and second virial coefficients of boron trifluoride

Fugacity coefficients and isothermal changes of enthalpy have been calculated and reported. The calculations cover a temperature range of 0° to 75°C. up to gas densities of 1.0 gram per cc. The generalized Benedict-Webb-Rubin constants evaluated from generalized PVT relations is found to predict the experimental data with an over-all absolute deviation of 3.1%. Second virial coefficients and potential energy parameters for Lennard-Jones (12-6) potential energy function are reported also.