1
|
Jara-Toro RA, Roncero O, Lique F. The temperature variation of the CH + + H reaction rate coefficients: a puzzle finally understood? Phys Chem Chem Phys 2024; 26:21370-21378. [PMID: 39099553 PMCID: PMC11323935 DOI: 10.1039/d4cp01902d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2024] [Accepted: 07/19/2024] [Indexed: 08/06/2024]
Abstract
CH+ was the first molecular ion identified in the interstellar medium and is found to be ubiquitous in interstellar clouds. However, its formation and destruction paths are not well understood, especially at low temperatures. A new theoretical approach based on the canonical variational transition state theory was used to study the H + CH+ reactive collisions. Rate coefficients for formation of C+ ions are calculated as a function of temperature. We considered the participation of a direct path and an indirect path in which the reactants should overcome an entropic barrier to form a van der Waals complex or pass through a CH2+ intermediate complex, respectively. We show that the contribution of both pathways to the formation of C+ has to be taken into account. The new reactive rate coefficients for the title reaction, complemented by reactive data for CH+/CH2+ in the H/H2/He mixture, have been used to simulate the corresponding kinetics experimentally measured using an Atomic Beam 22 Pole Trap apparatus at low temperature. A good agreement with the experimental findings was found at 50 K. At a lower temperature, the model overestimates the formation of C+. This shows that secondary reactions are not responsible for the weak C+ production in the experiments at such temperature. Then, we discuss the possible impact of non-adiabatic effects in the study of the H + CH+ reactive collisions and we found that such effects can be responsible for the decrease of the H + CH+ rate coefficients at low temperature. This study offers an explanation for the disagreement between H + CH+ theoretical and experimental rate coefficients which has been going on for 20 years and highlights the need for performing non-adiabatic studies for this simple chemical reaction.
Collapse
Affiliation(s)
- Rafael A Jara-Toro
- Univ Rennes, CNRS, IPR (Institut de Physique de Rennes), UMR 6251, F-35000 Rennes, France.
| | - Octavio Roncero
- Instituto de Física Fundamental (IFF-CSIC), C.S.I.C., Serrano 123, 28006 Madrid, Spain
| | - François Lique
- Univ Rennes, CNRS, IPR (Institut de Physique de Rennes), UMR 6251, F-35000 Rennes, France.
| |
Collapse
|
2
|
Minnette W, Hoy EP, Sand AM. The Use of Effective Core Potentials with Multiconfiguration Pair-Density Functional Theory. J Phys Chem A 2024; 128:6555-6565. [PMID: 39052857 PMCID: PMC11317981 DOI: 10.1021/acs.jpca.4c02666] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/24/2024] [Revised: 06/25/2024] [Accepted: 07/09/2024] [Indexed: 07/27/2024]
Abstract
The reliable and accurate prediction of chemical properties is a key goal in quantum chemistry. Transition-metal-containing complexes can often pose difficulties to quantum mechanical methods for multiple reasons, including many electron configurations contributing to the overall electronic description of the system and the large number of electrons significantly increasing the amount of computational resources required. Often, multiconfigurational electronic structure methods are employed for such systems, and the cost of these calculations can be reduced by the use of an effective core potential (ECP). In this work, we explore both theoretical considerations and performances of ECPs applied in the context of multiconfiguration pair-density functional theory (MC-PDFT). A mixed-basis set approach is used, using ECP basis sets for transition metals and all-electron basis sets for nonmetal atoms. We illustrate the effects that an ECP has on the key parameters used in the computation of MC-PDFT energies, and we explore how ECPs affect the prediction of physical observables for chemical systems. The dissociation curve for a metal dimer was explored, and ionization energies for transition metal-containing diatomic systems were computed and compared to experimental values. In general, we find that ECP approaches employed with MC-PDFT are able to predict ionization energies with improved accuracy compared to traditional Kohn-Sham density functional theory approaches.
Collapse
Affiliation(s)
- William
E. Minnette
- Department
of Chemistry and Biochemistry, Butler University, Indianapolis, Indiana 46208, United States
| | - Erik P. Hoy
- Department
of Chemistry, Rowan University, Glassboro, New Jersey 08028, United States
| | - Andrew M. Sand
- Department
of Chemistry and Biochemistry, Butler University, Indianapolis, Indiana 46208, United States
| |
Collapse
|
3
|
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
|
4
|
García de la Concepción J, Corchado JC, Cintas P, Babiano R. Norcaradiene-Cycloheptatriene Equilibrium: A Heavy-Atom Quantum Tunneling Case. J Org Chem 2024; 89:9336-9343. [PMID: 38888485 PMCID: PMC11232008 DOI: 10.1021/acs.joc.4c00464] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2024] [Revised: 05/31/2024] [Accepted: 06/07/2024] [Indexed: 06/20/2024]
Abstract
The equilibrium between norcaradiene and cycloheptatriene, which has captivated chemists for more than half a century, is revisited by state-of-the-art quantum chemical calculations. Our theoretical data significantly deviate from the experimental results (J. Am. Chem. Soc., 1981, 26, 7791-7792), especially at low temperatures, where isomerization is dominated by heavy-atom tunneling. This effect results in an extremely short half-life for norcaradiene, rendering it undetectable. This work sheds light on this equilibrium, updating the kinetic and thermodynamic data while also expanding the repertoire of organic reactions controlled by this exotic quantum effect.
Collapse
Affiliation(s)
- Juan García de la Concepción
- Departamento de Química Orgánica e Inorgánica, Facultad de Ciencias, and IACYS-Green Chemistry and Sustainable Development Unit, Universidad de Extremadura, 06006 Badajoz, Spain
| | - José C Corchado
- Departamento de Ingeniería Química y Química Física, Facultad de Ciencias, and ICCAEx, Universidad Extremadura, 06006 Badajoz, Spain
| | - Pedro Cintas
- Departamento de Química Orgánica e Inorgánica, Facultad de Ciencias, and IACYS-Green Chemistry and Sustainable Development Unit, Universidad de Extremadura, 06006 Badajoz, Spain
| | - Reyes Babiano
- Departamento de Química Orgánica e Inorgánica, Facultad de Ciencias, and IACYS-Green Chemistry and Sustainable Development Unit, Universidad de Extremadura, 06006 Badajoz, Spain
| |
Collapse
|
5
|
Mehta N, Martin JML. On the sensitivity of computed partial charges toward basis set and (exchange-)correlation treatment. J Comput Chem 2024; 45:1017-1032. [PMID: 38216516 DOI: 10.1002/jcc.27294] [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: 09/25/2023] [Revised: 12/05/2023] [Accepted: 12/10/2023] [Indexed: 01/14/2024]
Abstract
Partial charges are a central concept in general chemistry and chemical biology, yet dozens of different computational definitions exist. In prior work [Cho et al., ChemPhysChem 21, 688-696 (2020)], we showed that these can be reduced to at most three 'principal components of ionicity'. The present study addressed the dependence of computed partial charges q on 1-particle basis set and (for WFT methods) n -particle correlation treatment or (for DFT methods) exchange-correlation functional, for several representative partial charge definitions such as QTAIM, Hirshfeld, Hirshfeld-I, HLY (electrostatic), NPA, and GAPT. Our findings show that semi-empirical double hybrids can closely approach the CCSD(T) 'gold standard' for this property. In fact, owing to an error compensation in MP2, CCSD partial charges are further away from CCSD(T) than is MP2. The nonlocal correlation is important, especially when there is a substantial amount of nonlocal exchange. Employing range separation proves to be "mostly" not advantageous, while global hybrids perform optimally for 20%-30% Hartree-Fock exchange across all charge types. Basis set convergence analysis shows that an augmented triple-zeta heavy-aug-cc-pV(T+d)Z basis set or a partially augmented jun-cc-pV(T+d)Z basis set is sufficient for Hirshfeld, Hirshfeld-I, HLY, and GAPT charges. In contrast, QTAIM and NPA display slower basis set convergence. It is noteworthy that for both NPA and QTAIM, HF exhibits markedly slower basis set convergence than the correlation components of MP2 and CCSD. Triples corrections in CCSD(T), denoted as CCSD(T)-CCSD, exhibit even faster basis set convergence.
Collapse
Affiliation(s)
- Nisha Mehta
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Reḥovot, Israel
| | - Jan M L Martin
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Reḥovot, Israel
| |
Collapse
|
6
|
Gao M, He J, Tian L, Chen L, Shi S, Zhang C, Chen D. Comparative Study on the H-Abstraction Reactions of Isopropyl Acetate and Propyl Acetate by HO 2 and OH Radicals. J Phys Chem A 2024; 128:2825-2842. [PMID: 38563141 DOI: 10.1021/acs.jpca.4c00794] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/04/2024]
Abstract
Isopropyl acetate (IPA) and propyl acetate (PA) are recognized as promising biofuels suitable for applications as fuel additives and biodiesel models. The H-abstraction reactions with radicals stand out as the fundamental initiating reactions in the combustion kinetic models for IPA and PA. In the present work, the kinetic calculations of IPA and PA plus HO2 and OH radicals were investigated at M06-2X/cc-pVTZ//G4, M08-HX/maug-cc-pVTZ, and CCSD(T)/jul-cc-pVTZ levels. The thermodynamic calculations were obtained based on the G4 and CBS-APNO methods. Rate coefficients were calculated using both transition state theory and canonical variational transition state theory with tunneling correction at the temperature range of 250-2000 K. The total rate constants for the IPA + OH system were fitted as follows: k = 0.4674 × T3.927 exp(2128/T) (cm3 mol-1 s-1), and for the PA + OH system, the total rate constants were determined using the following equation: k = 0.0161 × T4.373 exp(2220/T) (cm3 mol-1 s-1). The rate coefficients of IPA + OH reactions determined based on the M08-HX/maug-cc-pVTZ level effectively replicate the experimental data, while H-abstraction rate coefficients of PA + OH by the CCSD(T)/jul-cc-pVTZ method accurately reproduce the experimental data. Refining the H-abstraction rate coefficients in the kinetic mechanism of PA, as proposed by Dayma et al. [Proc. Combust. Inst. 37 (2019) 429-436], has been achieved through incorporating the present calculated data, leading to the development of a revised mechanism. The validation of the updated mechanism against jet-stirred reactor data is presented, showcasing its effective performance in predicting JSR data.
Collapse
Affiliation(s)
- Mengjiao Gao
- College of Mathematics and Physics, Chengdu University of Technology, Chengdu, Sichuan 610059, China
| | - Jiuning He
- College of Mathematics and Physics, Chengdu University of Technology, Chengdu, Sichuan 610059, China
| | - Lilan Tian
- College of Mathematics and Physics, Chengdu University of Technology, Chengdu, Sichuan 610059, China
| | - Lei Chen
- College of Mathematics and Physics, Chengdu University of Technology, Chengdu, Sichuan 610059, China
| | - Shunping Shi
- College of Mathematics and Physics, Chengdu University of Technology, Chengdu, Sichuan 610059, China
| | - Changhua Zhang
- Institute of Atomic and Molecular Physics, Sichuan University, Chengdu 610065, China
| | - Deliang Chen
- School of Physics and Electronic, Guizhou Education University, Guiyang, Guizhou 550018, China
| |
Collapse
|
7
|
Cardia R, Dardenne N, Mula G, Pinna E, Rignanese GM, Charlier JC, Cappellini G. First-Principles Investigation of the Optical Properties of Eumelanin Protomolecules. J Phys Chem A 2023; 127:10797-10806. [PMID: 38109190 DOI: 10.1021/acs.jpca.3c04898] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2023]
Abstract
Using first-principles calculations, we investigate the absorption spectra (in the near-infrared, visible, and first UV range) of the two most probable eumelanin tetrameric molecules exhibiting either a linear open-chain or a cyclic porphyrine-like configuration. In order to simulate a realistic molecular system, an implicit solvent model is used in our calculations to mimic the effect of the solvated environment around the eumelanin molecule. Although the presence of solvent is found not to significantly affect the absorption pattern of both molecules, the onset of the spectra are shifted toward higher energies, especially for the linear tetramer. Interestingly, the absorption spectra and optical onsets of the two molecules differ significantly both in a vacuum and in ethanol. However, the two predicted spectra do not allow us to definitely discriminate between the two configurations when comparing the theoretical predictions with the available experimental spectrum. In addition, a mix of the two eumelanin configurations (close to fifty-fifty) leads to a maximum overlap between theoretical and experimental spectra. Consequently, this theoretical research shows that deeper insight can be gained using beyond DFT techniques on the real form of eumelanin protomolecules present in living systems as well as on their possible use in hybrid solar cells.
Collapse
Affiliation(s)
- Roberto Cardia
- Department of Physics, Università degli Studi di Cagliari, Cittadella Universitaria I-09042 Monserrato, Cagliari, Italy
| | - Nicolas Dardenne
- Institute of Condensed Matter and Nanosciences, Université catholique de Louvain (UCLouvain), B-1348 Louvain-la-Neuve, Belgium
| | - Guido Mula
- Department of Physics, Università degli Studi di Cagliari, Cittadella Universitaria I-09042 Monserrato, Cagliari, Italy
| | - Elisa Pinna
- Department of Physics, Università degli Studi di Cagliari, Cittadella Universitaria I-09042 Monserrato, Cagliari, Italy
| | - Gian-Marco Rignanese
- Institute of Condensed Matter and Nanosciences, Université catholique de Louvain (UCLouvain), B-1348 Louvain-la-Neuve, Belgium
- European Theoretical Spectroscopy Facility (ETSF)
| | - Jean-Christophe Charlier
- Institute of Condensed Matter and Nanosciences, Université catholique de Louvain (UCLouvain), B-1348 Louvain-la-Neuve, Belgium
- European Theoretical Spectroscopy Facility (ETSF)
| | - Giancarlo Cappellini
- Department of Physics, Università degli Studi di Cagliari, Cittadella Universitaria I-09042 Monserrato, Cagliari, Italy
- European Theoretical Spectroscopy Facility (ETSF)
| |
Collapse
|
8
|
Chebbi W, Derbel N, Alijah A, Cours T. UV-spectrum and photodecomposition of peroxynitrous acid in the troposphere. Phys Chem Chem Phys 2023; 26:123-129. [PMID: 38059643 DOI: 10.1039/d3cp04580c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/08/2023]
Abstract
The UV spectrum of peroxynitrous acid, HOONO, was computed at the B3LYP/AVTZ and MCSCF/AVTZ levels using the fewest switches surface hopping algorithm. Due to large-amplitude vibrational motions of this molecule, the maxima in the simulated spectra are displaced from the positions of vertical excitations. The three lowest excited electronic singlet states, which are all repulsive, can be reached by UV absorption. The photolysis products are determined, and the photolysis rate constant is provided for the first time. We found that near the tropopause the photolysis rate constant J ≈ 6 × 10-4 s-1, exceeds that for thermal decomposition by two orders of magnitude. The photolysis lifetime is about 30 minutes. Thus, photolysis is an important process and should be included in atmospheric models.
Collapse
Affiliation(s)
- Wiem Chebbi
- LSAMA, Laboratoire de Spectroscopie Atomique, Moléculaire et Applications, Department of Physics, University Tunis - El Manar, 1060 Tunis, Tunisia
- GSMA, Groupe de Spectrométrie Moléculaire et Atmosphérique, UMR CNRS 7331, University of Reims Champagne-Ardenne, 51100 Reims, France.
| | - Najoua Derbel
- LSAMA, Laboratoire de Spectroscopie Atomique, Moléculaire et Applications, Department of Physics, University Tunis - El Manar, 1060 Tunis, Tunisia
| | - Alexander Alijah
- GSMA, Groupe de Spectrométrie Moléculaire et Atmosphérique, UMR CNRS 7331, University of Reims Champagne-Ardenne, 51100 Reims, France.
| | - Thibaud Cours
- GSMA, Groupe de Spectrométrie Moléculaire et Atmosphérique, UMR CNRS 7331, University of Reims Champagne-Ardenne, 51100 Reims, France.
| |
Collapse
|
9
|
K Chandy S, Bowers S, Raghavachari K, Li LS. Structure Dependence of CO 2 Reduction Electrocatalyzed by Metal-Nanographene Complexes: A Computational Study. J Phys Chem A 2023; 127:8566-8573. [PMID: 37796447 DOI: 10.1021/acs.jpca.3c04564] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/06/2023]
Abstract
Improving the energy efficiency of electrocatalytic reduction of CO2 requires tuning of redox properties of electrocatalysts to match redox potentials of the substrate. Recently, we introduced nanographenes as ligands for metal complexes for such purposes by taking advantage of size-dependent properties of the conjugated systems. Here, we use computations to investigate the structure dependence of the electrocatalysis at Re(diimine)(CO)3Cl complexes with nanographene ligands that contain a polycyclic aromatic hydrocarbon moiety through a pyrazinyl linkage. We show that the reduction potentials of the complexes depend not only on conjugation size but also on shape and geometry of the ligands, revealing another parameter in tuning the redox properties of the electrocatalysts. In addition, our work reveals a compromise between reduction potentials and activation of this class of electrocatalysts.
Collapse
Affiliation(s)
- Sruthy K Chandy
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Scott Bowers
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Krishnan Raghavachari
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Liang-Shi Li
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| |
Collapse
|
10
|
Akher FB, Shu Y, Varga Z, Truhlar DG. Semiclassical Multistate Dynamics for Six Coupled 5A' States of O + O 2. J Chem Theory Comput 2023. [PMID: 37441750 DOI: 10.1021/acs.jctc.3c00517] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/15/2023]
Abstract
Dynamics simulations of high-energy O2-O collisions play an important role in simulating thermal energy content and heat flux in flows around hypersonic vehicles. To carry out such dynamics simulations efficiently requires accurate global potential energy surfaces and (in most algorithms) state couplings for many energetically accessible electronic states. The ability to treat collisions involving many coupled electronic states has been a challenge for decades. Very recently, a new diabatization method, the parametrically managed diabatization by deep neural network (PM-DDNN), has been developed. The PM-DDNN method uses a deep neural network architecture with an activation function parametrically dependent on input data to discover and fit the diabatic potential energy matrix (DPEM) as a function of geometry, and the adiabatic potential energy surfaces are obtained by diagonalization of a small matrix with analytic matrix elements. Here, we applied the PM-DDNN method to the six lowest-energy potential energy surfaces in the 5A' manifold of O3 to perform simultaneous diabatization and fitting; the data are obtained by extended multistate complete-active-space second-order perturbation theory. We then used the adiabatic surfaces for dynamics calculations with three methods: coherent switching with decay of mixing (CSDM), curvature-driven CSDM (κCSDM), and electronically curvature-driven CSDM (eκCSDM). The κCSDM calculations require only adiabatic potential energies and gradients. The three dynamical methods are in good agreement. We then calculated electronically nonadiabatic, electronically inelastic, and dissociative cross sections for seven initial collision energies, five initial vibrational levels, and four initial rotational levels. Trends in the electronically inelastic cross sections as functions of the initial collision energy and vibrational level were rationalized in terms of the coordinate ranges where the gaps between the second and third potential energy surfaces are small.
Collapse
Affiliation(s)
- Farideh Badichi Akher
- Department of Chemistry, Chemical Theory Center and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Yinan Shu
- Department of Chemistry, Chemical Theory Center and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Zoltan Varga
- Department of Chemistry, Chemical Theory Center and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
11
|
Rossomme E, Cunha LA, Li W, Chen K, McIsaac AR, Head-Gordon T, Head-Gordon M. The Good, the Bad, and the Ugly: Pseudopotential Inconsistency Errors in Molecular Applications of Density Functional Theory. J Chem Theory Comput 2023; 19:2827-2841. [PMID: 37156013 DOI: 10.1021/acs.jctc.3c00089] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/10/2023]
Abstract
The pseudopotential (PP) approximation is one of the most common techniques in computational chemistry. Despite its long history, the development of custom PPs has not tracked with the explosion of different density functional approximations (DFAs). As a result, the use of PPs with exchange/correlation models for which they were not developed is widespread, although this practice is known to be theoretically unsound. The extent of PP inconsistency errors (PPIEs) associated with this practice has not been systematically explored across the types of energy differences commonly evaluated in chemical applications. We evaluate PPIEs for a number of PPs and DFAs across 196 chemically relevant systems of both transition-metal and main-group elements, as represented by the W4-11, TMC34, and S22 data sets. Near the complete basis set limit, these PPs are found to cleanly approach all-electron (AE) results for noncovalent interactions but introduce root-mean-squared errors (RMSEs) upwards of 15 kcal mol-1 into predictions of covalent bond energies for a number of popular DFAs. We achieve significant improvements through the use of empirical atom- and DFA-specific PP corrections, indicating considerable systematicity of the PPIEs. The results of this work have implications for chemical modeling in both molecular contexts and for DFA design, which we discuss.
Collapse
Affiliation(s)
- Elliot Rossomme
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Leonardo A Cunha
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Wanlu Li
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Kaixuan Chen
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Alexandra R McIsaac
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Teresa Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
- Department of Bioengineering, University of California, Berkeley, California 94720, United States
- Department of Chemical and Biomolecular Engineering, University of California, Berkeley, California 94720, United States
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
12
|
Tao C, Yang J, Hong Q, Sun Q, Li J. Global and Full-Dimensional Potential Energy Surfaces of the N 2 + O 2 Reaction for Hyperthermal Collisions. J Phys Chem A 2023; 127:4027-4042. [PMID: 37128765 DOI: 10.1021/acs.jpca.3c01065] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/03/2023]
Abstract
The energy transfer, dissociations, and chemical reactions between O2 and N2 play an important role in the re-entry process of aircraft and many atmospheric, combustion, and plasma processes. Recently, Varga et al. (J. Chem. Phys., 2016, 144, 024310) developed a full-dimensional high-precision potential energy surface (PES) of the ground triplet electronic state for the O2 and N2 system based on ca. 55,000 data points, whose energies were calculated by multi-state complete-active-space second-order perturbation theory/minimally augmented correlation-consistent polarized valence triple-zeta electronic structure calculations plus dynamically scaled external correlation. The fitting function adopted the many-body expansion form with the four-body interactions fitted by the permutationally invariant polynomial in terms of bond-order functions of the six interatomic distances (MB-PIP). In this work, we refit the PES of the N2O2 system by two methods based on the same data set that was used by Varga et al. The first uses the permutation invariant polynomial-neural network (PIP-NN) method to fit the entire energy of the 55,000 data points. In the second approach, the PIP-NN method is used to fit only the four-body interaction component, a similar treatment in the MB-PIP method, and the resulting PES is named MB-PIP-NN. Then, the performances of these new PESs and the MB-PIP PES are comprehensively and systematically compared, such as comparisons of various scans, properties of stationary points, and dynamics simulations. Possible improvements for the PES of N2O2 are suggested. A more reliable PES of the system can be constructed in terms of data sampling range, electronic structure calculation level, and fitting method for high-temperature calculation and simulation in the future.
Collapse
Affiliation(s)
- Chun Tao
- School of Chemistry and Chemical Engineering and Chongqing Key Laboratory of Theoretical and Computational Chemistry, Chongqing University, Chongqing 401331, China
| | - Jiawei Yang
- School of Chemistry and Chemical Engineering and Chongqing Key Laboratory of Theoretical and Computational Chemistry, Chongqing University, Chongqing 401331, China
| | - Qizhen Hong
- State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy of Sciences, Beijing 100190, China
| | - Quanhua Sun
- State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy of Sciences, Beijing 100190, China
| | - Jun Li
- School of Chemistry and Chemical Engineering and Chongqing Key Laboratory of Theoretical and Computational Chemistry, Chongqing University, Chongqing 401331, China
| |
Collapse
|
13
|
DFT and TD-DFT study of hydrogen bonded complexes of aspartic acid and n water (n = 1 and 2). J Mol Model 2023; 29:94. [PMID: 36905452 DOI: 10.1007/s00894-023-05500-z] [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: 01/18/2023] [Accepted: 03/02/2023] [Indexed: 03/12/2023]
Abstract
CONTEXT Hydrogen bonds (HB) influence the conformational preferences of biomolecules and their optical and electronic properties. The directional interaction of molecules of water can be a prototype to understand the effects of HBs on biomolecules. Among the neurotransmitters (NT), L-aspartic acid (ASP) stands out due to its importance in health and as a precursor of several biomolecules. As it presents different functional groups and readily forms inter- and intramolecular HBs, ASP can be considered a prototype for understanding the behavior of NTs when interacting by HB with other substances. Although several theoretical studies have been performed in the past on isolated ASP and its formed complexes with water, both in gas and liquid phases, using DFT and TD-DFT formalisms, these works did not perform large basis set calculations or study electronic transitions of ASP-water complexes. We investigated the HB interactions in complexes of ASP and water molecules. The results show that the interactions between the carboxylic groups of ASP with water molecules, forming cyclic structures with two HBs, lead to more stable and less polar complexes than other conformers formed between water and the NH2 group. It was observed that there is a relationship between the deviation in the UV-Vis absorption band of the ASP and the interactions of water with the HOMO and LUMO orbitals with the stabilization/destabilization of the S1 state to the S0 of the complexes. However, in some cases, such as 1:1 complex ASP-W2, this analysis may be inaccurate due to small changes in ΔE. METHODS We studied the landscapes of the ground state surface of different conformers of isolated L-ASP and the L-ASP-(H2O)n complexes (n = 1 and 2) using the DFT formalism, with the B3LYP functional, and six different basis sets: 6-31 + + G(d,p), 6-311 + + G(d,p), D95 + + (d,p), D95V + + (d,p), cc-pVDZ, and, cc-pVTZ basis sets. The cc-pVTZ basis set provides the minimum energy of all conformers, and therefore, we performed the analysis with this basis set. We evaluated the stabilization of the ASP and complexes using the minimum ground state energy, corrected by the zero point energy and the interaction energy between the ASP and the water molecules. We also calculated the vertical electronic transitions S1 ← S0, and their properties using the TD-DFT formalism at B3LYP/cc-pVTZ level with the optimized geometries for S0 state with the same basis set. For the analysis of the vertical transitions of isolated ASP and the ASP-(H2O)n complexes, we calculated the electrostatic energy in the S0 and S1 states. We performed the calculations with the Gaussian 09 software package. We used the VMD software package to visualize the geometries and shapes of the molecule and complexes.
Collapse
|
14
|
Ning J, Truhlar DG. Electronic Excitation of ortho-Fluorothiophenol. J Phys Chem A 2023; 127:1469-1474. [PMID: 36749710 DOI: 10.1021/acs.jpca.2c08674] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
ortho-Fluorothiophenol (o-FTP) photodissociates through the well-known πσ* process. The fluorine atom of o-FTP introduces a feature in the photodissociation of o-FTP that does not occur in most other πσ* processes because the fluorine atom can form a hydrogen bond with the hydrogen atom of the SH group. Theoretical computations can serve as a good way to study these reactions because they usually proceed very quickly, and the current spectroscopies cannot probe the details of the processes as thoroughly as theory can. Here we use completely renormalized equation-of-motion coupled cluster theory with single and double excitations and a quasiperturbative treatment of connected triple excitations (CR-EOM-CCSD(T)) and quasidegenerate perturbation theory, in particular extended multistate complete-active-space second-order perturbation theory (XMS- CASPT2), to calculate the four lowest singlet states of o-FTP and hybrid density functional theory to optimize the geometries of the two lowest singlet states. We find that ten active electrons in nine active orbitals are sufficient to provide a good reference function for all four states. We find that the ground electronic state and the first excited singlet state both exhibit strongly bent hydrogen bonds. We also use density functional theory with the Tamm-Dancoff approximation and the SMD solvation model to successfully simulate the electronic spectrum of o-FTP in n-hexane solvent.
Collapse
Affiliation(s)
- Jiaxin Ning
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| |
Collapse
|
15
|
Properties of Gaseous Deprotonated L-Cysteine S-Sulfate Anion [cysS-SO 3] -: Intramolecular H-Bond Network, Electron Affinity, Chemically Active Site, and Vibrational Fingerprints. Int J Mol Sci 2023; 24:ijms24021682. [PMID: 36675196 PMCID: PMC9862062 DOI: 10.3390/ijms24021682] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2022] [Revised: 01/10/2023] [Accepted: 01/11/2023] [Indexed: 01/18/2023] Open
Abstract
L-cysteine S-sulfate, Cys-SSO3H, and their derivatives play essential roles in biological chemistry and pharmaceutical synthesis, yet their intrinsic molecular properties have not been studied to date. In this contribution, the deprotonated anion [cysS-SO3]- was introduced in the gas phase by electrospray and characterized by size-selected, cryogenic, negative ion photoelectron spectroscopy. The electron affinity of the [cysS-SO3]• radical was determined to be 4.95 ± 0.10 eV. In combination with theoretical calculations, it was found that the most stable structure of [cysS-SO3]- (S1) is stabilized via three intramolecular hydrogen bonds (HBs); i.e., one O-H⋯⋯N between the -COOH and -NH2 groups, and two N-H⋯⋯O HBs between -NH2 and -SO3, in which the amino group serves as both HB acceptor and donor. In addition, a nearly iso-energetic conformer (S2) with the formation of an O-H⋯⋯N-H⋯⋯O-S chain-type binding motif competes with S1 in the source. The most reactive site of the molecule susceptible for electrophilic attacks is the linkage S atom. Theoretically predicted infrared spectra indicate that O-H and N-H stretching modes are the fingerprint region (2800 to 3600 cm-1) to distinguish different isomers. The obtained information lays out a foundation to better understand the transformation and structure-reactivity correlation of Cys-SSO3H in biologic settings.
Collapse
|
16
|
Recio P, Alessandrini S, Vanuzzo G, Pannacci G, Baggioli A, Marchione D, Caracciolo A, Murray VJ, Casavecchia P, Balucani N, Cavallotti C, Puzzarini C, Barone V. Intersystem crossing in the entrance channel of the reaction of O( 3P) with pyridine. Nat Chem 2022; 14:1405-1412. [PMID: 36175514 DOI: 10.1038/s41557-022-01047-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2021] [Accepted: 08/25/2022] [Indexed: 01/04/2023]
Abstract
Two quantum effects can enable reactions to take place at energies below the barrier separating reactants from products: tunnelling and intersystem crossing between coupled potential energy surfaces. Here we show that intersystem crossing in the region between the pre-reactive complex and the reaction barrier can control the rate of bimolecular reactions for weakly coupled potential energy surfaces, even in the absence of heavy atoms. For O(3P) plus pyridine, a reaction relevant to combustion, astrochemistry and biochemistry, crossed-beam experiments indicate that the dominant products are pyrrole and CO, obtained through a spin-forbidden ring-contraction mechanism. The experimental findings are interpreted-by high-level quantum-chemical calculations and statistical non-adiabatic computations of branching fractions-in terms of an efficient intersystem crossing occurring before the high entrance barrier for O-atom addition to the N-atom lone pair. At low to moderate temperatures, the computed reaction rates prove to be dominated by intersystem crossing.
Collapse
Affiliation(s)
- Pedro Recio
- Dipartimento di Chimica, Biologia e Biotecnologie, Università degli Studi di Perugia, Perugia, Italy
| | - Silvia Alessandrini
- Scuola Normale Superiore, Pisa, Italy
- Dipartimento di Chimica 'Giacomo Ciamician', University of Bologna, Bologna, Italy
| | - Gianmarco Vanuzzo
- Dipartimento di Chimica, Biologia e Biotecnologie, Università degli Studi di Perugia, Perugia, Italy
| | - Giacomo Pannacci
- Dipartimento di Chimica, Biologia e Biotecnologie, Università degli Studi di Perugia, Perugia, Italy
| | - Alberto Baggioli
- Dipartimento di Chimica, Materiali e Ingegneria Chimica 'Giulio Natta', Politecnico di Milano, Milan, Italy
| | - Demian Marchione
- Dipartimento di Chimica, Biologia e Biotecnologie, Università degli Studi di Perugia, Perugia, Italy
| | - Adriana Caracciolo
- Dipartimento di Chimica, Biologia e Biotecnologie, Università degli Studi di Perugia, Perugia, Italy
- Department of Aerospace Engineering Sciences, University of Colorado, Boulder, CO, USA
| | - Vanessa J Murray
- Dipartimento di Chimica, Biologia e Biotecnologie, Università degli Studi di Perugia, Perugia, Italy
- Montana State University, Bozeman, MT, USA
| | - Piergiorgio Casavecchia
- Dipartimento di Chimica, Biologia e Biotecnologie, Università degli Studi di Perugia, Perugia, Italy
| | - Nadia Balucani
- Dipartimento di Chimica, Biologia e Biotecnologie, Università degli Studi di Perugia, Perugia, Italy.
| | - Carlo Cavallotti
- Dipartimento di Chimica, Materiali e Ingegneria Chimica 'Giulio Natta', Politecnico di Milano, Milan, Italy.
| | - Cristina Puzzarini
- Dipartimento di Chimica 'Giacomo Ciamician', University of Bologna, Bologna, Italy.
| | | |
Collapse
|
17
|
Mitchell EC, Scott TR, Bao JJ, Truhlar DG. Application of Multiconfiguration Pair-Density Functional Theory to the Diels–Alder Reaction. J Phys Chem A 2022; 126:8834-8843. [DOI: 10.1021/acs.jpca.2c06433] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
Affiliation(s)
- Erica C. Mitchell
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota − Twin Cities, Minneapolis, Minnesota 55455-0931, United States
- Department of Chemistry, University of Georgia, Athens, Georgia 30602, United States
| | - Thais R. Scott
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota − Twin Cities, Minneapolis, Minnesota 55455-0931, United States
| | - Jie J. Bao
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota − Twin Cities, Minneapolis, Minnesota 55455-0931, United States
| | - Donald G. Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota − Twin Cities, Minneapolis, Minnesota 55455-0931, United States
| |
Collapse
|
18
|
Chandy SK, Bowers SA, Yin M, Liu L, Raghavachari K, Li LS. Proton-Coupled, Low-Energy Pathway for Electrocatalytic CO 2 Reduction at Re(Diimine) Complexes with a Conjugated Pyrazinyl Moiety. Inorg Chem 2022; 61:17505-17514. [DOI: 10.1021/acs.inorgchem.2c02400] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Sruthy K. Chandy
- Department of Chemistry, Indiana University, Bloomington, Indiana47405, United States
| | - Scott A. Bowers
- Department of Chemistry, Indiana University, Bloomington, Indiana47405, United States
| | - Minyang Yin
- Department of Chemistry, Indiana University, Bloomington, Indiana47405, United States
| | - Lu Liu
- Department of Chemistry, Indiana University, Bloomington, Indiana47405, United States
| | - Krishnan Raghavachari
- Department of Chemistry, Indiana University, Bloomington, Indiana47405, United States
| | - Liang-shi Li
- Department of Chemistry, Indiana University, Bloomington, Indiana47405, United States
| |
Collapse
|
19
|
Nottoli G, Ballotta B, Rampino S. Local Charge-Displacement Analysis: Targeting Local Charge-Flows in Complex Intermolecular Interactions. J Chem Phys 2022; 157:084107. [DOI: 10.1063/5.0095142] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Charge-displacement (CD) analysis has recently proven to be a simple and powerful scheme for quantitatively analyzing the profile of the charge redistribution occurring upon intermolecular interactions along a given interaction axis. However, when two molecular fragments bind through complex interactions involving multiple concurrent charge flows, ordinary CD analysis is capable of providing only an averaged picture of the related charge-flow profiles and no detailed information on each of them. In this article, we combine CD analysis with a Hirshfeld partitioning of the molecular charge redistribution for a local analysis on focused portions of the molecule, allowing for a detailed characterization of one charge flow at a time. The resulting scheme - the local charge-displacement (LCD) analysis - is tested on the intriguing case of the dimethyl sulfide (DMS)-sulfur dioxide (SO2) complex, characterized by concurrent charge flows relating to a sulfur-sulfur homo-chalcogen interaction and a pair of hydrogen bonds. The LCD scheme is then applied to the analysis of multiple hydrogen bonding in the acetic acid dimer, of base-pairing interactions in DNA, and of ambifunctional hydrogen bonding in the ammonia-pyridine complex.
Collapse
|
20
|
Zheng P, Gan Z, Zhou C, Su P, Wu W. λ-DFVB(U): A hybrid density functional valence bond method based on unpaired electron density. J Chem Phys 2022; 156:204103. [DOI: 10.1063/5.0091592] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this paper, a hybrid density functional valence bond method based on unpaired electron density, called λ-DFVB(U), is presented, which is a combination of the valence bond self-consistent field (VBSCF) method and Kohn–Sham density functional theory. In λ-DFVB(U), the double-counting error of electron correlation is mitigated by a linear decomposition of the electron–electron interaction using a parameter λ, which is a function of an index based on the number of effectively unpaired electrons. In addition, λ-DFVB(U) is based on the approximation that correlation functionals in KS-DFT only cover dynamic correlation and exchange functionals mimic some amount of static correlation. Furthermore, effective spin densities constructed from unpaired density are used to address the symmetry dilemma problem in λ-DFVB(U). The method is applied to test calculations of atomization energies, atomic excitation energies, and reaction barriers. It is shown that the accuracy of λ-DFVB(U) is comparable to that of CASPT2, while its computational cost is approximately the same as VBSCF.
Collapse
Affiliation(s)
- Peikun Zheng
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, Department of Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| | - Zixi Gan
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, Department of Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| | - Chen Zhou
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, Department of Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| | - Peifeng Su
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, Department of Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| | - Wei Wu
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, Department of Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| |
Collapse
|
21
|
Ceselin G, Salta Z, Bloino J, Tasinato N, Barone V. Accurate Quantum Chemical Spectroscopic Characterization of Glycolic Acid: A Route Toward its Astrophysical Detection. J Phys Chem A 2022; 126:2373-2387. [PMID: 35384666 PMCID: PMC9036519 DOI: 10.1021/acs.jpca.2c01419] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2022] [Revised: 03/18/2022] [Indexed: 11/29/2022]
Abstract
The first step to shed light on the abiotic synthesis of biochemical building blocks, and their further evolution toward biological systems, is the detection of the relevant species in astronomical environments, including earthlike planets. To this end, the species of interest need to be accurately characterized from structural, energetic, and spectroscopic viewpoints. This task is particularly challenging when dealing with flexible systems, whose spectroscopic signature is ruled by the interplay of small- and large-amplitude motions (SAMs and LAMs, respectively) and is further tuned by the conformational equilibrium. In such instances, quantum chemical (QC) calculations represent an invaluable tool for assisting the interpretation of laboratory measurements or even observations. In the present work, the role of QC results is illustrated with reference to glycolic acid (CH2OHCOOH), a molecule involved in photosynthesis and plant respiration and a precursor of oxalate in humans, which has been detected in the Murchison meteorite but not yet in the interstellar medium or in planetary atmospheres. In particular, the equilibrium structure of the lowest-energy conformer is derived by employing the so-called semiexperimental approach. Then, accurate yet cost-effective QC calculations relying on composite post-Hartree-Fock schemes and hybrid coupled-cluster/density functional theory approaches are used to predict the structural and ro-vibrational spectroscopic properties of the different conformers within the framework of the second-order vibrational perturbation theory. A purposely tailored discrete variable representation anharmonic approach is used to treat the LAMs related to internal rotations. The computed spectroscopic data, particularly those in the infrared region, complement the available experimental investigations, thus enhancing the possibility of an astronomical detection of this molecule.
Collapse
Affiliation(s)
- Giorgia Ceselin
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126, Pisa, Italy
| | - Zoi Salta
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126, Pisa, Italy
| | - Julien Bloino
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126, Pisa, Italy
| | - Nicola Tasinato
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126, Pisa, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126, Pisa, Italy
| |
Collapse
|
22
|
Barclay AJ, McKellar ARW, Pietropolli Charmet A, Moazzen-Ahmadi N. Weakly-bound clusters of atmospheric molecules: infrared spectra and structural calculations of (CO 2) n-(CO) m-(N 2) p, ( n, m, p) = (2,1,0), (2,0,1), (1,2,0), (1,0,2), (1,1,1), (1,3,0), (1,0,3), (1,2,1), (1,1,2). Phys Chem Chem Phys 2022; 24:7231-7242. [PMID: 35266943 DOI: 10.1039/d2cp00263a] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Structural calculations and high-resolution infrared spectra are reported for trimers and tetramers containing CO2 together with CO and/or N2. Among the 9 clusters studied here, only (CO2)2-CO was previously observed by high-resolution spectroscopy. The spectra, which occur in the region of the ν3 fundamental of CO2 (≈2350 cm-1), were recorded using a tunable optical parametric oscillator source to probe a pulsed supersonic slit jet expansion. The trimers (CO2)2-CO and (CO2)2-N2 have structures in which the CO or N2 is aligned along the symmetry axis of a staggered side-by-side CO2 dimer unit. The observation of two fundamental bands for (CO2)2-CO and (CO2)2-N2 shows that this CO2 dimer unit is non-planar, unlike (CO2)2 itself. For the trimers CO2-(CO)2 and CO2-(N2)2, the CO or N2 monomers occupy equivalent positions in the 'equatorial plane' of the CO2, pointing toward its C atom. To form the tetramers CO2-(CO)3 and CO2-(N2)3, a third CO or N2 monomer is then added off to the 'side' of the first two. In the mixed tetramers CO2-(CO)2-N2 and CO2-CO-(N2)2, this 'side' position is taken by N2 and not CO. In addition to the fundamental bands, combination bands are also observed for (CO2)2-CO, CO2-(CO)2, and CO2-(N2)2, yielding some information about their low-frequency intermolecular vibrations.
Collapse
Affiliation(s)
- A J Barclay
- Department of Physics and Astronomy, University of Calgary, 2500 University Drive North West, Calgary, Alberta T2N 1N4, Canada.
| | - A R W McKellar
- National Research Council of Canada, Ottawa, Ontario K1A 0R6, Canada
| | - A Pietropolli Charmet
- Dipartimento di ScienzeMolecolari e Nanosistemi, Università Ca' Foscari Venezia, Via Torino 155, I-30172, Mestre, Venezia, Italy
| | - N Moazzen-Ahmadi
- Department of Physics and Astronomy, University of Calgary, 2500 University Drive North West, Calgary, Alberta T2N 1N4, Canada.
| |
Collapse
|
23
|
Li X, Spada L, Alessandrini S, Zheng Y, Lengsfeld KG, Grabow J, Feng G, Puzzarini C, Barone V. Stacked but not Stuck: Unveiling the Role of π→π* Interactions with the Help of the Benzofuran-Formaldehyde Complex. Angew Chem Int Ed Engl 2022; 61:e202113737. [PMID: 34697878 PMCID: PMC9298890 DOI: 10.1002/anie.202113737] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2021] [Indexed: 12/14/2022]
Abstract
The 1:1 benzofuran-formaldehyde complex has been chosen as model system for analyzing π→π* interactions in supramolecular organizations involving heteroaromatic rings and carbonyl groups. A joint "rotational spectroscopy-quantum chemistry" strategy unveiled the dominant role of π→π* interactions in tuning the intermolecular interactions of such adduct. The exploration of the intermolecular potential energy surface led to the identification of 14 low-energy minima, with 4 stacked isomers being more stable than those linked by hydrogen bond or lone-pair→π interactions. All energy minima are separated by loose transition states, thus suggesting an effective relaxation to the global minimum under the experimental conditions. This expectation has been confirmed by the experimental detection of only one species, which was unambiguously assigned owing to the computation of accurate spectroscopic parameters and the characterization of 11 isotopologues. The large number of isotopic species opened the way to the determination of the first semi-experimental equilibrium structure for a molecular complex of such a dimension.
Collapse
Affiliation(s)
- Xiaolong Li
- School of Chemistry and Chemical EngineeringChongqing UniversityDaxuecheng South Rd. 55Chongqing401331China
| | - Lorenzo Spada
- Scuola Normale SuperiorePiazza dei Cavalieri 756126PisaItaly
| | - Silvia Alessandrini
- Scuola Normale SuperiorePiazza dei Cavalieri 756126PisaItaly
- Dipartimento di Chimica “Giacomo Ciamician”University of BolognaVia F. Selmi 240126BolognaItaly
| | - Yang Zheng
- School of Chemistry and Chemical EngineeringChongqing UniversityDaxuecheng South Rd. 55Chongqing401331China
| | - Kevin Gregor Lengsfeld
- Institut für Physikalische Chemie and ElektrochemieGottfried Wilhelm Leibniz Universität HannoverCallinstrasse 3A30167HannoverGermany
| | - Jens‐Uwe Grabow
- Institut für Physikalische Chemie and ElektrochemieGottfried Wilhelm Leibniz Universität HannoverCallinstrasse 3A30167HannoverGermany
| | - Gang Feng
- School of Chemistry and Chemical EngineeringChongqing UniversityDaxuecheng South Rd. 55Chongqing401331China
| | - Cristina Puzzarini
- Dipartimento di Chimica “Giacomo Ciamician”University of BolognaVia F. Selmi 240126BolognaItaly
| | - Vincenzo Barone
- Scuola Normale SuperiorePiazza dei Cavalieri 756126PisaItaly
| |
Collapse
|
24
|
Li X, Spada L, Alessandrini S, Zheng Y, Lengsfeld KG, Grabow J, Feng G, Puzzarini C, Barone V. Gestapelt, nicht geklebt: Enthüllung der π→π*‐Wechselwirkung mithilfe des Benzofuran‐Formaldehyd‐Komplexes. Angew Chem Int Ed Engl 2022. [DOI: 10.1002/ange.202113737] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Affiliation(s)
- Xiaolong Li
- School of Chemistry and Chemical Engineering Chongqing University Daxuecheng South Rd. 55 Chongqing 401331 China
| | - Lorenzo Spada
- Scuola Normale Superiore Piazza dei Cavalieri 7 56126 Pisa Italien
| | - Silvia Alessandrini
- Scuola Normale Superiore Piazza dei Cavalieri 7 56126 Pisa Italien
- Dipartimento di Chimica “Giacomo Ciamician” Università di Bologna Via F. Selmi 2 40126 Bologna Italien
| | - Yang Zheng
- School of Chemistry and Chemical Engineering Chongqing University Daxuecheng South Rd. 55 Chongqing 401331 China
| | - Kevin Gregor Lengsfeld
- Institut für Physikalische Chemie und Elektrochemie Gottfried Wilhelm Leibniz Universität Hannover Callinstraße 3A 30167 Hannover Deutschland
| | - Jens‐Uwe Grabow
- Institut für Physikalische Chemie und Elektrochemie Gottfried Wilhelm Leibniz Universität Hannover Callinstraße 3A 30167 Hannover Deutschland
| | - Gang Feng
- School of Chemistry and Chemical Engineering Chongqing University Daxuecheng South Rd. 55 Chongqing 401331 China
| | - Cristina Puzzarini
- Dipartimento di Chimica “Giacomo Ciamician” Università di Bologna Via F. Selmi 2 40126 Bologna Italien
| | - Vincenzo Barone
- Scuola Normale Superiore Piazza dei Cavalieri 7 56126 Pisa Italien
| |
Collapse
|
25
|
Barone V, Fusè M, Pinto SMV, Tasinato N. A Computational Journey across Nitroxide Radicals: From Structure to Spectroscopic Properties and Beyond. Molecules 2021; 26:molecules26237404. [PMID: 34885980 PMCID: PMC8659111 DOI: 10.3390/molecules26237404] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2021] [Revised: 11/25/2021] [Accepted: 12/02/2021] [Indexed: 11/17/2022] Open
Abstract
Nitroxide radicals are characterized by a long-lived open-shell electronic ground state and are strongly sensitive to the chemical environment, thus representing ideal spin probes and spin labels for paramagnetic biomolecules and materials. However, the interpretation of spectroscopic parameters in structural and dynamic terms requires the aid of accurate quantum chemical computations. In this paper we validate a computational model rooted into double-hybrid functionals and second order vibrational perturbation theory. Then, we provide reference quantum chemical results for the structures, vibrational frequencies and other spectroscopic features of a large panel of nitroxides of current biological and/or technological interest.
Collapse
|
26
|
Varga Z, Truhlar DG. Potential energy surface for high-energy N + N 2 collisions. Phys Chem Chem Phys 2021; 23:26273-26284. [PMID: 34787127 DOI: 10.1039/d1cp04373k] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Potential energy surface calculations yield physical insight into the structure of intermediates and the dynamics of molecular collisions, and they are the first step toward molecular simulations that provide physical insight into energy transfer, reaction, and dissociation probabilities. The potential energy surface for high-energy collisions of N2 with N can be used for modeling chemical dynamics and energy transfer in atmospheric shock waves. Here we present an analytic ground-state. (4A'') potential energy surface for N3 that governs electronically adiabatic collisions of N2(1Σ+g) with N(4S). The fitted surface consists of a pairwise potential based on an accurate diatomic potential energy curve plus a connected permutationally invariant polynomial (PIP) in mixed-exponential-Gaussian bond order variables (MEGs) for the three-body part. The three-body fit is based on multireference complete active space second order perturbation theory (CASPT2) calculations. The quality of the quartet N3 fit is comparable to that for a previous fit of the NO2 potential. We characterize two local minima of N3, two tight transition structures, two van der Waals geometries, and the noncollinear reaction path for the symmetric exchange reaction. The nonreactive approach of an N atom to N2 along the perpendicular bisector is more repulsive than the collinear reproach, but plots of the force on the bond versus the potential energy at the distance of closest approach allow us to infer that vibrational energy transfer should occur much more readily in high-energy collinear collisions than in high-energy perpendicular-bisector collisions.
Collapse
Affiliation(s)
- Zoltan Varga
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA.
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA.
| |
Collapse
|
27
|
Wang P, Shu C, Ye H, Biczysko M. Structural and Energetic Properties of Amino Acids and Peptides Benchmarked by Accurate Theoretical and Experimental Data. J Phys Chem A 2021; 125:9826-9837. [PMID: 34752094 DOI: 10.1021/acs.jpca.1c06504] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
Structural, energetic, and spectroscopic data derived in this work aim at the setup of an "experimentally validated" database for amino acids and polypeptides conformers. First, the "cheap" composite scheme (ChS, CCSD(T)/(CBS+CV)MP2) is tested for evaluation of conformational energies of all eight stable conformers of glycine, by comparing to the more accurate CCSD(T)/CBS+CV computations (Phys. Chem. Chem. Phys. 2013, 15, 10094-10111 and J Mol. Model. 2020, 26, 129). The recently proposed jun-ChS (J. Chem. Theory and Comput. 2020, 16, 988-1006), employing the jun-cc-pVnZ basis set family for CCSD(T) computations and CBS extrapolation, yields conformational energies accurate to 0.2 kJ·mol-1, at reduced computational cost with respect to aug-ChS employing aug-cc-pVnZ basis sets. The jun-ChS composite scheme is further applied to derive conformational energies for three dipeptide analogues Ac-Gly-NH2, Ac-Ala-NH2, and Gly-Gly. Finally, dipeptide conformational energies and semiexperimental equilibrium rotational constants along with the CCSD(T)/(CBS+CV)MP2 structural parameters (J. Phys. Chem. Lett. 2014, 5, 534-540) stand as the reference for benchmarking of selected density functional methodologies. The double-hybrid functionals B2-PLYP-D3(BJ) and DSD-PBEP86, perform best for structural and energetic characterization of all dipeptide analogues. From hybrid functionals CAM-B3LYP-D3(BJ) and ωB97X-D3(BJ) represent promising methods applicable for larger peptide-based systems for which computations with double-hybrid functionals are not feasible.
Collapse
Affiliation(s)
- Ping Wang
- International Centre for Quantum and Molecular Structures, Physics Department, College of Science, Shanghai University, 99 Shangda Road, Shanghai 200444, China
| | - Chong Shu
- International Centre for Quantum and Molecular Structures, Physics Department, College of Science, Shanghai University, 99 Shangda Road, Shanghai 200444, China
| | - Hexu Ye
- International Centre for Quantum and Molecular Structures, Physics Department, College of Science, Shanghai University, 99 Shangda Road, Shanghai 200444, China
| | - Malgorzata Biczysko
- International Centre for Quantum and Molecular Structures, Physics Department, College of Science, Shanghai University, 99 Shangda Road, Shanghai 200444, China
| |
Collapse
|
28
|
Wang Y, Cheng X. State-to-State Transition Study of the Exchange Reaction for N( 4S) and O 2(X 3Σ g-) Collision by Quasi-Classical Trajectory. J Phys Chem A 2021; 125:9318-9326. [PMID: 34652155 DOI: 10.1021/acs.jpca.1c06386] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Based on the new 2A' and 4A' potential energy surfaces of NO2 fitted by Varga et al., we conducted a quasi-classical trajectory study on the N(4S) +O2(X3Σg- ) → NO(2Π) + O(3P) reaction, focusing on the high vibrational state up to ν = 25. For different rovibrational states of O2, within the relative translational energy (Ec) range of 0.1-30 eV, the total exchange cross section (ECS) is calculated, and it is found that the initial relative translational energy and vibration excitation have a significant effect on ECSs, while rotational excitation has little influence; the rate coefficient of the high rovibrational state of O2 molecules at high temperatures is studied, and it is found that when the vibrational level ν of O2 is in the range of 0-15, the value of log10 k(T, ν, j) with the vibrational level ν is almost linear, while when ν is greater than 15, it becomes gentle with the increase in ν. Finally, the state-to-state rate coefficients are calculated; our results supply the advantageous state-to-state process data in the NO2 system, and they are useful for further studying the related hypersonic gas flow at very high temperature.
Collapse
Affiliation(s)
- Yan Wang
- Institute of Atomic and Molecular Physics, Sichuan University, Chengdu 610065, China
| | - Xinlu Cheng
- Institute of Atomic and Molecular Physics and Key Laboratory of High Energy Density Physics and Technology of Ministry of Education, Sichuan University, Chengdu 610065, China
| |
Collapse
|
29
|
Sagresti L, Rampino S. Charge-Flow Profiles along Curvilinear Paths: A Flexible Scheme for the Analysis of Charge Displacement upon Intermolecular Interactions. Molecules 2021; 26:molecules26216409. [PMID: 34770822 PMCID: PMC8586930 DOI: 10.3390/molecules26216409] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2021] [Revised: 10/19/2021] [Accepted: 10/21/2021] [Indexed: 11/26/2022] Open
Abstract
The Charge-Displacement (CD) analysis has proven to be a powerful tool for a quantitative characterization of the electron-density flow occurring upon chemical bonding along a suitably chosen interaction axis. In several classes of interesting intermolecular interactions, however, an interaction axis cannot be straightforwardly defined, and the CD analysis loses consistency and usefulness. In this article, we propose a general, flexible reformulation of the CD analysis capable of providing a quantitative view of the charge displacement along custom curvilinear paths. The new scheme naturally reduces to ordinary CD analysis if the path is chosen to be a straight line. An implementation based on a discrete sampling of the electron densities and a Voronoi space partitioning is described and shown in action on two test cases of a metal-carbonyl and a pyridine-ammonia complex.
Collapse
Affiliation(s)
- Luca Sagresti
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy;
- Istituto Nazionale di Fisica Nucleare (INFN)—Sezione di Pisa, Largo Bruno Pontecorvo 3, 56127 Pisa, Italy
| | - Sergio Rampino
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy;
- Istituto Nazionale di Fisica Nucleare (INFN)—Sezione di Pisa, Largo Bruno Pontecorvo 3, 56127 Pisa, Italy
- Correspondence:
| |
Collapse
|
30
|
Sheng M, Silvestrini F, Biczysko M, Puzzarini C. Structural and Vibrational Properties of Amino Acids from Composite Schemes and Double-Hybrid DFT: Hydrogen Bonding in Serine as a Test Case. J Phys Chem A 2021; 125:9099-9114. [PMID: 34623165 DOI: 10.1021/acs.jpca.1c06993] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
The structures, relative stabilities, and vibrational wavenumbers of the two most stable conformers of serine, stabilized by the O-H···N, O-H···O═C and N-H···O-H intramolecular hydrogen bonds, have been evaluated by means of state-of-the-art composite schemes based on coupled-cluster (CC) theory. The so-called "cheap" composite approach (CCSD(T)/(CBS+CV)MP2) allowed determination of accurate equilibrium structures and harmonic vibrational wavenumbers, also pointing out significant corrections beyond the CCSD(T)/cc-pVTZ level. These accurate results stand as a reference for benchmarking selected hybrid and double-hybrid, dispersion-corrected DFT functionals. B2PLYP-D3 and DSDPBEP86 in conjunction with a triple-ζ basis set have been confirmed as effective methodologies for structural and spectroscopic studies of medium-sized flexible biomolecules, also showing intramolecular hydrogen bonding. These best performing double-hybrid functionals have been employed to simulate IR spectra by means of vibrational perturbation theory, also considering hybrid CC/DFT schemes. The best overall agreement with experiment, with mean absolute error of 8 cm-1, has been obtained by combining CCSD(T)/(CBS+CV)MP2 harmonic wavenumbers with B2PLYP-D3/maug-cc-pVTZ anharmonic corrections. Finally, a composite scheme entirely based on CCSD(T) calculations (CCSD(T)/CBS+CV) has been employed for energetics, further confirming that serine II is the most stable conformer, also when zero-point vibrational energy corrections are included.
Collapse
Affiliation(s)
- Mingzhu Sheng
- International Centre for Quantum and Molecular Structures, Physics Department, Shanghai University, 99 Shangda Road, Shanghai, 200444, China
| | - Filippo Silvestrini
- Department of Chemistry "Giacomo Ciamician", University of Bologna, Via F. Selmi 2, 40126 Bologna, Italy
| | - Malgorzata Biczysko
- International Centre for Quantum and Molecular Structures, Physics Department, Shanghai University, 99 Shangda Road, Shanghai, 200444, China
| | - Cristina Puzzarini
- Department of Chemistry "Giacomo Ciamician", University of Bologna, Via F. Selmi 2, 40126 Bologna, Italy
| |
Collapse
|
31
|
Li Y, Zhou B, Wang K, Zhang J, Sun W, Zhang L, Guo Y. Powerful Steroid-Based Chiral Selector for High-Throughput Enantiomeric Separation of α-Amino Acids Utilizing Ion Mobility-Mass Spectrometry. Anal Chem 2021; 93:13589-13596. [PMID: 34597017 DOI: 10.1021/acs.analchem.1c02691] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Stereospecific recognition of amino acids (AAs) plays a crucial role in chiral biomarker-based diagnosis and prognosis. Separation of AA enantiomers is a long and tedious task due to the requirement of AA derivatization prior to the chromatographic or electrophoretic steps which are also time-consuming. Here, a mass-tagged chiral selector named [d0]/[d5]-estradiol-3-benzoate-17β-chloroformate ([d0]/[d5]-17β-EBC) with high reactivity and good enantiomeric resolution in regard to AAs was developed. After a quick and easy chemical derivatization step of AAs using 17β-EBC as the single chiral selector before ion mobility-mass spectrometry analysis, good enantiomer separation was achieved for 19 chiral proteinogenic AAs in a single analytical run (∼2 s). A linear calibration curve of enantiomeric excess was also established using [d0]/[d5]-17β-EBC. It was demonstrated to be capable of determining enantiomeric ratios down to 0.5% in the nanomolar range. 17β-EBC was successfully applied to investigate the absolute configuration of AAs among peptide drugs and detect trace levels of d-AAs in complex biological samples. These results indicated that [d0]/[d5]-17β-EBC may contribute to entail a valuable step forward in peptide drug quality control and discovering chiral disease biomarkers.
Collapse
Affiliation(s)
- Yuling Li
- State Key Laboratory of Organometallic Chemistry and National Center for Organic Mass Spectrometry in Shanghai, Shanghai Institute of Organic Chemistry, University of Chinese Academy of Sciences, Chinese Academy of Sciences, Shanghai 200032, China
| | - Bowen Zhou
- State Key Laboratory of Organometallic Chemistry and National Center for Organic Mass Spectrometry in Shanghai, Shanghai Institute of Organic Chemistry, University of Chinese Academy of Sciences, Chinese Academy of Sciences, Shanghai 200032, China
| | - Keke Wang
- Shimadzu Research Laboratory (Shanghai) Co. Ltd., Shanghai 201206, China
| | - Jing Zhang
- State Key Laboratory of Organometallic Chemistry and National Center for Organic Mass Spectrometry in Shanghai, Shanghai Institute of Organic Chemistry, University of Chinese Academy of Sciences, Chinese Academy of Sciences, Shanghai 200032, China
| | - Wenjian Sun
- Shimadzu Research Laboratory (Shanghai) Co. Ltd., Shanghai 201206, China
| | - Li Zhang
- State Key Laboratory of Organometallic Chemistry and National Center for Organic Mass Spectrometry in Shanghai, Shanghai Institute of Organic Chemistry, University of Chinese Academy of Sciences, Chinese Academy of Sciences, Shanghai 200032, China
| | - Yinlong Guo
- State Key Laboratory of Organometallic Chemistry and National Center for Organic Mass Spectrometry in Shanghai, Shanghai Institute of Organic Chemistry, University of Chinese Academy of Sciences, Chinese Academy of Sciences, Shanghai 200032, China
| |
Collapse
|
32
|
Danilack AD, Mulvihill CR, Klippenstein SJ, Goldsmith CF. Diastereomers and Low-Temperature Oxidation. J Phys Chem A 2021; 125:8064-8073. [PMID: 34469163 DOI: 10.1021/acs.jpca.1c05635] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Diastereomers have historically been ignored when building kinetic mechanisms for combustion. Low-temperature oxidation kinetics, which continues to gain interest in both combustion and atmospheric communities, may be affected by the inclusion of diastereomers in radical chain-branching pathways. In this work, key intermediates and transition states lacking stereochemical specification in an existing diethyl ether low-temperature oxidation mechanism were replaced with their diastereomeric counterparts. Rate coefficients for reactions involving diastereomers were computed with ab initio transition state theory master equation calculations. The presence of diastereomers increased rate coefficients by factors of 1.2-1.6 across various temperatures and pressures. Ignition delay simulations incorporating these revised rate coefficients indicate that the diastereomers enhanced the overall reactivity of the mechanism by almost 15% and increased the peak ketohydroperoxide concentration by 30% in the negative temperature coefficient region at combustion-relevant pressures. These results provide an illustrative indication of the important role of stereomeric effects in oxidation kinetics.
Collapse
Affiliation(s)
- Aaron D Danilack
- Chemical and Environmental Engineering, School of Engineering, Brown University, Providence, Rhode Island 02912, United States
| | - Clayton R Mulvihill
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Stephen J Klippenstein
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - C Franklin Goldsmith
- Chemical and Environmental Engineering, School of Engineering, Brown University, Providence, Rhode Island 02912, United States
| |
Collapse
|
33
|
Di Foggia M, Tugnoli V, Ottani S, Dettin M, Zamuner A, Sanchez-Cortes S, Cesini D, Torreggiani A. SERS Investigation on Oligopeptides Used as Biomimetic Coatings for Medical Devices. Biomolecules 2021; 11:959. [PMID: 34209793 PMCID: PMC8301923 DOI: 10.3390/biom11070959] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2021] [Revised: 06/14/2021] [Accepted: 06/24/2021] [Indexed: 11/22/2022] Open
Abstract
The surface-enhanced Raman scattering (SERS) spectra of three amphiphilic oligopeptides derived from EAK16 (AEAEAKAK)2 were examined to study systematic amino acid substitution effects on the corresponding interaction with Ag colloidal nanoparticles. Such self-assembling molecular systems, known as "molecular Lego", are of particular interest for their uses in tissue engineering and as biomimetic coatings for medical devices because they can form insoluble macroscopic membranes under physiological conditions. Spectra were collected for both native and gamma-irradiated samples. Quantum mechanical data on two of the examined oligopeptides were also obtained to clarify the assignment of the prominent significative bands observed in the spectra. In general, the peptide-nanoparticles interaction occurs through the COO- groups, with the amide bond and the aliphatic chain close to the colloid surface. After gamma irradiation, mimicking a free oxidative radical attack, the SERS spectra of the biomaterials show that COO- groups still provide the main peptide-nanoparticle interactions. However, the spatial arrangement of the peptides is different, exhibiting a systematic decrease in the distance between aliphatic chains and colloid nanoparticles.
Collapse
Affiliation(s)
- Michele Di Foggia
- Dipartimento di Scienze Biomediche e Neuromotorie, Università di Bologna, 40126 Bologna, Italy;
| | - Vitaliano Tugnoli
- Dipartimento di Scienze Biomediche e Neuromotorie, Università di Bologna, 40126 Bologna, Italy;
| | - Stefano Ottani
- Istituto per la Sintesi Organica e la Fotoreattività, Consiglio Nazionale delle Ricerche (ISOF-CNR), 40129 Bologna, Italy; (S.O.); (A.T.)
| | - Monica Dettin
- Dipartimento di Ingegneria Industriale, Università di Padova, 35131 Padova, Italy; (M.D.); (A.Z.)
| | - Annj Zamuner
- Dipartimento di Ingegneria Industriale, Università di Padova, 35131 Padova, Italy; (M.D.); (A.Z.)
| | - Santiago Sanchez-Cortes
- Instituto de Estructura de la Materia, Consejo Superior de Investigaciones Cientificas (CSIC), 28006 Madrid, Spain;
| | - Daniele Cesini
- Dipartimento CNAF, Istituto Nazionale di Fisica Nucleare (INFN-CNAF), 40127 Bologna, Italy;
| | - Armida Torreggiani
- Istituto per la Sintesi Organica e la Fotoreattività, Consiglio Nazionale delle Ricerche (ISOF-CNR), 40129 Bologna, Italy; (S.O.); (A.T.)
| |
Collapse
|
34
|
Abstract
Understanding the mechanism of hydrazine oxidation reaction by OH radical along with the rate constants of all possible pathways leads to explain the fate of hydrazine in the atmosphere. In this article, the comprehensive mechanisms and kinetics of the hydrazine plus hydroxyl radical reaction have been investigated theoretically at different temperatures and pressures. To achieve the main goals, a series of high levels of quantum chemical calculations have been widely implemented in reliable channels of the H-abstraction, SN2, and addition/elimination reactions. The energy profile of all pathways accompanied by the molecular properties of the involved stationary points has been characterized at the MP2, M06-2X, and CCSD(T)/CBS levels. To estimate accurate barrier energies of the H-abstraction channels, large numbers of the CCSD (T) calculations in conjunction with various augmented basis sets have been implemented. The direct dynamic calculations have been carried out using the validated M06-2X/maug-cc-pVTZ level, and also by the CCSD(T) (energies) + MP2 (partition functions) level. The pressure-dependent rate constants of the barrierless pathways have been investigated by the strong collision approach. Therefore, the main behaviors of the N2H4 + OH reaction have been explored according to the influences of temperature and pressure on the computed rate coefficients within the well-behaved theoretical frameworks of the TST, VTST, and RRKM theories. It has been found that the H-abstraction mechanism (to form N2H3) is dominant relative to the SN2 reaction and OH-addition to the N center of N2H4 moiety (to form H2NOH + NH2). The computed high pressure limit rate constant of the main reaction pathway, k(298.15) = 7.31 × 10–11 cm3 molecule−1 s−1, has an excellent agreement with the experimental value (k (298.15) = (6.50 ± 1.3) × 10–11 cm3 molecule−1 s−1) recommended by Vaghjiani. Also, the atmospheric lifetime of hydrazine degradation by OH radicals has been demonstrated to be 32.80 to 1161.11 h at the altitudes of 0–50 km. Finally, the disagreement in the calculated rate constants between the previous theoretical study and experimental results has been rectified.
Collapse
|
35
|
Léon I, Tasinato N, Spada L, Alonso ER, Mata S, Balbi A, Puzzarini C, Alonso JL, Barone V. Looking for the Elusive Imine Tautomer of Creatinine: Different States of Aggregation Studied by Quantum Chemistry and Molecular Spectroscopy. Chempluschem 2021; 86:1374-1386. [PMID: 34255935 PMCID: PMC8519097 DOI: 10.1002/cplu.202100224] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2021] [Revised: 06/22/2021] [Indexed: 01/06/2023]
Abstract
New spectroscopic experiments and state-of-the-art quantum-chemical computations of creatinine in different aggregation states unequivocally unveiled a significant tuning of tautomeric equilibrium by the environment: from the exclusive presence of the amine tautomer in the solid state and aqueous solution to a mixture of amine and imine tautomers in the gas phase. Quantum-chemical calculations predict the amine species as the most stable tautomer by about 30 kJ mol-1 in condensed phases. On the contrary, moving to the isolated forms, both Z and E imine isomers become more stable by about 7 kJ mol-1 . Since the imine isomers and one amine tautomer are separated by significant energy barriers, all of them should be present in the gas phase. This prediction has indeed been confirmed by high-resolution rotational spectroscopy, which provides the first experimental characterization of the elusive imine tautomer. The interpretation of the complicated hyperfine structure of the rotational spectrum, originated by three 14 N nuclei, makes it possible to use the spectral signatures as a sort of fingerprint for each individual tautomer in the complex sample.
Collapse
Affiliation(s)
- Iker Léon
- Grupo de Espectroscopia Molecular (GEM), Edificio Quifima, Laboratorios de Espectroscopia y Bioespectroscopia Parque Cientifico UVa, Universidad de Valladolid, 47011, Valladolid, Spain
| | - Nicola Tasinato
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126, Pisa, Italy
| | - Lorenzo Spada
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126, Pisa, Italy.,Dipartimento di Chimica "Giacomo Ciamician", Alma Mater Studiorum -, Università di Bologna, Via Selmi 2, 40126, Bologna, Italy
| | - Elena R Alonso
- Grupo de Espectroscopia Molecular (GEM), Edificio Quifima, Laboratorios de Espectroscopia y Bioespectroscopia Parque Cientifico UVa, Universidad de Valladolid, 47011, Valladolid, Spain
| | - Santiago Mata
- Grupo de Espectroscopia Molecular (GEM), Edificio Quifima, Laboratorios de Espectroscopia y Bioespectroscopia Parque Cientifico UVa, Universidad de Valladolid, 47011, Valladolid, Spain
| | - Alice Balbi
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126, Pisa, Italy
| | - Cristina Puzzarini
- Dipartimento di Chimica "Giacomo Ciamician", Alma Mater Studiorum -, Università di Bologna, Via Selmi 2, 40126, Bologna, Italy
| | - Jose L Alonso
- Grupo de Espectroscopia Molecular (GEM), Edificio Quifima, Laboratorios de Espectroscopia y Bioespectroscopia Parque Cientifico UVa, Universidad de Valladolid, 47011, Valladolid, Spain
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126, Pisa, Italy
| |
Collapse
|
36
|
Li W, Li J, Ning H, Shang Y, Luo SN. Multistructural Variational Reaction Kinetics of the Simplest Unsaturated Methyl Ester: H-Abstraction from Methyl Acrylate by H, OH, CH 3, and HO 2 Radicals. J Phys Chem A 2021; 125:5103-5116. [PMID: 34082530 DOI: 10.1021/acs.jpca.1c01788] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The H-abstraction reaction kinetics of methyl acrylate (MA) + H/OH/CH3/HO2 radicals have been investigated theoretically in the present work. For these reactions, the reaction energies and barrier heights are first computed using several density functionals and compared to the coupled cluster CCSD(T)-F12/jun-cc-pVTZ benchmark calculations. The M062X/maug-cc-pVTZ method shows the best performance with the smallest mean unsigned deviation (MUD) of 0.42 kcal mol-1. Combined with the electronic structure calculations using the M062X/maug-cc-pVTZ method, the multistructural canonical variational transition-state theory (MS-CVT) with small-curvature tunneling (SCT) is employed to calculate the reaction rate constants at 500-2000 K. The variational effect is between 0.56 and 1.0, the multistructural torsional anharmonicity factor ranges from 0.004 to 4.57, and the tunneling coefficient is in the range of 0.5-4.70. Notably, given the existence of reactant complexes (RCs) between reactants and transition states for the reaction systems MA + OH/HO2, we further compare the rate constants under the low-pressure limit (LPL) kinetic model, which treats the reaction as a single-step process and neglects RCs, and the pre-equilibrium model, which takes RCs into account in the reaction and treats the reaction as a two-step process. The rate constants calculated by these two models are similar within the combustion temperature range, and apparent differences occur at lower temperatures. In addition, we determine the branching ratios as a function of temperature and find that the methyl site (S3) abstractions by OH and H radicals are dominant in the low- and high-temperature ranges, respectively. Moreover, we update the kinetic model with the calculated H-abstraction rate constants to simulate the ignition delay times of MA. The simulations of the updated model are in good agreement with experimental results. The accurate reaction kinetics determined in this work are useful for the understanding and prediction of consumption branching fractions and ignition properties of the unsaturated methyl esters.
Collapse
Affiliation(s)
- Wenrui Li
- School of Chemistry and Chemical Engineering & Chongqing Key Laboratory of Theoretical and Computational Chemistry, Chongqing University, Chongqing 401331, China.,Key Laboratory of Advanced Technologies of Materials, Ministry of Education, Southwest Jiaotong University, Chengdu, Sichuan 610031, China
| | - Jun Li
- School of Chemistry and Chemical Engineering & Chongqing Key Laboratory of Theoretical and Computational Chemistry, Chongqing University, Chongqing 401331, China
| | - Hongbo Ning
- Key Laboratory of Advanced Technologies of Materials, Ministry of Education, Southwest Jiaotong University, Chengdu, Sichuan 610031, China
| | - Yanlei Shang
- Key Laboratory of Advanced Technologies of Materials, Ministry of Education, Southwest Jiaotong University, Chengdu, Sichuan 610031, China
| | - Sheng-Nian Luo
- Key Laboratory of Advanced Technologies of Materials, Ministry of Education, Southwest Jiaotong University, Chengdu, Sichuan 610031, China
| |
Collapse
|
37
|
Paoloni L, Maris A. Interplay of Rotational and Pseudorotational Motions in Flexible Cyclic Molecules. J Phys Chem A 2021; 125:4098-4113. [PMID: 33973473 PMCID: PMC8279653 DOI: 10.1021/acs.jpca.1c01472] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2021] [Revised: 04/23/2021] [Indexed: 11/29/2022]
Abstract
Solutions to the time-independent nuclear Schrödinger equation associated with the pseudorotational motion of three flexible cyclic molecules are presented and discussed. Structural relaxations related to the pseudorotational motion are described as functions of a pseudorotation angle ϕ which is formulated according to the definition of ring-puckering coordinates originally proposed by Cremer and Pople ( J. Am. Chem. Soc. 1975, 97 (6), 1354-1358). In order to take into account the interplay between pseudorotational and rotational motions, the rovibrational Hamiltonian matrices are formulated for the rotational quantum numbers J = 0 and J = 1. The rovibrational Hamiltonian matrices are constructed and diagonalized using a Python program developed by the authors. Suitable algorithms for (i) the construction of one-dimensional cuts of potential energy surfaces along the pseudorotation angle ϕ and (ii) the assignment of the vibrorotational wave functions (which are needed for the automatic calculation of rotational transition energies J = 0 → J = 1) are described and discussed.
Collapse
Affiliation(s)
- Lorenzo Paoloni
- Dipartimento
di Fisica e Astronomia, Università
di Padova, via Marzolo 8, I-35131 Padova, Italy
| | - Assimo Maris
- Dipartimento
di Chimica G. Ciamician, Università
di Bologna, via Selmi 2, I-40126 Bologna, Italy
| |
Collapse
|
38
|
Varga Z, Liu Y, Li J, Paukku Y, Guo H, Truhlar DG. Potential energy surfaces for high-energy N + O 2 collisions. J Chem Phys 2021; 154:084304. [PMID: 33639765 DOI: 10.1063/5.0039771] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Potential energy surfaces for high-energy collisions between an oxygen molecule and a nitrogen atom are useful for modeling chemical dynamics in shock waves. In the present work, we present doublet, quartet, and sextet potential energy surfaces that are suitable for studying collisions of O2(3Σg -) with N(4S) in the electronically adiabatic approximation. Two sets of surfaces are developed, one using neural networks (NNs) with permutationally invariant polynomials (PIPs) and one with the least-squares many-body (MB) method, where a two-body part is an accurate diatomic potential and the three-body part is expressed with connected PIPs in mixed-exponential-Gaussian bond order variables (MEGs). We find, using the same dataset for both fits, that the fitting performance of the PIP-NN method is significantly better than that of the MB-PIP-MEG method, even though the MB-PIP-MEG fit uses a higher-order PIP than those used in previous MB-PIP-MEG fits of related systems (such as N4 and N2O2). However, the evaluation of the PIP-NN fit in trajectory calculations requires about 5 times more computer time than is required for the MB-PIP-MEG fit.
Collapse
Affiliation(s)
- Zoltan Varga
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| | - Yang Liu
- School of Chemistry and Chemical Engineering & Chongqing Key Laboratory of Theoretical and Computational Chemistry, Chongqing University, Chongqing 401331, China
| | - Jun Li
- School of Chemistry and Chemical Engineering & Chongqing Key Laboratory of Theoretical and Computational Chemistry, Chongqing University, Chongqing 401331, China
| | - Yuliya Paukku
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| |
Collapse
|
39
|
Barclay AJ, Pietropolli Charmet A, McKellar ARW, Moazzen-Ahmadi N. Exploring the next step in micro-solvation of CO in water: Infrared spectra and structural calculations of (H2O)4–CO and (D2O)4–CO. J Chem Phys 2021; 154:044310. [DOI: 10.1063/5.0038188] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023] Open
Affiliation(s)
- A. J. Barclay
- Department of Physics and Astronomy, University of Calgary, 2500 University Drive North West, Calgary, Alberta T2N 1N4, Canada
| | - A. Pietropolli Charmet
- Dipartimento di Scienze Molecolari e Nanosistemi, Università Ca’ Foscari Venezia, Via Torino 155, I-30172 Mestre, Venezia, Italy
| | - A. R. W. McKellar
- National Research Council of Canada, Ottawa, Ontario K1A 0R6, Canada
| | - N. Moazzen-Ahmadi
- Department of Physics and Astronomy, University of Calgary, 2500 University Drive North West, Calgary, Alberta T2N 1N4, Canada
| |
Collapse
|
40
|
Zheng P, Ji C, Ying F, Su P, Wu W. A Valence-Bond-Based Multiconfigurational Density Functional Theory: The λ-DFVB Method Revisited. Molecules 2021; 26:521. [PMID: 33498268 PMCID: PMC7863953 DOI: 10.3390/molecules26030521] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2021] [Revised: 01/17/2021] [Accepted: 01/18/2021] [Indexed: 11/29/2022] Open
Abstract
A recently developed valence-bond-based multireference density functional theory, named λ-DFVB, is revisited in this paper. λ-DFVB remedies the double-counting error of electron correlation by decomposing the electron-electron interactions into the wave function term and density functional term with a variable parameter λ. The λ value is defined as a function of the free valence index in our previous scheme, denoted as λ-DFVB(K) in this paper. Here we revisit the λ-DFVB method and present a new scheme based on natural orbital occupation numbers (NOONs) for parameter λ, named λ-DFVB(IS), to simplify the process of λ-DFVB calculation. In λ-DFVB(IS), the parameter λ is defined as a function of NOONs, which are straightforwardly determined from the many-electron wave function of the molecule. Furthermore, λ-DFVB(IS) does not involve further self-consistent field calculation after performing the valence bond self-consistent field (VBSCF) calculation, and thus, the computational effort in λ-DFVB(IS) is approximately the same as the VBSCF method, greatly reduced from λ-DFVB(K). The performance of λ-DFVB(IS) was investigated on a broader range of molecular properties, including equilibrium bond lengths and dissociation energies, atomization energies, atomic excitation energies, and chemical reaction barriers. The computational results show that λ-DFVB(IS) is more robust without losing accuracy and comparable in accuracy to high-level multireference wave function methods, such as CASPT2.
Collapse
Affiliation(s)
| | | | | | - Peifeng Su
- Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, The State Key Laboratory of Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China; (P.Z.); (C.J.); (F.Y.)
| | - Wei Wu
- Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, The State Key Laboratory of Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China; (P.Z.); (C.J.); (F.Y.)
| |
Collapse
|
41
|
Barone V, Ceselin G, Fusè M, Tasinato N. Accuracy Meets Interpretability for Computational Spectroscopy by Means of Hybrid and Double-Hybrid Functionals. Front Chem 2020; 8:584203. [PMID: 33195078 PMCID: PMC7645164 DOI: 10.3389/fchem.2020.584203] [Citation(s) in RCA: 41] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2020] [Accepted: 08/17/2020] [Indexed: 11/13/2022] Open
Abstract
Accuracy and interpretability are often seen as the devil and holy grail in computational spectroscopy and their reconciliation remains a primary research goal. In the last few decades, density functional theory has revolutionized the situation, paving the way to reliable yet effective models for medium size molecules, which could also be profitably used by non-specialists. In this contribution we will compare the results of some widely used hybrid and double hybrid functionals with the aim of defining the most suitable recipe for all the spectroscopic parameters of interest in rotational and vibrational spectroscopy, going beyond the rigid rotor/harmonic oscillator model. We will show that last-generation hybrid and double hybrid functionals in conjunction with partially augmented double- and triple-zeta basis sets can offer, in the framework of second order vibrational perturbation theory, a general, robust, and user-friendly tool with unprecedented accuracy for medium-size semi-rigid molecules.
Collapse
Affiliation(s)
- Vincenzo Barone
- SMART Laboratory, Scuola Normale Superiore di Pisa, Pisa, Italy
| | - Giorgia Ceselin
- SMART Laboratory, Scuola Normale Superiore di Pisa, Pisa, Italy
| | - Marco Fusè
- SMART Laboratory, Scuola Normale Superiore di Pisa, Pisa, Italy
| | - Nicola Tasinato
- SMART Laboratory, Scuola Normale Superiore di Pisa, Pisa, Italy
| |
Collapse
|
42
|
Alessandrini S, Dell'Isola V, Spada L, Barone V, Puzzarini C. A computational journey in the CH 2O 2S land: an accurate rotational and ro-vibrational analysis of the sulfene molecule and the O,S- and O,O-monothiocarbonic acids. Mol Phys 2020. [DOI: 10.1080/00268976.2020.1766707] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
Affiliation(s)
- S. Alessandrini
- Scuola Normale Superiore, Pisa, Italy
- Dipartimento di Chimica “Giacomo Ciamician”, Università di Bologna, Bologna, Italy
| | - V. Dell'Isola
- Dipartimento di Chimica “Giacomo Ciamician”, Università di Bologna, Bologna, Italy
| | - L. Spada
- Scuola Normale Superiore, Pisa, Italy
- Dipartimento di Chimica “Giacomo Ciamician”, Università di Bologna, Bologna, Italy
| | - V. Barone
- Scuola Normale Superiore, Pisa, Italy
| | - C. Puzzarini
- Dipartimento di Chimica “Giacomo Ciamician”, Università di Bologna, Bologna, Italy
| |
Collapse
|
43
|
Qin Z, Hou GL, Yang Z, Valiev M, Wang XB. Distonic radical anion species in cysteine oxidation processes. Phys Chem Chem Phys 2020; 22:17554-17558. [PMID: 32716467 DOI: 10.1039/d0cp02165b] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Abstract
Oxidation of cysteine residues constitutes an important regulatory mechanism in the function of biological systems. Much of this behavior is controlled by the specific chemical properties of the thiol side-chain group, where reactions with reactive oxygen species take place. Herein, we investigated the entire cysteine oxidation cycle Cys-SH → Cys-SOnH (n = 1, 2, and 3) using cryogenic negative ion photoelectron spectroscopy and quantum-chemical calculations. The conventional view of the first reversible oxidation step (n = 1) is associated with sulfenate species. Yet our results indicate that an alternative option exists in the form of a novel distonic radical anion, ˙OS-CH2CH(NH2)-COO-, with an unpaired electron on the thiol group and excess negative charge on the carboxylate group. Higher order oxidation states (n = 2 and 3) are thought to be associated with irreversible oxidative damage, and our results show that excess negative charge in those cases migrates to the -SOn- group. Furthermore, these species are stable towards 1e oxidation, as opposed to the n = 1 case that undergoes intra-molecular proton transfer. The molecular level insights reported in this work provide direct spectroscopic evidence of the unique chemical versatility of Cys-sulfenic acid (Cys-SOH) in post-translational modifications of protein systems.
Collapse
Affiliation(s)
- Zhengbo Qin
- Anhui Province Key Laboratory of Optoelectric Materials Science and Technology, Anhui Normal University, Wuhu, 241000, China.
| | | | | | | | | |
Collapse
|
44
|
Chan B. Fullerene Thermochemical Stability: Accurate Heats of Formation for Small Fullerenes, the Importance of Structural Deformation on Reactivity, and the Special Stability of C 60. J Phys Chem A 2020; 124:6688-6698. [PMID: 32786665 DOI: 10.1021/acs.jpca.0c04732] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
We have used quantum chemistry computations, in conjunction with isodesmic-type reactions, to obtain accurate heats of formation (HoFs) for the small fullerenes C20 (2358.2 ± 8.0 kJ mol-1), C24 (2566.2 ± 7.6), and the lowest-energy isomers of C32 (2461.1 ± 15.4), C42 (2629.0 ± 20.5), and C54 (2686.2 ± 25.3). As part of this endeavor, we have also obtained accurate HoFs for several medium-sized molecules, namely 216.6 ± 1.4 for fulvene, 375.5 ± 1.5 for pentalene, 670.8 ± 2.9 for acepentalene, and 262.7 ± 2.5 for acenaphthylene. We combine the energies of the small fullerenes and previously obtained energies for larger fullerenes (from C60 to C6000) into a full picture of fullerene thermochemical stability. In general, the per-carbon energies can be reasonably approximated by the "R+D" model that we have previously developed [Chan et al. J. Chem. Theory Comput. 2019, 15, 1255-1264], which takes into account Resonance and structural Deformation factors. In a case study on C54, we find that most of the high-deformation-energy atoms correspond to the sites of the C-Cl bond in the experimentally captured C54Cl8. In another case study, we find that C60 has the lowest value for the maximum local-deformation energy when compared with similar-sized fullerenes, which is consistent with its "special stability". These results are indicative of structural deformation playing an important role in the reactivity of fullerenes.
Collapse
Affiliation(s)
- Bun Chan
- Graduate School of Engineering, Nagasaki University, Bunkyo 1-14, Nagasaki-shi, Nagasaki 852-8521, Japan
| |
Collapse
|
45
|
Paoloni L, Fusè M, Baiardi A, Barone V. Interplay of Stereoelectronic and Vibrational Modulation Effects in Tuning the UPS Spectra of Unsaturated Hydrocarbon Cage Compounds. J Chem Theory Comput 2020; 16:5218-5226. [PMID: 32667793 PMCID: PMC8009518 DOI: 10.1021/acs.jctc.0c00645] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2020] [Indexed: 11/28/2022]
Abstract
The UPS spectra of six hydrocarbon cage compounds have been investigated by a Green-function approach in conjunction with a full harmonic treatment of vibrational modulation effects. The remarkable agreement with experimental results points out the reliability of the proposed computational approach and the strong interplay of stereoelectronic and vibrational effects in tuning the overall spectra.
Collapse
Affiliation(s)
- Lorenzo Paoloni
- Scuola
Normale Superiore, Piazza dei Cavalieri 7, 56125 Pisa, Italy
| | - Marco Fusè
- Scuola
Normale Superiore, Piazza dei Cavalieri 7, 56125 Pisa, Italy
| | - Alberto Baiardi
- Lab.
für Physikalische Chemie, Vladimir-Prelog-Weg 1-5/10, 8093 Zürich, Switzerland
| | - Vincenzo Barone
- Scuola
Normale Superiore, Piazza dei Cavalieri 7, 56125 Pisa, Italy
| |
Collapse
|
46
|
Ye H, Mendolicchio M, Kruse H, Puzzarini C, Biczysko M, Barone V. The challenging equilibrium structure of HSSH: Another success of the rotational spectroscopy / quantum chemistry synergism. J Mol Struct 2020. [DOI: 10.1016/j.molstruc.2020.127933] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
47
|
Baiano C, Lupi J, Tasinato N, Puzzarini C, Barone V. The Role of State-of-the-Art Quantum-Chemical Calculations in Astrochemistry: Formation Route and Spectroscopy of Ethanimine as a Paradigmatic Case. Molecules 2020; 25:E2873. [PMID: 32580443 PMCID: PMC7357107 DOI: 10.3390/molecules25122873] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2020] [Revised: 06/17/2020] [Accepted: 06/18/2020] [Indexed: 12/13/2022] Open
Abstract
The gas-phase formation and spectroscopic characteristics of ethanimine have been re-investigated as a paradigmatic case illustrating the accuracy of state-of-the-art quantum-chemical (QC) methodologies in the field of astrochemistry. According to our computations, the reaction between the amidogen, NH, and ethyl, C2H5, radicals is very fast, close to the gas-kinetics limit. Although the main reaction channel under conditions typical of the interstellar medium leads to methanimine and the methyl radical, the predicted amount of the two E,Z stereoisomers of ethanimine is around 10%. State-of-the-art QC and kinetic models lead to a [E-CH3CHNH]/[Z-CH3CHNH] ratio of ca. 1.4, slightly higher than the previous computations, but still far from the value determined from astronomical observations (ca. 3). An accurate computational characterization of the molecular structure, energetics, and spectroscopic properties of the E and Z isomers of ethanimine combined with millimeter-wave measurements up to 300 GHz, allows for predicting the rotational spectrum of both isomers up to 500 GHz, thus opening the way toward new astronomical observations.
Collapse
Affiliation(s)
- Carmen Baiano
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy; (C.B.); (J.L.); (N.T.)
| | - Jacopo Lupi
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy; (C.B.); (J.L.); (N.T.)
| | - Nicola Tasinato
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy; (C.B.); (J.L.); (N.T.)
| | - Cristina Puzzarini
- 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; (C.B.); (J.L.); (N.T.)
| |
Collapse
|
48
|
Puzzarini C, Spada L, Alessandrini S, Barone V. The challenge of non-covalent interactions: theory meets experiment for reconciling accuracy and interpretation. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2020; 32:343002. [PMID: 32203942 DOI: 10.1088/1361-648x/ab8253] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/27/2019] [Accepted: 03/23/2020] [Indexed: 06/10/2023]
Abstract
In the past decade, many gas-phase spectroscopic investigations have focused on the understanding of the nature of weak interactions in model systems. Despite the fact that non-covalent interactions play a key role in several biological and technological processes, their characterization and interpretation are still far from being satisfactory. In this connection, integrated experimental and computational investigations can play an invaluable role. Indeed, a number of different issues relevant to unraveling the properties of bulk or solvated systems can be addressed from experimental investigations on molecular complexes. Focusing on the interaction of biological model systems with solvent molecules (e.g., water), since the hydration of the biomolecules controls their structure and mechanism of action, the study of the molecular properties of hydrated systems containing a limited number of water molecules (microsolvation) is the basis for understanding the solvation process and how structure and reactivity vary from gas phase to solution. Although hydrogen bonding is probably the most widespread interaction in nature, other emerging classes, such as halogen, chalcogen and pnicogen interactions, have attracted much attention because of the role they play in different fields. Their understanding requires, first of all, the characterization of the directionality, strength, and nature of such interactions as well as a comprehensive analysis of their competition with other non-covalent bonds. In this review, it is shown how state-of-the-art quantum-chemical computations combined with rotational spectroscopy allow for fully characterizing intermolecular interactions taking place in molecular complexes from both structural and energetic points of view. The transition from bi-molecular complex to microsolvation and then to condensed phase is shortly addressed.
Collapse
Affiliation(s)
- Cristina Puzzarini
- Dipartimento di Chimica 'Giacomo Ciamician', Via F. Selmi 2, I-40126 Bologna, Italy
| | - Lorenzo Spada
- Dipartimento di Chimica 'Giacomo Ciamician', Via F. Selmi 2, I-40126 Bologna, Italy
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Silvia Alessandrini
- Dipartimento di Chimica 'Giacomo Ciamician', Via F. Selmi 2, I-40126 Bologna, Italy
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| |
Collapse
|
49
|
Toward accurate prediction of amino acid derivatives structure and energetics from DFT: glycine conformers and their interconversions. J Mol Model 2020; 26:129. [DOI: 10.1007/s00894-020-4342-7] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2019] [Accepted: 03/04/2020] [Indexed: 12/30/2022]
|
50
|
DFT meets the segmented polarization consistent basis sets: Performances in the computation of molecular structures, rotational and vibrational spectroscopic properties. J Mol Struct 2020. [DOI: 10.1016/j.molstruc.2020.127886] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|