1
|
Thorpe JH, Windom ZW, Bartlett RJ, Matthews DA. Factorized Quadruples and a Predictor of Higher-Level Correlation in Thermochemistry. J Phys Chem A 2024. [PMID: 39194317 DOI: 10.1021/acs.jpca.4c04460] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/29/2024]
Abstract
Coupled cluster theory has had a momentous impact on the ab initio prediction of molecular properties, and remains a staple ingratiate in high-accuracy thermochemical model chemistries. However, these methods require inclusion of at least some connected quadruple excitations, which generally scale at best as O ( N 9 ) with the number of basis functions. It is very difficult to predict, a priori, the effect correlation of past CCSD(T) on a given reaction energy. The purpose of this work is to examine cost-effective quadruple corrections based on the factorization theorem of the many-body perturbation theory that may address these challenges. We show that the O ( N 7 ) factorized CCSD(TQf) method introduces minimal error to predicted correlation and reaction energies as compared to the O ( N 9 ) CCSD(TQ). Further, we examine the performance of Goodson's continued fraction method in the estimation of CCSDT(Q)Λ contributions to reaction energies as well as a "new" method related to %TAE[(T)] that we refer to as a scaled perturbation estimator. We find that the scaled perturbation estimator based upon CCSD(TQf)/cc-pVDZ is capable of predicting CCSDT(Q)Λ/cc-pVDZ contributions to reaction energies with an average error of 0.07 kcal mol-1 and an L2D of 0.52 kcal mol-1 when applied to a test-suite of nearly 3000 reactions. This offers a means by which to reliably "ballpark" how important post-CCSD(T) contributions are to reaction energies while incurring no more than CCSD(T) formal cost and a little mental math.
Collapse
Affiliation(s)
- James H Thorpe
- Department of Chemistry, Southern Methodist University, Dallas, Texas 75275, United States
| | - Zachary W Windom
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, United States
| | - Rodney J Bartlett
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, United States
| | - Devin A Matthews
- Department of Chemistry, Southern Methodist University, Dallas, Texas 75275, United States
| |
Collapse
|
2
|
Schröder B, Rauhut G. From the Automated Calculation of Potential Energy Surfaces to Accurate Infrared Spectra. J Phys Chem Lett 2024; 15:3159-3169. [PMID: 38478898 PMCID: PMC10961845 DOI: 10.1021/acs.jpclett.4c00186] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2024] [Revised: 02/20/2024] [Accepted: 02/28/2024] [Indexed: 03/22/2024]
Abstract
Advances in the development of quantum chemical methods and progress in multicore architectures in computer science made the simulation of infrared spectra of isolated molecules competitive with respect to established experimental methods. Although it is mainly the multidimensional potential energy surface that controls the accuracy of these calculations, the subsequent vibrational structure calculations need to be carefully converged in order to yield accurate results. As both aspects need to be considered in a balanced way, we focus on approaches for molecules of up to 12-15 atoms with respect to both parts, which have been automated to some extent so that they can be employed in routine applications. Alternatives to machine learning will be discussed, which appear to be attractive, as long as local regions of the potential energy surface are sufficient. The automatization of these methods is still in its infancy, and the generalization to molecules with large amplitude motions or molecular clusters is far from trivial, but many systems relevant for astrophysical studies are already in reach.
Collapse
Affiliation(s)
- Benjamin Schröder
- Institute
of Physical Chemistry, University of Goettingen, Tammannstrasse 6, Göttingen 37077, Germany
| | - Guntram Rauhut
- Institute
for Theoretical Chemistry, University of
Stuttgart, Pfaffenwaldring 55, Stuttgart 70569, Germany
| |
Collapse
|
3
|
Vysotskiy VP, Filippi C, Ryde U. Scalar Relativistic All-Electron and Pseudopotential Ab Initio Study of a Minimal Nitrogenase [Fe(SH) 4H] - Model Employing Coupled-Cluster and Auxiliary-Field Quantum Monte Carlo Many-Body Methods. J Phys Chem A 2024; 128:1358-1374. [PMID: 38324717 PMCID: PMC10895656 DOI: 10.1021/acs.jpca.3c05808] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2023] [Revised: 01/03/2024] [Accepted: 01/03/2024] [Indexed: 02/09/2024]
Abstract
Nitrogenase is the only enzyme that can cleave the triple bond in N2, making nitrogen available to organisms. The detailed mechanism of this enzyme is currently not known, and computational studies are complicated by the fact that different density functional theory (DFT) methods give very different energetic results for calculations involving nitrogenase models. Recently, we designed a [Fe(SH)4H]- model with the fifth proton binding either to Fe or S to mimic different possible protonation states of the nitrogenase active site. We showed that the energy difference between these two isomers (ΔE) is hard to estimate with quantum-mechanical methods. Based on nonrelativistic single-reference coupled-cluster (CC) calculations, we estimated that the ΔE is 101 kJ/mol. In this study, we demonstrate that scalar relativistic effects play an important role and significantly affect ΔE. Our best revised single-reference CC estimates for ΔE are 85-91 kJ/mol, including energy corrections to account for contributions beyond triples, core-valence correlation, and basis-set incompleteness error. Among coupled-cluster approaches with approximate triples, the canonical CCSD(T) exhibits the largest error for this problem. Complementary to CC, we also used phaseless auxiliary-field quantum Monte Carlo calculations (ph-AFQMC). We show that with a Hartree-Fock (HF) trial wave function, ph-AFQMC reproduces the CC results within 5 ± 1 kJ/mol. With multi-Slater-determinant (MSD) trials, the results are 82-84 ± 2 kJ/mol, indicating that multireference effects may be rather modest. Among the DFT methods tested, τ-HCTH, r2SCAN with 10-13% HF exchange with and without dispersion, and O3LYP/O3LYP-D4, and B3LYP*/B3LYP*-D4 generally perform the best. The r2SCAN12 (with 12% HF exchange) functional mimics both the best reference MSD ph-AFQMC and CC ΔE results within 2 kJ/mol.
Collapse
Affiliation(s)
- Victor P. Vysotskiy
- Department
of Computational Chemistry, Lund University,
Chemical Centre, SE-221 00 Lund, Sweden
| | - Claudia Filippi
- MESA+
Institute for Nanotechnology, University
of Twente, P.O. Box 217, Enschede 7500 AE, Netherlands
| | - Ulf Ryde
- Department
of Computational Chemistry, Lund University,
Chemical Centre, SE-221 00 Lund, Sweden
| |
Collapse
|
4
|
Vysotskiy VP, Torbjörnsson M, Jiang H, Larsson ED, Cao L, Ryde U, Zhai H, Lee S, Chan GKL. Assessment of DFT functionals for a minimal nitrogenase [Fe(SH)4H]- model employing state-of-the-art ab initio methods. J Chem Phys 2023; 159:044106. [PMID: 37486046 DOI: 10.1063/5.0152611] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2023] [Accepted: 07/03/2023] [Indexed: 07/25/2023] Open
Abstract
We have designed a [Fe(SH)4H]- model with the fifth proton binding either to Fe or S. We show that the energy difference between these two isomers (∆E) is hard to estimate with quantum-mechanical (QM) methods. For example, different density functional theory (DFT) methods give ∆E estimates that vary by almost 140 kJ/mol, mainly depending on the amount of exact Hartree-Fock included (0%-54%). The model is so small that it can be treated by many high-level QM methods, including coupled-cluster (CC) and multiconfigurational perturbation theory approaches. With extrapolated CC series (up to fully connected coupled-cluster calculations with singles, doubles, and triples) and semistochastic heat-bath configuration interaction methods, we obtain results that seem to be converged to full configuration interaction results within 5 kJ/mol. Our best result for ∆E is 101 kJ/mol. With this reference, we show that M06 and B3LYP-D3 give the best results among 35 DFT methods tested for this system. Brueckner doubles coupled cluster with perturbaitve triples seems to be the most accurate coupled-cluster approach with approximate triples. CCSD(T) with Kohn-Sham orbitals gives results within 4-11 kJ/mol of the extrapolated CC results, depending on the DFT method. Single-reference CC calculations seem to be reasonably accurate (giving an error of ∼5 kJ/mol compared to multireference methods), even if the D1 diagnostic is quite high (0.25) for one of the two isomers.
Collapse
Affiliation(s)
- Victor P Vysotskiy
- Department of Computational Chemistry, Lund University, Chemical Centre, SE-221 00 Lund, Sweden
| | - Magne Torbjörnsson
- Department of Computational Chemistry, Lund University, Chemical Centre, SE-221 00 Lund, Sweden
| | - Hao Jiang
- Department of Computational Chemistry, Lund University, Chemical Centre, SE-221 00 Lund, Sweden
| | - Ernst D Larsson
- Department of Computational Chemistry, Lund University, Chemical Centre, SE-221 00 Lund, Sweden
| | - Lili Cao
- Department of Computational Chemistry, Lund University, Chemical Centre, SE-221 00 Lund, Sweden
| | - Ulf Ryde
- Department of Computational Chemistry, Lund University, Chemical Centre, SE-221 00 Lund, Sweden
| | - Huanchen Zhai
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Seunghoon Lee
- 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
| |
Collapse
|
5
|
Rocha CMR, Linnartz H. High-level ab initio quartic force fields and spectroscopic characterization of C 2N . Phys Chem Chem Phys 2021; 23:26227-26240. [PMID: 34787132 DOI: 10.1039/d1cp03505c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
While it is now well established that large carbon chain species and radiative electron attachment (REA) are key ingredients triggering interstellar anion chemistry, the role played by smaller molecular anions, for which REA appears to be an unlikely formation pathway, is as yet elusive. Advancing this research undoubtedly requires the knowledge (and modeling) of their astronomical abundances which, for the case of C2N-, is largely hindered by a lack of accurate spectroscopic signatures. In this work, we provide such data for both ground -CCN-(3Σ-) and low-lying c-CNC-(1A1) isomers and their singly-substituted isotopologues by means of state-of-the-art rovibrational quantum chemical techniques. Their quartic force fields are herein calibrated using a high-level composite energy scheme that accounts for extrapolations to both one-particle and (approximate) -particle basis set limits, in addition to relativistic effects, with the final forms being subsequently subject to nuclear motion calculations. Besides standard spectroscopic attributes, the full set of computed properties includes fine and hyperfine interaction constants and can be readily introduced as guesses in conventional experimental data reduction analyses through effective Hamiltonians. On the basis of benchmark calculations performed anew for a minimal test set of prototypical triatomics and limited (low-resolution) experimental data for -CCN-(3Σ-), the target accuracies are determined to be better than 0.1% of experiment for rotational constants and 0.3% for vibrational fundamentals. Apart from laboratory investigations, the results here presented are expected to also prompt future astronomical surveys on C2N-. To this end and using the theoretically-predicted spectroscopic constants, the rotational spectra of both -CCN-(3Σ-) and c-CNC-(1A1) are derived and their likely detectability in the interstellar medium is further explored in connection with working frequency ranges of powerful astronomical facilities. Our best theoretical estimate places c-CNC-(1A1) at about 15.3 kcal mol-1 above the ground-state -CCN-(3Σ-) species.
Collapse
Affiliation(s)
- C M R Rocha
- Laboratory for Astrophysics, Leiden Observatory, Leiden University, P.O. Box 9513, NL-2300 RA Leiden, The Netherlands.
| | - H Linnartz
- Laboratory for Astrophysics, Leiden Observatory, Leiden University, P.O. Box 9513, NL-2300 RA Leiden, The Netherlands.
| |
Collapse
|
6
|
Taamalli S, Pitoňák M, Dibble TS, Černušák I, Louis F. Theoretical Study of the Monohydration of Mercury Compounds of Atmospheric Interest. J Phys Chem A 2021; 125:5819-5828. [PMID: 34180661 DOI: 10.1021/acs.jpca.1c02772] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The structures, vibrational frequencies, and model IR spectra of the monohydrates of oxygenated mercury compounds (BrHgO, BrHgOH, BrHgOOH, BrHgNO2, BrHgONO, and HgOH) have been theoretically studied using the ωB97X-D/aug-cc-pVTZ level of theory. The ground state potential energy surface exhibits several stable structures of these monohydrates. The thermodynamic properties of the hydration reactions have been calculated at different levels of theory including DFT and coupled-cluster calculations DK-CCSD(T) with the ANO-RCC-Large basis sets. Standard enthalpies and Gibbs free energies of hydration were computed. The temperature dependence of ΔrG°(T) was evaluated for the most stable complexes over the temperature range 200-400 K. Thermodynamic data revealed that the highest fraction hydrated at 298 K and 100% relative humidity will be BrHgNO2-H2O at ∼5%.
Collapse
Affiliation(s)
- Sonia Taamalli
- CNRS, UMR 8522 - PC2A - PhysicoChimie des Processus de Combustion et de l'Atmosphère, Université Lille, 59000 Lille, France
| | - Michal Pitoňák
- Department of Physical and Theoretical Chemistry, Faculty of Natural Sciences, Comenius University in Bratislava, Ilkovičova 6, 84215 Bratislava, Slovakia.,Computing Center, Centre of Operations of the Slovak Academy of Sciences, Dúbravská cesta 9, 845 35 Bratislava, Slovakia
| | - Theodore S Dibble
- Department of Chemistry, State University of New York-College of Environmental Science and Forestry, 1 Forestry Drive, Syracuse, New York 13210, United States
| | - Ivan Černušák
- Department of Physical and Theoretical Chemistry, Faculty of Natural Sciences, Comenius University in Bratislava, Ilkovičova 6, 84215 Bratislava, Slovakia
| | - Florent Louis
- CNRS, UMR 8522 - PC2A - PhysicoChimie des Processus de Combustion et de l'Atmosphère, Université Lille, 59000 Lille, France
| |
Collapse
|
7
|
Feller D, Bross DH, Ruscic B. Enthalpy of Formation of C2H2O4 (Oxalic Acid) from High-Level Calculations and the Active Thermochemical Tables Approach. J Phys Chem A 2019; 123:3481-3496. [DOI: 10.1021/acs.jpca.8b12329] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- David Feller
- Department of Chemistry, Washington State University, Pullman, Washington 99164-4630, United States
| | - David H. Bross
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - Branko Ruscic
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
- Consortium for Advanced Science and Engineering, The University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
8
|
Schulz CE, Dutta AK, Izsák R, Pantazis DA. Systematic High-Accuracy Prediction of Electron Affinities for Biological Quinones. J Comput Chem 2018; 39:2439-2451. [PMID: 30281169 DOI: 10.1002/jcc.25570] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2018] [Revised: 08/06/2018] [Accepted: 08/07/2018] [Indexed: 11/07/2022]
Abstract
Quinones play vital roles as electron carriers in fundamental biological processes; therefore, the ability to accurately predict their electron affinities is crucial for understanding their properties and function. The increasing availability of cost-effective implementations of correlated wave function methods for both closed-shell and open-shell systems offers an alternative to density functional theory approaches that have traditionally dominated the field despite their shortcomings. Here, we define a benchmark set of quinones with experimentally available electron affinities and evaluate a range of electronic structure methods, setting a target accuracy of 0.1 eV. Among wave function methods, we test various implementations of coupled cluster (CC) theory, including local pair natural orbital (LPNO) approaches to canonical and parameterized CCSD, the domain-based DLPNO approximation, and the equations-of-motion approach for electron affinities, EA-EOM-CCSD. In addition, several variants of canonical, spin-component-scaled, orbital-optimized, and explicitly correlated (F12) Møller-Plesset perturbation theory are benchmarked. Achieving systematically the target level of accuracy is challenging and a composite scheme that combines canonical CCSD(T) with large basis set LPNO-based extrapolation of correlation energy proves to be the most accurate approach. Methods that offer comparable performance are the parameterized LPNO-pCCSD, the DLPNO-CCSD(T0 ), and the orbital optimized OO-SCS-MP2. Among DFT methods, viable practical alternatives are only the M06 and the double hybrids, but the latter should be employed with caution because of significant basis set sensitivity. A highly accurate yet cost-effective DLPNO-based coupled cluster approach is used to investigate the methoxy conformation effect on the electron affinities of ubiquinones found in photosynthetic bacterial reaction centers. © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Christine E Schulz
- Fakultät für Chemie und Biochemie, Ruhr-Universität Bochum, 44780, Bochum, Germany
- Max-Planck-Institut für Chemische Energiekonversion, Stiftstr. 34-36, 45470, Mülheim an der Ruhr, Germany
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470, Mülheim an der Ruhr, Germany
| | - Achintya Kumar Dutta
- Max-Planck-Institut für Chemische Energiekonversion, Stiftstr. 34-36, 45470, Mülheim an der Ruhr, Germany
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai 400076, India
| | - Róbert Izsák
- Max-Planck-Institut für Chemische Energiekonversion, Stiftstr. 34-36, 45470, Mülheim an der Ruhr, Germany
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470, Mülheim an der Ruhr, Germany
| | - Dimitrios A Pantazis
- Max-Planck-Institut für Chemische Energiekonversion, Stiftstr. 34-36, 45470, Mülheim an der Ruhr, Germany
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470, Mülheim an der Ruhr, Germany
| |
Collapse
|
9
|
Hu L, Chen K, Chen H. Modeling σ-Bond Activations by Nickel(0) Beyond Common Approximations: How Accurately Can We Describe Closed-Shell Oxidative Addition Reactions Mediated by Low-Valent Late 3d Transition Metal? J Chem Theory Comput 2017; 13:4841-4853. [PMID: 28881134 DOI: 10.1021/acs.jctc.7b00708] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Accurate modelings of reactions involving 3d transition metals (TMs) are very challenging to both ab initio and DFT approaches. To gain more knowledge in this field, we herein explored typical σ-bond activations of H-H, C-H, C-Cl, and C-C bonds promoted by nickel(0), a low-valent late 3d TM. For the key parameters of activation energy (ΔE‡) and reaction energy (ΔER) for these reactions, various issues related to the computational accuracy were systematically investigated. From the scrutiny of convergence issue with one-electron basis set, augmented (A) basis functions are found to be important, and the CCSD(T)/CBS level with complete basis set (CBS) limit extrapolation based on augmented double-ζ and triple-ζ basis pair (ADZ and ATZ), which produces deviations below 1 kcal/mol from the reference, is recommended for larger systems. As an alternative, the explicitly correlated F12 method can accelerate the basis set convergence further, especially after its CBS extrapolations. Thus, the CCSD(T)-F12/CBS(ADZ-ATZ) level with computational cost comparable to the conventional CCSD(T)/CBS(ADZ-ATZ) level, is found to reach the accuracy of the conventional CCSD(T)/A5Z level, which produces deviations below 0.5 kcal/mol from the reference, and is also highly recommendable. Scalar relativistic effects and 3s3p core-valence correlation are non-negligible for achieving chemical accuracy of around 1 kcal/mol. From the scrutiny of convergence issue with the N-electron basis set, in comparison with the reference CCSDTQ result, CCSD(T) is found to be able to calculate ΔE‡ quite accurately, which is not true for the ΔER calculations. Using highest-level CCSD(T) results of ΔE‡ in this work as references, we tested 18 DFT methods and found that PBE0 and CAM-B3LYP are among the three best performing functionals, irrespective of DFT empirical dispersion correction. With empirical dispersion correction included, ωB97XD is also recommendable due to its improved performance.
Collapse
Affiliation(s)
- Lianrui Hu
- Beijing National Laboratory for Molecular Sciences (BNLMS), CAS Key Laboratory of Photochemistry, CAS Research/Education Center for Excellence in Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences , Beijing 100190, China
| | - Kejuan Chen
- Beijing National Laboratory for Molecular Sciences (BNLMS), CAS Key Laboratory of Photochemistry, CAS Research/Education Center for Excellence in Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences , Beijing 100190, China
| | - Hui Chen
- Beijing National Laboratory for Molecular Sciences (BNLMS), CAS Key Laboratory of Photochemistry, CAS Research/Education Center for Excellence in Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences , Beijing 100190, China
| |
Collapse
|
10
|
Feller D, Bross DH, Ruscic B. Enthalpy of Formation of N2H4 (Hydrazine) Revisited. J Phys Chem A 2017; 121:6187-6198. [DOI: 10.1021/acs.jpca.7b06017] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- David Feller
- Department
of Chemistry Washington State University, Pullman, Washington 99164-4630, United States
| | - David H. Bross
- Chemical
Sciences and Engineering Division Argonne National Laboratory Argonne, Illinois 60439, United States
| | - Branko Ruscic
- Chemical
Sciences and Engineering Division Argonne National Laboratory Argonne, Illinois 60439, United States
- Computation
Institute, The University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
11
|
Feller D. Estimating the intrinsic limit of the Feller-Peterson-Dixon composite approach when applied to adiabatic ionization potentials in atoms and small molecules. J Chem Phys 2017; 147:034103. [DOI: 10.1063/1.4993625] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
12
|
The complete basis set Full-CI roto-vibrational spectroscopic constants of AlH, $$\hbox {AlH}^{+}$$ AlH + and $$\hbox {AlH}^{-}$$ AlH -. Theor Chem Acc 2016. [DOI: 10.1007/s00214-016-2026-2] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
13
|
Khanniche S, Louis F, Cantrel L, Černušák I. A Density Functional Theory and ab Initio Investigation of the Oxidation Reaction of CO by IO Radicals. J Phys Chem A 2016; 120:1737-49. [PMID: 26908233 DOI: 10.1021/acs.jpca.6b00047] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
To get an insight into the possible reactivity between iodine oxides and CO, a first step was to study the thermochemical properties and kinetic parameters of the reaction between IO and CO using theoretical chemistry tools. All stationary points involved were optimized using the Becke's three-parameter hybrid exchange functional coupled with the Lee-Yang-Parr nonlocal correlation functional (B3LYP) and the Møller-Plesset second-order perturbation theory (MP2). Single-point energy calculations were performed using the coupled cluster theory with the iterative inclusion of singles and doubles and the perturbative estimation for triple excitations (CCSD(T)) and the aug-cc-pVnZ (n = T, Q, and 5) basis sets on geometries previously optimized at the aug-cc-pVTZ level. The energetics was then recalculated using the one-component DK-CCSD(T) approach with the relativistic ANO basis sets. The spin-orbit coupling for the iodine containing species was calculated a posteriori using the restricted active space state interaction method in conjunction with the multiconfigurational perturbation theory (CASPT2/RASSI) employing the complete active space (CASSCF) wave function as the reference. The CCSD(T) energies were also corrected for BSSE for molecular complexes and refined with the extrapolation to CBS limit while the DK-CCSD(T) values were refined with the extrapolation to FCI. The exploration of the potential energy surface revealed a two-steps mechanism with a trans and a cis pathway. The rate constants for the direct and complex mechanism were computed as a function of temperature (250-2500 K) using the canonical transition state theory. The three-parameter Arrhenius expressions obtained for the direct and indirect mechanism at the DK-CCSD(T)-cf level of theory is 1.49 × 10(-17) × T(1.77) exp(-47.4 (kJ mol(-1))/RT).
Collapse
Affiliation(s)
- Sarah Khanniche
- University Lille, CNRS, UMR 8522-PC2A, PhysicoChimie des Processus de Combustion et de l'Atmosphère, F-59000 Lille, France.,Laboratoire de Recherche Commun IRSN-CNRS-Lille1 "Cinétique Chimique, Combustion, Réactivité″ (C3R), Cadarache, St Paul Lez Durance, 13115, France
| | - Florent Louis
- University Lille, CNRS, UMR 8522-PC2A, PhysicoChimie des Processus de Combustion et de l'Atmosphère, F-59000 Lille, France.,Laboratoire de Recherche Commun IRSN-CNRS-Lille1 "Cinétique Chimique, Combustion, Réactivité″ (C3R), Cadarache, St Paul Lez Durance, 13115, France
| | - Laurent Cantrel
- Institut de Radioprotection et de Sûreté Nucléaire (IRSN), PSN-RES, Cadarache, St Paul Lez Durance, 13115, France.,Laboratoire de Recherche Commun IRSN-CNRS-Lille1 "Cinétique Chimique, Combustion, Réactivité″ (C3R), Cadarache, St Paul Lez Durance, 13115, France
| | - Ivan Černušák
- Department of Physical and Theoretical Chemistry, Faculty of Natural Sciences, Comenius University in Bratislava , Mlynská dolina CH1, 84215 Bratislava, Slovakia
| |
Collapse
|
14
|
Šulková K, Cantrel L, Louis F. Gas-Phase Reactivity of Cesium-Containing Species by Quantum Chemistry. J Phys Chem A 2015; 119:9373-84. [PMID: 26237575 DOI: 10.1021/acs.jpca.5b05548] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Thermodynamics and kinetics of cesium species reactions have been studied by using high-level quantum chemical tools. A systematic theoretical study has been done to find suitable methodology for calculation of reliable thermodynamic properties, allowing us to determine bimolecular rate constants with appropriate kinetic theories of gas-phase reactions. Four different reactions have been studied in this work: CsO + H2 = CsOH + H (R1), Cs + HI = CsI + H (R2), CsI + H2O = CsOH + HI (R3), and CsI + OH = CsOH + I (R4). All reactions involve steam, hydrogen, and iodine in addition of cesium. Most of the reactions are fast and (R3) and (R4) proceed even without energetic barrier. In terms of chemical reactivity in the reactor coolant system (RCS) in the case of severe accident, it can be expected that there will be no kinetic limitations for main cesium species (CsOH and CsI) transported along the RCS. Cs chemical speciation inside the RCS should be governed by the thermodynamics.
Collapse
Affiliation(s)
- Katarína Šulková
- Institut de Radioprotection et de Sûreté Nucléaire (IRSN), PSN-RES, SAG, LETR, St Paul Lez Durance 13115, France
| | - Laurent Cantrel
- Institut de Radioprotection et de Sûreté Nucléaire (IRSN), PSN-RES, SAG, LETR, St Paul Lez Durance 13115, France.,Laboratoire de Recherche Commun IRSN-CNRS-Lille1 "Cinétique Chimique, Combustion, Réactivité" (C3R), Centre de Cadarache, BP3, St Paul Lez Durance 13115, France
| | - Florent Louis
- PhysicoChimie des Processus de Combustion et de l'Atmosphère (PC2A), UMR CNRS 8522, Université Lille 1 Sciences et Technologies, 59655 Villeneuve d'Ascq Cedex, France.,Laboratoire de Recherche Commun IRSN-CNRS-Lille1 "Cinétique Chimique, Combustion, Réactivité" (C3R), Centre de Cadarache, BP3, St Paul Lez Durance 13115, France
| |
Collapse
|
15
|
Cao J, Wang ZX, Gao LJ, Fu F. Water-assisted isomerization of the [H, C, N, O] system. J Mol Model 2015; 21:66. [DOI: 10.1007/s00894-015-2628-y] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2014] [Accepted: 02/15/2015] [Indexed: 10/23/2022]
|
16
|
Louis F. A Theoretical Study of the Kinetics of the Hydrogen Atom Abstraction Reactions from Cyclopropane by H, O (3
P), and Cl (2
P3/2
) Atoms and OH Radicals. INT J CHEM KINET 2015. [DOI: 10.1002/kin.20905] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Affiliation(s)
- Florent Louis
- PhysicoChimie des Processus de Combustion et de l'Atmosphère (PC2A); UMR 8522 CNRS/Lille1, Université Lille 1 Sciences et Technologies; 59655 Villeneuve d'Ascq Cedex France
| |
Collapse
|
17
|
Sudolská M, Louis F, Cernušák I. Reactivity of CHI3 with OH radicals: X-abstraction reaction pathways (X = H, I), atmospheric chemistry, and nuclear safety. J Phys Chem A 2014; 118:9512-20. [PMID: 25207959 DOI: 10.1021/jp5051832] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The X-abstraction (X = H, I) pathways in the reaction of CHI3 with OH radical, a possible iodoform removal process relevant to the Earth's atmosphere and conditions prevailing in the case of a nuclear accident, have been studied applying highly correlated ab initio quantum chemistry methods and canonical transition-state theory to obtain reaction energy profiles and rate constants. Geometry optimizations of reactants, products, molecular complexes, and transition states determined at the MP2/cc-pVTZ level of theory have been followed by DK-CCSD(T)/ANO-RCC single-point energy calculations. Further improvement of electronic energies has been achieved by applying spin-orbit coupling, corrections toward full configuration interaction, vibration contributions, and tunneling corrections. Calculated reaction enthalpies at 0 K are -108.2 and -5.1 kJ mol(-1) for the H- and I-abstraction pathways, respectively; the strongly exothermic H-abstraction pathway is energetically favored over the modestly exothermic I-abstraction one. The overall rate constant at 298 K based on our ab initio calculations is 4.90 × 10(-11) cm(3) molecule(-1) s(-1), with the I-abstraction pathway being the major channel over the temperature range of 250-2000 K. The CHI3 atmospheric lifetime with respect to the removal reaction with OH radical is predicted to be about 6 h, very short compared to that of other halomethanes.
Collapse
Affiliation(s)
- Mária Sudolská
- Department of Physical and Theoretical Chemistry, Faculty of Natural Sciences, Comenius University , Mlynská dolina CH1, 84215 Bratislava, Slovakia
| | | | | |
Collapse
|
18
|
|
19
|
Swart M, Bickelhaupt FM. Benchmark study on the smallest bimolecular nucleophilic substitution reaction: H⁻+CH₄-->CH₄+H⁻. Molecules 2013; 18:7726-38. [PMID: 23823873 PMCID: PMC6270058 DOI: 10.3390/molecules18077726] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2013] [Revised: 06/03/2013] [Accepted: 06/28/2013] [Indexed: 11/16/2022] Open
Abstract
We report here a benchmark study on the bimolecular nucleophilic substitution (S(N)2) reaction between hydride and methane, for which we have obtained reference energies at the coupled cluster toward full configuration-interaction limit (CC-cf/CBS). Several wavefunction (HF, MP2, coupled cluster) and density functional methods are compared for their reliability regarding these reference data.
Collapse
Affiliation(s)
- Marcel Swart
- Institució Catalana de Recerca i Estudis Avançats (ICREA), Pg. Lluís Companys 23, 08010 Barcelona, Spain
- Institut de Química Computacional i Catàlisi (IQCC) and Departament de Química, Universitat de Girona, Campus Montilivi, 17071 Girona, Spain
| | - F. Matthias Bickelhaupt
- Department of Theoretical Chemistry & Amsterdam Center for Multiscale Modeling, VU University, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands; E-Mail:
- Institute of Molecules and Materials, Radboud University Nijmegen, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands
| |
Collapse
|
20
|
Šulka M, Šulková K, Louis F, Neogrády P, Černušák I. A Theoretical Study of the X-Abstraction Reactions (X = H, Br, or I) from CH2IBr by OH Radicals: Implications for Atmospheric Chemistry. Z PHYS CHEM 2013. [DOI: 10.1524/zpch.2013.0391] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
21
|
Abstract
By analysing the properties of the electron density in the structurally simple perhalogenated ethanes, X3C-CY3 (X, Y = F, Cl), a previously overlooked non-covalent attraction between halogens attached to opposite carbon atoms is found. Quantum chemical calculations extrapolated towards the full solution of the Schrödinger equation reveal the complex nature of the interaction. When at least one of the halogens is a chlorine, the strength of the interaction is comparable to that of hydrogen bonds. Further analysis shows that the bond character is quite different from standard non-covalent halogen bonds and hydrogen bonds; no bond critical points are found between the halogens, and the σ-holes of the halogens are not utilised for bonding. Thus, the nature of the intramolecular halogen···halogen bonding studied here appears to be of an unusually strong van der Waals type.
Collapse
Affiliation(s)
- Mikael P Johansson
- Institut de Química Computational i Catàlisi and Departament de Química, Universitat de Girona, Campus Montilivi, ES-17071 Girona, Spain
| | | |
Collapse
|
22
|
Guerra M, Pinto R, Santos J, Paiva A. Towards the assignment of the REMPI spectrum of Ph2O using CIS and TD-DFT methods. Mol Phys 2013. [DOI: 10.1080/00268976.2013.783638] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
23
|
Šulková K, Šulka M, Louis F, Neogrády P. Atmospheric Reactivity of CH2ICl with OH Radicals: High-Level OVOS CCSD(T) Calculations for the X-Abstraction Pathways (X = H, Cl, or I). J Phys Chem A 2013; 117:771-82. [PMID: 23294186 DOI: 10.1021/jp312447x] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Affiliation(s)
- Katarína Šulková
- Department of Physical and Theoretical
Chemistry, Faculty of Natural Sciences, Comenius University, Mlynská dolina CH1, 84215 Bratislava, Slovakia
| | - Martin Šulka
- Department of Physical and Theoretical
Chemistry, Faculty of Natural Sciences, Comenius University, Mlynská dolina CH1, 84215 Bratislava, Slovakia
| | - Florent Louis
- PhysicoChimie des Processus
de Combustion et de l’Atmosphère (PC2A), UMR 8522 CNRS/Lille1, Université Lille 1 Sciences et Technologies,
Cité Scientifique, Bât. C11/C5, 59655 Villeneuve d’Ascq
Cedex, France
| | - Pavel Neogrády
- Department of Physical and Theoretical
Chemistry, Faculty of Natural Sciences, Comenius University, Mlynská dolina CH1, 84215 Bratislava, Slovakia
| |
Collapse
|
24
|
Feller D, Peterson KA, Dixon DA. Further benchmarks of a composite, convergent, statistically calibrated coupled-cluster-based approach for thermochemical and spectroscopic studies. Mol Phys 2012. [DOI: 10.1080/00268976.2012.684897] [Citation(s) in RCA: 122] [Impact Index Per Article: 10.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Affiliation(s)
- David Feller
- a Department of Chemistry , Washington State University , Pullman , Washington 99164-4630 , USA
| | - Kirk A. Peterson
- a Department of Chemistry , Washington State University , Pullman , Washington 99164-4630 , USA
| | - David A. Dixon
- b Chemistry Department , The University of Alabama , Shelby Hall, Tuscaloosa , Alabama 35487-0336 , USA
| |
Collapse
|
25
|
Dixon DA, Feller D, Peterson KA. A Practical Guide to Reliable First Principles Computational Thermochemistry Predictions Across the Periodic Table. ANNUAL REPORTS IN COMPUTATIONAL CHEMISTRY VOLUME 8 2012. [DOI: 10.1016/b978-0-444-59440-2.00001-6] [Citation(s) in RCA: 80] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
26
|
Feller D, Vasiliu M, Grant DJ, Dixon DA. Thermodynamic Properties of Arsenic Compounds and the Heat of Formation of the As Atom from High Level Electronic Structure Calculations. J Phys Chem A 2011; 115:14667-76. [DOI: 10.1021/jp209097x] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- David Feller
- Department of Chemistry, Washington State University, Pullman, Washington 99164-4630, United States
| | - Monica Vasiliu
- Department of Chemistry, The University of Alabama, Shelby Hall, Tuscaloosa, Alabama 35487-0336, United States
| | - Daniel J. Grant
- Department of Chemistry, The University of Alabama, Shelby Hall, Tuscaloosa, Alabama 35487-0336, United States
| | - David A. Dixon
- Department of Chemistry, The University of Alabama, Shelby Hall, Tuscaloosa, Alabama 35487-0336, United States
| |
Collapse
|
27
|
Mečiarová K, Šulka M, Canneaux S, Louis F, Černušák I. A theoretical study of the kinetics of the forward and reverse reactions HI+CH3=I+CH4. Chem Phys Lett 2011. [DOI: 10.1016/j.cplett.2011.10.029] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/16/2022]
|
28
|
Goodson DZ. Resummation methods. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2011. [DOI: 10.1002/wcms.92] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
|
29
|
Grant Hill J, Mitrushchenkov A, Yousaf KE, Peterson KA. Accurate ab initio ro-vibronic spectroscopy of the X̃2Π CCN radical using explicitly correlated methods. J Chem Phys 2011; 135:144309. [DOI: 10.1063/1.3647311] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
|
30
|
Ojeda OU, Çağın T. Hydrogen Bonding and Molecular Rearrangement in 1,3,5-Triamino-2,4,6-trinitrobenzene under Compression. J Phys Chem B 2011; 115:12085-93. [DOI: 10.1021/jp2007649] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Affiliation(s)
- Oscar U. Ojeda
- Laboratory of Computational Engineering of Nanomaterials, Artie McFerrin Department of Chemical Engineering, Texas A&M University, College Station, Texas 77843-3122, United States
| | - Tahir Çağın
- Laboratory of Computational Engineering of Nanomaterials, Artie McFerrin Department of Chemical Engineering, Texas A&M University, College Station, Texas 77843-3122, United States
| |
Collapse
|
31
|
Sancho-García JC. Stability of Hydrocarbons of the Polyhedrane Family: Convergence of ab Initio Calculations and Corresponding Assessment of DFT Main Approximations. J Chem Theory Comput 2011; 7:2761-5. [DOI: 10.1021/ct200198z] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Affiliation(s)
- J. C. Sancho-García
- Departamento de Química Física, University of Alicante, E-03080 Alicante, Spain
| |
Collapse
|
32
|
Atmospheric reactivity of CH3I and CH2I2 with OH radicals: A comparative study of the H- versus I-abstraction. COMPUT THEOR CHEM 2011. [DOI: 10.1016/j.theochem.2010.09.022] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
33
|
Feller D, Peterson KA, Dixon DA. Ab Initio Coupled Cluster Determination of the Heats of Formation of C2H2F2, C2F2, and C2F4. J Phys Chem A 2011; 115:1440-51. [DOI: 10.1021/jp111644h] [Citation(s) in RCA: 37] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- David Feller
- Department of Chemistry, Washington State University, Pullman, Washington 99164-4630, United States
| | - Kirk A. Peterson
- Department of Chemistry, Washington State University, Pullman, Washington 99164-4630, United States
| | - David A. Dixon
- Department of Chemistry, The University of Alabama, Shelby Hall, Box 870336, Tuscaloosa, Alabama 35487-0336, United States
| |
Collapse
|
34
|
Vasiliu M, Feller D, Gole JL, Dixon DA. Structures and Heats of Formation of Simple Alkaline Earth Metal Compounds: Fluorides, Chlorides, Oxides, and Hydroxides for Be, Mg, and Ca. J Phys Chem A 2010; 114:9349-58. [DOI: 10.1021/jp1050657] [Citation(s) in RCA: 39] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Monica Vasiliu
- Chemistry Department, Shelby Hall, The University of Alabama, Shelby Hall, Box 870336, Tuscaloosa, Alabama 35487-0336, Department of Chemistry, Washington State University, Pullman, Washington 99164-4630, and Schools of Physics and Mechanical Engineering, Georgia Institute of Technology, 837 State Street, NW, Atlanta, Georgia 30332-0430
| | - David Feller
- Chemistry Department, Shelby Hall, The University of Alabama, Shelby Hall, Box 870336, Tuscaloosa, Alabama 35487-0336, Department of Chemistry, Washington State University, Pullman, Washington 99164-4630, and Schools of Physics and Mechanical Engineering, Georgia Institute of Technology, 837 State Street, NW, Atlanta, Georgia 30332-0430
| | - James L. Gole
- Chemistry Department, Shelby Hall, The University of Alabama, Shelby Hall, Box 870336, Tuscaloosa, Alabama 35487-0336, Department of Chemistry, Washington State University, Pullman, Washington 99164-4630, and Schools of Physics and Mechanical Engineering, Georgia Institute of Technology, 837 State Street, NW, Atlanta, Georgia 30332-0430
| | - David A. Dixon
- Chemistry Department, Shelby Hall, The University of Alabama, Shelby Hall, Box 870336, Tuscaloosa, Alabama 35487-0336, Department of Chemistry, Washington State University, Pullman, Washington 99164-4630, and Schools of Physics and Mechanical Engineering, Georgia Institute of Technology, 837 State Street, NW, Atlanta, Georgia 30332-0430
| |
Collapse
|
35
|
Vasiliu M, Li S, Peterson KA, Feller D, Gole JL, Dixon DA. Structures and heats of formation of simple alkali metal compounds: hydrides, chlorides, fluorides, hydroxides, and oxides for Li, Na, and K. J Phys Chem A 2010; 114:4272-81. [PMID: 20201583 DOI: 10.1021/jp911735c] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
Geometry parameters, frequencies, heats of formation, and bond dissociation energies are predicted for simple alkali metal compounds (hydrides, chlorides, fluorides, hydroxides and oxides) of Li, Na, and K from coupled cluster theory [CCSD(T)] calculations including core-valence correlation with the aug-cc-pwCVnZ basis set (n = D, T, Q, and 5). To accurately calculate the heats of formation, the following additional correction were included: scalar relativistic effects, atomic spin-orbit effects, and vibrational zero-point energies. For calibration purposes, the properties of some of the lithium compounds were predicted with iterative triple and quadruple excitations via CCSDT and CCSDTQ. The calculated geometry parameters, frequencies, heats of formation, and bond dissociation energies were compared with all available experimental measurements and are in excellent agreement with high-quality experimental data. High-level calculations are required to correctly predict that K(2)O is linear and that the ground state of KO is (2)Sigma(+), not (2)Pi, as in LiO and NaO. This reliable and consistent set of calculated thermodynamic data is appropriate for use in combustion and atmospheric simulations.
Collapse
Affiliation(s)
- Monica Vasiliu
- Chemistry Department, Shelby Hall, The University of Alabama, Box 870336, Tuscaloosa, Alabama 35487-0336, USA
| | | | | | | | | | | |
Collapse
|
36
|
Feller D, Peterson KA, Dixon DA. Refined Theoretical Estimates of the Atomization Energies and Molecular Structures of Selected Small Oxygen Fluorides. J Phys Chem A 2009; 114:613-23. [DOI: 10.1021/jp908128g] [Citation(s) in RCA: 34] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- David Feller
- Department of Chemistry, Washington State University, Pullman, Washington 99164-4630
| | - Kirk A. Peterson
- Department of Chemistry, Washington State University, Pullman, Washington 99164-4630
| | - David A. Dixon
- Chemistry Department, The University of Alabama, Shelby Hall, Tuscaloosa, Alabama 35487-0336
| |
Collapse
|
37
|
Lukeš V, Šolc R, Lischka H, Kauffmann HF. Theoretical Study of the Relations between Structure and Photophysical Properties of Model Oligofluorenes with Central Keto Defect. J Phys Chem A 2009; 113:14141-9. [DOI: 10.1021/jp902658u] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Vladimír Lukeš
- Department of Chemical Physics, Slovak University of Technology, Radlinského 9, SK-812 37 Bratislava, Slovakia, Institute for Theoretical Chemistry, University of Vienna, Währingerstrasse 17, A-1090 Wien, Austria, and Institute of Physical Chemistry, University of Vienna, Währingerstrasse 42, A-1090 Vienna, Austria
| | - Roland Šolc
- Department of Chemical Physics, Slovak University of Technology, Radlinského 9, SK-812 37 Bratislava, Slovakia, Institute for Theoretical Chemistry, University of Vienna, Währingerstrasse 17, A-1090 Wien, Austria, and Institute of Physical Chemistry, University of Vienna, Währingerstrasse 42, A-1090 Vienna, Austria
| | - Hans Lischka
- Department of Chemical Physics, Slovak University of Technology, Radlinského 9, SK-812 37 Bratislava, Slovakia, Institute for Theoretical Chemistry, University of Vienna, Währingerstrasse 17, A-1090 Wien, Austria, and Institute of Physical Chemistry, University of Vienna, Währingerstrasse 42, A-1090 Vienna, Austria
| | - Harald Friedrich Kauffmann
- Department of Chemical Physics, Slovak University of Technology, Radlinského 9, SK-812 37 Bratislava, Slovakia, Institute for Theoretical Chemistry, University of Vienna, Währingerstrasse 17, A-1090 Wien, Austria, and Institute of Physical Chemistry, University of Vienna, Währingerstrasse 42, A-1090 Vienna, Austria
| |
Collapse
|
38
|
Feller D, Craig NC. High Level ab Initio Energies and Structures for the Rotamers of 1,3-Butadiene. J Phys Chem A 2009; 113:1601-7. [DOI: 10.1021/jp8095709] [Citation(s) in RCA: 37] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Affiliation(s)
- David Feller
- Department of Chemistry, Washington State University, Pullman, Washington 99164-4630, and Department of Chemistry and Biochemistry, Oberlin College, Oberlin, Ohio 44074
| | - Norman C. Craig
- Department of Chemistry, Washington State University, Pullman, Washington 99164-4630, and Department of Chemistry and Biochemistry, Oberlin College, Oberlin, Ohio 44074
| |
Collapse
|
39
|
A comparative study of modern and robust computational methods applied to π-complexes of moderate size: The case of the ethene/benzenium ion complex. Chem Phys Lett 2009. [DOI: 10.1016/j.cplett.2008.12.009] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
40
|
Feller D, Peterson KA, Dixon DA. A survey of factors contributing to accurate theoretical predictions of atomization energies and molecular structures. J Chem Phys 2008; 129:204105. [DOI: 10.1063/1.3008061] [Citation(s) in RCA: 312] [Impact Index Per Article: 19.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
41
|
Lukeš V, Šolc R, Barbatti M, Elstner M, Lischka H, Kauffmann HF. Torsional potentials and full-dimensional simulation of electronic absorption and fluorescence spectra of para-phenylene oligomers using the semiempirical self-consistent charge density-functional tight binding approach. J Chem Phys 2008; 129:164905. [DOI: 10.1063/1.2998523] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
42
|
Johansson MP, Olsen J. Torsional Barriers and Equilibrium Angle of Biphenyl: Reconciling Theory with Experiment. J Chem Theory Comput 2008; 4:1460-71. [DOI: 10.1021/ct800182e] [Citation(s) in RCA: 172] [Impact Index Per Article: 10.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Affiliation(s)
- Mikael P. Johansson
- Lundbeck Foundation Center for Theoretical Chemistry, Aarhus University, DK-8000 Århus C, Denmark, and Department of Chemistry, University of Helsinki, Helsinki, Finland
| | - Jeppe Olsen
- Lundbeck Foundation Center for Theoretical Chemistry, Aarhus University, DK-8000 Århus C, Denmark, and Department of Chemistry, University of Helsinki, Helsinki, Finland
| |
Collapse
|
43
|
Figgen D, Peterson KA, Stoll H. Energy-consistent relativistic pseudopotentials for the 4d elements: Atomic and molecular applications. J Chem Phys 2008; 128:034110. [DOI: 10.1063/1.2822992] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
44
|
Karton A, Taylor PR, Martin JML. Basis set convergence of post-CCSD contributions to molecular atomization energies. J Chem Phys 2007; 127:064104. [PMID: 17705585 DOI: 10.1063/1.2755751] [Citation(s) in RCA: 121] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Abstract
Basis set convergence of correlation effects on molecular atomization energies beyond the coupled cluster with singles and doubles (CCSD) approximation has been studied near the one-particle basis set limit. Quasiperturbative connected triple excitations, (T), converge more rapidly than L(-3) (where L is the highest angular momentum represented in the basis set), while higher-order connected triples, T3-(T), converge more slowly--empirically, proportional to L(-5/2). Quasiperturbative connected quadruple excitations, (Q), converge smoothly as proportional to L(-3) starting with the cc-pVTZ basis set, while the cc-pVDZ basis set causes overshooting of the contribution in highly polar systems. Higher-order connected quadruples display only weak, but somewhat erratic, basis set dependence. Connected quintuple excitations converge very rapidly with the basis set, to the point where even an unpolarized double-zeta basis set yields useful numbers. In cases where fully iterative coupled cluster up to connected quintuples (CCSDTQ5) calculations are not an option, CCSDTQ(5) (i.e., coupled cluster up to connected quadruples plus a quasiperturbative connected quintuples correction) cannot be relied upon in the presence of significant nondynamical correlation, whereas CCSDTQ(5)(Lambda) represents a viable alternative. Connected quadruples corrections to the core-valence contribution are thermochemically significant in some systems. We propose an additional variant of W4 theory [A. Karton et al., J. Chem. Phys. 125, 144108 (2006)], denoted W4.4 theory, which is shown to yield a rms deviation from experimental atomization energies (active thermochemical tables, ATcT) of only 0.05 kcal/mol for systems for which ATcT values are available. We conclude that "3sigma <or=1 kJ/mol" thermochemistry is feasible with current technology, but that the more ambitious goal of +/-10 cm(-1) accuracy is illusory, at least for atomization energies.
Collapse
Affiliation(s)
- Amir Karton
- Department of Organic Chemistry, Weizmann Institute of Science, IL-76100 Rehovot, Israel
| | | | | |
Collapse
|
45
|
Feller D, Peterson KA. Probing the limits of accuracy in electronic structure calculations: Is theory capable of results uniformly better than “chemical accuracy”? J Chem Phys 2007; 126:114105. [PMID: 17381194 DOI: 10.1063/1.2464112] [Citation(s) in RCA: 82] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023] Open
Abstract
Current limitations in electronic structure methods are discussed from the perspective of their potential to contribute to inherent uncertainties in predictions of molecular properties, with an emphasis on atomization energies (or heats of formation). The practical difficulties arising from attempts to achieve high accuracy are illustrated via two case studies: the carbon dimer (C2) and the hydroperoxyl radical (HO2). While the HO2 wave function is dominated by a single configuration, the carbon dimer involves considerable multiconfigurational character. In addition to these two molecules, statistical results will be presented for a much larger sample of molecules drawn from the Computational Results Database. The goal of this analysis will be to determine if a combination of coupled cluster theory with large 1-particle basis sets and careful incorporation of several computationally expensive smaller corrections can yield uniform agreement with experiment to better than "chemical accuracy" (+/-1 kcal/mol). In the case of HO2, the best current theoretical estimate of the zero-point-inclusive, spin-orbit corrected atomization energy (SigmaD0=166.0+/-0.3 kcal/mol) and the most recent Active Thermochemical Table (ATcT) value (165.97+/-0.06 kcal/mol) are in excellent agreement. For C2 the agreement is only slightly poorer, with theory (D0=143.7+/-0.3 kcal/mol) almost encompassing the most recent ATcT value (144.03+/-0.13 kcal/mol). For a larger collection of 68 molecules, a mean absolute deviation of 0.3 kcal/mol was found. The same high level of theory that produces good agreement for atomization energies also appears capable of predicting bond lengths to an accuracy of +/-0.001 A.
Collapse
Affiliation(s)
- David Feller
- Department of Chemistry, Washington State University, Pullman, Washington 99164-4630, USA.
| | | |
Collapse
|
46
|
Peterson KA, Shepler BC, Figgen D, Stoll H. On the spectroscopic and thermochemical properties of ClO, BrO, IO, and their anions. J Phys Chem A 2007; 110:13877-83. [PMID: 17181347 DOI: 10.1021/jp065887l] [Citation(s) in RCA: 610] [Impact Index Per Article: 35.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A coupled cluster composite approach has been used to accurately determine the spectroscopic constants, bond dissociation energies, and heats of formation for the X1(2)II(3/2) states of the halogen oxides ClO, BrO, and IO, as well as their negative ions ClO-, BrO-, and IO-. After determining the frozen core, complete basis set (CBS) limit CCSD(T) values, corrections were added for core-valence correlation, relativistic effects (scalar and spin-orbit), the pseudopotential approximation (BrO and IO), iterative connected triple excitations (CCSDT), and iterative quadruples (CCSDTQ). The final ab initio equilibrium bond lengths and harmonic frequencies for ClO and BrO differ from their accurate experimental values by an average of just 0.0005 A and 0.8 cm-1, respectively. The bond length of IO is overestimated by 0.0047 A, presumably due to an underestimation of molecular spin-orbit coupling effects. Spectroscopic constants for the spin-orbit excited X2(2)III(1/2) states are also reported for each species. The predicted bond lengths and harmonic frequencies for the closed-shell anions are expected to be accurate to within about 0.001 A and 2 cm-1, respectively. The dissociation energies of the radicals have been determined by both direct calculation and through use of negative ion thermochemical cycles, which made use of a small amount of accurate experimental data. The resulting values of D0, 63.5, 55.8, and 54.2 kcal/mol for ClO, BrO, and IO, respectively, are the most accurate ab initio values to date, and those for ClO and BrO differ from their experimental values by just 0.1 kcal/mol. These dissociation energies lead to heats of formation, DeltaH(f) (298 K), of 24.2 +/- 0.3, 29.6 +/- 0.4, and 29.9 +/- 0.6 kcal/mol for ClO, BrO, and IO, respectively. Also, the final calculated electron affinities are all within 0.2 kcal/mol of their experimental values. Improved pseudopotential parameters for the iodine atom are also reported, together with revised correlation consistent basis sets for this atom.
Collapse
Affiliation(s)
- Kirk A Peterson
- Department of Chemistry, Washington State University, Pullman, WA 99164-4630, USA.
| | | | | | | |
Collapse
|
47
|
Balabanov NB, Peterson KA. Basis set limit electronic excitation energies, ionization potentials, and electron affinities for the 3d transition metal atoms: Coupled cluster and multireference methods. J Chem Phys 2006; 125:074110. [PMID: 16942325 DOI: 10.1063/1.2335444] [Citation(s) in RCA: 207] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Recently developed correlation consistent basis sets for the first row transition metal elements Sc-Zn have been utilized to determine complete basis set (CBS) scalar relativistic electron affinities, ionization potentials, and 4s(2)3d(n-2)-4s(1)d(n-1) electronic excitation energies with single reference coupled cluster methods [CCSD(T), CCSDT, and CCSDTQ] and multireference configuration interaction with three reference spaces: 3d4s, 3d4s4p, and 3d4s4p3d'. The theoretical values calculated with the highest order coupled cluster techniques at the CBS limit, including extrapolations to full configuration interaction, are well within 1 kcal/mol of the corresponding experimental data. For the early transition metal elements (Sc-Mn) the internally contracted multireference averaged coupled pair functional method yielded excellent agreement with experiment; however, the atomic properties for the late transition metals (Mn-Zn) proved to be much more difficult to describe with this level of theory, even with the largest reference function of the present work.
Collapse
Affiliation(s)
- Nikolai B Balabanov
- Department of Chemistry, Washington State University, Pullman, Washington 99164-4630, USA.
| | | |
Collapse
|
48
|
Abstract
The convergence behavior of Møller-Plesset (MP) perturbation series is governed by the singularity structure of the energy, with the energy treated as a function of the perturbation parameter. Singularity locations, determined from quadratic approximant analysis of high-order series, are presented for a variety of atoms and small molecules. These results can be used as benchmarks for understanding the convergence of low-order methods such as MP4 and for developing and testing summation methods that model the singularity structure. The positions and types of singularities confirm previous qualitative predictions based on functional analysis of the Schrodinger equation.
Collapse
Affiliation(s)
- Alexey V Sergeev
- Department of Chemistry and Biochemistry, University of Massachusetts Dartmouth, North Dartmouth, Massachusetts 02747, USA
| | | |
Collapse
|
49
|
Feller D, Peterson KA, Crawford TD. Sources of error in electronic structure calculations on small chemical systems. J Chem Phys 2006; 124:054107. [PMID: 16468851 DOI: 10.1063/1.2137323] [Citation(s) in RCA: 153] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
Abstract
The sources of error in electronic structure calculations arising from the truncation of the one-particle and n-particle expansions are examined with very large correlation consistent basis sets, in some cases up through valence 10-zeta quality, and coupled cluster methods, up through connected quadruple excitations. A limited number of full configuration interaction corrections are also considered. For cases where full configuration interaction calculations were unavailable or prohibitively expensive, a continued fraction approximation was used. In addition, errors arising from corevalence and relativistic corrections are also probed for a number of small chemical systems. The accuracies of several formulas for estimating total energies and atomization energies in the complete basis set limit are compared in light of the present large basis set findings. In agreement with previous work, the CCSD(T) method is found to provide results that are closer to the CCSDTQ and full configuration-interaction results than the less approximate CCSDT method.
Collapse
Affiliation(s)
- David Feller
- Department of Chemistry, Washington State University, Pullman, Washington 99164-4630, USA.
| | | | | |
Collapse
|
50
|
Hou H, Deng L, Li J, Wang B. A Systematic Computational Study of the Reactions of HO2 with RO2: The HO2 + CH2ClO2, CHCl2O2, and CCl3O2 Reactions. J Phys Chem A 2005; 109:9299-309. [PMID: 16833272 DOI: 10.1021/jp052718c] [Citation(s) in RCA: 37] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Potential energy surfaces for the reactions of HO(2) with CH(2)ClO(2), CHCl(2)O(2), and CCl(3)O(2) have been calculated using coupled cluster theory and density functional theory (B3LYP). It is revealed that all the reactions take place on both singlet and triplet surfaces. Potential wells exist in the entrance channels for both surfaces. The reaction mechanism on the triplet surface is simple, including hydrogen abstraction and S(N)2-type displacement. The reaction mechanism on the singlet surface is more complicated. Interestingly, the corresponding transition states prefer to be 4-, 5-, or 7-member-ring structures. For the HO(2) + CH(2)ClO(2) reaction, there are two major product channels, viz., the formation of CH(2)ClOOH + O(2) via hydrogen abstraction on the triplet surface and the formation of CHClO + OH + HO(2) via a 5-member-ring transition state. Meanwhile, two O(3)-forming channels, namely, CH(2)O + HCl + O(3) and CH(2)ClOH + O(3) might be competitive at elevated temperatures. The HO(2) + CHCl(2)O(2) reaction has a mechanism similar to that of the HO(2) + CH(2)ClO(2) reaction. For the HO(2) + CCl(3)O(2) reaction, the formation of CCl(3)O(2)H + O(2) is the dominant channel. The Cl-substitution effect on the geometries, barriers, and heats of reaction is discussed. In addition, the unimolecular decomposition of the excited ROOH (e.g., CH(2)ClOOH, CHCl(2)OOH, and CCl(3)OOH) molecules has been investigated. The implication of the present mechanisms in atmospheric chemistry is discussed in comparison with the experimental measurements.
Collapse
Affiliation(s)
- Hua Hou
- College of Chemistry and Molecular Sciences, Wuhan University, Wuhan 430072, People's Republic of China
| | | | | | | |
Collapse
|