ONIOM study of dissociated hydrogen and water on ZnO surface

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Study of water and dimethylformamide interaction by computer simulation

Monte Carlo simulations of water-dimethylformamide (DMF) mixtures were performed in the isothermal and isobaric ensemble at 298.15 K and 1 atm. The intermolecular interaction energy was calculated using the classical 6-12 Lennard-Jones pairwise potential plus a Coulomb term. The TIP4P model was used for simulating water molecules, and a six-site model previously optimised by us was used to represent DMF. The potential energy for the water-DMF interaction was obtained via standard geometric combining rules using the original potential parameters for the pure liquids. The radial distribution functions calculated for water-DMF mixtures show well characterised hydrogen bonds between the oxygen site of DMF and hydrogen of water. A structureless correlation curve was observed for the interaction between the hydrogen site of the carbonyl group and the oxygen site of water. Hydration effects on the stabilisation of the DMF molecule in aqueous solution have been investigated using statistical perturbation theory. The results show that energetic changes involved in the hydration process are not strong enough to stabilise another configuration of DMF than the planar one.

The kinetic of mullite crystallization: Effect of water content

The kinetic of mullite crystallization from sol–gel method, with different water content, was investigated under non-isothermal conditions using DTA. The sols were obtained from Al(NO3)3.9H2O (ANN) and Si(OC2H5)4 (TEOS) mixtures by varying the water–alcohol content of the system. The crystalline phase changes were verified by X-ray diffraction (XRD). For a sample prepared using ethanol-based alkoxide solution (M0), only Al-poor mullite (p-mullite) crystallizes at 1000 °C; for the one synthesized with low water concentration (M6) Al-rich mullite (r-mullite) and spinel crystallize together; and for a sample prepared using a water-based alkoxide solution only spinel is formed. Thus, the variation of water contents during the synthesis caused great variations in the course of mullitization process. The average value of the apparent activation energy determined for p-mullite, r-mullite and spinel phase crystallization were found to be E = (899 ± 61) kJ mol−1, E = (1015 ± 272) kJ mol−1 and E = (980 ± 196) kJ mol−1, respectively. These results showed that sample M(0) was a monophasic gel, where aluminum and silicon atoms are mixed at a molecular level while sample M(100) was a diphasic gel, where silicon and aluminum atoms are distributed in a nanometric level. The fast reaction between TEOS and water molecules is responsible for this great difference in the sample's homogeneity. The kinetic model of the crystallization process was determined using Malek's procedure. It was established that the crystallization of p-mullite, r-mullite and spinel phase can be described by Šesták–Berggren autocatalytic model.

Host-guest interactions between xanthones and water: the role of O-H center dot center dot center dot O, C-H center dot center dot center dot O, and pi center dot center dot center dot pi contacts in the channel- and cage-type frameworks

In this article were studied two xanthone derivatives known as 1,5-dihydroxy-8-methoxyxanthone (I) and 1,3,7-trihydroxy-8-methoxyxanthone (II), which show one water molecule into their crystal structures. In xanthone I, there are water wires contributing to build up channel-like cavities along the c axis, whereas in xanthone II the water is surrounded by three xanthone molecules forming a cage-type structure. The geometries of I and II were optimized using the density functional theory method with B3LYP functional, and the results were compared with crystal structure. Both theoretical and experimental investigations reveal a concordance between structural parameters, with the xanthone core presenting an almost flat conformation and substituents adopting the more stable orientations. In the two compounds, the hydroxyl group linked at position 1 is involved in a resonance-assisted hydrogen bond with the carbonyl group. Besides, the supramolecular arrangement of the host/guest systems are stabilized mainly by classical intermolecular hydrogen bonds (O-H center dot center dot center dot O) involving xanthone-to-water and xanthone-to-xanthone. In addition, C-H center dot center dot center dot O weak hydrogen bonds, as well as pi-pi interactions play an important role to stabilize the crystal self-assembly of xanthones I and II. The results reported here underline the role of inclusion of water molecules and their different arrangement into the crystal structure of two xanthone host/guest systems.

Hydrogen bond and the resonance effect on the formamide-water complexes

