1
|
Beć KB, Grabska J, Hawranek JP, Huck CW. Carbonyl stretching band in amides as Lorentz oscillator. Insights into anharmonicity and local environment in the liquid phase from NIR and MIR spectra of N-methylformamide and di-N,N-methylformamide. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2025; 324:124954. [PMID: 39180970 DOI: 10.1016/j.saa.2024.124954] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/08/2024] [Revised: 08/02/2024] [Accepted: 08/08/2024] [Indexed: 08/27/2024]
Abstract
We investigated the anharmonicity and intermolecular interactions of N-methylformamide (NMF) and di-N,N-methylformamide (DMF) in the neat liquid phase with particular interest in the amide bands. The vibrational spectra, complex refractive index, and complex electric permittivity were determined in in the mid- (MIR) and near-infrared (NIR) regions (11,500-560 cm-1; 870-17857 nm). Dispersion analysis was based on the Classical Damped Harmonic Oscillator (CDHO) and simultaneous modelling of the real and imaginary components of the spectra. This data delivered insights into the vibrational energy dissipation and self-association in liquid amides. Identification of the MIR and NIR bands was based on anharmonic GVPT2//B3LYP/6-311++G(d,p) calculations. DMF and NMF follow distinct self-association, evidenced in the MIR fingerprint by the two components of the νCO, the analog of the Amide I band. These conclusions are supported by the structural information derived from the NIR spectra. Furthermore, the contribution of overtones and combination bands in the MIR spectra of amides was examined. The conclusions on molecular interactions and structural dynamics of NMF and DMF contribute to a deeper understanding of the effects of changes in the local environment of the amide group.
Collapse
Affiliation(s)
- Krzysztof B Beć
- Institute of Analytical Chemistry and Radiochemistry, University of Innsbruck, Innsbruck, Austria
| | - Justyna Grabska
- Institute of Analytical Chemistry and Radiochemistry, University of Innsbruck, Innsbruck, Austria.
| | | | - Christian W Huck
- Institute of Analytical Chemistry and Radiochemistry, University of Innsbruck, Innsbruck, Austria
| |
Collapse
|
2
|
Khanifaev J, Schrader T, Perlt E. Machine-learning to predict anharmonic frequencies: a study of models and transferability. Phys Chem Chem Phys 2024; 26:23495-23502. [PMID: 39222042 DOI: 10.1039/d4cp01789g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/04/2024]
Abstract
With more and more accurate electronic structure methods at hand, the inclusion of anharmonic effects in the post-processing of such data towards thermochemical properties is the next step. In this context, the description of anharmonicity has been an important topic of physical chemistry and chemical physics for a long time. In this study, anharmonic frequencies of various hydrogen-halides and halogenated hydrocarbon molecular clusters are calculated using harmonic as well as explicitly anharmonic methods, i.e., normal mode analysis and vibrational self-consistent field. Simple harmonic model based descriptors were used to predict anharmonic frequencies via multilinear regression and gradient boosting regression. Gradient boosting regression is capable of predicting reliable anharmonic data and even the simple multilinear regression model yields reasonable predictions that can account for mode-to-mode couplings. Moreover, the transferability to unseen chemical systems is assessed and it is confirmed that the machine-learned models can be applied to larger, unseen molecules.
Collapse
Affiliation(s)
| | - Tim Schrader
- Friedrich Schiller University Jena, Löbdergraben 32, 07743 Jena, Germany.
| | - Eva Perlt
- Friedrich Schiller University Jena, Löbdergraben 32, 07743 Jena, Germany.
| |
Collapse
|
3
|
Mendolicchio M, Barone V. Accurate Vibrational and Ro-Vibrational Contributions to the Properties of Large Molecules by a New Engine Employing Curvilinear Internal Coordinates and Vibrational Perturbation Theory to Second Order. J Chem Theory Comput 2024. [PMID: 39215708 DOI: 10.1021/acs.jctc.4c00857] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/04/2024]
Abstract
The unbiased comparison between theory and experiment requires approaches more sophisticated than the basic harmonic-oscillator rigid-rotor model, for taking into account vibrational averaging effects and ro-vibrational couplings in molecules of increasing size. Second-order vibrational perturbation theory based on curvilinear internal coordinates (ICs) offers a remarkable compromise between accuracy and computational cost, thanks to the reduction of mode-mode couplings with respect to their counterparts based on Cartesian coordinates. Therefore, we have developed, implemented, and validated a general engine employing ICs, which allows the accurate evaluation of vibrational averages and ro-vibrational couplings for molecules containing up to about 50 atoms beyond the harmonic approximation. After validation of the new tool for relatively small molecules, the effectiveness of ICs has been demonstrated for some flexible and/or quite large molecular bricks of life.
Collapse
|
4
|
Bloino J, Jähnigen S, Merten C. After 50 Years of Vibrational Circular Dichroism Spectroscopy: Challenges and Opportunities of Increasingly Accurate and Complex Experiments and Computations. J Phys Chem Lett 2024; 15:8813-8828. [PMID: 39167088 DOI: 10.1021/acs.jpclett.4c01700] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/23/2024]
Abstract
VCD research continues to thrive, driven by ongoing experimental and theoretical advances. Modern studies deal with increasingly complex samples featuring weak intermolecular interactions and shallow potential energy surfaces. Likewise, the combination of VCD measurements with, for instance, cryo-spectroscopic techniques has significantly increased their sensitivity. The extent to which such modern measurements enhance the informative value of VCD depends significantly on the quality of the theoretical models, which must adequately account for anharmonicity, solvation and molecular dynamics. We herein discuss how experimental advancements engage in a stimulating interplay with recent theoretical developments, pursuing either the static or the dynamic computational route. Both paths have their own strengths and limitations, each addressing fundamentally different problems. We give an outlook on future challenges of VCD research, including the possibility to combine static and dynamic approaches to obtain a full picture of the sample.
Collapse
Affiliation(s)
- Julien Bloino
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Sascha Jähnigen
- Freie Universität Berlin, Institut für Chemie und Biochemie, Arnimallee 22, 14195 Berlin, Germany
| | - Christian Merten
- Ruhr-Universität Bochum, Fakultät für Chemie und Biochemie, Universitätsstraße 150, 44801 Bochum, Germany
| |
Collapse
|
5
|
Díaz Soto LJ, Oliveira RR, Baptista L, da Silveira EF, Nascimento MAC. Energy and spectroscopic parameters of neutral and cations isomers of the C nH 2 (n = 2-6) families using high-level ab-initio approaches. J Comput Chem 2024. [PMID: 39177429 DOI: 10.1002/jcc.27485] [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: 05/10/2024] [Revised: 07/27/2024] [Accepted: 08/02/2024] [Indexed: 08/24/2024]
Abstract
Cationic species, previously detected from ion-induced desorption of solid methane by plasma desorption mass spectrometry (PDMS), and neutral species, are investigated using high-level ab-initio approaches. From a set of 25 cationic and 26 neutral structures belonging to CnH2 (n = 2-6) families, it was obtained the energy, rotational constants, harmonic vibrational frequency, charge distribution and excitation energies. The ZPVE-corrected energies, at CCSD(T)-F12; CCSD(T)-F12/RI/(cc-pVTZ-F12, cc-pVTZ-F12-CABS, cc-pVQZ/C) (n = 2-5) and CCSD(T)/cc-pVTZ (n = 6) levels, reveal that the topology of the most stable isomer vary with n and the charge. Out of 674 harmonic frequencies, those with maximum intensity are generally in the 3000-3500 cm-1 range. Analysis of 169 vertical transition energies calculated with the EOM-CCSD approach, suggest three C6H2 species as potential carriers of the diffuse interstellar bands (DIB). Systematic comparison of properties between neutral and cationic species can assist in the structural description of complex matrices.
Collapse
Affiliation(s)
- Lenin J Díaz Soto
- Instituto de Química, Universidade Federal do Rio de Janeiro, Cidade Universitária, CT, Rio de Janeiro, Brazil
| | - Ricardo R Oliveira
- Instituto de Química, Universidade Federal do Rio de Janeiro, Cidade Universitária, CT, Rio de Janeiro, Brazil
| | - Leonardo Baptista
- Departamento de Química e Ambiental, Faculdade de Tecnologia, Universidade do Estado do Rio de Janeiro, Resende, Brazil
| | - Enio F da Silveira
- Departamento de Física, Pontifícia Universidade Católica do Rio de Janeiro, Rio de Janeiro, Brazil
| | | |
Collapse
|
6
|
Góbi S, Keresztes B, Schneiker A, Tarczay G. Hydrogen-atom-assisted processes on thioacetamide in para-H 2 matrix - formation of thiol tautomers. Phys Chem Chem Phys 2024; 26:21589-21597. [PMID: 39083030 DOI: 10.1039/d4cp02400a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/15/2024]
Abstract
TA has been isolated in low-temperature para-H2 matrices and it has been exposed to H atoms. In accordance with previous experimental results, TA exclusively exists in its more stable thione tautomeric form in the freshly deposited matrix. However, upon H atom generation, the bands belonging to the precursor start decreasing with the simultaneous appearance of new bands. By comparing the position of these new peaks with earlier results obtained in para-H2, they can be undoubtedly attributed to the presence of the higher-energy thiol tautomeric forms of TA. No other products could be observed, except for NH3. Quantum-chemical computations have been invoked to understand the mechanism behind the observed thione-thiol tautomerization assisted by H atoms. Accordingly, tautomerization starts with barrierless H-atom addition on the S atom of the thione resulting in the formation of the intermediate 1-amino-1-mercapto-ethyl radical (H3C-Ċ(-SH)-NH2), which has been detected tentatively for the first time. The second step is a barrierless H-atom induced H-abstraction from the -NH2 moiety of the H3C-Ċ(-SH)-NH2 radical. A comprehensive mechanism for tautomerization is proposed based on the experimental and theoretical results. Although earlier studies showed the possibility of TA thione-thiol tautomerization in cryogenic matrices achieved by energetic UV irradiation, the present study points out that this can also take place through a barrierless pathway by simply exposing the TA molecules to H atoms. As such, this is the first evidence for the occurrence of such a reaction in a matrix-isolated environment. The current results also help us better understand the mystery behind thione-thiol tautomerization in low-temperature environments.
Collapse
Affiliation(s)
- Sándor Góbi
- Laboratory of Molecular Spectroscopy, Institute of Chemistry, ELTE Eötvös Loránd University, PO Box 32, H-1518 Budapest, Hungary.
- MTA-ELTE Lendület Laboratory Astrochemistry Research Group, Institute of Chemistry, ELTE Eötvös Loránd University, PO Box 32, H-1518 Budapest, Hungary
| | - Barbara Keresztes
- Laboratory of Molecular Spectroscopy, Institute of Chemistry, ELTE Eötvös Loránd University, PO Box 32, H-1518 Budapest, Hungary.
- Hevesy György PhD School of Chemistry, Institute of Chemistry, ELTE Eötvös Loránd University, PO Box 32, H-1518 Budapest, Hungary
| | - Anita Schneiker
- Laboratory of Molecular Spectroscopy, Institute of Chemistry, ELTE Eötvös Loránd University, PO Box 32, H-1518 Budapest, Hungary.
- Hevesy György PhD School of Chemistry, Institute of Chemistry, ELTE Eötvös Loránd University, PO Box 32, H-1518 Budapest, Hungary
| | - György Tarczay
- Laboratory of Molecular Spectroscopy, Institute of Chemistry, ELTE Eötvös Loránd University, PO Box 32, H-1518 Budapest, Hungary.
- MTA-ELTE Lendület Laboratory Astrochemistry Research Group, Institute of Chemistry, ELTE Eötvös Loránd University, PO Box 32, H-1518 Budapest, Hungary
- Centre for Astrophysics and Space Science, ELTE Eötvös Loránd University, PO Box 32, H-1518 Budapest, Hungary
| |
Collapse
|
7
|
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
|
8
|
Lopes Jesus AJ, Nunes CM, Ferreira GA, Keyvan K, Fausto R. Photochemical Generation and Characterization of C-Aminophenyl-Nitrilimines: Insights on Their Bond-Shift Isomers by Matrix-Isolation IR Spectroscopy and Density Functional Theory Calculations. Molecules 2024; 29:3497. [PMID: 39124902 PMCID: PMC11314218 DOI: 10.3390/molecules29153497] [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/11/2024] [Revised: 07/20/2024] [Accepted: 07/22/2024] [Indexed: 08/12/2024] Open
Abstract
The intriguing ability of C-phenyl-nitrilimine to co-exist as allenic and propargylic bond-shift isomers motivated us to investigate how substituents in the phenyl ring influence this behavior. Building on our previous work on the meta- and para-OH substitution, here we extended this investigation to explore the effect of the NH2 substitution. For this purpose, C-(4-aminophenyl)- and C-(3-aminophenyl)-nitrilimines were photogenerated in an argon matrix at 15 K by narrowband UV-light irradiation (λ = 230 nm) of 5-(4-aminophenyl)- and 5-(3-aminophenyl)-tetrazole, respectively. The produced nitrilimines were further photoisomerized to carbodiimides via 1H-diazirines by irradiations at longer wavelengths (λ = 380 or 330 nm). Combining IR spectroscopic measurements and DFT calculations, it was found that the para-NH2-substituted nitrilimine exists as a single isomeric structure with a predominant allenic character. In contrast, the meta-NH2-substituted nitrilimine co-exists as two bond-shift isomers characterized by propargylic and allenic structures. To gain further understanding of the effects of phenyl substitution on the bond-shift isomerism of the nitrilimine fragment, we compared geometric and charge distribution data derived from theoretical calculations performed for C-phenyl-nitrilimine with those performed for the derivatives resulting from NH2 (electron-donating group) and NO2 (electron-withdrawing group) phenyl substitutions.
Collapse
Affiliation(s)
- A. J. Lopes Jesus
- University of Coimbra, CQC-IMS, Faculty of Pharmacy, 3004-295 Coimbra, Portugal
| | - Cláudio M. Nunes
- University of Coimbra, CQC-IMS, Department of Chemistry, 3004-535 Coimbra, Portugal; (C.M.N.); (G.A.F.); (K.K.); (R.F.)
| | - Gil A. Ferreira
- University of Coimbra, CQC-IMS, Department of Chemistry, 3004-535 Coimbra, Portugal; (C.M.N.); (G.A.F.); (K.K.); (R.F.)
| | - Kiarash Keyvan
- University of Coimbra, CQC-IMS, Department of Chemistry, 3004-535 Coimbra, Portugal; (C.M.N.); (G.A.F.); (K.K.); (R.F.)
| | - R. Fausto
- University of Coimbra, CQC-IMS, Department of Chemistry, 3004-535 Coimbra, Portugal; (C.M.N.); (G.A.F.); (K.K.); (R.F.)
- Istanbul Kultur University, Faculty Sciences and Letters, Department of Physics, Bakirkoy, 34158 Istanbul, Turkey
| |
Collapse
|
9
|
Xie F, Mendolicchio M, Omarouayache W, Murugachandran SI, Lei J, Gou Q, Sanz ME, Barone V, Schnell M. Structural and Electronic Evolution of Ethanolamine upon Microhydration: Insights from Hyperfine Resolved Rotational Spectroscopy. Angew Chem Int Ed Engl 2024:e202408622. [PMID: 38982982 DOI: 10.1002/anie.202408622] [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: 05/07/2024] [Revised: 06/21/2024] [Accepted: 06/21/2024] [Indexed: 07/11/2024]
Abstract
Ethanolamine hydrates containing from one to seven water molecules were identified via rotational spectroscopy with the aid of accurate quantum chemical methods considering anharmonic vibrational corrections. Ethanolamine undergoes significant conformational changes upon hydration to form energetically favorable hydrogen bond networks. The final structures strongly resemble the pure (H2O)3-9 complexes reported before when replacing two water molecules by ethanolamine. The 14N nuclear quadrupole coupling constants of all the ethanolamine hydrates have been determined and show a remarkable correlation with the strength of hydrogen bonds involving the amino group. After addition of the seventh water molecule, both hydrogen atoms of the amino group actively contribute to hydrogen bond formation, reinforcing the network and introducing approximately 21-27 % ionicity towards the formation of protonated amine. These findings highlight the critical role of microhydration in altering the electronic environment of ethanolamine, enhancing our understanding of amine hydration dynamics.
Collapse
Affiliation(s)
- Fan Xie
- Deutsches Elektronen-Synchrotron DESY, Notkestr. 85, 22607, Hamburg, Germany
| | | | | | | | - Juncheng Lei
- Department of Chemistry, School of Chemistry and Chemical Engineering, Chongqing Key Laboratory of Theoretical and Computational Chemistry, Chongqing University, 401331, Chongqing, China
| | - Qian Gou
- Department of Chemistry, School of Chemistry and Chemical Engineering, Chongqing Key Laboratory of Theoretical and Computational Chemistry, Chongqing University, 401331, Chongqing, China
| | - M Eugenia Sanz
- Department of Chemistry, King's College London, London, SE1 1DB, U.K
| | | | - Melanie Schnell
- Deutsches Elektronen-Synchrotron DESY, Notkestr. 85, 22607, Hamburg, Germany
- Institut für Physikalische Chemie, Christian-Albrechts-Universität zu Kiel, Max-Eyth-Str. 1, 24118, Kiel, Germany
| |
Collapse
|
10
|
Alberini A, Fornaro T, García-Florentino C, Biczysko M, Poblacion I, Aramendia J, Madariaga JM, Poggiali G, Vicente-Retortillo Á, Benison KC, Siljeström S, Biancalani S, Lorenz C, Cloutis EA, Applin DM, Gómez F, Steele A, Wiens RC, Hand KP, Brucato JR. Investigating the stability of aromatic carboxylic acids in hydrated magnesium sulfate under UV irradiation to assist detection of organics on Mars. Sci Rep 2024; 14:15945. [PMID: 38987581 PMCID: PMC11237158 DOI: 10.1038/s41598-024-66669-8] [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: 05/28/2024] [Accepted: 07/03/2024] [Indexed: 07/12/2024] Open
Abstract
The Scanning Habitable Environments with Raman and Luminescence for Organics and Chemicals (SHERLOC) instrument onboard the Mars 2020 Perseverance rover detected so far some of the most intense fluorescence signals in association with sulfates analyzing abraded patches of rocks at Jezero crater, Mars. To assess the plausibility of an organic origin of these signals, it is key to understand if organics can survive exposure to ambient Martian UV after exposure by the Perseverance abrasion tool and prior to analysis by SHERLOC. In this work, we investigated the stability of organo-sulfate assemblages under Martian-like UV irradiation and we observed that the spectroscopic features of phthalic and mellitic acid embedded into hydrated magnesium sulfate do not change for UV exposures corresponding to at least 48 Martian sols and, thus, should still be detectable in fluorescence when the SHERLOC analysis takes place, thanks to the photoprotective properties of magnesium sulfate. In addition, different photoproduct bands diagnostic of the parent carboxylic acid molecules could be observed. The photoprotective behavior of hydrated magnesium sulfate corroborates the hypothesis that sulfates might have played a key role in the preservation of organics on Mars, and that the fluorescence signals detected by SHERLOC in association with sulfates could potentially arise from organic compounds.
Collapse
Affiliation(s)
- Andrew Alberini
- INAF- Astrophysical Observatory of Arcetri, L.go E. Fermi 5, 50125, Firenze, Italy.
- Department of Physics and Astronomy, University of Florence, Via Giovanni Sansone 1, Sesto Fiorentino, 50019, Florence, Italy.
| | - Teresa Fornaro
- INAF- Astrophysical Observatory of Arcetri, L.go E. Fermi 5, 50125, Firenze, Italy.
| | - Cristina García-Florentino
- INAF- Astrophysical Observatory of Arcetri, L.go E. Fermi 5, 50125, Firenze, Italy
- Department of Analytical Chemistry, University of the Basque Country UPV/EHU, 48080, Bilbao, Spain
| | - Malgorzata Biczysko
- College of Science, Shanghai University, 99 Shangda Road, Shanghai, 200444, China
| | - Iratxe Poblacion
- Department of Analytical Chemistry, University of the Basque Country UPV/EHU, 48080, Bilbao, Spain
| | - Julene Aramendia
- Department of Analytical Chemistry, University of the Basque Country UPV/EHU, 48080, Bilbao, Spain
| | - Juan Manuel Madariaga
- Department of Analytical Chemistry, University of the Basque Country UPV/EHU, 48080, Bilbao, Spain
| | - Giovanni Poggiali
- INAF- Astrophysical Observatory of Arcetri, L.go E. Fermi 5, 50125, Firenze, Italy
- LESIA - Observatoire de Paris, CNRS, Université Paris Cité, Université PSL, Sorbonne Université, 5 Place Jules Janssen, 92190, Meudon, France
| | | | - Kathleen C Benison
- Department of Geology and Geography, West Virginia University, Morgantown, WV, USA
| | | | - Sole Biancalani
- INAF- Astrophysical Observatory of Arcetri, L.go E. Fermi 5, 50125, Firenze, Italy
- Department of Physics, University of Trento, Via Sommarive 14, 38123, Povo, Italy
- Italian Space Angency (ASI), Viale del Politecnico Snc, 00133, Rome, Italy
- Department of Earth Sciences, University of Florence, Via G. La Pira 4, 50121, Florence, Italy
| | - Christian Lorenz
- INAF- Astrophysical Observatory of Arcetri, L.go E. Fermi 5, 50125, Firenze, Italy
- Department of Biology, University of Naples Federico II, Via Cinthia, 80126, Naples, Italy
| | - Edward A Cloutis
- Centre for Terrestrial and Planetary Exploration, University of Winnipeg, Winnipeg, MB, R3B 2E9, Canada
| | - Dan M Applin
- Centre for Terrestrial and Planetary Exploration, University of Winnipeg, Winnipeg, MB, R3B 2E9, Canada
| | - Felipe Gómez
- Centro de Astrobiología (CAB), CSIC-INTA, Torrejón de Ardoz, Spain
| | | | - Roger C Wiens
- Earth, Atmospheric, and Planetary Sciences, Purdue University, West Lafayette, IN, USA
| | - Kevin P Hand
- Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA, USA
| | - John R Brucato
- INAF- Astrophysical Observatory of Arcetri, L.go E. Fermi 5, 50125, Firenze, Italy
| |
Collapse
|
11
|
Chawananon S, Goubet M, Pirali O, Georges R, Roucou A, Hadj Said I, Senent ML, Cuisset A, Asselin P. High resolution far-infrared synchrotron spectroscopy of 2-furfural conformers: Fundamental and hot bands. J Chem Phys 2024; 161:014308. [PMID: 38958165 DOI: 10.1063/5.0213834] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2024] [Accepted: 06/16/2024] [Indexed: 07/04/2024] Open
Abstract
In the continuity of a previous jet-cooled rovibrational study of trans and cis conformers of 2-furfural in the mid-infrared region (700-1750 cm-1) [Chawananon et al., Molecules 28 (10), 4165 (2023)], the present work investigates the far-infrared spectroscopy of 2-furfural using a long path absorption cell coupled to a high-resolution Fourier transform spectrometer and synchrotron radiation at the AILES beamline of the SOLEIL synchrotron. Guided by anharmonic calculations, vibrational energy levels and excited-state rotational constants are sufficiently predictive for a complete assignment of all fundamental and combination bands up to 700 cm-1, as well as the rovibrational analysis of 4 (1) low-frequency modes of trans-(cis-)2-furfural. A global rovibrational simulation, including far-infrared rovibrational lines and microwave and millimeter-wave rotational lines assigned in a previous study [Motiyenko et al., J. Mol. Spectrosc., 244, 9 (2007)] provides a reliable set of ground- and excited-state rotational parameters involving ring torsion, bending, and ring puckering modes of 2-furfural. In a second step, a rovibrational analysis of several hot band sequences, mainly involving the lowest frequency ring CHO torsion mode, is carried out. Reliable values of some anharmonic coefficients are obtained experimentally and could serve as a benchmark for validating advanced anharmonic calculations related to these large amplitude motions of flexible molecules.
Collapse
Affiliation(s)
- Sathapana Chawananon
- Sorbonne Université, CNRS, MONARIS, UMR 8233, 4 Place Jussieu, Paris F-75005, France
| | - Manuel Goubet
- Université de Lille, CNRS, UMR 8523 - PhLAM - Physique des Lasers Atomes et Molécules, F-59000 Lille, France
| | - Olivier Pirali
- Université de Paris-Saclay, CNRS, Institut des Sciences Moléculaires d'Orsay, F-91405 Orsay, France
| | - Robert Georges
- Université de Rennes, CNRS, IPR (Institut de Physique de Rennes)-UMR 6251, F-35000 Rennes, France
| | - Anthony Roucou
- Université du Littoral Côte d'Opale, UR4493, LPCA, Laboratoire de Physico-Chimie de l'Atmosphère, F-59140 Dunkerque, France
| | - Ikram Hadj Said
- USTHB, Laboratory of Thermodynamics and Molecular Modeling, Faculty of Chemistry, BP32, El Alia, 16111 Bab Ezzouar, Algiers, Algeria
- Departamento de Química y Física Teóricas, Instituto de Estructura de la Materia, IEM-CSIC, Unidad Asociada GIFMAN, CSIC-UHU, Serrano121, Madrid 28006, Spain
| | - María Luisa Senent
- Departamento de Química y Física Teóricas, Instituto de Estructura de la Materia, IEM-CSIC, Unidad Asociada GIFMAN, CSIC-UHU, Serrano121, Madrid 28006, Spain
| | - Arnaud Cuisset
- Université du Littoral Côte d'Opale, UR4493, LPCA, Laboratoire de Physico-Chimie de l'Atmosphère, F-59140 Dunkerque, France
| | - Pierre Asselin
- Sorbonne Université, CNRS, MONARIS, UMR 8233, 4 Place Jussieu, Paris F-75005, France
| |
Collapse
|
12
|
Uribe L, Lazzari F, Di Grande S, Crisci L, Mendolicchio M, Barone V. Accurate structures and rotational constants of bicyclic monoterpenes at DFT cost by means of the bond-corrected Pisa composite scheme (BPCS). J Chem Phys 2024; 161:014307. [PMID: 38958160 DOI: 10.1063/5.0216384] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2024] [Accepted: 06/17/2024] [Indexed: 07/04/2024] Open
Abstract
The structural, conformational, and spectroscopic properties in the gas phase of 20 bicyclic monoterpenes and monoterpenoids have been analyzed by a new accurate, reduced-cost computational strategy. In detail, the revDSD-PBEP86 double-hybrid functional in conjunction with the D3BJ empirical dispersion corrections and a suitable triple-zeta basis set provides accurate geometrical parameters, whence equilibrium rotational constants, which are further improved by proper account of core-valence correlation. Average deviations within 0.1% between computed and experimental rotational constants are reached when taking into account the vibrational corrections obtained by the B3LYP functional in conjunction with a double-zeta basis set in the framework of second-order vibrational perturbation theory. In addition to their intrinsic interest, the studied terpenes further extend the panel of systems for which the proposed strategy has provided accurate results at density functional theory cost. Therefore, a very accurate yet robust and user-friendly tool is now available for systematic investigations of the role of stereo-electronic effects on the properties of large systems of current technological and/or biological interest by experimentally oriented researchers.
Collapse
Affiliation(s)
- Lina Uribe
- 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
| | - 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
| | - Luigi Crisci
- 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
|
13
|
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
|
14
|
Missaoui D, Brahem S, Najar F, Yazidi O, Senent ML. The Far-Infrared Spectrum of Methoxymethanol (CH 3-O-CH 2OH): A Theoretical Study. ACS EARTH & SPACE CHEMISTRY 2024; 8:1236-1245. [PMID: 38919855 PMCID: PMC11194847 DOI: 10.1021/acsearthspacechem.4c00053] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 03/08/2024] [Revised: 05/16/2024] [Accepted: 05/17/2024] [Indexed: 06/27/2024]
Abstract
Methoxymethanol (CH3OCH2OH), an oxygenated volatile organic compound of low stability detected in the interstellar medium, represents an example of nonrigid organic molecules showing various interacting and inseparable large-amplitude motions. The species discloses a relevant coupling among torsional modes, strong enough to prevent complete assignments using effective Hamiltonians of reduced dimensionality. Theoretical models for rotational spectroscopy can improve if they treat three vibrational coordinates together. In this paper, the nonrigid properties and the far-infrared region are analyzed using highly correlated ab initio methods and a three-dimensional vibrational model. The molecule displays two gauche-gauche (CGcg and CGcg') and one trans-gauche (Tcg) conformers, whose relative energies are very small (CGcg/CGcg'/Tcg = 0.0:641.5:792.7 cm-1). The minima are separated by relatively low barriers (1200-1500 cm-1), and the corresponding methyl torsional barriers V 3 are estimated to be 595.7, 829.0, and 683.7 cm-1, respectively. The ground vibrational state rotational constants of the most stable geometry have been computed to be A 0 = 17233.99 MHz, B 0 = 5572.58 MHz, and C 0 = 4815.55 MHz, at ΔA 0 = -3.96 MHz, ΔB 0 = 4.76 MHz, and ΔC 0 = 2.51 MHz from previous experimental data. Low-energy levels and their tunneling splittings are determined variationally up to 700 cm-1. The A/E splitting of the ground vibrational state was computed to be 0.003 cm-1, as was expected given the methyl torsional barrier (∼600 cm-1). The fundamental levels (100), (010), and (001) are predicted at 132.133 and 132.086 cm-1 (methyl torsion), 186.507 and 186.467 cm-1 (O-CH3 torsion), and 371.925 and 371.950 cm-1 (OH torsion), respectively.
Collapse
Affiliation(s)
- Dorsaf Missaoui
- Laboratoire
de Spectroscopie Atomique Moléculaire et Applications, Faculté
des Sciences de Tunis, Université
de Tunis El Manar, Tunis 2092, Tunisia
- Departamento
de Química y Física Teóricas, Instituto de Estructura de la Materia, IEM-CSIC, Serrano 121, Madrid 28006, Spain
- Unidad
Asociada GIFMAN, CSIC-UHU, Huelva 21071, Spain
| | - Sinda Brahem
- Laboratoire
de Spectroscopie Atomique Moléculaire et Applications, Faculté
des Sciences de Tunis, Université
de Tunis El Manar, Tunis 2092, Tunisia
- Departamento
de Química y Física Teóricas, Instituto de Estructura de la Materia, IEM-CSIC, Serrano 121, Madrid 28006, Spain
- Unidad
Asociada GIFMAN, CSIC-UHU, Huelva 21071, Spain
| | - Faouzi Najar
- Laboratoire
de Spectroscopie Atomique Moléculaire et Applications, Faculté
des Sciences de Tunis, Université
de Tunis El Manar, Tunis 2092, Tunisia
| | - Ounaies Yazidi
- Laboratoire
de Spectroscopie Atomique Moléculaire et Applications, Faculté
des Sciences de Tunis, Université
de Tunis El Manar, Tunis 2092, Tunisia
| | - María Luisa Senent
- Departamento
de Química y Física Teóricas, Instituto de Estructura de la Materia, IEM-CSIC, Serrano 121, Madrid 28006, Spain
- Unidad
Asociada GIFMAN, CSIC-UHU, Huelva 21071, Spain
| |
Collapse
|
15
|
Esposito VJ, Fortenberry RC, Boersma C, Allamandola LJ. Assigning the CH stretch overtone spectrum of benzene and naphthalene with extension to anthracene and tetracene using 2- and 3-quanta anharmonic quantum chemical computations. J Chem Phys 2024; 160:211101. [PMID: 38828805 DOI: 10.1063/5.0208597] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2024] [Accepted: 05/05/2024] [Indexed: 06/05/2024] Open
Abstract
The CH stretch overtone region (5750-6300 cm-1) of benzene and naphthalene is assigned herein using anharmonic quantum chemical computations, and the trend of how this extends to larger polycyclic aromatic hydrocarbons (PAHs) is established. The assignment of all experimental bands to specific vibrational states is performed for the first time. Resonance polyads and the inclusion of 3-quanta vibrational states are both needed to compute accurate vibrational frequencies with the proper density-of-states to match the experimental band shape. Hundreds of 3-quanta states produce the observed band structure in naphthalene, anthracene, and tetracene, and this number is expected to increase drastically for larger PAHs. The width and shape of the main peak are consistent from naphthalene to anthracene, necessitating further exploration of this trend to confirm whether it is representative of all PAHs in the CH stretch overtone region. Understanding observations of PAH sources in the 1-3 μm region from the NIRSpec instrument aboard JWST requires new computational data, and this study provides a benchmark and foundation for their computation.
Collapse
Affiliation(s)
- Vincent J Esposito
- Astrophysics Branch, NASA Ames Research Center, Moffett Field, California 94035, USA
| | - Ryan C Fortenberry
- Department of Chemistry and Biochemistry, University of Mississippi, University, Mississippi 38677-1848, USA
| | - Christiaan Boersma
- Astrophysics Branch, NASA Ames Research Center, Moffett Field, California 94035, USA
| | - Louis J Allamandola
- Astrophysics Branch, NASA Ames Research Center, Moffett Field, California 94035, USA
| |
Collapse
|
16
|
Zou M, Hassan Y, Roy TK, McCoy AB, Lester MI. Infrared spectroscopy of the syn-methyl-substituted Criegee intermediate: A combined experimental and theoretical study. J Chem Phys 2024; 160:204309. [PMID: 38818894 DOI: 10.1063/5.0210122] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2024] [Accepted: 05/06/2024] [Indexed: 06/01/2024] Open
Abstract
An IR-vacuum ultraviolet (VUV) ion-dip spectroscopy method is utilized to examine the IR spectrum of acetaldehyde oxide (CH3CHOO) in the overtone CH stretch (2νCH) spectral region. IR activation creates a depletion of the ground state population that reduces the VUV photoionization signal on the parent mass channel. IR activation of the more stable and populated syn-CH3CHOO conformer results in rapid unimolecular decay to OH + vinoxy products and makes the most significant contribution to the observed spectrum. The resultant IR-VUV ion-dip spectrum of CH3CHOO is similar to that obtained previously for syn-CH3CHOO using IR action spectroscopy with UV laser-induced fluorescence detection of OH products. The prominent IR features at 5984 and 6081 cm-1 are also observed using UV + VUV photoionization of OH products. Complementary theoretical calculations utilizing a general implementation of second-order vibrational perturbation theory provide new insights on the vibrational transitions that give rise to the experimental spectrum in the overtone CH stretch region. The introduction of physically motivated small shifts of the harmonic frequencies yields remarkably improved agreement between experiment and theory in the overtone CH stretch region. The prominent features are assigned as highly mixed states with contributions from two quanta of CH stretch and/or a combination of CH stretch with an overtone in mode 4. The generality of this approach is demonstrated by applying it to three different levels of electronic structure theory/basis sets, all of which provide spectra that are virtually indistinguishable despite showing large deviations prior to introducing the shifts to the harmonic frequencies.
Collapse
Affiliation(s)
- Meijun Zou
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104-6323, USA
| | - Yarra Hassan
- Department of Chemistry, University of Washington, Seattle, Washington 98195-1700, USA
| | - Tarun Kumar Roy
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104-6323, USA
| | - Anne B McCoy
- Department of Chemistry, University of Washington, Seattle, Washington 98195-1700, USA
| | - Marsha I Lester
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104-6323, USA
| |
Collapse
|
17
|
Achevski B, Pejov L. Molecular Hydrogen Acts as a Hydrogen Bond Proton Acceptor: From Protonated Betaine Tagging to the Weakest Hydrogen Bond. J Phys Chem A 2024; 128:3968-3981. [PMID: 38720275 DOI: 10.1021/acs.jpca.4c01331] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/24/2024]
Abstract
In an attempt to gain further insights into the intermolecular interactions implied by Rizzo's group's cautionary tale related to molecular tagging in infrared multiple photon dissociation (IRMPD) spectroscopy with molecular messengers [Masson, A. . J. Chem. Phys. 2015, 143, 104313], in the present study, we provide an in-depth analysis of the noncovalent interaction between the molecular hydrogen and protonated betaine molecule in the gas phase. We aim to shed some new light on the fundamental issues concerning the wide diapason of hydrogen-bonding-type intermolecular interactions, with a wide variety of proton acceptors. We demonstrate that in the course of tagging the protonated betaine with molecular hydrogen from the OH group side, it is the σ bond of molecular hydrogen that plays the role of hydrogen-bonding proton acceptor. The tagging thus induces a small yet significant red shift of the protonated betaine O-H stretching mode. We investigate the performance of a wide range of density functional theory (DFT) functionals for the calculation of anharmonic vibrational frequency shifts of the studied system, which are essential for the correct interpretation of the experimental IRMPD data. For an accurate prediction of the OH stretching frequency shifts, specifically designed functionals such as Handy's group HCTH/407 should be applied. The empirical dispersion correction enhances the systematic overestimation of the anharmonic frequency shift, characteristic of the most widely used DFT functionals. Combining the full-wave function approach with the charge field perturbation and natural bond orbital (NBO) deletion analyses, we demonstrate that the frequency shift in the OH-tagged structure is governed by the σHH → σ*OH intermolecular charge transfer. This interaction stabilizes the OH-tagged dimer as well, in contrast to the dipole-quadrupole electrostatic interaction energy term. Topological analysis of the electron density reveals the presence of an intermolecular bond critical point with a positive value of the density Laplacian very close to the lower limit for hydrogen bonds. NCI analyses demonstrate that the OH···H2 interaction is weaker than the intramolecular CH···O one within the protonated betaine molecule, with the through of reduced density gradient appearing at less negative sign(λ2)·ρ values. Analyzing the O-H stretching vibrational potential with the second-generation absolutely localized molecular orbitals energy decomposition analysis (ALMO-EDA 2) revealed that in the case of betaineH(+) tagged from the OH group side, the permanent electrostatics (ΔEelec), polarization (ΔEpol), and charge-transfer (ΔEct) contributions to the total intermolecular interaction energy contribute favorably to the weak hydrogen bond formation and to the red shift of the fundamental O-H stretching frequency, the ΔEct contribution being the most significant in the last context. The Pauli repulsion term, on the other hand, favors an O-H stretching frequency blue shift as a consequence of the vibrational confinement effects.
Collapse
Affiliation(s)
- Blagoj Achevski
- Institute of Chemistry, Faculty of Natural Sciences and Mathematics, SS Cyril and Methodius University, PO Box 162, 1000 Skopje, Republic of North Macedonia
- Faculty of Pharmacy, SS Cyril and Methodius University, 1000 Skopje, Republic of North Macedonia
| | - Ljupcho Pejov
- Institute of Chemistry, Faculty of Natural Sciences and Mathematics, SS Cyril and Methodius University, PO Box 162, 1000 Skopje, Republic of North Macedonia
- Department of Chemistry, Bioscience and Environmental Engineering, Faculty of Science and Technology, University of Stavanger, Mailbox 8600, 4036 Stavanger, Norway
- The Polytechnic School, Ira A. Fulton Schools of Engineering, Arizona State University, 6049 S Backus Mall Mesa, Arizona 85212, United States
| |
Collapse
|
18
|
Kontodimas V, Yaman M, Greis K, Lettow M, Pagel K, Marianski M. Reinvestigation of the internal glycan rearrangement of Lewis a and blood group type H1 epitopes. Phys Chem Chem Phys 2024; 26:14160-14170. [PMID: 38712976 PMCID: PMC11147448 DOI: 10.1039/d3cp04491b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/08/2024]
Abstract
Protonated ions of fucose-containing oligosaccharides are prone to undergo internal glycan rearrangement which results in chimeric fragments that obfuscate mass-spectrometric analysis. Lack of accessible tools that would facilitate systematic analysis of glycans in the gas phase limits our understanding of this phenomenon. In this work, we use density functional theory modeling to interpret cryogenic IR spectra of Lewis a and blood group type H1 trisaccharides and to establish whether these trisaccharides undergo the rearrangement during gas-phase analysis. Structurally unconstrained search reveals that none of the parent ions constitute a thermodynamic global minimum. In contrast, predicted collision cross sections and anharmonic IR spectra provide a good match to available experimental data which allowed us to conclude that fucose migration does not occur in these antigens. By comparing the predicted structures with those obtained for Lewis x and blood group type H2 epitopes, we demonstrate that the availability of the mobile proton and a large difference in the relative stability of the parent ions and rearrangement products constitute the prerequisites for the rearrangement reaction.
Collapse
Affiliation(s)
- Vasilis Kontodimas
- Department of Chemistry, Hunter College, The City University of New York, New York, NY 10065, USA.
| | - Murat Yaman
- Department of Chemistry, Hunter College, The City University of New York, New York, NY 10065, USA.
- Ph.D. Programs in Chemistry and Biochemistry, The Graduate Center, The City University of New York, New York, NY 10016, USA
| | - Kim Greis
- Fritz-Haber-Intitut der Max Planck Gesellschaft, 14195 Berlin, Germany
- Institute of Chemistry and Biochemistry, Freie Universität Berlin, 14195 Berlin, Germany
- Laboratory of Organic Chemistry, Department of Chemistry and Applied Biosciences, ETH Zürich, 8093 Zürich, Switzerland
| | - Maike Lettow
- Fritz-Haber-Intitut der Max Planck Gesellschaft, 14195 Berlin, Germany
- Institute of Chemistry and Biochemistry, Freie Universität Berlin, 14195 Berlin, Germany
| | - Kevin Pagel
- Fritz-Haber-Intitut der Max Planck Gesellschaft, 14195 Berlin, Germany
- Institute of Chemistry and Biochemistry, Freie Universität Berlin, 14195 Berlin, Germany
| | - Mateusz Marianski
- Department of Chemistry, Hunter College, The City University of New York, New York, NY 10065, USA.
- Ph.D. Programs in Chemistry and Biochemistry, The Graduate Center, The City University of New York, New York, NY 10016, USA
| |
Collapse
|
19
|
McIntosh O, García-Florentino C, Fornaro T, Marabello D, Alberini A, Siljeström S, Biczysko M, Szopa C, Brucato J. Undecanoic Acid and L-Phenylalanine in Vermiculite: Detection, Characterization, and UV Degradation Studies for Biosignature Identification on Mars. ASTROBIOLOGY 2024; 24:518-537. [PMID: 38669050 DOI: 10.1089/ast.2023.0088] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/22/2024]
Abstract
Solar radiation that arrives on the surface of Mars interacts with organic molecules present in the soil. The radiation can degrade or transform the organic matter and make the search for biosignatures on the planet's surface difficult. Therefore, samples to be analyzed by instruments on board Mars probes for molecular content should be selectively chosen to have the highest organic preservation content. To support the identification of organic molecules on Mars, the behavior under UV irradiation of two organic compounds, undecanoic acid and L-phenylalanine, in the presence of vermiculite and two chloride salts, NaCl and MgCl, was studied. The degradation of the molecule's bands was monitored through IR spectroscopy. Our results show that, while vermiculite acts as a photoprotective mineral with L-phenylalanine, it catalyzes the photodegradation of undecanoic acid molecules. On the other hand, both chloride salts studied decreased the degradation of both organic species acting as photoprotectors. While these results do not allow us to conclude on the preservation capabilities of vermiculite, they show that places where chloride salts are present could be good candidates for in situ analytic experiments on Mars due to their organic preservation capacity under UV radiation.
Collapse
Affiliation(s)
- Ophélie McIntosh
- LATMOS/IPSL, UVSQ Université Paris-Saclay, Sorbonne Université, CNRS, Guyancourt, France
- INAF - Astrophysical Observatory of Arcetri, Firenze, Italy
| | - Cristina García-Florentino
- INAF - Astrophysical Observatory of Arcetri, Firenze, Italy
- Department of Analytical Chemistry, University of the Basque Country UPV/EHU, Bilbao, Spain
| | - Teresa Fornaro
- INAF - Astrophysical Observatory of Arcetri, Firenze, Italy
| | - Domenica Marabello
- Department of Chemistry, University of Torino, Torino, Italy
- Interdepartmental Center for Crystallography, University of Torino, Torino, Italy
| | | | - Sandra Siljeström
- Department of Methodology, Textiles and Medical Technology, RISE Research Institutes of Sweden, Stockholm, Sweden
| | - Malgorzata Biczysko
- International Centre for Quantum and Molecular Structures, Physics Department, College of Science, Shanghai University, Shanghai, China
| | - Cyril Szopa
- LATMOS/IPSL, UVSQ Université Paris-Saclay, Sorbonne Université, CNRS, Guyancourt, France
| | - John Brucato
- INAF - Astrophysical Observatory of Arcetri, Firenze, Italy
| |
Collapse
|
20
|
Brahem S, Missaoui D, Yazidi O, Najar F, Senent ML. Theoretical structural and spectroscopic characterization of peroxyacetic acid (CH 3-CO-OOH): study of the far infrared region. Phys Chem Chem Phys 2024; 26:12600-12609. [PMID: 38597218 DOI: 10.1039/d3cp05783f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/11/2024]
Abstract
Peroxyacetic acid, a non-rigid oxygenated organic molecule which acts in the atmosphere as a reservoir of HOX and ROX radicals, is studied using highly correlated ab initio methods with the aim of its spectroscopic characterization in the gas phase. The study focuses on the far infrared region providing reliable rovibrational parameters such as energy levels and splittings. The molecule presents three conformers that inter-convert by internal rotation, drawing a potential energy surface of 12 minima. One of them shows prominent stability due to the formation of one weak intramolecular bond between the hydrogen atom of the hydroperoxy group and the oxygen atom of the carbonyl group. For the three minimum energy structures, rotational constants and centrifugal distortion constants are provided. It may be expected that the most stable conformer is the only one contributing to the spectral features in further measurements at low temperature. In this structure, the methyl torsional barrier has been found to be very low, V3 = 88.6 cm-1 producing a splitting of 2.262 cm-1 for the ground vibrational state. The study confirms that the ν20 torsional mode interacts strongly with the other two torsional modes ν19 and ν21, but slightly with the remaining vibrations. Then, a variational procedure in three dimensions allows the exploration of the low-frequency modes. The methyl torsional fundamental ν21 was found to be 49.1 cm-1 (Ai) and 33.4 cm-1 (E). The fundamentals of ν20 (C-O bond torsion) and ν19 (OH torsion) have been computed to be 216.7 cm-1 (A2) and 218.5 cm-1 (E) and 393.6 cm-1 (A2) and 394.1 cm-1. Since non-rigidity can have effects on the reactivity due to the conformer interconversion, and transitions involving low-lying levels can be observed with many spectroscopic techniques, this work can help kinetic studies and assignments of further spectroscopic studies needed for the detection in the gas phase of trace molecules.
Collapse
Affiliation(s)
- Sinda Brahem
- Laboratoire de Spectroscopie Atomique Moléculaire et Applications, Faculté des Sciences de Tunis, Université de Tunis El Manar, 2092, Tunisia
- Departamento de Química y Física Teóricas, Instituto de Estructura de la Materia, IEM-CSIC, Serrano 121, Madrid 28006, Spain
- Unidad Asociada GIFMAN, CSIC-UHU, Spain.
| | - Dorsaf Missaoui
- Laboratoire de Spectroscopie Atomique Moléculaire et Applications, Faculté des Sciences de Tunis, Université de Tunis El Manar, 2092, Tunisia
- Departamento de Química y Física Teóricas, Instituto de Estructura de la Materia, IEM-CSIC, Serrano 121, Madrid 28006, Spain
- Unidad Asociada GIFMAN, CSIC-UHU, Spain.
| | - Ounaies Yazidi
- Laboratoire de Spectroscopie Atomique Moléculaire et Applications, Faculté des Sciences de Tunis, Université de Tunis El Manar, 2092, Tunisia
| | - Faouzi Najar
- Laboratoire de Spectroscopie Atomique Moléculaire et Applications, Faculté des Sciences de Tunis, Université de Tunis El Manar, 2092, Tunisia
| | - María Luisa Senent
- Departamento de Química y Física Teóricas, Instituto de Estructura de la Materia, IEM-CSIC, Serrano 121, Madrid 28006, Spain
- Unidad Asociada GIFMAN, CSIC-UHU, Spain.
| |
Collapse
|
21
|
Ma J, Insausti A, Al-Jabiri MH, Carlson CD, Jäger W, Xu Y. Unlocking a new hydrogen-bonding marker: C-O bond shortening in vicinal diols revealed by rotational spectroscopy. J Chem Phys 2024; 160:154304. [PMID: 38624119 DOI: 10.1063/5.0203123] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2024] [Accepted: 03/28/2024] [Indexed: 04/17/2024] Open
Abstract
The conformational space of cis-1,2-cyclohexanediol, a model molecule for cyclic vicinal diols, was investigated using rotational spectroscopy and density functional theory calculations. Four low energy conformers within an energy window of 5 kJ mol-1 were identified computationally. A rotational spectrum of jet-cooled cis-1,2-cyclohexanediol was recorded with a chirped pulse Fourier transform microwave spectrometer. Two sets of rotational transitions were observed and could be assigned to conformers of cis-1,2-cyclohexanediol. The non-observation of other low energy conformers was explained by conformational conversion barrier height calculations and results from experimental spectra recorded with different carrier gases. Eight isotopologues, including those with 13C and 18O, of the lowest energy conformer were observed, allowing the determination of the semi-experimental equilibrium structure, reSE. Interestingly, the structural analysis revealed that the C-O bond length of the intramolecular hydrogen-bond donor is shorter than that of the acceptor. This appears to be a general characteristic of vicinal diols and can be used as a novel hydrogen-bond marker in such compounds.
Collapse
Affiliation(s)
- Jiarui Ma
- Department of Chemistry, University of Alberta, Edmonton, Alberta T6G 2G2, Canada
| | - Aran Insausti
- Department of Chemistry, University of Alberta, Edmonton, Alberta T6G 2G2, Canada
- Departamento de Química Física, Universidad del País Vasco (UPV/EHU) Barrio Sarriena s/n, Leioa, Spain
- Instituto Biofisika (UPV/EHU, CSIC) Barrio Sarriena s/n, Leioa, Spain
| | - Mohamad H Al-Jabiri
- Department of Chemistry, University of Alberta, Edmonton, Alberta T6G 2G2, Canada
| | - Colton D Carlson
- Department of Chemistry, University of Alberta, Edmonton, Alberta T6G 2G2, Canada
| | - Wolfgang Jäger
- Department of Chemistry, University of Alberta, Edmonton, Alberta T6G 2G2, Canada
| | - Yunjie Xu
- Department of Chemistry, University of Alberta, Edmonton, Alberta T6G 2G2, Canada
| |
Collapse
|
22
|
Wensink FJ, Smink CE, Stevenson BC, Steele RP, Bakker JM, Armentrout PB. IR spectroscopic characterization of [M,C,2H] + (M = Ru and Rh) products formed by reacting 4d transition metal cations with oxirane: Spectroscopic evidence for multireference character in RhCH 2. Phys Chem Chem Phys 2024; 26:11445-11458. [PMID: 38572552 PMCID: PMC11022548 DOI: 10.1039/d4cp00012a] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/02/2024] [Accepted: 03/13/2024] [Indexed: 04/05/2024]
Abstract
A combination of infrared multiple-photon dissociation (IRMPD) action spectroscopy and quantum chemical calculations was employed to investigate the [M,C,2H]+ (M = Ru and Rh) species. These ions were formed by reacting laser ablated M+ ions with oxirane (ethylene oxide, c-C2H4O) in a room-temperature ion trap. IRMPD spectra for the Ru species exhibit one major band and two side bands, whereas spectra for the Rh species contain more distinct bands. Comparison with density functional theory (DFT), coupled-cluster (CCSD), and equation-of-motion spin-flip CCSD (EOM-SF-CCSD) calculations allows assignment of the [M,C,2H]+ structures. For the spectrum of [Ru,C,2H]+, a combination of HRuCH+ and RuCH2+ structures reproduces the observed spectrum at all levels of theory. The well-resolved spectrum of [Rh,C,2H]+ could not be assigned unambiguously to any calculated structure using DFT approaches. The EOM-SF-CCSD calculations showed that the ground-state surface has multireference electronic character, and symmetric carbenes in both the 1A1 and 3A2 states are needed to reproduce the observed spectrum.
Collapse
Affiliation(s)
- Frank J Wensink
- Radboud University, Institute for Molecules and Materials, HFML-FELIX, Toernooiveld 7, 6525 ED Nijmegen, The Netherlands.
| | - Corry E Smink
- Radboud University, Institute for Molecules and Materials, HFML-FELIX, Toernooiveld 7, 6525 ED Nijmegen, The Netherlands.
| | - Brandon C Stevenson
- Department of Chemistry, University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112, USA.
| | - Ryan P Steele
- Department of Chemistry, University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112, USA.
| | - Joost M Bakker
- Radboud University, Institute for Molecules and Materials, HFML-FELIX, Toernooiveld 7, 6525 ED Nijmegen, The Netherlands.
| | - P B Armentrout
- Department of Chemistry, University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112, USA.
| |
Collapse
|
23
|
Fusè M, Mazzeo G, Longhi G, Abbate S, Yang Q, Bloino J. Scaling-up VPT2: A feasible route to include anharmonic correction on large molecules. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2024; 311:123969. [PMID: 38330757 DOI: 10.1016/j.saa.2024.123969] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/02/2023] [Revised: 01/18/2024] [Accepted: 01/24/2024] [Indexed: 02/10/2024]
Abstract
Vibrational analysis plays a crucial role in the investigation of molecular systems. Though methodologies like second-order vibrational perturbation theory (VPT2) have paved the way to more accurate simulations, the computational cost remains a difficult barrier to overcome when the molecular size increases. Building upon recent advances in the identification of resonances, we propose an approach making anharmonic simulations possible for large-size systems, typically unreachable by standard means. This relies on the fact that, often, only portions of the whole spectra are of actual interest. Therefore, the anharmonic corrections can be included selectively on subsets of normal modes directly related to the regions of interest. Starting from the VPT2 equations, we evaluate rigorously and systematically the impact of the truncated anharmonic treatment onto simulations. The limit and feasibility of the reduced-dimensionality approach are detailed, starting on a smaller model system. The methodology is then challenged on the IR absorption and vibrational circular dichroism spectra of an organometallic complex in three different spectral ranges.
Collapse
Affiliation(s)
- Marco Fusè
- Dipartimento di Medicina Molecolare e Traslazionale, Università di Brescia, Viale Europa 11, 25123, Brescia, Italy.
| | - Giuseppe Mazzeo
- Dipartimento di Medicina Molecolare e Traslazionale, Università di Brescia, Viale Europa 11, 25123, Brescia, Italy
| | - Giovanna Longhi
- Dipartimento di Medicina Molecolare e Traslazionale, Università di Brescia, Viale Europa 11, 25123, Brescia, Italy; Istituto Nazionale di Ottica (INO), CNR, Research Unit of Brescia, c/o CSMT, VIA Branze 45, 25123, Brescia, Italy
| | - Sergio Abbate
- Dipartimento di Medicina Molecolare e Traslazionale, Università di Brescia, Viale Europa 11, 25123, Brescia, Italy; Istituto Nazionale di Ottica (INO), CNR, Research Unit of Brescia, c/o CSMT, VIA Branze 45, 25123, Brescia, Italy
| | - Qin Yang
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo náměstí 542/2, 160 00, Prague, Czech Republic
| | - Julien Bloino
- Scuola Normale Superiore, Piazza dei Cavalieri, 56125, Pisa, Italy.
| |
Collapse
|
24
|
Mendolicchio M, Barone V. Unbiased Comparison between Theoretical and Experimental Molecular Structures and Properties: Toward an Accurate Reduced-Cost Evaluation of Vibrational Contributions. J Chem Theory Comput 2024; 20:2842-2857. [PMID: 38556752 DOI: 10.1021/acs.jctc.4c00023] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/02/2024]
Abstract
The tremendous development of hardware and software is constantly increasing the role of quantum chemical (QC) computations in the assignment and interpretation of experimental results. However, an unbiased comparison between theory and experiment requires the proper account of vibrational averaging effects. In particular, high-resolution spectra in the gas phase are now available for molecules containing up to about 50 atoms, which are too large for a brute-force approach with the available QC methods of sufficient accuracy. In the present paper, we introduce hybrid approaches, which allow the accurate evaluation of vibrational averaging effects for molecules of this size beyond the harmonic approximation, with special attention being devoted to rotational constants. After the validation of new tools for relatively small molecules, the β-estradiol hormone and a prototypical molecular motor have been considered to witness the feasibility of accurate computations for large molecules.
Collapse
|
25
|
Uribe L, Di Grande S, Crisci L, Lazzari F, Mendolicchio M, Barone V. Accurate Structures and Rotational Constants of Steroid Hormones at DFT Cost: Androsterone, Testosterone, Estrone, β-Estradiol, and Estriol. J Phys Chem A 2024; 128:2629-2642. [PMID: 38530336 DOI: 10.1021/acs.jpca.4c00573] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/27/2024]
Abstract
A comprehensive analysis of the structural, conformational, and spectroscopic properties in the gas phase has been performed for five prototypical steroid hormones, namely, androsterone, testosterone, estrone, β-estradiol, and estriol. The revDSD-PBEP86 double-hybrid functional in conjunction with the D3BJ empirical dispersion and a suitable triple-ζ basis set provides accurate conformational energies and equilibrium molecular structures, with the latter being further improved by proper account of core-valence correlation. Average deviations within 0.1% between computed and experimental ground state rotational constants are reached when adding to those equilibrium values vibrational corrections obtained at the cost of standard harmonic frequencies thanks to the use of a new computational tool. Together with the intrinsic interest of the studied hormones, the accuracy of the results obtained at DFT cost for molecules containing about 50 atoms paves the way toward the accurate investigations of other flexible bricks of life.
Collapse
Affiliation(s)
- Lina Uribe
- 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
| | - Luigi Crisci
- Scuola Normale Superiore di Pisa, Piazza dei Cavalieri 7, 56126 Pisa, 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
|
26
|
Khanifaev J, Schrader T, Perlt E. The effect of machine learning predicted anharmonic frequencies on thermodynamic properties of fluid hydrogen fluoride. J Chem Phys 2024; 160:124302. [PMID: 38516969 DOI: 10.1063/5.0195386] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/02/2024] [Accepted: 03/02/2024] [Indexed: 03/23/2024] Open
Abstract
Anharmonic effects play a crucial role in determining thermochemical properties of liquids and gases. For such extended phases, the inclusion of anharmonicity in reliable electronic structure methods is computationally extremely demanding, and hence, anharmonic effects are often lacking in thermochemical calculations. In this study, we apply the quantum cluster equilibrium method to transfer density functional theory calculations at the cluster level to the macroscopic, liquid, and gaseous phase of hydrogen fluoride. This allows us to include anharmonicity, either via vibrational self-consistent field calculations for smaller clusters or using a regression model for larger clusters. We obtain the structural composition of the fluid phases in terms of the population of different clusters as well as isobaric heat capacities as an example for thermodynamic properties. We study the role of anharmonicities for these analyses and observe that, in particular, the dominating structural motifs are rather sensitive to the anharmonicity in vibrational frequencies. The regression model proves to be a promising way to get access to anharmonic features, and the extension to more sophisticated machine-learning models is promising.
Collapse
Affiliation(s)
- Jamoliddin Khanifaev
- Otto Schott Institute of Materials Research, Friedrich Schiller University Jena, 07743 Jena, Germany
| | - Tim Schrader
- Otto Schott Institute of Materials Research, Friedrich Schiller University Jena, 07743 Jena, Germany
| | - Eva Perlt
- Otto Schott Institute of Materials Research, Friedrich Schiller University Jena, 07743 Jena, Germany
| |
Collapse
|
27
|
Crișan G, Stan Ș, Chiș V. Exploring Geometrical, Electronic and Spectroscopic Properties of 2-Nitroimidazole-Based Radiopharmaceuticals via Computational Chemistry Methods. Molecules 2024; 29:1505. [PMID: 38611785 PMCID: PMC11013577 DOI: 10.3390/molecules29071505] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2024] [Revised: 03/22/2024] [Accepted: 03/25/2024] [Indexed: 04/14/2024] Open
Abstract
Tumor hypoxia plays an important role in the clinical management and treatment planning of various cancers. The use of 2-nitroimidazole-based radiopharmaceuticals has been the most successful for positron emission tomography (PET) and single photon emission computed tomography (SPECT) imaging probes, offering noninvasive means to assess tumor hypoxia. In this study we performed detailed computational investigations of the most used compounds for PET imaging, focusing on those derived from 2-nitroimidazole: fluoromisonidazole (FMISO), fluoroazomycin arabinoside (FAZA), fluoroetanidazole (FETA), fluoroerythronitroimidazole (FETNIM) and 2-(2-nitroimidazol-1-yl)-N-(2,2,3,3,3-pentafluoropropyl)acetamide (EF5). Conformational analysis, structural parameters, vibrational IR and Raman properties (within both harmonic and anharmonic approximations), as well as the NMR shielding tensors and spin-spin coupling constants were obtained by density functional theory (DFT) calculations and then correlated with experimental findings, where available. Furthermore, time-dependent DFT computations reveal insight into the excited states of the compounds. Our results predict a significant change in the conformational landscape of most of the investigated compounds when transitioning from the gas phase to aqueous solution. According to computational data, the 2-nitroimidazole moiety determines to a large extent the spectroscopic properties of its derivatives. Due to the limited structural information available in the current literature for the investigated compounds, the findings presented herein deepen the current understanding of the electronic structures of these five radiopharmaceuticals.
Collapse
Affiliation(s)
- George Crișan
- Faculty of Physics, Babeș-Bolyai University, Str. M. Kogălniceanu 1, RO-400084 Cluj-Napoca, Romania; (G.C.); (Ș.S.)
- Department of Nuclear Medicine, County Clinical Hospital, Clinicilor 3-5, RO-400006 Cluj-Napoca, Romania
| | - Ștefan Stan
- Faculty of Physics, Babeș-Bolyai University, Str. M. Kogălniceanu 1, RO-400084 Cluj-Napoca, Romania; (G.C.); (Ș.S.)
| | - Vasile Chiș
- Faculty of Physics, Babeș-Bolyai University, Str. M. Kogălniceanu 1, RO-400084 Cluj-Napoca, Romania; (G.C.); (Ș.S.)
- Institute for Research, Development and Innovation in Applied Natural Sciences, Babeș-Bolyai University, Str. Fântânele 30, RO-400327 Cluj-Napoca, Romania
| |
Collapse
|
28
|
Esposito VJ, Ferrari P, Buma WJ, Fortenberry RC, Boersma C, Candian A, Tielens AGGM. The infrared absorption spectrum of phenylacetylene and its deuterated isotopologue in the mid- to far-IR. J Chem Phys 2024; 160:114312. [PMID: 38501470 DOI: 10.1063/5.0191404] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2023] [Accepted: 02/18/2024] [Indexed: 03/20/2024] Open
Abstract
Anharmonicity strongly influences the absorption and emission spectra of polycyclic aromatic hydrocarbon (PAH) molecules. Here, IR-UV ion-dip spectroscopy experiments together with detailed anharmonic computations reveal the presence of fundamental, overtone, as well as 2- and 3-quanta combination band transitions in the far- and mid-infrared absorption spectra of phenylacetylene and its singly deuterated isotopologue. Strong absorption features in the 400-900 cm-1 range originate from CH(D) in-plane and out-of-plane wags and bends, as well as bending motions including the C≡C and CH bonds of the acetylene substituent and the aromatic ring. For phenylacetylene, every absorption feature is assigned either directly or indirectly to a single or multiple vibrational mode(s). The measured spectrum is dense, broad, and structureless in many regions but well characterized by computations. Upon deuteration, large isotopic shifts are observed. At frequencies above 1500 cm-1 for d1-phenylacetylene, a one-to-one match is seen when comparing computations and experiments with all features assigned to combination bands and overtones. The C≡C stretch observed in phenylacetylene is not observed in d1-phenylacetylene due to a computed 40-fold drop in intensity. Overall, a careful treatment of anharmonicity that includes 2- and 3-quanta modes is found to be crucial to understand the rich details of the infrared spectrum of phenylacetylene. Based on these results, it can be expected that such an all-inclusive anharmonic treatment will also be key for unraveling the infrared spectra of PAHs in general.
Collapse
Affiliation(s)
- Vincent J Esposito
- NASA Ames Research Center, MS 245-6, Moffett Field, California 94035, USA
| | - Piero Ferrari
- Radboud University, Institute for Molecules and Materials, HFML-FELIX, 6525 ED Nijmegen, The Netherlands
| | - Wybren Jan Buma
- Radboud University, Institute for Molecules and Materials, HFML-FELIX, 6525 ED Nijmegen, The Netherlands
- Van't Hoff Institute for Molecular Sciences, University of Amsterdam, 1098 XH Amsterdam, The Netherlands
| | - Ryan C Fortenberry
- Department of Chemistry and Biochemistry, University of Mississippi, University, Mississippi 38677-1848, USA
| | - Christiaan Boersma
- NASA Ames Research Center, MS 245-6, Moffett Field, California 94035, USA
| | - Alessandra Candian
- Anton Pannekoek Institute for Astronomy, University of Amsterdam, 1098 XH Amsterdam, The Netherlands
| | - Alexander G G M Tielens
- Leiden Observatory, Leiden University, 2333 CA Leiden, The Netherlands
- Astronomy Department, University of Maryland, College Park, Maryland 20742-2421, USA
| |
Collapse
|
29
|
Hartwig B, Schnell M, Suhm MA, Obenchain DA. Weak hydrogen bonding to halogens and chirality communication in propanols: Raman and microwave spectroscopy benchmark theory. Phys Chem Chem Phys 2024; 26:9432-9452. [PMID: 38446207 DOI: 10.1039/d3cp04911f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/07/2024]
Abstract
Constitutional and conformational isomers of bromopropanol are vibrationally and rotationally characterised with parallels drawn to the structural chlorine analogues. A previous microwave spectroscopic study of the chloropropanols is re-examined and all systems are explored by Raman jet spectroscopy. For bromine, the entire nuclear quadrupole coupling tensors are accurately determined and compared to their chlorine counterparts. Tensor asymmetry parameters are determined and linked with the hydrogen bond strength as indicated by the downshift of the OH-stretching frequency. The spectroscopic constants derived from the observed transitions are used as benchmarks for a large variety of electronic structure methods followed by harmonic and anharmonic rovibrational treatments. The CCSD(T) electronic structure calculations provide the best performance, in particular once anharmonic and relativistic corrections are applied or implied. Standard DFT approaches vary substantially with respect to their systematic error cancellation across the investigated species, and cost-effective compromises for the different observables are proposed.
Collapse
Affiliation(s)
- Beppo Hartwig
- Institut für Physikalische Chemie, Tammannstrasse 6, 37077 Göttingen, Germany.
| | - Melanie Schnell
- Deutsches Elektronen-Synchrotron DESY, Notkestr. 85, 22607 Hamburg, Germany
- Institute of Physical Chemistry, Christian-Albrechts-Universität zu Kiel, Max-Eyth-Straße 1, 24118 Kiel, Germany
| | - Martin A Suhm
- Institut für Physikalische Chemie, Tammannstrasse 6, 37077 Göttingen, Germany.
| | - Daniel A Obenchain
- Institut für Physikalische Chemie, Tammannstrasse 6, 37077 Göttingen, Germany.
| |
Collapse
|
30
|
Käser S, Meuwly M. Numerical Accuracy Matters: Applications of Machine Learned Potential Energy Surfaces. J Phys Chem Lett 2024:3419-3424. [PMID: 38506827 DOI: 10.1021/acs.jpclett.3c03405] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/21/2024]
Abstract
The role of numerical accuracy in training and evaluating neural network-based potential energy surfaces is examined for different experimental observables. For observables that require third- and fourth-order derivatives of the potential energy with respect to Cartesian coordinates single-precision arithmetics as is typically used in ML-based approaches is insufficient and leads to roughness of the underlying PES as is explicitly demonstrated. Increasing the numerical accuracy to double-precision gives a smooth PES with higher-order derivatives that are numerically stable and yield meaningful anharmonic frequencies and tunneling splitting as is demonstrated for H2CO and malonaldehyde. For molecular dynamics simulations, which only require first-order derivatives, single-precision arithmetics appears to be sufficient, though.
Collapse
Affiliation(s)
- Silvan Käser
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| |
Collapse
|
31
|
Fusè M, Mazzeo G, Abbate S, Ruzziconi R, Bloino J, Longhi G. Mid-IR and CH stretching vibrational circular dichroism spectroscopy to distinguish various sources of chirality: The case of quinophaneoxazoline based ruthenium(II) complexes. Chirality 2024; 36:e23649. [PMID: 38409881 DOI: 10.1002/chir.23649] [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: 12/13/2023] [Revised: 01/13/2024] [Accepted: 01/15/2024] [Indexed: 02/28/2024]
Abstract
Five diastereomers of ruthenium(II) complexes based on quinolinophaneoxazoline ligands were investigated by vibrational circular dichroism (VCD) in the mid-IR and CH stretching regions. Diastereomers differ in three sources of chirality: the planar chirality of the quinolinophane moiety, the central chirality of an asymmetric carbon atom of the oxazoline ring, and the chirality of the ruthenium atom. VCD, allied to DFT calculations, has been found to be effective in disentangling the various forms of chirality. In particular, a VCD band is identified in the CH stretching region directly connected to the chirality of the metal. The analysis of the calculated VCD spectra is carried out by partitioning the complexes into fragments. The anharmonic analysis is also performed with a recently proposed reduced-dimensionality approach: such treatment is particularly important when examining spectroscopic regions highly perturbed by resonances, like the CH stretching region.
Collapse
Affiliation(s)
- Marco Fusè
- Dipartimento di Medicina Molecolare e Traslazionale (DMMT), Università di Brescia, Brescia, Italy
| | - Giuseppe Mazzeo
- Dipartimento di Medicina Molecolare e Traslazionale (DMMT), Università di Brescia, Brescia, Italy
| | - Sergio Abbate
- Dipartimento di Medicina Molecolare e Traslazionale (DMMT), Università di Brescia, Brescia, Italy
- Istituto Nazionale di Ottica (INO), CNR, Research Unit Brescia, Brescia, Italy
| | - Renzo Ruzziconi
- Dipartimento di Chimica Biologia e Biotecnologie, Università di Perugia, Perugia, Italy
| | - Julien Bloino
- Scuola Normale Superiore, Piazza dei Cavalieri, Pisa, Italy
| | - Giovanna Longhi
- Dipartimento di Medicina Molecolare e Traslazionale (DMMT), Università di Brescia, Brescia, Italy
- Istituto Nazionale di Ottica (INO), CNR, Research Unit Brescia, Brescia, Italy
| |
Collapse
|
32
|
Lazzari F, Mendolicchio M, Barone V. Accurate Geometries of Large Molecules by Integration of the Pisa Composite Scheme and the Templating Synthon Approach. J Phys Chem A 2024; 128:1385-1395. [PMID: 38347709 DOI: 10.1021/acs.jpca.3c08382] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/23/2024]
Abstract
An effective yet reliable computational workflow is proposed, which permits the computation of accurate geometrical structures for large flexible molecules at an affordable cost thanks to the integration of machine learning tools and DFT models together with reduced scaling computations of vibrational averaging effects. After validation of the different components of the overall strategy, a panel of molecules of biological interest have been analyzed. The results confirm that very accurate geometrical parameters can be obtained at reasonable cost for molecules including up to about 50 atoms, which are the largest ones for which comparison with high-resolution rotational spectra is possible. Since the whole computational workflow can be followed employing standard electronic structure codes, accurate results for large-sized molecules can be obtained at DFT cost also by nonspecialists.
Collapse
Affiliation(s)
- 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
|
33
|
Brás EM, Zimmermann C, Fausto R, Suhm MA. Benchmarking the anisotropy of nitroxyl radical solvation with IR spectroscopy. Phys Chem Chem Phys 2024; 26:5822-5829. [PMID: 38314587 DOI: 10.1039/d3cp05668f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2024]
Abstract
Two simple nitroxyl radicals, di-tert-butyl nitroxyl (DTBN) and 2,2,6,6-tetramethylpiperidinyloxyl (TEMPO) are solvated by one or two water, methanol, tert-butyl alcohol or phenol molecules. The resulting low temperature IR spectra of the vacuum-isolated microsolvates in the OH stretching range are assigned based on harmonic DFT predictions for closed shell solvent dimers and trimers and their offset from experiment, to minimise theory-guided assignment bias. Systematic conformational preferences for the first and second solvent molecule are observed, depending on the conformational rigidity of the radical. These assignments are collected into an experimental benchmark data set and used to assess the spectral predicting power of different DFT approaches. The goal is to find inexpensive computational methods which provide reliable spectral predictions for this poorly explored class of microsolvates.
Collapse
Affiliation(s)
- Elisa M Brás
- University of Göttingen, Institute of Physical Chemistry, Tammannstr. 6, 37077 Göttingen, Germany.
- CQC-IMS, Department of Chemistry, University of Coimbra, 3004-535 Coimbra, Portugal
| | - Charlotte Zimmermann
- University of Göttingen, Institute of Physical Chemistry, Tammannstr. 6, 37077 Göttingen, Germany.
| | - Rui Fausto
- CQC-IMS, Department of Chemistry, University of Coimbra, 3004-535 Coimbra, Portugal
- Faculty of Sciences and Letters, Department of Physics, Istanbul Kultur University, Ataköy Campus, Bakirköy 34156, Istanbul, Turkey
| | - Martin A Suhm
- University of Göttingen, Institute of Physical Chemistry, Tammannstr. 6, 37077 Göttingen, Germany.
| |
Collapse
|
34
|
Barone V. Quantum chemistry meets high-resolution spectroscopy for characterizing the molecular bricks of life in the gas-phase. Phys Chem Chem Phys 2024; 26:5802-5821. [PMID: 38099409 DOI: 10.1039/d3cp05169b] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/15/2024]
Abstract
Computation of accurate geometrical structures and spectroscopic properties of large flexible molecules in the gas-phase is tackled at an affordable cost using a general exploration/exploitation strategy. The most distinctive feature of the approach is the careful selection of different quantum chemical models for energies, geometries and vibrational frequencies with the aim of maximizing the accuracy of the overall description while retaining a reasonable cost for all the steps. In particular, a composite wave-function method is used for energies, whereas a double-hybrid functional (with the addition of core-valence correlation) is employed for geometries and harmonic frequencies and a cheaper hybrid functional for anharmonic contributions. A thorough benchmark based on a wide range of prototypical molecular bricks of life shows that the proposed strategy is close to the accuracy of state-of-the-art composite wave-function methods, and is applicable to much larger systems. A freely available web-utility post-processes the geometries optimized by standard electronic structure codes paving the way toward the accurate yet not prohibitively expensive study of medium- to large-sized molecules by experimentally-oriented researchers.
Collapse
Affiliation(s)
- Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy.
| |
Collapse
|
35
|
Dral PO, Ge F, Hou YF, Zheng P, Chen Y, Barbatti M, Isayev O, Wang C, Xue BX, Pinheiro Jr M, Su Y, Dai Y, Chen Y, Zhang L, Zhang S, Ullah A, Zhang Q, Ou Y. MLatom 3: A Platform for Machine Learning-Enhanced Computational Chemistry Simulations and Workflows. J Chem Theory Comput 2024; 20:1193-1213. [PMID: 38270978 PMCID: PMC10867807 DOI: 10.1021/acs.jctc.3c01203] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2023] [Revised: 12/29/2023] [Accepted: 01/03/2024] [Indexed: 01/26/2024]
Abstract
Machine learning (ML) is increasingly becoming a common tool in computational chemistry. At the same time, the rapid development of ML methods requires a flexible software framework for designing custom workflows. MLatom 3 is a program package designed to leverage the power of ML to enhance typical computational chemistry simulations and to create complex workflows. This open-source package provides plenty of choice to the users who can run simulations with the command-line options, input files, or with scripts using MLatom as a Python package, both on their computers and on the online XACS cloud computing service at XACScloud.com. Computational chemists can calculate energies and thermochemical properties, optimize geometries, run molecular and quantum dynamics, and simulate (ro)vibrational, one-photon UV/vis absorption, and two-photon absorption spectra with ML, quantum mechanical, and combined models. The users can choose from an extensive library of methods containing pretrained ML models and quantum mechanical approximations such as AIQM1 approaching coupled-cluster accuracy. The developers can build their own models using various ML algorithms. The great flexibility of MLatom is largely due to the extensive use of the interfaces to many state-of-the-art software packages and libraries.
Collapse
Affiliation(s)
- Pavlo O. Dral
- State
Key Laboratory of Physical Chemistry of Solid Surfaces, College of
Chemistry and Chemical Engineering, and Innovation Laboratory for
Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
- Fujian
Provincial Key Laboratory of Theoretical and Computational Chemistry, Xiamen, Fujian 361005, China
| | - Fuchun Ge
- State
Key Laboratory of Physical Chemistry of Solid Surfaces, College of
Chemistry and Chemical Engineering, and Innovation Laboratory for
Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
- Fujian
Provincial Key Laboratory of Theoretical and Computational Chemistry, Xiamen, Fujian 361005, China
| | - Yi-Fan Hou
- State
Key Laboratory of Physical Chemistry of Solid Surfaces, College of
Chemistry and Chemical Engineering, and Innovation Laboratory for
Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
- Fujian
Provincial Key Laboratory of Theoretical and Computational Chemistry, Xiamen, Fujian 361005, China
| | - Peikun Zheng
- State
Key Laboratory of Physical Chemistry of Solid Surfaces, College of
Chemistry and Chemical Engineering, and Innovation Laboratory for
Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
- Fujian
Provincial Key Laboratory of Theoretical and Computational Chemistry, Xiamen, Fujian 361005, China
| | - Yuxinxin Chen
- State
Key Laboratory of Physical Chemistry of Solid Surfaces, College of
Chemistry and Chemical Engineering, and Innovation Laboratory for
Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
- Fujian
Provincial Key Laboratory of Theoretical and Computational Chemistry, Xiamen, Fujian 361005, China
| | - Mario Barbatti
- Aix
Marseille University, CNRS, ICR, Marseille 13013, France
- Institut
Universitaire de France, Paris 75231, France
| | - Olexandr Isayev
- Department
of Chemistry, Carnegie Mellon University, Pittsburgh, Pennsylvania15213, United States
| | - Cheng Wang
- State
Key Laboratory of Physical Chemistry of Solid Surfaces, College of
Chemistry and Chemical Engineering, and Innovation Laboratory for
Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
- iChem, Xiamen University, Xiamen, Fujian 361005, China
| | - Bao-Xin Xue
- State
Key Laboratory of Physical Chemistry of Solid Surfaces, College of
Chemistry and Chemical Engineering, and Innovation Laboratory for
Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
- Fujian
Provincial Key Laboratory of Theoretical and Computational Chemistry, Xiamen, Fujian 361005, China
| | - Max Pinheiro Jr
- Aix
Marseille University, CNRS, ICR, Marseille 13013, France
| | - Yuming Su
- State
Key Laboratory of Physical Chemistry of Solid Surfaces, College of
Chemistry and Chemical Engineering, and Innovation Laboratory for
Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
- iChem, Xiamen University, Xiamen, Fujian 361005, China
| | - Yiheng Dai
- State
Key Laboratory of Physical Chemistry of Solid Surfaces, College of
Chemistry and Chemical Engineering, and Innovation Laboratory for
Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
- iChem, Xiamen University, Xiamen, Fujian 361005, China
| | - Yangtao Chen
- State
Key Laboratory of Physical Chemistry of Solid Surfaces, College of
Chemistry and Chemical Engineering, and Innovation Laboratory for
Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
- iChem, Xiamen University, Xiamen, Fujian 361005, China
| | - Lina Zhang
- State
Key Laboratory of Physical Chemistry of Solid Surfaces, College of
Chemistry and Chemical Engineering, and Innovation Laboratory for
Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
- Fujian
Provincial Key Laboratory of Theoretical and Computational Chemistry, Xiamen, Fujian 361005, China
| | - Shuang Zhang
- State
Key Laboratory of Physical Chemistry of Solid Surfaces, College of
Chemistry and Chemical Engineering, and Innovation Laboratory for
Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
- Fujian
Provincial Key Laboratory of Theoretical and Computational Chemistry, Xiamen, Fujian 361005, China
| | - Arif Ullah
- School
of Physics and Optoelectronic Engineering, Anhui University, Hefei230601, China
| | - Quanhao Zhang
- State
Key Laboratory of Physical Chemistry of Solid Surfaces, College of
Chemistry and Chemical Engineering, and Innovation Laboratory for
Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
- Fujian
Provincial Key Laboratory of Theoretical and Computational Chemistry, Xiamen, Fujian 361005, China
| | - Yanchi Ou
- State
Key Laboratory of Physical Chemistry of Solid Surfaces, College of
Chemistry and Chemical Engineering, and Innovation Laboratory for
Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
- Fujian
Provincial Key Laboratory of Theoretical and Computational Chemistry, Xiamen, Fujian 361005, China
| |
Collapse
|
36
|
Finney JM, McCoy AB. Correlations between the Structures and Spectra of Protonated Water Clusters. J Phys Chem A 2024; 128:868-879. [PMID: 38265889 DOI: 10.1021/acs.jpca.3c07338] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2024]
Abstract
Badger's rule-like correlations between OH stretching frequencies and intensities and the OH bond length are used to develop a spectral mapping procedure for studies of pure and protonated water clusters. This approach utilizes the vibrationally averaged OH bond lengths, which were obtained from diffusion Monte Carlo simulations that were performed using the general potential developed by Yu and Bowman. Good agreement is achieved between the spectra obtained using this approach and previously reported spectra for H+(H2O)n clusters, with n = 3, 4, and 5, as well as their perdeuterated analogues. The analysis of the spectra obtained by this spectral mapping approach supports previous work that assigned the spectrum of H+(H2O)6 to a mixture of Eigen and Zundel-like structures. Analysis of the calculated spectra also suggests a reassignment of the frequency of one of the transitions that involves the OH stretching vibration of the OH bonds in the hydronium core in the Eigen-like structure of H+(H2O)6 from 1917 cm-1 to roughly 2100 cm-1. For D+(D2O)6, comparison of the measured spectrum to those obtained by using the spectral mapping approach suggests that the carrier of the measured spectrum is one or more of the isomers of D+(D2O)6 that contain a four-membered ring and two flanking water molecules. While there are several candidate structures, the two flanking water molecules most likely form a chain that is bound to the hydronium core.
Collapse
Affiliation(s)
- Jacob M Finney
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Anne B McCoy
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
37
|
Hui WCH, Lemke KH. The sulfur dioxide-water complex: CCSD(T)/CBS anharmonic vibrational spectroscopy of stacked and hydrogen-bonded dimers. J Chem Phys 2024; 160:054307. [PMID: 38341706 DOI: 10.1063/5.0177077] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2023] [Accepted: 01/11/2024] [Indexed: 02/13/2024] Open
Abstract
This study examines the structures, energies, and IR vibrational spectra of the sulfur dioxide-water SO2(H2O) complexes by employing coupled cluster theory CCSD(T) with Dunning style correlation consistent type basis sets aug-cc-pV(n+d)Z (n = D, T, Q, 5). Complete basis set (CBS) extrapolations have been carried out to predict binding energies for two isomers of the SO2(H2O) complex: a stacked global minimum (1A) structure and a hydrogen-bonded local minimum (1B) structure. The CCSD(T)/CBS extrapolation predicts an intermolecular S-O distance rS⋯O = 2.827 Å for the stacked isomer, which is in excellent agreement with an experimental measurement of 2.824 Å [K. Matsumura et al., J. Chem. Phys., 91, 5887 (1989)]. The CCSD(T)/CBS binding energy for the stacked dimer 1A and hydrogen-bonded form 1B is De = -4.37 kcal/mol and De = -2.40 kcal/mol, respectively. This study also employs anharmonic VPT2 MP2/aug-cc-pV(n+d)Z level corrections to CCSD(T)/aug-cc-pV(n+d)Z vibrational frequencies in both forms of SO2(H2O). The anharmonic CCSD(T)/aug-cc-pV(Q+d)Z OH stretching frequencies in the stacked structure 1A are 3743 cm-1 (ν3) and 3647 cm-1 (ν1), and these align well with the recorded IR spectroscopic values of 3745 and 3643 cm-1, respectively [C. Wang et al., J. Phys. Chem. Lett., 13, 5654 (2022)]. If we combine CCSD(T)/aug-cc-pV(n+d)Z De values with VPT2 vibrational frequencies, we obtain a new CCSD(T)/aug-cc-pV(Q+d)Z anharmonic dissociation energy D0 = -3.48 kcal/mol for 1A and D0 = -1.74 kcal/mol for 1B. In summary, the results presented here demonstrate that the application of CCSD(T) calculations with aug-cc-pV(n+d)Z basis sets and CBS extrapolations is critical in probing the structure and IR spectroscopic properties of the sulfur dioxide-water complex.
Collapse
Affiliation(s)
- Wallace C H Hui
- Department of Earth Sciences, University of Hong Kong, Pokfulam Road, Hong Kong SAR
| | - Kono H Lemke
- Department of Earth Sciences, University of Hong Kong, Pokfulam Road, Hong Kong SAR
| |
Collapse
|
38
|
Gleeson R, Aggelund PA, Østergaard FC, Schaltz KF, Sauer SPA. Exploring Alternate Methods for the Calculation of High-Level Vibrational Corrections of NMR Spin-Spin Coupling Constants. J Chem Theory Comput 2024. [PMID: 38299500 DOI: 10.1021/acs.jctc.3c01223] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2024]
Abstract
Traditional nuclear magnetic resonance (NMR) calculations typically treat systems with a Born-Oppenheimer-derived electronic wave function that is solved for a fixed nuclear geometry. One can numerically account for this neglected nuclear motion by averaging over property values for all nuclear geometries with a vibrational wave function and adding this expectation value as a correction to an equilibrium geometry property value. Presented are benchmark coupled-cluster singles and doubles (CCSD) vibrational corrections to spin-spin coupling constants (SSCCs) computed at the level of vibrational second-order perturbation theory (VPT2) using the vibrational averaging driver of the CFOUR program. As CCSD calculations of vibrational corrections are very costly, cheaper electronic structure methods are explored via a newly developed Python vibrational averaging program within the Dalton Project. Namely, results obtained with the second-order polarization propagator approximation (SOPPA) and density functional theory (DFT) with the B3LYP and PBE0 exchange-correlation functionals are compared to the benchmark CCSD//CCSD(T) and experimental values. CCSD//CCSD(T) corrections are also combined with literature CC3 equilibrium geometry values to form the highest-order vibrationally corrected values available (i.e., CC3//CCSD(T) + CCSD//CCSD(T)). CCSD//CCSD(T) statistics showed favorable statistics in comparison to experimental values, albeit at an unfavorably high computational cost. A cheaper CCSD//CCSD(T) + B3LYP method showed quite similar mean absolute deviation (MAD) values as CCSD//CCSD(T), concluding that CCSD//CCSD(T) + B3LYP is optimal in terms of cost and accuracy. With reference to experimental values, a vibrational correction was not worth the cost for all of the other methods tested. Finally, deviation statistics showed that CC3//CCSD(T) + CCSD//CCSD(T) vibrational-corrected equilibrium values deteriorated in comparison to CCSD//CCSD(T) attributed to the use of a smaller basis set or lack of solvation effects for the CC3 equilibrium calculations.
Collapse
Affiliation(s)
- Ronan Gleeson
- Department of Chemistry, University of Copenhagen, 2100 Copenhagen, Denmark
| | - Patrick A Aggelund
- Department of Chemistry, University of Copenhagen, 2100 Copenhagen, Denmark
| | | | - Kasper F Schaltz
- Department of Chemistry, University of Copenhagen, 2100 Copenhagen, Denmark
| | - Stephan P A Sauer
- Department of Chemistry, University of Copenhagen, 2100 Copenhagen, Denmark
| |
Collapse
|
39
|
Zhu GZ, Lao G, Dickerson CE, Caram JR, Campbell WC, Alexandrova AN, Hudson ER. Extending the Large Molecule Limit: The Role of Fermi Resonance in Developing a Quantum Functional Group. J Phys Chem Lett 2024; 15:590-597. [PMID: 38198595 DOI: 10.1021/acs.jpclett.3c03177] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2024]
Abstract
Polyatomic molecules equipped with optical cycling centers (OCCs), enabling continuous photon scattering during optical excitation, are exciting candidates for advancing quantum information science. However, as these molecules grow in size and complexity, the interplay of complex vibronic couplings on optical cycling becomes a critical but relatively unexplored consideration. Here, we present an extensive exploration of Fermi resonances in large-scale OCC-containing molecules using high-resolution dispersed laser-induced fluorescence and excitation spectroscopy. These resonances manifest as vibrational coupling leading to intensity borrowing by combination bands near optically active harmonic bands, which require additional repumping lasers for effective optical cycling. To mitigate these effects, we explore altering the vibrational energy level spacing through substitutions on the phenyl ring or changes in the OCC itself. While the complete elimination of vibrational coupling in complex molecules remains challenging, our findings highlight significant mitigation possibilities, opening new avenues for optimizing optical cycling in large polyatomic molecules.
Collapse
Affiliation(s)
- Guo-Zhu Zhu
- Department of Physics and Astronomy, University of California, Los Angeles, California 90095, United States
| | - Guanming Lao
- Department of Physics and Astronomy, University of California, Los Angeles, California 90095, United States
| | - Claire E Dickerson
- Department of Chemistry and Biochemistry, University of California, Los Angeles, California 90095, United States
| | - Justin R Caram
- Department of Chemistry and Biochemistry, University of California, Los Angeles, California 90095, United States
- Center for Quantum Science and Engineering, University of California, Los Angeles, California 90095, United States
| | - Wesley C Campbell
- Department of Physics and Astronomy, University of California, Los Angeles, California 90095, United States
- Center for Quantum Science and Engineering, University of California, Los Angeles, California 90095, United States
- Challenge Institute for Quantum Computation, University of California, Los Angeles, California 90095, United States
| | - Anastassia N Alexandrova
- Department of Chemistry and Biochemistry, University of California, Los Angeles, California 90095, United States
- Center for Quantum Science and Engineering, University of California, Los Angeles, California 90095, United States
| | - Eric R Hudson
- Department of Physics and Astronomy, University of California, Los Angeles, California 90095, United States
- Center for Quantum Science and Engineering, University of California, Los Angeles, California 90095, United States
- Challenge Institute for Quantum Computation, University of California, Los Angeles, California 90095, United States
| |
Collapse
|
40
|
Zasimov PV, Volosatova AD, Góbi S, Keresztes B, Tyurin DA, Feldman VI, Tarczay G. Infrared spectroscopy of the α-hydroxyethyl radical isolated in cryogenic solid media. J Chem Phys 2024; 160:024308. [PMID: 38205854 DOI: 10.1063/5.0177189] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2023] [Accepted: 12/20/2023] [Indexed: 01/12/2024] Open
Abstract
The α-hydroxyethyl radical (CH3·CHOH, 2A) is a key intermediate in ethanol biochemistry, combustion, atmospheric chemistry, radiation chemistry, and astrochemistry. Experimental data on the vibrational spectrum of this radical are crucially important for reliable detection and understanding of the chemical dynamics of this species. This study represents the first detailed experimental report on the infrared absorption bands of the α-hydroxyethyl radical complemented by ab initio computations. The radical was generated in solid para-H2 and Xe matrices via the reactions of hydrogen atoms with matrix-isolated ethanol molecules and radiolysis of isolated ethanol molecules with x rays. The absorption bands with maxima at 3654.6, 3052.1, 1425.7, 1247.9, 1195.6 (1177.4), and 1048.4 cm-1, observed in para-H2 matrices appearing upon the H· atom reaction, were attributed to the OHstr, α-CHstr, CCstr, COstr + CCObend, COstr, and CCstr + CCObend vibrational modes of the CH3·CHOH radical, respectively. The absorption bands with the positions slightly red-shifted from those observed in para-H2 were detected in both the irradiated and post-irradiation annealed Xe matrices containing C2H5OH. The results of the experiments with the isotopically substituted ethanol molecules (CH3CD2OH and CD3CD2OH) and the quantum-chemical computations at the UCCSD(T)/L2a_3 level support the assignment. The photolysis with ultraviolet light (240-300 nm) results in the decay of the α-hydroxyethyl radical, yielding acetaldehyde and its isomer, vinyl alcohol. A comparison of the experimental and theoretical results suggests that the radical adopts the thermodynamically more stable anti-conformation in both matrices.
Collapse
Affiliation(s)
- Pavel V Zasimov
- MTA-ELTE Lendület Laboratory Astrochemistry Research Group, Institute of Chemistry, ELTE Eötvös Loránd University, P.O. Box 32, H-1518 Budapest, Hungary
- Department of Chemistry, Lomonosov Moscow State University, 119991 Moscow, Russia
| | - Anastasia D Volosatova
- MTA-ELTE Lendület Laboratory Astrochemistry Research Group, Institute of Chemistry, ELTE Eötvös Loránd University, P.O. Box 32, H-1518 Budapest, Hungary
- Department of Chemistry, Lomonosov Moscow State University, 119991 Moscow, Russia
| | - Sándor Góbi
- MTA-ELTE Lendület Laboratory Astrochemistry Research Group, Institute of Chemistry, ELTE Eötvös Loránd University, P.O. Box 32, H-1518 Budapest, Hungary
| | - Barbara Keresztes
- Laboratory of Molecular Spectroscopy, Institute of Chemistry, ELTE Eötvös Loránd University, P.O. Box 32, H-1518 Budapest, Hungary
- Hevesy György PhD School of Chemistry, Institute of Chemistry, ELTE Eötvös Loránd University, P.O. Box 32, H-1518 Budapest, Hungary
| | - Daniil A Tyurin
- Department of Chemistry, Lomonosov Moscow State University, 119991 Moscow, Russia
| | - Vladimir I Feldman
- Department of Chemistry, Lomonosov Moscow State University, 119991 Moscow, Russia
| | - György Tarczay
- MTA-ELTE Lendület Laboratory Astrochemistry Research Group, Institute of Chemistry, ELTE Eötvös Loránd University, P.O. Box 32, H-1518 Budapest, Hungary
- Laboratory of Molecular Spectroscopy, Institute of Chemistry, ELTE Eötvös Loránd University, P.O. Box 32, H-1518 Budapest, Hungary
- Centre for Astrophysics and Space Science, ELTE Eötvös Loránd University, P.O. Box 32, H-1518 Budapest, Hungary
| |
Collapse
|
41
|
Cong F, Cai L, Cheng J, Pu Z, Wang X. Beryllium Dimer Reactions with Acetonitrile: Formation of Strong Be-Be Bonds. Molecules 2023; 29:177. [PMID: 38202759 PMCID: PMC10779904 DOI: 10.3390/molecules29010177] [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: 11/11/2023] [Revised: 12/25/2023] [Accepted: 12/26/2023] [Indexed: 01/12/2024] Open
Abstract
Laser ablated Be atoms have been reacted with acetonitrile molecules in 4 K solid neon matrix. The diberyllium products BeBeNCCH3 and CNBeBeCH3 have been identified by D and 13C isotopic substitutions and quantum chemical calculations. The stabilization of the diberyllium species is rationalized from the formation of the real Be-Be single bonds with bond distances as 2.077 and 2.058 Å and binding energies as -27.1 and -77.2 kcal/mol calculated at CCSD (T)/aug-cc-pVTZ level of theory for BeBeNCCH3 and CNBeBeCH3, respectively. EDA-NOCV analysis described the interaction between Be2 and NC···CH3 fragments as Lewis "acid-base" interactions. In the complexes, the Be2 moiety carries positive charges which transfer from antibonding orbital of Be2 to the bonding fragments significantly strengthen the Be-Be bonds that are corroborated by AIM, LOL and NBO analyses. In addition, mono beryllium products BeNCCH3, CNBeCH3, HBeCH2CN and HBeNCCH2 have also been observed in our experiments.
Collapse
Affiliation(s)
- Fei Cong
- Shanghai Key Lab of Chemical Assessment and Sustainability, School of Chemical Science and Engineering, Tongji University, Shanghai 200092, China; (F.C.); (L.C.); (J.C.)
| | - Liyan Cai
- Shanghai Key Lab of Chemical Assessment and Sustainability, School of Chemical Science and Engineering, Tongji University, Shanghai 200092, China; (F.C.); (L.C.); (J.C.)
| | - Juanjuan Cheng
- Shanghai Key Lab of Chemical Assessment and Sustainability, School of Chemical Science and Engineering, Tongji University, Shanghai 200092, China; (F.C.); (L.C.); (J.C.)
| | - Zhen Pu
- China Academy of Engineering and Physics, Mianyang 621900, China
| | - Xuefeng Wang
- Shanghai Key Lab of Chemical Assessment and Sustainability, School of Chemical Science and Engineering, Tongji University, Shanghai 200092, China; (F.C.); (L.C.); (J.C.)
| |
Collapse
|
42
|
Zhang Y, Wang Y, Xu X, Chen Z, Yang Y. Vibrational Spectra of Highly Anharmonic Water Clusters: Molecular Dynamics and Harmonic Analysis Revisited with Constrained Nuclear-Electronic Orbital Methods. J Chem Theory Comput 2023; 19:9358-9368. [PMID: 38096546 DOI: 10.1021/acs.jctc.3c01037] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2023]
Abstract
Vibrational spectroscopy is widely used to gain insights into structural and dynamic properties of chemical, biological, and materials systems. Thus, an efficient and accurate method to simulate vibrational spectra is desired. In this paper, we justify and employ a microcanonical molecular simulation scheme to calculate the vibrational spectra of three challenging water clusters: the neutral water dimer (H4O2), the protonated water trimer (H7O3+), and the protonated water tetramer (H9O4+). We find that with the accurate description of quantum nuclear delocalization effects through the constrained nuclear-electronic orbital framework, including vibrational mode coupling effects through molecular dynamics simulations can additionally improve the vibrational spectrum calculations. In contrast, without the quantum nuclear delocalization picture, conventional ab initio molecular dynamics may even lead to less accurate results than harmonic analysis.
Collapse
Affiliation(s)
- Yuzhe Zhang
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Yiwen Wang
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Xi Xu
- Center for Advanced Materials Research, Beijing Normal University, Zhuhai 519087, China
| | - Zehua Chen
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Yang Yang
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| |
Collapse
|
43
|
Greis K, Kirschbaum C, Ober K, Taccone MI, Torres-Boy A, Meijer G, Pagel K, von Helden G. Infrared Spectroscopy of Fluorenyl Cations at Cryogenic Temperatures. J Phys Chem Lett 2023; 14:11313-11317. [PMID: 38064287 PMCID: PMC10749476 DOI: 10.1021/acs.jpclett.3c02928] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2023] [Revised: 12/01/2023] [Accepted: 12/05/2023] [Indexed: 12/22/2023]
Abstract
The notion of (anti)aromaticity is a successful concept in chemistry to explain the structure and stability of polycyclic hydrocarbons. Cyclopentadienyl and fluorenyl cations are among the most studied classical antiaromatic systems. In this work, fluorenyl cations are investigated by high-resolution gas-phase infrared spectroscopy in helium droplets. Bare fluorenyl cations are generated in the gas phase by electrospray ionization. After mass-to-charge selection, ions are captured in ultracold helium nanodroplets and probed by infrared spectroscopy using a widely tunable free-electron laser in the 600-1700 cm-1 range. The highly resolved cryogenic infrared spectra confirm, in combination with DFT computations, that all cations are present in their singlet states.
Collapse
Affiliation(s)
- Kim Greis
- Fritz
Haber Institute of the Max Planck Society, Berlin, Faradayweg 4-6, 14195 Berlin, Germany
- Institute
of Chemistry and Biochemistry, Freie Universität
Berlin, Altensteinstraße 23A, 14195 Berlin, Germany
| | - Carla Kirschbaum
- Fritz
Haber Institute of the Max Planck Society, Berlin, Faradayweg 4-6, 14195 Berlin, Germany
- Institute
of Chemistry and Biochemistry, Freie Universität
Berlin, Altensteinstraße 23A, 14195 Berlin, Germany
| | - Katja Ober
- Fritz
Haber Institute of the Max Planck Society, Berlin, Faradayweg 4-6, 14195 Berlin, Germany
| | - Martín I. Taccone
- Fritz
Haber Institute of the Max Planck Society, Berlin, Faradayweg 4-6, 14195 Berlin, Germany
| | - América
Y. Torres-Boy
- Fritz
Haber Institute of the Max Planck Society, Berlin, Faradayweg 4-6, 14195 Berlin, Germany
| | - Gerard Meijer
- Fritz
Haber Institute of the Max Planck Society, Berlin, Faradayweg 4-6, 14195 Berlin, Germany
| | - Kevin Pagel
- Fritz
Haber Institute of the Max Planck Society, Berlin, Faradayweg 4-6, 14195 Berlin, Germany
- Institute
of Chemistry and Biochemistry, Freie Universität
Berlin, Altensteinstraße 23A, 14195 Berlin, Germany
| | - Gert von Helden
- Fritz
Haber Institute of the Max Planck Society, Berlin, Faradayweg 4-6, 14195 Berlin, Germany
| |
Collapse
|
44
|
Barone V, Lazzari F. Hunting for Complex Organic Molecules in the Interstellar Medium: The Role of Accurate Low-Cost Theoretical Geometries and Rotational Constants. J Phys Chem A 2023; 127:10517-10527. [PMID: 38033327 PMCID: PMC10726368 DOI: 10.1021/acs.jpca.3c06649] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2023] [Revised: 11/07/2023] [Accepted: 11/08/2023] [Indexed: 12/02/2023]
Abstract
A new approach to computation at affordable cost of accurate geometrical structures and rotational constants for medium-sized molecules in the gas phase is further improved and applied to a large panel of interstellar complex organic molecules. The most distinctive feature of the new model is the effective inclusion of core-valence correlation and vibrational averaging effects in the framework of density functional theory (DFT). In particular, a double-hybrid functional in conjunction with a quadruple-ζ valence/triple-ζ polarization basis set is employed for geometry optimizations, whereas a cheaper hybrid functional in conjunction with a split-valence basis set is used for the evaluation of vibrational corrections. A thorough benchmark based on a wide range of prototypical systems shows that the new scheme approaches the accuracy of state-of-the-art wave function methods with the computational cost of the standard methods (DFT or MP2) routinely employed in the interpretation of microwave spectra. Since the whole computational workflow involves the postprocessing of the output of standard electronic structure codes by a new freely available web utility, the way is paved for the accurate yet not prohibitively expensive study of medium- to large-sized molecules also by nonspecialists.
Collapse
Affiliation(s)
- Vincenzo Barone
- Scuola Normale Superiore, di Pisa, Piazza dei Cavalieri 7, Pisa 56125, Italy
| | - Federico Lazzari
- Scuola Normale Superiore, di Pisa, Piazza dei Cavalieri 7, Pisa 56125, Italy
| |
Collapse
|
45
|
Melosso M, Alessandrini S, Spada L, Melli A, Wang X, Zheng Y, Duan C, Li J, Du W, Gou Q, Bizzocchi L, Dore L, Barone V, Puzzarini C. Rotational spectra and semi-experimental structures of furonitrile and its water cluster. Phys Chem Chem Phys 2023; 25:31281-31291. [PMID: 37955344 DOI: 10.1039/d3cp03984f] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2023]
Abstract
Rotational spectroscopy represents an invaluable tool for several applications: from the identification of new molecules in interstellar objects to the characterization of van der Waals complexes, but also for the determination of very accurate molecular structures and for conformational analyses. In this work, we used high-resolution rotational spectroscopic techniques in combination with high-level quantum-chemical calculations to address all these aspects for two isomers of cyanofuran, namely 2-furonitrile and 3-furonitrile. In particular, we have recorded and analyzed the rotational spectra of both of them from 6 to 320 GHz; rotational transitions belonging to several singly-substituted isotopologues have been identified as well. The rotational constants derived in this way have been used in conjunction with computed rotation-vibration interaction constants in order to derive a semi-experimental equilibrium structure for both isomers. Moreover, we observed the rotational spectra of four different intermolecular adducts formed by furonitrile and water, whose identification has been supported by a conformational analysis and a theoretical spectroscopic characterization. A semi-experimental determination of the intermolecular parameters has been achieved for all of them and the results have been compared with those obtained for the analogous system formed by benzonitrile and water.
Collapse
Affiliation(s)
- Mattia Melosso
- 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.
| | - Lorenzo Spada
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - 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
| | - Xiujuan Wang
- Department of Chemistry, School of Chemistry and Chemical Engineering, Chongqing University, Daxuecheng South Rd. 55, 401331, Chongqing, China.
| | - Yang Zheng
- Department of Chemistry, School of Chemistry and Chemical Engineering, Chongqing University, Daxuecheng South Rd. 55, 401331, Chongqing, China.
| | - Chunguo Duan
- Department of Chemistry, School of Chemistry and Chemical Engineering, Chongqing University, Daxuecheng South Rd. 55, 401331, Chongqing, China.
| | - Jiayi Li
- Department of Chemistry, School of Chemistry and Chemical Engineering, Chongqing University, Daxuecheng South Rd. 55, 401331, Chongqing, China.
| | - Weiping Du
- Department of Chemistry, School of Chemistry and Chemical Engineering, Chongqing University, Daxuecheng South Rd. 55, 401331, Chongqing, China.
| | - Qian Gou
- Department of Chemistry, School of Chemistry and Chemical Engineering, Chongqing University, Daxuecheng South Rd. 55, 401331, Chongqing, China.
| | - Luca Bizzocchi
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy.
| | - Luca Dore
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy.
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Cristina Puzzarini
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy.
| |
Collapse
|
46
|
Tran HK, Berkelbach TC. Vibrational heat-bath configuration interaction with semistochastic perturbation theory using harmonic oscillator or VSCF modals. J Chem Phys 2023; 159:194101. [PMID: 37965997 PMCID: PMC10653875 DOI: 10.1063/5.0172702] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2023] [Accepted: 10/23/2023] [Indexed: 11/16/2023] Open
Abstract
Vibrational heat-bath configuration interaction (VHCI)-a selected configuration interaction technique for vibrational structure theory-has recently been developed in two independent works [J. H. Fetherolf and T. C. Berkelbach, J. Chem. Phys. 154, 074104 (2021); A. U. Bhatty and K. R. Brorsen, Mol. Phys. 119, e1936250 (2021)], where it was shown to provide accuracy on par with the most accurate vibrational structure methods with a low computational cost. Here, we eliminate the memory bottleneck of the second-order perturbation theory correction using the same (semi)stochastic approach developed previously for electronic structure theory. This allows us to treat, in an unbiased manner, much larger perturbative spaces, which are necessary for high accuracy in large systems. Stochastic errors are easily controlled to be less than 1 cm-1. We also report two other developments: (i) we propose a new heat-bath criterion and an associated exact implicit sorting algorithm for potential energy surfaces expressible as a sum of products of one-dimensional potentials; (ii) we formulate VHCI to use a vibrational self-consistent field (VSCF) reference, as opposed to the harmonic oscillator reference configuration used in previous reports. Our tests are done with quartic and sextic force fields, for which we find that with VSCF, the minor improvements to accuracy are outweighed by the higher computational cost associated the matrix element evaluations. We expect VSCF-based VHCI to be important for more general potential representations, for which the harmonic oscillator basis function integrals are no longer analytic.
Collapse
Affiliation(s)
- Henry K. Tran
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | | |
Collapse
|
47
|
Crisci L, Di Grande S, Cavallotti C, Barone V. Toward an Accurate Black-Box Tool for the Kinetics of Gas-Phase Reactions Involving Barrier-less Elementary Steps. J Chem Theory Comput 2023; 19:7626-7639. [PMID: 37880932 PMCID: PMC10653117 DOI: 10.1021/acs.jctc.3c00857] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2023] [Revised: 10/06/2023] [Accepted: 10/13/2023] [Indexed: 10/27/2023]
Abstract
An enhanced computational protocol has been devised for the accurate characterization of gas-phase barrier-less reactions in the framework of the reaction-path (RP) and variable reaction coordinate variational transition-state theory. In particular, the synergistic combination of density functional theory and Monte Carlo sampling to optimize reactive fluxes led to a reliable yet effective computational workflow. A black-box strategy has been developed for selecting the most suited density functional with reference to a high-level one-dimensional reference potential. At the same time, different descriptions of hindered rotations are automatically selected, depending on the corresponding harmonic frequencies along the RP. The performance of the new tool is investigated by means of two prototypical reactions involving different degrees of static and dynamic correlation, namely, H2S + Cl and CH3 + CH3. The remarkable agreement of the computed kinetic parameters with the available experimental data confirms the accuracy and robustness of the proposed approach. Together with their intrinsic interest, these results also pave the way toward systematic investigations of gas-phase reactions involving barrier-less elementary steps by a reliable, user-friendly tool, which can be confidently used also by nonspecialists.
Collapse
Affiliation(s)
- Luigi Crisci
- Scuola
Normale Superiore di Pisa, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
- Department
of Chemical Sciences, University of Napoli
Federico II, Complesso Universitario di M.S. Angelo, via Cintia 21, I-80126 Napoli, Italy
| | - Silvia Di Grande
- Scuola
Normale Superiore di Pisa, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
- Scuola
Superiore Meridionale, Largo San Marcellino 10, I-80138 Napoli, Italy
| | - Carlo Cavallotti
- Department
of Chemistry, Materials and Chemical Engineering “G. Natta”, Politecnico di Milano, I-20131 Milano, Italy
| | - Vincenzo Barone
- Scuola
Normale Superiore di Pisa, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| |
Collapse
|
48
|
Chang JL, Chen HY, Huang YJ. Reassignment of the Photoelectron Spectrum of Methylketene Using a Hybrid Model of Harmonic and Anharmonic Oscillators to Compute Franck-Condon Factors. ACS OMEGA 2023; 8:40685-40694. [PMID: 37953835 PMCID: PMC10637627 DOI: 10.1021/acsomega.3c05750] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/05/2023] [Accepted: 10/02/2023] [Indexed: 11/14/2023]
Abstract
We constructed a hybrid model of harmonic and anharmonic oscillators to compute Franck-Condon factors and interpret the photoelectron spectrum of methylketene. The equilibrium structures of methylketene and its cation were optimized, and then, the harmonic and anharmonic vibrational frequencies were computed using the B3LYP, PBE0, APFD, and ωB97XD approaches of the density functional theory. The photoelectron spectrum of methylketene was simulated by computing the Franck-Condon factors with both the harmonic and hybrid models. The adiabatic ionization energy of methylketene was computed by using the CCSD(T) approach extrapolating to the complete basis set limit. The simulated photoelectron spectra are consistent with those from the experiment for both the harmonic and hybrid models. However, the error in band positions is reduced by using the hybrid model. The computed adiabatic ionization energies of methylketene are in agreement with the experiment, with the smallest error being 0.017 eV. Our interpretation based on the theoretical spectrum led to the reassignment of the experimental photoelectron spectrum of methylketene.
Collapse
Affiliation(s)
- Jia-Lin Chang
- Department of Science Education
and Application, National Taichung University
of Education, Taichung 403514, Taiwan,
Republic of China
| | - Hsiang-Yu Chen
- Department of Science Education
and Application, National Taichung University
of Education, Taichung 403514, Taiwan,
Republic of China
| | - Yun-Jhu Huang
- Department of Science Education
and Application, National Taichung University
of Education, Taichung 403514, Taiwan,
Republic of China
| |
Collapse
|
49
|
Perera SM, Aikawa T, Shaner SE, Moran SD, Wang L. Effects of the Intramolecular Group and Solvent on Vibrational Coupling Modes and Strengths of Fermi Resonances in Aryl Azides: A DFT Study of 4-Azidotoluene and 4-Azido- N-phenylmaleimide. J Phys Chem A 2023; 127:8911-8921. [PMID: 37819373 DOI: 10.1021/acs.jpca.3c06312] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/13/2023]
Abstract
The high transition dipole strength of the azide asymmetric stretch makes aryl azides good candidates as vibrational probes (VPs). However, aryl azides have complex absorption profiles due to Fermi resonances (FRs). Understanding the origin and the vibrational modes involved in FRs of aryl azides is critically important toward developing them as VPs for studies of protein structures and structural changes in response to their surroundings. As such, we studied vibrational couplings in 4-azidotoluene and 4-azido-N-phenylmaleimide in two solvents, N,N-dimethylacetamide and tetrahydrofuran, to explore the origin and the effects of intramolecular group and solvent on the FRs of aryl azides using density functional theory (DFT) calculations with the B3LYP functional and seven basis sets, 6-31G(d,p), 6-31+G(d,p), 6-31++G(d,p), 6-311G(d,p), 6-311+G(d,p), 6-311++G(d,p), and 6-311++G(df,pd). Two combination bands consisting of the azide symmetric stretch and another mode form strong FRs with the azide asymmetric stretch for both molecules. The FR profile was altered by replacing the methyl group with maleimide. Solvents change the relative peak position and intensity more significantly for 4-azido-N-phenylmaleimide, which makes it a more sensitive VP. Furthermore, the DFT results indicate that a comparison among the results from different basis sets can be used as a means to predict more reliable vibrational spectra.
Collapse
Affiliation(s)
- Sathya M Perera
- School of Chemical and Biomolecular Sciences, Southern Illinois University, Carbondale, Illinois 62901, United States
| | - Tenyu Aikawa
- School of Chemical and Biomolecular Sciences, Southern Illinois University, Carbondale, Illinois 62901, United States
| | - Sarah E Shaner
- Department of Chemistry and Physics, Southeast Missouri State University, Cape Girardeau, Missouri 63701, United States
| | - Sean D Moran
- School of Chemical and Biomolecular Sciences, Southern Illinois University, Carbondale, Illinois 62901, United States
| | - Lichang Wang
- School of Chemical and Biomolecular Sciences, Southern Illinois University, Carbondale, Illinois 62901, United States
| |
Collapse
|
50
|
Harold SE, Warf SL, Shields GC. Prebiotic dimer and trimer peptide formation in gas-phase atmospheric nanoclusters of water. Phys Chem Chem Phys 2023; 25:28517-28532. [PMID: 37847315 DOI: 10.1039/d3cp02915h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2023]
Abstract
Insight into the origin of prebiotic molecules is key to our understanding of how living systems evolved into the complex network of biological processes on Earth. By modelling diglycine and triglycine peptide formation in the prebiotic atmosphere, we provide a plausible pathway for peptide growth. By examining different transition states (TSs), we conclude that the formation of diglycine and triglycine in atmospheric nanoclusters of water in the prebiotic atmosphere kinetically favors peptide growth by an N-to-C synthesis of glycines through a trans conformation. Addition of water stabilizes the TS structures and lowers the Gibbs free activation energies. At temperatures that model the prebiotic atmosphere, the free energies of activation with a six water nanocluster as part of the TS are predicted to be 16 kcal mol-1 relative to the prereactive complex. Examination of the trans vs. cis six water transition states reveals that a homodromic water network that maximizes the acceptor/donor nature of the six waters is responsible for enhanced kinetic favorability of the trans N-to-C pathway. Compared to the non-hydrated trans TS, the trans six-water TS accelerates the reaction of diglycine and glycine to form triglycine by 13 orders of magnitude at 217 K. Nature uses the trans N-to-C pathway to synthesize proteins in the ribosome, and we note the similarities in hydrogen bond stabilization between the transition state for peptide synthesis in the ribosome and the transition states formed in nanoclusters of water in the same pathway. These results support the hypothesis that small oligomers formed in the prebiotic atmosphere and rained onto earth's surface.
Collapse
Affiliation(s)
- Shannon E Harold
- Department of Chemistry, Furman University, Greenville, South Carolina 29613, USA.
| | - Skyler L Warf
- Department of Chemistry, Furman University, Greenville, South Carolina 29613, USA.
| | - George C Shields
- Department of Chemistry, Furman University, Greenville, South Carolina 29613, USA.
| |
Collapse
|