1
|
Heßelmann A, Giner E, Reinhardt P, Knowles PJ, Werner HJ, Toulouse J. A density-fitting implementation of the density-based basis-set correction method. J Comput Chem 2024; 45:1247-1253. [PMID: 38348951 DOI: 10.1002/jcc.27325] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2023] [Revised: 12/13/2023] [Accepted: 01/26/2024] [Indexed: 04/19/2024]
Abstract
This work reports an efficient density-fitting implementation of the density-based basis-set correction (DBBSC) method in the MOLPRO software. This method consists in correcting the energy calculated by a wave-function method with a given basis set by an adapted basis-set correction density functional incorporating the short-range electron correlation effects missing in the basis set, resulting in an accelerated convergence to the complete-basis-set limit. Different basis-set correction density-functional approximations are explored and the complementary-auxiliary-basis-set single-excitation correction is added. The method is tested on a benchmark set of reaction energies at the second-order Møller-Plesset (MP2) level and a comparison with the explicitly correlated MP2-F12 method is provided. The results show that the DBBSC method greatly accelerates the basis convergence of MP2 reaction energies, without reaching the accuracy of the MP2-F12 method but with a lower computational cost.
Collapse
Affiliation(s)
- Andreas Heßelmann
- Institute for Theoretical Chemistry, University of Stuttgart, Stuttgart, Germany
| | - Emmanuel Giner
- Laboratoire de Chimie Théorique, Sorbonne Université and CNRS, Paris, France
| | - Peter Reinhardt
- Laboratoire de Chimie Théorique, Sorbonne Université and CNRS, Paris, France
| | | | - Hans-Joachim Werner
- Institute for Theoretical Chemistry, University of Stuttgart, Stuttgart, Germany
| | - Julien Toulouse
- Laboratoire de Chimie Théorique, Sorbonne Université and CNRS, Paris, France
- Institut Universitaire de France, Paris, France
| |
Collapse
|
2
|
Werner HJ, Hansen A. Accurate Calculation of Isomerization and Conformational Energies of Larger Molecules Using Explicitly Correlated Local Coupled Cluster Methods in Molpro and ORCA. J Chem Theory Comput 2023; 19:7007-7030. [PMID: 37486154 DOI: 10.1021/acs.jctc.3c00270] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/25/2023]
Abstract
An overview of the approximations in the explicitly correlated local coupled cluster methods PNO-LCCSD(T)-F12 in Molpro and DLPNO-CCSD(T)F12 in ORCA is given. Options to select the domains of projected atomic orbitals (PAOs), pair natural orbitals (PNOs), and triples natural orbitals (TNOs) in both programs are described and compared in detail. The two programs are applied to compute isomerization and conformational energies of the ISOL24 and ACONFL test sets, where the former is part of the GMTKN55 benchmark suite. Thorough studies of basis set effects are presented for selected systems. These revealed large intramolecular basis set superposition effects that make it practically impossible to reliably determine the complete basis set (CBS) limits without including explicitly correlated terms. The latter strongly reduce the basis set dependence and at the same time also errors caused by the local domain approximations. On the basis of these studies, the PNO-LCCSD(T)-F12 method is applied to determine new reference energies for the above-mentioned benchmark sets. We are confident that our results should agree within a few tenths of a kcal mol-1 with the (unknown) CCSD(T)/CBS values, which therefore allowed us to define computational settings for accurate explicitly correlated local coupled cluster methods with moderate computational effort. With these protocols, especially PNO-LCCSD(T)-F12b/AVTZ', reliable reference values for comprehensive benchmark sets can be generated efficiently. This can significantly advance the development and evaluation of the performance of approximate electronic structure methods, especially improved density functional approximations or machine learning approaches.
Collapse
Affiliation(s)
- Hans-Joachim Werner
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstrasse 4, D-53115 Bonn, Germany
| |
Collapse
|
3
|
Zhao Q, Vaddadi SM, Woulfe M, Ogunfowora LA, Garimella SS, Isayev O, Savoie BM. Comprehensive exploration of graphically defined reaction spaces. Sci Data 2023; 10:145. [PMID: 36935430 PMCID: PMC10025260 DOI: 10.1038/s41597-023-02043-z] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Grants] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2022] [Accepted: 02/27/2023] [Indexed: 03/21/2023] Open
Abstract
Existing reaction transition state (TS) databases are comparatively small and lack chemical diversity. Here, this data gap has been addressed using the concept of a graphically-defined model reaction to comprehensively characterize a reaction space associated with C, H, O, and N containing molecules with up to 10 heavy (non-hydrogen) atoms. The resulting dataset is composed of 176,992 organic reactions possessing at least one validated TS, activation energy, heat of reaction, reactant and product geometries, frequencies, and atom-mapping. For 33,032 reactions, more than one TS was discovered by conformational sampling, allowing conformational errors in TS prediction to be assessed. Data is supplied at the GFN2-xTB and B3LYP-D3/TZVP levels of theory. A subset of reactions were recalculated at the CCSD(T)-F12/cc-pVDZ-F12 and ωB97X-D2/def2-TZVP levels to establish relative errors. The resulting collection of reactions and properties are called the Reaction Graph Depth 1 (RGD1) dataset. RGD1 represents the largest and most chemically diverse TS dataset published to date and should find immediate use in developing novel machine learning models for predicting reaction properties.
Collapse
Affiliation(s)
- Qiyuan Zhao
- Davidson School of Chemical Engineering, Purdue University, West Lafayette, IN, 47906, USA
| | - Sai Mahit Vaddadi
- Davidson School of Chemical Engineering, Purdue University, West Lafayette, IN, 47906, USA
| | - Michael Woulfe
- Davidson School of Chemical Engineering, Purdue University, West Lafayette, IN, 47906, USA
| | - Lawal A Ogunfowora
- Department of Chemistry, Purdue University, West Lafayette, IN, 47906, USA
| | - Sanjay S Garimella
- Davidson School of Chemical Engineering, Purdue University, West Lafayette, IN, 47906, USA
| | - Olexandr Isayev
- Department of Chemistry, Carnegie Mellon University, Pittsburgh, PA, 15213, USA
| | - Brett M Savoie
- Davidson School of Chemical Engineering, Purdue University, West Lafayette, IN, 47906, USA.
| |
Collapse
|
4
|
Zhao Q, Savoie BM. Algorithmic Explorations of Unimolecular and Bimolecular Reaction Spaces. Angew Chem Int Ed Engl 2022; 61:e202210693. [PMID: 36074520 PMCID: PMC9827825 DOI: 10.1002/anie.202210693] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2022] [Indexed: 01/12/2023]
Abstract
Algorithmic reaction exploration based on transition state searches has already made inroads into many niche applications, but its potential as a general-purpose tool is still largely unrealized. Computational cost and the absence of benchmark problems involving larger molecules remain obstacles to further progress. Here an ultra-low cost exploration algorithm is implemented and used to explore the reactivity of unimolecular and bimolecular reactants, comprising a total of 581 reactions involving 51 distinct reactants. The algorithm discovers all established reaction pathways, where such comparisons are possible, while also revealing a much richer reactivity landscape, including lower barrier reaction pathways and a strong dependence of reaction conformation in the apparent barriers of the reported reactions. The diversity of these benchmarks illustrate that reaction exploration algorithms are approaching general-purpose capability.
Collapse
Affiliation(s)
- Qiyuan Zhao
- Davidson School of Chemical EngineeringPurdue UniversityWest LafayetteIN47906USA
| | - Brett M. Savoie
- Davidson School of Chemical EngineeringPurdue UniversityWest LafayetteIN47906USA
| |
Collapse
|
5
|
Spiekermann K, Pattanaik L, Green WH. High accuracy barrier heights, enthalpies, and rate coefficients for chemical reactions. Sci Data 2022; 9:417. [PMID: 35851390 PMCID: PMC9293986 DOI: 10.1038/s41597-022-01529-6] [Citation(s) in RCA: 14] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2022] [Accepted: 06/30/2022] [Indexed: 12/13/2022] Open
Abstract
Quantitative chemical reaction data, including activation energies and reaction rates, are crucial for developing detailed kinetic mechanisms and accurately predicting reaction outcomes. However, such data are often difficult to find, and high-quality datasets are especially rare. Here, we use CCSD(T)-F12a/cc-pVDZ-F12//ωB97X-D3/def2-TZVP to obtain high-quality single point calculations for nearly 22,000 unique stable species and transition states. We report the results from these quantum chemistry calculations and extract the barrier heights and reaction enthalpies to create a kinetics dataset of nearly 12,000 gas-phase reactions. These reactions involve H, C, N, and O, contain up to seven heavy atoms, and have cleaned atom-mapped SMILES. Our higher-accuracy coupled-cluster barrier heights differ significantly (RMSE of ∼5 kcal mol-1) relative to those calculated at ωB97X-D3/def2-TZVP. We also report accurate transition state theory rate coefficients [Formula: see text] between 300 K and 2000 K and the corresponding Arrhenius parameters for a subset of rigid reactions. We believe this data will accelerate development of automated and reliable methods for quantitative reaction prediction.
Collapse
Affiliation(s)
- Kevin Spiekermann
- Department of Chemical Engineering, Massachusetts Institute of Technology, 77 Massachusetts Ave, Cambridge, MA, 02139, USA
| | - Lagnajit Pattanaik
- Department of Chemical Engineering, Massachusetts Institute of Technology, 77 Massachusetts Ave, Cambridge, MA, 02139, USA
| | - William H Green
- Department of Chemical Engineering, Massachusetts Institute of Technology, 77 Massachusetts Ave, Cambridge, MA, 02139, USA.
| |
Collapse
|
6
|
Villot C, Ballesteros F, Wang D, Lao KU. Coupled Cluster Benchmarking of Large Noncovalent Complexes in L7 and S12L as Well as the C 60 Dimer, DNA-Ellipticine, and HIV-Indinavir. J Phys Chem A 2022; 126:4326-4341. [PMID: 35766331 DOI: 10.1021/acs.jpca.2c01421] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
In this work, we report the benchmark binding energies of the seven complexes within the L7 data set, six host-guest complexes from the S12L data set, a C60 dimer, the DNA-ellipticine intercalation complex, and the largest system of the study, the HIV-indinavir system, which contained 343 atoms or 139 heavy atoms. The high-quality values reported were obtained via a focal point method that relies on the canonical form of second-order Møller-Plesset theory and the domain-based local pair natural orbital scheme for the coupled cluster with single double and perturbative triple excitations [DLPNO-CCSD(T)] extrapolated to the complete basis set (CBS) limit. The results in this work not only corroborate but also improve upon some previous benchmark values for large noncovalent complexes albeit at a relatively steep cost. Although local CCSD(T) and the largely successful fixed-node diffusion Monte Carlo (FN-DMC) have been shown to generally agree for small- to medium-size systems, a discrepancy in their reported binding energy values arises for large complexes, where the magnitude of the disagreement is a definite cause for concern. For example, the largest deviation in the L7 data set was 2.8 kcal/mol (∼10%) on the low end in C3GC. Such a deviation only grows worse in the S12L set, which showed a difference of up to 10.4 kcal/mol (∼25%) by a conservative estimation in buckycatcher-C60. The DNA-ellipticine complex also generated a disagreement of 4.4 kcal/mol (∼10%) between both state-of-the-art methods. The disagreement between local CCSD(T) and FN-DMC in large noncovalent complexes shows that it is urgently needed to have the canonical CCSD(T), the Monte Carlo CCSD(T), or the full configuration interaction quantum Monte Carlo approaches available to large systems on the hundred-atom scale to solve this dilemma. In addition, the performances of cheaper popular computational methods were assessed for the studied complexes with respect to DLPNO-CCSD(T)/CBS. r2SCAN-3c, B97M-V, and PBE0+D4 work well in large noncovalent complexes in this work, and GFN2-xTB performs well in π-π stacking complexes. B97M-V is the most reliable computationally efficient approach to predicting noncovalent interactions for large complexes, being the only one to have binding errors within the so-called 1 kcal/mol "chemical accuracy". The benchmark interaction energies of these host-guest complexes, molecular materials, and biological systems with electronic and medicinal implications provide crucial reference data for the improvement of current and future lower-cost methods.
Collapse
Affiliation(s)
- Corentin Villot
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284 United States
| | - Francisco Ballesteros
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284 United States
| | - Danyang Wang
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284 United States
| | - Ka Un Lao
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284 United States
| |
Collapse
|
7
|
Spiekermann KA, Pattanaik L, Green WH. Fast Predictions of Reaction Barrier Heights: Toward Coupled-Cluster Accuracy. J Phys Chem A 2022; 126:3976-3986. [PMID: 35727075 DOI: 10.1021/acs.jpca.2c02614] [Citation(s) in RCA: 23] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Quantitative estimates of reaction barriers are essential for developing kinetic mechanisms and predicting reaction outcomes. However, the lack of experimental data and the steep scaling of accurate quantum calculations often hinder the ability to obtain reliable kinetic values. Here, we train a directed message passing neural network on nearly 24,000 diverse gas-phase reactions calculated at CCSD(T)-F12a/cc-pVDZ-F12//ωB97X-D3/def2-TZVP. Our model uses 75% fewer parameters than previous studies, an improved reaction representation, and proper data splits to accurately estimate performance on unseen reactions. Using information from only the reactant and product, our model quickly predicts barrier heights with a testing MAE of 2.6 kcal mol-1 relative to the coupled-cluster data, making it more accurate than a good density functional theory calculation. Furthermore, our results show that future modeling efforts to estimate reaction properties would significantly benefit from fine-tuning calibration using a transfer learning technique. We anticipate this model will accelerate and improve kinetic predictions for small molecule chemistry.
Collapse
Affiliation(s)
- Kevin A Spiekermann
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Lagnajit Pattanaik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - William H Green
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
8
|
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
|
9
|
Varandas AJC. Canonical and explicitly-correlated coupled cluster correlation energies of sub-kJ mol -1 accuracy via cost-effective hybrid-post-CBS extrapolation. Phys Chem Chem Phys 2021; 23:9571-9584. [PMID: 33885095 DOI: 10.1039/d1cp00357g] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Cost-effectiveness and accuracy are two basic pillars in electronic structure calculations. While cost-effectiveness enhances applicability, high accuracy is sustained when employing advanced computational tools. With the gold standard method of ab initio quantum chemistry at the focal point, canonical CCSD(T) and modern explicitly correlated CCSD(T)-F12 calculations are employed hand in hand to develop accurate hybrid post-CBS extrapolation schemes, which are validated using popular training sets involving a total of 130 molecules. By using raw valence-only calculations at CCSD(T)/VDZ and CCSD(T)/VQZ-F12 levels of theory, the novel scheme leads to the prediction of absolute energies that differ on average (-0.170 ± 0.224) kcal mol-1 from the highest affordable CCSD(T)-F12b/V(Q,5)Z-F12 extrapolations, but only (-0.048 ± 0.228) kcal mol-1 from the post-CBS extrapolated values based on CBS(D,T), CBS(D,Q) and CBS(T,Q) energies. From the cost-effectiveness standpoint, the approach is a kind of pseudo one-point extrapolation scheme since its cost is basically that of the highest-rung raw energy where it is based. Variants that imply no additional cost are also discussed, emerging h-pCBS(dt,dq)ab as the most effective. The approach can also be used with PNO-based local correlation methods that gained popularity due to allowing coupled-cluster calculations even for large molecules at reduced computational cost, namely local PNO-CCSD(T) and PNO-CCSD(T)-F12b. To gauge the approach performance, both the hydrogen molecule and the O-C2H5 torsion path of ethyl-methyl-ether, an extra molecule here considered with presupposed existence in astrophysical objects, are also studied. Additionally, the nonbonding interactions in the A24 test set are revisited per se. The results show that the title approach may be useful in high-accuracy quantum chemistry, with further improvements requiring the inclusion of contributions beyond the theory here employed such as the ones due to relativistic and nonadiabatic effects.
Collapse
Affiliation(s)
- A J C Varandas
- School of Physics and Physical Engineering, Qufu Normal University, 273165 Qufu, China
| |
Collapse
|
10
|
Sun Y, Somers KP, Wang QD, Farrell C, Curran HJ. Hindered rotor benchmarks for the transition states of free radical additions to unsaturated hydrocarbons. Phys Chem Chem Phys 2020; 22:27241-27254. [PMID: 33226373 DOI: 10.1039/d0cp04194g] [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/21/2022]
Abstract
The hindered internal rotors of 32 transition states (TSs) formed through four free radicals, namely methyl, vinyl, ethyl, methoxy (ĊH3, Ċ2H3, Ċ2H5, CH3) additions to acetylene, ethylene, allene, propyne, and propene (C2H2/C2H4/C3H4-a/C3H4-p/C3H6) are studied. To validate the uncertainties of rate constants that stem from the use of different electronic structure methods to treat hindered rotors, the rotations of the newly formed C-C and/or C-O rotors in the transition states are calculated using commonly used DFT methods (B3LYP, M06-2X, ωB97X-D and B2PLYP-D3 with two Pople basis sets (6-31+G(d,p), 6-311++G(d,p)) and cc-pVTZ). The hindrance potential energies V(χ) calculated using the M06-2X/6-311++G(d,p) method are benchmarked at the CCSD(T), CCSD(T)-F12, DLPNO-CCSD(T) levels of theory with cc-pVTZ-F12 and cc-pVXZ (X = T, Q) basis sets and are extrapolated to the complete basis set (CBS) limit. The DLPNO-CCSD(T)/CBS method is proven to reproduce the CCSD(T)/CBS energies within 0.5 kJ mol-1 and this method is selected as the benchmark for all of the rotors in this study. Rotational constants B(χ) are computed for each method based on the optimized geometries for the hindrance potential via the I(2,3) approximation. Thereafter, the V(χ) and B(χ) values are used to compute hindered internal rotation partition functions, QHR, as a function of temperature. The uncertainties in the V(χ), B(χ) and QHR calculations stem from the use of different DFT methods for the internal rotor treatment are discussed for these newly formed rotors. For rotors formed by + C2 alkenes/alkynes, the V(χ) and QHR values calculated using DFT methods are compared with the DLPNO-CCSD(T)/CBS results and analysed according to reaction types. Based on comparisons of the DFT methods with the benchmarking method, reliable DFT methods are recommended for the treatment of internal rotors for different reaction types considering both accuracy and computational cost. This work, to the authors' knowledge, is the first to systematically benchmark hindrance potentials which can be used to estimate uncertainties in theoretically derived rate constants arising from the choice of different electronic structure methods.
Collapse
Affiliation(s)
- Yanjin Sun
- Combustion Chemistry Centre, School of Chemistry, Ryan Institute, MaREI, National University of Ireland, Galway, Ireland.
| | | | | | | | | |
Collapse
|
11
|
Sorathia K, Tew DP. Basis set extrapolation in pair natural orbital theories. J Chem Phys 2020; 153:174112. [PMID: 33167642 DOI: 10.1063/5.0022077] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present the results of a benchmark study of the effect of Pair Natural Orbital (PNO) truncation errors on the performance of basis set extrapolation. We find that reliable conclusions from the application of Helgaker's extrapolation method are only obtained when using tight PNO thresholds of at least 10-7. The use of looser thresholds introduces a significant risk of observing a false basis set convergence and underestimating the residual basis set errors. We propose an alternative extrapolation approach based on the PNO truncation level that only requires a single basis set and show that it is a viable alternative to hierarchical basis set extrapolation methods.
Collapse
Affiliation(s)
- Kesha Sorathia
- Max Planck Institute for Solid State Research, Heisenbergstrasse 1, 70569 Stuttgart, GermanyUniversity of Oxford, South Parks Road, Oxford OX1 3QZ, United Kingdom
| | - David P Tew
- Max Planck Institute for Solid State Research, Heisenbergstrasse 1, 70569 Stuttgart, Germany
| |
Collapse
|
12
|
Oliveira VP, Kraka E, Machado FBC. Pushing 3c–4e Bonds to the Limit: A Coupled Cluster Study of Stepwise Fluorination of First-Row Atoms. Inorg Chem 2019; 58:14777-14789. [DOI: 10.1021/acs.inorgchem.9b02458] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Vytor P. Oliveira
- Instituto Tecnológico de Aeronáutica (ITA), Departamento de Química, São José dos Campos, 12228-900 São Paulo, Brazil
| | - Elfi Kraka
- Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| | - Francisco B. C. Machado
- Instituto Tecnológico de Aeronáutica (ITA), Departamento de Química, São José dos Campos, 12228-900 São Paulo, Brazil
| |
Collapse
|
13
|
Shaw RA, Hill JG. A Linear-Scaling Method for Noncovalent Interactions: An Efficient Combination of Absolutely Localized Molecular Orbitals and a Local Random Phase Approximation Approach. J Chem Theory Comput 2019; 15:5352-5369. [PMID: 31465215 DOI: 10.1021/acs.jctc.9b00615] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A novel method for the accurate and efficient calculation of interaction energies in weakly bound complexes composed of a large number of molecules is presented. The new ALMO+RPAd method circumvents the prohibitive scaling of coupled cluster singles and doubles while still providing similar accuracy across a diverse range of weakly bound chemical systems. Linear-scaling procedures for the Fock build are given utilizing absolutely localized molecular orbitals (ALMOs), resulting in the a priori exclusion of basis set superposition errors. A bespoke data structure and algorithm using density fitting are described, leading to linear scaling for the storage and computation of the two-electron integrals. Electron correlation is included through a new, linear-scaling pairwise local random phase approximation approach, including exchange interactions, and decomposed into purely dispersive excitations (RPAxd). Collectively, these allow meaningful decomposition of the interaction energy into physically distinct contributions: electrostatic, polarization, charge transfer, and dispersion. Comparison with symmetry-adapted perturbation theory shows good qualitative agreement. Tests on various dimers and the S66 benchmark set demonstrate results within 0.5 kcal mol-1 of coupled cluster singles and doubles results. On a large cluster of water molecules, we achieve calculations involving over 3500 orbital and 12,000 auxiliary basis functions in under 10 min on a single CPU core.
Collapse
Affiliation(s)
- Robert A Shaw
- Department of Chemistry , University of Sheffield , Sheffield S3 7HF , United Kingdom
| | - J Grant Hill
- Department of Chemistry , University of Sheffield , Sheffield S3 7HF , United Kingdom
| |
Collapse
|
14
|
Riera M, Lambros E, Nguyen TT, Götz AW, Paesani F. Low-order many-body interactions determine the local structure of liquid water. Chem Sci 2019; 10:8211-8218. [PMID: 32133122 PMCID: PMC6927411 DOI: 10.1039/c9sc03291f] [Citation(s) in RCA: 29] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/04/2019] [Accepted: 07/21/2019] [Indexed: 12/30/2022] Open
Abstract
Two-body and three-body energies, modulated by higher-body terms and nuclear quantum effects, determine the structure of liquid water and require sub-chemical accuracy that is achieved by the MB-pol model but not by existing DFT functionals.
Despite its apparent simplicity, water displays unique behavior across the phase diagram which is strictly related to the ability of the water molecules to form dense, yet dynamic, hydrogen-bond networks that continually fluctuate in time and space. The competition between different local hydrogen-bonding environments has been hypothesized as a possible origin of the anomalous properties of liquid water. Through a systematic application of the many-body expansion of the total energy, we demonstrate that the local structure of liquid water at room temperature is determined by a delicate balance between two-body and three-body energies, which is further modulated by higher-order many-body effects. Besides providing fundamental insights into the structure of liquid water, this analysis also emphasizes that a correct representation of two-body and three-body energies requires sub-chemical accuracy that is nowadays only achieved by many-body models rigorously derived from the many-body expansion of the total energy, which thus hold great promise for shedding light on the molecular origin of the anomalous behavior of liquid water.
Collapse
Affiliation(s)
- Marc Riera
- Department of Chemistry and Biochemistry , University of California , San Diego , La Jolla , California 92093 , USA
| | - Eleftherios Lambros
- Department of Chemistry and Biochemistry , University of California , San Diego , La Jolla , California 92093 , USA
| | - Thuong T Nguyen
- Department of Chemistry and Biochemistry , University of California , San Diego , La Jolla , California 92093 , USA
| | - Andreas W Götz
- San Diego Supercomputer Center , University of California , San Diego , La Jolla , California 92093 , USA
| | - Francesco Paesani
- Department of Chemistry and Biochemistry , University of California , San Diego , La Jolla , California 92093 , USA.,San Diego Supercomputer Center , University of California , San Diego , La Jolla , California 92093 , USA.,Materials Science and Engineering , University of California , San Diego , La Jolla , California 92093 , USA .
| |
Collapse
|
15
|
Höfener S, Schieschke N, Klopper W, Köhn A. The extended explicitly-correlated second-order approximate coupled-cluster singles and doubles ansatz suitable for response theory. J Chem Phys 2019; 150:184110. [DOI: 10.1063/1.5094434] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Affiliation(s)
- Sebastian Höfener
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), P.O. Box 6980, D-76049 Karlsruhe, Germany
| | - Nils Schieschke
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), P.O. Box 6980, D-76049 Karlsruhe, Germany
| | - Wim Klopper
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), P.O. Box 6980, D-76049 Karlsruhe, Germany
| | - Andreas Köhn
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| |
Collapse
|
16
|
Mattalia JM, Nava P. C-C Bond Breaking in Addition-Elimination Reactions on Nitriles. European J Org Chem 2019. [DOI: 10.1002/ejoc.201900119] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
| | - Paola Nava
- CNRS, Centrale Marseille, iSm2, Marseille; Aix-Marseille Univ; France
| |
Collapse
|
17
|
Ma Q, Werner H. Explicitly correlated local coupled‐cluster methods using pair natural orbitals. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2018. [DOI: 10.1002/wcms.1371] [Citation(s) in RCA: 100] [Impact Index Per Article: 16.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Affiliation(s)
- Qianli Ma
- Institute for Theoretical ChemistryUniversity of StuttgartStuttgartGermany
| | | |
Collapse
|
18
|
Győrffy W, Werner HJ. Analytical energy gradients for explicitly correlated wave functions. II. Explicitly correlated coupled cluster singles and doubles with perturbative triples corrections: CCSD(T)-F12. J Chem Phys 2018; 148:114104. [DOI: 10.1063/1.5020436] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Affiliation(s)
- Werner Győrffy
- 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
|
19
|
Varandas AJC. Straightening the Hierarchical Staircase for Basis Set Extrapolations: A Low-Cost Approach to High-Accuracy Computational Chemistry. Annu Rev Phys Chem 2018; 69:177-203. [PMID: 29394151 DOI: 10.1146/annurev-physchem-050317-021148] [Citation(s) in RCA: 52] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Because the one-electron basis set limit is difficult to reach in correlated post-Hartree-Fock ab initio calculations, the low-cost route of using methods that extrapolate to the estimated basis set limit attracts immediate interest. The situation is somewhat more satisfactory at the Hartree-Fock level because numerical calculation of the energy is often affordable at nearly converged basis set levels. Still, extrapolation schemes for the Hartree-Fock energy are addressed here, although the focus is on the more slowly convergent and computationally demanding correlation energy. Because they are frequently based on the gold-standard coupled-cluster theory with single, double, and perturbative triple excitations [CCSD(T)], correlated calculations are often affordable only with the smallest basis sets, and hence single-level extrapolations from one raw energy could attain maximum usefulness. This possibility is examined. Whenever possible, this review uses raw data from second-order Møller-Plesset perturbation theory, as well as CCSD, CCSD(T), and multireference configuration interaction methods. Inescapably, the emphasis is on work done by the author's research group. Certain issues in need of further research or review are pinpointed.
Collapse
Affiliation(s)
- António J C Varandas
- Coimbra Chemistry Center and Department of Chemistry, University of Coimbra, Coimbra 3004-535, Portugal;
| |
Collapse
|
20
|
|
21
|
Ma Q, Werner HJ. Scalable Electron Correlation Methods. 5. Parallel Perturbative Triples Correction for Explicitly Correlated Local Coupled Cluster with Pair Natural Orbitals. J Chem Theory Comput 2017; 14:198-215. [DOI: 10.1021/acs.jctc.7b01141] [Citation(s) in RCA: 59] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
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
|
22
|
Győrffy W, Knizia G, Werner HJ. Analytical energy gradients for explicitly correlated wave functions. I. Explicitly correlated second-order Møller-Plesset perturbation theory. J Chem Phys 2017; 147:214101. [DOI: 10.1063/1.5003065] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023] Open
Affiliation(s)
- Werner Győrffy
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - Gerald Knizia
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
- Department of Chemistry, Pennsylvania State University, 401A Chemistry Building, University Park, Pennsylvania 16802, USA
| | - Hans-Joachim Werner
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| |
Collapse
|
23
|
Ma Q, Schwilk M, Köppl C, Werner HJ. Scalable Electron Correlation Methods. 4. Parallel Explicitly Correlated Local Coupled Cluster with Pair Natural Orbitals (PNO-LCCSD-F12). J Chem Theory Comput 2017; 13:4871-4896. [DOI: 10.1021/acs.jctc.7b00799] [Citation(s) in RCA: 65] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Qianli Ma
- Institut für Theoretische
Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - Max Schwilk
- Institut für Theoretische
Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - Christoph Köppl
- 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
|
24
|
Schwilk M, Ma Q, Köppl C, Werner HJ. Scalable Electron Correlation Methods. 3. Efficient and Accurate Parallel Local Coupled Cluster with Pair Natural Orbitals (PNO-LCCSD). J Chem Theory Comput 2017; 13:3650-3675. [DOI: 10.1021/acs.jctc.7b00554] [Citation(s) in RCA: 100] [Impact Index Per Article: 14.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Affiliation(s)
- Max Schwilk
- Institut für Theoretische
Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - Qianli Ma
- Institut für Theoretische
Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - Christoph Köppl
- 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
|
25
|
Shaw RA, Hill JG. Approaching the Hartree–Fock Limit through the Complementary Auxiliary Basis Set Singles Correction and Auxiliary Basis Sets. J Chem Theory Comput 2017; 13:1691-1698. [DOI: 10.1021/acs.jctc.7b00140] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Robert A. Shaw
- Department of Chemistry, University of Sheffield, Sheffield S3 7HF, U.K
| | - J. Grant Hill
- Department of Chemistry, University of Sheffield, Sheffield S3 7HF, U.K
| |
Collapse
|
26
|
Reddy SK, Straight SC, Bajaj P, Huy Pham C, Riera M, Moberg DR, Morales MA, Knight C, Götz AW, Paesani F. On the accuracy of the MB-pol many-body potential for water: Interaction energies, vibrational frequencies, and classical thermodynamic and dynamical properties from clusters to liquid water and ice. J Chem Phys 2016; 145:194504. [DOI: 10.1063/1.4967719] [Citation(s) in RCA: 161] [Impact Index Per Article: 20.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Affiliation(s)
- Sandeep K. Reddy
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| | - Shelby C. Straight
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| | - Pushp Bajaj
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| | - C. Huy Pham
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| | - Marc Riera
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| | - Daniel R. Moberg
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| | - Miguel A. Morales
- Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California 94550, USA
| | - Chris Knight
- Leadership Computing Facility, Argonne National Laboratory, 9700 South Cass Avenue, Argonne, Illinois 60439, USA
| | - Andreas W. Götz
- San Diego Supercomputer Center, University of California, San Diego, La Jolla, California 92093, USA
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| |
Collapse
|
27
|
Tew DP. Explicitly correlated coupled-cluster theory with Brueckner orbitals. J Chem Phys 2016; 145:074103. [DOI: 10.1063/1.4960655] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- David P. Tew
- School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| |
Collapse
|
28
|
On the performance of various hierarchized bases in extrapolating the correlation energy to the complete basis set limit. Chem Phys Lett 2015. [DOI: 10.1016/j.cplett.2015.10.064] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
29
|
Medders GR, Götz AW, Morales MA, Bajaj P, Paesani F. On the representation of many-body interactions in water. J Chem Phys 2015; 143:104102. [DOI: 10.1063/1.4930194] [Citation(s) in RCA: 78] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Gregory R. Medders
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| | - Andreas W. Götz
- San Diego Supercomputer Center, University of California, San Diego, La Jolla, California 92093, USA
| | - Miguel A. Morales
- Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California 94550, USA
| | - Pushp Bajaj
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| |
Collapse
|
30
|
Pansini F, Varandas A. Toward a unified single-parameter extrapolation scheme for the correlation energy: Systems formed by atoms of hydrogen through neon. Chem Phys Lett 2015. [DOI: 10.1016/j.cplett.2015.04.052] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
31
|
Hehn AS, Tew DP, Klopper W. Explicitly correlated ring-coupled-cluster-doubles theory. J Chem Phys 2015; 142:194106. [DOI: 10.1063/1.4921256] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
32
|
Lao KU, Schäffer R, Jansen G, Herbert JM. Accurate description of intermolecular interactions involving ions using symmetry-adapted perturbation theory. J Chem Theory Comput 2015; 11:2473-86. [PMID: 26575547 DOI: 10.1021/ct5010593] [Citation(s) in RCA: 80] [Impact Index Per Article: 8.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Three new data sets for intermolecular interactions, AHB21 for anion-neutral dimers, CHB6 for cation-neutral dimers, and IL16 for ion pairs, are assembled here, with complete-basis CCSD(T) results for each. These benchmarks are then used to evaluate the accuracy of the single-exchange approximation that is used for exchange energies in symmetry-adapted perturbation theory (SAPT), and the accuracy of SAPT based on wave function and density-functional descriptions of the monomers is evaluated. High-level SAPT calculations afford poor results for these data sets, and this includes the recently proposed "gold", "silver", and "bronze standards" of SAPT, namely, SAPT2+(3)-δMP2/aug-cc-pVTZ, SAPT2+/aug-cc-pVDZ, and sSAPT0/jun-cc-pVDZ, respectively [ Parker , T. M. , et al. , J. Chem. Phys. 2014 , 140 , 094106 ]. Especially poor results are obtained for symmetric shared-proton systems of the form X(-)···H(+)···X(-), for X = F, Cl, or OH. For the anionic data set, the SAPT2+(CCD)-δMP2/aug-cc-pVTZ method exhibits the best performance, with a mean absolute error (MAE) of 0.3 kcal/mol and a maximum error of 0.7 kcal/mol. For the cationic data set, the highest-level SAPT method, SAPT2+3-δMP2/aug-cc-pVQZ, outperforms the rest of the SAPT methods, with a MAE of 0.2 kcal/mol and a maximum error of 0.4 kcal/mol. For the ion-pair data set, the SAPT2+3-δMP2/aug-cc-pVTZ performs the best among all SAPT methods with a MAE of 0.3 kcal/mol and a maximum error of 0.9 kcal/mol. Overall, SAPT2+3-δMP2/aug-cc-pVTZ affords a small and balanced MAE (<0.5 kcal/mol) for all three data sets, with an overall MAE of 0.4 kcal/mol. Despite the breakdown of perturbation theory for ionic systems at short-range, SAPT can still be saved given two corrections: a "δHF" correction, which requires a supermolecular Hartree-Fock calculation to incorporate polarization effects beyond second order, and a "δMP2" correction, which requires a supermolecular MP2 calculation to account for higher-order induction-dispersion coupling. These corrections serve to remove artifacts introduced by the single exchange approximation in the exchange-induction and exchange-dispersion interactions, and obviate the need for ad hoc scaling of the first- and second-order exchange energies. Finally, some density-functional and MP2-based electronic structure methods are assessed as well, and we find that the best density-functional method for computing binding energies in these data sets is B97M-V/aug-cc-pVTZ, which affords a MAE of 0.4 kcal/mol, whereas complete-basis MP2 affords an MAE of 0.3 kcal/mol.
Collapse
Affiliation(s)
- Ka Un Lao
- Department of Chemistry and Biochemistry, The Ohio State University , Columbus, Ohio 43210, United States
| | - Rainer Schäffer
- Fakultät für Chemie, Universität Duisburg-Essen , 45117 Essen, Germany
| | - Georg Jansen
- Fakultät für Chemie, Universität Duisburg-Essen , 45117 Essen, Germany
| | - John M Herbert
- Department of Chemistry and Biochemistry, The Ohio State University , Columbus, Ohio 43210, United States
| |
Collapse
|
33
|
Shallcross DE, Leather KE, Bacak A, Xiao P, Lee EPF, Ng M, Mok DKW, Dyke JM, Hossaini R, Chipperfield MP, Khan MAH, Percival CJ. Reaction between CH3O2 and BrO Radicals: A New Source of Upper Troposphere Lower Stratosphere Hydroxyl Radicals. J Phys Chem A 2015; 119:4618-32. [DOI: 10.1021/jp5108203] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
| | - Kimberley E. Leather
- School
of Earth, Atmospheric and Environmental Sciences, The University of Manchester, Williamson Building, Oxford Road, Manchester M13 9PL, U.K
| | - Asan Bacak
- School
of Earth, Atmospheric and Environmental Sciences, The University of Manchester, Williamson Building, Oxford Road, Manchester M13 9PL, U.K
| | - Ping Xiao
- School
of Chemistry, University of Bristol, Bristol BS8 1TS, U.K
| | - Edmond P. F. Lee
- School
of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, U.K
- Department
of Applied Biology and Chemical Technology, The Hong Kong Polytechnic University, Hung Hom, Kowloon, Hong Kong
| | - Maggie Ng
- Department
of Applied Biology and Chemical Technology, The Hong Kong Polytechnic University, Hung Hom, Kowloon, Hong Kong
| | - Daniel K. W. Mok
- Department
of Applied Biology and Chemical Technology, The Hong Kong Polytechnic University, Hung Hom, Kowloon, Hong Kong
| | - John M. Dyke
- School
of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, U.K
| | - Ryan Hossaini
- School
of Earth and Environment, University of Leeds, Leeds LS2 9JT, U.K
| | | | - M. Anwar H. Khan
- School
of Chemistry, University of Bristol, Bristol BS8 1TS, U.K
| | - Carl J. Percival
- School
of Earth, Atmospheric and Environmental Sciences, The University of Manchester, Williamson Building, Oxford Road, Manchester M13 9PL, U.K
| |
Collapse
|
34
|
Temelso B, Renner CR, Shields GC. Importance and Reliability of Small Basis Set CCSD(T) Corrections to MP2 Binding and Relative Energies of Water Clusters. J Chem Theory Comput 2015; 11:1439-48. [DOI: 10.1021/ct500944v] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Berhane Temelso
- Dean’s
Office, College of Arts and Sciences, and Department of Chemistry, Bucknell University, Lewisburg, Pennsylvania 17837, United States
| | - Carla R. Renner
- Dean’s
Office, College of Arts and Sciences, and Department of Chemistry, Bucknell University, Lewisburg, Pennsylvania 17837, United States
| | - George C. Shields
- Dean’s
Office, College of Arts and Sciences, and Department of Chemistry, Bucknell University, Lewisburg, Pennsylvania 17837, United States
| |
Collapse
|
35
|
Class CA, Aguilera-Iparraguirre J, Green WH. A kinetic and thermochemical database for organic sulfur and oxygen compounds. Phys Chem Chem Phys 2015; 17:13625-39. [DOI: 10.1039/c4cp05631k] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Potential energy surfaces and reaction kinetics were calculated for reactions involving sulfur and oxygen, which are potentially relevant in combustion and desulfurization chemistry.
Collapse
Affiliation(s)
- Caleb A. Class
- Department of Chemical Engineering
- Massachusetts Institute of Technology
- Cambridge
- USA
| | | | - William H. Green
- Department of Chemical Engineering
- Massachusetts Institute of Technology
- Cambridge
- USA
| |
Collapse
|
36
|
Lao KU, Herbert JM. Accurate and Efficient Quantum Chemistry Calculations for Noncovalent Interactions in Many-Body Systems: The XSAPT Family of Methods. J Phys Chem A 2014; 119:235-52. [DOI: 10.1021/jp5098603] [Citation(s) in RCA: 72] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Ka Un Lao
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| | - John M. Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| |
Collapse
|
37
|
Grüneis A, Shepherd JJ, Alavi A, Tew DP, Booth GH. Explicitly correlated plane waves: Accelerating convergence in periodic wavefunction expansions. J Chem Phys 2013; 139:084112. [DOI: 10.1063/1.4818753] [Citation(s) in RCA: 56] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
38
|
Hollman DS, Wilke JJ, Schaefer HF. Explicitly correlated atomic orbital basis second order Møller–Plesset theory. J Chem Phys 2013; 138:064107. [DOI: 10.1063/1.4790582] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
39
|
Zhang J, Valeev EF. Prediction of Reaction Barriers and Thermochemical Properties with Explicitly Correlated Coupled-Cluster Methods: A Basis Set Assessment. J Chem Theory Comput 2012; 8:3175-86. [DOI: 10.1021/ct3005547] [Citation(s) in RCA: 62] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Jinmei Zhang
- Department of Chemistry, Virginia Tech, Blacksburg,
Virginia 24061, United States
| | - Edward F. Valeev
- Department of Chemistry, Virginia Tech, Blacksburg,
Virginia 24061, United States
| |
Collapse
|
40
|
Husar DE, Temelso B, Ashworth AL, Shields GC. Hydration of the Bisulfate Ion: Atmospheric Implications. J Phys Chem A 2012; 116:5151-63. [DOI: 10.1021/jp300717j] [Citation(s) in RCA: 56] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Devon E. Husar
- Dean’s Office, College of
Arts and Sciences,
and Department of Chemistry, Bucknell University, Lewisburg, Pennsylvania 17837, United States
| | - Berhane Temelso
- Dean’s Office, College of
Arts and Sciences,
and Department of Chemistry, Bucknell University, Lewisburg, Pennsylvania 17837, United States
| | - Alexa L. Ashworth
- Dean’s Office, College of
Arts and Sciences,
and Department of Chemistry, Bucknell University, Lewisburg, Pennsylvania 17837, United States
| | - George C. Shields
- Dean’s Office, College of
Arts and Sciences,
and Department of Chemistry, Bucknell University, Lewisburg, Pennsylvania 17837, United States
| |
Collapse
|
41
|
Perez C, Muckle MT, Zaleski DP, Seifert NA, Temelso B, Shields GC, Kisiel Z, Pate BH. Structures of Cage, Prism, and Book Isomers of Water Hexamer from Broadband Rotational Spectroscopy. Science 2012; 336:897-901. [DOI: 10.1126/science.1220574] [Citation(s) in RCA: 335] [Impact Index Per Article: 27.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/02/2022]
|
42
|
Ryding MJ, Ruusuvuori K, Andersson PU, Zatula AS, McGrath MJ, Kurtén T, Ortega IK, Vehkamäki H, Uggerud E. Structural Rearrangements and Magic Numbers in Reactions between Pyridine-Containing Water Clusters and Ammonia. J Phys Chem A 2012; 116:4902-8. [DOI: 10.1021/jp3021326] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Affiliation(s)
- Mauritz J. Ryding
- Mass Spectrometry Laboratory
and Centre of Theoretical and Computational Chemistry, Department
of Chemistry, University of Oslo, P.O.
Box 1033 Blindern, NO-0315 Oslo, Norway
| | - Kai Ruusuvuori
- Division of Atmospheric Sciences,
Department of Physics, University of Helsinki, P.O. Box 64, FI-00014 University of Helsinki, Helsinki, Finland
| | | | - Alexey S. Zatula
- Mass Spectrometry Laboratory
and Centre of Theoretical and Computational Chemistry, Department
of Chemistry, University of Oslo, P.O.
Box 1033 Blindern, NO-0315 Oslo, Norway
| | - Matthew J. McGrath
- Department of Biophysics,
Graduate
School of Science, Kyoto University, Kyoto
606-8502, Japan
| | - Theo Kurtén
- Division of Atmospheric Sciences,
Department of Physics, University of Helsinki, P.O. Box 64, FI-00014 University of Helsinki, Helsinki, Finland
| | - Ismael K. Ortega
- Division of Atmospheric Sciences,
Department of Physics, University of Helsinki, P.O. Box 64, FI-00014 University of Helsinki, Helsinki, Finland
| | - Hanna Vehkamäki
- Division of Atmospheric Sciences,
Department of Physics, University of Helsinki, P.O. Box 64, FI-00014 University of Helsinki, Helsinki, Finland
| | - Einar Uggerud
- Mass Spectrometry Laboratory
and Centre of Theoretical and Computational Chemistry, Department
of Chemistry, University of Oslo, P.O.
Box 1033 Blindern, NO-0315 Oslo, Norway
| |
Collapse
|
43
|
Haunschild R, Klopper W. New accurate reference energies for the G2/97 test set. J Chem Phys 2012; 136:164102. [DOI: 10.1063/1.4704796] [Citation(s) in RCA: 40] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
44
|
Haunschild R, Mao S, Mukherjee D, Klopper W. A universal explicit electron correlation correction applied to Mukherjee’s multi-reference perturbation theory. Chem Phys Lett 2012. [DOI: 10.1016/j.cplett.2012.02.020] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
|
45
|
Temelso B, Morrell TE, Shields RM, Allodi MA, Wood EK, Kirschner KN, Castonguay TC, Archer KA, Shields GC. Quantum Mechanical Study of Sulfuric Acid Hydration: Atmospheric Implications. J Phys Chem A 2012; 116:2209-24. [DOI: 10.1021/jp2119026] [Citation(s) in RCA: 98] [Impact Index Per Article: 8.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Berhane Temelso
- Dean’s Office, College
of Arts and Sciences, and Department of Chemistry, Bucknell University, Lewisburg, Pennsylvania
17837, United States
| | - Thomas E. Morrell
- Dean’s Office, College
of Arts and Sciences, and Department of Chemistry, Bucknell University, Lewisburg, Pennsylvania
17837, United States
| | - Robert M. Shields
- Dean’s Office, College
of Arts and Sciences, and Department of Chemistry, Bucknell University, Lewisburg, Pennsylvania
17837, United States
| | - Marco A. Allodi
- Dean’s Office, College
of Arts and Sciences, and Department of Chemistry, Bucknell University, Lewisburg, Pennsylvania
17837, United States
| | - Elena K. Wood
- Dean’s Office, College
of Arts and Sciences, and Department of Chemistry, Bucknell University, Lewisburg, Pennsylvania
17837, United States
| | - Karl N. Kirschner
- Department of Simulation
Engineering, Fraunhofer-Institute for Algorithms and Scientific Computing (SCAI), Schloss Birlinghoven,
53754 Sankt Augustin, Germany
| | - Thomas C. Castonguay
- Department of Chemistry, Iona College, New Rochelle, New York
10801, United States
| | - Kaye A. Archer
- Dean’s Office, College
of Arts and Sciences, and Department of Chemistry, Bucknell University, Lewisburg, Pennsylvania
17837, United States
| | - George C. Shields
- Dean’s Office, College
of Arts and Sciences, and Department of Chemistry, Bucknell University, Lewisburg, Pennsylvania
17837, United States
| |
Collapse
|
46
|
Haunschild R, Klopper W. Theoretical reference values for the AE6 and BH6 test sets from explicitly correlated coupled-cluster theory. Theor Chem Acc 2012. [DOI: 10.1007/s00214-012-1112-3] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/14/2022]
|
47
|
Hättig C, Klopper W, Köhn A, Tew DP. Explicitly Correlated Electrons in Molecules. Chem Rev 2011; 112:4-74. [DOI: 10.1021/cr200168z] [Citation(s) in RCA: 401] [Impact Index Per Article: 30.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Christof Hättig
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, D-44780 Bochum, Germany
| | - Wim Klopper
- Abteilung für Theoretische Chemie, Institut für Physikalische Chemie, Karlsruher Institut für Technologie, KIT-Campus Süd, Postfach 6980, D-76049 Karlsruhe, Germany
| | - Andreas Köhn
- Institut für Physikalische Chemie, Johannes Gutenberg-Universität Mainz, D-55099 Mainz, Germany
| | - David P. Tew
- School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| |
Collapse
|
48
|
Tew DP, Helmich B, Hättig C. Local explicitly correlated second-order Møller–Plesset perturbation theory with pair natural orbitals. J Chem Phys 2011; 135:074107. [DOI: 10.1063/1.3624370] [Citation(s) in RCA: 76] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023] Open
|
49
|
Bachorz RA, Bischoff FA, Glöss A, Hättig C, Höfener S, Klopper W, Tew DP. The MP2-F12 method in the Turbomole program package. J Comput Chem 2011; 32:2492-513. [PMID: 21590779 DOI: 10.1002/jcc.21825] [Citation(s) in RCA: 91] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2010] [Revised: 03/25/2011] [Accepted: 03/27/2011] [Indexed: 01/21/2023]
Abstract
A detailed description of the explicitly correlated second-order Møller-Plesset perturbation theory (MP2-F12) method, as implemented in the TURBOMOLE program package, is presented. The TURBOMOLE implementation makes use of density fitting, which greatly reduces the prefactor for integral evaluation. Methods are available for the treatment of ground states of open- and closed-shell species, using unrestricted as well as restricted (open-shell) Hartree-Fock reference determinants. Various methodological choices and approximations are discussed. The performance of the TURBOMOLE implementation is illustrated by example calculations of the molecules leflunomide, prednisone, methotrexate, ethylenedioxytetrafulvalene, and a cluster model for the adsorption of methanol on the zeolite H-ZSM-5. Various basis sets are used, including the correlation-consistent basis sets specially optimized for explicitly correlated calculations (cc-pVXZ-F12).
Collapse
Affiliation(s)
- Rafał A Bachorz
- Center for Functional Nanostructures (CFN, Karlsruhe Institute of Technology, Karlsruhe, Germany
| | | | | | | | | | | | | |
Collapse
|
50
|
Shiozaki T, Knizia G, Werner HJ. Explicitly correlated multireference configuration interaction: MRCI-F12. J Chem Phys 2011; 134:034113. [DOI: 10.1063/1.3528720] [Citation(s) in RCA: 198] [Impact Index Per Article: 15.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023] Open
|