The interaction of formamide and the two transition states of its amide group rotation with one, two, or three water molecules was studied in vacuum. Great differences between the electronic structure of formamide in its most stable form and the electronic structure of the transition states were noticed. Intermolecular interactions were intense, especially in the cases where the solvent interacted with the amide and the carbonyl groups simultaneously. In the transition states, the interaction between the lone pair of nitrogen and the water molecule becomes important. With the aid of the natural bond orbitals, natural resonance theory, and electron localization function (ELF) analyses an increase in the resonance of planar formamide with the addition of successive water molecules was observed. Such observation suggests that the hydrogen bonds in the formamidewater complexes may have some covalent character. These results are also supported by the quantitative ELF analyses. (C) 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012

A simple analysis of the influence of the solvent-induced electronic polarization on the N-15 magnetic shielding of pyridine in water

Electronic polarization induced by the interaction of a reference molecule with a liquid environment is expected to affect the magnetic shielding constants. Understanding this effect using realistic theoretical models is important for proper use of nuclear magnetic resonance in molecular characterization. In this work, we consider the pyridine molecule in water as a model system to briefly investigate this aspect. Thus, Monte Carlo simulations and quantum mechanics calculations based on the B3LYP/6-311++G (d,p) are used to analyze different aspects of the solvent effects on the N-15 magnetic shielding constant of pyridine in water. This includes in special the geometry relaxation and the electronic polarization of the solute by the solvent. The polarization effect is found to be very important, but, as expected for pyridine, the geometry relaxation contribution is essentially negligible. Using an average electrostatic model of the solvent, the magnetic shielding constant is calculated as -58.7 ppm, in good agreement with the experimental value of -56.3 ppm. The explicit inclusion of hydrogen-bonded water molecules embedded in the electrostatic field of the remaining solvent molecules gives the value of -61.8 ppm.

Atomistic Study of Water Confined in Silica.

In this work, we have used a combined of atomistic simulation methods to explore the effects of confinement of water molecules between silica surfaces. Firstly, the mechanical properties of water severe confined (~3A) between two silica alpha-quartz was determined based on first principles calculations within the density functional theory (DFT). Simulated annealing methods were employed due to the complex potential energry surface, and the difficulties to avoid local minima. Our results suggest that much of the stiffness of the material (46%) remains, even after the insertion of a water monolayer in the silica. Secondly, in order to access typical time scales for confined systems, classical molecular dynamics was used to determine the dynamical properties of water confined in silica cylindrical pores, with diameters varying from 10 to 40A. in this case we have varied the passivation of the silica surface, from 13% to 100% of SiOH, and the other terminations being SiOH2 and SiOH3, the distribution of the different terminations was obtained with a Monte Carlo simulation. The simulations indicates a lowering of the diffusion coefficientes as the diameter decreases, due to the structuration of hydrogen bonds of water molecules; we have also obtained the density profiles of the confined water and the interfacial tension.

Between the scales: Water from different perspectives

In der vorliegenden Arbeit werden verschiedene Wassermodelle in sogenannten Multiskalen-Computersimulationen mit zwei Auflösungen untersucht, in atomistischer Auflösung und in einer vergröberten Auflösung, die als "coarse-grained" bezeichnet wird. In der atomistischen Auflösung wird ein Wassermolekül, entsprechend seiner chemischen Struktur, durch drei Atome beschrieben, im Gegensatz dazu wird ein Molekül in der coarse-grained Auflösung durch eine Kugel dargestellt.rnrnDie coarse-grained Modelle, die in dieser Arbeit vorgestellt werden, werden mit verschiedenen coarse-graining Methoden entwickelt. Hierbei kommen hauptsächlich die "iterative Boltzmann Inversion" und die "iterative Monte Carlo Inversion" zum Einsatz. Beides sind struktur-basierte Ansätze, die darauf abzielen bestimmte strukturelle Eigenschaften, wie etwa die Paarverteilungsfunktionen, des zugrundeliegenden atomistischen Systems zu reproduzieren. Zur automatisierten Anwendung dieser Methoden wurde das Softwarepaket "Versatile Object-oriented Toolkit for Coarse-Graining Applications" (VOTCA) entwickelt.rnrnEs wird untersucht, in welchem Maße coarse-grained Modelle mehrere Eigenschaftenrndes zugrundeliegenden atomistischen Modells gleichzeitig reproduzieren können, z.B. thermodynamische Eigenschaften wie Druck und Kompressibilität oder strukturelle Eigenschaften, die nicht zur Modellbildung verwendet wurden, z.B. das tetraedrische Packungsverhalten, welches für viele spezielle Eigenschaft von Wasser verantwortlich ist.rnrnMit Hilfe des "Adaptive Resolution Schemes" werden beide Auflösungen in einer Simulation kombiniert. Dabei profitiert man von den Vorteilen beider Modelle:rnVon der detaillierten Darstellung eines räumlich kleinen Bereichs in atomistischer Auflösung und von der rechnerischen Effizienz des coarse-grained Modells, die den Bereich simulierbarer Zeit- und Längenskalen vergrössert.rnrnIn diesen Simulationen kann der Einfluss des Wasserstoffbrückenbindungsnetzwerks auf die Hydration von Fullerenen untersucht werden. Es zeigt sich, dass die Struktur der Wassermoleküle an der Oberfläche hauptsächlich von der Art der Wechselwirkung zwischen dem Fulleren und Wasser und weniger von dem Wasserstoffbrückenbindungsnetzwerk dominiert wird.rn

