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.

Structural Studies Of Analgesics And Their Interactions .1. Crystal And Molecular-Structure Of Antipyrine

The complex crystallizes in the space group P21/c with four formula units in a unit cell of dimensionsa= 12.747, b= 7.416, c= 17.894 A and/3= 90.2 °. The structure has been solved by the symbolic addition procedure using three dimensional photographic data and refined to an R value of 0.079 for 2019 observed reflexions. The pyramidal nature of the two hetero nitrogen atoms in the antipyrine molecule is inter:nediate between that observed in free antipyrine and in some of its metal complexes. The molecule is more polar than that in crystals of free antipyrine but less so compared with that in metal complexes. In the salicylic acid molecule, the hydroxyl group forms an internal hydrogen bond with one of the oxygen atoms in the carboxyl group. The association between the salicylic acid and the antipyrine molecules is achieved through an intermolecular hydrogen bond with the other carboxyl oxygen atom in the salicylic acid molecule as the proton donor and the carboxyl oxygen atom of the antipyrine molecule as the acceptor

Crystalline ethane-1,2-diol does not have intra-molecular hydrogen bonding: Experimental and theoretical charge density studies

The X-ray structure and electron density distribution of ethane-1,2-diol (ethylene glycol), obtained at a resolution extending to 1.00 Å−1 in sin θ/λ (data completion = 100% at 100 K) by in situ cryocrystallization technique is reported. The diol is in the gauche (g′Gt) conformation with the crystal structure stabilised by a network of inter-molecular hydrogen bonds. In addition to the well-recognized O–H···O hydrogen bonds there is topological evidence for C–H···O inter-molecular interactions. There is no experimental electron density based topological evidence for the occurrence of an intra-molecular hydrogen bond. The O···H spacing is not, vert, similar0.45 Å greater than in the gas-phase with an O–H···O angle close to 90°, calling into question the general assumption that the gauche conformation of ethane-1,2-diol is stabilised by the intra-molecular oxygen–hydrogen interaction.

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.

Molecular Basis for the Role of Staphylococcus aureus Penicillin Binding Protein 4 in Antimicrobial Resistance

Penicillin binding proteins (PBPs) are membrane-associated proteins that catalyze the final step of murein biosynthesis. These proteins function as either transpeptidases or carboxypeptidases and in a few cases demonstrate transglycosylase activity. Both transpeptidase and carboxypeptidase activities of PBPs occur at the D-Ala-D-Ala terminus of a murein precursor containing a disaccharide pentapeptide comprising N-acetyl-glucosamine and N-acetyl-muramic acid-L-Ala-D-Glu-L-Lys-D-Ala-D-Ala. beta-Lactam antibiotics inhibit these enzymes by competing with the pentapeptide precursor for binding to the active site of the enzyme. Here we describe the crystal structure, biochemical characteristics, and expression profile of PBP4, a low-molecular-mass PBP from Staphylococcus aureus strain COL. The crystal structures of PBP4-antibiotic complexes reported here were determined by molecular replacement, using the atomic coordinates deposited by the New York Structural Genomics Consortium. While the pbp4 gene is not essential for the viability of S. aureus, the knockout phenotype of this gene is characterized by a marked reduction in cross-linked muropeptide and increased vancomycin resistance. Unlike other PBPs, we note that expression of PBP4 was not substantially altered under different experimental conditions, nor did it change across representative hospital- or community-associated strains of S. aureus that were examined. In vitro data on purified recombinant S. aureus PBP4 suggest that it is a beta-lactamase and is not trapped as an acyl intermediate with beta-lactam antibiotics. Put together, the expression analysis and biochemical features of PBP4 provide a framework for understanding the function of this protein in S. aureus and its role in antimicrobial resistance.

Two-Phase Thermodynamic Model for Efficient and Accurate Absolute Entropy of Water from Molecular Dynamics Simulations

Presented here is the two-phase thermodynamic (2PT) model for the calculation of energy and entropy of molecular fluids from the trajectory of molecular dynamics (MD) simulations. In this method, the density of state (DoS) functions (including the normal modes of translation, rotation, and intramolecular vibration motions) are determined from the Fourier transform of the corresponding velocity autocorrelation functions. A fluidicity parameter (f), extracted from the thermodynamic state of the system derived from the same MD, is used to partition the translation and rotation modes into a diffusive, gas-like component (with 3Nf degrees of freedom) and a nondiffusive, solid-like component. The thermodynamic properties, including the absolute value of entropy, are then obtained by applying quantum statistics to the solid component and applying hard sphere/rigid rotor thermodynamics to the gas component. The 2PT method produces exact thermodynamic properties of the system in two limiting states: the nondiffusive solid state (where the fluidicity is zero) and the ideal gas state (where the fluidicity becomes unity). We examine the 2PT entropy for various water models (F3C, SPC, SPC/E, TIP3P, and TIP4P-Ew) at ambient conditions and find good agreement with literature results obtained based on other simulation techniques. We also validate the entropy of water in the liquid and vapor phases along the vapor-liquid equilibrium curve from the triple point to the critical point. We show that this method produces converged liquid phase entropy in tens of picoseconds, making it an efficient means for extracting thermodynamic properties from MD simulations.

