Molecular dynamics simulations of silicate and borate glasses and melts : structure, diffusion dynamics and vibrational properties

Molecular dynamics simulations of silicate and borate glasses and melts: Structure, diffusion dynamics and vibrational properties. In this work computer simulations of the model glass formers SiO2 and B2O3 are presented, using the techniques of classical molecular dynamics (MD) simulations and quantum mechanical calculations, based on density functional theory (DFT). The latter limits the system size to about 100−200 atoms. SiO2 and B2O3 are the two most important network formers for industrial applications of oxide glasses. Glass samples are generated by means of a quench from the melt with classical MD simulations and a subsequent structural relaxation with DFT forces. In addition, full ab initio quenches are carried out with a significantly faster cooling rate. In principle, the structural properties are in good agreement with experimental results from neutron and X-ray scattering, in all cases. A special focus is on the study of vibrational properties, as they give access to low-temperature thermodynamic properties. The vibrational spectra are calculated by the so-called ”frozen phonon” method. In all cases, the DFT curves show an acceptable agreement with experimental results of inelastic neutron scattering. In case of the model glass former B2O3, a new classical interaction potential is parametrized, based on the liquid trajectory of an ab initio MD simulation at 2300 K. In this course, a structural fitting routine is used. The inclusion of 3-body angular interactions leads to a significantly improved agreement of the liquid properties of the classical MD and ab initio MD simulations. However, the generated glass structures, in all cases, show a significantly lower fraction of 3-membered planar boroxol rings as predicted by experimental results (f=60%-80%). The largest boroxol ring fraction of f=15±5% is observed in the full ab initio quenches from 2300 K. In case of SiO2, the glass structures after the quantum mechanical relaxation are the basis for calculations of the linear thermal expansion coefficient αL(T), employing the quasi-harmonic approximation. The striking observation is a change change of sign of αL(T) going along with a temperature range of negative αL(T) at low temperatures, which is in good agreement with experimental results.

Looking for silicene: studies of silicon deposition on metallic and semiconductor substrates

Nel presente lavoro espongo i risultati degli esperimenti svolti durante la mia internship all’Institut des NanoSciences de Paris (INSP), presso l’Università Pierre et Marie Curie (Paris VI), nel team "Phisico-Chimie et Dynamique des Surfaces", sotto la supervisione del Dott. Geoffroy Prévot. L’elaborato è stato redatto e in- tegrato sotto la guida del Dott. Pasquini, del dipartimento di Fisica e Astronomia dell’Università di Bologna. La tesi s’inserisce nel campo di ricerca del silicene, i.e. l’allotropo bidimensionale del silicio. Il cosidetto free-standing silicene è stato predetto teoricamente nel 2009 utilizzando calcoli di Density Functional Theory, e da allora ha stimolato un’intensa ricerca per la sua realizzazione sperimentale. La sua struttura elettronica lo rende particolarmente adatto per eventuali appli- cazioni tecnologiche e sperimentali, mentre lo studio delle sue proprietà è di grande interesse per la scienza di base. Nel capitolo 1 presento innanzitutto la struttura del silicene e le proprietà previste dagli studi pubblicati nella letteratura scientifica. In seguito espongo alcuni dei risultati sperimentali ottenuti negli ultimi anni, in quanto utili per un paragone con i risultati ottenuti durante l’internship. Nel capitolo 2 presento le tecniche sperimentali che ho utilizzato per effettuare le misure. Molto tempo è stato investito per ottenere una certa dimistichezza con gli apparati in modo da svolgere gli esperimenti in maniera autonoma. Il capitolo 3 è dedicato alla discussione e analisi dei risultati delle misure, che sono presentati in relazione ad alcune considerazioni esposte nel primo capitolo. Infine le conclusioni riassumono brevemente quanto ottenuto dall’analisi dati. A partire da queste considerazioni propongo alcuni esperimenti che potrebbero ulteriormente contribuire alla ricerca del silicene. I risultati ottenuti su Ag(111) sono contenuti in un articolo accettato da Physical Review B.

Further Quantum Mechanical Evidence that Difluorotoluene does not Hydrogen Bond

