Azo-hydrazone tautomerism in protonated aminoazobenzenes: Resonance Raman spectroscopy and quantum-chemical calculations

The protonation effect on the vibrational and electronic spectra of 4-aminoazobenzene and 4-(dimethylamino)azobenzene was investigated by resonance Raman spectroscopy, and the results were discussed on the basis of quantum-chemical calculations. Although this class of molecular systems has been investigated in the past concerning the azo-hydrazone tautomerism, the present work is the first to use CASSCF/CASPT2 calculations to unveil the structure of both tautomers as well the nature of the molecular orbitals involved in chromophoric moieties responsible for the resonance Raman enhancement patterns. More specifically both the resonance Raman and theoretical results show clearly that in the neutral species, the charge transfer transition involves mainly the azo moiety, whereas in the protonated forms there is a great difference, depending on the tautomer. In fact, for the azo tautomer the transition is similar to that observed in the corresponding neutral species, whereas in the hydrazone tautomer such a transition is much more delocalized due to the contribution of the quinoid structure. The characterization of protonated species and the understanding of the tautomerization mechanism are crucial for controlling molecular properties depending on the polarity and pH of the medium.

Structural and optical properties of CaTiO3 perovskite-based materials obtained by microwave-assisted hydrothermal synthesis: An experimental and theoretical insight

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Spectroscopic studies of the Eu3+ and Er3+ ions in the fluorozirconate LaZr2F11 matrix

An investigation by optical spectroscopy of the Eu3 + and Er3 + active ions in the crystallized fluorozirconate matrix LaZr2F11 is presented. The D-5(1) --> F-7(0-5) emission lines of Eu3 + are used to extract the F-7(0-5) energy scheme and the observed extinctions permit the deduction of irreducible representations (IRREPS) associated with corresponding sub-levels in the D-2 symmetry. The crystal field analysis was carried out on a 387 x 387 basis set, comprising the F-7, D-5(1,2,3) F-5(1,2), (5)G(1,2,3) and P-3(1,2,3,4,5,6) terms of the Eu-3 (+) 4f(6) configuration. The deviation and rms are 6.8 and 7.9 cm (-1), respectively for 38 levels and ten parameters. The experimental crystal field parameters are in good agreement with the ab-initio ones. Moreover, the relative intensities of the D-5(0) --> F-7(2,3,4) emissions are well reproduced by an 'ab-initio' calculation, except for three lines. The Er3 + ions introduced in LaZr2F11, microcrystals also lie in an unique crystallographic site. A total of 31 energy levels were recorded and the crystal field analysis led to 6.6 and 7.8 cm (-1) for the deviation and rms, respectively, for nine variable parameters taken into account. The experimental CF parameters for Er3 + and Eu3 + are very similar, which seems to show that the host lattice contracts around the smaller Er3 + ion. The informations given by both Eu3 + and Eu3 + optical probes in LaZr2F11 are very consistent with the structure previously determined for the isotypic PrZr2F11 fluoride. (C) 2000 Elsevier B.V. B.V. All rights reserved.

Investigation on the structural properties in Er-doped PbTiO3 compounds: A correlation between experimental and theoretical results

The aim of this work is to investigate the structural properties of the (Pb1-xErx)TiO3 (PET) powders, with x varying from 0.01 to 0.08, prepared by the polymeric precursor method. Combined X-ray diffraction, Raman spectroscopy and ab initio calculation reveal a diffuse phase-transition of a tetragonal to a cubic phase. The crystalline models built allowed to calculate electronic properties and to analyze the behavior of the doping element in the structure of the material, which are consistent with the experimental results that indicate the beginning of phase-transition from tetragonal to cubic. (C) 2007 Elsevier B.V. All rights reserved.

Determination of the Θ+ parity from γn → K -K+n

It is demonstrated that measurements of photon asymmetry in the γn → K-K+n reaction, can most likely determine the parity of the newly discovered Θ+ pentaquark. We predict that if the parity of Θ+ is positive, the photon asymmetry is significantly positive; if the parity is negative, the photon asymmetry is significantly negative. If the background contribution is large, the photon asymmetry may become very small in magnitude, thereby making it difficult to distinguish between the positive and negative parity results. However, even in this case, a combined analysis of the (K+n) invariant mass distribution and photon asymmetry should allow a determination of the parity of Θ+. © 2004 Published by Elsevier B.V.

Towards an insight on photodamage in hair fibre by UV-light: An experimental and theoretical study

