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.

Long-range synthon Aufbau modules (LSAM) in crystal structures: systematic changes in C6H6-nFn (0 <= n <= 6) fluorobenzenes

We discuss the assembly of a three-dimensional molecular crystal in terms of short-range supramolecular synthons that spontaneously organize themselves according to Aufbau principles into long-range geometries characteristic of the molecules themselves. For this purpose we have examined the systematic changes in the known crystal structures of a family of fluorobenzenes, C6H6-nFn, where 0 <= n <= 6. Crystal assembly is initiated by forming long-range synthon Aufbau modules (LSAM) that carry the imprint of the synthons. For example, when 1 <= n <= 5 the short-range synthons use H center dot center dot center dot F interactions to form the LSAMs. In the n = 0 and n = 6 compounds, the synthons are H center dot center dot center dot C and F center dot center dot center dot C interactions, respectively. The LSAMs are usually one-dimensional. In this study we show that these 1D LSAMs assemble into 2D quasi-hexagonal close-packed layers. The 3D crystal structure is obtained from the various kinds of close-packing known for these 2D layers. The final stages of this 1D -> 2D -> 3D assembly seem to be more influenced by the packing of LSAMs than by any other factor. In these final stages, there may not be so much influence exerted by the stronger short-range synthons. We discuss the evolution of these fluorobenzene crystal structures in terms of putative LSAMs and the purely geometric relationships between the n and (6 - n) compounds that can thus be expected. Such particle-hole pairs show structural similarities. Our discussion is quantified by the interpretation of intermolecular distances in terms of atomic sizes and with qualitative predictions of magnetic model systems.

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.

Long-range synthon Aufbau modules (LSAM) in crystal structures: systematic changes in C6H6-nFn (0 <= n <= 6) fluorobenzenes

We discuss the assembly of a three-dimensional molecular crystal in terms of short-range supramolecular synthons that spontaneously organize themselves according to Aufbau principles into long-range geometries characteristic of the molecules themselves. For this purpose we have examined the systematic changes in the known crystal structures of a family of fluorobenzenes, C6H6-nFn, where 0 <= n <= 6. Crystal assembly is initiated by forming long-range synthon Aufbau modules (LSAM) that carry the imprint of the synthons. For example, when 1 <= n <= 5 the short-range synthons use H center dot center dot center dot F interactions to form the LSAMs. In the n = 0 and n = 6 compounds, the synthons are H center dot center dot center dot C and F center dot center dot center dot C interactions, respectively. The LSAMs are usually one-dimensional. In this study we show that these 1D LSAMs assemble into 2D quasi-hexagonal close-packed layers. The 3D crystal structure is obtained from the various kinds of close-packing known for these 2D layers. The final stages of this 1D -> 2D -> 3D assembly seem to be more influenced by the packing of LSAMs than by any other factor. In these final stages, there may not be so much influence exerted by the stronger short-range synthons. We discuss the evolution of these fluorobenzene crystal structures in terms of putative LSAMs and the purely geometric relationships between the n and (6 - n) compounds that can thus be expected. Such particle-hole pairs show structural similarities. Our discussion is quantified by the interpretation of intermolecular distances in terms of atomic sizes and with qualitative predictions of magnetic model systems.

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.

Type I seesaw mechanism for quasi-degenerate neutrinos

We discuss symmetries and scenarios leading to quasi-degenerate neutrinos in type I seesaw models. The existence of degeneracy in the present approach is not linked to any specific structure for the Dirac neutrino Yukawa coupling matrix y(D) and holds in general. Basic input is the application of the minimal flavour violation principle to the leptonic sector. Generalizing this principle, we assume that the structure of the right-handed neutrino mass matrix is determined by y(D) and the charged lepton Yukawa coupling matrix y(l) in an effective theory invariant under specific groups G(F) contained in the full symmetry group of the kinetic energy terms. G(F) invariance also leads to specific structure for the departure from degeneracy. The neutrino mass matrix (with degenerate mass m(0)) resulting after seesaw mechanism has a simple form Mv approximate to m(0)(I - py(l)y(l)(T)) in one particular scenario based on supersymmetry. This form is shown tolead to correct description of neutrino masses and mixing angles. The thermal leptogenesis after inclusion of flavour effects can account for the observed baryon asymmetry of the universe within the present scenario. Rates for lepton flavour violating processes can occur at observable levels in the supersymmetric version of the scenario. (c) 2010 Elsevier B.V. All rights reserved.

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.

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

Effect of coordinated ions on structure and flexibility of parallel G-quadruplexes: A molecular dynamics study

ingle tract guanine residues can associate to form stable parallel quadruplex structures in the presence of certain cations. Nanosecond scale molecular dynamics simulations have been performed on fully solvated fibre model of parallel d(G(7)) quadruplex structures with Na+ or K+ ions coordinated in the cavity formed by the O6 atoms of the guanine bases. The AMBER 4.1 force field and Particle Mesh Ewald technique for electrostatic interactions have been used in all simulations. There quadruplex structures are stable during the simulation, with the middle four base tetrads showing root mean square deviation values between 0.5 to 0.8 Angstrom from the initial structure as well the high resolution crystal structure. Even in the absence of any coordinated ion in the initial structure, the G-quadruplex structure remains intact throughout the simulation. During the 1.1 ns MD simulation, one Nai counter ion from the solvent as well as several water molecules enter the central cavity to occupy the empty coordination sites within the parallel quadruplex and help stabilize the structure. Hydrogen bonding pattern depends on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate cations of different sizes. in the absence of any coordinated ion. due to strong mutual repulsion, O6 atoms within G-tetrad are forced farther apart from each other, which leads to a considerably different hydrogen bonding scheme within the G-tetrads and very favourable interaction energy between the guanine bases constituting a G-tetrad. However, a coordinated ion between G-tetrads provides extra stacking energy for the G-tetrads and makes the quadruplex structure more rigid. Na+ ions, within the quadruplex cavity, are more mobile than coordinated K+ ions. A number of hydrogen bonded water molecules are observed within the grooves of all quadruplex structures.