Elemetary processes of radiation damage in organic molecules of biological interest

It was observed in the ‘80s that the radiation damage on biological systems strongly depends on processes occurring at the microscopic level, involving the elementary constituents of biological cells. Since then, lot of attention has been paid to study elementary processes of photo- and ion-chemistry of isolated organic molecule of biological interest. This work fits in this framework and aims to study the radiation damage mechanisms induced by different types of radiations on simple halogenated biomolecules used as radiosensitizers in radiotherapy. The research is focused on the photofragmentation of halogenated pyrimidine molecules (5Br-pyrimidine, 2Br-pyrimidine and 2Cl-pyrimidine) in the VUV range and on the 12C4+ ion-impact fragmentation of the 5Br-uracil and its homogeneous and hydrated clusters. Although halogen substituted pyrimidines have similar structure to the pyrimidine molecule, their photodissociation dynamics is quite different. These targets have been chosen with the purpose of investigating the effect of the specific halogen atom and site of halogenation on the fragmentation dynamics. Theoretical and experimental studies have highlighted that the site of halogenation and the type of halogen atom, lead either to the preferential breaking of the pyrimidinic ring or to the release of halogen/hydrogen radicals. The two processes can subsequently trigger different mechanisms of biological damage. To understand the effect of the environment on the fragmentation dynamic of the single molecule, the ion-induced fragmentation of homogenous and hydrated clusters of 5Br-uracil have been studied and compared to similar studies on the isolated molecule. The results show that the “protective effect” of the environment on the single molecule hold in the homogeneous clusters, but not in the hydrated clusters, where several hydrated fragments have been observed. This indicates that the presence of water molecules can inhibit some fragmentation channels and promote the keto-enol tautomerization, which is very important in the mutagenesis of the DNA.

Accurate Predictions of Water Cluster Formation, (H2O)n=2-10

An efficient mixed molecular dynamics/quantum mechanics model has been applied to the water cluster system. The use of the MP2 method and correlation consistent basis sets, with appropriate correction for BSSE, allows for the accurate calculation of electronic and free energies for the formation of clusters of 2−10 water molecules. This approach reveals new low energy conformers for (H2O)n=7,9,10. The water heptamer conformers comprise five different structural motifs ranging from a three-dimensional prism to a quasi-planar book structure. A prism-like structure is favored energetically at low temperatures, but a chair-like structure is the global Gibbs free energy minimum past 200 K. The water nonamers exhibit less complexity with all the low energy structures shaped like a prism. The decamer has 30 conformers that are within 2 kcal/mol of the Gibbs free energy minimum structure at 298 K. These structures are categorized into four conformer classes, and a pentagonal prism is the most stable structure from 0 to 320 K. Results can be used as benchmark values for empirical water models and density functionals, and the method can be applied to larger water clusters.

Quantum-Mechanical Investigation of Large Water Clusters