Synopsis: Objectives: In this research, an experimental and theoretical study was conducted to design a photodegradation mechanism of the amino acid tryptophan (Trp) in hair fibres. Methods: For the experimental research, Caucasian hair fibres were exposed to several different solar radiation simulation periods. Then, Trp and its photoproducts (N-formylkynurenine and kynurenine) were assayed by excitation and emission spectroscopic analysis. Results: For the theoretical study, reactions involved in the photodegradation of Trp were evaluated by high-level quantum mechanical calculations in a density functional theory (DFT) framework which indicate a probable Trp degradation mechanism with a minimum expended energy pathway. Conclusion: The biochemistry concerning these reactions is essentially important for a biological system where the degradation of Trp occurs. © 2013 John Wiley & Sons Ltd.

Propriedades ferroelétricas, microestruturais e ópticas dos materiais cerâmicos Ba0,5 Sr0,5 (Ti1-ySny) O3

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

MOLECULAR DYNAMICS STUDIES OF SIMPLE MODEL FLUIDS AND WATER CONFINED IN CARBON NANOTUBE

Molecular Dynamics (MD) simulation is one of the most important computational techniques with broad applications in physics, chemistry, chemical engineering, materials design and biological science. Traditional computational chemistry refers to quantum calculations based on solving Schrodinger equations. Later developed Density Functional Theory (DFT) based on solving Kohn-Sham equations became the more popular ab initio calculation technique which could deal with ~1000 atoms by explicitly considering electron interactions. In contrast, MD simulation based on solving classical mechanics equations of motion is a totally different technique in the field of computational chemistry. Electron interactions were implicitly included in the empirical atom-based potential functions and the system size to be investigated can be extended to ~106 atoms. The thermodynamic properties of model fluids are mainly determined by macroscopic quantities, like temperature, pressure, density. The quantum effects on thermodynamic properties like melting point, surface tension are not dominant. In this work, we mainly investigated the melting point, surface tension (liquid-vapor and liquid-solid) of model fluids including Lennard-Jones model, Stockmayer model and a couple of water models (TIP4P/Ew, TIP5P/Ew) by means of MD simulation. In addition, some new structures of water confined in carbon nanotube were discovered and transport behaviors of water and ions through nano-channels were also revealed.

Size-selective investigations of spectral and emission properties of small particles and nanotubes