Calculations were run on the methylated DNA base pairs adenine:thymine and adenine:difluorotoluene to further investigate the hydrogen-bonding properties of difluorotoluene (F). Geometries were optimized using hybrid density functional theory. Single-point calculations at the MP2(full) level were performed to obtain more rigorous energies. The functional counterpoise method was used to correct for the basis set superposition error (BSSE), and the interaction energies were also corrected for fragment relaxation. These corrections brought the B3LYP and MP2 interaction energies into excellent agreement. In the gas phase, the Gibbs free energies calculated at the B3LYP and MP2 levels of theory predict that A and T will spontaneously form an A:T pair while A:F spontaneously dissociates into A and F. Solvation effects on the pairing of the bases were explored using implicit solvent models for water and chloroform. In aqueous solution, both A:T and A:F are predicted to dissociate into their component monomers. Semiempirical calculations were performed on small sections of B-form DNA containing the two pairs, and the results provide support for the concept that base stacking is more important than hydrogen bonding for the stability of the A:F pair within a DNA helix.

Efficient and Accurate Characterization of the Bergman Cyclization for an Expanded Substructure of Esperamicin A1

Incorporation of enediynes into anticancer drugs remains an intriguing yet elusive strategy for the design of therapeutically active agents. Density functional theory was used to locate reactants, products, and transition states along the Bergman cyclization pathways connecting enediynes to reactive para-biradicals. Sum method correction to low-level calculations confirmed B3LYP/6-31G(d,p) as the method of choice in investigating enediynes. Herein described as MI:Sum, calculated reaction enthalpies differed from experiment by an average of 2.1 kcal·mol−1 (mean unsigned error). A combination of strain energy released across the reaction coordinate and the critical intramolecular distance between reacting diynes explains reactivity differences. Where experimental and calculated barrier heights are in disagreement, higher level multireference treatment of the enediynes confirms lower level estimates. Previous work concerning the chemically reactive fragment of esperamcin, MTC, is expanded to our model system MTC2.

Quantum Mechanical Study of Sulfuric Acid Hydration: Atmospheric Implications

The role of the binary nucleation of sulfuric acid in aerosol formation and its implications for global warming is one of the fundamental unsettled questions in atmospheric chemistry. We have investigated the thermodynamics of sulfuric acid hydration using ab initio quantum mechanical methods. For H2SO4(H2O)n where n = 1–6, we used a scheme combining molecular dynamics configurational sampling with high-level ab initio calculations to locate the global and many low lying local minima for each cluster size. For each isomer, we extrapolated the Møller–Plesset perturbation theory (MP2) energies to their complete basis set (CBS) limit and added finite temperature corrections within the rigid-rotor-harmonic-oscillator (RRHO) model using scaled harmonic vibrational frequencies. We found that ionic pair (HSO4–·H3O+)(H2O)n−1clusters are competitive with the neutral (H2SO4)(H2O)n clusters for n ≥ 3 and are more stable than neutral clusters for n ≥ 4 depending on the temperature. The Boltzmann averaged Gibbs free energies for the formation of H2SO4(H2O)n clusters are favorable in colder regions of the troposphere (T = 216.65–273.15 K) for n = 1–6, but the formation of clusters with n ≥ 5 is not favorable at higher (T > 273.15 K) temperatures. Our results suggest the critical cluster of a binary H2SO4–H2O system must contain more than one H2SO4 and are in concert with recent findings(1) that the role of binary nucleation is small at ambient conditions, but significant at colder regions of the troposphere. Overall, the results support the idea that binary nucleation of sulfuric acid and water cannot account for nucleation of sulfuric acid in the lower troposphere.

Au@Hg nanoalloy formation through direct amalgamation: Structural, spectroscopic, and computational evidence for slow nanoscale diffusion

Dynamic core-shell nanoparticles have received increasing attention in recent years. This paper presents a detailed study of Au-Hg nanoalloys, whose composing elements show a large difference in cohesive energy. A simple method to prepare Au@Hg particles with precise control over the composition up to 15 atom% mercury is introduced, based on reacting a citrate stabilized gold sol with elemental mercury. Transmission electron microscopy shows an increase of particle size with increasing mercury content and, together with X-ray powder diffraction, points towards the presence of a core-shell structure with a gold core surrounded by an Au-Hg solid solution layer. The amalgamation process is described by pseudo-zero-order reaction kinetics, which indicates slow dissolution of mercury in water as the rate determining step, followed by fast scavenging by nanoparticles in solution. Once adsorbed at the surface, slow diffusion of Hg into the particle lattice occurs, to a depth of ca. 3 nm, independent of Hg concentration. Discrete dipole approximation calculations relate the UV-vis spectra to the microscopic details of the nanoalloy structure. Segregation energies and metal distribution in the nanoalloys were modeled by density functional theory calculations. The results indicate slow metal interdiffusion at the nanoscale, which has important implications for synthetic methods aimed at core-shell particles.