The PM3 quantum-mechanical method has been used to study large water clusters ranging from 8 to 42 water molecules. These large clusters are built from smaller building blocks. The building blocks include cyclic tetramers, pentamers, octamers, and a pentagonal dodecahedron cage. The correlations between the strain energy resulting from bending of the hydrogen bonds formed by different cluster motifs and the number of waters involved in the cluster are discussed. The PM3 results are compared with TIP4P potential and ab initio results. The number of net hydrogen bonds per water increases with the cluster size. This places a limit on the size of clusters that would fit the Benson model of liquid water. Many of the 20-mer clusters fit the Benson model well. Calculations of the ion cluster (H20)4o(H30+)2 reveal that the m/e ratio obtainable by mass spectrometry experiments can uniquely indicate the conformation of the 20 water pentagonal dodecahedron cage present in the larger clusters.

Do Hydroxyl Radical−Water Clusters, OH(H2O)n, n = 1−5, Exist in the Atmosphere?

It has been speculated that the presence of OH(H2O)n clusters in the troposphere could have significant effects on the solar absorption balance and the reactivity of the hydroxyl radical. We have used the G3 and G3B3 model chemistries to model the structures and predict the frequencies of hydroxyl radical/water clusters containing one to five water molecules. The reaction between hydroxyl radical clusters and methane was examined as a function of water cluster size to gain an understanding of how cluster size affects the hydroxyl radical reactivity.

Prediction of Accurate Anharmonic Experimental Vibrational Frequencies for Water Clusters, (H2O)n, n = 2−5

Accurate anharmonic experimental vibrational frequencies for water clusters consisting of 2−5 water molecules have been predicted on the basis of comparing different methods with MP2/aug-cc-pVTZ calculated and experimental anharmonic frequencies. The combination of using HF/6-31G* scaled frequencies for intramolecular modes and anharmonic frequencies for intermolecular modes gives excellent agreement with experiment for the water dimer and trimer and are as good as the expensive anharmonic MP2 calculations. The water trimer, the cyclic Ci and S4 tetramers, and the cyclic pentamer all have unique peaks in the infrared spectrum between 500 and 800 cm-1 and between 3400 and 3700 cm-1. Under the right experimental conditions these different clusters can be uniquely identified using high-resolution IR spectroscopy.

Comparison of Model Chemistry and Density Functional Theory Thermochemical Predictions with Experiment for Formation of Ionic Clusters of the Ammonium Cation Complexed with Water and Ammonia; Atmospheric Implications

The G2, G3, CBS-QB3, and CBS-APNO model chemistry methods and the B3LYP, B3P86, mPW1PW, and PBE1PBE density functional theory (DFT) methods have been used to calculate ΔH° and ΔG° values for ionic clusters of the ammonium ion complexed with water and ammonia. Results for the clusters NH4+(NH3)n and NH4+(H2O)n, where n = 1−4, are reported in this paper and compared against experimental values. Agreement with the experimental values for ΔH° and ΔG° for formation of NH4+(NH3)n clusters is excellent. Comparison between experiment and theory for formation of the NH4+(H2O)n clusters is quite good considering the uncertainty in the experimental values. The four DFT methods yield excellent agreement with experiment and the model chemistry methods when the aug-cc-pVTZ basis set is used for energetic calculations and the 6-31G* basis set is used for geometries and frequencies. On the basis of these results, we predict that all ions in the lower troposphere will be saturated with at least one complete first hydration shell of water molecules.

Accurate Predictions of Water Cluster Formation, (H2O)n=2−10

An efficient mixed molecular dynamics/quantum mechanics model has been applied to the water cluster system. The use of the MP2 method and correlation consistent basis sets, with appropriate correction for BSSE, allows for the accurate calculation of electronic and free energies for the formation of clusters of 2−10 water molecules. This approach reveals new low energy conformers for (H2O)n=7,9,10. The water heptamer conformers comprise five different structural motifs ranging from a three-dimensional prism to a quasi-planar book structure. A prism-like structure is favored energetically at low temperatures, but a chair-like structure is the global Gibbs free energy minimum past 200 K. The water nonamers exhibit less complexity with all the low energy structures shaped like a prism. The decamer has 30 conformers that are within 2 kcal/mol of the Gibbs free energy minimum structure at 298 K. These structures are categorized into four conformer classes, and a pentagonal prism is the most stable structure from 0 to 320 K. Results can be used as benchmark values for empirical water models and density functionals, and the method can be applied to larger water clusters.