In dieser Arbeit wurde die Elektronenemission von Nanopartikeln auf Oberflächen mittels spektroskopischen Photoelektronenmikroskopie untersucht. Speziell wurden metallische Nanocluster untersucht, als selbstorganisierte Ensembles auf Silizium oder Glassubstraten, sowie ferner ein Metall-Chalcogenid (MoS2) Nanoröhren-Prototyp auf Silizium. Der Hauptteil der Untersuchungen war auf die Wechselwirkung von fs-Laserstrahlung mit den Nanopartikeln konzentriert. Die Energie der Lichtquanten war kleiner als die Austrittsarbeit der untersuchten Proben, so dass Ein-Photonen-Photoemission ausgeschlossen werden konnte. Unsere Untersuchungen zeigten, dass ausgehend von einem kontinuierlichen Metallfilm bis hin zu Clusterfilmen ein anderer Emissionsmechanismus konkurrierend zur Multiphotonen-Photoemission auftritt und für kleine Cluster zu dominieren beginnt. Die Natur dieses neuen Mechanismus` wurde durch verschiedenartige Experimente untersucht. Der Übergang von einem kontinuierlichen zu einem Nanopartikelfilm ist begleitet von einer Zunahme des Emissionsstroms von mehr als eine Größenordnung. Die Photoemissions-Intensität wächst mit abnehmender zeitlicher Breite des Laserpulses, aber diese Abhängigkeit wird weniger steil mit sinkender Partikelgröße. Die experimentellen Resultate wurden durch verschiedene Elektronenemissions-Mechanismen erklärt, z.B. Multiphotonen-Photoemission (nPPE), thermionische Emission und thermisch unterstützte nPPE sowie optische Feldemission. Der erste Mechanismus überwiegt für kontinuierliche Filme und Partikel mit Größen oberhalb von mehreren zehn Nanometern, der zweite und dritte für Filme von Nanopartikeln von einer Größe von wenigen Nanometern. Die mikrospektroskopischen Messungen bestätigten den 2PPE-Emissionsmechanismus von dünnen Silberfilmen bei „blauer“ Laseranregung (hν=375-425nm). Das Einsetzen des Ferminiveaus ist relativ scharf und verschiebt sich um 2hν, wenn die Quantenenergie erhöht wird, wogegen es bei „roter“ Laseranregung (hν=750-850nm) deutlich verbreitert ist. Es zeigte sich, dass mit zunehmender Laserleistung die Ausbeute von niederenergetischen Elektronen schwächer zunimmt als die Ausbeute von höherenergetischen Elektronen nahe der Fermikante in einem Spektrum. Das ist ein klarer Hinweis auf eine Koexistenz verschiedener Emissionsmechanismen in einem Spektrum. Um die Größenabhängigkeit des Emissionsverhaltens theoretisch zu verstehen, wurde ein statistischer Zugang zur Lichtabsorption kleiner Metallpartikel abgeleitet und diskutiert. Die Elektronenemissionseigenschaften bei Laseranregung wurden in zusätzlichen Untersuchungen mit einer anderen Anregungsart verglichen, der Passage eines Tunnelstroms durch einen Metall-Clusterfilm nahe der Perkolationsschwelle. Die elektrischen und Emissionseigenschaften von stromtragenden Silberclusterfilmen, welche in einer schmalen Lücke (5-25 µm Breite) zwischen Silberkontakten auf einem Isolator hergestellt wurden, wurden zum ersten Mal mit einem Emissions-Elektronenmikroskop (EEM) untersucht. Die Elektronenemission beginnt im nicht-Ohmschen Bereich der Leitungsstrom-Spannungskurve des Clusterfilms. Wir untersuchten das Verhalten eines einzigen Emissionszentrums im EEM. Es zeigte sich, dass die Emissionszentren in einem stromleitenden Silberclusterfilm Punktquellen für Elektronen sind, welche hohe Emissions-Stromdichten (mehr als 100 A/cm2) tragen können. Die Breite der Energieverteilung der Elektronen von einem einzelnen Emissionszentrum wurde auf etwa 0.5-0.6 eV abgeschätzt. Als Emissionsmechanismus wird die thermionische Emission von dem „steady-state“ heißen Elektronengas in stromdurchflossenen metallischen Partikeln vorgeschlagen. Größenselektierte, einzelne auf Si-Substraten deponierte MoS2-Nanoröhren wurden mit einer Flugzeit-basierten Zweiphotonen-Photoemissions-Spektromikroskopie untersucht. Die Nanoröhren-Spektren wiesen bei fs-Laser Anregung eine erstaunlich hohe Emissionsintensität auf, deutlich höher als die SiOx Substratoberfläche. Dagegen waren die Röhren unsichtbar bei VUV-Anregung bei hν=21.2 eV. Eine ab-initio-Rechnung für einen MoS2-Slab erklärt die hohe Intensität durch eine hohe Dichte freier intermediärer Zustände beim Zweiphotonen-Übergang bei hν=3.1 eV.

Investigation of heterogeneous catalysts for the synthesis of fine chemicals

In order to match the more stringent environmental regulations, heterogenization of traditional homogeneous processes is one of the main challenges of the modern chemical industry. Great results have been achieved in the fields of petrochemicals and base chemicals, whereas in fine chemical industry most of the synthetic procedures are based on multistep processes catalyzed by homogeneous catalysts mainly used in stoichiometric amounts. In the fine chemicals manufacture not so much efforts have been devoted to the investigation of suitable solid catalysts for the development of greener processes, then this sector represent a very attractive field of research. In this context, the present work deals with the extensive investigation of the possibility to heterogenize existing processes, in particular two different classes of reactions have been studied: alkylation of aromatic and heteroaromatic compounds and selective oxidation of aromatic alcohols. Traditional solid acid catalysts, such as zeolites, clays and alumina have been tested in the gas phase alkylation of 1,2-methylendioxybenzene, core building block of many drugs, pesticides and fragrances. The observed reactivity were clarified through a deep FTIR investigation complemented by ab initio calculation. The same catalysts were tested in the gas phase isopropylation of thiophene with the aim of clearly attribute the role of the reaction parameters in the reaction proceeding and verify the possibility to enhance the selectivity of one of the two possible isomers. Finally various Au/CeO2 catalysts were tested in the synthesis of benzaldehyde and piperonal, two aldehydes largely employed in the manufacture of fine chemical products, through liquid phase oxidation of the corresponding alcohols in very mild conditions.

Hadronic matrix elements in lattice QCD

The lattice formulation of Quantum ChromoDynamics (QCD) has become a reliable tool providing an ab initio calculation of low-energy quantities. Despite numerous successes, systematic uncertainties, such as discretisation effects, finite-size effects, and contaminations from excited states, are inherent in any lattice calculation. Simulations with controlled systematic uncertainties and close to the physical pion mass have become state-of-the-art. We present such a calculation for various hadronic matrix elements using non-perturbatively O(a)-improved Wilson fermions with two dynamical light quark flavours. The main topics covered in this thesis are the axial charge of the nucleon, the electro-magnetic form factors of the nucleon, and the leading hadronic contributions to the anomalous magnetic moment of the muon. Lattice simulations typically tend to underestimate the axial charge of the nucleon by 5 − 10%. We show that including excited state contaminations using the summed operator insertion method leads to agreement with the experimentally determined value. Further studies of systematic uncertainties reveal only small discretisation effects. For the electro-magnetic form factors of the nucleon, we see a similar contamination from excited states as for the axial charge. The electro-magnetic radii, extracted from a dipole fit to the momentum dependence of the form factors, show no indication of finite-size or cutoff effects. If we include excited states using the summed operator insertion method, we achieve better agreement with the radii from phenomenology. The anomalous magnetic moment of the muon can be measured and predicted to very high precision. The theoretical prediction of the anomalous magnetic moment receives contribution from strong, weak, and electro-magnetic interactions, where the hadronic contributions dominate the uncertainties. A persistent 3σ tension between the experimental determination and the theoretical calculation is found, which is considered to be an indication for physics beyond the Standard Model. We present a calculation of the connected part of the hadronic vacuum polarisation using lattice QCD. Partially twisted boundary conditions lead to a significant improvement of the vacuum polarisation in the region of small momentum transfer, which is crucial in the extraction of the hadronic vacuum polarisation.

Photodissociation dynamics of tert-butylnitrite following excitation to the S(1) and S(2) states. A study by velocity-map ion-imaging and 3D-REMPI spectroscopy

Excitation of tert-butylnitrite into the first and second UV absorption bands leads to efficient dissociation into the fragment radicals NO and tert-butoxy in their electronic ground states (2)Π and (2)E, respectively. Velocity distributions and angular anisotropies for the NO fragment in several hundred rotational and vibrational quantum states were obtained by velocity-map imaging and the recently developed 3D-REMPI method. Excitation into the well resolved vibronic progression bands (k = 0, 1, 2) of the NO stretch mode in the S(1) ← S(0) transition produces NO fragments mostly in the vibrational state with v = k, with smaller fractions in v = k - 1 and v = k - 2. It is concluded that dissociation occurs on the purely repulsive PES of S(1) without barrier. All velocity distributions from photolysis via the S(1)(nπ*) state are monomodal and show high negative anisotropy (β ≈ -1). The rotational distributions peak near j = 30.5 irrespective of the vibronic state S(1)(k) excited and the vibrational state v of the NO fragment. On average 46% of the excess energy is converted to kinetic energy, 23% and 31% remain as internal energy in the NO fragment and the t-BuO radical, respectively. Photolysis via excitation into the S(2) ← S(0) transition at 227 nm yields NO fragments with about equal populations in v = 0 and v = 1. The rotational distributions have a single maximum near j = 59.5. The velocity distributions are monomodal with positive anisotropy β ≈ 0.8. The average fractions of the excess energy distributed into translation, internal energy of NO, and internal energy of t-BuO are 39%, 23%, and 38%, respectively. In all cases ∼8500 cm(-1) of energy remain in the internal degrees of freedom of the t-BuO fragment. This is mostly assigned to rotational energy. An ab initio calculation of the dynamic reaction path shows that not only the NO fragment but also the t-BuO fragment gain large angular momentum during dissociation on the purely repulsive potential energy surface of S(2).

Peptide nucleic acids: Mechanism of tight binding and application in artificial nuclease

Peptide nucleic acids (PNA) are mimics of nucleic acids with a peptidic backbone. Duplexes and triplexes formed between PNA and DNA or RNA possess remarkable thermal stability, they are resistant to nuclease cleavage and can better discriminate mismatches. Understanding the mechanism for the tight binding between PNA and oligonucleotides is important for the design and development of better PNA-based drugs.^ We have performed molecular dynamics (MD) simulations of 8-mer PNA/DNA duplex and two analogous duplexes with chiral modification of PNA strand (D- or L-Alanine modification). MD simulations were performed with explicit water and Na$\sp{+}$ counter ions. The 1.5-ns simulations were carried out with AMBER using periodic boundary and particle mesh Ewald summation. The point charges for PNA monomers were derived from fitting electrostatic potentials, obtained from ab initio calculation, to atomic centers using RESP. Derived charges reveal significantly altered charge distribution on the PNA bases and predict the Watson-Crick H-bonds involving PNA to be stronger. Results from NMR studies investigating H-bond interactions between DNA-DNA and DNA-PNA base pairs in non-polar environment are consistent with this prediction. MD simulations demonstrated that the PNA strand is more flexible than the DNA strand in the same duplex. That this flexibility might be important for the duplex stability is tested by introducing modification into the PNA backbones. Results from MD simulation revealed dramatically altered structures for the modified PNA-DNA duplexes. Consistent with previous NMR results, we also found no intrachain hydrogen bonds between O7$\sp\prime$ and N1$\sp\prime$ of the neighboring residues in our MD study. Our study reveals that in addition to the lack of charge repulsion, stronger Watson-Crick hydrogen bonds together with flexible backbone are important factors for the enhanced stability of the PNA-DNA duplex.^ In a related study, we have developed an application of Gly-Gly-His-(Gly)$\sb3$-PNA conjugate as an artificial nuclease. We were able to demonstrate cleavage of single stranded DNA at a single site upon Ni(II) binding to Gly-Gly-His tripeptide and activation of nuclease with monoperoxyphthalic acid. ^

Vibrational dynamics and structure of natural feldspars

The present study gives a contribution to the knowledge on the Na-feldspar and plagioclases, extending the database of the Raman spectra of plagioclases with different chemical compositions and structural orders. This information may be used for the future planetary explorations by “rovers”, for the investigation of ceramics nanocrystal materials and for the mineralogical phase identification in sediments. Na-feldspar and plagioclase solid solution have been investigated by Raman spectroscopy in order to determine the relationships between the vibrational changes and the plagioclase crystal chemistry and structure. We focused on the Raman micro-spectroscopy technique, being a non-destructive method, suited for contactless analysis with high spatial resolution. Chemical and structural analyses have been performed on natural samples to test the usefulness of Raman spectroscopy as a tool in the study of the pressure-induced structural deformations, the disordering processes due to change in the Al-Si distribution in the tetrahedral sites and, finally, in the determination of the anorthitic content (Anx) in plagioclase minerals. All the predicted 39 Ag Raman active modes have been identified and assigned to specific patterns of atomic vibrational motion. A detailed comparison between experimental and computed Raman spectra has been performed and previous assignments have been revised, solving some discrepancies reported in recent literature. The ab initio calculation at the hybrid HF/DFT level with the WC1LYP Hamiltonian has proven to give excellent agreement between calculated and experimentally measured Raman wavenumbers and intensities in triclinic minerals. A short digression on the 36 infrared active modes of Na-feldspar has been done too. The identification of all 39 computed Raman modes in the experimentally measured spectra of the fully ordered Na-feldspar, known as low albite, along with the detailed description of each vibrational mode, has been essential to extend the comparative analysis to the high pressure and high temperature structural forms of albite, which reflect the physical–chemical conditions of the hosting rocks. The understanding of feldspar structure response to pressure and temperature is crucial in order to constrain crustal behaviour. The compressional behaviour of the Na-feldspar has been investigated for the first time by Raman spectroscopy. The absence of phase transitions and the occurrence of two secondary compression mechanisms acting at different pressures have been confirmed. Moreover, Raman data suggest that the internal structural changes are confined to a small pressure interval, localized around 6 GPa, not spread out from 4 to 8 GPa as suggested by previous X-rays studies on elasticity. The dominant compression mechanisms act via tetrahedral tilting, while the T-O bond lengths remain nearly constant at moderate compressional regimes. At the spectroscopic level, this leads to the strong pressure dependencies of T-O-T bending modes, as found for the four modes at 478, 508, 578 and 815 cm-1. The Al-Si distribution in the tetrahedral sites affects also the Raman spectrum of Na-feldspar. In particular, peak broadening is more sensitive than peak position to changes in the degree of order. Raman spectroscopy is found to be a good probe for local ordering, in particular being sensitive to the first annealing steps, when the macroscopic order parameter is still high. Even though Raman data are scattered and there are outliers in the estimated values of the degree of order, the average peak linewidths of the Na-feldspar characteristic doublet band, labelled here as υa and υb, as a function of the order parameter Qod show interesting trends: both peak linewidths linearly increase until saturation. From Qod values lower than 0.6, peak broadening is no more affected by the Al-Si distribution. Moreover, the disordering process is found to be heterogeneous. SC-XRD and Raman data have suggested an inter-crystalline inhomogeneity of the samples, i.e., the presence of regions with different defect density on the micrometric scale. Finally, the influence of Ca-Na substitution in the plagioclase Raman spectra has been investigated. Raman spectra have been collected on a series of well characterized natural, low structural plagioclases. The variations of the Raman modes as a function of the chemical composition and the structural order have been determined. The number of the observed Raman bands at each composition gives information about the unit-cell symmetry: moving away from the C1 structures, the number of the Raman bands enhances, as the number of formula units in the unit cell increases. The modification from an “albite-like” Raman spectrum to a more “anorthite-like” spectrum occurs from sample An78 onwards, which coincides with the appearance of c reflections in the diffraction patterns of the samples. The evolution of the Raman bands υa and υb displays two changes in slope at ~An45 and ~An75: the first one occurs between e2 and e1 plagioclases, the latter separates e1 and I1 plagioclases with only b reflections in their diffraction patterns from I1 and P1 samples having b and c reflections too. The first variation represents exactly the e2→e1 phase transitions, whereas the second one corresponds in good approximation to the C1→I1 transition, which has been determined at ~An70 by previous works. The I1→P1 phase transition in the anorthite-rich side of the solid solution is not highlighted in the collected Raman spectra. Variations in peak broadening provide insights into the behaviour of the order parameter on a local scale, suggesting an increase in the structural disorder within the solid solution, as the structures have to incorporate more Al atoms to balance the change from monovalent to divalent cations. All the information acquired on these natural plagioclases has been used to produce a protocol able to give a preliminary estimation of the chemical composition of an unknown plagioclase from its Raman spectrum. Two calibration curves, one for albite-rich plagioclases and the other one for the anorthite-rich plagioclases, have been proposed by relating the peak linewidth of the most intense Raman band υa and the An content. It has been pointed out that the dependence of the composition from the linewidth can be obtained only for low structural plagioclases with a degree of order not far away from the references. The proposed tool has been tested on three mineralogical samples, two of meteoric origin and one of volcanic origin. Chemical compositions by Raman spectroscopy compare well, within an error of about 10%, with those obtained by elemental techniques. Further analyses on plagioclases with unknown composition will be necessary to validate the suggested method and introduce it as routine tool for the determination of the chemical composition from Raman data in planetary missions.

Applying many-body expansion aiming to make ab in initio molecular dynamics more efficient

In this present work we present a methodology that aims to apply the many-body expansion to decrease the computational cost of ab initio molecular dynamics, keeping acceptable accuracy on the results. We implemented this methodology in a program which we called ManBo. In the many-body expansion approach, we partitioned the total energy E of the system in contributions of one body, two bodies, three bodies, etc., until the contribution of the Nth body [1-3]: E = E1 + E2 + E3 + …EN. The E1 term is the sum of the internal energy of the molecules; the term E2 is the energy due to interaction between all pairs of molecules; E3 is the energy due to interaction between all trios of molecules; and so on. In Manbo we chose to truncate the expansion in the contribution of two or three bodies, both for the calculation of the energy and for the calculation of the atomic forces. In order to partially include the many-body interactions neglected when we truncate the expansion, we can include an electrostatic embedding in the electronic structure calculations, instead of considering the monomers, pairs and trios as isolated molecules in space. In simulations we made we chose to simulate water molecules, and use the Gaussian 09 as external program to calculate the atomic forces and energy of the system, as well as reference program for analyzing the accuracy of the results obtained with the ManBo. The results show that the use of the many-body expansion seems to be an interesting approach for reducing the still prohibitive computational cost of ab initio molecular dynamics. The errors introduced on atomic forces in applying such methodology are very small. The inclusion of an embedding electrostatic seems to be a good solution for improving the results with only a small increase in simulation time. As we increase the level of calculation, the simulation time of ManBo tends to largely decrease in relation to a conventional BOMD simulation of Gaussian, due to better scalability of the methodology presented. References [1] E. E. Dahlke and D. G. Truhlar; J. Chem. Theory Comput., 3, 46 (2007). [2] E. E. Dahlke and D. G. Truhlar; J. Chem. Theory Comput., 4, 1 (2008). [3] R. Rivelino, P. Chaudhuri and S. Canuto; J. Chem. Phys., 118, 10593 (2003).