Electronic structures of electron donor-acceptor complexes: results from ultraviolet photoelectron spectroscopy and molecular orbital calculations

HeI photoelectron spectra of 1:1 electron donor-acceptor complexes are discussed in the light of molecular orbital calculations. The complexes discussed include those formed by BH3, BF3 and SO2. Some systematics have been found in the ionization energy shifts of the complexes compared to the free components and these are related to the strength of the donor-acceptor bond. Hel spectra of hydrogen bonded complexes are discussed in comparison with results from MO calculations. Limitations of such studies as well as scope for further investigations are indicated.

Denoising signals corrupted by chaotic noise

We consider the problem of signal estimation where the observed time series is modeled as y(i) = x(i) + s(i) with {x(i)} being an orbit of a chaotic self-map on a compact subset of R-d and {s(i)} a sequence in R-d converging to zero. This model is motivated by experimental results in the literature where the ocean ambient noise and the ocean clutter are found to be chaotic. Making use of observations up to time n, we propose an estimate of s(i) for i < n and show that it approaches s(i) as n -> infinity for typical asymptotic behaviors of orbits. (C) 2010 Elsevier B.V. All rights reserved.

Estimation of signals in colored non gaussian noise based on Gaussian Mixture models

Non-Gaussianity of signals/noise often results in significant performance degradation for systems, which are designed using the Gaussian assumption. So non-Gaussian signals/noise require a different modelling and processing approach. In this paper, we discuss a new Bayesian estimation technique for non-Gaussian signals corrupted by colored non Gaussian noise. The method is based on using zero mean finite Gaussian Mixture Models (GMMs) for signal and noise. The estimation is done using an adaptive non-causal nonlinear filtering technique. The method involves deriving an estimator in terms of the GMM parameters, which are in turn estimated using the EM algorithm. The proposed filter is of finite length and offers computational feasibility. The simulations show that the proposed method gives a significant improvement compared to the linear filter for a wide variety of noise conditions, including impulsive noise. We also claim that the estimation of signal using the correlation with past and future samples leads to reduced mean squared error as compared to signal estimation based on past samples only.

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.

Comparative study of filter-bank mean-energy distance for automated segmentation of speech signals

This paper describes a method of automated segmentation of speech assuming the signal is continuously time varying rather than the traditional short time stationary model. It has been shown that this representation gives comparable if not marginally better results than the other techniques for automated segmentation. A formulation of the 'Bach' (music semitonal) frequency scale filter-bank is proposed. A comparative study has been made of the performances using Mel, Bark and Bach scale filter banks considering this model. The preliminary results show up to 80 % matches within 20 ms of the manually segmented data, without any information of the content of the text and without any language dependence. 'Bach' filters are seen to marginally outperform the other filters.

Am unexpected entry into the silver(I) chemistry of diphosphazane ligands: Crystal and molecular structures of [Ag{Ph2PN(Pr-i)PPh2}(CF3SO3)](2), [Ag{Ph2PN(Pr-i)PPh(OC6H3Me2-2,6)}(2)]PF6, and [Ag-3(mu-Cl)(2){Ph2PN(R)PPh2}(3)]PF6 (R = H, Pr-i)

Mononuclear, binuclear and trinuclear silver(l) complexes were obtained unexpectedly while probing the reactivity of diphosphazane ligands of the type X2PN(Pr-i)PXY towards the ruthenium-based precursor Ru(bipy)(2)Cl-2 center dot 2H(2)O, in the presence of a silver salt as a chloride scavenger. Subsequently, the reactions of AgX [X = Cl, NO3 or CF3SO3] with Ph2PN(R)PPh(Y) [R = H, Y = Ph; R = Pr-i, Y = Ph or OC6H3Me2-2,6] in a 1: 1 or 1:2 molar ratio have been investigated. Mononuclear or binuclear Ag(I) complexes containing either chelating or bridging diphosphazane ligands are obtained. Trinuclear silver(l) complexes are accessible by the treatment of diphosphazane ligands, Ph2PN(R)PPh2 [R = H, Pr-i] with AgCl using piperidine as the solvent. In the presence of a suitable chloride donor species, the mononuclear and binuclear complexes of Ph2PN(Pr-i)PPh2 are transformed slowly to the trinuclear complex [Ag-3(mu-Cl)(2){Ph2PN(Pr-i)PPh2}(3)]X, over a period 20 h. The structures of representative complexes have been confirmed by X-ray crystallography and the salient structural features are discussed