51
|
Accurate estimation of singlet-triplet gap of strongly correlated systems by CCSD(T) method using improved orbitals. COMPUT THEOR CHEM 2021. [DOI: 10.1016/j.comptc.2021.113326] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
52
|
Ciborowski SM, Liu G, Blankenhorn M, Harris RM, Marshall MA, Zhu Z, Bowen KH, Peterson KA. The electron affinity of the uranium atom. J Chem Phys 2021; 154:224307. [PMID: 34241193 DOI: 10.1063/5.0046315] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The results of a combined experimental and computational study of the uranium atom are presented with the aim of determining its electron affinity. Experimentally, the electron affinity of uranium was measured via negative ion photoelectron spectroscopy of the uranium atomic anion, U-. Computationally, the electron affinities of both thorium and uranium were calculated by conducting relativistic coupled-cluster and multi-reference configuration interaction calculations. The experimentally determined value of the electron affinity of the uranium atom was determined to be 0.309 ± 0.025 eV. The computationally predicted electron affinity of uranium based on composite coupled cluster calculations and full four-component spin-orbit coupling was found to be 0.232 eV. Predominately due to a better convergence of the coupled cluster sequence for Th and Th-, the final calculated electron affinity of Th, 0.565 eV, was in much better agreement with the accurate experimental value of 0.608 eV. In both cases, the ground state of the anion corresponds to electron attachment to the 6d orbital.
Collapse
Affiliation(s)
- Sandra M Ciborowski
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Gaoxiang Liu
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Moritz Blankenhorn
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Rachel M Harris
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Mary A Marshall
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Zhaoguo Zhu
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Kit H Bowen
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Kirk A Peterson
- Department of Chemistry, Washington State University, Pullman, Washington 99162, USA
| |
Collapse
|
53
|
Chilkuri VG, Neese F. Comparison of Many-Particle Representations for Selected Configuration Interaction: II. Numerical Benchmark Calculations. J Chem Theory Comput 2021; 17:2868-2885. [PMID: 33886300 PMCID: PMC8279407 DOI: 10.1021/acs.jctc.1c00081] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
The present work
is the second part in our three-part series on
the comparison of many-particle representations for the selected configuration
interaction (CI) method. In this work, we present benchmark calculations
based on our selected CI program called the iterative configuration
expansion (ICE) that is inspired by the CIPSI method of Malrieu and
co-workers (Malrieu1973, 58, ( (12), ), 5745−5759). We describe the main parameters
that enter in this algorithm and perform benchmark calculations on
a set of 21 small molecules and compare ground state energies with
full configuration interaction (FCI) results (FCI21 test set). The
focus is the comparison of the performance of three different types
of many-particle basis functions (MPBFs): (1) individual Slater determinants
(DETS), (2) individual spin-adapted configuration state functions
(CSFs), and (3) all CSFs of a given total spin that can be generated
from spatial configurations (CFGs). An analysis of the cost of the
calculation in terms of the number of wavefunction parameters and
the energy error is evaluated for the DET-, CFG-, and CSF-based ICE.
The main differences for the three many-particle basis representations
show up in the number of wavefunction parameters and the rate of convergence
toward the FCI limit with the thresholds of the ICE. Next, we analyze
the best way to extrapolate the ICE energies toward the FCI results
as a function of the thresholds. The efficiency of the extrapolation
is investigated relative to the FCI21 test set as well as near FCI
calculations on three moderately sized hydrocarbon molecules CH4, C2H4, and C4H6. Finally, we comment on the size-inconsistency error for the three
many-particle representations and compare it with the error in the
total energy. The implication for selected CI implementations with
any of the three many-particle representations is discussed.
Collapse
Affiliation(s)
- Vijay Gopal Chilkuri
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, Mülheim an der Ruhr 45470, Germany
| | - Frank Neese
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, Mülheim an der Ruhr 45470, Germany
| |
Collapse
|
54
|
Stoychev GL, Auer AA, Gauss J, Neese F. DLPNO-MP2 second derivatives for the computation of polarizabilities and NMR shieldings. J Chem Phys 2021; 154:164110. [PMID: 33940835 DOI: 10.1063/5.0047125] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a derivation and efficient implementation of the formally complete analytic second derivatives for the domain-based local pair natural orbital second order Møller-Plesset perturbation theory (MP2) method, applicable to electric or magnetic field-response properties but not yet to harmonic frequencies. We also discuss the occurrence and avoidance of numerical instability issues related to singular linear equation systems and near linear dependences in the projected atomic orbital domains. A series of benchmark calculations on medium-sized systems is performed to assess the effect of the local approximation on calculated nuclear magnetic resonance shieldings and the static dipole polarizabilities. Relative deviations from the resolution of the identity-based MP2 (RI-MP2) reference for both properties are below 0.5% with the default truncation thresholds. For large systems, our implementation achieves quadratic effective scaling, is more efficient than RI-MP2 starting at 280 correlated electrons, and is never more than 5-20 times slower than the equivalent Hartree-Fock property calculation. The largest calculation performed here was on the vancomycin molecule with 176 atoms, 542 correlated electrons, and 4700 basis functions and took 3.3 days on 12 central processing unit cores.
Collapse
Affiliation(s)
- Georgi L Stoychev
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Alexander A Auer
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Jürgen Gauss
- Department Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, 55128 Mainz, Germany
| | - Frank Neese
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
55
|
Matthews DA. A critical analysis of least-squares tensor hypercontraction applied to MP3. J Chem Phys 2021; 154:134102. [PMID: 33832252 DOI: 10.1063/5.0038764] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The least-squares tensor hypercontraction (LS-THC) approach is a promising method of reducing the high polynomial scaling of wavefunction methods, for example, those based on many-body perturbation theory or coupled cluster. Here, we focus on LS-THC-MP3 and identify four variants with differing errors and efficiency characteristics. The performance of LS-THC-MP3 is analyzed for regular test systems with up to 40 first-row atoms. We also analyze the size-extensivity/size-consistency and grid- and basis set dependence of LS-THC-MP3. Overall, the errors observed are favorably small in comparison with standard density fitting, and a more streamlined method of generating grids via pruning is suggested. A practical crossover (the point at which LS-THC-MP3 is cheaper than the canonical method) is achieved around 240 correlated electrons. Despite several drawbacks of LS-THC that have been identified: an initial non-linearity of error when increasing system size, poor description of angular correlation, and a potentially large increase in error with the basis set size, the results show that LS-THC has significant potential for practical application to MP3 and other wavefunction methods.
Collapse
|
56
|
Szabó PB, Csóka J, Kállay M, Nagy PR. Linear-Scaling Open-Shell MP2 Approach: Algorithm, Benchmarks, and Large-Scale Applications. J Chem Theory Comput 2021; 17:2886-2905. [PMID: 33819030 PMCID: PMC8154337 DOI: 10.1021/acs.jctc.1c00093] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
A linear-scaling
local second-order Møller–Plesset
(MP2) method is presented for high-spin open-shell molecules based
on restricted open-shell (RO) reference functions. The open-shell
local MP2 (LMP2) approach inherits the iteration- and redundancy-free
formulation and the completely integral-direct, OpenMP-parallel, and
memory and disk use economic algorithms of our closed-shell LMP2 implementation.
By utilizing restricted local molecular orbitals for the demanding
integral transformation step and by introducing a novel long-range
spin-polarization approximation, the computational cost of RO-LMP2
approaches that of closed-shell LMP2. Extensive benchmarks were performed
for reactions of radicals, ionization potentials, as well as spin-state
splittings of carbenes and transition-metal complexes. Compared to
the conventional MP2 reference for systems of up to 175 atoms, local
errors of at most 0.1 kcal/mol were found, which are well below the
intrinsic accuracy of MP2. RO-LMP2 computations are presented for
challenging protein models of up to 601 atoms and 11 000 basis
functions, which involve either spin states of a complexed iron ion
or a highly delocalized singly occupied orbital. The corresponding
runtimes of 9–15 h obtained with a single, many-core CPU demonstrate
that MP2, as well as spin-scaled MP2 and double-hybrid density functional
methods, become widely accessible for open-shell systems of unprecedented
size and complexity.
Collapse
Affiliation(s)
- P Bernát Szabó
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - József Csóka
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Mihály Kállay
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Péter R Nagy
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| |
Collapse
|
57
|
Véril M, Scemama A, Caffarel M, Lipparini F, Boggio‐Pasqua M, Jacquemin D, Loos P. QUESTDB
: A database of highly accurate excitation energies for the electronic structure community. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2021. [DOI: 10.1002/wcms.1517] [Citation(s) in RCA: 31] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
- Mickaël Véril
- Laboratoire de Chimie et Physique Quantiques Université de Toulouse, CNRS, UPS Toulouse France
| | - Anthony Scemama
- Laboratoire de Chimie et Physique Quantiques Université de Toulouse, CNRS, UPS Toulouse France
| | - Michel Caffarel
- Laboratoire de Chimie et Physique Quantiques Université de Toulouse, CNRS, UPS Toulouse France
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale University of Pisa Pisa Italy
| | - Martial Boggio‐Pasqua
- Laboratoire de Chimie et Physique Quantiques Université de Toulouse, CNRS, UPS Toulouse France
| | | | - Pierre‐François Loos
- Laboratoire de Chimie et Physique Quantiques Université de Toulouse, CNRS, UPS Toulouse France
| |
Collapse
|
58
|
Ni Z, Guo Y, Neese F, Li W, Li S. Cluster-in-Molecule Local Correlation Method with an Accurate Distant Pair Correction for Large Systems. J Chem Theory Comput 2021; 17:756-766. [PMID: 33410327 DOI: 10.1021/acs.jctc.0c00831] [Citation(s) in RCA: 31] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The cluster-in-molecule (CIM) local correlation approach with an accurate distant pair correlation energy correction is presented. For large systems, the inclusion of distant pair correlation energies is essential for the accurate prediction of absolute correlation energies and relative energies. Here, we propose a simple and efficient scheme for evaluating the distant pair correlation energy correction for the CIM approaches. The corrections can be readily extracted from electron correlation calculations of clusters with almost no additional effort. Benchmark calculations show that the improved CIM approach can recover more than 99.94% of the correlation energy calculated by the parent method. By combining the CIM approach with the domain-based local pair natural orbital (DLPNO) local correlation approach, we have provided accurate binding energies at the CIM-DLPNO-CCSD(T) level for a test set consisting of eight weakly bound complexes ranging in size from 200 to 1027 atoms. With these results as the reference data, the accuracy and applicability of other electron correlation methods and a few density functional methods for large systems have been assessed.
Collapse
Affiliation(s)
- Zhigang Ni
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, China.,College of Materials, Chemistry and Chemical Engineering, Hangzhou Normal University, Hangzhou 311121, China
| | - Yang Guo
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao 266237, China
| | - Frank Neese
- Max Planck Institut für Kohlenforschung, Kaiser-Wilhelm Platz 1, D-45470 Mülheim an der Ruhr, Germany.,FAccTs GmbH, Rolandstr. 67, 50677 Köln, Germany
| | - Wei Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, China
| | - Shuhua Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, China
| |
Collapse
|
59
|
Franke PR, Stanton JF, Douberly GE. How to VPT2: Accurate and Intuitive Simulations of CH Stretching Infrared Spectra Using VPT2+K with Large Effective Hamiltonian Resonance Treatments. J Phys Chem A 2021; 125:1301-1324. [PMID: 33506678 DOI: 10.1021/acs.jpca.0c09526] [Citation(s) in RCA: 65] [Impact Index Per Article: 21.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
This article primarily discusses the utility of vibrational perturbation theory for the prediction of X-H stretching vibrations with particular focus on the specific variant, second-order vibrational perturbation theory with resonances (VPT2+K). It is written as a tutorial, reprinting most important formulas and providing numerous simple examples. It discusses the philosophy and practical considerations behind vibrational simulations with VPT2+K, including but not limited to computational method selection, cost-saving approximations, approaches to evaluating intensity, resonance identification, and effective Hamiltonian structure. Particular attention is given to resonance treatments, beginning with simple Fermi dyads and gradually progressing to arbitrarily large polyads that describe both Fermi and Darling-Dennison resonances. VPT2+K combined with large effective Hamiltonians is shown to be a reliable framework for modeling the complicated CH stretching spectra of alkenes. An error is also corrected in the published analytic formula for the VPT2 transition moment between the vibrational ground state and triply excited states.
Collapse
Affiliation(s)
- Peter R Franke
- Department of Chemistry, University of Florida, Gainesville, Florida 32611, United States.,Department of Chemistry, University of Georgia, Athens, Georgia 30602, United States
| | - John F Stanton
- Department of Chemistry, University of Florida, Gainesville, Florida 32611, United States
| | - Gary E Douberly
- Department of Chemistry, University of Georgia, Athens, Georgia 30602, United States
| |
Collapse
|
60
|
Ma Q, Werner HJ. Scalable Electron Correlation Methods. 8. Explicitly Correlated Open-Shell Coupled-Cluster with Pair Natural Orbitals PNO-RCCSD(T)-F12 and PNO-UCCSD(T)-F12. J Chem Theory Comput 2021; 17:902-926. [PMID: 33405921 DOI: 10.1021/acs.jctc.0c01129] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Abstract
We present explicitly correlated open-shell pair natural orbital local coupled-cluster methods, PNO-RCCSD(T)-F12 and PNO-UCCSD(T)-F12. The methods are extensions of our previously reported PNO-R/UCCSD methods (J. Chem. Theory Comput., 2020, 16, 3135-3151, https://pubs.acs.org/doi/10.1021/acs.jctc.0c00192) with additions of explicit correlation and perturbative triples corrections. The explicit correlation treatment follows the spin-orbital CCSD-F12b theory using Ansatz 3*A, which is found to yield comparable or better basis set convergence than the more rigorous Ansatz 3C in computed ionization potentials and reaction energies using double- to quaduple-ζ basis sets. The perturbative triples correction is adapted from the spin-orbital (T) theory to use triples natural orbitals (TNOs). To address the coupling due to off-diagonal Fock matrix elements, the local triples amplitudes are iteratively solved using small domains of TNOs, and a semicanonical (T0) domain correction with larger domains is applied to reduce the domain errors. The performance of the methods is demonstrated through benchmark calculations on ionization potentials, radical stabilization energies, reaction energies of fragmentations and rearrangements in radical cations, and spin-state energy differences of iron complexes. For a few test sets where canonical calculations are feasible, PNO-RCCSD(T)-F12 results agree with the canonical ones to within 0.4 kcal mol-1, and this maximum error is reduced to below 0.2 kcal mol-1 when large local domains are used. For larger systems, results using different thresholds for the local approximations are compared to demonstrate that 1 kcal mol-1 level of accuracy can be achieved using our default settings. For a couple of difficult cases, it is demonstrated that the errors from individual approximations are only a fraction of 1 kcal mol-1, and the overall accuracy of the method does not rely on error compensations. In contrast to canonical calculations, the use of spin-orbitals does not lead to a significant increase of computational time and memory usage in the most expensive steps of PNO-R/UCCSD(T)-F12 calculations. The only exception is the iterative solution of the (T) amplitudes, which can be avoided without significant errors by using a perturbative treatment of the off-diagonal coupling, known as (T1) approximation. For most systems, even the semicanonical approximation (T0) leads only to small errors in relative energies. Our program is well parallelized and capable of computing accurate correlation energies for molecules with 100-200 atoms using augmented triple-ζ basis sets in less than a day of elapsed time on a small computer cluster.
Collapse
Affiliation(s)
- Qianli Ma
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - Hans-Joachim Werner
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| |
Collapse
|
61
|
Adjieufack AI, Liégeois V, Mbouombouo Ndassa I, Champagne B. Topological investigation of the reaction mechanism of glycerol carbonate decomposition by bond evolution theory. RSC Adv 2021; 11:10083-10093. [PMID: 35423535 PMCID: PMC8695523 DOI: 10.1039/d0ra09755a] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2020] [Accepted: 02/26/2021] [Indexed: 11/30/2022] Open
Abstract
The reaction mechanisms of the decomposition of glycerol carbonate have been investigated at the density functional theory level within the bond evolution theory. The four reaction pathways yield to 3-hydroxypropanal (TS1), glycidol (TS2a and TS2b), and 4-methylene-1,3-dioxolan-2-one (TS3). The study reveals non-concerted processes with the same number (four) of structural stability domains for each reaction pathway. For the two decarboxylation mechanisms, the two first steps are similar. They correspond to the cleavage of two single CO bonds to the detriment of the increased population of the lone pairs of two O atoms. These are followed, along TS1, by the transformation of a CO single bond into a double bond together with a proton transfer to create a CH bond. For TS2a and TS2b, the last step is a cyclization by CO bond formation. For the TS3 pathway, the first stage consists in the cleavage of a CH bond and the transfer of its electron population to both a proton and a C atom, the second step corresponds to the formation of an OH bond, and the last one describes the formation of a CC double bond. Moreover, the analysis of the energies, enthalpies, and free enthalpies of reaction and of activation leads to the conclusion that 3-hydroxypropanal is both the thermodynamic and kinetic product, independent of the method of calculation. ELFs of glycerol carbonate and of its kinetic and thermodynamic decomposition product, 3-hydroxypropanal (+CO2).![]()
Collapse
Affiliation(s)
- Abel Idrice Adjieufack
- Physical and Theoretical Chemistry Laboratory
- Cameroon
- Computational Chemistry Laboratory
- High Teacher Training College
- Cameroon
| | - Vincent Liégeois
- Laboratory of Theoretical Chemistry and Namur Institute of Structured Matter (NISM)
- University of Namur
- B-5000 Namur
- Belgium
| | | | - Benoît Champagne
- Laboratory of Theoretical Chemistry and Namur Institute of Structured Matter (NISM)
- University of Namur
- B-5000 Namur
- Belgium
| |
Collapse
|
62
|
Tew DP. Principal domains in F12 explicitly correlated theory. ADVANCES IN QUANTUM CHEMISTRY 2021. [DOI: 10.1016/bs.aiq.2021.06.001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
63
|
Calvin JA, Peng C, Rishi V, Kumar A, Valeev EF. Many-Body Quantum Chemistry on Massively Parallel Computers. Chem Rev 2020; 121:1203-1231. [DOI: 10.1021/acs.chemrev.0c00006] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Affiliation(s)
- Justus A. Calvin
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Chong Peng
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Varun Rishi
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Ashutosh Kumar
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Edward F. Valeev
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| |
Collapse
|
64
|
Tajti A, Szalay PG, Kochanov R, Tyuterev VG. Diagonal Born-Oppenheimer corrections to the ground electronic state potential energy surfaces of ozone: improvement of ab initio vibrational band centers for the 16O 3, 17O 3 and 18O 3 isotopologues. Phys Chem Chem Phys 2020; 22:24257-24269. [PMID: 33089270 DOI: 10.1039/d0cp02457k] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Mass-dependent diagonal Born-Oppenheimer corrections (DBOCs) to the ab initio electronic ground state potential energy surface for the main 16O3 isotopologue and for homogeneous isotopic substitutions 17O3 and 18O3 of the ozone molecule are reported for the first time. The system being of strongly multiconfigurational character, multireference configuration interaction wave function ansatz with different complete active spaces was used. The reliable DBOC calculations with the targeted accuracy were possible to carry out up to about half of the dissociation threshold D0. The comparison with the experimental band centers shows a significant improvement of the accuracy with respect to the best Born-Oppenheimer (BO) ab initio calculations reducing the total root-mean-squares (calculated-observed) deviations by about a factor of two. For the set of 16O3 vibrations up to five bending and four stretching quanta, the mean (calculated-observed) deviations drop down from 0.7 cm-1 (BO) to about 0.1 cm-1, with the most pronounced improvement seen for bending states and for mixed bending-stretching polyads. In the case of bending band centers directly observed under high spectral resolutions, the errors are reduced by more than an order of magnitude down to 0.02 cm-1 from the observed levels, approaching nearly experimental accuracy. A similar improvement for heavy isotopologues shows that the reported DBOC corrections almost remove the systematic BO errors in vibrational levels below D0/2, though the scatter increases towards higher energies. The possible reasons for this finding, as well as remaining issues are discussed in detail. The reported results provide an encouraging accuracy validation for the multireference methods of the ab initio theory. New sets of ab initio vibrational states can be used for improving effective spectroscopic models for analyses of the observed high-resolution spectra, particularly in the cases of accidental resonances with "dark" states requiring accurate theoretical predictions.
Collapse
Affiliation(s)
- Attila Tajti
- ELTE Eötvös Loránd University, Institute of Chemistry, Laboratory of Theoretical Chemistry, P. O. Box 32, H-1518, Budapest 112, Hungary.
| | | | | | | |
Collapse
|
65
|
Domenichini G, von Rudorff GF, von Lilienfeld OA. Effects of perturbation order and basis set on alchemical predictions. J Chem Phys 2020; 153:144118. [PMID: 33086815 DOI: 10.1063/5.0023590] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
Alchemical perturbation density functional theory has been shown to be an efficient and computationally inexpensive way to explore chemical compound space. We investigate approximations made, in terms of atomic basis sets and the perturbation order, introduce an electron-density based estimate of errors of the alchemical prediction, and propose a correction for effects due to basis set incompleteness. Our numerical analysis of potential energy estimates, and resulting binding curves, is based on coupled-cluster single double (CCSD) reference results and is limited to all neutral diatomics with 14 electrons (AlH⋯NN). The method predicts binding energy, equilibrium distance, and vibrational frequencies of neighboring out-of-sample diatomics with near CCSD quality using perturbations up to the fifth order. We also discuss simultaneous alchemical mutations at multiple sites in benzene.
Collapse
Affiliation(s)
- Giorgio Domenichini
- Institute of Physical Chemistry and National Center for Computational Design and Discovery of Novel Materials (MARVEL), Department of Chemistry, University of Basel, 4056 Basel, Switzerland
| | - Guido Falk von Rudorff
- Institute of Physical Chemistry and National Center for Computational Design and Discovery of Novel Materials (MARVEL), Department of Chemistry, University of Basel, 4056 Basel, Switzerland
| | - O Anatole von Lilienfeld
- Institute of Physical Chemistry and National Center for Computational Design and Discovery of Novel Materials (MARVEL), Department of Chemistry, University of Basel, 4056 Basel, Switzerland
| |
Collapse
|
66
|
Mallick S, Roy B, Kumar P. A comparison of DLPNO-CCSD(T) and CCSD(T) method for the determination of the energetics of hydrogen atom transfer reactions. COMPUT THEOR CHEM 2020. [DOI: 10.1016/j.comptc.2020.112934] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
|
67
|
Vitale E, Alavi A, Kats D. FCIQMC-Tailored Distinguishable Cluster Approach. J Chem Theory Comput 2020; 16:5621-5634. [PMID: 32786911 PMCID: PMC7482318 DOI: 10.1021/acs.jctc.0c00470] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2020] [Indexed: 11/28/2022]
Abstract
The tailored approach is applied to the distinguishable cluster method together with a stochastic FCI solver (FCIQMC). It is demonstrated that the new method is more accurate than the corresponding tailored coupled cluster and the pure distinguishable cluster methods. An F12 correction for tailored methods and FCIQMC is introduced, which drastically improves the basis set convergence. A new black-box approach to define the active space using the natural orbitals from the distinguishable cluster is evaluated and found to be a convenient alternative to the usual CASSCF approach.
Collapse
Affiliation(s)
- Eugenio Vitale
- Max Planck Institute for
Solid State Research, Heisenbergstrasse 1, 70569 Stuttgart, Germany
| | - Ali Alavi
- Max Planck Institute for
Solid State Research, Heisenbergstrasse 1, 70569 Stuttgart, Germany
| | - Daniel Kats
- Max Planck Institute for
Solid State Research, Heisenbergstrasse 1, 70569 Stuttgart, Germany
| |
Collapse
|
68
|
Burton HGA, Thom AJW. Reaching Full Correlation through Nonorthogonal Configuration Interaction: A Second-Order Perturbative Approach. J Chem Theory Comput 2020; 16:5586-5600. [DOI: 10.1021/acs.jctc.0c00468] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Hugh G. A. Burton
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, U.K
| | - Alex J. W. Thom
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, U.K
| |
Collapse
|
69
|
Veccham SP, Head-Gordon M. Density Functionals for Hydrogen Storage: Defining the H2Bind275 Test Set with Ab Initio Benchmarks and Assessment of 55 Functionals. J Chem Theory Comput 2020; 16:4963-4982. [PMID: 32603109 DOI: 10.1021/acs.jctc.0c00292] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Efficient and high-capacity storage materials are indispensable for a hydrogen-based economy. In silico tools can accelerate the process of discovery of new adsorbent materials with optimal hydrogen adsorption enthalpies. Density functional theory is well-poised to become a very useful tool for enabling high-throughput screening of potential materials. In this work, we have identified density functional approximations that provide good performance for hydrogen binding applications following a two-pronged approach. First, we have compiled a data set (H2Bind275) that comprehensively represents the hydrogen binding problem capturing the chemical and mechanistic diversity in the binding sites encountered in hydrogen storage materials. We have also computed reference interaction energies for this data set using coupled-cluster theory. Second, we have assessed the performance of 55 density functional approximations for predicting H2 interaction energies and have identified two hybrid density functionals (ωB97X-V and ωB97M-V), two double hybrid density functionals (DSD-PBEPBE-D3(BJ) and PBE0-DH), and one semilocal density functional (B97M-V) as the best performing ones. We have recommended the addition of empirical dispersion corrections to systematically underbinding density functionals such as revPBE, BLYP, and B3LYP for improvements in performance at negligible additional cost. We have also recommended the usage of the def2-TZVPP basis set as it represents a good compromise between accuracy and cost, limiting the finite basis set errors to less than 1 kJ/mol.
Collapse
Affiliation(s)
- Srimukh Prasad Veccham
- Department of Chemistry, University of California, Berkeley, California 94720, United States.,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, United States.,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
70
|
Small DW. Remarkable Accuracy of an O( N6) Perturbative Correction to Opposite-Spin CCSD: Are Triples Necessary for Chemical Accuracy in Coupled Cluster? J Chem Theory Comput 2020; 16:4014-4020. [PMID: 32551642 DOI: 10.1021/acs.jctc.0c00244] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The focus of this work is OS-CCSD-SPT(2), which is a second-order similarity transformed perturbation theory correction to opposite spin coupled cluster singles doubles, where in the latter the same-spin amplitudes are removed and the opposite-spin ones are solved self-consistently. OS-CCSD-SPT(2) is free of empirical parameters, has an instrinsic scaling of O(N6), and makes no use of triples. We demonstrate that, for non-multireference molecules, OS-CCSD-SPT(2) produces relative energies whose accuracy is significantly higher than what is generally expected of a triples-free model. For example, using PBE0 orbitals in the reference, OS-CCSD-SPT(2) exhibits a mean absolute deviation (MAD) of 1.13 kcal/mol with respect to CCSD(2F) benchmark values for the non-multireference subset of W4-08 atomization energies (cf. a MAD > 6.5 kcal/mol for CCSD) and a MAD of 0.68 kcal/mol for the energies of reactions generated from the W4-08 molecules. These MADs are reduced to 0.61 and 0.63 kcal/mol, respectively, by a simple one-parameter spin-component scaling of the OS-CCSD-SPT(2) same-spin correlation energy. OS-CCSD is also naturally amenable to higher order corrections: the associated third-order correction, OS-CCSD-SPT(3), which does involve connected triples and quadruples, exhibits a MAD of 0.44 kcal/mol for the same atomization-energy benchmark.
Collapse
Affiliation(s)
- David W Small
- Molecular Graphics and Computation Facility, College of Chemistry, University of California, Berkeley, California 94720, United States
| |
Collapse
|
71
|
Karton A. Effective basis set extrapolations for CCSDT, CCSDT(Q), and CCSDTQ correlation energies. J Chem Phys 2020; 153:024102. [PMID: 32668917 DOI: 10.1063/5.0011674] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
It is well established that extrapolating the coupled-cluster single double triple [CCSD and (T)] correlation energies using empirically motivated extrapolation exponents can accelerate the basis set convergence. Here, we consider the extrapolation of coupled-cluster expansion terms beyond the CCSD(T) level to the complete basis set (CBS) limit. We obtain reference CCSDT-CCSD(T) [T3-(T)], CCSDT(Q)-CCSDT [(Q)], and CCSDTQ-CCSDT(Q) [T4-(Q)] contributions from cc-pV{5,6}Z extrapolations for a diverse set of 16 first- and second-row systems. We use these basis-set limit results to fit extrapolation exponents in conjunction with the cc-pV{D,T}Z, cc-pV{T,Q}Z, and cc-pV{Q,5}Z basis set pairs. The optimal extrapolation exponents result in noticeable improvements in performance (relative to α = 3.0) in conjunction with the cc-pV{T,Q}Z basis set pair; however, smaller improvements are obtained for the other basis sets. These results confirm that the basis sets and basis set extrapolations used for obtaining post-CCSD(T) components in composite thermochemical theories such as Weizmann-4 and HEAT are sufficiently close to the CBS limit for attaining sub-kJ/mole accuracy. The fitted extrapolation exponents demonstrate that the T3-(T) correlation component converges more slowly to the CBS limit than the (Q) and T4 terms. A systematic investigation of the effect of diffuse functions shows that it diminishes (i) in the order T3-(T) > (Q) > T4-(Q) and (ii) with the size of the basis set. Importantly, we find that diffuse functions tend to systematically reduce the T3-(T) contribution but systematically increases the (Q) contribution. Thus, the use of the cc-pVnZ basis sets benefits from a certain degree of error cancellation between these two components.
Collapse
Affiliation(s)
- Amir Karton
- School of Molecular Sciences, The University of Western Australia, Perth, WA 6009, Australia
| |
Collapse
|
72
|
Schurkus HF, Chan GKL, Chen DT, Cheng HP, Stanton JF. Theoretical prediction of magnetic exchange coupling constants from broken-symmetry coupled cluster calculations. J Chem Phys 2020; 152:234115. [PMID: 32571049 DOI: 10.1063/1.5144696] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Exchange coupling constants (J) are fundamental to the understanding of spin spectra of magnetic systems. Here, we investigate the broken-symmetry (BS) approaches of Noodleman and Yamaguchi in conjunction with coupled cluster (CC) methods to obtain exchange couplings. J values calculated from CC in this fashion converge smoothly toward the full configuration interaction result with increasing level of CC excitation. We compare this BS-CC scheme to the complementary equation-of-motion CC approach on a selection of bridged molecular cases and give results from a few other methodologies for context.
Collapse
Affiliation(s)
- Henry F Schurkus
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Garnet Kin-Lic Chan
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Dian-Teng Chen
- Department of Physics, University of Florida, Gainesville, Florida 32611, USA
| | - Hai-Ping Cheng
- Department of Physics, University of Florida, Gainesville, Florida 32611, USA
| | - John F Stanton
- Department of Chemistry, University of Florida, Gainesville, Florida 32611, USA
| |
Collapse
|
73
|
Neese F, Wennmohs F, Becker U, Riplinger C. The ORCA quantum chemistry program package. J Chem Phys 2020; 152:224108. [DOI: 10.1063/5.0004608] [Citation(s) in RCA: 697] [Impact Index Per Article: 174.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Frank Neese
- Max Planck Institut für Kohlenforschung, Kaiser-Wilhelm Platz 1, D-45470 Mülheim an der Ruhr, Germany
- FAccTs GmbH, Rolandstr. 67, 50677 Köln, Germany
| | - Frank Wennmohs
- Max Planck Institut für Kohlenforschung, Kaiser-Wilhelm Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Ute Becker
- Max Planck Institut für Kohlenforschung, Kaiser-Wilhelm Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | | |
Collapse
|
74
|
Matthews DA, Cheng L, Harding ME, Lipparini F, Stopkowicz S, Jagau TC, Szalay PG, Gauss J, Stanton JF. Coupled-cluster techniques for computational chemistry: The CFOUR program package. J Chem Phys 2020; 152:214108. [DOI: 10.1063/5.0004837] [Citation(s) in RCA: 214] [Impact Index Per Article: 53.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Affiliation(s)
- Devin A. Matthews
- Department of Chemistry, Southern Methodist University, Dallas, Texas 75275, USA
| | - Lan Cheng
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Michael E. Harding
- Institut für Physikalische Chemie, Karlsruher Institut für Technologie (KIT), Kaiserstr. 12, D-76131 Karlsruhe, Germany
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, I-56124 Pisa, Italy
| | - Stella Stopkowicz
- Department Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, D-55128 Mainz, Germany
| | - Thomas-C. Jagau
- Department of Chemistry, University of Munich (LMU), Butenandtstr. 5-13, D-81377 Munich, Germany
| | - Péter G. Szalay
- ELTE Eötvös Loránd University, Institute of Chemistry, Pázmány Péter sétány 1/A, H-1117 Budapest, Hungary
| | - Jürgen Gauss
- Department Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, D-55128 Mainz, Germany
| | - John F. Stanton
- Quantum Theory Project, Departments of Chemistry and Physics, University of Florida, Gainesville, Florida 32611, USA
| |
Collapse
|
75
|
Margócsy Á, Szabados Á. Ring coupled cluster doubles at the multireference level. J Chem Phys 2020; 152:204114. [PMID: 32486660 DOI: 10.1063/5.0005075] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A ring approximation within an internally contracted multireference (MR) Coupled Cluster (CC) framework is worked out and tested. Derivation of equations utilizes MR based, generalized normal ordering and the corresponding generalized Wick-theorem (MR-GWT). Contractions among cluster operators are avoided by adopting a normal ordered exponential ansatz. The original version of the MR ring CC doubles (MR-rCCD) equations [Á. Szabados and Á. Margócsy, Mol. Phys. 115, 2731 (2017)] is rectified in two aspects. On the one hand, over-completeness of double excitations is treated by relying on the concept of frames. On the other hand, restriction on the maximal cumulant rank is lifted from two to four. This is found essential for obtaining reliable correlation corrections to the energy. The MR function underlying the approach is provided by the Generalized Valence Bond (GVB) model. The pair structure of the reference ensures a fragment structure of GVB cumulants. This represents a benefit when evaluating cumulant contractions appearing as a consequence of MR-GWT. In particular, cumulant involving terms remain less expensive than their traditional, pair-contracted counterpart, facilitating an O(N6) eventual scaling of the proposed MR-rCCD method. Pilot applications are presented for covalent bond breaking, deprotonation energies, and torsional potentials.
Collapse
Affiliation(s)
- Á Margócsy
- Laboratory of Theoretical Chemistry, Institute of Chemistry, Faculty of Science, ELTE Eötvös Loránd University, Budapest, Hungary
| | - Á Szabados
- Laboratory of Theoretical Chemistry, Institute of Chemistry, Faculty of Science, ELTE Eötvös Loránd University, Budapest, Hungary
| |
Collapse
|
76
|
Saue T, Bast R, Gomes ASP, Jensen HJA, Visscher L, Aucar IA, Di Remigio R, Dyall KG, Eliav E, Fasshauer E, Fleig T, Halbert L, Hedegård ED, Helmich-Paris B, Iliaš M, Jacob CR, Knecht S, Laerdahl JK, Vidal ML, Nayak MK, Olejniczak M, Olsen JMH, Pernpointner M, Senjean B, Shee A, Sunaga A, van Stralen JNP. The DIRAC code for relativistic molecular calculations. J Chem Phys 2020; 152:204104. [PMID: 32486677 DOI: 10.1063/5.0004844] [Citation(s) in RCA: 118] [Impact Index Per Article: 29.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
DIRAC is a freely distributed general-purpose program system for one-, two-, and four-component relativistic molecular calculations at the level of Hartree-Fock, Kohn-Sham (including range-separated theory), multiconfigurational self-consistent-field, multireference configuration interaction, electron propagator, and various flavors of coupled cluster theory. At the self-consistent-field level, a highly original scheme, based on quaternion algebra, is implemented for the treatment of both spatial and time reversal symmetry. DIRAC features a very general module for the calculation of molecular properties that to a large extent may be defined by the user and further analyzed through a powerful visualization module. It allows for the inclusion of environmental effects through three different classes of increasingly sophisticated embedding approaches: the implicit solvation polarizable continuum model, the explicit polarizable embedding model, and the frozen density embedding model.
Collapse
Affiliation(s)
- Trond Saue
- Laboratoire de Chimie et Physique Quantique, UMR 5626 CNRS-Université Toulouse III-Paul Sabatier, 118 Route de Narbonne, F-31062 Toulouse, France
| | - Radovan Bast
- Department of Information Technology, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - André Severo Pereira Gomes
- Université de Lille, CNRS, UMR 8523-PhLAM-Physique des Lasers, Atomes et Molécules, F-59000 Lille, France
| | - Hans Jørgen Aa Jensen
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| | - Lucas Visscher
- Department of Chemistry and Pharmaceutical Sciences, Vrije Universiteit Amsterdam, NL-1081HV Amsterdam, The Netherlands
| | - Ignacio Agustín Aucar
- Instituto de Modelado e Innovación Tecnológica, CONICET, and Departamento de Física-Facultad de Ciencias Exactas y Naturales, UNNE, Avda. Libertad 5460, W3404AAS Corrientes, Argentina
| | - Roberto Di Remigio
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Kenneth G Dyall
- Dirac Solutions, 10527 NW Lost Park Drive, Portland, Oregon 97229, USA
| | - Ephraim Eliav
- School of Chemistry, Tel Aviv University, Ramat Aviv, Tel Aviv 69978, Israel
| | - Elke Fasshauer
- Department of Physics and Astronomy, Aarhus University, Ny Munkegade 120, 8000 Aarhus, Denmark
| | - Timo Fleig
- Laboratoire de Chimie et Physique Quantique, UMR 5626 CNRS-Université Toulouse III-Paul Sabatier, 118 Route de Narbonne, F-31062 Toulouse, France
| | - Loïc Halbert
- Université de Lille, CNRS, UMR 8523-PhLAM-Physique des Lasers, Atomes et Molécules, F-59000 Lille, France
| | - Erik Donovan Hedegård
- Division of Theoretical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Benjamin Helmich-Paris
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Miroslav Iliaš
- Department of Chemistry, Faculty of Natural Sciences, Matej Bel University, Tajovského 40, 974 01 Banská Bystrica, Slovakia
| | - Christoph R Jacob
- Technische Universität Braunschweig, Institute of Physical and Theoretical Chemistry, Gaußstr. 17, 38106 Braunschweig, Germany
| | - Stefan Knecht
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Jon K Laerdahl
- Department of Microbiology, Oslo University Hospital, Oslo, Norway
| | - Marta L Vidal
- Department of Chemistry, Technical University of Denmark, 2800 Kgs. Lyngby, Denmark
| | - Malaya K Nayak
- Theoretical Chemistry Section, Bhabha Atomic Research Centre, Trombay, Mumbai 400085, India
| | - Małgorzata Olejniczak
- Centre of New Technologies, University of Warsaw, S. Banacha 2c, 02-097 Warsaw, Poland
| | - Jógvan Magnus Haugaard Olsen
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | | | - Bruno Senjean
- Department of Chemistry and Pharmaceutical Sciences, Vrije Universiteit Amsterdam, NL-1081HV Amsterdam, The Netherlands
| | - Avijit Shee
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Ayaki Sunaga
- Department of Chemistry, Tokyo Metropolitan University, 1-1 Minami-Osawa, Hachioji-city, Tokyo 192-0397, Japan
| | - Joost N P van Stralen
- Department of Chemistry and Pharmaceutical Sciences, Vrije Universiteit Amsterdam, NL-1081HV Amsterdam, The Netherlands
| |
Collapse
|
77
|
de Jongh T, Besemer M, Shuai Q, Karman T, van der Avoird A, Groenenboom GC, van de Meerakker SYT. Imaging the onset of the resonance regime in low-energy NO-He collisions. Science 2020; 368:626-630. [DOI: 10.1126/science.aba3990] [Citation(s) in RCA: 50] [Impact Index Per Article: 12.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2019] [Accepted: 03/19/2020] [Indexed: 11/02/2022]
Affiliation(s)
- Tim de Jongh
- Institute for Molecules and Materials, Radboud University, 6525 AJ Nijmegen, Netherlands
| | - Matthieu Besemer
- Institute for Molecules and Materials, Radboud University, 6525 AJ Nijmegen, Netherlands
| | - Quan Shuai
- Institute for Molecules and Materials, Radboud University, 6525 AJ Nijmegen, Netherlands
| | - Tijs Karman
- Institute for Theoretical Atomic Molecular and Optical Physics, Center For Astrophysics, Harvard and Smithsonian, Cambridge, MA 02138, USA
| | - Ad van der Avoird
- Institute for Molecules and Materials, Radboud University, 6525 AJ Nijmegen, Netherlands
| | - Gerrit C. Groenenboom
- Institute for Molecules and Materials, Radboud University, 6525 AJ Nijmegen, Netherlands
| | | |
Collapse
|
78
|
Loos PF, Scemama A, Boggio-Pasqua M, Jacquemin D. Mountaineering Strategy to Excited States: Highly Accurate Energies and Benchmarks for Exotic Molecules and Radicals. J Chem Theory Comput 2020; 16:3720-3736. [DOI: 10.1021/acs.jctc.0c00227] [Citation(s) in RCA: 37] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Affiliation(s)
- Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques, CNRS, UPS, Université de Toulouse, 31000 Toulouse, France
| | - Anthony Scemama
- Laboratoire de Chimie et Physique Quantiques, CNRS, UPS, Université de Toulouse, 31000 Toulouse, France
| | - Martial Boggio-Pasqua
- Laboratoire de Chimie et Physique Quantiques, CNRS, UPS, Université de Toulouse, 31000 Toulouse, France
| | - Denis Jacquemin
- CEISAM UMR 6230, CNRS, Université de Nantes, F-44000 Nantes, France
| |
Collapse
|
79
|
Franke PR, Duncan MA, Douberly GE. Infrared photodissociation spectroscopy and anharmonic vibrational study of the HO 4 + molecular ion. J Chem Phys 2020; 152:174309. [PMID: 32384862 DOI: 10.1063/5.0005975] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Molecular cations of HO4 + and DO4 + are produced in a supersonic expansion. They are mass-selected, and infrared photodissociation spectra of these species are measured with the aid of argon-tagging. Although previous theoretical studies have modeled these systems as proton-bound dimers of molecular oxygen, infrared spectra have free OH stretching bands, suggesting other isomeric structures. As a consequence, we undertook extensive computational studies. Our conformer search used a composite method based on an economical combination of single- and multi-reference theories. Several conformers were located on the quintet, triplet, and singlet surfaces, spanning in energy of only a few thousand wavenumbers. Most of the singlet and triplet conformers have pronounced multiconfigurational character. Previously unidentified covalent-like structures (H-O-O-O-O) on the singlet and triplet surfaces likely represent the global minima. In our experiments, HO4 + is formed in a relatively hot environment, and similar experiments have been shown capable of producing multiple conformers in low-lying electronic states. None of the predicted HO4 + isomers can be ruled out a priori based on energetic arguments. We interpret our argon-tagged spectra with Second-Order Vibrational Perturbation Theory with Resonances (VPT2+K). The presence of one or more covalent-like isomers is the only reasonable explanation for the spectral features observed.
Collapse
Affiliation(s)
- Peter R Franke
- Department of Chemistry, University of Georgia, Athens, Georgia 30602, USA
| | - Michael A Duncan
- Department of Chemistry, University of Georgia, Athens, Georgia 30602, USA
| | - Gary E Douberly
- Department of Chemistry, University of Georgia, Athens, Georgia 30602, USA
| |
Collapse
|
80
|
Yu F, Wang Y. Dual‐hybrid direct random phase approximation and second‐order screened exchange with nonlocal van der Waals correlations for noncovalent interactions. J Comput Chem 2020; 41:1018-1025. [DOI: 10.1002/jcc.26149] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2019] [Revised: 01/05/2020] [Accepted: 01/06/2020] [Indexed: 11/09/2022]
Affiliation(s)
- Feng Yu
- Department of Physics, School of ScienceXi'an Technological University Xi'an Shaanxi China
| | - Yaoting Wang
- Department of Physics, School of ScienceXi'an Technological University Xi'an Shaanxi China
| |
Collapse
|
81
|
Werner HJ, Knowles PJ, Manby FR, Black JA, Doll K, Heßelmann A, Kats D, Köhn A, Korona T, Kreplin DA, Ma Q, Miller TF, Mitrushchenkov A, Peterson KA, Polyak I, Rauhut G, Sibaev M. The Molpro quantum chemistry package. J Chem Phys 2020; 152:144107. [PMID: 32295355 DOI: 10.1063/5.0005081] [Citation(s) in RCA: 502] [Impact Index Per Article: 125.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
Molpro is a general purpose quantum chemistry software package with a long development history. It was originally focused on accurate wavefunction calculations for small molecules but now has many additional distinctive capabilities that include, inter alia, local correlation approximations combined with explicit correlation, highly efficient implementations of single-reference correlation methods, robust and efficient multireference methods for large molecules, projection embedding, and anharmonic vibrational spectra. In addition to conventional input-file specification of calculations, Molpro calculations can now be specified and analyzed via a new graphical user interface and through a Python framework.
Collapse
Affiliation(s)
- Hans-Joachim Werner
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Peter J Knowles
- School of Chemistry, Cardiff University, Main Building, Park Place, Cardiff CF10 3AT, United Kingdom
| | - Frederick R Manby
- School of Chemistry, University of Bristol, Cantock's Close, Bristol BS8 1TS, United Kingdom
| | - Joshua A Black
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Klaus Doll
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Andreas Heßelmann
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Daniel Kats
- Max-Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Andreas Köhn
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Tatiana Korona
- Faculty of Chemistry, University of Warsaw, L. Pasteura 1 St., 02-093 Warsaw, Poland
| | - David A Kreplin
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Qianli Ma
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Thomas F Miller
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | | | - Kirk A Peterson
- Department of Chemistry, Washington State University, Pullman, Washington 99164-4630, USA
| | - Iakov Polyak
- School of Chemistry, Cardiff University, Main Building, Park Place, Cardiff CF10 3AT, United Kingdom
| | - Guntram Rauhut
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Marat Sibaev
- School of Chemistry, Cardiff University, Main Building, Park Place, Cardiff CF10 3AT, United Kingdom
| |
Collapse
|
82
|
Stein CJ, Reiher M. Semiclassical Dispersion Corrections Efficiently Improve Multiconfigurational Theory with Short-Range Density-Functional Dynamic Correlation. J Phys Chem A 2020; 124:2834-2841. [PMID: 32186877 DOI: 10.1021/acs.jpca.0c02130] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Abstract
Multiconfigurational wave functions are known to describe the electronic structure across a Born-Oppenheimer surface qualitatively correct. However, for quantitative reaction energies, dynamic correlation originating from the many configurations involving excitations out of the restricted orbital space, the active space, must be considered. Standard procedures involve approximations that eventually limit the ultimate accuracy achievable (most prominently, multireference perturbation theory). At the same time, the computational cost increases dramatically due to the necessity to obtain higher-order reduced density matrices. It is this disproportion that leads us here to propose an MC-srDFT-D hybrid approach of semiclassical dispersion (D) corrections to cover long-range dynamic correlation in a multiconfigurational (MC) wave function theory, which includes short-range (sr) dynamic correlation by density functional theory (DFT) without double counting. We demonstrate that the reliability of this approach is very good (at negligible cost), especially when considering that standard second-order multireference perturbation theory usually overestimates dispersion interactions.
Collapse
Affiliation(s)
- Christopher J Stein
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Markus Reiher
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
83
|
Zheng X, Liu J, Doumy G, Young L, Cheng L. Hetero-site Double Core Ionization Energies with Sub-electronvolt Accuracy from Delta-Coupled-Cluster Calculations. J Phys Chem A 2020; 124:4413-4426. [DOI: 10.1021/acs.jpca.0c00901] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Xuechen Zheng
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Junzi Liu
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Gilles Doumy
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Linda Young
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Lan Cheng
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, United States
| |
Collapse
|
84
|
Sajid A, Ford MJ, Reimers JR. Single-photon emitters in hexagonal boron nitride: a review of progress. REPORTS ON PROGRESS IN PHYSICS. PHYSICAL SOCIETY (GREAT BRITAIN) 2020; 83:044501. [PMID: 31846956 DOI: 10.1088/1361-6633/ab6310] [Citation(s) in RCA: 41] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/13/2023]
Abstract
This report summarizes progress made in understanding properties such as zero-phonon-line energies, emission and absorption polarizations, electron-phonon couplings, strain tuning and hyperfine coupling of single photon emitters in hexagonal boron nitride. The primary aims of this research are to discover the chemical nature of the emitting centres and to facilitate deployment in device applications. Critical analyses of the experimental literature and data interpretation, as well as theoretical approaches used to predict properties, are made. In particular, computational and theoretical limitations and challenges are discussed, with a range of suggestions made to overcome these limitations, striving to achieve realistic predictions concerning the nature of emitting centers. A symbiotic relationship is required in which calculations focus on properties that can easily be measured, whilst experiments deliver results in a form facilitating mass-produced calculations.
Collapse
Affiliation(s)
- A Sajid
- University of Technology Sydney, School of Mathematical and Physical Sciences, Ultimo, New South Wales 2007, Australia. CAMD, Department of Physics, Technical University of Denmark, 2800 Kgs. Lyngby, Denmark. Department of Physics, GC University Faisalabad, Allama Iqbal Road, 38000 Faisalabad, Pakistan. Author to whom any correspondence should be addressed
| | | | | |
Collapse
|
85
|
Chou CW, Collopy AL, Kurz C, Lin Y, Harding ME, Plessow PN, Fortier T, Diddams S, Leibfried D, Leibrandt DR. Frequency-comb spectroscopy on pure quantum states of a single molecular ion. Science 2020; 367:1458-1461. [PMID: 32217722 PMCID: PMC10652508 DOI: 10.1126/science.aba3628] [Citation(s) in RCA: 26] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/27/2019] [Accepted: 03/04/2020] [Indexed: 01/21/2023]
Abstract
Spectroscopy is a powerful tool for studying molecules and is commonly performed on large thermal molecular ensembles that are perturbed by motional shifts and interactions with the environment and one another, resulting in convoluted spectra and limited resolution. Here, we use quantum-logic techniques to prepare a trapped molecular ion in a single quantum state, drive terahertz rotational transitions with an optical frequency comb, and read out the final state nondestructively, leaving the molecule ready for further manipulation. We can resolve rotational transitions to 11 significant digits and derive the rotational constant of 40CaH+ to be B R = 142 501 777.9(1.7) kilohertz. Our approach is suited for a wide range of molecular ions, including polyatomics and species relevant for tests of fundamental physics, chemistry, and astrophysics.
Collapse
Affiliation(s)
- C W Chou
- Time and Frequency Division, National Institute of Standards and Technology, Boulder, CO 80305, USA.
| | - A L Collopy
- Time and Frequency Division, National Institute of Standards and Technology, Boulder, CO 80305, USA
| | - C Kurz
- Time and Frequency Division, National Institute of Standards and Technology, Boulder, CO 80305, USA
| | - Y Lin
- CAS Key Laboratory of Microscale Magnetic Resonance and Department of Modern Physics, University of Science and Technology of China, Hefei 230026, China
- Time and Frequency Division, National Institute of Standards and Technology, Boulder, CO 80305, USA
- Hefei National Laboratory for Physical Sciences at the Microscale, University of Science and Technology of China, Hefei 230026, China
- Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei 230026, China
| | - M E Harding
- Institute of Nanotechnology, Karlsruhe Institute of Technology (KIT), 76021 Karlsruhe, Germany
| | - P N Plessow
- Institute of Catalysis Research and Technology, Karlsruhe Institute of Technology (KIT), 76021 Karlsruhe, Germany
| | - T Fortier
- Time and Frequency Division, National Institute of Standards and Technology, Boulder, CO 80305, USA
- Department of Physics, University of Colorado, Boulder, CO 80309, USA
| | - S Diddams
- Time and Frequency Division, National Institute of Standards and Technology, Boulder, CO 80305, USA
- Department of Physics, University of Colorado, Boulder, CO 80309, USA
| | - D Leibfried
- Time and Frequency Division, National Institute of Standards and Technology, Boulder, CO 80305, USA
- Department of Physics, University of Colorado, Boulder, CO 80309, USA
| | - D R Leibrandt
- Time and Frequency Division, National Institute of Standards and Technology, Boulder, CO 80305, USA
- Department of Physics, University of Colorado, Boulder, CO 80309, USA
| |
Collapse
|
86
|
Weidman JD, Turney JM, Schaefer HF. Energetics and mechanisms for the acetonyl radical + O 2 reaction: An important system for atmospheric and combustion chemistry. J Chem Phys 2020; 152:114301. [PMID: 32199416 DOI: 10.1063/1.5141859] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Abstract
The acetonyl radical (•CH2COCH3) is relevant to atmospheric and combustion chemistry due to its prevalence in many important reaction mechanisms. One such reaction mechanism is the decomposition of Criegee intermediates in the atmosphere that can produce acetonyl radical and OH. In order to understand the fate of the acetonyl radical in these environments and to create more accurate kinetics models, we have examined the reaction system of the acetonyl radical with O2 using highly reliable theoretical methods. Structures were optimized using coupled cluster theory with singles, doubles, and perturbative triples [CCSD(T)] with an atomic natural orbital (ANO0) basis set. Energetics were computed to chemical accuracy using the focal point approach involving perturbative treatment of quadruple excitations [CCSDT(Q)] and basis sets as large as cc-pV5Z. The addition of O2 to the acetonyl radical produces the acetonylperoxy radical, and multireference computations on this reaction suggest it to be barrierless. No submerged pathways were found for the unimolecular isomerization of the acetonylperoxy radical. Besides dissociation to reactants, the lowest energy pathway available for the acetonylperoxy radical is a 1-5 H shift from the methyl group to the peroxy group through a transition state that is 3.3 kcal mol-1 higher in energy than acetonyl radical + O2. The ultimate products from this pathway are the enol tautomer of the acetonyl radical along with O2. Multiple pathways that lead to OH formation are considered; however, all of these pathways are predicted to be energetically inaccessible, except at high temperatures.
Collapse
Affiliation(s)
- Jared D Weidman
- Center for Computational Quantum Chemistry, University of Georgia, Athens, Georgia 30602, USA
| | - Justin M Turney
- Center for Computational Quantum Chemistry, University of Georgia, Athens, Georgia 30602, USA
| | - Henry F Schaefer
- Center for Computational Quantum Chemistry, University of Georgia, Athens, Georgia 30602, USA
| |
Collapse
|
87
|
Wu CH, Magers DB, Harding LB, Klippenstein SJ, Allen WD. Reaction Profiles and Kinetics for Radical-Radical Hydrogen Abstraction via Multireference Coupled Cluster Theory. J Chem Theory Comput 2020; 16:1511-1525. [PMID: 32073856 DOI: 10.1021/acs.jctc.9b00966] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Radical-radical abstractions in hydrocarbon oxidation chemistry are disproportionation reactions that are generally exothermic with little or no barrier yet are underappreciated and poorly studied. Such challenging multireference electronic structure problems are tackled here using the recently developed state-specific multireference coupled cluster methods Mk-MRCCSD and Mk-MRCCSD(T), as well as the companion perturbation theory Mk-MRPT2 and the established MRCISD, MRCISD+Q, and CASPT2 approaches. Reaction paths are investigated for five prototypes involving radical-radical hydrogen abstraction: H + BeH → H2+ Be, H + NH2 → H2 + NH, CH3 + C2H5 → CH4 + C2H4, H + C2H5 → H2 + C2H4, and H + HCO → H2 + CO. Full configuration interaction (FCI) benchmark computations for the H + BeH, H + NH2, and H + HCO reactions prove that Mk-MRCCSD(T) provides superior accuracy for the interaction energies in the entrance channel, with mean absolute errors less than 0.3 kcal mol-1 and percentage deviations less than 10% over the fragment separations of relevance to kinetics. To facilitate combustion studies, energetics for the CH3 + C2H5, H + C2H5, and H + HCO reactions were computed at each level of theory with correlation-consistent basis sets (cc-pVXZ, X = T, Q, 5) and extrapolated to the complete basis set (CBS) limit. These CBS energies were coupled with CASPT2 projected vibrational frequencies along a minimum energy path to obtain rate constants for these three reactions. The rigorous Mk-MRCCSD(T)/CBS results demonstrate unequivocally that these three reactions proceed with no barrier in the entrance channel, contrary to some earlier predictions. Mk-MRCCSD(T) also reveals that the economical CASPT2 method performs well for large interfragment separations but may deteriorate substantially at shorter distances.
Collapse
Affiliation(s)
- Chia-Hua Wu
- Center for Computational Quantum Chemistry and Department of Chemistry, University of Georgia, Athens, Georgia 30602, United States
| | - D Brandon Magers
- Center for Computational Quantum Chemistry and Department of Chemistry, University of Georgia, Athens, Georgia 30602, United States.,Department of Chemistry and Physics, Belhaven University, Jackson, Mississippi 39202, United States
| | - Lawrence B Harding
- 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
| | - Wesley D Allen
- Center for Computational Quantum Chemistry and Department of Chemistry, University of Georgia, Athens, Georgia 30602, United States
| |
Collapse
|
88
|
Battey SR, Bross DH, Peterson KA, Persinger TD, VanGundy RA, Heaven MC. Spectroscopic and theoretical studies of UN and UN . J Chem Phys 2020; 152:094302. [PMID: 33480743 DOI: 10.1063/1.5144299] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The low-energy electronic states of UN and UN+ have been examined using high-level electronic structure calculations and two-color photoionization techniques. The experimental measurements provided an accurate ionization energy for UN (IE = 50 802 ± 5 cm-1). Spectra for UN+ yielded ro-vibrational constants and established that the ground state has the electronic angular momentum projection Ω = 4. Ab initio calculations were carried out using the spin-orbit state interacting approach with the complete active space second-order perturbation theory method. A series of correlation consistent basis sets were used in conjunction with small-core relativistic pseudopotentials on U to extrapolate to the complete basis set limits. The results for UN correctly obtained an Ω = 3.5 ground state and demonstrated a high density of configurationally related excited states with closely similar ro-vibrational constants. Similar results were obtained for UN+, with reduced complexity owing to the smaller number of outer-shell electrons. The calculated IE for UN was in excellent agreement with the measured value. Improved values for the dissociation energies of UN and UN+, as well as their heats of formation, were obtained using the Feller-Peterson-Dixon composite thermochemistry method, including corrections up through coupled cluster singles, doubles, triples and quadruples. An analysis of the ab initio results from the perspective of the ligand field theory shows that the patterns of electronic states for both UN and UN+ can be understood in terms of the underlying energy level structure of the atomic metal ion.
Collapse
Affiliation(s)
- S R Battey
- Department of Chemistry, Washington State University, Pullman, Washington 99164-4630, USA
| | - D H Bross
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, USA
| | - K A Peterson
- Department of Chemistry, Washington State University, Pullman, Washington 99164-4630, USA
| | - T D Persinger
- Department of Chemistry, Emory University, Atlanta, Georgia 30322, USA
| | - R A VanGundy
- Department of Chemistry, Emory University, Atlanta, Georgia 30322, USA
| | - M C Heaven
- Department of Chemistry, Emory University, Atlanta, Georgia 30322, USA
| |
Collapse
|
89
|
Armentrout PB, Peterson KA. Guided Ion Beam and Quantum Chemical Investigation of the Thermochemistry of Thorium Dioxide Cations: Thermodynamic Evidence for Participation of f Orbitals in Bonding. Inorg Chem 2020; 59:3118-3131. [PMID: 32083480 DOI: 10.1021/acs.inorgchem.9b03488] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Kinetic energy dependent reactions of ThO+ with O2 are studied using a guided ion beam tandem mass spectrometer. The formation of ThO2+ in the reaction of ThO+ with O2 is observed to be slightly endothermic and also exhibits two obvious features in the cross section. These kinetic energy dependent cross sections were modeled to determine a 0 K bond dissociation energy of D0(OTh+-O) = 4.94 ± 0.06 eV. This value is slightly larger but within experimental uncertainty of less precise previously reported experimental values. The higher energy feature in the ThO2+ cross section was also analyzed and suggests formation of an excited state of the product ion lying 3.1 ± 0.2 eV above the ground state. Additionally, the thermochemistry of ThO2+ was explored by quantum chemical calculations, including a full Feller-Peterson-Dixon (FPD) composite approach with correlation contributions up to CCSDT(Q) and four-component spin-orbit corrections, as well as more approximate CCSD(T) calculations including semiempirical estimates of spin-orbit energy contributions. The FPD approach predicts D0(OTh+-O) = 4.87 ± 0.04 eV, in good agreement with the experimental value. Analogous FPD results for ThO+, ThO, and ThO2 are also presented, including ionization energies for both ThO and ThO2. The ThO2+ bond energy is larger than those of its transition metal congeners, TiO2+ and ZrO2+, which can be attributed partially to an actinide contraction, but also to contributions from the participation of f orbitals on thorium that are unavailable to the transition metal systems.
Collapse
Affiliation(s)
- P B Armentrout
- Department of Chemistry, University of Utah, Salt Lake City, Utah 84112-0850, United States
| | - Kirk A Peterson
- Department of Chemistry, Washington State University, Pullman, Washington 99164-4630, United States
| |
Collapse
|
90
|
Hansen AS, Baardsen G, Rebolini E, Maschio L, Pedersen TB. Representation of the virtual space in extended systems – a correlation energy convergence study. Mol Phys 2020. [DOI: 10.1080/00268976.2020.1733118] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
Affiliation(s)
- A. S. Hansen
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, Oslo, Norway
| | - G. Baardsen
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, Oslo, Norway
| | - E. Rebolini
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, Oslo, Norway
- Institut Laue Langevin, Grenoble, France
| | - L. Maschio
- Dipartimento di Chimica, Universitá di Torino, Torino, Italy
| | - T. B. Pedersen
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, Oslo, Norway
| |
Collapse
|
91
|
Kállay M, Nagy PR, Mester D, Rolik Z, Samu G, Csontos J, Csóka J, Szabó PB, Gyevi-Nagy L, Hégely B, Ladjánszki I, Szegedy L, Ladóczki B, Petrov K, Farkas M, Mezei PD, Ganyecz Á. The MRCC program system: Accurate quantum chemistry from water to proteins. J Chem Phys 2020; 152:074107. [DOI: 10.1063/1.5142048] [Citation(s) in RCA: 152] [Impact Index Per Article: 38.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
Affiliation(s)
- Mihály Kállay
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Péter R. Nagy
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Dávid Mester
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Zoltán Rolik
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Gyula Samu
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - József Csontos
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - József Csóka
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - P. Bernát Szabó
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - László Gyevi-Nagy
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Bence Hégely
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - István Ladjánszki
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Lóránt Szegedy
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Bence Ladóczki
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Klára Petrov
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Máté Farkas
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Pál D. Mezei
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Ádám Ganyecz
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| |
Collapse
|
92
|
Maranzana A, Tonachini G. Multireference Study of the H 2COO (Criegee Intermediate) + O 3 Addition: A Reaction of Possible Tropospheric Interest. J Phys Chem A 2020; 124:1112-1120. [DOI: 10.1021/acs.jpca.9b11430] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Andrea Maranzana
- Università di Torino, Dipartimento di Chimica, Corso Massimo D’Azeglio, 48, I-10125 Torino, Italy
| | - Glauco Tonachini
- Università di Torino, Dipartimento di Chimica, Corso Massimo D’Azeglio, 48, I-10125 Torino, Italy
| |
Collapse
|
93
|
Chen JL, Sun T, Wang YB, Wang W. Toward a less costly but accurate calculation of the CCSD(T)/CBS noncovalent interaction energy. J Comput Chem 2020; 41:1252-1260. [PMID: 32045021 DOI: 10.1002/jcc.26171] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2019] [Revised: 01/12/2020] [Accepted: 02/01/2020] [Indexed: 01/14/2023]
Abstract
The popular method of calculating the noncovalent interaction energies at the coupled-cluster single-, double-, and perturbative triple-excitations [CCSD(T)] theory level in the complete basis set (CBS) limit was to add a CCSD(T) correction term to the CBS second-order Møller-Plesset perturbation theory (MP2). The CCSD(T) correction term is the difference between the CCSD(T) and MP2 interaction energies evaluated in a medium basis set. However, the CCSD(T) calculations with the medium basis sets are still very expensive for systems with more than 30 atoms. Comparatively, the domain-based local pair natural orbital coupled-cluster method [DLPNO-CCSD(T)] can be applied to large systems with over 1,000 atoms. Considering both the computational accuracy and efficiency, in this work, we propose a new scheme to calculate the CCSD(T)/CBS interaction energies. In this scheme, the MP2/CBS term keeps intact and the CCSD(T) correction term is replaced by a DLPNO-CCSD(T) correction term which is the difference between the DLPNO-CCSD(T) and DLPNO-MP2 interaction energies evaluated in a medium basis set. The interaction energies of the noncovalent systems in the S22, HSG, HBC6, NBC10, and S66 databases were recalculated employing this new scheme. The consistent and tight settings of the truncation parameters for DLPNO-CCSD(T) and DLPNO-MP2 in this noncanonical CCSD(T)/CBS calculations lead to the maximum absolute deviation and root-mean-square deviation from the canonical CCSD(T)/CBS interaction energies of less than or equal to 0.28 kcal/mol and 0.09 kcal/mol, respectively. The high accuracy and low cost of this new computational scheme make it an excellent candidate for the study of large noncovalent systems.
Collapse
Affiliation(s)
- Jiu-Li Chen
- Department of Chemistry, and Key Laboratory of Guizhou High Performance Computational Chemistry, Guizhou University, Guiyang, China
| | - Tao Sun
- Department of Chemistry, and Key Laboratory of Guizhou High Performance Computational Chemistry, Guizhou University, Guiyang, China
| | - Yi-Bo Wang
- Department of Chemistry, and Key Laboratory of Guizhou High Performance Computational Chemistry, Guizhou University, Guiyang, China
| | - Weizhou Wang
- College of Chemistry and Chemical Engineering, and Henan Key Laboratory of Function-Oriented Porous Materials, Luoyang Normal University, Luoyang, China
| |
Collapse
|
94
|
Förster A, Franchini M, van Lenthe E, Visscher L. A Quadratic Pair Atomic Resolution of the Identity Based SOS-AO-MP2 Algorithm Using Slater Type Orbitals. J Chem Theory Comput 2020; 16:875-891. [PMID: 31930915 PMCID: PMC7027358 DOI: 10.1021/acs.jctc.9b00854] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2019] [Indexed: 01/04/2023]
Abstract
We report a production level implementation of pair atomic resolution of the identity (PARI) based second-order Møller-Plesset perturbation theory (MP2) in the Slater type orbital (STO) based Amsterdam Density Functional (ADF) code. As demonstrated by systematic benchmarks, dimerization and isomerization energies obtained with our code using STO basis sets of triple-ζ-quality show mean absolute deviations from Gaussian type orbital, canonical, basis set limit extrapolated, global density fitting (DF)-MP2 results of less than 1 kcal/mol. Furthermore, we introduce a quadratic scaling atomic orbital based spin-opposite-scaled (SOS)-MP2 approach with a very small prefactor. Due to a worst-case scaling of [Formula: see text], our implementation is very fast already for small systems and shows an exceptionally early crossover to canonical SOS-PARI-MP2. We report computational wall time results for linear as well as for realistic three-dimensional molecules and show that triple-ζ quality calculations on molecules of several hundreds of atoms are only a matter of a few hours on a single compute node, the bottleneck of the computations being the SCF rather than the post-SCF energy correction.
Collapse
Affiliation(s)
- Arno Förster
- Theoretical Chemistry, Vrije
Universiteit, De Boelelaan 1083, NL-1081 HV Amsterdam, The
Netherlands
| | - Mirko Franchini
- Theoretical Chemistry, Vrije
Universiteit, De Boelelaan 1083, NL-1081 HV Amsterdam, The
Netherlands
- Scientific Computing & Modelling
NV, De Boelelaan 1083, NL-1081 HV Amsterdam, The
Netherlands
| | - Erik van Lenthe
- Scientific Computing & Modelling
NV, De Boelelaan 1083, NL-1081 HV Amsterdam, The
Netherlands
| | - Lucas Visscher
- Theoretical Chemistry, Vrije
Universiteit, De Boelelaan 1083, NL-1081 HV Amsterdam, The
Netherlands
| |
Collapse
|
95
|
Sommerfeld T, Davis MC. Excluded-volume descriptors for dipole-bound anions: Amine N-oxides as a test case. J Chem Phys 2020; 152:054102. [PMID: 32035463 DOI: 10.1063/1.5142624] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Dipole-bound anions can be theoretically characterized at three fundamentally different levels. The highest are ab initio calculations, which themselves range from fairly approximate, say, Koopmans's Theorem (KT) or second-order Møller-Plesset perturbation theory, to highly sophisticated, say, the electron affinity equation-of-motion couple-cluster with single, double, and perturbative triple substitutions, which rivals experiments in reliability. The next level down is represented by one-electron model Hamiltonians. Again, one-electron model Hamiltonians can be fairly approximate, especially if the molecular system is modeled by a simple point-dipole and point-polarizable site; however, very reliable models have been developed for specific systems, for example, water clusters. At the lowest level, one can qualitatively explain trends in classes of dipole-bound anions in terms of the dipole moment, μ, the polarizability, α, and the so-called excluded volume, Vx. This project aims at the qualitative level. While the dipole moment and the polarizability possess clear-cut definitions, the excluded volume must-similar to all molecular volumes-remain a rather vaguely defined term, and so far, we are unaware of any quantitative definition in the literature. Here, we introduce and investigate three descriptors for Vx. To this end, we first establish a dataset with consistent ab initio results for 25 amine N-oxides structures. Then, we demonstrate that the descriptors are indeed able to explain trends for sets of isomers and conformers and investigate to what extent the descriptors are able to predict electron binding energy of dipole-bound states using simple quantitative structure-property relationship-like models. It turns out that μ and Vx provide a reasonably accurate prediction of the electrostatic part of the electron bind energy (the KT value) and that the polarizability α provides an acceptable prediction of the electron correlation contribution.
Collapse
Affiliation(s)
- Thomas Sommerfeld
- Department of Chemistry and Physics, Southeastern Louisiana University, SLU 10878, Hammond, Louisiana 70402, USA
| | - Megan C Davis
- Department of Chemistry and Physics, Southeastern Louisiana University, SLU 10878, Hammond, Louisiana 70402, USA
| |
Collapse
|
96
|
Loos PF, Lipparini F, Boggio-Pasqua M, Scemama A, Jacquemin D. A Mountaineering Strategy to Excited States: Highly Accurate Energies and Benchmarks for Medium Sized Molecules. J Chem Theory Comput 2020; 16:1711-1741. [PMID: 31986042 DOI: 10.1021/acs.jctc.9b01216] [Citation(s) in RCA: 107] [Impact Index Per Article: 26.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Following our previous work focusing on compounds containing up to 3 non-hydrogen atoms [J. Chem. Theory Comput. 2018, 14, 4360-4379], we present here highly accurate vertical transition energies obtained for 27 molecules encompassing 4, 5, and 6 non-hydrogen atoms: acetone, acrolein, benzene, butadiene, cyanoacetylene, cyanoformaldehyde, cyanogen, cyclopentadiene, cyclopropenone, cyclopropenethione, diacetylene, furan, glyoxal, imidazole, isobutene, methylenecyclopropene, propynal, pyrazine, pyridazine, pyridine, pyrimidine, pyrrole, tetrazine, thioacetone, thiophene, thiopropynal, and triazine. To obtain these energies, we use equation-of-motion/linear-response coupled cluster theory up to the highest technically possible excitation order for these systems (CC3, EOM-CCSDT, and EOM-CCSDTQ) and selected configuration interaction (SCI) calculations (with tens of millions of determinants in the reference space), as well as the multiconfigurational n-electron valence state perturbation theory (NEVPT2) method. All these approaches are applied in combination with diffuse-containing atomic basis sets. For all transitions, we report at least CC3/aug-cc-pVQZ vertical excitation energies as well as CC3/aug-cc-pVTZ oscillator strengths for each dipole-allowed transition. We show that CC3 almost systematically delivers transition energies in agreement with higher-level methods with a typical deviation of ±0.04 eV, except for transitions with a dominant double excitation character where the error is much larger. The present contribution gathers a large, diverse, and accurate set of more than 200 highly accurate transition energies for states of various natures (valence, Rydberg, singlet, triplet, n → π*, π → π*, ...). We use this series of theoretical best estimates to benchmark a series of popular methods for excited state calculations: CIS(D), ADC(2), CC2, STEOM-CCSD, EOM-CCSD, CCSDR(3), CCSDT-3, CC3, and NEVPT2. The results of these benchmarks are compared to the available literature data.
Collapse
Affiliation(s)
- Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques, CNRS et Université Toulouse III - Paul Sabatier, 118 route de Narbonne, 31062 Toulouse, France
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, Via Moruzzi 3, 56124 Pisa, Italy
| | - Martial Boggio-Pasqua
- Laboratoire de Chimie et Physique Quantiques, CNRS et Université Toulouse III - Paul Sabatier, 118 route de Narbonne, 31062 Toulouse, France
| | - Anthony Scemama
- Laboratoire de Chimie et Physique Quantiques, CNRS et Université Toulouse III - Paul Sabatier, 118 route de Narbonne, 31062 Toulouse, France
| | - Denis Jacquemin
- CEISAM Lab, UMR 6230, Université de Nantes, CNRS, F-44000 Nantes, France
| |
Collapse
|
97
|
Job N, Karton A, Thirumoorthy K, Cooksy AL, Thimmakondu VS. Theoretical Studies of SiC 4H 2 Isomers Delineate Three Low-Lying Silylidenes Are Missing in the Laboratory. J Phys Chem A 2020; 124:987-1002. [PMID: 31904236 DOI: 10.1021/acs.jpca.9b11742] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Eleven isomers of SiC4H2 lying within 50 kcal mol-1 have been theoretically investigated using density functional theory and high-level coupled-cluster methods. Among them, four isomers, 1-ethynyl-3-silacycloprop-1(2)-en-3-ylidene (1), diethynylsilylidene (2), 1-sila-1,2,3,4-pentatetraenylidene (4), and 1,3-butadiynylsilylidene (5), have already been identified in the laboratory. The current investigation reports three low-lying (<1 eV) silylidenes [2-methylenesilabicyclo[1.1.0]but-1(3)-en-4-ylidene (3), 4-sila-2-methylenebicyclo[1.1.0]but-1(3)-en-4-ylidene (6), and 3-ethynyl-1-silapropadienylidene (7)] and three high-lying (>1 eV) silylidenes [2-sila-(didehydrovinylidene)cyclopropene (8), an isomer with a planar tetracoordinate carbon (ptC) atom (10), and 1-ethynyl-1-silapropadienylidene (11)], which remain elusive in the laboratory to date. Isomer 9 also contains a ptC atom, which turned out to be a transition state at all levels. Though all isomers are polar (μ ≠ 0), rotational spectrum is available only for 4. Using matrix isolation, three isomers (1, 2, and 5) have been trapped in the laboratory at 10 K. Considering the astrochemical relevance of silicon-carbide clusters in the interstellar medium, the current theoretical data demand new molecular spectroscopic studies on SiC4H2. Surprisingly, unlike the isovalent C5H2 isomers, where the bent carbenes are yet to be identified in the laboratory, the bent silylidenes (2 and 5) have been trapped in the case of SiC4H2. In both the cases, molecules with transannular C-C and/or Si-C bonds remain elusive, though they lie in the low-lying region. Using suitable precursors, whether these peculiar geometries (especially 3 and 6) would be identified or not in the laboratory needs to be addressed by molecular spectroscopists. The present investigation documents structural and spectroscopic information of SiC4H2 isomers, which may compliment future molecular spectroscopic observations including radioastronomical searches.
Collapse
Affiliation(s)
- Nisha Job
- Department of Chemistry, School of Advanced Sciences , Vellore Institute of Technology , Vellore 632014 , Tamil Nadu , India
| | - Amir Karton
- School of Molecular Sciences , The University of Western Australia , Perth , Western Australia 6009 , Australia
| | - Krishnan Thirumoorthy
- Department of Chemistry, School of Advanced Sciences , Vellore Institute of Technology , Vellore 632014 , Tamil Nadu , India
| | - Andrew L Cooksy
- Department of Chemistry and Biochemistry , San Diego State University , San Diego , California 92182-1030 , United States
| | - Venkatesan S Thimmakondu
- Department of Chemistry and Biochemistry , San Diego State University , San Diego , California 92182-1030 , United States
| |
Collapse
|
98
|
Annaberdiyev A, Melton CA, Bennett MC, Wang G, Mitas L. Accurate Atomic Correlation and Total Energies for Correlation Consistent Effective Core Potentials. J Chem Theory Comput 2020; 16:1482-1502. [DOI: 10.1021/acs.jctc.9b00962] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Abdulgani Annaberdiyev
- Department of Physics, North Carolina State University, Raleigh, North Carolina 27695-8202, United States
| | - Cody A. Melton
- Department of Physics, North Carolina State University, Raleigh, North Carolina 27695-8202, United States
- Sandia National Laboratories, Albuquerque, New Mexico 87123, United States
| | - M. Chandler Bennett
- Department of Physics, North Carolina State University, Raleigh, North Carolina 27695-8202, United States
| | - Guangming Wang
- Department of Physics, North Carolina State University, Raleigh, North Carolina 27695-8202, United States
| | - Lubos Mitas
- Department of Physics, North Carolina State University, Raleigh, North Carolina 27695-8202, United States
| |
Collapse
|
99
|
Smith CD, Karton A. Kinetics and Thermodynamics of Reactions Involving Criegee Intermediates: An Assessment of Density Functional Theory and Ab Initio Methods Through Comparison with CCSDT(Q)/CBS Data. J Comput Chem 2020; 41:328-339. [PMID: 31750964 DOI: 10.1002/jcc.26106] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2019] [Revised: 10/21/2019] [Accepted: 10/22/2019] [Indexed: 12/21/2022]
Abstract
Reactions involving Criegee intermediates (CIs, R1 R2 COO) are important in atmospheric ozonolysis models. In recent years, density functional theory (DFT) and CCSD(T)-based ab initio methods are increasingly being used for modeling reaction profiles involving CIs. We obtain highly accurate CCSDT(Q)/CBS reaction energies and barrier heights for ring-closing reactions involving atmospherically important CIs (R1 /R2 = H, Me, OH, OMe, F, CN, cyclopropene, ethylene, acetaldehyde, and acrolein). We use this benchmark data to evaluate the performance of DFT, double-hybrid DFT (DHDFT), and ab initio methods for the kinetics and thermodynamics of these reactions. We find that reaction energies are more challenging for approximate theoretical procedures than barrier heights. Overall, taking both reaction energies and barrier heights into account, only one of the 58 considered DFT methods (the meta-GGA MN12-L) attains near chemical accuracy, with root-mean-square deviations (RMSDs) of 3.5 (barrier heights) and 4.7 (reaction energies) kJ mol-1 . Therefore, MN12-L is recommended for investigations where CCSD(T)-based methods are not computationally feasible. For reaction barrier heights performance does not strictly follow Jacob's Ladder, for example, DHDFT methods do not perform better than conventional DFT methods. Of the ab initio methods, the cost-effective CCSD(T)/CBS(MP2) approach gives the best performance for both reaction energies and barrier heights, with RMSDs of 1.7 and 1.4 kJ mol-1 , respectively. All the considered Gaussian-n methods show good performance with RMSDs below the threshold of chemical accuracy for both reaction energies and barrier heights, where G4(MP2) shows the best overall performance with RMSDs of 2.9 and 1.5 kJ mol-1 , respectively. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Cameron D Smith
- School of Molecular Sciences, The University of Western Australia, Perth, Western Australia, 6009, Australia
| | - Amir Karton
- School of Molecular Sciences, The University of Western Australia, Perth, Western Australia, 6009, Australia
| |
Collapse
|
100
|
Zhang T, Wen M, Zeng Z, Lu Y, Wang Y, Wang W, Shao X, Wang Z, Makroni L. Effect of NH 3 and HCOOH on the H 2O 2 + HO → HO 2 + H 2O reaction in the troposphere: competition between the one-step and stepwise mechanisms. RSC Adv 2020; 10:9093-9102. [PMID: 35496523 PMCID: PMC9050117 DOI: 10.1039/d0ra00024h] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/02/2020] [Accepted: 02/13/2020] [Indexed: 11/21/2022] Open
Abstract
The H2O2 + HO → HO2 + H2O reaction is an important reservoir for both radicals of HO and HO2 catalyzing the destruction of O3. Here, this reaction assisted by NH3 and HCOOH catalysts was explored using the CCSD(T)-F12a/cc-pVDZ-F12//M06-2X/aug-cc-pVTZ method and canonical variational transition state theory with small curvature tunneling. Two possible sets of mechanisms, (i) one-step routes and (ii) stepwise processes, are possible. Our results show that in the presence of both NH3 and HCOOH catalysts under relevant atmospheric temperature, mechanism (i) is favored both energetically and kinetically than the corresponding mechanism (ii). At 298 K, the relative rate for mechanism (i) in the presence of NH3 (10, 2900 ppbv) and HCOOH (10 ppbv) is respectively 3–5 and 2–4 orders of magnitude lower than that of the water-catalyzed reaction. This is due to a comparatively lower concentration of NH3 and HCOOH than H2O which indicates the positive water effect under atmospheric conditions. Although NH3 and HCOOH catalysts play a negligible role in the reservoir for both radicals of HO and HO2 catalyzing the destruction of O3, the current study provides a comprehensive example of how acidic and basic catalysts assisted the gas-phase reactions. The H2O2 + HO → HO2 + H2O reaction is an important reservoir for both radicals of HO and HO2 catalyzing the destruction of O3.![]()
Collapse
Affiliation(s)
- Tianlei Zhang
- Shaanxi Key Laboratory of Catalysis
- School of Chemical & Environment Science
- Shaanxi University of Technology
- Hanzhong
- P. R. China
| | - Mingjie Wen
- Shaanxi Key Laboratory of Catalysis
- School of Chemical & Environment Science
- Shaanxi University of Technology
- Hanzhong
- P. R. China
| | - Zhaopeng Zeng
- Shaanxi Key Laboratory of Catalysis
- School of Chemical & Environment Science
- Shaanxi University of Technology
- Hanzhong
- P. R. China
| | - Yousong Lu
- Key Laboratory for Macromolecular Science of Shaanxi Province
- School of Chemistry & Chemical Engineering
- Shaanxi Normal University
- Xi'an
- P. R. China
| | - Yan Wang
- Shaanxi Key Laboratory of Catalysis
- School of Chemical & Environment Science
- Shaanxi University of Technology
- Hanzhong
- P. R. China
| | - Wei Wang
- Shaanxi Key Laboratory of Catalysis
- School of Chemical & Environment Science
- Shaanxi University of Technology
- Hanzhong
- P. R. China
| | - Xianzhao Shao
- Shaanxi Key Laboratory of Catalysis
- School of Chemical & Environment Science
- Shaanxi University of Technology
- Hanzhong
- P. R. China
| | - Zhiyin Wang
- Shaanxi Key Laboratory of Catalysis
- School of Chemical & Environment Science
- Shaanxi University of Technology
- Hanzhong
- P. R. China
| | - Lily Makroni
- Key Laboratory for Macromolecular Science of Shaanxi Province
- School of Chemistry & Chemical Engineering
- Shaanxi Normal University
- Xi'an
- P. R. China
| |
Collapse
|