Vincenzo Barone

Anharmonic theoretical simulations of infrared spectra of halogenated organic compounds

The recent implementation of the computation of infrared (IR) intensities beyond the double-harmonic approximation [J. Bloino and V. Barone, J. Chem. Phys. 136, 124108 (2012)] paved the route to routine calculations of infrared spectra for a wide set of molecular systems. Halogenated organic compounds represent an interesting class of molecules, from both an atmospheric and computational point of view, due to the peculiar chemical features related to the halogen atoms. In this work, we simulate the IR spectra of eight halogenated molecules (CH2F2, CHBrF2, CH2DBr, CF3Br, CH2CHF, CF2CFCl, cis-CHFCHBr, cis-CHFCHI), using two common hybrid and double-hybrid density functionals in conjunction with both double- and triple-ζ quality basis sets (SNSD and cc-pVTZ) as well as employing the coupled-cluster theory with basis sets of at least triple-ζ quality. Finally, we compare our results with available experimental spectra, with the aim of checking the accuracy and the performances of the computational approaches.

Gas-phase formation of the prebiotic molecule formamide: insights from new quantum computations

New insights into the formation of interstellar formamide, a species of great relevance in prebiotic chemistry, are provided by electronic structure and kinetic calculations for the reaction NH2 + H2CO -> NH2CHO + H. Contrarily to what previously suggested, this reaction is essentially barrierless and can, therefore, occur under the low temperature conditions of interstellar objects thus providing a facile formation route of formamide. The rate coefficient parameters for the reaction channel leading to NH2CHO + H have been calculated to be A = 2.6x10^{-12} cm^3 s^{-1}, beta = -2.1 and gamma = 26.9 K in the range of temperatures 10-300 K. Including these new kinetic data in a refined astrochemical model, we show that the proposed mechanism can well reproduce the abundances of formamide observed in two very different interstellar objects: the cold envelope of the Sun-like protostar IRAS16293-2422 and the molecular shock L1157-B2. Therefore, the major conclusion of this Letter is that there is no need to invoke grain-surface chemistry to explain the presence of formamide provided that its precursors, NH2 and H2CO, are available in the gas-phase.

ACCURATE SPECTROSCOPIC CHARACTERIZATION OF OXIRANE: A VALUABLE ROUTE TO ITS IDENTIFICATION IN TITAN'S ATMOSPHERE AND THE ASSIGNMENT OF UNIDENTIFIED INFRARED BANDS

In an effort to provide an accurate spectroscopic characterization of oxirane, state-of-the-art computational methods and approaches have been employed to determine highly accurate fundamental vibrational frequencies and rotational parameters. Available experimental data were used to assess the reliability of our computations, and an accuracy on average of 10 cm-1 for fundamental transitions as well as overtones and combination bands has been pointed out. Moving to rotational spectroscopy, relative discrepancies of 0.1%, 2%-3%, and 3%-4% were observed for rotational, quartic, and sextic centrifugal-distortion constants, respectively. We are therefore confident that the highly accurate spectroscopic data provided herein can be useful for identification of oxirane in Titans atmosphere and the assignment of unidentified infrared bands. Since oxirane was already observed in the interstellar medium and some astronomical objects are characterized by very high D/H ratios, we also considered the accurate determination of the spectroscopic parameters for the mono-deuterated species, oxirane-d1. For the latter, an empirical scaling procedure allowed us to improve our computed data and to provide predictions for rotational transitions with a relative accuracy of about 0.02% (i.e., an uncertainty of about 40 MHz for a transition lying at 200 GHz).

CYANOMETHANIMINE ISOMERS IN COLD INTERSTELLAR CLOUDS: INSIGHTS FROM ELECTRONIC STRUCTURE AND KINETIC CALCULATIONS