Conformationally Controlled Electron Delocalization in n-Type Rods: Synthesis, Structure, and Optical, Electrochemical, and Spectroelectrochemical Properties of Dicyanocyclophanes

A series of dicyanobiphenyl-cyclophanes 1-6 with various pi-backbone conformations and characteristic n-type semiconductor properties is presented. Their synthesis, optical, structural, electrochemical, spectroelectrochemical, and packing properties are investigated. The X-ray crystal structures of all n-type rods allow the systematic correlation of structural features with physical properties. In addition, the results are supported by quantum mechanical calculations based on density functional theory. A two-step reduction process is observed for all n-type rods, in which the first step is reversible. The potential gap between the reduction processes depends linearly on the cos(2) value of the torsion angle phi between the pi-systems. Similarly, optical absorption spectroscopy shows that the vertical excitation energy of the conjugation band correlates with the cos(2) value of the torsion angle phi. These correlations demonstrate that the fixed intramolecular torsion angle phi is the dominant factor determining the extent of electron delocalization in these model compounds, and that the angle phi measured in the solid-state structure is a good proxy for the molecular conformation in solution. Spectroelectrochemical investigations demonstrate that conformational rigidity is maintained even in the radical anion form. In particular, the absorption bands corresponding to the SOMO-LUMO+i transitions are shifted bathochromically, whereas the absorption bands corresponding to the HOMO-SOMO transition are shifted hypsochromically with increasing torsion angle phi.

Intermolecular Clamping by Hydrogen Bonds: 2-Pyridone center dot NH3

A combined spectroscopic and ab initio theoretical study of the doubly hydrogen-bonded complex of 2-pyridone (2PY) with NH3 has been performed. The S-1 <- S-0 spectrum extends up to approximate to 1200 cm(-1) above the 0(0)(0) band, close to twice the range observed for 2PY. The S-1 state nonradiative decay for vibrations above approximate to 300 cm(-1) in the NH3 complex is dramatically slowed down relative to bare 2PY. Also, the Delta v=2,4,... overtone bands of the v(1)' and v(2)' out-of-plane vibrations that dominate the low-energy spectral region of 2PY are much weaker or missing for 2PY center dot NH3, which implies that the bridging (2PY)NH center dot center dot center dot NH3 and H2NH center dot center dot center dot O=C H-bonds clamp the 2PY at a planar geometry in the S-1 state. The mass-resolved UV vibronic spectra of jet-cooled 2PY center dot NH3 and its H/D mixed isotopomers are measured using two-color resonant two-photon ionization spectroscopy. The S-0 and S-1 equilibrium structures and normal-mode frequencies are calculated by density functional (B3LYP) and correlated ab initio methods (MP2 and approximate second-order coupled-cluster, CC2). The S-1 <- S-0 vibronic assignments are based on configuration interaction singles (CIS) and CC2 calculations. A doubly H-bonded bridged structure of C-S symmetry is predicted, in agreement with that of Held and Pratt [J. Am. Chem. Soc. 1993, 115, 9718]. While the B3LYP and MP2 calculated rotational constants are in very good agreement with experiment, the calculated H2NH center dot center dot center dot O=C H-bond distance is approximate to 0.7 angstrom shorter than that derived by Held and Pratt. On the other hand, this underlines their observation that ammonia can act as a strong H-bond donor when built into an H-bonded bridge. The CC2 calculations predict the H2NH center dot center dot center dot O distance to increase by 0.2 angstrom upon S-1 <- S-0 electronic excitation, while the (2PY)NH center dot center dot center dot NH3 H-bond remains nearly unchanged. Thus, the expansion of the doubly H-bonded bridge in the excited state is asymmetric and almost wholly due to the weakening of the interaction of ammonia with the keto acceptor group.

