1
|
Xu R, Jiang Z, Yang Q, Bloino J, Biczysko M. Harmonic and anharmonic vibrational computations for biomolecular building blocks: Benchmarking DFT and basis sets by theoretical and experimental IR spectrum of glycine conformers. J Comput Chem 2024; 45:1846-1869. [PMID: 38682874 DOI: 10.1002/jcc.27377] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2024] [Revised: 04/01/2024] [Accepted: 04/02/2024] [Indexed: 05/01/2024]
Abstract
Advanced vibrational spectroscopic experiments have reached a level of sophistication that can only be matched by numerical simulations in order to provide an unequivocal analysis, a crucial step to understand the structure-function relationship of biomolecules. While density functional theory (DFT) has become the standard method when targeting medium-size or larger systems, the problem of its reliability and accuracy are well-known and have been abundantly documented. To establish a reliable computational protocol, especially when accuracy is critical, a tailored benchmark is usually required. This is generally done over a short list of known candidates, with the basis set often fixed a priori. In this work, we present a systematic study of the performance of DFT-based hybrid and double-hybrid functionals in the prediction of vibrational energies and infrared intensities at the harmonic level and beyond, considering anharmonic effects through vibrational perturbation theory at the second order. The study is performed for the six-lowest energy glycine conformers, utilizing available "state-of-the-art" accurate theoretical and experimental data as reference. Focusing on the most intense fundamental vibrations in the mid-infrared range of glycine conformers, the role of the basis sets is also investigated considering the balance between computational cost and accuracy. Targeting larger systems, a broad range of hybrid schemes with different computational costs is also tested.
Collapse
Affiliation(s)
- Ruiqin Xu
- Department of Physics, College of Sciences, Shanghai University, Shanghai, China
| | | | - Qin Yang
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Science, Prague, Czechia
| | - Julien Bloino
- Classe di Scienze, Scuola Normale Superiore, Pisa, Italy
| | - Malgorzata Biczysko
- Department of Physics, College of Sciences, Shanghai University, Shanghai, China
| |
Collapse
|
2
|
Puzzarini C, Alessandrini S. Carbamic acid and its dimer: A computational study. J Comput Chem 2024. [PMID: 38970400 DOI: 10.1002/jcc.27442] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2024] [Revised: 05/15/2024] [Accepted: 05/17/2024] [Indexed: 07/08/2024]
Abstract
A recent work by Marks et al. on the formation of carbamic acid in NH3 $$ {}_3 $$ -CO2 $$ {}_2 $$ interstellar ices pointed out its stability in the gas phase and the concomitant production of its dimer. Prompted by these results and the lack of information on these species, we have performed an accurate structural, energetic and spectroscopic investigation of carbamic acid and its dimer. For the former, the structural and spectroscopic characterization employed composite schemes based on coupled cluster (CC) calculations that account for the extrapolation to the complete basis set limit and core correlation effects. A first important outcome is the definitive confirmation of the nonplanarity of carbamic acid, then followed by an accurate estimate of its rotational and vibrational spectroscopy parameters. As far as the carbamic acid dimer is concerned, the investigation started from the identification of its most stable forms. For them, structure and vibrational properties have been evaluated using density functional theory, while a composite scheme rooted in CC theory has been employed for the energetic characterization. Our results allowed us to provide a better interpretation of the feature observed in the recent experiment mentioned above.
Collapse
Affiliation(s)
- Cristina Puzzarini
- Department of Chemistry "Giacomo Ciamician", University of Bologna, Bologna, Italy
| | - Silvia Alessandrini
- Department of Chemistry "Giacomo Ciamician", University of Bologna, Bologna, Italy
| |
Collapse
|
3
|
Nonne M, Melosso M, Tonolo F, Bizzocchi L, Alessandrini S, Guillemin JC, Dore L, Puzzarini C. Tracing Prebiotic Molecules: Rotational Spectroscopy of Deuterated Glycolaldehyde and ( Z)-1,2-Ethenediol. J Phys Chem A 2024; 128:4850-4858. [PMID: 38842131 DOI: 10.1021/acs.jpca.4c02533] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2024]
Abstract
Glycolaldehyde, an important prebiotic molecule, along with its monodeuterated species and its higher energy tautomer, (Z)-1,2-ethenediol, has been detected in the interstellar medium. Although the elemental D/H ratio in the universe is only ∼1.6 × 10-5, the deuterium relative abundance in interstellar molecules might be by far larger than this. As such, it provides a remarkable and almost unique diagnostic tool. In particular, it might help elucidate the reaction mechanisms that lead to the formation of the so-called complex organic molecules. It is therefore crucial to extend the census of the interstellar deuterated molecules. To this aim, in this work, we present for the first time a spectroscopic investigation of the rotational spectra of the CHDOD-CHO bideuterated variant of glycolaldehyde and of mono- and bideuterated species of (Z)-1,2-ethenediol (CHOD═CHOD, CHOD═CHOH, and CHOH═CHOD rotamers). For each species, more than a hundred transitions have been assigned. Their analysis led to the accurate determination of all rotational constants as well as quartic and sextic centrifugal distortion terms, thus providing spectroscopic line catalogs suitable for supporting astronomical searches. In addition, the rotational constants of the bideuterated glycolaldehyde isotopologue studied in this work allowed us to improve the semiexperimental equilibrium structure determination for this molecule.
Collapse
Affiliation(s)
- Michela Nonne
- Scuola Superiore Meridionale, Largo San Marcellino 10, 80138 Naples, Italy
| | - Mattia Melosso
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy
| | - Francesca Tonolo
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy
| | - Luca Bizzocchi
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy
| | - Silvia Alessandrini
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy
| | - Jean-Claude Guillemin
- Univ Rennes, Ecole Nationale Supérieure de Chimie de Rennes, CNRS, ISCR-UMR6226, F-35000 Rennes, France
| | - Luca Dore
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy
| | - Cristina Puzzarini
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy
| |
Collapse
|
4
|
Mandelli G, Corneo L, Aieta C. Coupled Cluster Semiclassical Estimates of Experimental Reaction Rates: The Interconversion of Glycine Conformer VIp to Ip. J Phys Chem Lett 2023; 14:9996-10002. [PMID: 37906174 PMCID: PMC10641884 DOI: 10.1021/acs.jpclett.3c02560] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2023] [Revised: 10/20/2023] [Accepted: 10/23/2023] [Indexed: 11/02/2023]
Abstract
We apply the full-dimensional Semiclassical Transition State Theory (SCTST) to estimate the rate constant of glycine molecule interconversion between the VIp and Ip conformers. We have reached an electronic structure accuracy up to the explicitly correlated Coupled Cluster method (CCSD(T)-F12b/cc-pVDZ-F12) thanks to our parallel implementation. The reaction has been experimentally investigated in the literature and is known to proceed by quantum mechanical tunneling. The SCTST rates improve over other theoretical methods, and our results align with the experimental measurements, thus confirming the accuracy of the fully coupled anharmonic semiclassical tunneling treatment, providing that the level of electronic structure theory gives a reliable estimate of the reaction barrier height and shape. The comparison with experimental half-life times supports the validity of SCTST for glycine VIp-Ip conformer conversion in the cryogenic temperature range, where this theory is usually not considered applicable due to the onset of the deep tunneling regime.
Collapse
Affiliation(s)
- Giacomo Mandelli
- Dipartimento di Chimica, Università degli Studi di Milano, via C. Golgi 19, Milano 20133, Italy
| | - Luca Corneo
- Dipartimento di Chimica, Università degli Studi di Milano, via C. Golgi 19, Milano 20133, Italy
| | - Chiara Aieta
- Dipartimento di Chimica, Università degli Studi di Milano, via C. Golgi 19, Milano 20133, Italy
| |
Collapse
|
5
|
Di Grande S, Kállay M, Barone V. Accurate thermochemistry at affordable cost by means of an improved version of the JunChS-F12 model chemistry. J Comput Chem 2023; 44:2149-2157. [PMID: 37432050 DOI: 10.1002/jcc.27187] [Citation(s) in RCA: 9] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2023] [Accepted: 06/23/2023] [Indexed: 07/12/2023]
Abstract
The junChS-F12 composite method has been improved by means of the latest implementation of the CCSD(F12*)(T+) ansatz and validated for the thermochemistry of molecules containing atoms of the first three rows of the periodic table. A thorough benchmark showed that this model, in conjunction with cost-effective revDSD-PBEP86-D3(BJ) reference geometries, offers an optimal compromise between accuracy and computational cost. If improved geometries are sought, the most effective option is to add MP2-F12 core-valence correlation corrections to CCSD(T)-F12b/jun-cc-pVTZ geometries without the need of performing any extrapolation to the complete basis set limit. In the same vein, CCSD(T)-F12b/jun-cc-pVTZ harmonic frequencies are remarkably accurate without any additional contribution. Pilot applications to noncovalent intermolecular interactions, conformational landscapes, and tautomeric equilibria confirm the effectiveness and reliability of the model.
Collapse
Affiliation(s)
- S Di Grande
- Classe di Scienze, Scuola Normale Superiore di Pisa, Pisa, Italy
- Scuola Superiore Meridionale, Napoli, Italy
| | - M Kállay
- Department of Physical Chemistry and Materials Science, Faculty of Chemical Technology and Biotechnology, Budapest University of Technology and Economics, Budapest, Hungary
- ELKH-BME Quantum Chemistry Research Group, Budapest, Hungary
- MTA-BME Lendület Quantum Chemistry Research Group, Budapest, Hungary
| | - V Barone
- Classe di Scienze, Scuola Normale Superiore di Pisa, Pisa, Italy
| |
Collapse
|
6
|
Barone V, Uribe Grajales LM, Di Grande S, Lazzari F, Mendolicchio M. DFT Meets Wave-Function Methods for Accurate Structures and Rotational Constants of Histidine, Tryptophan, and Proline. J Phys Chem A 2023; 127:7534-7543. [PMID: 37665117 PMCID: PMC10510395 DOI: 10.1021/acs.jpca.3c04227] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2023] [Revised: 08/12/2023] [Indexed: 09/05/2023]
Abstract
A new computational strategy has been applied to the conformational and spectroscopic properties in the gas phase of amino acids with very distinctive features, ranging from different tautomeric forms (histidine) to ring puckering (proline), and heteroaromatic structures with non-equivalent rings (tryptophan). The integration of modern double-hybrid functionals and wave-function composite methods has allowed us to obtain accurate results for a large panel of conformers with reasonable computer times. The remarkable agreement between computations and microwave experiments allows an unbiased interpretation of the latter in terms of stereoelectronic effects.
Collapse
Affiliation(s)
- Vincenzo Barone
- Scuola
Normale Superiore di Pisa, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Lina Marcela Uribe Grajales
- Scuola
Normale Superiore di Pisa, Piazza dei Cavalieri 7, 56126 Pisa, Italy
- Scuola
Superiore Meridionale, Largo San Marcellino 10, 80138 Napoli, Italy
| | - Silvia Di Grande
- Scuola
Normale Superiore di Pisa, Piazza dei Cavalieri 7, 56126 Pisa, Italy
- Scuola
Superiore Meridionale, Largo San Marcellino 10, 80138 Napoli, Italy
| | - Federico Lazzari
- Scuola
Normale Superiore di Pisa, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Marco Mendolicchio
- Scuola
Normale Superiore di Pisa, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
7
|
Nelson PM, Glick ZL, Sherrill CD. Approximating large-basis coupled-cluster theory vibrational frequencies using focal-point approximations. J Chem Phys 2023; 159:094104. [PMID: 37655773 DOI: 10.1063/5.0168608] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2023] [Accepted: 08/09/2023] [Indexed: 09/02/2023] Open
Abstract
The focal-point approximation can be used to estimate a high-accuracy, slow quantum chemistry computation by combining several lower-accuracy, faster computations. We examine the performance of focal-point methods by combining second-order Møller-Plesset perturbation theory (MP2) with coupled-cluster singles, doubles, and perturbative triples [CCSD(T)] for the calculation of harmonic frequencies and that of fundamental frequencies using second-order vibrational perturbation theory (VPT2). In contrast to standard CCSD(T), the focal-point CCSD(T) method approaches the complete basis set (CBS) limit with only triple-ζ basis sets for the coupled-cluster portion of the computation. The predicted harmonic and fundamental frequencies were compared with the experimental values for a set of 20 molecules containing up to six atoms. The focal-point method combining CCSD(T)/aug-cc-pV(T + d)Z with CBS-extrapolated MP2 has mean absolute errors vs experiment of only 7.3 cm-1 for the fundamental frequencies, which are essentially the same as the mean absolute error for CCSD(T) extrapolated to the CBS limit using the aug-cc-pV(Q + d)Z and aug-cc-pV(5 + d)Z basis sets. However, for H2O, the focal-point procedure requires only 3% of the computation time as the extrapolated CCSD(T) result, and the cost savings will grow for larger molecules.
Collapse
Affiliation(s)
- Philip M Nelson
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Zachary L Glick
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| |
Collapse
|
8
|
Barone V. Accurate structures and spectroscopic parameters of α,α-dialkylated α-amino acids in the gas-phase: a joint venture of DFT and wave-function composite methods. Phys Chem Chem Phys 2023; 25:22768-22774. [PMID: 37591810 DOI: 10.1039/d3cp02503a] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/19/2023]
Abstract
Accurate computations of structural, conformational and spectroscopic properties in the gas phase have been performed for two α,α-dialkylated α-amino acids, namely aminoisobutyric acid and cyclopropylglycine. Thanks to the integration of modern double hybrid functionals and wave-function methods, several low-energy structures of the title molecules could be analyzed employing standard computer resources. The computed features of all the most stable conformers of the target amino acids closely match the corresponding spectroscopic parameters issued from microwave spectroscopic studies in the gas-phase. Together with their intrinsic interest, the accuracy of the results obtained with reasonable computer times paves the way for accurate investigations of other flexible bricks of life.
Collapse
Affiliation(s)
- Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy.
| |
Collapse
|
9
|
Barone V, Fusè M. Accurate Structures and Spectroscopic Parameters of Phenylalanine and Tyrosine in the Gas Phase: A Joint Venture of DFT and Composite Wave-Function Methods. J Phys Chem A 2023; 127:3648-3657. [PMID: 37052318 PMCID: PMC10150396 DOI: 10.1021/acs.jpca.3c01174] [Citation(s) in RCA: 11] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/14/2023]
Abstract
A general strategy for the accurate computation of conformational and spectroscopic properties of flexible molecules in the gas phase is applied to two representative proteinogenic amino acids with aromatic side chains, namely, phenylalanine and tyrosine. The main features of all the most stable conformers predicted by this computational strategy closely match those of the species detected in microwave and infrared experiments. Together with their intrinsic interest, the accuracy of the results obtained with reasonable computer times paves the route for accurate investigations of other flexible bricks of life.
Collapse
Affiliation(s)
- Vincenzo Barone
- Scuola Normale Superiore di Pisa, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Marco Fusè
- DMMT-sede Europa, Università di Brescia, Viale Europa 11, 25121 Brescia, Italy
| |
Collapse
|
10
|
Alessandrini S, Melosso M, Rivilla VM, Bizzocchi L, Puzzarini C. Computational Protocol for the Identification of Candidates for Radioastronomical Detection and Its Application to the C3H3NO Family of Isomers. Molecules 2023; 28:molecules28073226. [PMID: 37049990 PMCID: PMC10096335 DOI: 10.3390/molecules28073226] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2023] [Revised: 03/30/2023] [Accepted: 03/31/2023] [Indexed: 04/08/2023] Open
Abstract
The C3H3NO family of isomers is relevant in astrochemistry, even though its members are still elusive in the interstellar medium. To identify the best candidate for astronomical detection within this family, we developed a new computational protocol based on the minimum-energy principle. This approach aims to identify the most stable isomer of the family and consists of three steps. The first step is an extensive investigation that characterizes the vast number of compounds having the C3H3NO chemical formula, employing density functional theory for this purpose. The second step is an energy refinement, which is used to select isomers and relies on coupled cluster theory. The last step is a structural improvement with a final energy refinement that provides improved energies and a large set of accurate spectroscopic parameters for all isomers lying within 30 kJ mol−1 above the most stable one. According to this protocol, vinylisocyanate is the most stable isomer, followed by oxazole, which is about 5 kJ mol−1 higher in energy. The other stable species are pyruvonitrile, cyanoacetaldehyde, and cyanovinylalcohol. For all of these species, new computed rotational and vibrational spectroscopic data are reported, which complement those already available in the literature or fill current gaps.
Collapse
Affiliation(s)
- Silvia Alessandrini
- Dipartimento di Chimica “Giacomo Ciamician”, Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy
| | - Mattia Melosso
- Dipartimento di Chimica “Giacomo Ciamician”, Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy
| | - Víctor M. Rivilla
- Centro de Astrobiología (CSIC-INTA), Ctra. de Ajalvir km. 4, Torrejón de Ardoz, 28850 Madrid, Spain
| | - Luca Bizzocchi
- Dipartimento di Chimica “Giacomo Ciamician”, Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy
| | - Cristina Puzzarini
- Dipartimento di Chimica “Giacomo Ciamician”, Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy
| |
Collapse
|
11
|
Nacsa AB, Kígyósi M, Czakó G. Protonation of serine: conformers, proton affinities and gas-phase basicities at the "gold standard" and beyond. Phys Chem Chem Phys 2023; 25:8891-8902. [PMID: 36916632 DOI: 10.1039/d3cp00612c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/09/2023]
Abstract
The potential energy surfaces (PESs) of serine and its protonated counterparts are investigated to determine the structures of the minima. A total of 95 neutral serine, 15 N-(amino-) and 46 O-(carbonyl-)protonated serine conformers are found. Their relative energies, geometries and harmonic vibrational frequencies are determined at the MP2/aug-cc-pVDZ level of theory. To obtain highly accurate thermodynamic values, further computations are performed: the ten conformers with the lowest relative energies from each molecule type (neutral, N- and O-protonated) are further optimized using the explicitly correlated CCSD(T)-F12a/cc-pVDZ-F12 method (for neutral serine, harmonic vibrational frequencies were also computed). In addition, auxiliary corrections were determined: basis-set effects up to CCSD(T)-F12b/cc-pVQZ-F12, electron correlation effects up to CCSDT(Q), core correlation and second-order Douglas-Kroll relativistic effects along with zero-point energy contributions. Two important thermodynamic parameters (at 298.15 K), proton affinity (PA)/gas-phase basicity (GB) are calculated considering the two different protonation sites: 218.05 ± 0.2/209.86 ± 0.6 kcal mol-1 and 205.87 ± 0.2/196.36 ± 0.3 kcal mol-1 for the amino and carbonyl sites, respectively. The uncertainty of the determined values was approximated based on various sources including auxiliary corrections, basis-set effects, harmonic vibrational frequencies.
Collapse
Affiliation(s)
- András B Nacsa
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary.
| | - Máté Kígyósi
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary.
| | - Gábor Czakó
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary.
| |
Collapse
|
12
|
Jiang N, Melosso M, Alessandrini S, Bizzocchi L, Martin-Drumel MA, Pirali O, Puzzarini C. Insights into the molecular structure and infrared spectrum of the prebiotic species aminoacetonitrile. Phys Chem Chem Phys 2023; 25:4754-4763. [PMID: 36691972 DOI: 10.1039/d2cp05179f] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
Aminoacetonitrile is an interstellar molecule with a prominent prebiotic role, already detected in the chemically-rich molecular cloud Sagittarius B2(N) and postulated to be present in the atmosphere of the largest Saturn's moon, Titan. To further support its observation in such remote environments and laboratory experiments aimed at improving our understanding of interstellar chemistry, we report a thorough spectroscopic and structural characterization of aminoacetonitrile. Equilibrium geometry, fundamental bands as well as spectroscopic and molecular parameters have been accurately computed by exploiting a composite scheme rooted in the coupled-cluster theory that accounts for the extrapolation to the complete basis set limit and core-correlation effects. In addition, a semi-experimental approach that combines ground-state rotational constants for different isotopic species and calculated vibrational corrections has been employed for the structure determination. From the experimental side, we report the analysis of the three strongest fundamental bands of aminoacetonitrile observed between 500 and 1000 cm-1 in high-resolution infrared spectra. More generally, all computed band positions are in excellent agreement with the present and previous experiments. The only exception is the ν15 band, for which we provide a revision of the experimental assignment, now in good agreement with theory.
Collapse
Affiliation(s)
- Ningjing Jiang
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy.
| | - Mattia Melosso
- Scuola Superiore Meridionale, Largo San Marcellino 10, 80138 Naples, Italy.
| | - Silvia Alessandrini
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy. .,Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Luca Bizzocchi
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy.
| | | | - Olivier Pirali
- Université Paris-Saclay, CNRS, Institut des Sciences Moléculaires d'Orsay, 91405 Orsay, France.,SOLEIL Synchrotron, AILES beamline, l'Orme des Merisiers, 91190 Saint-Aubin, Gif-sur-Yvette, France
| | - Cristina Puzzarini
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy.
| |
Collapse
|
13
|
Gruber B, Czakó G. High-level ab initio mapping of the multiple H-abstraction pathways of the OH + glycine reaction. Phys Chem Chem Phys 2023; 25:5271-5281. [PMID: 36723222 DOI: 10.1039/d2cp03049g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
We perform a systematic search in the transition-state (TS) and product-channel complex (MIN) regions of the multi-channel OH + glycine → H2O + H2N-CH-COOH (CH)/HN-CH2-COOH (NH)/H2N-CH2-COO (COOH) reactions. Geometry optimizations reveal {7, 3, 3} CH-TS, {2, 2, 2} CH-MIN, {17, 10, 5} NH-TS, {35, 19, 19} NH-MIN, and {6, 5, 5} COOH-TS conformers at the {MP2/3-21G, MP2/aug-cc-pVDZ, CCSD(T)-F12b/aug-cc-pVDZ} levels of theory as well as 2 additional CH-TSs based on chemical intuition. The benchmark relative energies of the TS, MIN, and product conformers are obtained by considering basis set effects up to aug-cc-pVQZ using the explicitly-correlated CCSD(T)-F12b method, post-(T) correlation up to CCSDT(Q), core correlation, scalar relativistic effects, spin-orbit coupling, and zero-point energy corrections. All the CH [ΔEe(ΔH0) = -38.54(-38.61) kcal mol-1], NH [ΔEe(ΔH0) = -16.72(-17.98) kcal mol-1], and COOH [ΔEe = -4.98 kcal mol-1] reactions are exothermic and proceed via shallow, usually negative, classical(adiabatic) barriers of -0.37(-0.95), -1.91(-2.48), and 1.02(-0.57) kcal mol-1, respectively. In the entrance channel MRCI/aug-cc-pVTZ computations reveal several complexes with reactive(non-reactive) arrangements and binding energies of 1.0, 1.6, 3.3, (5.2 and 5.9) kcal mol-1, stabilized by CH⋯OH, NH⋯OH, COOH⋯OH, (OH⋯OC and OH⋯N) hydrogen bonds, respectively.
Collapse
Affiliation(s)
- Balázs Gruber
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary.
| | - Gábor Czakó
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary.
| |
Collapse
|
14
|
Franke PR, Stanton JF. Rotamers of Methanediol: Composite Ab Initio Predictions of Structures, Frequencies, and Rovibrational Constants. J Phys Chem A 2023; 127:924-937. [PMID: 36657011 DOI: 10.1021/acs.jpca.2c06686] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
Geminal diols are known to be important intermediates in atmospheric ozonolysis and the aerosol cycle. Recently, the simplest member of this class, methanediol, was interrogated in the gas phase with infrared spectroscopy. To aid in future spectroscopic investigations of methanediol, including in the interstellar medium, we report fundamental frequencies and rovibrational constants for the two rotamers of this molecule using ab initio composite methods along with vibrational perturbation theory. Sensitivity of the predictions to the level of theory and the treatment of anharmonic resonances are carefully assessed. The OH stretching harmonic frequencies of both rotamers are particularly sensitive to the level of theory. The CH stretches of the Cs rotamer are sensitive to the treatment of anharmonic resonances with VPT2-based effective Hamiltonian models. Equilibrium bond distances and harmonic frequencies are converged conservatively to within 0.0005 Å and 3 cm-1, respectively. The effect of tunneling on the rotational constants is investigated with a 2D variational calculation, based on a relaxed hydroxyl torsional potential energy surface. Tunneling is found to be negligible in the lower energy C2 rotamer but should modify the rotational constants of the Cs rotamer on the order of MHz, giving rise to rotational line splittings of the same order. The rovibrational constants of the Cs rotamer are dominated by hydroxyl torsional effects, and here we see evidence for the breakdown of vibrational perturbation theory.
Collapse
Affiliation(s)
- Peter R Franke
- Department of Chemistry, University of Florida, Gainesville, Florida32611, United States
| | - John F Stanton
- Department of Chemistry, University of Florida, Gainesville, Florida32611, United States
| |
Collapse
|
15
|
Barone V, Fusè M, Lazzari F, Mancini G. Benchmark Structures and Conformational Landscapes of Amino Acids in the Gas Phase: A Joint Venture of Machine Learning, Quantum Chemistry, and Rotational Spectroscopy. J Chem Theory Comput 2023; 19:1243-1260. [PMID: 36731119 PMCID: PMC9979611 DOI: 10.1021/acs.jctc.2c01143] [Citation(s) in RCA: 17] [Impact Index Per Article: 17.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
The accurate characterization of prototypical bricks of life can strongly benefit from the integration of high resolution spectroscopy and quantum mechanical computations. We have selected a number of representative amino acids (glycine, alanine, serine, cysteine, threonine, aspartic acid and asparagine) to validate a new computational setup rooted in quantum-chemical computations of increasing accuracy guided by machine learning tools. Together with low-lying energy minima, the barriers ruling their interconversion are evaluated in order to unravel possible fast relaxation paths. Vibrational and thermal effects are also included in order to estimate relative free energies at the temperature of interest in the experiment. The spectroscopic parameters of all the most stable conformers predicted by this computational strategy, which do not have low-energy relaxation paths available, closely match those of the species detected in microwave experiments. Together with their intrinsic interest, these accurate results represent ideal benchmarks for more approximate methods.
Collapse
Affiliation(s)
- Vincenzo Barone
- Scuola
Normale Superiore di Pisa, piazza dei Cavalieri 7, 56126 Pisa, Italy,
| | - Marco Fusè
- DMMT-sede
Europa, Universitá di Brescia, viale Europa 11, 25121 Brescia, Italy
| | - Federico Lazzari
- Scuola
Normale Superiore di Pisa, piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Giordano Mancini
- Scuola
Normale Superiore di Pisa, piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
16
|
Fayaz A, Banik S, Kanchan Roy T. The importance of electron correlations on vibrational anharmonicities and potential energy surfaces. COMPUT THEOR CHEM 2023. [DOI: 10.1016/j.comptc.2023.114059] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
|
17
|
Puzzarini C, Stanton JF. Connections between the accuracy of rotational constants and equilibrium molecular structures. Phys Chem Chem Phys 2023; 25:1421-1429. [PMID: 36562443 DOI: 10.1039/d2cp04706c] [Citation(s) in RCA: 27] [Impact Index Per Article: 27.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Rotational spectroscopy is the technique of choice for investigating molecular structures in the gas phase. Indeed, rotational constants are strongly connected to the geometry of the molecular system under consideration. Therefore, they are powerful tools for assessing the accuracy that quantum chemical approaches can reach in structural determinations. In this review article, it is shown how it is possible to measure the accuracy of a computed equilibrium geometry based on the comparison of rotational constants. But, it is also addressed what accuracy is required by computations for providing molecular structures and thus rotational constants that are useful to experiment. Quantum chemical methodologies for obtaining the "0.1% accuracy" for rotational constants are reviewed for systems ranging in size from small molecules to small polycyclic aromatic hydrocarbons. This accuracy for systems containing two dozen or so atoms opens the way towards future applications such as the accurate characterization of non-covalent interactions, which play a key role in several biological and technological processes.
Collapse
Affiliation(s)
- Cristina Puzzarini
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, via F. Selmi 2, 40126, Bologna, Italy.
| | - John F Stanton
- Department of Chemistry, University of Florida, Gainesville, FL, 32611, USA.
| |
Collapse
|
18
|
Barone V, Di Grande S, Puzzarini C. Toward Accurate yet Effective Computations of Rotational Spectroscopy Parameters for Biomolecule Building Blocks. MOLECULES (BASEL, SWITZERLAND) 2023; 28:molecules28020913. [PMID: 36677970 PMCID: PMC9863398 DOI: 10.3390/molecules28020913] [Citation(s) in RCA: 9] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 01/01/2023] [Revised: 01/01/2023] [Accepted: 01/12/2023] [Indexed: 01/19/2023]
Abstract
The interplay of high-resolution rotational spectroscopy and quantum-chemical computations plays an invaluable role in the investigation of biomolecule building blocks in the gas phase. However, quantum-chemical methods suffer from unfavorable scaling with the dimension of the system under consideration. While a complete characterization of flexible systems requires an elaborate multi-step strategy, in this work, we demonstrate that the accuracy obtained by quantum-chemical composite approaches in the prediction of rotational spectroscopy parameters can be approached by a model based on density functional theory. Glycine and serine are employed to demonstrate that, despite its limited cost, such a model is able to predict rotational constants with an accuracy of 0.3% or better, thus paving the way toward the accurate characterization of larger flexible building blocks of biomolecules.
Collapse
Affiliation(s)
- Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-50126 Pisa, Italy
- Correspondence: (V.B.); (C.P.)
| | - Silvia Di Grande
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-50126 Pisa, Italy
- Scuola Superiore Meridionale, Largo San Marcellino 10, I-80138 Napoli, Italy
| | - Cristina Puzzarini
- Rotational and Computational Spectroscopy Lab, Department of Chemistry “Giacomo Ciamician”, University of Bologna, Via F. Selmi 2, I-40126 Bologna, Italy
- Correspondence: (V.B.); (C.P.)
| |
Collapse
|
19
|
Nacsa AB, Czakó G. Benchmark Ab Initio Determination of the Conformers, Proton Affinities, and Gas-Phase Basicities of Cysteine. J Phys Chem A 2022; 126:9667-9679. [PMID: 36524999 PMCID: PMC9806835 DOI: 10.1021/acs.jpca.2c07035] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
A systematic conformational mapping combined with literature data leads to 85 stable neutral cysteine conformers. The implementation of the same mapping process for the protonated counterparts reveals 21 N-(amino-), 64 O-(carbonyl-), and 37 S-(thiol-)protonated cysteine conformers. Their relative energies and harmonic vibrational frequencies are given at the MP2/aug-cc-pVDZ level of theory. Further benchmark ab initio computations are performed for the 10 lowest-lying neutral and protonated amino acid conformers (for each type) such as CCSD(T)-F12a/cc-pVDZ-F12 geometry optimizations (and frequency computations for cysteine) as well as auxiliary correction computations of the basis set effects up to CCSD(T)-F12b/cc-pVQZ-F12, electron correlation effects up to CCSDT(Q), core correlation effects, second-order Douglass-Kroll relativistic effects, and zero-point energy contributions. Boltzmann-averaged 0 (298.15) K proton affinity and [298.15 K gas-phase basicity] values of cysteine are predicted to be 214.96 (216.39) [208.21], 201.83 (203.55) [194.16], and 193.31 (194.74) [186.40] kcal/mol for N-, O-, and S-protonation, respectively, also considering the previously described auxiliary corrections.
Collapse
|
20
|
Melli A, Melosso M, Bizzocchi L, Alessandrini S, Jiang N, Tonolo F, Boi S, Castellan G, Sapienza C, Guillemin JC, Dore L, Puzzarini C. Rotational Spectra of Unsaturated Carbon Chains Produced by Pyrolysis: The Case of Propadienone, Cyanovinylacetylene, and Allenylacetylene. J Phys Chem A 2022; 126:6210-6220. [PMID: 36044202 PMCID: PMC9483987 DOI: 10.1021/acs.jpca.2c05018] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
Several interstellar molecules are highly reactive unsaturated
carbon chains, which are unstable under terrestrial conditions. Laboratory
studies in support of their detection in space thus face the issue
of how to produce these species and how to correctly model their rotational
energy levels. In this work, we introduce a general approach for producing
and investigating unsaturated carbon chains by means of selected test
cases. We report a comprehensive theoretical/experimental spectroscopic
characterization of three species, namely, propadienone, cyanovinylacetylene,
and allenylacetylene, all of them being produced by means of flash
vacuum pyrolysis of a suitable precursor. For each species, quantum-chemical
calculations have been carried out with the aim of obtaining accurate
predictions of the missing spectroscopic information required to guide
spectral analysis and assignment. Rotational spectra of the title
molecules have been investigated up to 400 GHz by using a frequency-modulation
millimeter-/submillimeter-wave spectrometer, thus significantly extending
spectral predictions over a wide range of frequency and quantum numbers.
A comparison between our results and those available in the literature
points out the clear need of the reported laboratory measurements
at higher frequencies for setting up accurate line catalogs for astronomical
searches.
Collapse
Affiliation(s)
- Alessio Melli
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy.,Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Mattia Melosso
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy.,Scuola Superiore Meridionale, Largo San Marcellino 10, 80138 Naples, Italy
| | - Luca Bizzocchi
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy
| | - Silvia Alessandrini
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy.,Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Ningjing Jiang
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy
| | - Francesca Tonolo
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy.,Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Salvatore Boi
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy
| | - Giorgia Castellan
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy
| | - Carlotta Sapienza
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy
| | - Jean-Claude Guillemin
- Univ Rennes, Ecole Nationale Supérieure de Chimie de Rennes, CNRS, ISCR-UMR6226, F-35000 Rennes, France
| | - Luca Dore
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy
| | - Cristina Puzzarini
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy
| |
Collapse
|
21
|
Tasinato N, Pietropolli Charmet A, Ceselin G, Salta Z, Stoppa P. In Vitro and In Silico Vibrational-Rotational Spectroscopic Characterization of the Next-Generation Refrigerant HFO-1123. J Phys Chem A 2022; 126:5328-5342. [PMID: 35930010 PMCID: PMC9393866 DOI: 10.1021/acs.jpca.2c04680] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/04/2022] [Revised: 07/25/2022] [Indexed: 11/30/2022]
Abstract
Very short-lived substances have recently been proposed as replacements for hydrofluorocarbons (HFCs), in turn being used in place of ozone-depleting substances, in refrigerant applications. In this respect, hydro-fluoro-olefins (HFOs) are attracting particular interest because, due to their reduced global warming potential, they are supposed to be environmentally friendlier. Notwithstanding this feature, they represent a new class of compounds whose spectroscopic properties and reactivity need to be characterized to allow their atmospheric monitoring and to understand their environmental fate. In the present work, the structural, vibrational, and ro-vibrational properties of trifluorothene (HFO-1123, F2C = CHF) are studied by state-of-the-art quantum chemical calculations. The equilibrium molecular structure has an expected error within 2 mÅ and 0.2° for bond lengths and angles, respectively. This represents the first step toward the computation of highly accurate rotational constants for both the ground and first excited fundamental vibrational levels, which reproduce the available experimental data well within 0.1%. Centrifugal distortion parameters and vibrational-rotational coupling terms are computed as well and used to solve some conflicting experimental results. Simulation of the vibrational transition frequencies and intensities beyond the double harmonic approximation and up to three quanta of vibrational excitation provides insights into the couplings ruling the vibrational dynamics and guides the characterization of the gas-phase infrared spectrum experimentally recorded in the range of 200-5000 cm-1. The full characterization of the IR features is completed with the experimental determination of the absorption cross sections over the 400-5000 cm-1 region from which the radiative forcing and global warming potential of HFO-1123 are derived.
Collapse
Affiliation(s)
- Nicola Tasinato
- Scuola
Normale Superiore, SMART Laboratory, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Andrea Pietropolli Charmet
- Dipartimento
di Scienze Molecolari e Nanosistemi, Università
Ca’ Foscari Venezia, Via Torino 155, I-30172 Mestre, Italy
| | - Giorgia Ceselin
- Scuola
Normale Superiore, SMART Laboratory, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Zoi Salta
- Scuola
Normale Superiore, SMART Laboratory, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Paolo Stoppa
- Dipartimento
di Scienze Molecolari e Nanosistemi, Università
Ca’ Foscari Venezia, Via Torino 155, I-30172 Mestre, Italy
| |
Collapse
|
22
|
León I, Fusè M, Alonso ER, Mata S, Mancini G, Puzzarini C, Alonso JL, Barone V. Unbiased disentanglement of conformational baths with the help of microwave spectroscopy, quantum chemistry and artificial intelligence: the puzzling case of homocysteine. J Chem Phys 2022; 157:074107. [DOI: 10.1063/5.0102841] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
An integrated experimental-computational strategy for the accurate characterization of the conformational landscape of flexible biomolecule building blocks is proposed. This is based on the combination of rotational spectroscopy with quantum-chemical computations guided by artificial intelligence tools. The first step of the strategy is the conformer search and relative stability evaluation performed by means of an evolutionary algorithm. In this step, last generation semiempirical methods are exploited together with hybrid and double-hybrid density functionals. Next, the barriers ruling the interconversion between the low-lying conformers are evaluated in order to unravel possible fast relaxation paths. The relative stabilities and spectroscopic parameters of the ``surviving' conformers are then refined using state-of-the-art composite schemes. The reliability of the computational procedure is further improved by the inclusion of vibrational and thermal effects. The final step of the strategy is the comparison between experiment and theory without any ad hoc adjustment, which allows an unbiased assignment of the spectroscopic features in terms of different conformers and their spectroscopic parameters. The proposed approach has been tested and validated for homocysteine, a highly flexible non-proteinogenic alpha-amino acid. The synergism of the integrated strategy allowed the characterization of five conformers stabilized by bifurcated N-H-O=C hydrogen bonds, together with an additional conformer involving a more conventional HNH-O hydrogen bond. The stability order estimated from the experimental intensities as well as the number and type of conformers observed in the gas phase are in full agreement with the theoretical predictions. Analogously, a good match has been found for the spectroscopic parameters.
Collapse
Affiliation(s)
- Iker León
- Universidad de Valladolid - Campus Miguel Delibes, Spain
| | | | - Elena R. Alonso
- Química Física y Química Inorgánica, Universidad de Valladolid, Spain
| | - Santiago Mata
- Grupo de Espectroscopia Molecular (GEM). Edificio Quifima. Laboratorios de Espectroscopia y Bioespectroscopia. Unidad Asociada CSIC, Parque Científico Uva, Universidad de Valladolid, Spain
| | | | | | - Jose L. Alonso
- Grupo de Espectroscopia Molecular (GEM). Edificio Quifima. Laboratorios de Espectroscopia y Bioespectroscopia, Universidad de Valladolid Departamento Química Física y Química Inorgánica, Spain
| | | |
Collapse
|
23
|
Ceselin G, Salta Z, Bloino J, Tasinato N, Barone V. Accurate Quantum Chemical Spectroscopic Characterization of Glycolic Acid: A Route Toward its Astrophysical Detection. J Phys Chem A 2022; 126:2373-2387. [PMID: 35384666 PMCID: PMC9036519 DOI: 10.1021/acs.jpca.2c01419] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2022] [Revised: 03/18/2022] [Indexed: 11/29/2022]
Abstract
The first step to shed light on the abiotic synthesis of biochemical building blocks, and their further evolution toward biological systems, is the detection of the relevant species in astronomical environments, including earthlike planets. To this end, the species of interest need to be accurately characterized from structural, energetic, and spectroscopic viewpoints. This task is particularly challenging when dealing with flexible systems, whose spectroscopic signature is ruled by the interplay of small- and large-amplitude motions (SAMs and LAMs, respectively) and is further tuned by the conformational equilibrium. In such instances, quantum chemical (QC) calculations represent an invaluable tool for assisting the interpretation of laboratory measurements or even observations. In the present work, the role of QC results is illustrated with reference to glycolic acid (CH2OHCOOH), a molecule involved in photosynthesis and plant respiration and a precursor of oxalate in humans, which has been detected in the Murchison meteorite but not yet in the interstellar medium or in planetary atmospheres. In particular, the equilibrium structure of the lowest-energy conformer is derived by employing the so-called semiexperimental approach. Then, accurate yet cost-effective QC calculations relying on composite post-Hartree-Fock schemes and hybrid coupled-cluster/density functional theory approaches are used to predict the structural and ro-vibrational spectroscopic properties of the different conformers within the framework of the second-order vibrational perturbation theory. A purposely tailored discrete variable representation anharmonic approach is used to treat the LAMs related to internal rotations. The computed spectroscopic data, particularly those in the infrared region, complement the available experimental investigations, thus enhancing the possibility of an astronomical detection of this molecule.
Collapse
Affiliation(s)
- Giorgia Ceselin
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126, Pisa, Italy
| | - Zoi Salta
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126, Pisa, Italy
| | - Julien Bloino
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126, Pisa, Italy
| | - Nicola Tasinato
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126, Pisa, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126, Pisa, Italy
| |
Collapse
|
24
|
Dékány AÁ, Czakó G. Benchmark ab initio proton affinity and gas-phase basicity of α-alanine based on coupled-cluster theory and statistical mechanics. J Comput Chem 2022; 43:19-28. [PMID: 34676890 DOI: 10.1002/jcc.26767] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2021] [Revised: 09/16/2021] [Accepted: 09/27/2021] [Indexed: 11/06/2022]
Abstract
We determine the proton affinity (PA) and gas-phase basicity (GB) of amino acid α-alanine at a chemically accurate level by performing explicitly-correlated CCSD(T)-F12b/aug-cc-pVDZ geometry optimizations and normal mode vibrational frequency calculations as well as CCSD(T)-F12b/aug-cc-pVTZ energy computations at the possible neutral and protonated geometries. Temperature effects at 298.15 K considering translational, rotational, and vibrational enthalpy and entropy corrections are obtained via standard statistical mechanics utilizing the molecular geometries and the harmonic vibrational energy levels. Both the amino nitrogen (N) and the carbonyl oxygen (O) atoms are proven to be potential protonation sites and a systematic conformational search reveals 3 N- and 9 O-protonated conformers in the 0.00-7.88 and 25.43-30.43 kcal/mol energy ranges at 0 K, respectively. The final computed PA and GB values at (0)298.15 K in case of N-protonation are (214.47)216.80 and 207.07 kcal/mol, respectively, whereas the corresponding values for O-protonation are (189.04)190.63 and 182.31 kcal/mol. The results of the benchmark high-level coupled-cluster computations are utilized to assess the accuracy of several lower-level cost-effective methods such as MP2 and density functional theory with various functionals (SOGGA11-X, M06-2X, PBE0, B3LYP, M06, TPSS).
Collapse
Affiliation(s)
- Attila Á Dékány
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Szeged, Hungary
| | - Gábor Czakó
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Szeged, Hungary
| |
Collapse
|
25
|
Wang P, Shu C, Ye H, Biczysko M. Structural and Energetic Properties of Amino Acids and Peptides Benchmarked by Accurate Theoretical and Experimental Data. J Phys Chem A 2021; 125:9826-9837. [PMID: 34752094 DOI: 10.1021/acs.jpca.1c06504] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
Structural, energetic, and spectroscopic data derived in this work aim at the setup of an "experimentally validated" database for amino acids and polypeptides conformers. First, the "cheap" composite scheme (ChS, CCSD(T)/(CBS+CV)MP2) is tested for evaluation of conformational energies of all eight stable conformers of glycine, by comparing to the more accurate CCSD(T)/CBS+CV computations (Phys. Chem. Chem. Phys. 2013, 15, 10094-10111 and J Mol. Model. 2020, 26, 129). The recently proposed jun-ChS (J. Chem. Theory and Comput. 2020, 16, 988-1006), employing the jun-cc-pVnZ basis set family for CCSD(T) computations and CBS extrapolation, yields conformational energies accurate to 0.2 kJ·mol-1, at reduced computational cost with respect to aug-ChS employing aug-cc-pVnZ basis sets. The jun-ChS composite scheme is further applied to derive conformational energies for three dipeptide analogues Ac-Gly-NH2, Ac-Ala-NH2, and Gly-Gly. Finally, dipeptide conformational energies and semiexperimental equilibrium rotational constants along with the CCSD(T)/(CBS+CV)MP2 structural parameters (J. Phys. Chem. Lett. 2014, 5, 534-540) stand as the reference for benchmarking of selected density functional methodologies. The double-hybrid functionals B2-PLYP-D3(BJ) and DSD-PBEP86, perform best for structural and energetic characterization of all dipeptide analogues. From hybrid functionals CAM-B3LYP-D3(BJ) and ωB97X-D3(BJ) represent promising methods applicable for larger peptide-based systems for which computations with double-hybrid functionals are not feasible.
Collapse
Affiliation(s)
- Ping Wang
- International Centre for Quantum and Molecular Structures, Physics Department, College of Science, Shanghai University, 99 Shangda Road, Shanghai 200444, China
| | - Chong Shu
- International Centre for Quantum and Molecular Structures, Physics Department, College of Science, Shanghai University, 99 Shangda Road, Shanghai 200444, China
| | - Hexu Ye
- International Centre for Quantum and Molecular Structures, Physics Department, College of Science, Shanghai University, 99 Shangda Road, Shanghai 200444, China
| | - Malgorzata Biczysko
- International Centre for Quantum and Molecular Structures, Physics Department, College of Science, Shanghai University, 99 Shangda Road, Shanghai 200444, China
| |
Collapse
|
26
|
Lupi J, Alessandrini S, Puzzarini C, Barone V. junChS and junChS-F12 Models: Parameter-free Efficient yet Accurate Composite Schemes for Energies and Structures of Noncovalent Complexes. J Chem Theory Comput 2021; 17:6974-6992. [PMID: 34677974 DOI: 10.1021/acs.jctc.1c00869] [Citation(s) in RCA: 25] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A recently developed model chemistry (denoted as junChS [Alessandrini, S.; et al. J. Chem. Theory Comput. 2020, 16, 988-1006]) has been extended to the employment of explicitly correlated (F12) methods. This led us to propose a family of effective, reliable, and parameter-free schemes for the computation of accurate interaction energies of molecular complexes ruled by noncovalent interactions. A thorough benchmark based on a wide range of interactions showed that the so-called junChS-F12 model, which employs cost-effective revDSD-PBEP86-D3(BJ) reference geometries, has an improved performance with respect to its conventional counterpart and outperforms well-known model chemistries. Without employing any empirical parameter and at an affordable computational cost, junChS-F12 reaches subchemical accuracy. Accurate characterizations of molecular complexes are usually limited to energetics. To take a step forward, the conventional and F12 composite schemes developed for interaction energies have been extended to structural determinations. A benchmark study demonstrated that the most effective option is to add MP2-F12 core-valence correlation corrections to fc-CCSD(T)-F12/jun-cc-pVTZ geometries without the need of recovering the basis set superposition error and the extrapolation to the complete basis set.
Collapse
Affiliation(s)
- Jacopo Lupi
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Silvia Alessandrini
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy.,Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, I-40126 Bologna, Italy
| | - Cristina Puzzarini
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, I-40126 Bologna, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| |
Collapse
|
27
|
Sheng M, Silvestrini F, Biczysko M, Puzzarini C. Structural and Vibrational Properties of Amino Acids from Composite Schemes and Double-Hybrid DFT: Hydrogen Bonding in Serine as a Test Case. J Phys Chem A 2021; 125:9099-9114. [PMID: 34623165 DOI: 10.1021/acs.jpca.1c06993] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
The structures, relative stabilities, and vibrational wavenumbers of the two most stable conformers of serine, stabilized by the O-H···N, O-H···O═C and N-H···O-H intramolecular hydrogen bonds, have been evaluated by means of state-of-the-art composite schemes based on coupled-cluster (CC) theory. The so-called "cheap" composite approach (CCSD(T)/(CBS+CV)MP2) allowed determination of accurate equilibrium structures and harmonic vibrational wavenumbers, also pointing out significant corrections beyond the CCSD(T)/cc-pVTZ level. These accurate results stand as a reference for benchmarking selected hybrid and double-hybrid, dispersion-corrected DFT functionals. B2PLYP-D3 and DSDPBEP86 in conjunction with a triple-ζ basis set have been confirmed as effective methodologies for structural and spectroscopic studies of medium-sized flexible biomolecules, also showing intramolecular hydrogen bonding. These best performing double-hybrid functionals have been employed to simulate IR spectra by means of vibrational perturbation theory, also considering hybrid CC/DFT schemes. The best overall agreement with experiment, with mean absolute error of 8 cm-1, has been obtained by combining CCSD(T)/(CBS+CV)MP2 harmonic wavenumbers with B2PLYP-D3/maug-cc-pVTZ anharmonic corrections. Finally, a composite scheme entirely based on CCSD(T) calculations (CCSD(T)/CBS+CV) has been employed for energetics, further confirming that serine II is the most stable conformer, also when zero-point vibrational energy corrections are included.
Collapse
Affiliation(s)
- Mingzhu Sheng
- International Centre for Quantum and Molecular Structures, Physics Department, Shanghai University, 99 Shangda Road, Shanghai, 200444, China
| | - Filippo Silvestrini
- Department of Chemistry "Giacomo Ciamician", University of Bologna, Via F. Selmi 2, 40126 Bologna, Italy
| | - Malgorzata Biczysko
- International Centre for Quantum and Molecular Structures, Physics Department, Shanghai University, 99 Shangda Road, Shanghai, 200444, China
| | - Cristina Puzzarini
- Department of Chemistry "Giacomo Ciamician", University of Bologna, Via F. Selmi 2, 40126 Bologna, Italy
| |
Collapse
|
28
|
Ozaki Y, Beć KB, Morisawa Y, Yamamoto S, Tanabe I, Huck CW, Hofer TS. Advances, challenges and perspectives of quantum chemical approaches in molecular spectroscopy of the condensed phase. Chem Soc Rev 2021; 50:10917-10954. [PMID: 34382961 DOI: 10.1039/d0cs01602k] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The purpose of this review is to demonstrate advances, challenges and perspectives of quantum chemical approaches in molecular spectroscopy of the condensed phase. Molecular spectroscopy, particularly vibrational spectroscopy and electronic spectroscopy, has been used extensively for a wide range of areas of chemical sciences and materials science as well as nano- and biosciences because it provides valuable information about structure, functions, and reactions of molecules. In the meantime, quantum chemical approaches play crucial roles in the spectral analysis. They also yield important knowledge about molecular and electronic structures as well as electronic transitions. The combination of spectroscopic approaches and quantum chemical calculations is a powerful tool for science, in general. Thus, our article, which treats various spectroscopy and quantum chemical approaches, should have strong implications in the wider scientific community. This review covers a wide area of molecular spectroscopy from far-ultraviolet (FUV, 120-200 nm) to far-infrared (FIR, 400-10 cm-1)/terahertz and Raman spectroscopy. As quantum chemical approaches, we introduce several anharmonic approaches such as vibrational self-consistent field (VSCF) and the combination of periodic harmonic calculations with anharmonic corrections based on finite models, grid-based techniques like the Numerov approach, the Cartesian coordinate tensor transfer (CCT) method, Symmetry-Adapted Cluster Configuration-Interaction (SAC-CI), and the ZINDO (Semi-empirical calculations at Zerner's Intermediate Neglect of Differential Overlap). One can use anharmonic approaches and grid-based approaches for both infrared (IR) and near-infrared (NIR) spectroscopy, while CCT methods are employed for Raman, Raman optical activity (ROA), FIR/terahertz and low-frequency Raman spectroscopy. Therefore, this review overviews cross relations between molecular spectroscopy and quantum chemical approaches, and provides various kinds of close-reality advanced spectral simulation for condensed phases.
Collapse
Affiliation(s)
- Yukihiro Ozaki
- School of Biological and Environmental Sciences, Kwansei Gakuin University, Sanda, Hyogo 669-1337, Japan. and Toyota Physical and Chemical Research Institute, Yokomichi, Nagakute, Aichi 480-1192, Japan
| | - Krzysztof B Beć
- Institute of Analytical Chemistry and Radiochemistry, University of Innsbruck, Innrain 80/82, 6020 Innsbruck, Austria
| | - Yusuke Morisawa
- Department of Chemistry, School of Science and Engineering, Kindai University, Kowakae, Higashi-Osaka, Osaka 577-8502, Japan
| | - Shigeki Yamamoto
- Department of Chemistry, Graduate School of Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Ichiro Tanabe
- Department of Materials Engineering Science, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Christian W Huck
- Institute of Analytical Chemistry and Radiochemistry, University of Innsbruck, Innrain 80/82, 6020 Innsbruck, Austria
| | - Thomas S Hofer
- Institute of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innrain 80-82, A6020 Innsbruck, Austria
| |
Collapse
|
29
|
Rani N. Rotational and Vibrational Signatures of Astrophyscially relevant Gas-Phase Stereo-isomeric Species of Proteinogenic Amino acid Leucine. LIFE SCIENCES IN SPACE RESEARCH 2021; 30:29-38. [PMID: 34281662 DOI: 10.1016/j.lssr.2021.05.003] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/17/2021] [Revised: 05/16/2021] [Accepted: 05/17/2021] [Indexed: 06/13/2023]
Abstract
The search for life-supporting molecules in outer space is an ever-growing endeavour. Towards this, the computational chemistry supporting the astronomical spectroscopic observations is becoming a valuable tool to unravel the complex chemical network in interstellar medium (ISM). In the present work, quantum-mechanical computations, accounting for anharmonic effects, are performed to obtain the rotational and vibrational line-data for the gas-phase conformers of proteinogenic amino acid Leucine and its isomeric species predicted to be involved in its stereoinversion under the extreme environment of ISM. These species exhibit diverse chemistry including branched skeleton and zwitterionic ammonium ylides. A few of the species have significantly high dipole moment, which can act as tracer for the conformers of Leucine having low dipole moment. Besides this, the species, which are terrestrially less stable, can be of significant importance to the astronomers. Notably, the spectral database generated in this work can assist in the detection of proteinogenic Leucine and its isomeric species in different regions of ISM.
Collapse
Affiliation(s)
- Namrata Rani
- Quantum Chemistry Group, Department of Chemistry & Centre of Advanced Studies in Chemistry, Panjab University, Chandigarh-160014, India
| |
Collapse
|
30
|
Abstract
A systematic conformational search reveals three N- (amino) and eight O- (carbonyl) protonated glycine conformers with benchmark equilibrium(adiabatic) relative energies in the 0.00-7.51(0.00-7.37) and 25.91-31.61(24.45-30.28) kcal mol-1 ranges, respectively. Benchmark ab initio structures of the glycine conformers and its protonated species are obtained at the CCSD(T)-F12b/aug-cc-pVTZ level of theory and the relative energy computations consider basis-set effects up to aug-cc-pVQZ with CCSD(T)-F12b, electron correlation up to CCSDT(Q), core correlation corrections, scalar relativistic effects, and zero-point energy contributions. The best predictions for Boltzmann-averaged 0(298.15) K proton affinities and [298.15 K gas-phase basicities] of glycine are 211.00(212.43)[204.75] and 186.38(187.64)[180.21] kcal mol-1 for N- and O-protonation, respectively, in excellent agreement with experiments.
Collapse
Affiliation(s)
- András B Nacsa
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary.
| | | |
Collapse
|
31
|
Qu C, Houston PL, Conte R, Nandi A, Bowman JM. MULTIMODE Calculations of Vibrational Spectroscopy and 1d Interconformer Tunneling Dynamics in Glycine Using a Full-Dimensional Potential Energy Surface. J Phys Chem A 2021; 125:5346-5354. [PMID: 34110169 DOI: 10.1021/acs.jpca.1c03738] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A full-dimensional, permutationally invariant polynomial potential energy surface for glycine recently reported (R. Conte et al., J. Chem. Phys. 2020, 153, 244301) is used with the code MULTIMODE to determine the IR absorption spectra for Conformers I and II using a new separable dipole moment function. The calculated spectra agree well with the experimental ones. The full-dimensional nature of the potential allows us also to examine dynamical results, such as tunneling rates. Remarkably, using a one-dimensional path based on the potential energy surface to estimate the tunneling rate from Conformer VI to Conformer I, good agreement is found with the recent experimental measurement. Finally a brief comparison of our potential energy surface with a recently reported sGDML one is made.
Collapse
Affiliation(s)
- Chen Qu
- Department of Chemistry & Biochemistry, University of Maryland, College Park, Maryland 20742, United States
| | - Paul L Houston
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States.,Department of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | - Riccardo Conte
- Dipartimento di Chimica, Università degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| | - Apurba Nandi
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Joel M Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
32
|
Barone V, Alessandrini S, Biczysko M, Cheeseman JR, Clary DC, McCoy AB, DiRisio RJ, Neese F, Melosso M, Puzzarini C. Computational molecular spectroscopy. ACTA ACUST UNITED AC 2021. [DOI: 10.1038/s43586-021-00034-1] [Citation(s) in RCA: 23] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
|
33
|
Barone V, Puzzarini C. Looking for the bricks of the life in the interstellar medium: The fascinating world of astrochemistry. EPJ WEB OF CONFERENCES 2020. [DOI: 10.1051/epjconf/202024600021] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
The discovery in the interstellar medium of molecules showing a certain degree of complexity, and in particular those with a prebiotic character, has attracted great interest. A complex chemistry takes place in space, but the processes that lead to the production of molecular species are a matter of intense discussion, the knowledge still being at a rather primitive stage. Debate on the origins of interstellar molecules has been further stimulated by the identification of biomolecular building blocks, such as nucleobases and amino acids, in meteorites and comets. Since many of the molecules found in space play a role in the chemistry of life, the issue of their molecular genesis and evolution might be related to the profound question of the origin of life itself. Understanding the underlying chemical processes, including the production, reactions and destruction of compounds, requires the concomitant study of spectroscopy, gas-phase reactivity, and heterogeneous processes on dust-grains. The aim of this contribution is to provide a general view of a complex and multifaceted challenge, while focusing on the role played by molecular spectroscopy and quantum-chemical computations. In particular, the derivation of the molecular spectroscopic features and the investigation of gas-phase formation routes of prebiotic species in the interstellar medium are addressed from a computational point of view.
Collapse
|
34
|
Gawrilow M, Suhm MA. 2-Methoxyethanol: harmonic tricks, anharmonic challenges and chirality-sensitive chain aggregation. Phys Chem Chem Phys 2020; 22:15303-15311. [PMID: 32626860 DOI: 10.1039/d0cp02488k] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The Raman spectrum of jet-cooled, nearly monoconformational 2-methoxyethanol reveals a strong, aggregation-sensitive resonance in the backbone stretching region and pronounced chirality-dependent dimerisation effects in the OH stretching region. These effects are disentangled by a dual-detection scheme which simultaneously probes the two spectral regions. The performance of mass-scaled harmonic wavenumber predictions and of vibrational perturbation theory in reproducing these spectral features is explored.
Collapse
Affiliation(s)
- Maxim Gawrilow
- Institut für Physikalische Chemie, Georg-August-Universität Göttingen, Tammannstr. 6, 37077 Göttingen, Germany.
| | | |
Collapse
|
35
|
Ye H, Mendolicchio M, Kruse H, Puzzarini C, Biczysko M, Barone V. The challenging equilibrium structure of HSSH: Another success of the rotational spectroscopy / quantum chemistry synergism. J Mol Struct 2020. [DOI: 10.1016/j.molstruc.2020.127933] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
36
|
Orján EM, Nacsa AB, Czakó G. Conformers of dehydrogenated glycine isomers. J Comput Chem 2020; 41:2001-2014. [DOI: 10.1002/jcc.26375] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/24/2020] [Revised: 06/05/2020] [Accepted: 06/10/2020] [Indexed: 12/20/2022]
Affiliation(s)
- Erik M. Orján
- MTA‐SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials ScienceInstitute of Chemistry, University of Szeged Szeged Hungary
| | - András B. Nacsa
- MTA‐SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials ScienceInstitute of Chemistry, University of Szeged Szeged Hungary
| | - Gábor Czakó
- MTA‐SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials ScienceInstitute of Chemistry, University of Szeged Szeged Hungary
| |
Collapse
|
37
|
Puzzarini C, Spada L, Alessandrini S, Barone V. The challenge of non-covalent interactions: theory meets experiment for reconciling accuracy and interpretation. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2020; 32:343002. [PMID: 32203942 DOI: 10.1088/1361-648x/ab8253] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/27/2019] [Accepted: 03/23/2020] [Indexed: 06/10/2023]
Abstract
In the past decade, many gas-phase spectroscopic investigations have focused on the understanding of the nature of weak interactions in model systems. Despite the fact that non-covalent interactions play a key role in several biological and technological processes, their characterization and interpretation are still far from being satisfactory. In this connection, integrated experimental and computational investigations can play an invaluable role. Indeed, a number of different issues relevant to unraveling the properties of bulk or solvated systems can be addressed from experimental investigations on molecular complexes. Focusing on the interaction of biological model systems with solvent molecules (e.g., water), since the hydration of the biomolecules controls their structure and mechanism of action, the study of the molecular properties of hydrated systems containing a limited number of water molecules (microsolvation) is the basis for understanding the solvation process and how structure and reactivity vary from gas phase to solution. Although hydrogen bonding is probably the most widespread interaction in nature, other emerging classes, such as halogen, chalcogen and pnicogen interactions, have attracted much attention because of the role they play in different fields. Their understanding requires, first of all, the characterization of the directionality, strength, and nature of such interactions as well as a comprehensive analysis of their competition with other non-covalent bonds. In this review, it is shown how state-of-the-art quantum-chemical computations combined with rotational spectroscopy allow for fully characterizing intermolecular interactions taking place in molecular complexes from both structural and energetic points of view. The transition from bi-molecular complex to microsolvation and then to condensed phase is shortly addressed.
Collapse
Affiliation(s)
- Cristina Puzzarini
- Dipartimento di Chimica 'Giacomo Ciamician', Via F. Selmi 2, I-40126 Bologna, Italy
| | - Lorenzo Spada
- Dipartimento di Chimica 'Giacomo Ciamician', Via F. Selmi 2, I-40126 Bologna, Italy
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Silvia Alessandrini
- Dipartimento di Chimica 'Giacomo Ciamician', Via F. Selmi 2, I-40126 Bologna, Italy
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| |
Collapse
|
38
|
Toward accurate prediction of amino acid derivatives structure and energetics from DFT: glycine conformers and their interconversions. J Mol Model 2020; 26:129. [DOI: 10.1007/s00894-020-4342-7] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2019] [Accepted: 03/04/2020] [Indexed: 12/30/2022]
|
39
|
Experimental and computational analysis of 1-(4-chloro-3-nitrophenyl)-3-(3,4-dichlorophenyl)thiourea. J Mol Struct 2020. [DOI: 10.1016/j.molstruc.2019.127587] [Citation(s) in RCA: 34] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
|
40
|
A never-ending story in the sky: The secrets of chemical evolution. Phys Life Rev 2020; 32:59-94. [DOI: 10.1016/j.plrev.2019.07.001] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2019] [Revised: 06/26/2019] [Accepted: 07/02/2019] [Indexed: 01/13/2023]
|
41
|
Lee KLK, McCarthy M. Bayesian Analysis of Theoretical Rotational Constants from Low-Cost Electronic Structure Methods. J Phys Chem A 2020; 124:898-910. [DOI: 10.1021/acs.jpca.9b09982] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Kin Long Kelvin Lee
- Harvard & Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, Massachusetts 02138, United States
| | - Michael McCarthy
- Harvard & Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, Massachusetts 02138, United States
| |
Collapse
|
42
|
Puzzarini C, Barone V. The challenging playground of astrochemistry: an integrated rotational spectroscopy - quantum chemistry strategy. Phys Chem Chem Phys 2020; 22:6507-6523. [PMID: 32163090 DOI: 10.1039/d0cp00561d] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
While it is now well demonstrated that the interstellar medium (ISM) is characterized by a diverse and complex chemistry, a significant number of features in radioastronomical spectra are still unassigned and call for new laboratory efforts, which are increasingly based on integrated experimental and computational strategies. In parallel, the identification of an increasing number of molecules containing more than five atoms and at least one carbon atom (the so-called "interstellar" complex organic molecules), which can play a relevant role in the chemistry of life, raises the additional issue of how these species can be produced in the typical harsh conditions of the ISM. On these grounds, this perspective aims to present an integrated rotational spectroscopy - quantum chemistry approach for supporting radioastronomical observations and a computational strategy for contributing to the elucidation of chemical reactivity in the interstellar space.
Collapse
Affiliation(s)
- Cristina Puzzarini
- Dipartimento di Chimica "Giacomo Ciamician", University of Bologna, via F. Selmi 2, I-40126 Bologna, Italy.
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, Pisa, I-56126, Italy
| |
Collapse
|
43
|
Alessandrini S, Barone V, Puzzarini C. Extension of the “Cheap” Composite Approach to Noncovalent Interactions: The jun-ChS Scheme. J Chem Theory Comput 2019; 16:988-1006. [DOI: 10.1021/acs.jctc.9b01037] [Citation(s) in RCA: 39] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Silvia Alessandrini
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
- Dipartimento di Chimica “Giacomo Ciamician”, Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Cristina Puzzarini
- Dipartimento di Chimica “Giacomo Ciamician”, Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy
| |
Collapse
|
44
|
Myhre RH, Coriani S, Koch H. X-ray and UV Spectra of Glycine within Coupled Cluster Linear Response Theory. J Phys Chem A 2019; 123:9701-9711. [PMID: 31549830 DOI: 10.1021/acs.jpca.9b06590] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The coupled cluster models CCSD and CC3 are used to investigate the (core) excited states and ionization energies of glycine in the gas phase. Excited states and ionization energies in the UV spectral range are calculated using a standard coupled cluster linear response, while core-level excited states and ionization potentials are calculated using the core-valence separation approximation. The temperature dependence from different conformers is also assessed.
Collapse
Affiliation(s)
- Rolf H Myhre
- Department of Chemistry , Norwegian University of Science and Technology, NTNU , 7491 Trondheim , Norway.,Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry , University of Oslo , 0315 Oslo , Norway
| | - Sonia Coriani
- DTU Chemistry , Technical University of Denmark , DK-2800 Kongens Lyngby , Denmark.,Aarhus Institute of Advanced Studies , Aarhus University , DK-8000 Aarhus C , Denmark
| | - Henrik Koch
- Department of Chemistry , Norwegian University of Science and Technology, NTNU , 7491 Trondheim , Norway.,Scuola Normale Superiore , Piazza dei Cavalieri 7 , 56126 Pisa , Italy
| |
Collapse
|
45
|
The most stable isomer of H2C4-(OCS)2 van der Waals complex: Theory and experiment agree on a structure with C2 symmetry. Chem Phys Lett 2019. [DOI: 10.1016/j.cplett.2019.136610] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
|
46
|
Saladino R. From chemical complexity to origin of life: Comment on "A never-ending story in the sky: The secrets of chemical evolution" by C. Puzzarini and V. Barone. Phys Life Rev 2019; 32:111-113. [PMID: 31451434 DOI: 10.1016/j.plrev.2019.08.014] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2019] [Accepted: 08/19/2019] [Indexed: 10/26/2022]
Affiliation(s)
- Raffaele Saladino
- Department of Ecological and Biological Sciences, Via S. Camillo de Lellis, University of Tuscia, 01100, Viterbo, Italy.
| |
Collapse
|
47
|
Puzzarini C, Bloino J, Tasinato N, Barone V. Accuracy and Interpretability: The Devil and the Holy Grail. New Routes across Old Boundaries in Computational Spectroscopy. Chem Rev 2019; 119:8131-8191. [DOI: 10.1021/acs.chemrev.9b00007] [Citation(s) in RCA: 114] [Impact Index Per Article: 22.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Affiliation(s)
- Cristina Puzzarini
- Dipartimento di Chimica “Giacomo Ciamician”, Università di Bologna, Via F. Selmi 2, I-40126 Bologna, Italy
| | - Julien Bloino
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Nicola Tasinato
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| |
Collapse
|
48
|
Zhang H, Krupa J, Wierzejewska M, Biczysko M. The role of dispersion and anharmonic corrections in conformational analysis of flexible molecules: the allyl group rotamerization of matrix isolated safrole. Phys Chem Chem Phys 2019; 21:8352-8364. [PMID: 30958495 DOI: 10.1039/c9cp00926d] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Conformational changes of the monomeric safrole (5-(2-propenyl)-1,3-benzodioxole) isolated in low temperature xenon matrices were induced thermally or using narrow-band UV radiation. The rotation of the allyl group taking place in the studied matrices was followed by FTIR spectroscopy. Safrole represents a challenging example of a flexible molecule highlighting the importance of dispersion interactions and anharmonic effects in the structural, spectroscopic and energetic analysis. Structures of the safrole conformers, their energetics and infrared spectra have been calculated using various computational methods ranging from density functional theory (DFT) to coupled cluster (CC). The best theoretical results were obtained by integrating CCSD(T) energies including complete basis set extrapolation and core-valence corrections with B2PLYP-D3 equilibrium structures and hybrid B2PLYP-D3/B3LYP-D3 anharmonic computations for IR spectra and thermodynamics.
Collapse
Affiliation(s)
- Hongli Zhang
- International Centre for Quantum and Molecular Structures, Department of Physics, College of Sciences, Shanghai University, 99 Shangda Road, Shanghai, 200444, China.
| | | | | | | |
Collapse
|
49
|
Puzzarini C, Tasinato N, Bloino J, Spada L, Barone V. State-of-the-art computation of the rotational and IR spectra of the methyl-cyclopropyl cation: hints on its detection in space. Phys Chem Chem Phys 2019; 21:3431-3439. [PMID: 30110028 DOI: 10.1039/c8cp04629h] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Recent measurements by the Cassini Ion Neutral Mass Spectrometer demonstrated the presence of numerous carbocations in Titan's upper atmosphere. In [Ali et al., Planet. Space Sci., 2013, 87, 96], an analysis of these measurements revealed the formation of the three-membered cyclopropenyl cation and its methyl derivatives. As a starting point of a future coordinated effort of laboratory experiments, quantum-chemical calculations, and astronomical observations, in the present work the molecular structure and spectroscopic properties of the methyl-cyclopropenyl cation have been investigated by means of state-of-the-art computational approaches in order to simulate its rotational and infrared spectra. Rotational parameters have been predicted with an expected accuracy better than 0.1% for rotational constants and on the order of 1-2% for centrifugal-distortion terms. As for the infrared spectrum, despite the challenge of a large amplitude motion, fundamental transitions have been computed to a good accuracy, i.e., the uncertainties are expected to be smaller than 5-10 wavenumbers.
Collapse
Affiliation(s)
- Cristina Puzzarini
- Dipartimento di Chimica "Giacomo Ciamician", University of Bologna, via F. Selmi 2, I-40126 Bologna, Italy.
| | | | | | | | | |
Collapse
|
50
|
Chandramouli B, Del Galdo S, Fusè M, Barone V, Mancini G. Two-level stochastic search of low-energy conformers for molecular spectroscopy: implementation and validation of MM and QM models. Phys Chem Chem Phys 2019; 21:19921-19934. [DOI: 10.1039/c9cp03557e] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
The search for stationary points in the molecular potential energy surfaces (PES) is a problem of increasing relevance in molecular sciences especially for large, flexible systems featuring several large-amplitude internal motions.
Collapse
Affiliation(s)
| | | | | | - Vincenzo Barone
- Scuola Normale Superiore
- 56126 Pisa
- Italy
- Istituto Nazionale di Fisica Nucleare (INFN)
- Sezione di Pisa
| | - Giordano Mancini
- Scuola Normale Superiore
- 56126 Pisa
- Italy
- Istituto Nazionale di Fisica Nucleare (INFN)
- Sezione di Pisa
| |
Collapse
|