1
|
Tsushima S, Kretzschmar J, Doi H, Okuwaki K, Kaneko M, Mochizuki Y, Takao K. Towards tailoring hydrophobic interaction with uranyl(VI) oxygen for C-H activation. Chem Commun (Camb) 2024; 60:4769-4772. [PMID: 38563824 DOI: 10.1039/d4cc01030b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/04/2024]
Abstract
Bovine serum albumin (BSA) has a uranyl(VI) binding hotspot where uranium is tightly bound by three carboxylates. Uranyl oxygen is "soaked" into the hydrophobic core of BSA. Isopropyl hydrogen of Val is trapped near UO22+ and upon photoexcitation, C-H bond cleavage is initiated. A unique hydrophobic contact with "yl"-oxygen, as observed here, can be used to induce C-H activation.
Collapse
Affiliation(s)
- Satoru Tsushima
- Institute of Resource Ecology, Helmholtz-Zentrum Dresden-Rossendorf (HZDR), Dresden, 01328, Germany.
- International Research Frontiers Initiative (IRFI), Institute of Innovative Research, Tokyo Institute of Technology, Tokyo, 152-8550, Japan
| | - Jérôme Kretzschmar
- Institute of Resource Ecology, Helmholtz-Zentrum Dresden-Rossendorf (HZDR), Dresden, 01328, Germany.
| | - Hideo Doi
- Department of Chemistry and Research Center for Smart Molecules, Rikkyo University, Tokyo, 171-8501, Japan
| | - Koji Okuwaki
- Department of Chemistry and Research Center for Smart Molecules, Rikkyo University, Tokyo, 171-8501, Japan
| | - Masashi Kaneko
- Department of Chemistry, Osaka University, Osaka, 560-0043, Japan
| | - Yuji Mochizuki
- Department of Chemistry, Osaka University, Osaka, 560-0043, Japan
- Institute of Industrial Science, The University of Tokyo, Tokyo, 153-8505, Japan
| | - Koichiro Takao
- Laboratory for Zero-Carbon Energy, Institute of Innovative Research, Tokyo Institute of Technology, Tokyo, 152-8550, Japan
| |
Collapse
|
2
|
Masumian E, Boese AD. Benchmarking Swaths of Intermolecular Interaction Components with Symmetry-Adapted Perturbation Theory. J Chem Theory Comput 2024; 20:30-48. [PMID: 38117939 PMCID: PMC10782453 DOI: 10.1021/acs.jctc.3c00801] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2023] [Revised: 11/19/2023] [Accepted: 11/22/2023] [Indexed: 12/22/2023]
Abstract
A benchmark database for interaction energy components of various noncovalent interactions (NCIs) along their dissociation curve is one of the essential needs in theoretical chemistry, especially for the development of force fields and machine-learning methods. We utilize DFT-SAPT or SAPT(DFT) as one of the most accurate methods to generate an extensive stock of the energy components, including dispersion energies extrapolated to the complete basis set limit (CBS). Precise analyses of the created data, and benchmarking the total interaction energies against the best available CCSD(T)/CBS values, reveal different aspects of the methodology and the nature of NCIs. For example, error cancellation effects between the S2 approximation and nonexact xc-potentials occur, and large charge transfer energies in some systems, including heavy atoms, can explain the lower accuracy of DFT-SAPT. This method is perfect for neutral complexes containing light nonmetals, while other systems with heavier atoms should be treated carefully. In the last part, a representative data set for all NCIs is extracted from the original data.
Collapse
Affiliation(s)
- Ehsan Masumian
- Physical and Theoretical Chemistry,
Department of Chemistry, University of Graz, 8010 Graz, Austria
| | - A. Daniel Boese
- Physical and Theoretical Chemistry,
Department of Chemistry, University of Graz, 8010 Graz, Austria
| |
Collapse
|
3
|
Wong J, Ganoe B, Liu X, Neudecker T, Lee J, Liang J, Wang Z, Li J, Rettig A, Head-Gordon T, Head-Gordon M. An in-silico NMR laboratory for nuclear magnetic shieldings computed via finite fields: Exploring nucleus-specific renormalizations of MP2 and MP3. J Chem Phys 2023; 158:164116. [PMID: 37114707 PMCID: PMC10148725 DOI: 10.1063/5.0145130] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/03/2023] [Accepted: 04/10/2023] [Indexed: 04/29/2023] Open
Abstract
We developed and implemented a method-independent, fully numerical, finite difference approach to calculating nuclear magnetic resonance shieldings, using gauge-including atomic orbitals. The resulting capability can be used to explore non-standard methods, given only the energy as a function of finite-applied magnetic fields and nuclear spins. For example, standard second-order Møller-Plesset theory (MP2) has well-known efficacy for 1H and 13C shieldings and known limitations for other nuclei such as 15N and 17O. It is, therefore, interesting to seek methods that offer good accuracy for 15N and 17O shieldings without greatly increased compute costs, as well as exploring whether such methods can further improve 1H and 13C shieldings. Using a small molecule test set of 28 species, we assessed two alternatives: κ regularized MP2 (κ-MP2), which provides energy-dependent damping of large amplitudes, and MP2.X, which includes a variable fraction, X, of third-order correlation (MP3). The aug-cc-pVTZ basis was used, and coupled cluster with singles and doubles and perturbative triples [CCSD(T)] results were taken as reference values. Our κ-MP2 results reveal significant improvements over MP2 for 13C and 15N, with the optimal κ value being element-specific. κ-MP2 with κ = 2 offers a 30% rms error reduction over MP2. For 15N, κ-MP2 with κ = 1.1 provides a 90% error reduction vs MP2 and a 60% error reduction vs CCSD. On the other hand, MP2.X with a scaling factor of 0.6 outperformed CCSD for all heavy nuclei. These results can be understood as providing renormalization of doubles amplitudes to partially account for neglected triple and higher substitutions and offer promising opportunities for future applications.
Collapse
Affiliation(s)
- Jonathan Wong
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Brad Ganoe
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Xiao Liu
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Tim Neudecker
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Joonho Lee
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Jiashu Liang
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Zhe Wang
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Jie Li
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Adam Rettig
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | | | | |
Collapse
|
4
|
Xie Y, Glick ZL, Sherrill CD. Assessment of three-body dispersion models against coupled-cluster benchmarks for crystalline benzene, carbon dioxide, and triazine. J Chem Phys 2023; 158:094110. [PMID: 36889937 DOI: 10.1063/5.0143712] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/09/2023] Open
Abstract
To study the contribution of three-body dispersion to crystal lattice energies, we compute the three-body contributions to the lattice energies for crystalline benzene, carbon dioxide, and triazine using various computational methods. We show that these contributions converge quickly as the intermolecular distances between the monomers grow. In particular, the smallest value among the three pairwise intermonomer closest-contact distances, Rmin, shows a strong correlation with the three-body contribution to the lattice energy, and, here, the largest of the closest-contact distances, Rmax, serves as a cutoff criterion to limit the number of trimers to be considered. We considered all trimers up to Rmax=15Å. The trimers with Rmin<4Å contribute 90.4%, 90.6%, and 93.9% of the total three-body contributions for crystalline benzene, carbon dioxide, and triazine, respectively, for the coupled-cluster singles, doubles, and perturbative triples [CCSD(T)] method. For trimers with Rmin>4Å, the second-order Møller-Plesset perturbation theory (MP2) supplemented with the Axilrod-Teller-Muto (ATM) three-body dispersion correction reproduces the CCSD(T) values for the cumulative three-body contributions with errors of less than 0.1 kJ mol-1. Moreover, three-body contributions are converged within 0.15 kJ mol-1 by Rmax=10Å. From these results, it appears that in molecular crystals where dispersion dominates the three-body contribution to the lattice energy, the trimers with Rmin>4Å can be computed with the MP2+ATM method to reduce the computational cost, and those with Rmax>10Å appear to be basically negligible.
Collapse
Affiliation(s)
- Yi Xie
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Zachary L Glick
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| |
Collapse
|
5
|
Aðalsteinsson HM, Bjornsson R. Ionization energies of metallocenes: a coupled cluster study of cobaltocene. Phys Chem Chem Phys 2023; 25:4570-4587. [PMID: 36723003 DOI: 10.1039/d2cp04715b] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
Open-shell transition metal chemistry presents challenges to contemporary electronic structure methods, based on either density functional or wavefunction theory. While CCSD(T) is the well-trusted gold standard for maingroup thermochemistry, the accuracy and robustness of the method is less clear for open-shell transition metal chemistry, requiring benchmarking of CCSD(T)-based protocols against either higher-level theory or experiment. Ionization energies (IEs) of metallocenes provide an interesting test case with metallocenes being common redox reagents as well as playing roles as redox mediators and cocatalysts in redox catalysis. Using highly accurate ZEKE-MATI experimental measurements of gas phase adiabatic (5.3275 ± 0.0006 eV) and vertical (5.4424 ± 0.0006 eV) ionization energies of cobaltocene, we systematically assessed the accuracy of the local coupled-cluster method DLPNO-CCSD(T) with respect to geometry, reference determinant, basis set size and extrapolation schemes, PNO cut-off and extrapolation, local triples approximation, relativistic effects and core-valence correlation. We show that PNO errors are controllable via the recently introduced PNO extrapolation schemes and that the expensive iterative triples (T1) contribution can be made more manageable by calculating it as a smaller-basis/smaller PNO-cutoff correction. The reference determinant turns out to be a critical aspect in these calculations with the HF determinant resulting in large DLPNO-CCSD(T) errors, likely due to the qualitatively flawed molecular orbital spectrum. The BP86 functional on the other hand was found to provide reference orbitals giving small DLPNO-CCSD(T) errors, likely due to more realistic orbitals as suggested by the more consistent MO spectrum compared to HF. A protocol including complete basis set extrapolations with correlation-consistent basis sets, complete PNO space extrapolations, iterative triples- and core-valence correlation corrections was found to give errors of -0.07 eV and -0.03 eV for adiabatic- and vertical-IE of cobaltocene, respectively, giving close to chemical accuracy for both properties. A computationally efficient DLPNO-CCSD(T) protocol was devised and tested against adiabatic ionization energies of 6 different metallocenes (V, Cr, Mn, Fe, Co, Ni). For the other metallocenes, the iterative triples (T1) and PNO extrapolation contributions turn out to be even more important. The results give errors close to the experimental uncertainty, similar to recent auxiliary-field quantum Monte Carlo results. The quality of the reference determinant orbitals is identified as the main source of uncertainty in CCSD(T) calculations of metallocenes.
Collapse
Affiliation(s)
| | - Ragnar Bjornsson
- Science Institute, University of Iceland, 107 Reykjavik, Iceland.,Univ Grenoble Alpes, CNRS, CEA, IRIG, Laboratoire de Chimie et Biologie des Métaux, 17 Rue des Martyrs, F-38054 Grenoble Cedex, France.
| |
Collapse
|
6
|
Sargent CT, Metcalf DP, Glick ZL, Borca CH, Sherrill CD. Benchmarking two-body contributions to crystal lattice energies and a range-dependent assessment of approximate methods. J Chem Phys 2023; 158:054112. [PMID: 36754814 DOI: 10.1063/5.0141872] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023] Open
Abstract
Using the many-body expansion to predict crystal lattice energies (CLEs), a pleasantly parallel process, allows for flexibility in the choice of theoretical methods. Benchmark-level two-body contributions to CLEs of 23 molecular crystals have been computed using interaction energies of dimers with minimum inter-monomer separations (i.e., closest contact distances) up to 30 Å. In a search for ways to reduce the computational expense of calculating accurate CLEs, we have computed these two-body contributions with 15 different quantum chemical levels of theory and compared these energies to those computed with coupled-cluster in the complete basis set (CBS) limit. Interaction energies of the more distant dimers are easier to compute accurately and several of the methods tested are suitable as replacements for coupled-cluster through perturbative triples for all but the closest dimers. For our dataset, sub-kJ mol-1 accuracy can be obtained when calculating two-body interaction energies of dimers with separations shorter than 4 Å with coupled-cluster with single, double, and perturbative triple excitations/CBS and dimers with separations longer than 4 Å with MP2.5/aug-cc-pVDZ, among other schemes, reducing the number of dimers to be computed with coupled-cluster by as much as 98%.
Collapse
Affiliation(s)
- Caroline T Sargent
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Derek P Metcalf
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Zachary L Glick
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Carlos H Borca
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| |
Collapse
|
7
|
Förster A. Assessment of the Second-Order Statically Screened Exchange Correction to the Random Phase Approximation for Correlation Energies. J Chem Theory Comput 2022; 18:5948-5965. [PMID: 36150190 PMCID: PMC9558381 DOI: 10.1021/acs.jctc.2c00366] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
With increasing interelectronic distance, the screening
of the
electron–electron interaction by the presence of other electrons
becomes the dominant source of electron correlation. This effect is
described by the random phase approximation (RPA) which is therefore
a promising method for the calculation of weak interactions. The success
of the RPA relies on the cancellation of errors, which can be traced
back to the violation of the crossing symmetry of the 4-point vertex,
leading to strongly overestimated total correlation energies. By the
addition of second-order screened exchange (SOSEX) to the correlation
energy, this issue is substantially reduced. In the adiabatic connection
(AC) SOSEX formalism, one of the two electron–electron interaction
lines in the second-order exchange term is dynamically screened (SOSEX(W, vc)). A
related SOSEX expression in which both electron–electron interaction
lines are statically screened (SOSEX(W(0), W(0))) is obtained from the G3W2 contribution to the electronic self-energy. In contrast to SOSEX(W, vc), the
evaluation of this correlation energy expression does not require
an expensive numerical frequency integration and is therefore advantageous
from a computational perspective. We compare the accuracy of the statically
screened variant to RPA and RPA+SOSEX(W, vc) for a wide range of chemical
reactions. While both methods fail for barrier heights, SOSEX(W(0), W(0)) agrees very well with SOSEX(W, vc) for
charged excitations and noncovalent interactions where they lead to
major improvements over RPA.
Collapse
Affiliation(s)
- Arno Förster
- Theoretical Chemistry, Vrije Universiteit, De Boelelaan 1083, NL-1081 HV, Amsterdam, The Netherlands
| |
Collapse
|
8
|
Behnle S, Richter R, Völkl L, Idzko P, Förstner A, Bozkaya U, Fink RF. Accurate Property Prediction by Second Order Perturbation Theory: The REMP and OO-REMP Hybrids. J Chem Phys 2022; 157:104111. [DOI: 10.1063/5.0105628] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The prediction of molecular properties such as equilibrium structures or vibrationalwavenumbers is a routine task in computational chemistry. If very high accuracy is required, however, the use of computationally demanding ab initio wavefunction methods is mandatory. We present property calculations utilizing the REMP and OO-REMP hybrid perturbation theories showing that with the latter approach, very accurate results are obtained at second order in perturbation theory. Specifically, equilibrium structures and harmonic vibrational wavenumbers as well as dipole moments of closed and open shell molecules were calculated and compared to the best available experimental results or very accurate calculations.OO-REMP is capable of predicting bond lengths of small closed and open shell molecules with an accuracy of 0.2 pm and 0.5 pm, respectively, often within the range of experimental uncertainty. Equilibrium harmonic vibrational wavenumbers are predicted with an accuracy better than 20 cm−1 . Dipole moments of small closed and open shell molecules are reproduced with a relative error of less than 3 %. Across all investigated properties it turns out that a 20 %:80 % MP:RE mixing ratio consistently provides the best results. This is in line with our previous findings featuring closed and open shell reaction energies.
Collapse
Affiliation(s)
- Stefan Behnle
- Fachbereich II Chemie, Eberhard Karls Universität Tübingen Fachbereich II Chemie, Germany
| | - Robert Richter
- Eberhard Karls Universitat Tubingen Fachbereich II Chemie, Germany
| | - Luca Völkl
- Eberhard Karls Universitat Tubingen Fachbereich II Chemie, Germany
| | - Paul Idzko
- Eberhard Karls Universitat Tubingen Fachbereich II Chemie, Germany
| | - André Förstner
- Eberhard Karls Universitat Tubingen Fachbereich II Chemie, Germany
| | - Uğur Bozkaya
- Department of Chemistry, Hacettepe University, Turkey
| | - Reinhold F Fink
- Institute of Physical and Theoretical Chemistry, Eberhard Karls Universität Tübingen Fachbereich II Chemie, Germany
| |
Collapse
|
9
|
Huang HH, Wang YS, Chao SD. A Minimum Quantum Chemistry CCSD(T)/CBS Data Set of Dimeric Interaction Energies for Small Organic Functional Groups: Heterodimers. ACS OMEGA 2022; 7:20059-20080. [PMID: 35722020 PMCID: PMC9201891 DOI: 10.1021/acsomega.2c01888] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/28/2022] [Accepted: 05/04/2022] [Indexed: 06/15/2023]
Abstract
We extend our previous quantum chemistry calculations of interaction energies for 31 homodimers of small organic functional groups (the SOFG-31 data set) by including 239 heterodimers with monomers selected within the SOFG-31 data set, thus resulting in the SOFG-31+239 data set. The minimum-level theoretical scheme contains (1) the basis set superposition error corrected supermolecule (BSSE-SM) approach for intermolecular interactions; (2) the second-order Møller-Plesset perturbation theory (MP2) with the Dunning's aug-cc-pVXZ (X = D, T, Q) basis sets for the geometry optimization and correlation energy calculations; and (3) the single-point energy calculations with the coupled cluster with single, double, and perturbative triple excitations method at the complete basis set limit [CCSD(T)/CBS] using the well-tested extrapolation methods for the MP2 energy calibrations. In addition, we have performed a parallel series of energy decomposition calculations based on the symmetry adapted perturbation theory (SAPT) in order to gain chemical insights. That the above procedure cannot be further reduced has been proven to be very crucial for constructing reliable data sets of interaction energies. The calculated CCSD(T)/CBS interaction energy data can serve as a benchmark for testing or training less accurate but more efficient calculation methods, such as the electronic density functional theory. As an application, we employ a segmental SAPT model previously developed for the SOFG-31 data set to predict binding energies of large heterodimer complexes. These model energy "quanta" can be used in coarse-grained molecular dynamics simulations by avoiding large-scale calculations.
Collapse
Affiliation(s)
- Hsing-Hsiang Huang
- Institute
of Applied Mechanics, National Taiwan University, Taipei 10617, Taiwan R.O.C.
| | - Yi-Siang Wang
- School
of Chemistry & Biochemistry, Georgia
Institute of Technology, Atlanta, Georgia 30332, United States
| | - Sheng D. Chao
- Institute
of Applied Mechanics, National Taiwan University, Taipei 10617, Taiwan R.O.C.
| |
Collapse
|
10
|
Boggio-Pasqua M, Jacquemin DM, Loos PF. Benchmarking CASPT3 Vertical Excitation Energies. J Chem Phys 2022; 157:014103. [DOI: 10.1063/5.0095887] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Based on 280 reference vertical transition energies of various natures (singlet, triplet, valence, Rydberg, n → π∗, π → π∗, and double excitations) extracted from the QUEST database, we assess the accuracy of third-order multireference perturbation theory, CASPT3, in the context of molecular excited states. When one applies the disputable ionization- potential-electron-affinity (IPEA) shift, we show that CASPT3 provides a similar accuracy as its second-order counterpart, CASPT2, with the same mean absolute error of 0.11 eV. However, as already reported, we also observe that the accuracy of CASPT3 is almost insensitive to the IPEA shift, irrespective of the transition type and system size, with a small reduction of the mean absolute error to 0.09 eV when the IPEA shift is switched off.
Collapse
Affiliation(s)
| | - Denis M. Jacquemin
- Chimie Et Interdisciplinarité, Synthèse, Analyse, Modélisation, University of Nantes, France
| | | |
Collapse
|
11
|
Ehlert S, Grimme S, Hansen A. Conformational Energy Benchmark for Longer n-Alkane Chains. J Phys Chem A 2022; 126:3521-3535. [PMID: 35616628 DOI: 10.1021/acs.jpca.2c02439] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present the first benchmark set focusing on the conformational energies of highly flexible, long n-alkane chains, termed ACONFL. Unbranched alkanes are ubiquitous building blocks in nature, so the goal is to be able to calculate their properties most accurately to improve the modeling of, e.g., complex (biological) systems. Very accurate DLPNO-CCSD(T1)/CBS reference values are provided, which allow for a statistical meaningful evaluation of even the best available density functional methods. The performance of established and modern (dispersion corrected) density functionals is comprehensively assessed. The recently introduced r2SCAN-V functional shows excellent performance, similar to efficient composite DFT methods like B97-3c and r2SCAN-3c, which provide an even better cost-accuracy ratio, while almost reaching the accuracy of much more computationally demanding hybrid or double hybrid functionals with large QZ AO basis sets. In addition, we investigated the performance of common wave function methods, where MP2/CBS surprisingly performs worse compared to the simple D4 dispersion corrected Hartree-Fock. Furthermore, we investigate the performance of several semiempirical and force field methods, which are commonly used for the generation of conformational ensembles in multilevel workflows or in large scale molecular dynamics studies. Outstanding performance is obtained by the recently introduced general force field, GFN-FF, while other commonly applied methods like the universal force field yield large errors. We recommend the ACONFL as a helpful benchmark set for parametrization of new semiempirical or force field methods and machine learning potentials as well as a meaningful validation set for newly developed DFT or dispersion methods.
Collapse
Affiliation(s)
- Sebastian Ehlert
- Mulliken Center for Theoretical Chemistry, Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstrasse 4, 53115 Bonn, Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstrasse 4, 53115 Bonn, Germany
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstrasse 4, 53115 Bonn, Germany
| |
Collapse
|
12
|
Sirianni DA, Zhu X, Sitkoff DF, Cheney DL, Sherrill CD. The influence of a solvent environment on direct non-covalent interactions between two molecules: A symmetry-adapted perturbation theory study of polarization tuning of π-π interactions by water. J Chem Phys 2022; 156:194306. [PMID: 35597646 DOI: 10.1063/5.0087302] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
High-level quantum chemical computations have provided significant insight into the fundamental physical nature of non-covalent interactions. These studies have focused primarily on gas-phase computations of small van der Waals dimers; however, these interactions frequently take place in complex chemical environments, such as proteins, solutions, or solids. To better understand how the chemical environment affects non-covalent interactions, we have undertaken a quantum chemical study of π-π interactions in an aqueous solution, as exemplified by T-shaped benzene dimers surrounded by 28 or 50 explicit water molecules. We report interaction energies (IEs) using second-order Møller-Plesset perturbation theory, and we apply the intramolecular and functional-group partitioning extensions of symmetry-adapted perturbation theory (ISAPT and F-SAPT, respectively) to analyze how the solvent molecules tune the π-π interactions of the solute. For complexes containing neutral monomers, even 50 explicit waters (constituting a first and partial second solvation shell) change total SAPT IEs between the two solute molecules by only tenths of a kcal mol-1, while significant changes of up to 3 kcal mol-1 of the electrostatic component are seen for the cationic pyridinium-benzene dimer. This difference between charged and neutral solutes is attributed to large non-additive three-body interactions within solvated ion-containing complexes. Overall, except for charged solutes, our quantum computations indicate that nearby solvent molecules cause very little "tuning" of the direct solute-solute interactions. This indicates that differences in binding energies between the gas phase and solution phase are primarily indirect effects of the competition between solute-solute and solute-solvent interactions.
Collapse
Affiliation(s)
- Dominic A Sirianni
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Xiao Zhu
- Molecular Structure and Design, Bristol Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - Doree F Sitkoff
- Molecular Structure and Design, Bristol Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - Daniel L Cheney
- Molecular Structure and Design, Bristol Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| |
Collapse
|
13
|
Bauer MM, Dempwolff AL, Rehn DR, Dreuw A. Exploring the accuracy and usefulness of semi-empirically scaled ADC schemes by blending second and third order terms. J Chem Phys 2022; 156:144101. [DOI: 10.1063/5.0087173] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Different approaches to mixed-order algebraic-diagrammatic construction (ADC) schemes are investigated. The performance of two different strategies for scaling third-order contributions to the ADC secular matrix is evaluated. Both considered schemes employ a single tuning parameter and conserve general properties inherent to all ADC methods such as hermiticity and size-consistency.The first approach, scaled-matrix ADC[(2)+x(3)], scales all contributions first occurring in ADC(3) equally and leads to an improvement of the accuracy of excitation energies compared to ADC(3) for x=0.4-0.5. However, with respect to excited state dipole moments this method provides lower accuracy than ADC(3). The second scaling approach, MP[(1)+x(2)]-ISR(3), scales the second order contributions of the ground-state wavefunction and derives a rigorous ADC scheme via the intermediate state representation (ISR) formalism. Although the error in excitation energies is not improved, this method provides insight into the relevance of the individual terms of the ADC(3) matrix and indicates, that the MP(2) wavefunction is indeed the optimal reference wavefunction for deriving a third-order single-reference ADC scheme.
Collapse
Affiliation(s)
- Marco M. Bauer
- Ruprecht Karls University Heidelberg Interdisciplinary Center for Scientific Computing, Germany
| | - Adrian Ludwig Dempwolff
- Ruprecht Karls University Heidelberg Interdisciplinary Center for Scientific Computing, Germany
| | | | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Ruprecht Karls University Heidelberg Interdisciplinary Center for Scientific Computing, Germany
| |
Collapse
|
14
|
Behnle S, Fink RF. UREMP, RO-REMP, and OO-REMP: Hybrid perturbation theories for open-shell electronic structure calculations. J Chem Phys 2022; 156:124103. [DOI: 10.1063/5.0081285] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
An accurate description of the electron correlation energy in closed- and open-shell molecules is shown to be obtained by a second-order perturbation theory (PT) termed REMP. REMP is a hybrid of the Retaining the Excitation degree (RE) and the Møller–Plesset (MP) PTs. It performs particularly encouragingly in an orbital-optimized variant (OO-REMP) where the reference wavefunction is given by an unrestricted Slater determinant whose spin orbitals are varied such that the total energy becomes a minimum. While the approach generally behaves less satisfactorily with unrestricted Hartree–Fock references, reasonable performance is observed for restricted Hartree–Fock and restricted open-shell Hartree–Fock references. Inclusion of single excitations to OO-REMP is investigated and found—as in similar investigations—to be dissatisfying as it deteriorates performance. For the non-multireference subset of the accurate W4-11 benchmark set of Karton et al. [Chem. Phys. Lett. 510, 165–178 (2011)], OO-REMP predicts most atomization and reaction energies with chemical accuracy (1 kcal mol−1) if complete-basis-set extrapolation with augmented and core-polarized basis sets is used. For the W4-11 related test-sets, the error estimates obtained with the OO-REMP method approach those of coupled-cluster with singles, doubles and perturbative triples [CCSD(T)] within 20%–35%. The best performance of OO-REMP is found for a mixing ratio of 20%:80% MP:RE, which is essentially independent of whether radical stabilization energies, barrier heights, or reaction energies are investigated. Orbital optimization is shown to improve the REMP approach for both closed and open shell cases and outperforms coupled-cluster theory with singles and doubles (CCSD), spin-component scaled Møller-Plesset theory at second order (SCS-MP2), and density functionals, including double hybrids in all the cases considered.
Collapse
Affiliation(s)
- Stefan Behnle
- Institute for Physical and Theoretical Chemistry, Eberhard Karls University Tübingen, Auf der Morgenstelle 18, 72076 Tübingen, Germany
| | - Reinhold F. Fink
- Institute for Physical and Theoretical Chemistry, Eberhard Karls University Tübingen, Auf der Morgenstelle 18, 72076 Tübingen, Germany
| |
Collapse
|
15
|
Rozhenko AB, Horbenko YS, Kyrylchuk AA, Zarudnitskii EV, Mykhaylychenko SS, Shermolovich YG, Grafov AV. Stable Carbenes as Structural Components of Partially Saturated Sulfur-Containing Heterocycles. Molecules 2022; 27:molecules27051458. [PMID: 35268558 PMCID: PMC8911733 DOI: 10.3390/molecules27051458] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/24/2021] [Revised: 02/12/2022] [Accepted: 02/18/2022] [Indexed: 11/16/2022] Open
Abstract
Recently, an unusual elongation of the C-S bond was observed experimentally for some sulfur-containing heterocycles. Using a superior ab initio (SCS-MP2/cc-pVTZ) level of theory, we showed that the phenomenon can be explained by a contribution of a donor-acceptor adduct of a carbene with an unsaturated ligand. One may achieve further elongation of the C-S bond, eventually turning it to a coordinate one, by increasing the stability of each part of the system as, e.g., in the utmost case of spiro adducts with Arduengo carbenes. The effect of carbene stability was quantified by employing the isodesmic reactions of carbene exchange.
Collapse
Affiliation(s)
- Alexander B. Rozhenko
- Institute of Organic Chemistry, National Academy of Sciences, Murmanska Str. 5, 02094 Kyiv, Ukraine; (A.B.R.); (Y.S.H.); (A.A.K.); (E.V.Z.); (S.S.M.); (Y.G.S.)
- Igor Sikorsky Kyiv Polytechnic Institute, National Technical University of Ukraine, Prosp. Peremohy 37, 03056 Kyiv, Ukraine
| | - Yuliia S. Horbenko
- Institute of Organic Chemistry, National Academy of Sciences, Murmanska Str. 5, 02094 Kyiv, Ukraine; (A.B.R.); (Y.S.H.); (A.A.K.); (E.V.Z.); (S.S.M.); (Y.G.S.)
- Igor Sikorsky Kyiv Polytechnic Institute, National Technical University of Ukraine, Prosp. Peremohy 37, 03056 Kyiv, Ukraine
| | - Andrii A. Kyrylchuk
- Institute of Organic Chemistry, National Academy of Sciences, Murmanska Str. 5, 02094 Kyiv, Ukraine; (A.B.R.); (Y.S.H.); (A.A.K.); (E.V.Z.); (S.S.M.); (Y.G.S.)
| | - Evgenij V. Zarudnitskii
- Institute of Organic Chemistry, National Academy of Sciences, Murmanska Str. 5, 02094 Kyiv, Ukraine; (A.B.R.); (Y.S.H.); (A.A.K.); (E.V.Z.); (S.S.M.); (Y.G.S.)
| | - Sergiy S. Mykhaylychenko
- Institute of Organic Chemistry, National Academy of Sciences, Murmanska Str. 5, 02094 Kyiv, Ukraine; (A.B.R.); (Y.S.H.); (A.A.K.); (E.V.Z.); (S.S.M.); (Y.G.S.)
| | - Yuriy G. Shermolovich
- Institute of Organic Chemistry, National Academy of Sciences, Murmanska Str. 5, 02094 Kyiv, Ukraine; (A.B.R.); (Y.S.H.); (A.A.K.); (E.V.Z.); (S.S.M.); (Y.G.S.)
| | - Andriy V. Grafov
- Department of Chemistry, University of Helsinki, A.I.Virtasen Aukio 1, 00560 Helsinki, Finland
- Correspondence:
| |
Collapse
|
16
|
Krstić M, Fink K, Sharapa DI. The Adsorption of Small Molecules on the Copper Paddle-Wheel: Influence of the Multi-Reference Ground State. MOLECULES (BASEL, SWITZERLAND) 2022; 27:molecules27030912. [PMID: 35164179 PMCID: PMC8840508 DOI: 10.3390/molecules27030912] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/31/2021] [Revised: 01/26/2022] [Accepted: 01/26/2022] [Indexed: 11/16/2022]
Abstract
We report a theoretical study of the adsorption of a set of small molecules (C2H2, CO, CO2, O2, H2O, CH3OH, C2H5OH) on the metal centers of the “copper paddle-wheel”—a key structural motif of many MOFs. A systematic comparison between DFT of different rungs, single-reference post-HF methods (MP2, SOS–MP2, MP3, DLPNO–CCSD(T)), and multi-reference approaches (CASSCF, DCD–CAS(2), NEVPT2) is performed in order to find a methodology that correctly describes the complicated electronic structure of paddle-wheel structure together with a reasonable description of non-covalent interactions. Apart from comparison with literature data (experimental values wherever possible), benchmark calculations with DLPNO–MR–CCSD were also performed. Despite tested methods show qualitative agreement in the majority of cases, we showed and discussed reasons for quantitative differences as well as more fundamental problems of specific cases.
Collapse
Affiliation(s)
- Marjan Krstić
- Institute for Theoretical Solid State Physics (TFP), Karlsruhe Institute of Technology (KIT), Wolfgang-Gaede-Str. 1, 76131 Karlsruhe, Germany;
| | - Karin Fink
- Institute of Nanotechnology (INT), Karlsruhe Institute of Technology (KIT), Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany;
| | - Dmitry I. Sharapa
- Institute of Catalysis Research and Technology (IKFT), Karlsruhe Institute of Technology (KIT), Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany
- Correspondence:
| |
Collapse
|
17
|
Ermiş B, Ekinci E, Bozkaya U. State-Of-The-Art Computations of Vertical Electron Affinities with the Extended Koopmans' Theorem Integrated with the CCSD(T) Method. J Chem Theory Comput 2021; 17:7648-7656. [PMID: 34724787 DOI: 10.1021/acs.jctc.1c00938] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Accurate computation of electron affinities (EAs), within 0.10 eV, is one of the most challenging problems in modern computational quantum chemistry. The extended Koopmans' theorem (EKT) enables direct computations of electron affinities (EAs) from any level of the theory. In this research, the EKT approach based on the coupled-cluster singles and doubles with perturbative triples [CCSD(T)] method is applied to computations of EAs for the first time. For efficiency, the density-fitting (DF) technique is used for two-electron integrals. Further, the EKT-CCSD(T) method is applied to three test sets of atoms and closed- and open-shell molecules, denoted A16, C10, and O33, respectively, for comparison with the experimental electron affinities. For the A16, C10, and O33 sets, the EKT-CCSD(T) approach, along with the aug-cc-pV5Z basis set, provide mean absolute errors (MAEs) of 0.05, 0.08, and 0.09 eV, respectively. Hence, our results demonstrate that high-accuracy computations of EAs can be achieved with the EKT-CCSD(T) approach. Further, when the EKT-CCSD(T) approach is not computationally affordable, the EKT-MP2.5, EKT-LCCD, and EKT-CCSD methods can be considered, and their results are also reasonably accurate. The huge advantage of the EKT method for the computation of IPs is that it comes for free in an analytic gradient computation. Hence, one needs neither separate computations for neutral and ionic species, as in the case of common approaches, nor additional efforts to obtain IPs, as in the case of equation-of-motion approaches. Overall, we believe that the present research may open new avenues in EA computations.
Collapse
Affiliation(s)
- Betül Ermiş
- Department of Chemistry, Hacettepe University, Ankara 06800, Turkey
| | - Emel Ekinci
- Department of Chemistry, Hacettepe University, Ankara 06800, Turkey
| | - Uğur Bozkaya
- Department of Chemistry, Hacettepe University, Ankara 06800, Turkey
| |
Collapse
|
18
|
Santra G, Semidalas E, Martin JML. Surprisingly Good Performance of XYG3 Family Functionals Using a Scaled KS-MP3 Correlation. J Phys Chem Lett 2021; 12:9368-9376. [PMID: 34550706 DOI: 10.1021/acs.jpclett.1c02838] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/16/2023]
Abstract
By adding a GLPT3 (third-order Görling-Levy perturbation theory, or KS-MP3) term E3 to the XYG7 form for a double hybrid, we are able to bring down WTMAD2 (weighted total mean absolute deviation) for the very large and chemically diverse GMTKN55 benchmark to an unprecedented 1.17 kcal/mol, competitive with much costlier composite wave function ab initio approaches. Intriguingly, (a) the introduction of E3 makes an empirical dispersion correction redundant; (b) generalized gradient approximation (GGA) or meta-GGA semilocal correlation functionals offer no advantage over the local density approximation (LDA) in this framework; (c) if a dispersion correction is retained, then simple Slater exchange leads to no significant loss in accuracy. It is possible to create a six-parameter functional with WTMAD2 = 1.42 that has no post-LDA density functional theory components and no dispersion correction in the final energy.
Collapse
Affiliation(s)
- Golokesh Santra
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| | - Emmanouil Semidalas
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| | - Jan M L Martin
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| |
Collapse
|
19
|
Loipersberger M, Bertels LW, Lee J, Head-Gordon M. Exploring the Limits of Second- and Third-Order Møller-Plesset Perturbation Theories for Noncovalent Interactions: Revisiting MP2.5 and Assessing the Importance of Regularization and Reference Orbitals. J Chem Theory Comput 2021; 17:5582-5599. [PMID: 34382394 PMCID: PMC9948597 DOI: 10.1021/acs.jctc.1c00469] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
This work systematically assesses the influence of reference orbitals, regularization, and scaling on the performance of second- and third-order Møller-Plesset perturbation theory wave function methods for noncovalent interactions (NCIs). Testing on 19 data sets (A24, DS14, HB15, HSG, S22, X40, HW30, NC15, S66, AlkBind12, CO2Nitrogen16, HB49, Ionic43, TA13, XB18, Bauza30, CT20, XB51, and Orel26rad) covers a wide range of different NCIs including hydrogen bonding, dispersion, and halogen bonding. Inclusion of potential energy surfaces from different hydrogen bonds and dispersion-bound complexes gauges accuracy for nonequilibrium geometries. Fifteen methods are tested. In notation where nonstandard choices of orbitals are denoted as methods:orbitals, these are MP2, κ-MP2, SCS-MP2, OOMP2, κ-OOMP2, MP3, MP2.5, MP3:OOMP2, MP2.5:OOMP2, MP3:κ-OOMP2, MP2.5:κ-OOMP2, κ-MP3:κ-OOMP2, κ-MP2.5:κ-OOMP2, MP3:ωB97X-V, and MP2.5:ωB97X-V. Furthermore, we compare these methods to the ωB97M-V and B3LYP-D3 density functionals, as well as CCSD. We find that the κ-regularization (κ = 1.45 au was used throughout) improves the energetics in almost all data sets for both MP2 (in 17 out of 19 data sets) and OOMP2 (16 out of 19). The improvement is significant (e.g., the root-mean-square deviation (RMSD) for the S66 data set is 0.29 kcal/mol for κ-OOMP2 versus 0.67 kcal/mol for MP2) and for interactions between stable closed-shell molecules, not strongly dependent on the reference orbitals. Scaled MP3 (with a factor of 0.5) using κ-OOMP2 reference orbitals (MP2.5:κ-OOMP2) provides significantly more accurate results for NCIs across all data sets with noniterative O(N6) scaling (S66 data set RMSD: 0.10 kcal/mol). Across the entire data set of 356 points, the improvement over standard MP2.5 is approximately a factor of 2: RMSD for MP3:κ-OOMP2 is 0.25 vs 0.50 kcal/mol for MP2.5. The use of high-quality density functional reference orbitals (ωB97X-V) also significantly improves the results of MP2.5 for NCI over a Hartree-Fock orbital reference. All our assessments and conclusions are based on the use of the medium-sized aug-cc-pVTZ basis to yield results that are directly compared against complete basis set limit reference values.
Collapse
Affiliation(s)
| | - Luke W. Bertels
- Department of Chemistry, University of California, Berkeley, California 94720, USA,Present Address: Department of Chemistry, Virginia Tech, Blacksburg, VA 24061, USA
| | - Joonho Lee
- Department of Chemistry, University of California, Berkeley, California 94720, USA,Present Address: Department of Chemistry, Columbia University, NY
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, USA,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
20
|
Ahirwar MB, Gurav ND, Gadre SR, Deshmukh MM. Molecular Tailoring Approach for Estimating Individual Intermolecular Interaction Energies in Benzene Clusters. J Phys Chem A 2021; 125:6131-6140. [PMID: 34251827 DOI: 10.1021/acs.jpca.1c03907] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
There is no general method available for the estimation of individual intermolecular interaction energies in weakly bound molecular clusters, and such studies are limited only to the dimer. Recently, we proposed a molecular tailoring approach-based method for the estimation of individual O-H···O hydrogen bond energies in water clusters. In the present work, we extend the applicability of this method for estimating the individual intermolecular interaction energies in benzene clusters, which are expected to be small. The basis set superposition error (BSSE)-corrected individual intermolecular interaction energies in linear (LN) benzene clusters, LN-(Bz)n n = 3-7, were calculated to be in the range from -1.75 to -2.33 kcal/mol with the cooperativity contribution falling between 0.05 and 0.20 kcal/mol, calculated at the MP2.5/aug-cc-pVDZ level of theory. In the case of non-linear (NLN) benzene clusters, NLN-(Bz)n n = 3-5, the BSSE-corrected individual intermolecular interaction energies exhibit a wider range from -1.16 to -2.55 kcal/mol with cooperativity contribution in the range from 0.02 to -0.61 kcal/mol. The accuracy of these estimated values was validated by adding the sum of interaction energies to the sum of monomer energies. These estimated molecular energies of clusters were compared with their actual calculated values. The small difference (<0.3 kcal/mol) in these two values suggests that our estimated individual intermolecular interaction energies in benzene clusters are quite reliable.
Collapse
Affiliation(s)
- Mini Bharati Ahirwar
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya (A Central University), Sagar 470003, India
| | - Nalini D Gurav
- Department of Scientific Computing, Modelling and Simulation, Savitribai Phule Pune University, Pune 411 007, India
| | - Shridhar R Gadre
- Department of Scientific Computing, Modelling and Simulation, Savitribai Phule Pune University, Pune 411 007, India
| | - Milind M Deshmukh
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya (A Central University), Sagar 470003, India
| |
Collapse
|
21
|
Behnle S, Fink RF. OO-REMP: Approaching Chemical Accuracy with Second-Order Perturbation Theory. J Chem Theory Comput 2021; 17:3259-3266. [PMID: 34006110 DOI: 10.1021/acs.jctc.1c00280] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We present a perturbation theory (PT) providing second-order energies that reproduce main group chemistry benchmark sets for reaction energies, barrier heights, and atomization energies with mean absolute deviations below 1 kcal mol-1. The PT is defined as a constrained mixture of the unperturbed Hamiltonians of the Retaining the Excitation degree (RE) and the Møller-Plesset (MP) PTs. The orbitals of the reference wave function, a single unrestricted Slater determinant, are iteratively optimized to minimize the total energy. For all benchmark sets, good and near optimal performance of OO-REMP was observed for an unperturbed Hamiltonian consisting of 25% MP and 75% RE contributions.
Collapse
Affiliation(s)
- Stefan Behnle
- Eberhard Karls Universität, Auf der Morgenstelle 18, 72076 Tübingen, Germany
| | - Reinhold F Fink
- Eberhard Karls Universität, Auf der Morgenstelle 18, 72076 Tübingen, Germany
| |
Collapse
|
22
|
Ye HZ, Tran HK, Van Voorhis T. Accurate Electronic Excitation Energies in Full-Valence Active Space via Bootstrap Embedding. J Chem Theory Comput 2021; 17:3335-3347. [PMID: 33957050 DOI: 10.1021/acs.jctc.0c01221] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Fragment embedding has been widely used to circumvent the high computational scaling of using accurate electron correlation methods to describe the electronic ground states of molecules and materials. However, similar applications that utilize fragment embedding to treat electronic excited states are comparably less reported in the literature. The challenge here is twofold. First, most fragment embedding methods are most effective when the property of interest is local, but the change of the wave function upon excitation is nonlocal in general. Second, even for local excitations, an accurate estimate of, for example, the excitation energy can still be challenging owing to the need for a balanced treatment of both the ground and the excited states. In this work, we show that bootstrap embedding (BE), a fragment embedding method developed recently by our group, is promising toward describing general electronic excitations. Numerical simulations show that the excitation energies in full-valence active space (FVAS) can be well-estimated by BE to an error of ∼0.05 eV using relatively small fragments, for both local excitations and the excitations of some large dye molecules that exhibit strong charge-transfer characters. We hence anticipate BE to be a promising solution to accurately describing the excited states of large chemical systems.
Collapse
Affiliation(s)
- Hong-Zhou Ye
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Henry K Tran
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Troy Van Voorhis
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
23
|
Karton A, Martin JML. Prototypical π-π dimers re-examined by means of high-level CCSDT(Q) composite ab initio methods. J Chem Phys 2021; 154:124117. [PMID: 33810692 DOI: 10.1063/5.0043046] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023] Open
Abstract
The benzene-ethene and parallel-displaced (PD) benzene-benzene dimers are the most fundamental systems involving π-π stacking interactions. Several high-level ab initio investigations calculated the binding energies of these dimers using the coupled-cluster with singles, doubles, and quasi-perturbative triple excitations [CCSD(T)] method at the complete basis set [CBS] limit using various approaches such as reduced virtual orbital spaces and/or MP2-based basis set corrections. Here, we obtain CCSDT(Q) binding energies using a Weizmann-3-type approach. In particular, we extrapolate the self-consistent field (SCF), CCSD, and (T) components using large heavy-atom augmented Gaussian basis sets [namely, SCF/jul-cc-pV{5,6}Z, CCSD/jul-cc-pV{Q,5}Z, and (T)/jul-cc-pV{T,Q}Z]. We consider post-CCSD(T) contributions up to CCSDT(Q), inner-shell, scalar-relativistic, and Born-Oppenheimer corrections. Overall, our best relativistic, all-electron CCSDT(Q) binding energies are ∆Ee,all,rel = 1.234 (benzene-ethene) and 2.550 (benzene-benzene PD), ∆H0 = 0.949 (benzene-ethene) and 2.310 (benzene-benzene PD), and ∆H298 = 0.130 (benzene-ethene) and 1.461 (benzene-benzene PD) kcal mol-1. Important conclusions are reached regarding the basis set convergence of the SCF, CCSD, (T), and post-CCSD(T) components. Explicitly correlated calculations are used as a sanity check on the conventional binding energies. Overall, post-CCSD(T) contributions are destabilizing by 0.028 (benzene-ethene) and 0.058 (benzene-benzene) kcal mol-1, and thus, they cannot be neglected if sub-chemical accuracy is sought (i.e., errors below 0.1 kcal mol-1). CCSD(T)/aug-cc-pwCVTZ core-valence corrections increase the binding energies by 0.018 (benzene-ethene) and 0.027 (benzene-benzene PD) kcal mol-1. Scalar-relativistic and diagonal Born-Oppenheimer corrections are negligibly small. We use our best CCSDT(Q) binding energies to evaluate the performance of MP2-based, CCSD-based, and lower-cost composite ab initio procedures for obtaining these challenging π-π stacking binding energies.
Collapse
Affiliation(s)
- Amir Karton
- School of Molecular Sciences, The University of Western Australia, Perth, WA 6009, Australia
| | - Jan M L Martin
- Department of Organic Chemistry, Weizmann Institute of Science, 76100 Reḥovot, Israel
| |
Collapse
|
24
|
Véril M, Scemama A, Caffarel M, Lipparini F, Boggio‐Pasqua M, Jacquemin D, Loos P. QUESTDB
: A database of highly accurate excitation energies for the electronic structure community. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2021. [DOI: 10.1002/wcms.1517] [Citation(s) in RCA: 31] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
- Mickaël Véril
- Laboratoire de Chimie et Physique Quantiques Université de Toulouse, CNRS, UPS Toulouse France
| | - Anthony Scemama
- Laboratoire de Chimie et Physique Quantiques Université de Toulouse, CNRS, UPS Toulouse France
| | - Michel Caffarel
- Laboratoire de Chimie et Physique Quantiques Université de Toulouse, CNRS, UPS Toulouse France
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale University of Pisa Pisa Italy
| | - Martial Boggio‐Pasqua
- Laboratoire de Chimie et Physique Quantiques Université de Toulouse, CNRS, UPS Toulouse France
| | | | - Pierre‐François Loos
- Laboratoire de Chimie et Physique Quantiques Université de Toulouse, CNRS, UPS Toulouse France
| |
Collapse
|
25
|
Donchev AG, Taube AG, Decolvenaere E, Hargus C, McGibbon RT, Law KH, Gregersen BA, Li JL, Palmo K, Siva K, Bergdorf M, Klepeis JL, Shaw DE. Quantum chemical benchmark databases of gold-standard dimer interaction energies. Sci Data 2021; 8:55. [PMID: 33568655 PMCID: PMC7876112 DOI: 10.1038/s41597-021-00833-x] [Citation(s) in RCA: 26] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2020] [Accepted: 12/14/2020] [Indexed: 12/11/2022] Open
Abstract
Advances in computational chemistry create an ongoing need for larger and higher-quality datasets that characterize noncovalent molecular interactions. We present three benchmark collections of quantum mechanical data, covering approximately 3,700 distinct types of interacting molecule pairs. The first collection, which we refer to as DES370K, contains interaction energies for more than 370,000 dimer geometries. These were computed using the coupled-cluster method with single, double, and perturbative triple excitations [CCSD(T)], which is widely regarded as the gold-standard method in electronic structure theory. Our second benchmark collection, a core representative subset of DES370K called DES15K, is intended for more computationally demanding applications of the data. Finally, DES5M, our third collection, comprises interaction energies for nearly 5,000,000 dimer geometries; these were calculated using SNS-MP2, a machine learning approach that provides results with accuracy comparable to that of our coupled-cluster training data. These datasets may prove useful in the development of density functionals, empirically corrected wavefunction-based approaches, semi-empirical methods, force fields, and models trained using machine learning methods.
Collapse
Affiliation(s)
| | | | | | - Cory Hargus
- D. E. Shaw Research, New York, NY, 10036, USA
| | | | - Ka-Hei Law
- D. E. Shaw Research, New York, NY, 10036, USA
| | | | - Je-Luen Li
- D. E. Shaw Research, New York, NY, 10036, USA
| | - Kim Palmo
- D. E. Shaw Research, New York, NY, 10036, USA
| | | | | | | | - David E Shaw
- D. E. Shaw Research, New York, NY, 10036, USA. .,Department of Biochemistry and Molecular Biophysics, Columbia University, New York, NY, 10032, USA.
| |
Collapse
|
26
|
Modrzejewski M, Yourdkhani S, Śmiga S, Klimeš J. Random-Phase Approximation in Many-Body Noncovalent Systems: Methane in a Dodecahedral Water Cage. J Chem Theory Comput 2021; 17:804-817. [PMID: 33445879 DOI: 10.1021/acs.jctc.0c00966] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
The many-body expansion (MBE) of energies of molecular clusters or solids offers a way to detect and analyze errors of theoretical methods that could go unnoticed if only the total energy of the system was considered. In this regard, the interaction between the methane molecule and its enclosing dodecahedral water cage, CH4···(H2O)20, is a stringent test for approximate methods, including density functional theory (DFT) approximations. Hybrid and semilocal DFT approximations behave erratically for this system, with three- and four-body nonadditive terms having neither the correct sign nor magnitude. Here, we analyze to what extent these qualitative errors in different MBE contributions are conveyed to post-Kohn-Sham random-phase approximation (RPA), which uses approximate Kohn-Sham orbitals as its input. The results reveal a correlation between the quality of the DFT input states and the RPA results. Moreover, the renormalized singles energy (RSE) corrections play a crucial role in all orders of the many-body expansion. For dimers, RSE corrects the RPA underbinding for every tested Kohn-Sham model: generalized-gradient approximation (GGA), meta-GGA, (meta-)GGA hybrids, as well as the optimized effective potential at the correlated level. Remarkably, the inclusion of singles in RPA can also correct the wrong signs of three- and four-body nonadditive energies as well as mitigate the excessive higher-order contributions to the many-body expansion. The RPA errors are dominated by the contributions of compact clusters. As a workable method for large systems, we propose to replace those compact contributions with CCSD(T) energies and to sum up the remaining many-body contributions up to infinity with supermolecular or periodic RPA. As a demonstration of this approach, we show that for RPA(PBE0)+RSE it suffices to apply CCSD(T) to dimers and 30 compact, hydrogen-bonded trimers to get the methane-water cage interaction energy to within 1.6% of the reference value.
Collapse
Affiliation(s)
- Marcin Modrzejewski
- Faculty of Chemistry, University of Warsaw, 02-093 Warsaw, Pasteura 1, Poland.,Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, CZ-12116 Prague 2, Czech Republic
| | - Sirous Yourdkhani
- Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, CZ-12116 Prague 2, Czech Republic
| | - Szymon Śmiga
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University, Grudziądzka 5, 87-100 Toruń, Poland
| | - Jiří Klimeš
- Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, CZ-12116 Prague 2, Czech Republic
| |
Collapse
|
27
|
Akisawa K, Hatada R, Okuwaki K, Mochizuki Y, Fukuzawa K, Komeiji Y, Tanaka S. Interaction analyses of SARS-CoV-2 spike protein based on fragment molecular orbital calculations. RSC Adv 2021; 11:3272-3279. [PMID: 35424290 PMCID: PMC8694004 DOI: 10.1039/d0ra09555a] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2020] [Accepted: 01/06/2021] [Indexed: 12/13/2022] Open
Abstract
At the stage of SARS-CoV-2 infection in human cells, the spike protein consisting of three chains, A, B, and C, with a total of 3300 residues plays a key role, and thus its structural properties and the binding nature of receptor proteins to host human cells or neutralizing antibodies has attracted considerable interest. Here, we report on interaction analyses of the spike protein in both closed (PDB-ID: 6VXX) and open (6VYB) structures, based on large-scale fragment molecular orbital (FMO) calculations at the level of up to the fourth-order Møller–Plesset perturbation with singles, doubles, and quadruples (MP4(SDQ)). Inter-chain interaction energies were evaluated for both structures, and a mutual comparison indicated considerable losses of stabilization energies in the open structure, especially in the receptor binding domain (RBD) of chain-B. The role of charged residues in inter-chain interactions was illuminated as well. By two separate calculations for the RBD complexes with angiotensin-converting enzyme 2 (ACE2) (6M0J) and B38 Fab antibody (7BZ5), it was found that the binding with ACE2 or antibody partially compensated for this stabilization loss of RBD. Visualized IFIE results seen from chain-B of spike protein.![]()
Collapse
Affiliation(s)
- Kazuki Akisawa
- Department of Chemistry and Research Center for Smart Molecules
- Faculty of Science
- Rikkyo University
- Toshima-ku
- Japan
| | - Ryo Hatada
- Department of Chemistry and Research Center for Smart Molecules
- Faculty of Science
- Rikkyo University
- Toshima-ku
- Japan
| | - Koji Okuwaki
- Department of Chemistry and Research Center for Smart Molecules
- Faculty of Science
- Rikkyo University
- Toshima-ku
- Japan
| | - Yuji Mochizuki
- Department of Chemistry and Research Center for Smart Molecules
- Faculty of Science
- Rikkyo University
- Toshima-ku
- Japan
| | - Kaori Fukuzawa
- Institute of Industrial Science
- The University of Tokyo
- Meguro-ku
- Japan
- School of Pharmacy and Pharmaceutical Sciences
| | - Yuto Komeiji
- Health and Medical Research Institute
- AIST
- Tsukuba
- Japan
| | - Shigenori Tanaka
- Graduate School of System Informatics
- Department of Computational Science
- Kobe University
- Kobe 657-8501
- Japan
| |
Collapse
|
28
|
Rettig A, Hait D, Bertels LW, Head-Gordon M. Third-Order Møller-Plesset Theory Made More Useful? The Role of Density Functional Theory Orbitals. J Chem Theory Comput 2020; 16:7473-7489. [PMID: 33161713 DOI: 10.1021/acs.jctc.0c00986] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The practical utility of Møller-Plesset (MP) perturbation theory is severely constrained by the use of Hartree-Fock (HF) orbitals. It has recently been shown that the use of regularized orbital-optimized MP2 orbitals and scaling of MP3 energy could lead to a significant reduction in MP3 error [Bertels, L. W.; J. Phys. Chem. Lett. 2019, 10, 4170 4176]. In this work, we examine whether density functional theory (DFT)-optimized orbitals can be similarly employed to improve the performance of MP theory at both the MP2 and MP3 levels. We find that the use of DFT orbitals leads to significantly improved performance for prediction of thermochemistry, barrier heights, noncovalent interactions, and dipole moments relative to the standard HF-based MP theory. Indeed, MP3 (with or without scaling) with DFT orbitals is found to surpass the accuracy of coupled-cluster singles and doubles (CCSD) for several data sets. We also found that the results are not particularly functional sensitive in most cases (although range-separated hybrid functionals with low delocalization error perform the best). MP3 based on DFT orbitals thus appears to be an efficient, noniterative O(N6) scaling wave-function approach for single-reference electronic structure computations. Scaled MP2 with DFT orbitals is also found to be quite accurate in many cases, although modern double hybrid functionals are likely to be considerably more accurate.
Collapse
Affiliation(s)
- Adam Rettig
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States.,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Diptarka Hait
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States.,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Luke W Bertels
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States.,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
29
|
Ladóczki B, Uejima M, Ten-no SL. Third-order Epstein–Nesbet perturbative correction to the initiator approximation of configuration space quantum Monte Carlo. J Chem Phys 2020; 153:114112. [DOI: 10.1063/5.0022101] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Bence Ladóczki
- Graduate School of System Informatics, Kobe University, Nada-ku, Kobe 657-8501, Japan
| | - Motoyuki Uejima
- Graduate School of Science, Technology, and Innovation, Kobe University, Nada-ku, Kobe 657-8501, Japan
| | - Seiichiro L. Ten-no
- Graduate School of System Informatics, Kobe University, Nada-ku, Kobe 657-8501, Japan
- Graduate School of Science, Technology, and Innovation, Kobe University, Nada-ku, Kobe 657-8501, Japan
| |
Collapse
|
30
|
Řezáč J. Non-Covalent Interactions Atlas Benchmark Data Sets 2: Hydrogen Bonding in an Extended Chemical Space. J Chem Theory Comput 2020; 16:6305-6316. [DOI: 10.1021/acs.jctc.0c00715] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jan Řezáč
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, 166 10 Prague, Czech Republic
| |
Collapse
|
31
|
Suliman S, Pitoňák M, Cernusak I, Louis F. On the applicability of the MP2.5 approximation for open-shell systems. Case study of atmospheric reactivity. COMPUT THEOR CHEM 2020. [DOI: 10.1016/j.comptc.2020.112901] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
32
|
Zhang D, Truhlar DG. Spin Splitting Energy of Transition Metals: A New, More Affordable Wave Function Benchmark Method and Its Use to Test Density Functional Theory. J Chem Theory Comput 2020; 16:4416-4428. [PMID: 32525690 DOI: 10.1021/acs.jctc.0c00518] [Citation(s) in RCA: 33] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Accurately predicting the spin splitting energy of chemical species is important for understanding their reactivity and magnetic properties, but it is very challenging, especially for molecules containing transition metals. One impediment to progress is the scarcity of accurate benchmark data. Here we report a set of calculations designed to yield reliable benchmarks for simple transition-metal complexes that can be used to test density functional methods that are affordable for large systems of more practical interest. Various wave function methods are tested against experiment for Fe2+, Fe3+, and Co3+, including CASSCF, CASPT2, CASPT3, MRCISD, MRCISD+Q, ACPF, AQCC, CCSD(T), and CASPT2/CCSD(T) and also a new method called CASPT2.5, which is performed by taking the average of the CASPT2 and CASPT3 energies. We find that MRCISD+Q, ACPF, and AQCC require smaller active spaces for good accuracy than are required by CASPT2 and CASPT3, and this aspect may be important for calculations on larger molecules; here we find that CASPT2.5 extrapolated to a complete basis set is the most suitable method-in terms of computational cost and in terms of accuracy on monatomic systems-and therefore we chose this method for molecular benchmarks. Then Kohn-Sham density functional calculations with 60 exchange-correlation functionals are tested for FeF2, FeCl2, and CoF2. We find that MN15-L, M06-SX, and revM06 have very good agreement with CASPT2.5 benchmarks in terms of both the spin splitting energy and the optimized geometry for each spin state. In addition, we recommend def2-TZVP as the most suitable basis set to perform density functional calculations for molecular spin splitting energies; extra polarization functions in the basis set do not help to increase the accuracy of the spin splitting energy in KS calculations.
Collapse
Affiliation(s)
- Dayou Zhang
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
33
|
Smith DGA, Burns LA, Simmonett AC, Parrish RM, Schieber MC, Galvelis R, Kraus P, Kruse H, Di Remigio R, Alenaizan A, James AM, Lehtola S, Misiewicz JP, Scheurer M, Shaw RA, Schriber JB, Xie Y, Glick ZL, Sirianni DA, O’Brien JS, Waldrop JM, Kumar A, Hohenstein EG, Pritchard BP, Brooks BR, Schaefer HF, Sokolov AY, Patkowski K, DePrince AE, Bozkaya U, King RA, Evangelista FA, Turney JM, Crawford TD, Sherrill CD. Psi4 1.4: Open-source software for high-throughput quantum chemistry. J Chem Phys 2020; 152:184108. [PMID: 32414239 PMCID: PMC7228781 DOI: 10.1063/5.0006002] [Citation(s) in RCA: 337] [Impact Index Per Article: 84.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2020] [Accepted: 04/12/2020] [Indexed: 12/13/2022] Open
Abstract
PSI4 is a free and open-source ab initio electronic structure program providing implementations of Hartree-Fock, density functional theory, many-body perturbation theory, configuration interaction, density cumulant theory, symmetry-adapted perturbation theory, and coupled-cluster theory. Most of the methods are quite efficient, thanks to density fitting and multi-core parallelism. The program is a hybrid of C++ and Python, and calculations may be run with very simple text files or using the Python API, facilitating post-processing and complex workflows; method developers also have access to most of PSI4's core functionalities via Python. Job specification may be passed using The Molecular Sciences Software Institute (MolSSI) QCSCHEMA data format, facilitating interoperability. A rewrite of our top-level computation driver, and concomitant adoption of the MolSSI QCARCHIVE INFRASTRUCTURE project, makes the latest version of PSI4 well suited to distributed computation of large numbers of independent tasks. The project has fostered the development of independent software components that may be reused in other quantum chemistry programs.
Collapse
Affiliation(s)
| | - Lori A. Burns
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| | - Andrew C. Simmonett
- National Institutes of Health – National Heart,
Lung and Blood Institute, Laboratory of Computational Biology, Bethesda,
Maryland 20892, USA
| | - Robert M. Parrish
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| | - Matthew C. Schieber
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| | | | - Peter Kraus
- School of Molecular and Life Sciences, Curtin
University, Kent St., Bentley, Perth, Western Australia 6102,
Australia
| | - Holger Kruse
- Institute of Biophysics of the Czech Academy of
Sciences, Královopolská 135, 612 65 Brno, Czech
Republic
| | - Roberto Di Remigio
- Department of Chemistry, Centre for Theoretical
and Computational Chemistry, UiT, The Arctic University of Norway, N-9037
Tromsø, Norway
| | - Asem Alenaizan
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| | - Andrew M. James
- Department of Chemistry, Virginia
Tech, Blacksburg, Virginia 24061, USA
| | - Susi Lehtola
- Department of Chemistry, University of
Helsinki, P.O. Box 55 (A. I. Virtasen aukio 1), FI-00014 Helsinki,
Finland
| | - Jonathon P. Misiewicz
- Center for Computational Quantum Chemistry,
University of Georgia, Athens, Georgia 30602, USA
| | - Maximilian Scheurer
- Interdisciplinary Center for Scientific
Computing, Heidelberg University, D-69120 Heidelberg,
Germany
| | - Robert A. Shaw
- ARC Centre of Excellence in Exciton Science,
School of Science, RMIT University, Melbourne, VIC 3000,
Australia
| | - Jeffrey B. Schriber
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| | - Yi Xie
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| | - Zachary L. Glick
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| | - Dominic A. Sirianni
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| | - Joseph Senan O’Brien
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| | - Jonathan M. Waldrop
- Department of Chemistry and Biochemistry, Auburn
University, Auburn, Alabama 36849, USA
| | - Ashutosh Kumar
- Department of Chemistry, Virginia
Tech, Blacksburg, Virginia 24061, USA
| | - Edward G. Hohenstein
- SLAC National Accelerator Laboratory, Stanford
PULSE Institute, Menlo Park, California 94025,
USA
| | | | - Bernard R. Brooks
- National Institutes of Health – National Heart,
Lung and Blood Institute, Laboratory of Computational Biology, Bethesda,
Maryland 20892, USA
| | - Henry F. Schaefer
- Center for Computational Quantum Chemistry,
University of Georgia, Athens, Georgia 30602, USA
| | - Alexander Yu. Sokolov
- Department of Chemistry and Biochemistry, The
Ohio State University, Columbus, Ohio 43210, USA
| | - Konrad Patkowski
- Department of Chemistry and Biochemistry, Auburn
University, Auburn, Alabama 36849, USA
| | - A. Eugene DePrince
- Department of Chemistry and Biochemistry,
Florida State University, Tallahassee, Florida 32306-4390,
USA
| | - Uğur Bozkaya
- Department of Chemistry, Hacettepe
University, Ankara 06800, Turkey
| | - Rollin A. King
- Department of Chemistry, Bethel
University, St. Paul, Minnesota 55112, USA
| | | | - Justin M. Turney
- Center for Computational Quantum Chemistry,
University of Georgia, Athens, Georgia 30602, USA
| | | | - C. David Sherrill
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| |
Collapse
|
34
|
Yu F, Wang Y. Dual‐hybrid direct random phase approximation and second‐order screened exchange with nonlocal van der Waals correlations for noncovalent interactions. J Comput Chem 2020; 41:1018-1025. [DOI: 10.1002/jcc.26149] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2019] [Revised: 01/05/2020] [Accepted: 01/06/2020] [Indexed: 11/09/2022]
Affiliation(s)
- Feng Yu
- Department of Physics, School of ScienceXi'an Technological University Xi'an Shaanxi China
| | - Yaoting Wang
- Department of Physics, School of ScienceXi'an Technological University Xi'an Shaanxi China
| |
Collapse
|
35
|
Warden CE, Smith DGA, Burns LA, Bozkaya U, Sherrill CD. Efficient and automated computation of accurate molecular geometries using focal-point approximations to large-basis coupled-cluster theory. J Chem Phys 2020; 152:124109. [DOI: 10.1063/5.0004863] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Constance E. Warden
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Daniel G. A. Smith
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Lori A. Burns
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Uğur Bozkaya
- Department of Chemistry, Hacettepe University, Ankara 06800, Turkey
| | - C. David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| |
Collapse
|
36
|
Nguyen BD, Chen GP, Agee MM, Burow AM, Tang MP, Furche F. Divergence of Many-Body Perturbation Theory for Noncovalent Interactions of Large Molecules. J Chem Theory Comput 2020; 16:2258-2273. [DOI: 10.1021/acs.jctc.9b01176] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Affiliation(s)
- Brian D. Nguyen
- University of California, Irvine, Department of Chemistry, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Guo P. Chen
- University of California, Irvine, Department of Chemistry, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Matthew M. Agee
- University of California, Irvine, Department of Chemistry, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Asbjörn M. Burow
- University of California, Irvine, Department of Chemistry, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Matthew P. Tang
- University of California, Irvine, Department of Chemistry, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Filipp Furche
- University of California, Irvine, Department of Chemistry, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| |
Collapse
|
37
|
Chen J, Chan B, Shao Y, Ho J. How accurate are approximate quantum chemical methods at modelling solute-solvent interactions in solvated clusters? Phys Chem Chem Phys 2020; 22:3855-3866. [PMID: 32022044 PMCID: PMC7394230 DOI: 10.1039/c9cp06792b] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
In this paper, the performance of a wide range of DFT methods is assessed for the calculation of interaction energies of thermal clusters of a solute in water. Three different charge states (neutral, proton transfer transition state and zwitterion) of glycine were solvated by 1 to 40 water molecules as sampled from molecular dynamics simulations. While some ab initio composite methods that employ insufficiently large basis sets incurred significant errors even for a cluster containing only 5 water molecules relative to the W1X-2 benchmark, the DLPNO-CCSD(T)/CBS and DSD-PBEP86 (triple zeta basis set) levels of theory predicted very accurate interaction energies. These levels of theory were used to benchmark the performance of 16 density functionals from different rungs of Jacob's Ladder. Of the Rung 4 functionals examined, the ωB97M-V and ωB97X-V functionals stood out for predicting absolute interaction energies in 40-water clusters with mean absolute deviations (MAD) ∼4 kJ mol-1. The B3LYP-D3(BJ) functional performed exceptionally well with a MAD ∼1.7 kJ mol-1 and is the overall best performing method. Calculations of relative interaction energies allow for cancellation of systematic errors, including basis set truncation and superposition errors, and the ωB97M-V and B3LYP-D3(BJ) double zeta basis set calculations yielded relative interaction energies that are within ∼3 kJ mol-1 of the benchmark. The ONIOM approximation provides another strategy for accelerating the calculation of accurate absolute interaction energies provided that the calculations have converged with respect to the size of the "high-level-layer".
Collapse
Affiliation(s)
- Junbo Chen
- School of Chemistry, University of New South Wales, Sydney, NSW 2052, Australia.
| | - Bun Chan
- Graduate School of Engineering, Nagasaki University, Bunkyo-Machi 1-14, Nagasaki 852-8521, Japan.
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, USA
| | - Junming Ho
- School of Chemistry, University of New South Wales, Sydney, NSW 2052, Australia.
| |
Collapse
|
38
|
Sasanuma Y, Tanaka S. Molecular Design of Aromatic Polythionoesters. ACS OMEGA 2020; 5:3016-3029. [PMID: 32095725 PMCID: PMC7033963 DOI: 10.1021/acsomega.9b04111] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/02/2019] [Accepted: 01/10/2020] [Indexed: 06/10/2023]
Abstract
As an example of molecular design of new polymers, structures and properties of poly(ethylene thionoterephthalate) (PET[S2]) and the related polymers have been predicted from calculations of ab initio molecular orbital (MO) theory, rotational isomeric state (RIS) scheme, and periodic density functional theory (DFT). The MO calculations were confirmed by NMR experiments and introduced to the RIS scheme for PET[S2] to yield its configurational properties, which are compared herein with those of analogous polyester, polythioester, and polydithioester. Configurational properties of randomly thiono-substituted poly(ethylene terephthalate) (PET), PET[S z O1-z ], were also evaluated as a function of sulfidity (z). On the assumption that the crystal of PET[S2] can be expressed as an isomorphic replacement of the PET crystal, the crystal structure was optimized by a periodic DFT simulation and its Young's moduli in the a-, b-, and c-axis directions were, respectively, evaluated to be Ea = 0.94(7.20) GPa, E b = 19.58(22.26) GPa, and E c = 142.1(182.4) GPa, where the parenthetic values are those of the PET crystal. There is a possibility that properties of PET[S z O1-z ] will be controlled between those of PET and PET[S2] by adjusting the sulfidity. The potential practical applications of the polythionoesters are also discussed herein. By purely theoretical computations, the structures and properties of the not-yet synthesized polymers were predicted quantitatively; that is, the theoretical molecular design of new polymers has been achieved.
Collapse
|
39
|
Loos PF, Jacquemin D. Is ADC(3) as Accurate as CC3 for Valence and Rydberg Transition Energies? J Phys Chem Lett 2020; 11:974-980. [PMID: 31913639 DOI: 10.1021/acs.jpclett.9b03652] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
The search for new models rapidly delivering accurate excited-state energies and properties is one of the most active research lines of theoretical chemistry. Along with these developments, the performance of known methods is constantly reassessed on the basis of new benchmark values. In this Letter, we show that the third-order algebraic diagrammatic construction, ADC(3), does not yield transition energies of the same quality as the third-order coupled cluster method, CC3. This is demonstrated by extensive comparisons with several hundred high-quality vertical transition energies obtained with FCI, CCSDTQ, and CCSDT. Direct comparisons with experimental 0-0 energies of small- and medium-size molecules support the same conclusion, which holds for both valence and Rydberg transitions. Considering these results, we introduce a composite approach, ADC(2.5), which consists of averaging the ADC(2) and ADC(3) excitation energies. Although ADC(2.5) does not match the CC3 accuracy, it significantly improves the ADC(3) results, especially for vertical energies.
Collapse
Affiliation(s)
- Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), CNRS , Université de Toulouse , 31077 Toulouse , France
| | - Denis Jacquemin
- Laboratoire CEISAM UMR UN-CNRS 6230 , Université de Nantes , F-44000 Nantes , France
| |
Collapse
|
40
|
Theoretical Description of R-X⋯NH 3 Halogen Bond Complexes: Effect of the R Group on the Complex Stability and Sigma-Hole Electron Depletion. Molecules 2020; 25:molecules25030530. [PMID: 31991810 PMCID: PMC7037998 DOI: 10.3390/molecules25030530] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2019] [Revised: 01/14/2020] [Accepted: 01/17/2020] [Indexed: 11/16/2022] Open
Abstract
In the present work, a number of R–X⋯NH3 (X = Cl, Br, and I) halogen bonded systems were theoretical studied by means of DFT calculations performed at the ωB97XD/6-31+G(d,p) level of theory in order to get insights on the effect of the electron-donating or electron-withdrawing character of the different R substituent groups (R = halogen, methyl, partially fluorinated methyl, perfluoro-methyl, ethyl, vinyl, and acetyl) on the stability of the halogen bond. The results indicate that the relative stability of the halogen bond follows the Cl < Br < I trend considering the same R substituent whereas the more electron-withdrawing character of the R substituent the more stable the halogen bond. Refinement of the latter results, performed at the MP2/6-31+G(d,p) level showed that the DFT and the MP2 binding energies correlate remarkably well, suggesting that the Grimme’s type dispersion-corrected functional produces reasonable structural and energetic features of halogen bond systems. DFT results were also observed to agree with more refined calculations performed at the CCSD(T) level. In a further stage, a more thorough analysis of the R–Br⋯NH3 complexes was performed by means of a novel electron localization/delocalization tool, defined in terms of an Information Theory, IT, based quantity obtained from the conditional pair density. For the latter, our in-house developed C++/CUDA program, called KLD (acronym of Kullback–Leibler divergence), was employed. KLD results mapped onto the one-electron density plotted at a 0.04 a.u. isovalue, showed that (i) as expected, the localized electron depletion of the Br sigma-hole is largely affected by the electron-withdrawing character of the R substituent group and (ii) the R–X bond is significantly polarized due to the presence of the NH3 molecule in the complexes. The afore-mentioned constitutes a clear indication of the dominant character of electrostatics on the stabilization of halogen bonds in agreement with a number of studies reported in the main literature. Finally, the cooperative effects on the [Br—CN]n system (n = 1–8) was evaluated at the MP2/6-31+G(d,p) level, where it was observed that an increase of about ~14.2% on the complex stability is obtained when going from n = 2 to n = 8. The latter results were corroborated by the analysis of the changes on the Fermi-hole localization pattern on the halogen bond zones, which suggests an also important contribution of the electron correlation in the stabilization of these systems.
Collapse
|
41
|
Cabaleiro-Lago EM, Rodríguez-Otero J. Curvature and size effects hinder halogen bonds with extended π systems. Phys Chem Chem Phys 2020; 22:21988-22002. [DOI: 10.1039/d0cp03466e] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The curvature of aromatic systems strengthens the interaction by the concave face while it weakens by the convex one. Parallel structures are favoured over halogen bonded ones.
Collapse
Affiliation(s)
- Enrique M. Cabaleiro-Lago
- Departamento de Química Física
- Facultade de Ciencias
- Universidade de Santiago de Compostela
- Campus de Lugo
- Av. Alfonso X El Sabio
| | - Jesús Rodríguez-Otero
- Departamento de Química Física
- Facultade de Química
- Universidade de Santiago de Compostela
- 15782 Santiago de Compostela
- Spain
| |
Collapse
|
42
|
MOCHIZUKI Y, SAKAKURA K, WATANABE H, OKUWAKI K, KATO K, WATANABE N, OKIYAMA Y, FUKUZAWA K, NAKANO T. Development Status of ABINIT-MP in 2020. JOURNAL OF COMPUTER CHEMISTRY-JAPAN 2020. [DOI: 10.2477/jccj.2021-0015] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Affiliation(s)
- Yuji MOCHIZUKI
- Department of Chemistry, Rikkyo University (Nishi-Ikebukuro 3-34-1, Toshima-ku, Tokyo, 171-8501, Japan)
- Institute of Industrial Science, The University of Tokyo (Komaba 4-6-1, Meguro-ku, Tokyo, 153-8505, Japan)
| | - Kota SAKAKURA
- Foundation for Computational Science (Minatojima Minamicho 7-1-28 Computational Science Center Building 1F, Chuo-ku, Kobe, 650-0047, Japan)
| | - Hiromasa WATANABE
- HPC Systems Inc. (Kaigan 3-9-15 LOOP-X 8F, Minato-ku, Tokyo, 108-0022, Japan)
| | - Koji OKUWAKI
- Department of Chemistry, Rikkyo University (Nishi-Ikebukuro 3-34-1, Toshima-ku, Tokyo, 171-8501, Japan)
| | - Koichiro KATO
- Department of Applied Chemistry, Graduate School of Engineering, Kyushu University (744 Moto-oka, Nishi-ku, Fukuoka, 819-0395, Japan)
| | - Naoki WATANABE
- Mizuho Information & Research Institute, Inc. (Kanda-Nishikicho 2-3, Chiyoda-ku, Tokyo, 101-8443, Japan)
| | - Yoshio OKIYAMA
- National Institute of Health Sciences (Tonomachi 3-25-26, Kawasaki-ku, Kawasaki-shi Kanagawa, 210-951, Japan)
| | - Kaori FUKUZAWA
- Institute of Industrial Science, The University of Tokyo (Komaba 4-6-1, Meguro-ku, Tokyo, 153-8505, Japan)
- School of Pharmacy and Pharmaceutical Sciences, Hoshi University (Ebara 2-4-41, Shinagawa-ku, Tokyo, 142-8501, Japan)
| | - Tatsuya NAKANO
- National Institute of Health Sciences (Tonomachi 3-25-26, Kawasaki-ku, Kawasaki-shi Kanagawa, 210-951, Japan)
| |
Collapse
|
43
|
Bakowies D. Estimating Systematic Error and Uncertainty in Ab Initio Thermochemistry: II. ATOMIC(hc) Enthalpies of Formation for a Large Set of Hydrocarbons. J Chem Theory Comput 2019; 16:399-426. [DOI: 10.1021/acs.jctc.9b00974] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Dirk Bakowies
- Institute of Physical Chemistry, Department of Chemistry, University of Basel, Klingelbergstr. 80, CH 4056 Basel, Switzerland
| |
Collapse
|
44
|
Heßelmann A. Intermolecular interaction energies from fourth order many-body perturbation theory. Impact of individual electron correlation contributions. J Chem Phys 2019; 151:114105. [DOI: 10.1063/1.5112178] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Andreas Heßelmann
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91058 Erlangen, Germany
| |
Collapse
|
45
|
Bertels LW, Lee J, Head-Gordon M. Third-Order Møller-Plesset Perturbation Theory Made Useful? Choice of Orbitals and Scaling Greatly Improves Accuracy for Thermochemistry, Kinetics, and Intermolecular Interactions. J Phys Chem Lett 2019; 10:4170-4176. [PMID: 31259560 DOI: 10.1021/acs.jpclett.9b01641] [Citation(s) in RCA: 29] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
We develop and test methods that include second- and third-order perturbation theory (MP3) using orbitals obtained from regularized orbital-optimized second-order perturbation theory, κ-OOMP2, denoted as MP3:κ-OOMP2. Testing MP3:κ-OOMP2 shows RMS errors that are 1.7-5 times smaller than those of MP3 across 7 data sets. To do still better, empirical training of the scaling factors for the second- and third-order correlation energies and the regularization parameter on one of those data sets led to an unregularized scaled (c2 = 1.0; c3 = 0.8) denoted as MP2.8:κ-OOMP2. MP2.8:κ-OOMP2 yields significant additional improvement over MP3:κ-OOMP2 in 4 of 6 test data sets on thermochemistry, kinetics, and noncovalent interactions. Remarkably, these two methods outperform coupled cluster with singles and doubles in 5 of the 7 data sets considered, at greatly reduced cost (no O(N6) iterations).
Collapse
Affiliation(s)
- Luke W Bertels
- Department of Chemistry , University of California , Berkeley , California 94720 , United States
| | - Joonho Lee
- Department of Chemistry , University of California , Berkeley , California 94720 , United States
| | - Martin Head-Gordon
- Department of Chemistry , University of California , Berkeley , California 94720 , United States
- Chemical Sciences Division , Lawrence Berkeley National Laboratory , Berkeley , California 94720 , United States
| |
Collapse
|
46
|
Tian T, Sun X, Weiske T, Cai Y, Geng C, Li J, Schwarz H. Reassessment of the Mechanisms of Thermal C-H Bond Activation of Methane by Cationic Magnesium Oxides: A Critical Evaluation of the Suitability of Different Density Functionals. Chemphyschem 2019; 20:1812-1821. [PMID: 31120181 DOI: 10.1002/cphc.201900508] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2019] [Revised: 05/22/2019] [Indexed: 12/24/2022]
Abstract
The mechanisms of the thermal reactions of the two iconic magnesium oxide cations MgO.+ and Mg2 O2 .+ with methane have been re-evaluated at the CCSD(T)/CBS//CCSD/def2-TZVP level of theory. For the reaction of MgO.+ with CH4 , only the classical hydrogen-atom transfer (HAT) was found; in contrast, for the Mg2 O2 .+ /CH4 couple, both HAT and proton-coupled electron-transfer (PCET) exist as mechanistic variants. In order to evaluate the suitability of density functional theory (DFT) methods, the reactions were computed by using 27 density functionals. The results obtained demonstrate that the various DFT methods often deliver rather different results for both geometric and energetic features. As to the prediction of the apparent barriers, pure functionals give the largest mean absolute errors. BMK, ωB97XD, and the double-hybrid functional mPW2PLYP were confirmed to come closest to the results provided by CCSD(T)/CBS. Thus, mechanistic conclusions based on a single DFT method should be viewed with great caution. In summary, this study may assist in the selection of a suitable quantum chemical method to unravel the mechanistic details of C-H bond activation by charged metal oxides.
Collapse
Affiliation(s)
- Tian Tian
- Institute of Theoretical Chemistry, Jilin University, Changchun, 130023, People's Republic of China
| | - Xiaoli Sun
- Institute of Theoretical Chemistry, Jilin University, Changchun, 130023, People's Republic of China
| | - Thomas Weiske
- Institut für Chemie, Technische Universität Berlin, Straße des 17. Juni 115, 10623, Berlin, Germany
| | - Yuxi Cai
- Institute of Theoretical Chemistry, Jilin University, Changchun, 130023, People's Republic of China
| | - Caiyun Geng
- Institut für Chemie, Technische Universität Berlin, Straße des 17. Juni 115, 10623, Berlin, Germany
| | - Jilai Li
- Institute of Theoretical Chemistry, Jilin University, Changchun, 130023, People's Republic of China.,Institut für Chemie, Technische Universität Berlin, Straße des 17. Juni 115, 10623, Berlin, Germany
| | - Helmut Schwarz
- Institut für Chemie, Technische Universität Berlin, Straße des 17. Juni 115, 10623, Berlin, Germany
| |
Collapse
|
47
|
Cacelli I, Lipparini F, Greff da Silveira L, Jacobs M, Livotto PR, Prampolini G. Accurate interaction energies by spin component scaled Möller-Plesset second order perturbation theory calculations with optimized basis sets (SCS-MP2mod): Development and application to aromatic heterocycles. J Chem Phys 2019; 150:234113. [DOI: 10.1063/1.5094288] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Affiliation(s)
- Ivo Cacelli
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, I-56124 Pisa, Italy
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR), Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, I-56124 Pisa, Italy
| | - Leandro Greff da Silveira
- Instituto de Química, Universidade Federal do Rio Grande do Sul, Avenida Bento Gonçalves 9500, CEP 91501-970 Porto Alegre, Brazil
| | - Matheus Jacobs
- Institut für Physik, Humboldt-Universität zu Berlin, Newtonstr. 15, 12489 Berlin, Germany
- IRIS Adelrshof, Humboldt-Universität zu Berlin, Zum Großen Windkanal 6, 12489 Berlin, Germany
| | - Paolo Roberto Livotto
- Instituto de Química, Universidade Federal do Rio Grande do Sul, Avenida Bento Gonçalves 9500, CEP 91501-970 Porto Alegre, Brazil
| | - Giacomo Prampolini
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR), Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| |
Collapse
|
48
|
Kruse H, Banáš P, Šponer J. Investigations of Stacked DNA Base-Pair Steps: Highly Accurate Stacking Interaction Energies, Energy Decomposition, and Many-Body Stacking Effects. J Chem Theory Comput 2018; 15:95-115. [DOI: 10.1021/acs.jctc.8b00643] [Citation(s) in RCA: 43] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Holger Kruse
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Pavel Banáš
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University, 17 Listopadu 12, 77146 Olomouc, Czech Republic
| | - Jiřı́ Šponer
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University, 17 Listopadu 12, 77146 Olomouc, Czech Republic
| |
Collapse
|
49
|
Abstract
Noncovalent interactions in the single crystal of 3,6-Cl2-closo-1,2-P2B10H8 and in the crystal of closo-1,7-P2B10Cl10•toluene were analyzed by means of quantum chemical computations. The crystal packing in the second crystal was dominated by numerous B-Cl···Cl-B dihalogen and strong B-P···π pnictogen bonds, the latter of which were characterized by a small length of 3.08 Å and a large interaction energy value, exceeding −10 kcal mol−1.
Collapse
|
50
|
Ishikawa T, Sakakura K, Mochizuki Y. RI-MP3 calculations of biomolecules based on the fragment molecular orbital method. J Comput Chem 2018; 39:1970-1978. [PMID: 30277590 DOI: 10.1002/jcc.25368] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2018] [Revised: 05/04/2018] [Accepted: 05/09/2018] [Indexed: 12/24/2022]
Abstract
In this study, the third-order Møller-Plesset perturbation (MP3) theory using the resolution of the identity (RI) approximation was combined with the fragment molecular orbital (FMO) method to efficiently calculate a high-order electron correlation energy of biomolecular systems. We developed a new algorithm for the RI-MP3 calculation, which can be used with the FMO scheme. After test calculations using a small molecule, the FMO-RI-MP3 calculations were performed for two biomolecular systems comprising a protein and a ligand. The computational cost of these calculations was only around 5 and 4 times higher than those of the FMO-RHF calculations. The error associated with the RI approximation was around 2.0% of the third-order correlation contribution to the total energy. However, the RI approximation error in the interaction energy between the protein and ligand molecule was insignificantly small, which reflected the negligible error in the inter fragment interaction energy. © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Takeshi Ishikawa
- Department of Molecular Microbiology and Immunology, Graduate School of Biomedical Sciences, Nagasaki University, 1-12-4 Sakamoto, Nagasaki 852-8523, Japan
| | - Kota Sakakura
- 1st Government and Public Solutions Division, NEC Corporation, 7-1, Shiba 5-chome, Minato-ku, Tokyo, 108-8001, Japan
| | - Yuji Mochizuki
- Department of Chemistry and Research Center for Smart Molecules, Faculty of Science, Rikkyo University, 3-34-1 Nishi-ikebukuro, Toshima-ku, Tokyo, 171-8501, Japan.,Institute of Industrial Science, The University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo, 153-8505, Japan
| |
Collapse
|