The shoving model for the glass-former LiCl center dot 6H(2)O: A molecular dynamics simulation study

Molecular dynamics (MD) simulations of LiCl center dot 6H(2)O Showed that the diffusion coefficient D, and also I lie structural relaxation time , follow a power law at high temperatures, D(-1) proportional to (T - T(0))(-mu), with the same experimental parameters for viscosity (T(0) = 207 K, mu = 2.08). Decoupling between D and occurs at T(x) similar to 1.1 T(0). High frequency acoustic excitations for the LiCl center dot 6H(2)O model were obtained by the calculation of time correlation functions of mass current fluctuations. The temperature dependence of the instantaneous shear modulus, G,(T), was considered in the shoving model for supercooled liquids [J.C. Dyre, T. Christensen, N.B. Olsen, J. Non-Cryst. Solids 352 (2006) 4635] resulting in a linear relationship log (D(-1)) vs. G root T. (C) 2009 Elsevier B.V. All rights reserved.

Prigogine-Defay Ratio for an Ionic Glass-Former: Molecular Dynamics Simulations

The pressure dependence of the glass-transition temperature, T(g)(P), of the ionic glass-former 2Ca(NO(3))(2) center dot 3KNO(3), CKN, has been obtained by molecular dynamics (MD) simulations The liquid-glass difference of thermal expansivity, Delta alpha, heat capacity, Delta C(p), and isothermal compressibility, Delta kappa, have been calculated as a function of pressure. It has been found that the Ehrenfest relation dT(g)/dP = TV Delta alpha/Delta C(p) predicts the pressure dependence of T, but the other Ehrenfest relation, dT(g)/dP = Delta kappa/Delta alpha, does not. Consequently, the Prigogine-Defay ratio, Pi = Delta C(p)Delta kappa/TV Delta alpha(2), is Pi similar to 1.2 at low pressures, but increases 1 order of magnitude at high pressures. The pressure dependence of the Prigogine-Defay ratio is interpreted in light of recent explanations for the finding Pi > 1.

Alkoxy Chain Effect on the Viscosity of a Quaternary Ammonium Ionic Liquid: Molecular Dynamics Simulations

The viscosity of ionic liquids based on quaternary ammonium cations is reduced when one of the alkyl chains is replaced by an alkoxy chain (Zhou et al. Chem. Eur. J. 2005, 11, 752.). A microscopic picture of the role played by the ether function in decreasing the viscosity of quaternary ammonium ionic liquids is provided here by molecular dynamics (MD) simulations. A model for the ionic liquid N-ethyl-N,N-dimethyl-N-(2-methoxyethyl)ammonium bis(trifluoromethanesulfonyl)imide, MOENM(2)E TFSI, is compared to the tetraalky-lammonium counterpart. The alkoxy derivative has lower viscosity, higher ionic diffusion coefficients, and higher conductivity than the tetraalkyl system at the same density and temperature. A clear signature of the ether function on the liquid structure is observed in cation-cation correlations, but not in anion-anion or anion-cation correlations. In both the alkyl and the alkoxy ionic liquids, there is aggregation of long chains of neighboring cations within micelle-like structures. The MD simulations indicate that the less effective assembly between the more flexible alkoxy chains, in comparison to alkyl chains, is the structural reason for higher ionic mobility in MOENM(2)E TFSI.

Formation and decomposition of N-alkyinaphthalimides: experimental evidences and ab initio description of the reaction pathways