New insights into the formation of interstellar cyanomethanimine, a species of great relevance in prebiotic chemistry, are provided by electronic structure and kinetic calculations for the reaction CN + CH2 = NH. This reaction is a facile formation route of Z,E-C-cyanomethanimine, even under the extreme conditions of density and temperature typical of cold interstellar clouds. E-C-cyanomethanimine has been recently identified in Sgr B2(N) in the Green Bank Telescope (GBT) PRIMOS survey by P. Zaleski et al. and no efficient formation routes have been envisaged so far. The rate coefficient expression for the reaction channel leading to the observed isomer E-C-cyanomethanimine is 3.15 × 10-10 × (T/300)0.152 × e(−0.0948/T). According to the present study, the more stable Z-C-cyanomethanimine isomer is formed with a slightly larger yield (4.59 × 10−10 × (T/300)0.153 × e(−0.0871/T). As the detection of E-isomer is favored due to its larger dipole moment, the missing detection of the Z-isomer can be due to the sensitivity limit of the GBT PRIMOS survey and the detection of the Z-isomer should be attempted with more sensitive instrumentation. The CN + CH2 = NH reaction can also play a role in the chemistry of the upper atmosphere of Titan where the cyanomethanimine products can contribute to the buildup of the observed nitrogen-rich organic aerosols that cover the moon.

Reassessment of the Thermodynamic, Kinetic, and Spectroscopic Features of Cyanomethanimine Derivatives: A Full Anharmonic Perturbative Treatment.

Herein we report a full thermodynamic and vibrational investigation of C-cyanomethanimine isomers rooted into the Density Functional Theory (DFT) and the second-order vibrational perturbation theory (VPT2). We show that an anharmonic treatment affects dramatically the vibrational behavior of the molecules, especially thanks to the inclusion of interaction terms between the various modes. Furthermore, the equilibrium constant between the isomers, as well as the rate constant, have been obtained at both harmonic and anharmonic levels showing, as expected, slight but non-negligible differences. To support our investigation, dispersion effects have been employed.

Methods for Calculating Partition Functions of Molecules Involving Large Amplitude and/or Anharmonic Motions.

We present a method for calculating partition functions taking into account anharmonic contributions for systems involving both small-amplitude vibrations and hindered rotations. The Wang-Landau scheme is used in the first case, while two alternative schemes are used for hindered rotation based on imaginary time propagation and fitting of the exact energy levels as a function of quantum number. These two schemes are shown to be complementary in their ranges of applicability (in terms of the torsional rotational constant and the relevant potential). Partition functions for four different molecules are calculated and compared to simpler ones obtained using a harmonic model.

Accurate Spectroscopic Characterization of Protonated Oxirane: A Potential Prebiotic Species in Titan's Atmosphere

An accurate spectroscopic characterization of protonated oxirane has been carried out by means of state-of-the-art computational methods and approaches. The calculated spectroscopic parameters from our recent computational investigation of oxirane together with the corresponding experimental data available were used to assess the accuracy of our predicted rotational and IR spectra of protonated oxirane. We found an accuracy of about 10 cm-1 for vibrational transitions (fundamentals as well as overtones and combination bands) and, in relative terms, of 0.1% for rotational transitions. We are therefore confident that the spectroscopic data provided herein are a valuable support for the detection of protonated oxirane not only in Titans atmosphere but also in the interstellar medium.

A new Gaussian MCTDH program: implementation and validation on the levels of the water and glycine molecules.

We report the main features of a new general implementation of the Gaussian Multi-Configuration Time-Dependent Hartree model. The code allows effective computations of time-dependent phenomena, including calculation of vibronic spectra (in one or more electronic states), relative state populations, etc. Moreover, by expressing the Dirac-Frenkel variational principle in terms of an effective Hamiltonian, we are able to provide a new reliable estimate of the representation error. After validating the code on simple one-dimensional systems, we analyze the harmonic and anharmonic vibrational spectra of water and glycine showing that reliable and converged energy levels can be obtained with reasonable computing resources. The data obtained on water and glycine are compared with results of previous calculations using the vibrational second-order perturbation theory method. Additional features and perspectives are also shortly discussed.

Nonadiabatic photodynamics of phenol on a realistic potential energy surface by a novel multilayer Gaussian MCTDH program

We report the main features of a new implementation of the Gaussian Multi-Configuration Time-Dependent Hartree (G-MCTDH) model. The code allows effective computations of time-dependent phenomena, including calculation of vibronic spectra (in one or more electronic states), relative state populations etc., with the possibility of a multilayer formulation. We have validated the code on the diabatic surfaces recently published by Truhlar and coworkers to study the nonadiabatic photodynamics of phenol. Using an Ehrenfest-like, single-nuclear-configuration (but in a fully quantum formalism) model we calculate the optical spectrum and relative state populations of the system as a function of time.