N-H center dot center dot center dot pi hydrogen-bonding and large-amplitude tipping vibrations in jet-cooled pyrrole-benzene

The N-H center dot center dot center dot pi hydrogen bond is an important intermolecular interaction in many biological systems. We have investigated the infrared (IR) and ultraviolet (UV) spectra of the supersonic-jet cooled complex of pyrrole with benzene and benzene-d(6) (Pyr center dot Bz, Pyr center dot Bz-d(6)). DFT-D density functional, SCS-MP2 and SCS-CC2 calculations predict a T-shaped and (almost) C(s) symmetric structure with an N-H center dot center dot center dot pi hydrogen bond to the benzene ring. The pyrrole is tipped by omega(S(0)) = +/- 13 degrees relative to the surface normal of Bz. The N center dot center dot center dot ring distance is 3.13 angstrom. In the S(1) excited state, SCS-CC2 calculations predict an increased tipping angle omega(S(1)) = +/- 21 degrees. The IR depletion spectra support the T-shaped geometry: The NH stretch is redshifted by -59 cm(-1), relative to the "free" NH stretch of pyrrole at 3531 cm(-1), indicating a moderately strong N-H center dot center dot center dot pi interaction. The interaction is weaker than in the (Pyr)(2) dimer, where the NH donor shift is -87 cm(-1) [Dauster et al., Phys. Chem. Chem. Phys., 2008, 10, 2827]. The IR C-H stretch frequencies and intensities of the Bz subunit are very similar to those of the acceptor in the (Bz)(2) dimer, confirming that Bz acts as the acceptor. While the S(1) <- S(0) electronic origin of Bz is forbidden and is not observable in the gas-phase, the UV spectrum of Pyr center dot Bz in the same region exhibits a weak 0(0)(0) band that is red-shifted by 58 cm(-1) relative to that of Bz (38 086 cm(-1)). The origin appears due to symmetry-breaking of the p-electron system of Bz by the asymmetric pyrrole NH center dot center dot center dot pi hydrogen bond. This contrasts with (Bz)(2), which does not exhibit a 0(0)(0) band. The Bz moiety in Pyr center dot Bz exhibits a 6a(0)(1) band at 0(0)(0) + 518 cm(-1) that is about 20x more intense than the origin band. The symmetry breaking by the NH center dot center dot center dot pi hydrogen bond splits the degeneracy of the v(6)(e(2g)) vibration, giving rise to 6a' and 6b' sub-bands that are spaced by similar to 6 cm(-1). Both the 0(0)(0) and 6(0)(1) bands of Pyr center dot Bz carry a progression in the low-frequency (10 cm(-1)) excited-state tipping vibration omega', in agreement with the change of the omega tipping angle predicted by SCS-MP2 and SCS-CC2 calculations.

A tetrathiafulvalene-functionalized naphthalene diimide: synthesis, electrochemical and photophysical properties

A tetrathiafulvalene donor has been attached to the naphthalene diimide core via a rigid bridge affording a new planar molecular dyad. Its electronic properties have been studied experimentally by the combination of electrochemistry and UV-vis-NIR spectroscopy. Various electronic excited charge-transfer states are generated in different oxidation states, leading to almost full absorption in the visible to near-IR region with high extinction coefficients. The observed electronic properties are explained on the basis of density-functional-theory. In particular, the oxidized radical species show a strong tendency to undergo aggregation, in which the long-distance attractive interactions overcome the electrostatic repulsions. (C) 2011 Elsevier Ltd. All rights reserved.

A tetrathiafulvalene-functionalized schiff base macrocycle: synthesis, electrochemical, and photophysical properties

Highly selective formation of 2+2 macrocycle 1 from 2,5-bis(3-formyl-2-hydroxyphenyl)-1,3,4-oxadiazole and a diamine-functionalized tetrathiafulvalene derivative is reported. Its electronic properties have been studied experimentally by the combination of electrochemistry and UV-vis-NIR spectroscopy. Particularly, its largely extended pi-conjugation renders this novel macrocycle simultaneously a good multielectron donor and a strong chromophore, which is rationalized on the basis of density functional theory. (C) 2011 Elsevier Ltd. All rights reserved.

Excited state properties of formamide in water solution: An ab initio study

