1
|
Heßelmann A, Giner E, Reinhardt P, Knowles PJ, Werner HJ, Toulouse J. A density-fitting implementation of the density-based basis-set correction method. J Comput Chem 2024; 45:1247-1253. [PMID: 38348951 DOI: 10.1002/jcc.27325] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2023] [Revised: 12/13/2023] [Accepted: 01/26/2024] [Indexed: 04/19/2024]
Abstract
This work reports an efficient density-fitting implementation of the density-based basis-set correction (DBBSC) method in the MOLPRO software. This method consists in correcting the energy calculated by a wave-function method with a given basis set by an adapted basis-set correction density functional incorporating the short-range electron correlation effects missing in the basis set, resulting in an accelerated convergence to the complete-basis-set limit. Different basis-set correction density-functional approximations are explored and the complementary-auxiliary-basis-set single-excitation correction is added. The method is tested on a benchmark set of reaction energies at the second-order Møller-Plesset (MP2) level and a comparison with the explicitly correlated MP2-F12 method is provided. The results show that the DBBSC method greatly accelerates the basis convergence of MP2 reaction energies, without reaching the accuracy of the MP2-F12 method but with a lower computational cost.
Collapse
Affiliation(s)
- Andreas Heßelmann
- Institute for Theoretical Chemistry, University of Stuttgart, Stuttgart, Germany
| | - Emmanuel Giner
- Laboratoire de Chimie Théorique, Sorbonne Université and CNRS, Paris, France
| | - Peter Reinhardt
- Laboratoire de Chimie Théorique, Sorbonne Université and CNRS, Paris, France
| | | | - Hans-Joachim Werner
- Institute for Theoretical Chemistry, University of Stuttgart, Stuttgart, Germany
| | - Julien Toulouse
- Laboratoire de Chimie Théorique, Sorbonne Université and CNRS, Paris, France
- Institut Universitaire de France, Paris, France
| |
Collapse
|
2
|
Graf D, Thom AJW. Corrected density functional theory and the random phase approximation: Improved accuracy at little extra cost. J Chem Phys 2023; 159:174106. [PMID: 37921249 DOI: 10.1063/5.0168569] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2023] [Accepted: 10/16/2023] [Indexed: 11/04/2023] Open
Abstract
We recently introduced an efficient methodology to perform density-corrected Hartree-Fock density functional theory [DC(HF)-DFT] calculations and an extension to it we called "corrected" HF DFT [C(HF)-DFT] [Graf and Thom, J. Chem. Theory Comput. 19 5427-5438 (2023)]. In this work, we take a further step and combine C(HF)-DFT, augmented with a straightforward orbital energy correction, with the random phase approximation (RPA). We refer to the resulting methodology as corrected HF RPA [C(HF)-RPA]. We evaluate the proposed methodology across various RPA methods: direct RPA (dRPA), RPA with an approximate exchange kernel, and RPA with second-order screened exchange. C(HF)-dRPA demonstrates very promising performance; for RPA with exchange methods, on the other hand, we often find over-corrections.
Collapse
Affiliation(s)
- Daniel Graf
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, England
| | - Alex J W Thom
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, England
| |
Collapse
|
3
|
Werner HJ, Hansen A. Accurate Calculation of Isomerization and Conformational Energies of Larger Molecules Using Explicitly Correlated Local Coupled Cluster Methods in Molpro and ORCA. J Chem Theory Comput 2023; 19:7007-7030. [PMID: 37486154 DOI: 10.1021/acs.jctc.3c00270] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/25/2023]
Abstract
An overview of the approximations in the explicitly correlated local coupled cluster methods PNO-LCCSD(T)-F12 in Molpro and DLPNO-CCSD(T)F12 in ORCA is given. Options to select the domains of projected atomic orbitals (PAOs), pair natural orbitals (PNOs), and triples natural orbitals (TNOs) in both programs are described and compared in detail. The two programs are applied to compute isomerization and conformational energies of the ISOL24 and ACONFL test sets, where the former is part of the GMTKN55 benchmark suite. Thorough studies of basis set effects are presented for selected systems. These revealed large intramolecular basis set superposition effects that make it practically impossible to reliably determine the complete basis set (CBS) limits without including explicitly correlated terms. The latter strongly reduce the basis set dependence and at the same time also errors caused by the local domain approximations. On the basis of these studies, the PNO-LCCSD(T)-F12 method is applied to determine new reference energies for the above-mentioned benchmark sets. We are confident that our results should agree within a few tenths of a kcal mol-1 with the (unknown) CCSD(T)/CBS values, which therefore allowed us to define computational settings for accurate explicitly correlated local coupled cluster methods with moderate computational effort. With these protocols, especially PNO-LCCSD(T)-F12b/AVTZ', reliable reference values for comprehensive benchmark sets can be generated efficiently. This can significantly advance the development and evaluation of the performance of approximate electronic structure methods, especially improved density functional approximations or machine learning approaches.
Collapse
Affiliation(s)
- Hans-Joachim Werner
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstrasse 4, D-53115 Bonn, Germany
| |
Collapse
|
4
|
Graf D, Thom AJW. Simple and Efficient Route toward Improved Energetics within the Framework of Density-Corrected Density Functional Theory. J Chem Theory Comput 2023; 19:5427-5438. [PMID: 37525457 PMCID: PMC10448722 DOI: 10.1021/acs.jctc.3c00441] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2023] [Indexed: 08/02/2023]
Abstract
The crucial step in density-corrected Hartree-Fock density functional theory (DC(HF)-DFT) is to decide whether the density produced by the density functional for a specific calculation is erroneous and, hence, should be replaced by, in this case, the HF density. We introduce an indicator, based on the difference in noninteracting kinetic energies between DFT and HF calculations, to determine when the HF density is the better option. Our kinetic energy indicator directly compares the self-consistent density of the analyzed functional with the HF density, is size-intensive, reliable, and most importantly highly efficient. Moreover, we present a procedure that makes best use of the computed quantities necessary for DC(HF)-DFT by additionally evaluating a related hybrid functional and, in that way, not only "corrects" the density but also the functional itself; we call that procedure corrected Hartree-Fock density functional theory (C(HF)-DFT).
Collapse
Affiliation(s)
- Daniel Graf
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, United Kingdom
| | - Alex J. W. Thom
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
5
|
Corzo HH, Hillers-Bendtsen AE, Barnes A, Zamani AY, Pawłowski F, Olsen J, Jørgensen P, Mikkelsen KV, Bykov D. Corrigendum: Coupled cluster theory on modern heterogeneous supercomputers. Front Chem 2023; 11:1256510. [PMID: 37654900 PMCID: PMC10466216 DOI: 10.3389/fchem.2023.1256510] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2023] [Accepted: 07/11/2023] [Indexed: 09/02/2023] Open
Abstract
[This corrects the article DOI: 10.3389/fchem.2023.1154526.].
Collapse
Affiliation(s)
| | | | | | - Abdulrahman Y. Zamani
- Department of Chemistry and Biochemistry and Center for Chemical Computation and Theory, University of California, Merced, CA, United States
| | - Filip Pawłowski
- Department of Chemistry and Biochemistry, Auburn University, Auburn, AL, United States
| | - Jeppe Olsen
- Department of Chemistry, Aarhus University, Aarhus, Denmark
| | - Poul Jørgensen
- Department of Chemistry, Aarhus University, Aarhus, Denmark
| | - Kurt V. Mikkelsen
- Department of Chemistry, University of Copenhagen, Copenhagen, Denmark
| | - Dmytro Bykov
- Oak Ridge National Laboratory, Oak Ridge, TN, United States
| |
Collapse
|
6
|
Corzo HH, Hillers-Bendtsen AE, Barnes A, Zamani AY, Pawłowski F, Olsen J, Jørgensen P, Mikkelsen KV, Bykov D. Coupled cluster theory on modern heterogeneous supercomputers. Front Chem 2023; 11:1154526. [PMID: 37388945 PMCID: PMC10303140 DOI: 10.3389/fchem.2023.1154526] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2023] [Accepted: 05/11/2023] [Indexed: 07/01/2023] Open
Abstract
This study examines the computational challenges in elucidating intricate chemical systems, particularly through ab-initio methodologies. This work highlights the Divide-Expand-Consolidate (DEC) approach for coupled cluster (CC) theory-a linear-scaling, massively parallel framework-as a viable solution. Detailed scrutiny of the DEC framework reveals its extensive applicability for large chemical systems, yet it also acknowledges inherent limitations. To mitigate these constraints, the cluster perturbation theory is presented as an effective remedy. Attention is then directed towards the CPS (D-3) model, explicitly derived from a CC singles parent and a doubles auxiliary excitation space, for computing excitation energies. The reviewed new algorithms for the CPS (D-3) method efficiently capitalize on multiple nodes and graphical processing units, expediting heavy tensor contractions. As a result, CPS (D-3) emerges as a scalable, rapid, and precise solution for computing molecular properties in large molecular systems, marking it an efficient contender to conventional CC models.
Collapse
Affiliation(s)
| | | | | | - Abdulrahman Y. Zamani
- Department of Chemistry and Biochemistry and Center for Chemical Computation and Theory, University of California, Merced, CA, United States
| | - Filip Pawłowski
- Department of Chemistry and Biochemistry, Auburn University, Auburn, AL, United States
| | - Jeppe Olsen
- Department of Chemistry, Aarhus University, Aarhus, Denmark
| | - Poul Jørgensen
- Department of Chemistry, Aarhus University, Aarhus, Denmark
| | - Kurt V. Mikkelsen
- Department of Chemistry, University of Copenhagen, Copenhagen, Denmark
| | - Dmytro Bykov
- Oak Ridge National Laboratory, Oak Ridge, TN, United States
| |
Collapse
|
7
|
Nash HW, Shaw RA, Hill JG. Correlation consistent auxiliary basis sets in density fitting Hartree-Fock: The atoms sodium through argon revisited. J Comput Chem 2023; 44:1119-1128. [PMID: 36636897 PMCID: PMC10947126 DOI: 10.1002/jcc.27069] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2022] [Revised: 11/29/2022] [Accepted: 12/14/2022] [Indexed: 01/14/2023]
Abstract
We present a series of auxiliary basis sets, for the elements Na to Ar, for use in density-fitted Hartree-Fock calculations with the correlation consistent cc-pV(n + d)Z orbital basis sets. Benchmarking on total molecular energies, reaction energies and the spectroscopic constants of the SO molecule demonstrate that the new sets address the deficiencies of using existing auxiliary sets in combination with these orbital basis sets. We also report auxiliary basis sets for Na and Mg matched to cc-pVnZ, along with recommendations for pairing auxiliary basis sets to the cc-pVnZ-F12 basis sets for Hartree-Fock calculations.
Collapse
Affiliation(s)
- Harry W. Nash
- Department of ChemistryUniversity of SheffieldSheffieldUK
| | - Robert A. Shaw
- Department of ChemistryUniversity of SheffieldSheffieldUK
| | | |
Collapse
|
8
|
Hesselmann A, Werner HJ, Knowles PJ. Thermochemical evaluation of adaptive and fixed density functional theory quadrature schemes. J Chem Phys 2022; 157:234106. [PMID: 36550055 DOI: 10.1063/5.0119622] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
A systematic study is made of the accuracy and efficiency of a number of existing quadrature schemes for molecular Kohn-Sham Density-Functional Theory (DFT) using 408 molecules and 254 chemical reactions. Included are the fixed SG-x (x = 0-3) grids of Gill et al., Dasgupta, and Herbert, the 3-zone grids of Treutler and Ahlrichs, a fixed five-zone grid implemented in Molpro, and a new adaptive grid scheme. While all methods provide a systematic reduction of errors upon extension of the grid sizes, significant differences are observed in the accuracies for similar grid sizes with various approaches. For the tests in this work, the SG-x fixed grids are less suitable to achieve high accuracies in the DFT integration, while our new adaptive grid performed best among the schemes studied in this work. The extra computational time to generate the adaptive grid scales linearly with molecular size and is negligible compared with the time needed for the self-consistent field iterations for large molecules. A comparison of the grid accuracies using various density functionals shows that meta-GGA functionals need larger integration grids than GGA functionals to reach the same degree of accuracy, confirming previous investigations of the numerical stability of meta-GGA functionals. On the other hand, the grid integration errors are almost independent of the basis set, and the basis set errors are mostly much larger than the errors caused by the numerical integrations, even when using the smallest grids tested in this work.
Collapse
Affiliation(s)
- Andreas Hesselmann
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - Hans-Joachim Werner
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - Peter J Knowles
- School of Chemistry, Cardiff University, Main Building, Park Place, Cardiff CF10 3AT, United Kingdom
| |
Collapse
|
9
|
Mehta N, Martin JML. Reduced-Scaling Double Hybrid Density Functional Theory with Rapid Basis Set Convergence through Localized Pair Natural Orbital F12. J Phys Chem Lett 2022; 13:9332-9338. [PMID: 36178852 PMCID: PMC9575149 DOI: 10.1021/acs.jpclett.2c02620] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2022] [Accepted: 09/26/2022] [Indexed: 06/16/2023]
Abstract
Following earlier work [Mehta, N.; Martin, J. M. L. J. Chem. Theory Comput.2022, 10.1021/acs.jctc.2c00426] that showed how the slow basis set convergence of the double hybrid density functional theory can be obviated by the use of F12 explicit correlation in the GLPT2 step (second order Görling-Levy perturbation theory), we demonstrate here for the very large and chemically diverse GMTKN55 benchmark suite that the CPU time scaling of this step can be reduced (asymptotically linearized) using the localized pair natural orbital (PNO-L) approximation at negligible cost in accuracy.
Collapse
|
10
|
Mehta N, Martin JML. Explicitly Correlated Double-Hybrid DFT: A Comprehensive Analysis of the Basis Set Convergence on the GMTKN55 Database. J Chem Theory Comput 2022; 18:5978-5991. [PMID: 36099641 PMCID: PMC9558368 DOI: 10.1021/acs.jctc.2c00426] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
Double-hybrid density functional theory (DHDFT) offers
a pathway
to accuracy approaching composite wavefunction approaches such as
G4 theory. However, the Görling–Levy second-order perturbation
theory (GLPT2) term causes them to partially inherit the slow ∝L–3 (with L the maximum
angular momentum) basis set convergence of correlated wavefunction
methods. This could potentially be remedied by introducing F12 explicit
correlation: we investigate the basis set convergence of both DHDFT
and DHDFT-F12 (where GLPT2 is replaced by GLPT2-F12) for the large
and chemically diverse general main-group thermochemistry, kinetics,
and noncovalent interactions (GMTKN55) benchmark suite. The B2GP-PLYP-D3(BJ)
and revDSD-PBEP86-D4 DHDFs are investigated as test cases, together
with orbital basis sets as large as aug-cc-pV5Z and F12 basis sets
as large as cc-pVQZ-F12. We show that F12 greatly accelerates basis
set convergence of DHDFs, to the point that even the modest cc-pVDZ-F12
basis set is closer to the basis set limit than cc-pV(Q+d)Z or def2-QZVPPD
in orbital-based approaches, and in fact comparable in quality to
cc-pV(5+d)Z. Somewhat surprisingly, aug-cc-pVDZ-F12 is not required
even for the anionic subsets. In conclusion, DHDF-F12/VDZ-F12 eliminates
concerns about basis set convergence in both the development and applications
of double-hybrid functionals. Mass storage and I/O bottlenecks for
larger systems can be circumvented by localized pair natural orbital
approximations, which also exhibit much gentler system size scaling.
Collapse
Affiliation(s)
- Nisha Mehta
- 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
|
11
|
Bensberg M, Neugebauer J. Orbital Pair Selection for Relative Energies in the Domain-Based Local Pair Natural Orbital Coupled-Cluster Method. J Chem Phys 2022; 157:064102. [DOI: 10.1063/5.0100010] [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
For the accurate computation of relative energies, domain-based local pair natural orbital coupled-cluster [DLPNO-CCSD(T0)] has become increasingly popular. Even though DLPNO-CCSD(T0) shows a formally linear scaling of the computational effort with the system size, accurate predictions of relative energies remain costly. Therefore, multi-level approaches are attractive that focus the available computational resources on a minor part of the molecular system, e.g., a reaction center, where changes in the correlation energy are expected to be the largest. We present a pair-selected multi-level DLPNO-CCSD(T0) ansatz that automatically partitions the orbital pairs according to their contribution to the overall correlation energy change in a chemical reaction. To this end, the localized orbitals are mapped between structures in the reaction; all pair energies are approximated through computationally efficient semi-canonical second-order Møller--Plesser perturbation theory, and the orbital pairs for which the pair energies change significantly are identified. This multi-level approach is significantly more robust than our previously suggested, orbital selection-based multi-level DLPNO-CCSD(T0) ansatz [ J. Chem. Phys. 2021, 155, 224102] for reactions showing only small changes in the occupied orbitals. At the same time, it is even more efficient without added input complexity or accuracy loss compared to the full DLPNO-CCSD(T0) calculation. We demonstrate the accuracy of the multi-level approach for a total of 128 chemical reactions and potential energy curves of weakly interacting complexes from the S66x8 benchmark set.
Collapse
Affiliation(s)
- Moritz Bensberg
- Westfälische Wilhelms-Universität Münster Fachbereich 12 Chemie und Pharmazie, Germany
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut, Westfälische Wilhelms-Universität Münster Fachbereich 12 Chemie und Pharmazie, Germany
| |
Collapse
|
12
|
Prasad VK, Otero-de-la-Roza A, DiLabio GA. Small-Basis Set Density-Functional Theory Methods Corrected with Atom-Centered Potentials. J Chem Theory Comput 2022; 18:2913-2930. [PMID: 35412817 DOI: 10.1021/acs.jctc.2c00036] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Density functional theory (DFT) is currently the most popular method for modeling noncovalent interactions and thermochemistry. The accurate calculation of noncovalent interaction energies, reaction energies, and barrier heights requires choosing an appropriate functional and, typically, a relatively large basis set. Deficiencies of the density-functional approximation and the use of a limited basis set are the leading sources of error in the calculation of noncovalent and thermochemical properties in molecular systems. In this article, we present three new DFT methods based on the BLYP, M06-2X, and CAM-B3LYP functionals in combination with the 6-31G* basis set and corrected with atom-centered potentials (ACPs). ACPs are one-electron potentials that have the same form as effective-core potentials, except they do not replace any electrons. The ACPs developed in this work are used to generate energy corrections to the underlying DFT/basis-set method such that the errors in predicted chemical properties are minimized while maintaining the low computational cost of the parent methods. ACPs were developed for the elements H, B, C, N, O, F, Si, P, S, and Cl. The ACP parameters were determined using an extensive training set of 118655 data points, mostly of complete basis set coupled-cluster level quality. The target molecular properties for the ACP-corrected methods include noncovalent interaction energies, molecular conformational energies, reaction energies, barrier heights, and bond separation energies. The ACPs were tested first on the training set and then on a validation set of 42567 additional data points. We show that the ACP-corrected methods can predict the target molecular properties with accuracy close to complete basis set wavefunction theory methods, but at a computational cost of double-ζ DFT methods. This makes the new BLYP/6-31G*-ACP, M06-2X/6-31G*-ACP, and CAM-B3LYP/6-31G*-ACP methods uniquely suited to the calculation of noncovalent, thermochemical, and kinetic properties in large molecular systems.
Collapse
Affiliation(s)
- Viki Kumar Prasad
- Department of Chemistry, University of British Columbia, Okanagan, 3247 University Way, Kelowna, British Columbia V1V 1V7, Canada
| | - Alberto Otero-de-la-Roza
- Departamento de Química Física y Analítica, Facultad de Química, Universidad de Oviedo, MALTA Consolider Team, Oviedo E-33006, Spain
| | - Gino A DiLabio
- Department of Chemistry, University of British Columbia, Okanagan, 3247 University Way, Kelowna, British Columbia V1V 1V7, Canada
| |
Collapse
|
13
|
Seeger ZL, Izgorodina EI. A DLPNO-CCSD(T) benchmarking study of intermolecular interactions of ionic liquids. J Comput Chem 2022; 43:106-120. [PMID: 34687062 DOI: 10.1002/jcc.26776] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2021] [Revised: 09/14/2021] [Accepted: 09/15/2021] [Indexed: 11/11/2022]
Abstract
The accuracy of correlation energy recovered by coupled cluster single-, double-, and perturbative triple-excitations, CCSD(T), has led to the method being considered the gold standard of computational chemistry. The application of CCSD(T) has been limited to medium-sized molecular systems due to its steep scaling with molecular size. The recent development of alternative domain-based local pair natural orbital coupled-cluster method, DLPNO-CCSD(T), has significantly broadened the range of chemical systems to which CCSD(T) level calculations can be applied. Condensed systems such as ionic liquids (ILs) have a large contribution from London dispersion forces of up to 150 kJ mol-1 in large-scale clusters. Ionic liquids show appreciable charge transfer effects that result in the increased valence orbital delocalization over the entire ionic network, raising the question whether the application of methods based on localized orbitals is reliable for these semi-Coulombic materials. Here the performance of DLPNO-CCSD(T) is validated for the prediction of correlation interaction energies of two data sets incorporating single-ion pairs of protic and aprotic ILs. DLPNO-CCSD(T) produced results within chemical accuracy with tight parameter settings and a non-iterative treatment of triple excitations. To achieve spectroscopic accuracy of 1 kJ mol-1 , especially for hydrogen-bonded ILs and those containing halides, the DLPNO settings had to be increased by two orders of magnitude and include the iterative treatment of triple excitations, resulting in a 2.5-fold increase in computational cost. Two new sets of parameters are put forward to produce the performance of DLPNO-CCSD(T) within chemical and spectroscopic accuracy.
Collapse
Affiliation(s)
- Zoe L Seeger
- School of Chemistry, Monash University, Clayton, Victoria, Australia
| | | |
Collapse
|
14
|
Bensberg M, Neugebauer J. Direct orbital selection within the domain-based local pair natural orbital coupled-cluster method. J Chem Phys 2021; 155:224102. [PMID: 34911318 DOI: 10.1063/5.0071347] [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/14/2022] Open
Abstract
Domain-based local pair natural orbital coupled cluster (DLPNO-CC) has become increasingly popular to calculate relative energies (e.g., reaction energies and reaction barriers). It can be applied within a multi-level DLPNO-CC-in-DLPNO-CC ansatz to reduce the computational cost and focus the available computational resources on a specific subset of the occupied orbitals. We demonstrate how this multi-level DLPNO-CC ansatz can be combined with our direct orbital selection (DOS) approach [M. Bensberg and J. Neugebauer, J. Chem. Phys. 150, 214106 (2019)] to automatically select orbital sets for any multi-level calculation. We find that the parameters for the DOS procedure can be chosen conservatively such that they are transferable between reactions. The resulting automatic multi-level DLPNO-CC method requires no user input and is extremely robust and accurate. The computational cost is easily reduced by a factor of 3 without sacrificing accuracy. We demonstrate the accuracy of the method for a total of 61 reactions containing up to 174 atoms and use it to predict the relative stability of conformers of a Ru-based catalyst.
Collapse
Affiliation(s)
- Moritz Bensberg
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 36, 48149 Münster, Germany
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 36, 48149 Münster, Germany
| |
Collapse
|
15
|
Verma P, Huntington L, Coons MP, Kawashima Y, Yamazaki T, Zaribafiyan A. Scaling up electronic structure calculations on quantum computers: The frozen natural orbital based method of increments. J Chem Phys 2021; 155:034110. [DOI: 10.1063/5.0054647] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
Affiliation(s)
- Prakash Verma
- 1QB Information Technologies (1QBit), 200-1285 W Pender St., Vancouver, British Columbia V6E 4B1, Canada
| | - Lee Huntington
- 1QB Information Technologies (1QBit), 200-1285 W Pender St., Vancouver, British Columbia V6E 4B1, Canada
| | - Marc P. Coons
- Dow, Core R&D, Chemical Science, 1776 Building, Midland, Michigan 48674, USA
| | - Yukio Kawashima
- 1QB Information Technologies (1QBit), 200-1285 W Pender St., Vancouver, British Columbia V6E 4B1, Canada
| | - Takeshi Yamazaki
- 1QB Information Technologies (1QBit), 200-1285 W Pender St., Vancouver, British Columbia V6E 4B1, Canada
| | - Arman Zaribafiyan
- 1QB Information Technologies (1QBit), 200-1285 W Pender St., Vancouver, British Columbia V6E 4B1, Canada
| |
Collapse
|
16
|
Lipp A, Badir SO, Dykstra R, Gutierrez O, Molander GA. Catalyst-Free Decarbonylative Trifluoromethylthiolation Enabled by Electron Donor-Acceptor Complex Photoactivation. Adv Synth Catal 2021; 363:3507-3520. [PMID: 35273472 PMCID: PMC8903066 DOI: 10.1002/adsc.202100469] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2021] [Indexed: 08/06/2023]
Abstract
A catalyst- and additive-free decarbonylative trifluoromethylthiolation of aldehyde feedstocks has been developed. This operationally simple, scalable, and open-to-air transformation is driven by the selective photoexcitation of electron donor-acceptor (EDA) complexes, stemming from the association of 1,4-dihydropyridines (donor) with N-(trifluoromethylthio)phthalimide (acceptor), to trigger intermolecular single-electron transfer events under ambient- and visible light-promoted conditions. Extension to other electron acceptors enables the synthesis of thiocyanates and thioesters, as well as the difunctionalization of [1.1.1] propellane. The mechanistic intricacies of this photochemical paradigm are elucidated through a combination of experimental efforts and high-level quantum mechanical calculations [dispersion-corrected (U)DFT, DLPNO-CCSD(T), and TD-DFT]. This comprehensive study highlights the necessity for EDA complexation for efficient alkyl radical generation. Computation of subsequent ground state pathways reveals that SH2 addition of the alkyl radical to the intermediate radical EDA complex is extremely exergonic and results in a charge transfer event from the dihydropyridine donor to the N-(trifluoromethylthio)phthalimide acceptor of the EDA complex. Experimental and computational results further suggest that product formation also occurs via SH2 reaction of alkyl radicals with 1,2-bis(trifluoromethyl)disulfane, generated in-situ through combination of thiyl radicals.
Collapse
Affiliation(s)
- Alexander Lipp
- Roy and Diana Vagelos Laboratories, Department of Chemistry, University of Pennsylvania, 231 South 34th Street, Philadelphia, Pennsylvania 19104-6323, United States
| | - Shorouk O Badir
- Roy and Diana Vagelos Laboratories, Department of Chemistry, University of Pennsylvania, 231 South 34th Street, Philadelphia, Pennsylvania 19104-6323, United States
| | - Ryan Dykstra
- Department of Chemistry and Biochemistry, University of Maryland, College Park, Maryland 20742, United States
| | - Osvaldo Gutierrez
- Department of Chemistry and Biochemistry, University of Maryland, College Park, Maryland 20742, United States
| | - Gary A Molander
- Roy and Diana Vagelos Laboratories, Department of Chemistry, University of Pennsylvania, 231 South 34th Street, Philadelphia, Pennsylvania 19104-6323, United States
| |
Collapse
|
17
|
Ma Q, Werner HJ. Scalable Electron Correlation Methods. 8. Explicitly Correlated Open-Shell Coupled-Cluster with Pair Natural Orbitals PNO-RCCSD(T)-F12 and PNO-UCCSD(T)-F12. J Chem Theory Comput 2021; 17:902-926. [PMID: 33405921 DOI: 10.1021/acs.jctc.0c01129] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Abstract
We present explicitly correlated open-shell pair natural orbital local coupled-cluster methods, PNO-RCCSD(T)-F12 and PNO-UCCSD(T)-F12. The methods are extensions of our previously reported PNO-R/UCCSD methods (J. Chem. Theory Comput., 2020, 16, 3135-3151, https://pubs.acs.org/doi/10.1021/acs.jctc.0c00192) with additions of explicit correlation and perturbative triples corrections. The explicit correlation treatment follows the spin-orbital CCSD-F12b theory using Ansatz 3*A, which is found to yield comparable or better basis set convergence than the more rigorous Ansatz 3C in computed ionization potentials and reaction energies using double- to quaduple-ζ basis sets. The perturbative triples correction is adapted from the spin-orbital (T) theory to use triples natural orbitals (TNOs). To address the coupling due to off-diagonal Fock matrix elements, the local triples amplitudes are iteratively solved using small domains of TNOs, and a semicanonical (T0) domain correction with larger domains is applied to reduce the domain errors. The performance of the methods is demonstrated through benchmark calculations on ionization potentials, radical stabilization energies, reaction energies of fragmentations and rearrangements in radical cations, and spin-state energy differences of iron complexes. For a few test sets where canonical calculations are feasible, PNO-RCCSD(T)-F12 results agree with the canonical ones to within 0.4 kcal mol-1, and this maximum error is reduced to below 0.2 kcal mol-1 when large local domains are used. For larger systems, results using different thresholds for the local approximations are compared to demonstrate that 1 kcal mol-1 level of accuracy can be achieved using our default settings. For a couple of difficult cases, it is demonstrated that the errors from individual approximations are only a fraction of 1 kcal mol-1, and the overall accuracy of the method does not rely on error compensations. In contrast to canonical calculations, the use of spin-orbitals does not lead to a significant increase of computational time and memory usage in the most expensive steps of PNO-R/UCCSD(T)-F12 calculations. The only exception is the iterative solution of the (T) amplitudes, which can be avoided without significant errors by using a perturbative treatment of the off-diagonal coupling, known as (T1) approximation. For most systems, even the semicanonical approximation (T0) leads only to small errors in relative energies. Our program is well parallelized and capable of computing accurate correlation energies for molecules with 100-200 atoms using augmented triple-ζ basis sets in less than a day of elapsed time on a small computer cluster.
Collapse
Affiliation(s)
- Qianli Ma
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - Hans-Joachim Werner
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| |
Collapse
|
18
|
Mehta N, Goerigk L. Assessing the Applicability of the Geometric Counterpoise Correction in B2PLYP/Double-ζ Calculations for Thermochemistry, Kinetics, and Noncovalent Interactions. Aust J Chem 2021. [DOI: 10.1071/ch21133] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
Abstract
We present a proof-of-concept study of the suitability of Kruse and Grimme’s geometric counterpoise correction (gCP) for basis set superposition errors (BSSEs) in double-hybrid density functional calculations with a double-ζ basis set. The gCP approach only requires geometrical information as an input and no orbital/density information is needed. Therefore, this correction is practically free of any additional cost. gCP is trained against the Boys and Bernardi counterpoise correction across a set of 528 noncovalently bound dimers. We investigate the suitability of the approach for the B2PLYP/def2-SVP level of theory, and reveal error compensation effects—missing London dispersion and the BSSE—associated with B2PLYP/def2-SVP calculations, and present B2PLYP-gCP-D3(BJ)/def2-SVP with the reparametrised DFT-D3(BJ) and gCP corrections as a more balanced alternative. Benchmarking results on the S66x8 benchmark set for noncovalent interactions and the GMTKN55 database for main-group thermochemistry, kinetics, and noncovalent interactions show a statistical improvement of the B2PLYP-gCP-D3(BJ) scheme over plain B2PLYP and B2PLYP-D3(BJ). B2PLYP-D3(BJ) shows significant overestimation of interaction energies, barrier heights with larger deviations from the reference values, and wrong relative stabilities in conformers, all of which can be associated with BSSE. We find that the gCP-corrected method represents a significant improvement over B2PLYP-D3(BJ), particularly for intramolecular noncovalent interactions. These findings encourage future developments of efficient double-hybrid DFT strategies that can be applied when double-hybrid calculations with large basis sets are not feasible due to system size.
Collapse
|
19
|
Williams ZM, Wiles TC, Manby FR. Accurate Hybrid Density Functionals with UW12 Correlation. J Chem Theory Comput 2020; 16:6176-6194. [PMID: 32820918 DOI: 10.1021/acs.jctc.0c00442] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
In previous work, we suggested a single-parameter hybrid functional containing a novel correlation contribution based on the Unsöld approximation, UW12. This model resembles the explicitly correlated part of MP2-F12 theory and can be written as an explicit formula in terms of the single-particle reduced density matrix. Here, we further investigate hybrid functionals containing UW12 correlation and in particular look at functionals with a large fraction of exact exchange to reduce the self-interaction error. We suggest two new hybrid functionals B-LYP-osUW12 and fB-LYP-osUW12. On the test sets we use, our best hybrid functional overall (B-LYP-osUW12) is of similar accuracy to the best double hybrids considered while eliminating the need for virtual orbitals.
Collapse
Affiliation(s)
- Zack M Williams
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, U.K
| | - Timothy C Wiles
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, U.K
| | - Frederick R Manby
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, U.K
| |
Collapse
|
20
|
Frank MS, Schmitz G, Hättig C. Implementation of the iterative triples model CC3 for excitation energies using pair natural orbitals and Laplace transformation techniques. J Chem Phys 2020; 153:034109. [PMID: 32716174 DOI: 10.1063/5.0012597] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a pair natural orbital (PNO)-based implementation of CC3 excitation energies, which extends our previously published state-specific PNO ansatz for the solution of the excited state eigenvalue problem to methods including connected triple excitations. A thorough analysis of the equations for the excited state triples amplitudes is presented from which we derive a suitable state-specific triple natural orbital basis for the excited state triples amplitudes, which performs equally well for local and non-local excitations. The accuracy of the implementation is evaluated using a large and diverse test set. We find that for states with small contributions from double excitations, a T0 approximation to PNO-CC3 yields accurate results with a mean absolute error (MAE) for TPNO = 10-7 in the range of 0.02 eV. However, for states with larger double excitation contributions, the T0 approximation is found to yield significantly less accurate results, while the Laplace-transformed variant of PNO-CC3 shows a uniform accuracy for singly and doubly excited states (MAE and maximum error of 0.01 eV and 0.07 eV for TPNO = 10-7, respectively). Finally, we apply PNO-CC3 to the calculation of the first excited state of berenil at a S1 minimum geometry, which is shown to be close to a conical intersection. This calculation in the aug-cc-pVTZ basis set (more than 1300 basis functions) is the largest calculation ever performed with CC3 on excitation energies.
Collapse
Affiliation(s)
- Marius S Frank
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, D-44801 Bochum, Germany
| | - Gunnar Schmitz
- Deparment of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark
| | - Christof Hättig
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, D-44801 Bochum, Germany
| |
Collapse
|
21
|
Förster A, Visscher L. Double hybrid DFT calculations with Slater type orbitals. J Comput Chem 2020; 41:1660-1684. [PMID: 32297682 PMCID: PMC7317772 DOI: 10.1002/jcc.26209] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2020] [Revised: 03/31/2020] [Accepted: 04/01/2020] [Indexed: 12/14/2022]
Abstract
On a comprehensive database with 1,644 datapoints, covering several aspects of main-group as well as of transition metal chemistry, we assess the performance of 60 density functional approximations (DFA), among them 36 double hybrids (DH). All calculations are performed using a Slater type orbital (STO) basis set of triple-ζ (TZ) quality and the highly efficient pair atomic resolution of the identity approach for the exchange- and Coulomb-term of the KS matrix (PARI-K and PARI-J, respectively) and for the evaluation of the MP2 energy correction (PARI-MP2). Employing the quadratic scaling SOS-AO-PARI-MP2 algorithm, DHs based on the spin-opposite-scaled (SOS) MP2 approximation are benchmarked against a database of large molecules. We evaluate the accuracy of STO/PARI calculations for B3LYP as well as for the DH B2GP-PLYP and show that the combined basis set and PARI-error is comparable to the one obtained using the well-known def2-TZVPP Gaussian-type basis set in conjunction with global density fitting. While quadruple-ζ (QZ) calculations are currently not feasible for PARI-MP2 due to numerical issues, we show that, on the TZ level, Jacob's ladder for classifying DFAs is reproduced. However, while the best DHs are more accurate than the best hybrids, the improvements are less pronounced than the ones commonly found on the QZ level. For conformers of organic molecules and noncovalent interactions where very high accuracy is required for qualitatively correct results, DHs provide only small improvements over hybrids, while they still excel in thermochemistry, kinetics, transition metal chemistry and the description of strained organic systems.
Collapse
Affiliation(s)
- Arno Förster
- Theoretical ChemistryVrije UniversiteitAmsterdamThe Netherlands
| | - Lucas Visscher
- Theoretical ChemistryVrije UniversiteitAmsterdamThe Netherlands
| |
Collapse
|
22
|
Sirirak J, Lawan N, Van der Kamp MW, Harvey JN, Mulholland AJ. Benchmarking quantum mechanical methods for calculating reaction energies of reactions catalyzed by enzymes. PEERJ PHYSICAL CHEMISTRY 2020. [DOI: 10.7717/peerj-pchem.8] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
Abstract
To assess the accuracy of different quantum mechanical methods for biochemical modeling, the reaction energies of 20 small model reactions (chosen to represent chemical steps catalyzed by commonly studied enzymes) were calculated. The methods tested included several popular Density Functional Theory (DFT) functionals, second-order Møller Plesset perturbation theory (MP2) and its spin-component scaled variant (SCS-MP2), and coupled cluster singles and doubles and perturbative triples (CCSD(T)). Different basis sets were tested. CCSD(T)/aug-cc-pVTZ results for all 20 reactions were used to benchmark the other methods. It was found that MP2 and SCS-MP2 reaction energy calculation results are similar in quality to CCSD(T) (mean absolute error (MAE) of 1.2 and 1.3 kcal mol−1, respectively). MP2 calculations gave a large error in one case, and are more subject to basis set effects, so in general SCS-MP2 calculations are a good choice when CCSD(T) calculations are not feasible. Results with different DFT functionals were of reasonably good quality (MAEs of 2.5–5.1 kcal mol−1), whereas popular semi-empirical methods (AM1, PM3, SCC-DFTB) gave much larger errors (MAEs of 11.6–14.6 kcal mol−1). These results should be useful in guiding methodological choices and assessing the accuracy of QM/MM calculations on enzyme-catalyzed reactions.
Collapse
Affiliation(s)
- Jitnapa Sirirak
- Department of Chemistry, Faculty of Science, Silpakorn University, Nakhon Pathom, Thailand
| | - Narin Lawan
- Department of Chemistry, Faculty of Science, Chiang Mai University, Chiang Mai, Thailand
| | | | | | - Adrian J. Mulholland
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, United Kingdom
| |
Collapse
|
23
|
Balasubramani SG, Chen GP, Coriani S, Diedenhofen M, Frank MS, Franzke YJ, Furche F, Grotjahn R, Harding ME, Hättig C, Hellweg A, Helmich-Paris B, Holzer C, Huniar U, Kaupp M, Marefat Khah A, Karbalaei Khani S, Müller T, Mack F, Nguyen BD, Parker SM, Perlt E, Rappoport D, Reiter K, Roy S, Rückert M, Schmitz G, Sierka M, Tapavicza E, Tew DP, van Wüllen C, Voora VK, Weigend F, Wodyński A, Yu JM. TURBOMOLE: Modular program suite for ab initio quantum-chemical and condensed-matter simulations. J Chem Phys 2020; 152:184107. [PMID: 32414256 PMCID: PMC7228783 DOI: 10.1063/5.0004635] [Citation(s) in RCA: 523] [Impact Index Per Article: 130.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2020] [Accepted: 04/07/2020] [Indexed: 01/30/2023] Open
Abstract
TURBOMOLE is a collaborative, multi-national software development project aiming to provide highly efficient and stable computational tools for quantum chemical simulations of molecules, clusters, periodic systems, and solutions. The TURBOMOLE software suite is optimized for widely available, inexpensive, and resource-efficient hardware such as multi-core workstations and small computer clusters. TURBOMOLE specializes in electronic structure methods with outstanding accuracy-cost ratio, such as density functional theory including local hybrids and the random phase approximation (RPA), GW-Bethe-Salpeter methods, second-order Møller-Plesset theory, and explicitly correlated coupled-cluster methods. TURBOMOLE is based on Gaussian basis sets and has been pivotal for the development of many fast and low-scaling algorithms in the past three decades, such as integral-direct methods, fast multipole methods, the resolution-of-the-identity approximation, imaginary frequency integration, Laplace transform, and pair natural orbital methods. This review focuses on recent additions to TURBOMOLE's functionality, including excited-state methods, RPA and Green's function methods, relativistic approaches, high-order molecular properties, solvation effects, and periodic systems. A variety of illustrative applications along with accuracy and timing data are discussed. Moreover, available interfaces to users as well as other software are summarized. TURBOMOLE's current licensing, distribution, and support model are discussed, and an overview of TURBOMOLE's development workflow is provided. Challenges such as communication and outreach, software infrastructure, and funding are highlighted.
Collapse
Affiliation(s)
- Sree Ganesh Balasubramani
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| | - Guo P Chen
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| | - Sonia Coriani
- DTU Chemistry, Technical University of Denmark, Kemitorvet Build. 207, DK-2800 Kongens Lyngby, Denmark
| | - Michael Diedenhofen
- Dassault Systèmes Deutschland GmbH, Imbacher Weg 46, 51379 Leverkusen, Germany
| | - Marius S Frank
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44801 Bochum, Germany
| | - Yannick J Franzke
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), KIT Campus South, P.O. Box 6980, 76049 Karlsruhe, Germany
| | - Filipp Furche
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| | - Robin Grotjahn
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, 10623 Berlin, Germany
| | | | - Christof Hättig
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44801 Bochum, Germany
| | - Arnim Hellweg
- Dassault Systèmes Deutschland GmbH, Imbacher Weg 46, 51379 Leverkusen, Germany
| | - Benjamin Helmich-Paris
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Christof Holzer
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), KIT Campus South, P.O. Box 6980, 76049 Karlsruhe, Germany
| | - Uwe Huniar
- Dassault Systèmes Deutschland GmbH, Imbacher Weg 46, 51379 Leverkusen, Germany
| | - Martin Kaupp
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, 10623 Berlin, Germany
| | - Alireza Marefat Khah
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44801 Bochum, Germany
| | | | - Thomas Müller
- Forschungszentrum Jülich, Jülich Supercomputer Centre, Wilhelm-Jonen Straße, 52425 Jülich, Germany
| | - Fabian Mack
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), KIT Campus South, P.O. Box 6980, 76049 Karlsruhe, Germany
| | - Brian D Nguyen
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| | - Shane M Parker
- Department of Chemistry, Case Western Reserve University, 10900 Euclid Avenue, Cleveland, Ohio 44106, USA
| | - Eva Perlt
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| | - Dmitrij Rappoport
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, USA
| | - Kevin Reiter
- Institute of Nanotechnology, Karlsruhe Institute of Technology (KIT), KIT Campus North, P.O. Box 3640, 76021 Karlsruhe, Germany
| | - Saswata Roy
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| | - Matthias Rückert
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44801 Bochum, Germany
| | - Gunnar Schmitz
- Department of Chemistry, Aarhus Universitet, Langelandsgade 140, DK-8000 Aarhus, Denmark
| | - Marek Sierka
- TURBOMOLE GmbH, Litzenhardtstraße 19, 76135 Karlsruhe, Germany
| | - Enrico Tapavicza
- Department of Chemistry and Biochemistry, California State University, Long Beach, 1250 Bellflower Boulevard, Long Beach, California 90840, USA
| | - David P Tew
- Max Planck Institute for Solid State Research, Heisenbergstaße 1, 70569 Stuttgart, Germany
| | - Christoph van Wüllen
- Fachbereich Chemie and Forschungszentrum OPTIMAS, Technische Universität Kaiserslautern, Erwin-Schrödinger-Staße 52, 67663 Kaiserslautern, Germany
| | - Vamsee K Voora
- Department of Chemical Sciences, Tata Institute of Fundamental Research, Homi Bhabha Road, Colaba, Mumbai 400005, India
| | - Florian Weigend
- Institute of Nanotechnology, Karlsruhe Institute of Technology (KIT), KIT Campus North, P.O. Box 3640, 76021 Karlsruhe, Germany
| | - Artur Wodyński
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, 10623 Berlin, Germany
| | - Jason M Yu
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| |
Collapse
|
24
|
Ma Q, Werner HJ. Scalable Electron Correlation Methods. 7. Local Open-Shell Coupled-Cluster Methods Using Pair Natural Orbitals: PNO-RCCSD and PNO-UCCSD. J Chem Theory Comput 2020; 16:3135-3151. [PMID: 32275428 DOI: 10.1021/acs.jctc.0c00192] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present well-parallelized local implementations of high-spin open-shell coupled cluster methods with single and double excitations (CCSD) using pair natural orbitals (PNOs). The methods are based on the spin-orbital coupled cluster theory using restricted open-shell Hartree-Fock (ROHF) reference functions. Two variants, namely, PNO-UCCSD and PNO-RCCSD are implemented and compared. In PNO-UCCSD, the coupled cluster amplitudes are spin-unrestricted, while in PNO-RCCSD the linear terms are spin-adapted by a spin-projection approach as described in J. Chem. Phys. 1993, 99, 5219-5227. Near linear scaling of the computational cost with the number of correlated electrons is achieved by applying domain and pair approximations. The PNOs are spin-independent and obtained using a semicanonical spin-restricted MP2 approximation with large domains of projected atomic orbitals (PAOs). The pair approximations of our previously described closed-shell PNO-LCCSD method are carefully revised so that they are compatible to the UCCSD theory, and PNO-UCCSD or PNO-RCCSD calculations for closed-shell molecules yield exactly the same results as corresponding spin-free closed-shell PNO-LCCSD calculations. The convergence of the results with respect to the thresholds and options that control the domain and pair approximations is demonstrated. It is found that large domains are required for the single excitations in open-shell calculations in order to obtain converged results. In general, the errors of relative energies caused by the local approximations can be reduced to below 1 kcal mol-1, even for difficult cases. Presently, PNO-RCCSD and PNO-UCCSD calculations for molecules with 100-200 atoms and augmented triple-ζ basis sets can be carried out in a few hours of elapsed time using ∼100 CPU cores. In addition, the program is also capable of performing distinguishable cluster (PNO-RDCSD and PNO-UDCSC) calculations. The present work is a critical step in developing fully local open-shell PNO-RCCSD(T)-F12 methods.
Collapse
Affiliation(s)
- Qianli Ma
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - Hans-Joachim Werner
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| |
Collapse
|
25
|
Werner HJ, Knowles PJ, Manby FR, Black JA, Doll K, Heßelmann A, Kats D, Köhn A, Korona T, Kreplin DA, Ma Q, Miller TF, Mitrushchenkov A, Peterson KA, Polyak I, Rauhut G, Sibaev M. The Molpro quantum chemistry package. J Chem Phys 2020; 152:144107. [PMID: 32295355 DOI: 10.1063/5.0005081] [Citation(s) in RCA: 487] [Impact Index Per Article: 121.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
Molpro is a general purpose quantum chemistry software package with a long development history. It was originally focused on accurate wavefunction calculations for small molecules but now has many additional distinctive capabilities that include, inter alia, local correlation approximations combined with explicit correlation, highly efficient implementations of single-reference correlation methods, robust and efficient multireference methods for large molecules, projection embedding, and anharmonic vibrational spectra. In addition to conventional input-file specification of calculations, Molpro calculations can now be specified and analyzed via a new graphical user interface and through a Python framework.
Collapse
Affiliation(s)
- Hans-Joachim Werner
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Peter J Knowles
- School of Chemistry, Cardiff University, Main Building, Park Place, Cardiff CF10 3AT, United Kingdom
| | - Frederick R Manby
- School of Chemistry, University of Bristol, Cantock's Close, Bristol BS8 1TS, United Kingdom
| | - Joshua A Black
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Klaus Doll
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Andreas Heßelmann
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Daniel Kats
- Max-Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Andreas Köhn
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Tatiana Korona
- Faculty of Chemistry, University of Warsaw, L. Pasteura 1 St., 02-093 Warsaw, Poland
| | - David A Kreplin
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Qianli Ma
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Thomas F Miller
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | | | - Kirk A Peterson
- Department of Chemistry, Washington State University, Pullman, Washington 99164-4630, USA
| | - Iakov Polyak
- School of Chemistry, Cardiff University, Main Building, Park Place, Cardiff CF10 3AT, United Kingdom
| | - Guntram Rauhut
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Marat Sibaev
- School of Chemistry, Cardiff University, Main Building, Park Place, Cardiff CF10 3AT, United Kingdom
| |
Collapse
|
26
|
Wappett DA, Goerigk L. Toward a Quantum-Chemical Benchmark Set for Enzymatically Catalyzed Reactions: Important Steps and Insights. J Phys Chem A 2019; 123:7057-7074. [DOI: 10.1021/acs.jpca.9b05088] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
| | - Lars Goerigk
- School of Chemistry, The University of Melbourne, Victoria 3010, Australia
| |
Collapse
|
27
|
Mihm TN, McIsaac AR, Shepherd JJ. An optimized twist angle to find the twist-averaged correlation energy applied to the uniform electron gas. J Chem Phys 2019; 150:191101. [PMID: 31117769 DOI: 10.1063/1.5091445] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023] Open
Abstract
We explore an alternative to twist averaging in order to obtain more cost-effective and accurate extrapolations to the thermodynamic limit (TDL) for coupled cluster doubles (CCD) calculations. We seek a single twist angle to perform calculations at, instead of integrating over many random points or a grid. We introduce the concept of connectivity, a quantity derived from the nonzero four-index integrals in an MP2 calculation. This allows us to find a special twist angle that provides appropriate connectivity in the energy equation, which yields results comparable to full twist averaging. This special twist angle effectively makes the finite electron number CCD calculation represent the TDL more accurately, reducing the cost of twist-averaged CCD over Ns twist angles from Ns CCD calculations to Ns MP2 calculations plus one CCD calculation.
Collapse
Affiliation(s)
- Tina N Mihm
- Department of Chemistry, University of Iowa, Iowa City, Iowa 52242-1002, USA
| | - Alexandra R McIsaac
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - James J Shepherd
- Department of Chemistry, University of Iowa, Iowa City, Iowa 52242-1002, USA
| |
Collapse
|
28
|
Pinski P, Neese F. Analytical gradient for the domain-based local pair natural orbital second order Møller-Plesset perturbation theory method (DLPNO-MP2). J Chem Phys 2019; 150:164102. [DOI: 10.1063/1.5086544] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Affiliation(s)
- Peter Pinski
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Frank Neese
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
29
|
Morgante P, Peverati R. ACCDB: A collection of chemistry databases for broad computational purposes. J Comput Chem 2018; 40:839-848. [DOI: 10.1002/jcc.25761] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/05/2018] [Revised: 10/09/2018] [Accepted: 11/11/2018] [Indexed: 01/04/2023]
Affiliation(s)
- Pierpaolo Morgante
- Chemistry Program; Florida Institute of Technology, 150 W. University Blvd.; Melbourne Florida, 32901
| | - Roberto Peverati
- Chemistry Program; Florida Institute of Technology, 150 W. University Blvd.; Melbourne Florida, 32901
| |
Collapse
|
30
|
Kozłowska J, Schwilk M, Roztoczyńska A, Bartkowiak W. Assessment of DFT for endohedral complexes' dipole moment: PNO-LCCSD-F12 as a reference method. Phys Chem Chem Phys 2018; 20:29374-29388. [PMID: 30451255 DOI: 10.1039/c8cp05928d] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We present a systematic evaluation of the performance of a wide range of exchange-correlation functionals and related dispersion correction schemes for the computation of dipole moments of endohedral complexes, formed through the encapsulation of an AB molecule (AB = LiF, HCl) inside carbon nanotubes (CNTs) of different diameter. The consistency and accuracy of (i) generalized gradient approximation, (ii) meta GGA, (iii) global hybrid, and (iv) range-separated hybrid density functionals are assessed. In total, 37 density functionals are tested. The results obtained using the highly accurate pair natural orbitals based explicitly correlated local coupled cluster singles doubles (PNO-LCCSD-F12) method of Werner and co-workers [Schwilk et al., J. Chem. Theory Comput., 2017, 13, 3650; Ma et al., J. Chem. Theory Comput., 2017, 13, 4871] with the aug-cc-pVTZ basis set serve as a reference. The static electric dipole moment is computed via the finite field response or, when possible, as the expectation value of the dipole operator. Among others, it is shown that functionals belonging to the class of range-separated hybrids, provide results closest to the coupled cluster reference data. In particular, the ωB97X as well as the M11 functional may be considered as a promising choice for computing electric properties of noncovalent endohedral complexes. On the other hand, the worst performance was found for the functionals which do not include the Hartree-Fock exchange. The analysis of both the coupled cluster and the DFT results indicates a strong coupling of dispersion and polarization that may also explain why lower level DFT methods, as well as Hartree-Fock and MP2, cannot yield dipole moments beyond a qualitative agreement with the higher order reference data. Interestingly, the much smaller and less systematically constructed basis sets of Pople of moderate size provide results of accuracy at least comparable with the extended Dunning's aug-cc-pVTZ basis set.
Collapse
Affiliation(s)
- Justyna Kozłowska
- Department of Physical and Quantum Chemistry, Wrocław University of Science and Technology, Wybrzeże Wyspiańskiego 27, PL-50370 Wrocław, Poland.
| | | | | | | |
Collapse
|
31
|
Dornbach M, Werner HJ. Analytical energy gradients for local second-order Møller-Plesset perturbation theory using intrinsic bond orbitals. Mol Phys 2018. [DOI: 10.1080/00268976.2018.1537529] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Affiliation(s)
- Mark Dornbach
- Institut für Theoretische Chemie, Universität Stuttgart, Stuttgart, Germany
| | | |
Collapse
|
32
|
Najibi A, Goerigk L. The Nonlocal Kernel in van der Waals Density Functionals as an Additive Correction: An Extensive Analysis with Special Emphasis on the B97M-V and ωB97M-V Approaches. J Chem Theory Comput 2018; 14:5725-5738. [DOI: 10.1021/acs.jctc.8b00842] [Citation(s) in RCA: 108] [Impact Index Per Article: 18.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Asim Najibi
- School of Chemistry, The University of Melbourne, Parkville, Victoria 3010, Australia
| | - Lars Goerigk
- School of Chemistry, The University of Melbourne, Parkville, Victoria 3010, Australia
| |
Collapse
|
33
|
Mehta N, Casanova-Páez M, Goerigk L. Semi-empirical or non-empirical double-hybrid density functionals: which are more robust? Phys Chem Chem Phys 2018; 20:23175-23194. [PMID: 30062343 DOI: 10.1039/c8cp03852j] [Citation(s) in RCA: 89] [Impact Index Per Article: 14.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The development of non-empirical double-hybrid density functionals (DHDFs) is a very active research area with the number of approaches in this field having increased rapidly. At the same time, there is a lack of published work that provides a fair assessment and comparison between non-empirical and semi-empirical DHDFs on an equal footing. Herein, we close this gap and present a thorough analysis of both classes of DHDFs on the large GMTKN55 benchmark database for general main-group thermochemistry, kinetics, and noncovalent interactions [Goerigk et al., Phys. Chem. Chem. Phys., 2017, 19, 32184-32215]. In total, 115 variations of dispersion-corrected and -uncorrected DHDFs are tested, which will be condensed to an in-depth assessment of 31 methods: 19 non-empirical and 12 semi-empirical DHDFs. As such, our study represents the largest DHDF study ever conducted and can serve as an important benchmark informing method developers and users alike. Our results show that semi-empirical DHDFs are the most robust density functional approximations and more reliable and accurate than non-empirical ones. In fact, some non-empirical approaches are even outperformed by hybrid approaches or even dispersion-corrected and -uncorrected MP2 and SCS-MP2. SOS0-PBE0-2-D3(BJ) is the only exception and the only non-empirical DHDF that we can safely recommend for general applicability. However, it is still outperformed by six semi-empirical DHDFs, of which we would like to particularly recommend the following five: ωB97X-2-D3(BJ), DSD-BLYP-D3(BJ), DSD-PBEP86-D3(BJ), B2NC-PLYP-D3(BJ), and B2GPPLYP-D3(BJ). Our findings seriously question current trends in the field and they highlight that novel strategies have to be found in order to outperform the currently best density functional theory methods on the market. We hope that our study can function as an important cornerstone inspiring such a change of direction in the field.
Collapse
Affiliation(s)
- Nisha Mehta
- School of Chemistry, The University of Melbourne, Parkville, Australia.
| | | | | |
Collapse
|
34
|
Ma Q, Werner H. Explicitly correlated local coupled‐cluster methods using pair natural orbitals. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2018. [DOI: 10.1002/wcms.1371] [Citation(s) in RCA: 100] [Impact Index Per Article: 16.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Affiliation(s)
- Qianli Ma
- Institute for Theoretical ChemistryUniversity of StuttgartStuttgartGermany
| | | |
Collapse
|
35
|
Fiedler B, Himmel D, Krossing I, Friedrich J. More Stable Template Localization for an Incremental Focal-Point Approach—Implementation and Application to the Intramolecular Decomposition of Tris-perfluoro- tert-butoxyalane. J Chem Theory Comput 2018; 14:557-571. [DOI: 10.1021/acs.jctc.7b00707] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Benjamin Fiedler
- Institut
für Chemie, Technische Universität Chemnitz, 09111 Chemnitz, Germany
| | - Daniel Himmel
- Institut
für Anorganische und Analytische Chemie, Albert-Ludwigs-Universität Freiburg, 79104 Freiburg, Germany
| | - Ingo Krossing
- Institut
für Anorganische und Analytische Chemie, Albert-Ludwigs-Universität Freiburg, 79104 Freiburg, Germany
| | - Joachim Friedrich
- Institut
für Chemie, Technische Universität Chemnitz, 09111 Chemnitz, Germany
| |
Collapse
|
36
|
Gould T. ‘Diet GMTKN55’ offers accelerated benchmarking through a representative subset approach. Phys Chem Chem Phys 2018; 20:27735-27739. [DOI: 10.1039/c8cp05554h] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The GMTKN55 benchmarking protocol allows comprehensive analysis and ranking of density functional approximations with diverse chemical behaviours. This work reports diet versions of GMTKN55 which reproduce key properties of the full protocol at substantially reduced numerical cost. ‘Diet GMTKN55’ can thus be used for benchmarking expensive methods, or in combination with solid state benchmarks.
Collapse
Affiliation(s)
- Tim Gould
- Qld Micro- and Nanotechnology Centre, Griffith University
- Nathan
- Australia
| |
Collapse
|
37
|
Ma Q, Werner HJ. Scalable Electron Correlation Methods. 5. Parallel Perturbative Triples Correction for Explicitly Correlated Local Coupled Cluster with Pair Natural Orbitals. J Chem Theory Comput 2017; 14:198-215. [DOI: 10.1021/acs.jctc.7b01141] [Citation(s) in RCA: 59] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Qianli Ma
- Institut für Theoretische
Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - Hans-Joachim Werner
- Institut für Theoretische
Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| |
Collapse
|
38
|
Fiedler B, Schmitz G, Hättig C, Friedrich J. Combining Accuracy and Efficiency: An Incremental Focal-Point Method Based on Pair Natural Orbitals. J Chem Theory Comput 2017; 13:6023-6042. [PMID: 29045786 DOI: 10.1021/acs.jctc.7b00654] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
In this work, we present a new pair natural orbitals (PNO)-based incremental scheme to calculate CCSD(T) and CCSD(T0) reaction, interaction, and binding energies. We perform an extensive analysis, which shows small incremental errors similar to previous non-PNO calculations. Furthermore, slight PNO errors are obtained by using TPNO = TTNO with appropriate values of 10-7 to 10-8 for reactions and 10-8 for interaction or binding energies. The combination with the efficient MP2 focal-point approach yields chemical accuracy relative to the complete basis-set (CBS) limit. In this method, small basis sets (cc-pVDZ, def2-TZVP) for the CCSD(T) part are sufficient in case of reactions or interactions, while some larger ones (e.g., (aug)-cc-pVTZ) are necessary for molecular clusters. For these larger basis sets, we show the very high efficiency of our scheme. We obtain not only tremendous decreases of the wall times (i.e., factors >102) due to the parallelization of the increment calculations as well as of the total times due to the application of PNOs (i.e., compared to the normal incremental scheme) but also smaller total times with respect to the standard PNO method. That way, our new method features a perfect applicability by combining an excellent accuracy with a very high efficiency as well as the accessibility to larger systems due to the separation of the full computation into several small increments.
Collapse
Affiliation(s)
- Benjamin Fiedler
- Institut für Chemie, Technische Universität Chemnitz , 09111 Chemnitz, Germany
| | - Gunnar Schmitz
- Institut for Kemi, Aarhus Universitet , 8000 Aarhus C, Denmark
| | - Christof Hättig
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum , 44801 Bochum, Germany
| | - Joachim Friedrich
- Institut für Chemie, Technische Universität Chemnitz , 09111 Chemnitz, Germany
| |
Collapse
|
39
|
Ma Q, Schwilk M, Köppl C, Werner HJ. Scalable Electron Correlation Methods. 4. Parallel Explicitly Correlated Local Coupled Cluster with Pair Natural Orbitals (PNO-LCCSD-F12). J Chem Theory Comput 2017; 13:4871-4896. [DOI: 10.1021/acs.jctc.7b00799] [Citation(s) in RCA: 65] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Qianli Ma
- Institut für Theoretische
Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - Max Schwilk
- Institut für Theoretische
Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - Christoph Köppl
- Institut für Theoretische
Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - Hans-Joachim Werner
- Institut für Theoretische
Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| |
Collapse
|
40
|
Schwilk M, Ma Q, Köppl C, Werner HJ. Scalable Electron Correlation Methods. 3. Efficient and Accurate Parallel Local Coupled Cluster with Pair Natural Orbitals (PNO-LCCSD). J Chem Theory Comput 2017; 13:3650-3675. [DOI: 10.1021/acs.jctc.7b00554] [Citation(s) in RCA: 100] [Impact Index Per Article: 14.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Affiliation(s)
- Max Schwilk
- Institut für Theoretische
Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - Qianli Ma
- Institut für Theoretische
Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - Christoph Köppl
- Institut für Theoretische
Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - Hans-Joachim Werner
- Institut für Theoretische
Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| |
Collapse
|
41
|
Kjærgaard T, Baudin P, Bykov D, Kristensen K, Jørgensen P. The divide–expand–consolidate coupled cluster scheme. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2017. [DOI: 10.1002/wcms.1319] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
| | - Pablo Baudin
- Department of ChemistryAarhus UniversityAarhusDenmark
| | - Dmytro Bykov
- Department of ChemistryAarhus UniversityAarhusDenmark
| | | | | |
Collapse
|
42
|
Calbo J, Sancho-García JC, Ortí E, Aragó J. DLPNO-CCSD(T) scaled methods for the accurate treatment of large supramolecular complexes. J Comput Chem 2017; 38:1869-1878. [PMID: 28558123 DOI: 10.1002/jcc.24835] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2017] [Revised: 04/10/2017] [Accepted: 04/25/2017] [Indexed: 01/19/2023]
Abstract
In this work, we present scaled variants of the DLPNO-CCSD(T) method, dubbed as (LS)DLPNO-CCSD(T) and (NS)DLPNO-CCSD(T), to obtain accurate interaction energies in supramolecular complexes governed by noncovalent interactions. The novel scaled schemes are based on the linear combination of the DLPNO-CCSD(T) correlation energies calculated with the standard (LoosePNO and NormalPNO) and modified (Loose2PNO and Normal2PNO) DLPNO-CCSD(T) accuracy levels. The scaled DLPNO-CCSD(T) variants provide nearly TightPNO accuracy, which is essential for the quantification of weak noncovalent interactions, with a noticeable saving in computational cost. Importantly, the accuracy of the proposed schemes is preserved irrespective of the nature and strength of the supramolecular interaction. The (LS)DLPNO-CCSD(T) and (NS)DLPNO-CCSD(T) protocols have been used to study in depth the role of the CH-π versus π-π interactions in the supramolecular complex formed by the electron-donor truxene-tetrathiafulvalene (truxTTF) and the electron-acceptor hemifullerene (C30 H12 ). (NS)DLPNO-CCSD(T)/CBS calculations clearly reveal the higher stability of staggered (dominated by CH-π interactions) versus bowl-in-bowl (dominated by π-π interactions) arrangements in the truxTTF•C30 H12 heterodimer. Hemifullerene and similar carbon-based buckybowls are therefore expected to self-assemble with donor compounds in a richer way other than the typical concave-convex π-π arrangement found in fullerene-based aggregates. © 2017 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Joaquín Calbo
- Instituto de Ciencia Molecular, Universidad de Valencia, Paterna, E-46980, Spain
| | - Juan C Sancho-García
- Departamento de Química Física, Universidad de Alicante, Alicante, E-03080, Spain
| | - Enrique Ortí
- Instituto de Ciencia Molecular, Universidad de Valencia, Paterna, E-46980, Spain
| | - Juan Aragó
- Instituto de Ciencia Molecular, Universidad de Valencia, Paterna, E-46980, Spain
| |
Collapse
|
43
|
Schmitz G, Hättig C. Accuracy of Explicitly Correlated Local PNO-CCSD(T). J Chem Theory Comput 2017; 13:2623-2633. [DOI: 10.1021/acs.jctc.7b00180] [Citation(s) in RCA: 33] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Gunnar Schmitz
- Department
of Chemistry, Aarhus Universitet, Aarhus 8000, Denmark
| | - Christof Hättig
- Arbeitsgruppe
Quantenchemie, Ruhr-Universität, Bochum 44780, Germany
| |
Collapse
|
44
|
Goerigk L, Hansen A, Bauer C, Ehrlich S, Najibi A, Grimme S. A look at the density functional theory zoo with the advanced GMTKN55 database for general main group thermochemistry, kinetics and noncovalent interactions. Phys Chem Chem Phys 2017; 19:32184-32215. [DOI: 10.1039/c7cp04913g] [Citation(s) in RCA: 854] [Impact Index Per Article: 122.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
We present the updated and extended GMTKN55 benchmark database for more accurate and extensive energetic evaluation of density functionals and other electronic structure methods with detailed guidelines for method users.
Collapse
Affiliation(s)
- Lars Goerigk
- School of Chemistry
- The University of Melbourne
- Parkville
- Australia
| | - Andreas Hansen
- Universität Bonn
- Mulliken Center for Theoretical Chemistry
- Bonn
- Germany
| | - Christoph Bauer
- Universität Bonn
- Mulliken Center for Theoretical Chemistry
- Bonn
- Germany
| | - Stephan Ehrlich
- Universität Bonn
- Mulliken Center for Theoretical Chemistry
- Bonn
- Germany
| | - Asim Najibi
- School of Chemistry
- The University of Melbourne
- Parkville
- Australia
| | - Stefan Grimme
- Universität Bonn
- Mulliken Center for Theoretical Chemistry
- Bonn
- Germany
| |
Collapse
|
45
|
Marianski M, Supady A, Ingram T, Schneider M, Baldauf C. Assessing the Accuracy of Across-the-Scale Methods for Predicting Carbohydrate Conformational Energies for the Examples of Glucose and α-Maltose. J Chem Theory Comput 2016; 12:6157-6168. [DOI: 10.1021/acs.jctc.6b00876] [Citation(s) in RCA: 73] [Impact Index Per Article: 9.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Mateusz Marianski
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany
| | - Adriana Supady
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany
| | - Teresa Ingram
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany
| | - Markus Schneider
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany
| | - Carsten Baldauf
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany
| |
Collapse
|
46
|
Pavošević F, Pinski P, Riplinger C, Neese F, Valeev EF. SparseMaps--A systematic infrastructure for reduced-scaling electronic structure methods. IV. Linear-scaling second-order explicitly correlated energy with pair natural orbitals. J Chem Phys 2016; 144:144109. [PMID: 27083710 DOI: 10.1063/1.4945444] [Citation(s) in RCA: 78] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a formulation of the explicitly correlated second-order Møller-Plesset (MP2-F12) energy in which all nontrivial post-mean-field steps are formulated with linear computational complexity in system size. The two key ideas are the use of pair-natural orbitals for compact representation of wave function amplitudes and the use of domain approximation to impose the block sparsity. This development utilizes the concepts for sparse representation of tensors described in the context of the domain based local pair-natural orbital-MP2 (DLPNO-MP2) method by us recently [Pinski et al., J. Chem. Phys. 143, 034108 (2015)]. Novel developments reported here include the use of domains not only for the projected atomic orbitals, but also for the complementary auxiliary basis set (CABS) used to approximate the three- and four-electron integrals of the F12 theory, and a simplification of the standard B intermediate of the F12 theory that avoids computation of four-index two-electron integrals that involve two CABS indices. For quasi-1-dimensional systems (n-alkanes), the ON DLPNO-MP2-F12 method becomes less expensive than the conventional ON(5) MP2-F12 for n between 10 and 15, for double- and triple-zeta basis sets; for the largest alkane, C200H402, in def2-TZVP basis, the observed computational complexity is N(∼1.6), largely due to the cubic cost of computing the mean-field operators. The method reproduces the canonical MP2-F12 energy with high precision: 99.9% of the canonical correlation energy is recovered with the default truncation parameters. Although its cost is significantly higher than that of DLPNO-MP2 method, the cost increase is compensated by the great reduction of the basis set error due to explicit correlation.
Collapse
Affiliation(s)
- Fabijan Pavošević
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Peter Pinski
- Max Planch Institute for Chemical Energy Conversion, Stiftstr. 34-36, D-45470 Mülheim an der Ruhr, Germany
| | - Christoph Riplinger
- Max Planch Institute for Chemical Energy Conversion, Stiftstr. 34-36, D-45470 Mülheim an der Ruhr, Germany
| | - Frank Neese
- Max Planch Institute for Chemical Energy Conversion, Stiftstr. 34-36, D-45470 Mülheim an der Ruhr, Germany
| | - Edward F Valeev
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| |
Collapse
|
47
|
Shepherd JJ. Communication: Convergence of many-body wave-function expansions using a plane-wave basis in the thermodynamic limit. J Chem Phys 2016; 145:031104. [DOI: 10.1063/1.4958461] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- James J. Shepherd
- Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Ave., Cambridge, Massachusetts 02139, USA
| |
Collapse
|
48
|
Fiedler B, Coriani S, Friedrich J. Molecular Dipole Moments within the Incremental Scheme Using the Domain-Specific Basis-Set Approach. J Chem Theory Comput 2016; 12:3040-52. [PMID: 27300371 DOI: 10.1021/acs.jctc.6b00076] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
We present the first implementation of the fully automated incremental scheme for CCSD unrelaxed dipole moments using the domain-specific basis-set approach. Truncation parameters are varied, and the accuracy of the method is statistically analyzed for a test set of 20 molecules. The local approximations introduce small errors at second order and negligible ones at third order. For a third-order incremental CCSD expansion with a CC2 error correction, a cc-pVDZ/SV domain-specific basis set (tmain = 3.5 Bohr), and the truncation parameter f = 30 Bohr, we obtain a mean error of 0.00 mau (-0.20 mau) and a standard deviation of 1.95 mau (2.17 mau) for the total dipole moments (Cartesian components of the dipole vectors). By analyzing incremental CCSD energies, we demonstrate that the MP2 and CC2 error correction schemes are an exclusive correction for the domain-specific basis-set error. Our implementation of the incremental scheme provides fully automated computations of highly accurate dipole moments at reduced computational cost and is fully parallelized in terms of the calculation of the increments. Therefore, one can utilize the incremental scheme, on the same hardware, to extend the basis set in comparison to standard CCSD and thus obtain a better total accuracy.
Collapse
Affiliation(s)
- Benjamin Fiedler
- Institute for Chemistry, Technische Universität Chemnitz , Straße der Nationen 62, D-09111 Chemnitz, Germany
| | - Sonia Coriani
- Dipartimento di Scienze Chimiche e Farmaceutiche, Università degli Studi di Trieste , Via L. Giorgieri 1, I-34127 Trieste, Italy.,Aarhus Institute of Advanced Studies, Aarhus University , Høegh-Guldbergs Gade 6B, DK-8000 Aarhus C, Denmark
| | - Joachim Friedrich
- Institute for Chemistry, Technische Universität Chemnitz , Straße der Nationen 62, D-09111 Chemnitz, Germany
| |
Collapse
|
49
|
Friedrich J, Fiedler B. Accurate calculation of binding energies for molecular clusters – Assessment of different models. Chem Phys 2016. [DOI: 10.1016/j.chemphys.2016.02.022] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
50
|
Wang YM, Hättig C, Reine S, Valeev E, Kjærgaard T, Kristensen K. Explicitly correlated second-order Møller-Plesset perturbation theory in a Divide-Expand-Consolidate (DEC) context. J Chem Phys 2016; 144:204112. [DOI: 10.1063/1.4951696] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yang Min Wang
- qLEAP Center for Theoretical Chemistry, Department of Chemistry, Aarhus University, Langelandsgade 140, DK-8000 Aarhus C, Denmark
| | - Christof Hättig
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, D-44780 Bochum, Germany
| | - Simen Reine
- Centre for Theoretical and Computational Chemistry, Department of Chemistry, University of Oslo, P.O. Box 1033, N-1315 Blindern, Norway
| | - Edward Valeev
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Thomas Kjærgaard
- qLEAP Center for Theoretical Chemistry, Department of Chemistry, Aarhus University, Langelandsgade 140, DK-8000 Aarhus C, Denmark
| | - Kasper Kristensen
- qLEAP Center for Theoretical Chemistry, Department of Chemistry, Aarhus University, Langelandsgade 140, DK-8000 Aarhus C, Denmark
| |
Collapse
|