The kinetics of hydrolysis of 1,8-N-butyl-naphthalimide (1,8-NBN) to 1,8-N-butyl-naphthalamide (1,8-NBAmide) and of 2,3-N-butyl-naphthalimide (2,3-NBN) to 2,3-N-butyl-naphthalamide (2,3-NBAmide), as well as the formation of the respective anhydrides from the amides were investigated in a wide acidity range. 1,8-NBN equilibrates with 1,8-NBAmide in mild alkali. Under the same conditions 2,3-NBN quantitatively yields 2,3-NBAmide. Over a wide range of acidities the reactions of the 1,8- and 2,3-N-butyl-naphthalamides (or imides) yield similar products but with widely different rates and at distinct pH`s. Anhydride formation in acid was demonstrated for 1,8-NBAmide. The reactions mechanisms were rationalized in the manifold pathways of ab initio calculations. The differences in rates and pH ranges in the reactions of the 1,8- and 2,3-N-butyl-naphthalamides were attributed to differences in the stability of the tetrahedral intermediates in alkali as well as the relative stabilities of the five and six-membered ring intermediates. The rate of carboxylic acid assisted 1,8-N-Butyl-naphthalamide hydrolysis is one of the largest described for amide hydrolysis models. Copyright (C) 2010 John Wiley & Sons, Ltd.

Docking and molecular dynamics simulation of quinone compounds with trypanocidal activity

In this work, two different docking programs were used, AutoDock and FlexX, which use different types of scoring functions and searching methods. The docking poses of all quinone compounds studied stayed in the same region in the trypanothione reductase. This region is a hydrophobic pocket near to Phe396, Pro398 and Leu399 amino acid residues. The compounds studied displays a higher affinity in trypanothione reductase (TR) than glutathione reductase (GR), since only two out of 28 quinone compounds presented more favorable docking energy in the site of human enzyme. The interaction of quinone compounds with the TR enzyme is in agreement with other studies, which showed different binding sites from the ones formed by cysteines 52 and 58. To verify the results obtained by docking, we carried out a molecular dynamics simulation with the compounds that presented the highest and lowest docking energies. The results showed that the root mean square deviation (RMSD) between the initial and final pose were very small. In addition, the hydrogen bond pattern was conserved along the simulation. In the parasite enzyme, the amino acid residues Leu399, Met400 and Lys402 are replaced in the human enzyme by Met406, Tyr407 and Ala409, respectively. In view of the fact that Leu399 is an amino acid of the Z site, this difference could be explored to design selective inhibitors of TR.

Distribution of neutral prilocaine in a phospholipid bilayer: Insights from molecular dynamics simulations

In this work, we report a 20-ns constant pressure molecular dynamics simulation of prilocaine (PLC), in amine-amide local anesthetic, in a hydrated liquid crystal bilayer of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine. The partition of PLC induces the lateral expansion of the bilayer and a concomitant contraction in its thickness. PLC molecules are preferentially found in the hydrophobic acyl chains region, with a maximum probability at similar to 12 angstrom from the center of the bilayer (between the C(4) and C(5) methylene groups). A decrease in the acyl chain segmental order parameter, vertical bar S-CD vertical bar, compared to neat bilayers, is found, in good agreement with experimental H-2-NMR studies. The decrease in vertical bar S-CD vertical bar induced by PLC is attributed to a larger accessible volume per lipid in the acyl chain region. (C) 2008 Wiley Periodicals, Inc.

Preferential location of prilocaine and etidocaine in phospholipid bilayers: A molecular dynamics study

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Cálculos AB initio de propriedades estruturais, eletrônicas e ópticas de colmpostos CaXo3 (X=Si, Ge, Sn)

In this work we present a study of structural, electronic and optical properties, at ambient conditions, of CaSiO3, CaGeO3 and CaSnO3 crystals, all of them a member of Ca-perovskite class. To each one, we have performed density functional theory ab initio calculations within LDA and GGA approximations of the structural parameters, geometry optimization, unit cell volume, density, angles and interatomic length, band structure, carriers effective masses, total and partial density of states, dielectric function, refractive index, optical absorption, reflectivity, optical conductivity and loss function. A result comparative procedure was done between LDA and GGA calculations, a exception to CaSiO3 where only LDA calculation was performed, due high computational cost that its low symmetry crystalline structure imposed. The Ca-perovskite bibliography have shown the absence of electronic structure calculations about this materials, justifying the present work

Sistemas Nanoestruturados: Heteroestruturas Quasi-Periódicas de Nitretos e Cálculos Ab initio em Polimorfos CaCO3

The physical properties and the excitations spectrum in oxides and semiconductors materials are presented in this work, whose the first part presents a study on the confinement of optical phonons in artificial systems based on III-V nitrides, grown in periodic and quasiperiodic forms. The second part of this work describes the Ab initio calculations which were carried out to obtain the optoeletronic properties of Calcium Oxide (CaO) and Calcium Carbonate (CaCO3) crystals. For periodic and quasi-periodic superlattices, we present some dynamical properties related to confined optical phonons (bulk and surface), obtained through simple theories, such as the dielectric continuous model, and using techniques such as the transfer-matrix method. The localization character of confined optical phonon modes, the magnitude of the bands in the spectrum and the power laws of these structures are presented as functions of the generation number of sequence. The ab initio calculations have been carried out using the CASTEP software (Cambridge Total Sequential Energy Package), and they were based on ultrasoft-like pseudopotentials and Density Functional Theory (DFT). Two di®erent geometry optimizations have been e®ectuated for CaO crystals and CaCO3 polymorphs, according to LDA (local density approximation) and GGA (generalized gradient approximation) approaches, determining several properties, e. g. lattice parameters, bond length, electrons density, energy band structures, electrons density of states, e®ective masses and optical properties, such as dielectric constant, absorption, re°ectivity, conductivity and refractive index. Those results were employed to investigate the confinement of excitons in spherical Si@CaCO3 and CaCO3@SiO2 quantum dots and in calcium carbonate nanoparticles, and were also employed in investigations of the photoluminescence spectra of CaCO3 crystal

Structural phase transition study of FePt alloys using ab initio calculation

The FePt alloy undergoes the cubic to tetragonal lattice transformation in the ferromagnetic state. We calculated the electronic structure for both cubic and tetragonal structures using the FPLAPW method with APW + lo. Comparing the density of states of the cubic and tetragonal structures, it is expected that the lattice transformation is caused by the band Jahn-Teller effect. (C) 2009 Elsevier B.V. All rights reserved.

Estudo da adsorção de hidrogênio e sulfeto na superfície de paládio: aspectos experimentais (eletroquímica) e teóricos (ab initio e Teoria do Funcional da Densidade)

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Estudo ab-initio da a-alanina em meio aquoso

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

An ab initio study of oxygen vacancies and doping process of Nb and Cr atoms on TiO2 (110) surface models

We theoretically investigated how the formation of oxygen vacancies and the addition of niobium and chromium atoms as dopants modify the varistor properties of TiO2. The calculations were carried out at the HF level using a contracted basis set, developed by Huzinaga et al.. to represent the atomic centers on the (110) surface for the large (TiO2)(15) cluster model. The change of the values for the net atomic charges and band gap after oxygen vacancy formation and the presence of dopants in the lattice are analyzed and discussed. It is shown that the formation of oxygen vacancies decreases the band gap while an opposite effect is found when dopants are located in the reduced surface. The theoretical results are compared with available experimental data. A plausible explanation of the varistor behavior of this system is proposed. (C) 1997 John Wiley & Sons, Inc.