1
|
Alessandrini S, Melosso M, Bizzocchi L, Barone V, Puzzarini C. The Semiexperimental Approach at Work: Equilibrium Structure of Radical Species. J Phys Chem A 2024. [PMID: 38991181 DOI: 10.1021/acs.jpca.4c01758] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/13/2024]
Abstract
The so-called semiexperimental (SE) approach is a powerful technique for obtaining highly accurate equilibrium structures for isolated systems. This Featured Article describes its extension to open-shell species, thus providing the first systematic investigation on radical equilibrium geometries to be used for benchmarking purposes. The small yet significant database obtained demonstrates that there is no reduction in accuracy when moving from closed-shell species to radicals. We also provide an extension of the applicability of the SE approach to medium-/large-sized radicals by exploiting the so-called "Lego-brick" approach, which is based on the assumption that a molecular system can be seen as formed by smaller fragments for which the SE equilibrium structure is available. In this Featured Article we show that this model can be successfully applied also to open-shell species.
Collapse
Affiliation(s)
- Silvia Alessandrini
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, I-40126 Bologna, Italy
| | - Mattia Melosso
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, I-40126 Bologna, Italy
| | - Luca Bizzocchi
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, I-40126 Bologna, Italy
| | | | - Cristina Puzzarini
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, I-40126 Bologna, Italy
| |
Collapse
|
2
|
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
|
3
|
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
|
4
|
Fortenberry RC. Picking up Good Vibrations through Quartic Force Fields and Vibrational Perturbation Theory. J Phys Chem Lett 2024; 15:6528-6537. [PMID: 38875074 DOI: 10.1021/acs.jpclett.4c01089] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/16/2024]
Abstract
Quartic force fields (QFFs) define sparse potential energy surfaces (compared to semiglobal surfaces) that are the cheapest and easiest means of computing anharmonic vibrational frequencies, especially when utilized with second-order vibrational perturbation theory (VPT2). However, flat and shallow potential surfaces are exceedingly difficult for QFFs to treat through a combination of numerical noise in the often numerically computed derivatives and in competing energy factors in the composite energies often utilized to provide high-level spectroscopic predictions. While some of these issues can be alleviated with analytic derivatives, hybrid QFFs, and intelligent choices in coordinate systems, the best practice is for predicting good molecular vibrations via QFFs is to understand what they cannot do, and this manuscript documents such cases where QFFs may fail.
Collapse
Affiliation(s)
- Ryan C Fortenberry
- Department of Chemistry & Biochemistry, University of Mississippi, University, Mississippi 38677-1848, United States
| |
Collapse
|
5
|
Di Grande S, Barone V. Toward Accurate Quantum Chemical Methods for Molecules of Increasing Dimension: The New Family of Pisa Composite Schemes. J Phys Chem A 2024; 128:4886-4900. [PMID: 38847454 DOI: 10.1021/acs.jpca.4c01673] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/21/2024]
Abstract
The new versions of the Pisa composite scheme introduced in the present paper are based on 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, the computation of accurate electronic energies has been further improved introducing more reliable complete basis set extrapolations and estimation of core-valence correlation, together with improved basis sets for third-row atoms. Furthermore, the reduced-cost frozen natural orbital (FNO) model has been introduced and validated for large molecules. Accurate molecular structures can be obtained avoiding complete basis set extrapolation and evaluating core-valence correlation at the MP2 level. Unfortunately, analytical gradients are not available for the FNO version of the model. Therefore, for large molecules, an accurate reduced-cost alternative is offered by evaluation of valence contributions with a double-hybrid functional in conjunction with the same MP2 contribution for core-valence correlation or by means of a one-parameter approximation. The same double-hybrid functional and basis set are employed to evaluate zero-point energies and partition functions. After the validation of the new models for small systems, a panel of molecular bricks of life has been used to analyze their performances for problems of current fundamental or technological interest. The fully black-box implementation of the computational workflow paves the way toward the accurate yet not prohibitively expensive study of medium- to large-sized molecules also by experimentally oriented researchers.
Collapse
Affiliation(s)
- 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
| | | |
Collapse
|
6
|
Puzzarini C, Ye H, Alessandrini S. Isomerism of CH 2 SO : Accurate structural, energetic, and spectroscopic characterization. J Comput Chem 2024; 45:777-786. [PMID: 38116807 DOI: 10.1002/jcc.27283] [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: 10/17/2023] [Revised: 11/24/2023] [Accepted: 11/28/2023] [Indexed: 12/21/2023]
Abstract
A recent work [Ye et al. Mon. Not. R. Astron. Soc. 2023, 525, 1158] on the gas-phase formation of t-HC(O)SH, already detected in the interstellar medium, pointed out that the trans form of HC(S)OH is a potential candidate for astronomical observations. Prompted by these results, theCH 2 SO family of isomers has been investigated from an energetic point of view using a double-hybrid density functional in combination with a partially augmented triple-zeta basis set. This preliminary study showed that the most stable species of the family are the cis and trans forms of HC(O)SH and HC(S)OH. For their structural and spectroscopic characterization, a composite scheme based on coupled cluster (CC) calculations that incorporates up to the quadruple excitations and accounts for the extrapolation to the complete basis set limit and core correlation effects has been employed. This approach opens to the prediction of rotational constants with an accuracy of 0.1%. A hybrid scheme, based on harmonic frequencies computed using the CC singles, doubles and a perturbative treatment of triples method (CCSD(T)) in conjunction with a quadruple-zeta basis set, allowed us to obtain fundamental vibrational frequencies with a mean absolute error of about 1%.
Collapse
Affiliation(s)
- Cristina Puzzarini
- Department of Chemistry "Giacomo Ciamician", University of Bologna, Bologna, Italy
| | - Hexu Ye
- Department of Chemistry "Giacomo Ciamician", University of Bologna, Bologna, Italy
| | - Silvia Alessandrini
- Department of Chemistry "Giacomo Ciamician", University of Bologna, Bologna, Italy
| |
Collapse
|
7
|
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
|
8
|
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
|
9
|
Fortenberry RC. Quantum Chemistry and Astrochemistry: A Match Made in the Heavens. J Phys Chem A 2024; 128:1555-1565. [PMID: 38381079 DOI: 10.1021/acs.jpca.3c07601] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/22/2024]
Abstract
Quantum chemistry can uniquely answer astrochemical questions that no other technique can provide. Computations can be parallelized, automated, and left to run continuously providing exceptional molecular throughput that cannot be done through experimentation. Additionally, the granularity of the individual computations that are required of potential energy surfaces, reaction mechanism pathways, or other quantum chemically derived observables produces a unique mosaic that make up the larger whole. These pieces can be dissected for their individual contributions or evaluated in an ad hoc fashion for each of their roles in generating the larger whole. No other scientific approach is capable of reporting such fine-grained insights. Quantum chemistry also works from a bottom-up approach in providing properties directly from the desired molecule instead of a top-down perspective as required of experiment where molecules have to be linked to observed phenomena. Furthermore, modern quantum chemistry is well within the range of "chemical accuracy" and is approaching "spectroscopic accuracy." As such, the seemingly difficult questions asked by astrochemistry that would not be asked initially for any other application require quantum chemical reference data. While the results of quantum chemical computations are needed to interpret astrochemical observation, modeling, or laboratory experimentation, such hard questions, regardless of the original need to answer them, produce unique solutions. While questions in astrochemistry often require novel developments in and implementations of quantum chemistry as outlined herein, the applications of these solutions will stretch beyond astrochemistry and may yet impact fields much closer to Earth.
Collapse
Affiliation(s)
- Ryan C Fortenberry
- Department of Chemistry & Biochemistry, University of Mississippi, Oxford, Mississippi 38677-1848, United States
| |
Collapse
|
10
|
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: 0] [Impact Index Per Article: 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
|
11
|
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: 6] [Impact Index Per Article: 6.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
|
12
|
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
|
13
|
Puzzarini C, Linguerri R, Hochlaf M. Insights into the Molecular Structure and Spectroscopic Properties of HONCO: An Accurate Ab Initio Study. J Phys Chem A 2023; 127:9502-9512. [PMID: 37922399 PMCID: PMC10658623 DOI: 10.1021/acs.jpca.3c05741] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2023] [Revised: 10/06/2023] [Accepted: 10/09/2023] [Indexed: 11/05/2023]
Abstract
In an effort to provide the first accurate structural and spectroscopic characterization of the quasi-linear chain HONCO in its electronic ground state, state-of-the-art computational approaches mainly based on coupled-cluster (CC) theory have been employed. Equilibrium geometries have been calculated by means of a composite scheme based on CC calculations that incorporates up to the quadruple excitations and accounts for the extrapolation to the complete basis set limit and core correlation effects. This approach is proven to provide molecular structures with an accuracy better than 0.001 Å and 0.05° for bond lengths and angles, respectively. Incorporation of vibrational effects permits this level of theory to predict rotational constants with an estimated accuracy of 0.1% or better. Vibrational fundamental bands have been evaluated by means of a hybrid scheme based on harmonic frequencies computed using the CC singles, doubles, and a perturbative treatment of the triples method (CCSD(T)) in conjunction with a quadruple-ζ basis set, with all electrons being correlated, and anharmonic corrections from CCSD(T) calculations using a triple-ζ basis set, within the frozen-core approximation. Such a hybrid approach allowed us to obtain fundamental frequencies with a mean absolute error of about 1%. To complete the spectroscopic characterization, vertical electronic excitation energies have been calculated for the lowest singlet and triplet states using the internally contracted multireference configuration interaction (MRCI) method. Computations show that HONCO dissociates into OH + NCO upon the absorption of UV-vis light. In conclusion, we are confident that the highly accurate spectroscopic data provided herein can be useful for guiding future experimental investigations and supporting the characterization of this molecule in atmospheric and astrophysical media, as well as in combustion.
Collapse
Affiliation(s)
- Cristina Puzzarini
- Dipartimento
di Chimica “Giacomo Ciamician″, Università di Bologna, Via Selmi 2, 40126 Bologna, Italy
| | - Roberto Linguerri
- Université
Gustave Eiffel, COSYS/IMSE, 5 Bd Descartes, 77454 Champs sur Marne, France
| | - Majdi Hochlaf
- Université
Gustave Eiffel, COSYS/IMSE, 5 Bd Descartes, 77454 Champs sur Marne, France
| |
Collapse
|
14
|
Gao T, Lei J, Zou S, Wang C, Xu X, Gou Q. Conformational equilibria and interaction preference in the complex of isoprene-maleic anhydride. Phys Chem Chem Phys 2023; 25:27798-27804. [PMID: 37814807 DOI: 10.1039/d3cp03712f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/11/2023]
Abstract
The rotational spectrum of the isoprene-maleic anhydride complex has been investigated by pulsed jet Fourier transform microwave spectroscopy and interpreted with complementary quantum chemical calculations. Theoretical predictions have yielded four plausible isomers, all residing within an energy window of 12 kJ mol-1. However, two distinct isomers characterized by a π-π stacked configuration have been experimentally observed in pulsed jets, which have differed in the orientation of isoprene over maleic anhydride. The relative population ratio of the two detected isomers has been estimated to be NI/NII ≈ 3/1 from rigorous measurements of the relative intensity on a set of μc-type transitions. Remarkably, this study underscores the pivotal role played by the interaction between the CC bonding orbital (π) of isoprene and the CC antibonding orbital (π*) of maleic anhydride in stabilizing the target complex.
Collapse
Affiliation(s)
- Tianyue Gao
- Department of Chemistry, School of Chemistry and Chemical Engineering, Chongqing University, Daxuecheng South Rd. 55, 401331 Chongqing, China.
| | - Juncheng Lei
- Department of Chemistry, School of Chemistry and Chemical Engineering, Chongqing University, Daxuecheng South Rd. 55, 401331 Chongqing, China.
| | - Siyu Zou
- Department of Chemistry, School of Chemistry and Chemical Engineering, Chongqing University, Daxuecheng South Rd. 55, 401331 Chongqing, China.
| | - Chenxu Wang
- Department of Chemistry, School of Chemistry and Chemical Engineering, Chongqing University, Daxuecheng South Rd. 55, 401331 Chongqing, China.
| | - Xuefang Xu
- Department of Chemistry, School of Chemistry and Chemical Engineering, Chongqing University, Daxuecheng South Rd. 55, 401331 Chongqing, China.
- Chongqing Key Laboratory of Theoretical and Computational Chemistry, 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.
- Chongqing Key Laboratory of Theoretical and Computational Chemistry, Chongqing University, Daxuecheng South Rd. 55, 401331, Chongqing, China
| |
Collapse
|
15
|
Puzzarini C, Alessandrini S, Bizzocchi L, Melosso M. Hunting for interstellar molecules: rotational spectra of reactive species. Faraday Discuss 2023; 245:309-326. [PMID: 37318358 DOI: 10.1039/d3fd00052d] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/21/2023]
Abstract
Interstellar molecules are often highly reactive species, which are unstable under terrestrial conditions, such as radicals, ions and unsaturated carbon chains. Their detection in space is usually based on the astronomical observation of their rotational fingerprints. However, laboratory investigations have to face the issue of efficiently producing these molecules and preserving them during rotational spectroscopy measurements. A general approach for producing and investigating unstable/reactive species is presented by means of selected case-study molecules. The overall strategy starts from quantum-chemical calculations that aim at obtaining accurate predictions of the missing spectroscopic information required to guide spectral analysis and assignment. Rotational spectra of these species are then recorded by exploiting the approach mentioned above, and their subsequent analysis leads to accurate spectroscopic parameters. These are then used for setting up accurate line catalogs for astronomical searches.
Collapse
Affiliation(s)
- Cristina Puzzarini
- ROT&Comp Lab, Department of Chemistry "Giacomo Ciamician", University of Bologna, Via F. Selmi 2, I-40126 Bologna, Italy.
| | - Silvia Alessandrini
- ROT&Comp Lab, Department of Chemistry "Giacomo Ciamician", University of Bologna, Via F. Selmi 2, I-40126 Bologna, Italy.
| | - Luca Bizzocchi
- ROT&Comp Lab, Department of Chemistry "Giacomo Ciamician", University of Bologna, Via F. Selmi 2, I-40126 Bologna, Italy.
| | - Mattia Melosso
- ROT&Comp Lab, Department of Chemistry "Giacomo Ciamician", University of Bologna, Via F. Selmi 2, I-40126 Bologna, Italy.
| |
Collapse
|
16
|
Thimmakondu VS, Karton A. CCSD(T) Rotational Constants for Highly Challenging C 5H 2 Isomers-A Comparison between Theory and Experiment. Molecules 2023; 28:6537. [PMID: 37764314 PMCID: PMC10537648 DOI: 10.3390/molecules28186537] [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/30/2023] [Revised: 09/06/2023] [Accepted: 09/07/2023] [Indexed: 09/29/2023] Open
Abstract
We evaluate the accuracy of CCSD(T) and density functional theory (DFT) methods for the calculation of equilibrium rotational constants (Ae, Be, and Ce) for four experimentally detected low-lying C5H2 isomers (ethynylcyclopropenylidene (2), pentatetraenylidene (3), ethynylpropadienylidene (5), and 2-cyclopropen-1-ylidenethenylidene (8)). The calculated rotational constants are compared to semi-experimental rotational constants obtained by converting the vibrationally averaged experimental rotational constants (A0, B0, and C0) to equilibrium values by subtracting the vibrational contributions (calculated at the B3LYP/jun-cc-pVTZ level of the theory). The considered isomers are closed-shell carbenes, with cumulene, acetylene, or strained cyclopropene moieties, and are therefore highly challenging from an electronic structure point of view. We consider both frozen-core and all-electron CCSD(T) calculations, as well as a range of DFT methods. We find that calculating the equilibrium rotational constants of these C5H2 isomers is a difficult task, even at the CCSD(T) level. For example, at the all-electron CCSD(T)/cc-pwCVTZ level of the theory, we obtain percentage errors ≤0.4% (Ce of isomer 3, Be and Ce of isomer 5, and Be of isomer 8) and 0.9-1.5% (Be and Ce of isomer 2, Ae of isomer 5, and Ce of isomer 8), whereas for the Ae rotational constant of isomers 2 and 8 and Be rotational constant of isomer 3, high percentage errors above 3% are obtained. These results highlight the challenges associated with calculating accurate rotational constants for isomers with highly challenging electronic structures, which is further complicated by the need to convert vibrationally averaged experimental rotational constants to equilibrium values. We use our best CCSD(T) rotational constants (namely, ae-CCSD(T)/cc-pwCVTZ for isomers 2 and 5, and ae-CCSD(T)/cc-pCVQZ for isomers 3 and 8) to evaluate the performance of DFT methods across the rungs of Jacob's Ladder. We find that the considered pure functionals (BLYP-D3BJ, PBE-D3BJ, and TPSS-D3BJ) perform significantly better than the global and range-separated hybrid functionals. The double-hybrid DSD-PBEP86-D3BJ method shows the best overall performance, with percentage errors below 0.5% in nearly all cases.
Collapse
Affiliation(s)
- Venkatesan S Thimmakondu
- Department of Chemistry and Biochemistry, San Diego State University, San Diego, CA 92182-1030, USA
| | - Amir Karton
- School of Science and Technology, University of New England, Armidale, NSW 2351, Australia
| |
Collapse
|
17
|
Barone V. Accurate structures and spectroscopic parameters of α,α-dialkylated α-amino acids in the gas-phase: a joint venture of DFT and wave-function composite methods. Phys Chem Chem Phys 2023; 25:22768-22774. [PMID: 37591810 DOI: 10.1039/d3cp02503a] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/19/2023]
Abstract
Accurate computations of structural, conformational and spectroscopic properties in the gas phase have been performed for two α,α-dialkylated α-amino acids, namely aminoisobutyric acid and cyclopropylglycine. Thanks to the integration of modern double hybrid functionals and wave-function methods, several low-energy structures of the title molecules could be analyzed employing standard computer resources. The computed features of all the most stable conformers of the target amino acids closely match the corresponding spectroscopic parameters issued from microwave spectroscopic studies in the gas-phase. Together with their intrinsic interest, the accuracy of the results obtained with reasonable computer times paves the way for accurate investigations of other flexible bricks of life.
Collapse
Affiliation(s)
- Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy.
| |
Collapse
|
18
|
Barone V, Di Grande S, Lazzari F, Mendolicchio M. Accurate Structures and Spectroscopic Parameters of Guanine Tautomers in the Gas Phase by the Pisa Conventional and Explicitly Correlated Composite Schemes (PCS and PCS-F12). J Phys Chem A 2023; 127:6771-6778. [PMID: 37535450 PMCID: PMC10440789 DOI: 10.1021/acs.jpca.3c03999] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2023] [Revised: 07/21/2023] [Indexed: 08/05/2023]
Abstract
A general strategy for the accurate computation of structural and spectroscopic properties of biomolecule building blocks in the gas phase is proposed and validated for tautomeric equilibria. The main features of the new model are the inclusion of core-valence correlation in geometry optimizations by a double hybrid functional and the systematic use of wave-function composite methods in conjunction with cc-pVnZ-F12 basis sets with separate extrapolation of MP2 and post-MP2 contributions. The resulting Pisa composite scheme employing conventional (PCS) or explicitly correlated (PCS-F12) approaches is applied to the challenging problem of guanine tautomers in the gas phase. The results are in remarkable agreement with the experimental structures, relative stabilities, and spectroscopic signatures of different tautomers. The accuracy of the results obtained at reasonable cost by means of black-box parameter-free approaches paves the way toward systematic investigations of other molecular bricks of life also by non-specialists.
Collapse
Affiliation(s)
- Vincenzo Barone
- Scuola
Normale Superiore, Piazza dei Cavalieri 7, Pisa 56126, Italy
| | - Silvia Di Grande
- Scuola
Normale Superiore, Piazza dei Cavalieri 7, Pisa 56126, Italy
- Scuola
Superiore Meridionale, Largo San Marcellino 10, Napoli 80138, Italy
| | - Federico Lazzari
- Scuola
Normale Superiore, Piazza dei Cavalieri 7, Pisa 56126, Italy
| | | |
Collapse
|
19
|
Barone V. DFT Meets Wave-Function Composite Methods for Characterizing Cytosine Tautomers in the Gas Phase. J Chem Theory Comput 2023; 19:4970-4981. [PMID: 37479680 PMCID: PMC10413851 DOI: 10.1021/acs.jctc.3c00465] [Citation(s) in RCA: 10] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2023] [Indexed: 07/23/2023]
Abstract
A general strategy for the accurate computation of structural and spectroscopic properties of biomolecule building blocks in the gas phase has been further improved and validated with a special reference to tautomeric equilibria. The main improvements concern the use of the cc-pVTZ-F12 basis set in both DFT and CCSD(T)-F12 computations, the inclusion of core-valence correlation in geometry optimizations by double hybrid functionals, and the use of the cc-pVQZ-F12 basis set for complete basis set extrapolation at the MP2-F12 level. The resulting model chemistry is applied to the challenging problem of cytosine tautomers in the gas phase. The results are in remarkable agreement with experiment concerning both rotational and vibrational spectroscopic parameters and permit their unbiased interpretation in terms of structural and thermochemical features. Together with the intrinsic interest of the studied molecule, the accuracy of the results obtained at reasonable cost without any empirical parameter suggests that the proposed composite method can be profitably employed for accurate investigations of other molecular bricks of life.
Collapse
Affiliation(s)
- Vincenzo Barone
- Scuola Normale Superiore
di Pisa, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
20
|
Abstract
pbqff is an open-source program for fully automating the production of quartic force fields (QFFs) and their corresponding anharmonic spectroscopic data. Rather than being a monolithic piece of code, it consists of several key modules including a generic interface to quantum chemistry codes and, notably, queuing systems; a molecular point group symmetry library; an internal-to-Cartesian coordinate conversion module; a module for the ordinary least-squares fitting of potential energy surfaces; and an improved second-order rotational and vibrational perturbation theory package for asymmetric and symmetric tops that handles type-1 and -2 Fermi resonances, Fermi resonance polyads, and Coriolis resonances. All of these pieces are written in Rust, a modern, safe, and performant programming language that has much to offer for scientific programming. This work introduces pbqff and its surrounding ecosystem, in addition to reporting new anharmonic vibrational data for c-(C)C3H2 and describing how the components of pbqff can be leveraged in other projects.
Collapse
Affiliation(s)
- Brent R Westbrook
- Department of Chemistry & Biochemistry, University of Mississippi, University, Mississippi 38677-1848, United States
| | - Ryan C Fortenberry
- Department of Chemistry & Biochemistry, University of Mississippi, University, Mississippi 38677-1848, United States
| |
Collapse
|
21
|
Hadki HE, Koziol KJ, Kabbaj OK, Komiha N, Kleiner I, Nguyen HVL. The Microwave Rotational Electric Resonance (RER) Spectrum of Benzothiazole. Molecules 2023; 28:molecules28083419. [PMID: 37110653 PMCID: PMC10146593 DOI: 10.3390/molecules28083419] [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: 03/14/2023] [Revised: 04/04/2023] [Accepted: 04/06/2023] [Indexed: 04/29/2023] Open
Abstract
The microwave spectra of benzothiazole were measured in the frequency range 2-26.5 GHz using a pulsed molecular jet Fourier transform microwave spectrometer. Hyperfine splittings arising from the quadrupole coupling of the 14N nucleus were fully resolved and analyzed simultaneously with the rotational frequencies. In total, 194 and 92 hyperfine components of the main species and the 34S isotopologue, respectively, were measured and fitted to measurement accuracy using a semi-rigid rotor model supplemented by a Hamiltonian accounting for the 14N nuclear quadrupole coupling effect. Highly accurate rotational constants, centrifugal distortion constants, and 14N nuclear quadrupole coupling constants were deduced. A large number of method and basis set combinations were used to optimize the molecular geometry of benzothiazole, and the calculated rotational constants were compared with the experimentally determined constants in the course of a benchmarking effort. The similar value of the χcc quadrupole coupling constant when compared to other thiazole derivatives indicates only very small changes of the electronic environment at the nitrogen nucleus in these compounds. The small negative inertial defect of -0.056 uÅ2 hints that low-frequency out-of-plane vibrations are present in benzothiazole, similar to the observation for some other planar aromatic molecules.
Collapse
Affiliation(s)
- Hamza El Hadki
- Laboratory of Spectroscopy, Molecular Modeling, Materials, Nanomaterials, Water, and Environment, Faculty of Sciences, Mohammed V University, Av Ibn Battouta, Rabat B.P. 1014, Morocco
- Univ Paris Est Creteil and Université Paris Cité, CNRS, LISA, 94010 Créteil, France
| | - Kenneth J Koziol
- Univ Paris Est Creteil and Université Paris Cité, CNRS, LISA, 94010 Créteil, France
| | - Oum Keltoum Kabbaj
- Laboratory of Spectroscopy, Molecular Modeling, Materials, Nanomaterials, Water, and Environment, Faculty of Sciences, Mohammed V University, Av Ibn Battouta, Rabat B.P. 1014, Morocco
| | - Najia Komiha
- Laboratory of Spectroscopy, Molecular Modeling, Materials, Nanomaterials, Water, and Environment, Faculty of Sciences, Mohammed V University, Av Ibn Battouta, Rabat B.P. 1014, Morocco
| | - Isabelle Kleiner
- Université Paris Cité and Univ Paris Est Creteil, CNRS, LISA, 75013 Paris, France
| | - Ha Vinh Lam Nguyen
- Univ Paris Est Creteil and Université Paris Cité, CNRS, LISA, 94010 Créteil, France
- Institut Universitaire de France (IUF), 75231 Paris, France
| |
Collapse
|
22
|
Barone V, Fusè M. Accurate Structures and Spectroscopic Parameters of Phenylalanine and Tyrosine in the Gas Phase: A Joint Venture of DFT and Composite Wave-Function Methods. J Phys Chem A 2023; 127:3648-3657. [PMID: 37052318 PMCID: PMC10150396 DOI: 10.1021/acs.jpca.3c01174] [Citation(s) in RCA: 11] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/14/2023]
Abstract
A general strategy for the accurate computation of conformational and spectroscopic properties of flexible molecules in the gas phase is applied to two representative proteinogenic amino acids with aromatic side chains, namely, phenylalanine and tyrosine. The main features of all the most stable conformers predicted by this computational strategy closely match those of the species detected in microwave and infrared experiments. Together with their intrinsic interest, the accuracy of the results obtained with reasonable computer times paves the route for accurate investigations of other flexible bricks of life.
Collapse
Affiliation(s)
- Vincenzo Barone
- Scuola Normale Superiore di Pisa, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Marco Fusè
- DMMT-sede Europa, Università di Brescia, Viale Europa 11, 25121 Brescia, Italy
| |
Collapse
|
23
|
Alessandrini S, Melosso M, Rivilla VM, Bizzocchi L, Puzzarini C. Computational Protocol for the Identification of Candidates for Radioastronomical Detection and Its Application to the C3H3NO Family of Isomers. Molecules 2023; 28:molecules28073226. [PMID: 37049990 PMCID: PMC10096335 DOI: 10.3390/molecules28073226] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2023] [Revised: 03/30/2023] [Accepted: 03/31/2023] [Indexed: 04/08/2023] Open
Abstract
The C3H3NO family of isomers is relevant in astrochemistry, even though its members are still elusive in the interstellar medium. To identify the best candidate for astronomical detection within this family, we developed a new computational protocol based on the minimum-energy principle. This approach aims to identify the most stable isomer of the family and consists of three steps. The first step is an extensive investigation that characterizes the vast number of compounds having the C3H3NO chemical formula, employing density functional theory for this purpose. The second step is an energy refinement, which is used to select isomers and relies on coupled cluster theory. The last step is a structural improvement with a final energy refinement that provides improved energies and a large set of accurate spectroscopic parameters for all isomers lying within 30 kJ mol−1 above the most stable one. According to this protocol, vinylisocyanate is the most stable isomer, followed by oxazole, which is about 5 kJ mol−1 higher in energy. The other stable species are pyruvonitrile, cyanoacetaldehyde, and cyanovinylalcohol. For all of these species, new computed rotational and vibrational spectroscopic data are reported, which complement those already available in the literature or fill current gaps.
Collapse
Affiliation(s)
- Silvia Alessandrini
- Dipartimento di Chimica “Giacomo Ciamician”, Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy
| | - Mattia Melosso
- Dipartimento di Chimica “Giacomo Ciamician”, Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy
| | - Víctor M. Rivilla
- Centro de Astrobiología (CSIC-INTA), Ctra. de Ajalvir km. 4, Torrejón de Ardoz, 28850 Madrid, Spain
| | - Luca Bizzocchi
- Dipartimento di Chimica “Giacomo Ciamician”, Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy
| | - Cristina Puzzarini
- Dipartimento di Chimica “Giacomo Ciamician”, Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy
| |
Collapse
|
24
|
Sutikdja LW, Nguyen HVL, Jelisavac D, Stahl W, Mouhib H. Benchmarking quantum chemical methods for accurate gas-phase structure predictions of carbonyl compounds: the case of ethyl butyrate. Phys Chem Chem Phys 2023; 25:7688-7696. [PMID: 36857713 PMCID: PMC10015624 DOI: 10.1039/d2cp05774c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/11/2023]
Abstract
High-resolution spectroscopy techniques play a pivotal role to validate and efficiently benchmark available methods from quantum chemistry. In this work, we analyzed the microwave spectrum of ethyl butyrate within the scope of a systematic investigation to benchmark state-of-the-art exchange-correlation functionals and ab initio methods, to accurately predict the lowest energy conformers of carbonyl compounds in their isolated state. Under experimental conditions, we observed two distinct conformers, one of Cs and one of C1 symmetry. As reported earlier in the cases of some ethyl and methyl alkynoates, structural optimizations of the most abundant conformer that exhibits a C1 symmetry proved extremely challenging for several quantum chemical levels. To probe the sensitivity of different methods and basis sets, we use the identified soft-degree of freedom in proximity to the carbonyl group as an order parameter. The results of our study provide useful insight for spectroscopists to select an adapted method for structure prediction of carbonyl compounds based on their available computational resources, suggesting a reasonable trade-off between accuracy and CPU cost. At the same time, our observations and the resulting sets of highly accurate experimental constants from high-resolution spectroscopy experiments give an appeal to theoretical groups to look further into this seemingly simple family of chemical compounds, which may prove useful for the further development and parametrization of theoretical methods in computational chemistry.
Collapse
Affiliation(s)
- Lilian W Sutikdja
- Institute of Physical Chemistry, RWTH Aachen University, Landoltweg 2, D-52074, Aachen, Germany
| | - Ha Vinh Lam Nguyen
- Univ Paris Est Creteil and Université Paris Cité, CNRS, LISA, F-94010, Créteil, France. .,Institut Universitaire de France (IUF), F-75231, Paris cedex 05, France
| | - Dragan Jelisavac
- Institute of Physical Chemistry, RWTH Aachen University, Landoltweg 2, D-52074, Aachen, Germany
| | - Wolfgang Stahl
- Institute of Physical Chemistry, RWTH Aachen University, Landoltweg 2, D-52074, Aachen, Germany
| | - Halima Mouhib
- Department of Computer Science, VU Bioinformatics, Vrije Universiteit Amsterdam, De Boelelaan 1111, 1081 HV, Amsterdam, The Netherlands.
| |
Collapse
|
25
|
Barone V, Fusè M, Lazzari F, Mancini G. Benchmark Structures and Conformational Landscapes of Amino Acids in the Gas Phase: A Joint Venture of Machine Learning, Quantum Chemistry, and Rotational Spectroscopy. J Chem Theory Comput 2023; 19:1243-1260. [PMID: 36731119 PMCID: PMC9979611 DOI: 10.1021/acs.jctc.2c01143] [Citation(s) in RCA: 16] [Impact Index Per Article: 16.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
The accurate characterization of prototypical bricks of life can strongly benefit from the integration of high resolution spectroscopy and quantum mechanical computations. We have selected a number of representative amino acids (glycine, alanine, serine, cysteine, threonine, aspartic acid and asparagine) to validate a new computational setup rooted in quantum-chemical computations of increasing accuracy guided by machine learning tools. Together with low-lying energy minima, the barriers ruling their interconversion are evaluated in order to unravel possible fast relaxation paths. Vibrational and thermal effects are also included in order to estimate relative free energies at the temperature of interest in the experiment. The spectroscopic parameters of all the most stable conformers predicted by this computational strategy, which do not have low-energy relaxation paths available, closely match those of the species detected in microwave experiments. Together with their intrinsic interest, these accurate results represent ideal benchmarks for more approximate methods.
Collapse
Affiliation(s)
- Vincenzo Barone
- Scuola
Normale Superiore di Pisa, piazza dei Cavalieri 7, 56126 Pisa, Italy,
| | - Marco Fusè
- DMMT-sede
Europa, Universitá di Brescia, viale Europa 11, 25121 Brescia, Italy
| | - Federico Lazzari
- Scuola
Normale Superiore di Pisa, piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Giordano Mancini
- Scuola
Normale Superiore di Pisa, piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
26
|
Barone V, Di Grande S, Puzzarini C. Toward Accurate yet Effective Computations of Rotational Spectroscopy Parameters for Biomolecule Building Blocks. MOLECULES (BASEL, SWITZERLAND) 2023; 28:molecules28020913. [PMID: 36677970 PMCID: PMC9863398 DOI: 10.3390/molecules28020913] [Citation(s) in RCA: 9] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 01/01/2023] [Revised: 01/01/2023] [Accepted: 01/12/2023] [Indexed: 01/19/2023]
Abstract
The interplay of high-resolution rotational spectroscopy and quantum-chemical computations plays an invaluable role in the investigation of biomolecule building blocks in the gas phase. However, quantum-chemical methods suffer from unfavorable scaling with the dimension of the system under consideration. While a complete characterization of flexible systems requires an elaborate multi-step strategy, in this work, we demonstrate that the accuracy obtained by quantum-chemical composite approaches in the prediction of rotational spectroscopy parameters can be approached by a model based on density functional theory. Glycine and serine are employed to demonstrate that, despite its limited cost, such a model is able to predict rotational constants with an accuracy of 0.3% or better, thus paving the way toward the accurate characterization of larger flexible building blocks of biomolecules.
Collapse
Affiliation(s)
- Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-50126 Pisa, Italy
- Correspondence: (V.B.); (C.P.)
| | - Silvia Di Grande
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-50126 Pisa, Italy
- Scuola Superiore Meridionale, Largo San Marcellino 10, I-80138 Napoli, Italy
| | - Cristina Puzzarini
- Rotational and Computational Spectroscopy Lab, Department of Chemistry “Giacomo Ciamician”, University of Bologna, Via F. Selmi 2, I-40126 Bologna, Italy
- Correspondence: (V.B.); (C.P.)
| |
Collapse
|