We present ab initio quantum calculation of the optical properties of formamide in vapor phase and in water solution. We employ time dependent density functional theory for the isolated molecule and many-body perturbation theory methods for the system in solution. An average over several molecular dynamics snapshots is performed to take into account the disorder of the liquid. We find that the excited stateproperties of the gas-phase formamide are strongly modified by the presence of the water solvent: the geometry of the molecule is distorted and the electronic and optical properties are severely modified. The important interaction among the formamide and the water molecules forces us to use fully quantum methods for the calculation of the excited stateproperties of this system. The excitonic wave function is localized both on the solute and on part of the solvent.

Structural Analysis of alpha-Fetoprotein (AFP)-like Peptides with Anti-Breast-Cancer Properties

The abundance of alpha-fetoprotein (AFP), a natural protein produced by the fetal yolk sac during pregnancy, correlates with lower incidence of estrogen receptor positive (ER+) breast cancer. The pharmacophore region of AFP has been narrowed down to a four amino acid (AA) region in the third domain of the 591 AA peptide. Our computational study focuses on a 4-mer segment consisting of the amino acids threonine-proline-valine-asparagine (TPVN). We have run replica exchange molecular dynamics (REMD) simulations and used 120 configurational snapshots from the total trajectory as starting configurations for quantum chemical calculations. We optimized structures using semiempirical (PM3, PM6, PM6-D2, PM6-H2, PM6-DH+, PM6-DH2) and density functional methods (TPSS, PBE0, M06-2X). By comparing the accuracy of these methods against RI-MP2 benchmarks, we devised a protocol for calculating the lowest energy conformers of these peptides accurately and efficiently. This protocol screens out high-energy conformers using lower levels of theory and outlines a general method for predicting small peptide structures.

A Computational Analysis of the Collision-Induced Dissociation Mechanisms of the Dipeptide Glycine-Serine Under Mass Spectrometry

Collision-induced dissociation (CID) of peptides using tandem mass spectrometry (MS) has been used to determine the identity of peptides and other large biological molecules. Mass spectrometry (MS) is a useful tool for determining the identity of molecules based on their interaction with electromagnetic fields. If coupled with another method like infrared (IR) vibrational spectroscopy, MS can provide structural information, but in its own right, MS can only provide the mass-to-charge (m/z) ratio of the fragments produced, which may not be enough information to determine the mechanism of the collision-induced dissociation (CID) of the molecule. In this case, theoretical calculations provide a useful companion for MS data and yield clues about the energetics of the dissociation. In this study, negative ion electrospray tandem MS was used to study the CID of the deprotonated dipeptide glycine-serine (Gly-Ser). Though negative ion MS is not as popular a choice as positive ion MS, studies by Bowie et al. show that it yields unique clues about molecular structure which complement positive ion spectroscopy, such as characteristic fragmentations like the loss of formaldehyde from the serine residue.2 The increase in the collision energy in the mass spectrometer alters the flexibility of the dipeptide backbone, enabling isomerizations (reactions not resulting in a fragment loss) and dissociations to take place. The mechanism of the CID of Gly-Ser was studied using two computational methods, B3LYP/6-311+G* and M06-2X/6-311++G**. The main pathway for molecular dissociation was analyzed in 5 conformers in an attempt to verify the initial mechanism proposed by Dr. James Swan after examination of the MS data. The results suggest that the loss of formaldehyde from serine, which Bowie et al. indicates is a characteristic of the presence of serine in a protein residue, is an endothermic reaction that is made possible by the conversion of the translational energy of the ion into internal energy as the ion collides with the inert collision gas. It has also been determined that the M06-2X functional¿s improved description of medium and long-range correlation makes it more effective than the B3LYP functional at finding elusive transition states. M06-2X also more accurately predicts the energy of those transition states than does B3LYP. A second CID mechanism, which passes through intermediates with the same m/z ratio as the main pathway for molecular dissociation, but different structures, including a diketopiperazine intermediate, was also studied. This pathway for molecular dissociation was analyzed with 3 conformers and the M06-2X functional, due to its previously determined effectiveness. The results suggest that the latter pathway, which meets the same intermediate masses as the first mechanism, is lower in overall energy and therefore a more likely pathway of dissociation than the first mechanism.