1
|
Mester D, Kállay M. Basis Set Limit of CCSD(T) Energies: Explicit Correlation Versus Density-Based Basis-Set Correction. J Chem Theory Comput 2023; 19:8210-8222. [PMID: 37950703 PMCID: PMC10688194 DOI: 10.1021/acs.jctc.3c00979] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2023] [Revised: 10/21/2023] [Accepted: 10/23/2023] [Indexed: 11/13/2023]
Abstract
A thorough comparison is carried out for explicitly correlated and density-based basis-set correction approaches, which were primarily developed to mitigate the basis-set incompleteness error of wave function methods. An efficient implementation of the density-based scheme is also presented, utilizing the density-fitting approximation. The performance of these approaches is comprehensively tested for the second-order Møller-Plesset (MP2), coupled-cluster singles and doubles (CCSD), and CCSD with perturbative triples [CCSD(T)] methods with respect to the corresponding complete basis set references. It is demonstrated that the density-based correction together with complementary auxiliary basis set (CABS)-corrected Hartree-Fock energies is highly robust and effectively reduces the error of the standard approaches; however, it does not outperform the corresponding explicitly correlated methods. Nevertheless, what still makes the density-corrected CCSD and CCSD(T) methods competitive is that their computational costs are roughly half of those of the corresponding explicitly correlated variants. Additionally, an incremental approach for standard CCSD and CCSD(T) is introduced. In this simple scheme, the total energies are corrected with the CABS correction and explicitly correlated MP2 contributions. As demonstrated, the resulting methods yield surprisingly good results, below 1 kcal/mol for thermochemical properties even with a double-ζ basis, while their computational expenses are practically identical to those of the density-based basis-set correction approaches.
Collapse
Affiliation(s)
- Dávid Mester
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- HUN-REN-BME
Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- MTA-BME
Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| | - Mihály Kállay
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- HUN-REN-BME
Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- MTA-BME
Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| |
Collapse
|
2
|
Szabó PB, Csóka J, Kállay M, Nagy PR. Linear-Scaling Local Natural Orbital CCSD(T) Approach for Open-Shell Systems: Algorithms, Benchmarks, and Large-Scale Applications. J Chem Theory Comput 2023; 19:8166-8188. [PMID: 37921429 PMCID: PMC10687875 DOI: 10.1021/acs.jctc.3c00881] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2023] [Revised: 10/05/2023] [Accepted: 10/16/2023] [Indexed: 11/04/2023]
Abstract
The extension of the highly optimized local natural orbital (LNO) coupled cluster (CC) with single-, double-, and perturbative triple excitations [LNO-CCSD(T)] method is presented for high-spin open-shell molecules based on restricted open-shell references. The techniques enabling the outstanding efficiency of the closed-shell LNO-CCSD(T) variant are adopted, including the iteration- and redundancy-free second-order Møller-Plesset and (T) formulations as well as the integral-direct, memory- and disk use-economic, and OpenMP-parallel algorithms. For large molecules, the efficiency of our open-shell LNO-CCSD(T) method approaches that of its closed-shell parent method due to the application of restricted orbital sets for demanding integral transformations and a novel approximation for higher-order long-range spin-polarization effects. The accuracy of open-shell LNO-CCSD(T) is extensively tested for radicals and reactions thereof, ionization processes, as well as spin-state splittings, and transition-metal compounds. At the size range where the canonical CCSD(T) reference is accessible (up to 20-30 atoms), the average open-shell LNO-CCSD(T) correlation energies are found to be 99.9 to 99.95% accurate, which translates into average absolute deviations of a few tenths of kcal/mol in the investigated energy differences already with the default settings. For more extensive molecules, the local errors may grow, but they can be estimated and decreased via affordable systematic convergence studies. This enables the accurate modeling of large systems with complex electronic structures, as illustrated on open-shell organic radicals and transition-metal complexes of up to 179 atoms as well as on challenging biochemical systems, including up to 601 atoms and 11,000 basis functions. While the protein models involve difficulties for local approximations, such as the spin states of a bounded iron ion or an extremely delocalized singly occupied orbital, the corresponding single-node LNO-CCSD(T) computations were feasible in a matter of days with 10s to 100 GB of memory use. Therefore, the new LNO-CCSD(T) implementation enables highly accurate computations for open-shell systems of unprecedented size and complexity with widely accessible hardware.
Collapse
Affiliation(s)
- P. Bernát Szabó
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Műegyetem rkp. 3, H-1111 Budapest, Hungary
| | - József Csóka
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Műegyetem rkp. 3, H-1111 Budapest, Hungary
- HUN-REN-BME
Quantum Chemistry Research Group, Műegyetem rkp. 3, H-1111 Budapest, Hungary
- MTA-BME
Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3, H-1111 Budapest, Hungary
| | - Mihály Kállay
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Műegyetem rkp. 3, H-1111 Budapest, Hungary
- HUN-REN-BME
Quantum Chemistry Research Group, Műegyetem rkp. 3, H-1111 Budapest, Hungary
- MTA-BME
Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3, H-1111 Budapest, Hungary
| | - Péter R. Nagy
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Műegyetem rkp. 3, H-1111 Budapest, Hungary
- HUN-REN-BME
Quantum Chemistry Research Group, Műegyetem rkp. 3, H-1111 Budapest, Hungary
- MTA-BME
Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3, H-1111 Budapest, Hungary
| |
Collapse
|
3
|
Ten-No SL. Nonunitary projective transcorrelation theory inspired by the F12 ansatz. J Chem Phys 2023; 159:171103. [PMID: 37921247 DOI: 10.1063/5.0175337] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2023] [Accepted: 10/16/2023] [Indexed: 11/04/2023] Open
Abstract
An alternative nonunitary transcorrelation, inspired by the F12 ansatz, is investigated. In contrast to the Jastrow transcorrelation of Boys-Handy, the effective Hamiltonian of this projective transcorrelation features: 1. a series terminating formally at four-body interactions. 2. no spin-contamination within the non-relativistic framework. 3. simultaneous satisfaction of the singlet and triplet first-order cusp conditions. 4. arbitrary choices of pairs for correlation including frozen-core approximations. We discuss the connection between the projective transcorrelation and F12 theory with applications to small molecules, to show that the cusp conditions play an important role to reduce the uncertainty arising from the nonunitary transformation.
Collapse
Affiliation(s)
- Seiichiro L Ten-No
- Graduate School of System Informatics, Kobe University, Nada-ku, Kobe 657-8501, Japan
| |
Collapse
|
4
|
Valeev EF, Harrison RJ, Holmes AA, Peterson CC, Penchoff DA. Direct Determination of Optimal Real-Space Orbitals for Correlated Electronic Structure of Molecules. J Chem Theory Comput 2023; 19:7230-7241. [PMID: 37791808 DOI: 10.1021/acs.jctc.3c00732] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/05/2023]
Abstract
We demonstrate how to determine numerically nearly exact orthonormal orbitals that are optimal for the evaluation of the energy of arbitrary (correlated) states of atoms and molecules by minimization of the energy Lagrangian. Orbitals are expressed in real space using a multiresolution spectral element basis that is refined adaptively to achieve the user-specified target precision while avoiding the ill-conditioning issues that plague AO basis set expansions traditionally used for correlated models of molecular electronic structure. For light atoms, the orbital solver, in conjunction with a variational electronic structure model [selected Configuration Interaction (CI)] provides energies of comparable precision to a state-of-the-art atomic CI solver. The computed electronic energies of atoms and molecules are significantly more accurate than the counterparts obtained with the orbital sets of the same rank expanded in Gaussian AO bases, and can be determined even when linear dependence issues preclude the use of the AO bases. It is feasible to optimize more than 100 fully correlated numerical orbitals on a single computer node, and significant room exists for additional improvement. These findings suggest that real-space orbital representations might be the preferred alternative to AO representations for high-end models of correlated electronic states of molecules and materials.
Collapse
Affiliation(s)
- Edward F Valeev
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Robert J Harrison
- Department of Applied Mathematics & Statistics, Stony Brook University, Stony Brook, New York 11794, United States
| | - Adam A Holmes
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Charles C Peterson
- Office of Advanced Research Computing, University of California, Los Angeles, Los Angeles, California 90095, United States
| | - Deborah A Penchoff
- UT Innovative Computing Laboratory, University of Tennessee, Knoxville, Tennessee 37996, United States
| |
Collapse
|
5
|
Thorpe JH, Feller D, Bross DH, Ruscic B, Stanton JF. Sub 20 cm -1 computational prediction of the CH bond energy - a case of systematic error in computational thermochemistry. Phys Chem Chem Phys 2023; 25:21162-21172. [PMID: 36200428 DOI: 10.1039/d2cp03964h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
The bond dissociation energy of methylidyne, D0(CH), is studied using an improved version of the High-Accuracy Extrapolated ab initio Thermochemistry (HEAT) approach as well as the Feller-Peterson-Dixon (FPD) model chemistry. These calculations, which include basis sets up to nonuple (aug-cc-pCV9Z) quality, are expected to be capable of providing results substantially more accurate than the ca. 1 kJ mol-1 level that is characteristic of standard high-accuracy protocols for computational thermochemistry. The calculated 0 K CH bond energy (27 954 ± 15 cm-1 for HEAT and 27 956 ± 15 cm-1 for FPD), along with equivalent treatments of the CH ionization energy and the CH+ dissociation energy (85 829 ± 15 cm-1 and 32 946 ± 15 cm-1, respectively), were compared to the existing benchmarks from Active Thermochemical Tables (ATcT), uncovering an unexpected difference for D0(CH). This has prompted a detailed reexamination of the provenance of the corresponding ATcT benchmark, allowing the discovery and subsequent correction of a systematic error present in several published high-level calculations, ultimately yielding an amended ATcT benchmark for D0(CH). Finally, the current theoretical results were added to the ATcT Thermochemical Network, producing refined ATcT estimates of 27 957.3 ± 6.0 cm-1 for D0(CH), 32 946.7 ± 0.6 cm-1 for D0(CH+), and 85 831.0 ± 6.0 cm-1 for IE(CH).
Collapse
Affiliation(s)
- James H Thorpe
- The Quantum Theory Project, Department of Chemistry, University of Florida, Gainesville, Florida, 32611, USA.
| | - David Feller
- Washington State University, Pullman, Washington 99164-4630, USA
- University of Alabama, Tuscaloosa, Alabama 35487-0336, USA
| | - David H Bross
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Argonne, Illinois 60439, USA.
| | - Branko Ruscic
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Argonne, Illinois 60439, USA.
| | - John F Stanton
- The Quantum Theory Project, Department of Chemistry, University of Florida, Gainesville, Florida, 32611, USA.
| |
Collapse
|
6
|
Neese F. The SHARK integral generation and digestion system. J Comput Chem 2023; 44:381-396. [PMID: 35678278 DOI: 10.1002/jcc.26942] [Citation(s) in RCA: 17] [Impact Index Per Article: 17.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/18/2022] [Revised: 04/22/2022] [Accepted: 04/28/2022] [Indexed: 01/03/2023]
Abstract
In this paper, the SHARK integral generation and digestion engine is described. In essence, SHARK is based on a reformulation of the popular McMurchie/Davidson approach to molecular integrals. This reformulation leads to an efficient algorithm that is driven by BLAS level 3 operations. The algorithm is particularly efficient for high angular momentum basis functions (up to L = 7 is available by default, but the algorithm is programmed for arbitrary angular momenta). SHARK features a significant number of specific programming constructs that are designed to greatly simplify the workflow in quantum chemical program development and avoid undesirable code duplication to the largest possible extent. SHARK can handle segmented, generally and partially generally contracted basis sets. It can be used to generate a host of one- and two-electron integrals over various kernels including, two-, three-, and four-index repulsion integrals, integrals over Gauge Including Atomic Orbitals (GIAOs), relativistic integrals and integrals featuring a finite nucleus model. SHARK provides routines to evaluate Fock like matrices, generate integral transformations and related tasks. SHARK is the essential engine inside the ORCA package that drives essentially all tasks that are related to integrals over basis functions in version ORCA 5.0 and higher. Since the core of SHARK is based on low-level basic linear algebra (BLAS) operations, it is expected to not only perform well on present day but also on future hardware provided that the hardware manufacturer provides a properly optimized BLAS library for matrix and vector operations. Representative timings and comparisons to the Libint library used by ORCA are reported for Intel i9 and Apple M1 max processors.
Collapse
Affiliation(s)
- Frank Neese
- Department of Molecular Theory and Spectroscopy, Max Planck Institut für Kohlenforschung, Mülheim an der Ruhr, Germany
| |
Collapse
|
7
|
Kállay M, Horváth RA, Gyevi-Nagy L, Nagy PR. Basis Set Limit CCSD(T) Energies for Extended Molecules via a Reduced-Cost Explicitly Correlated Approach. J Chem Theory Comput 2022; 19:174-189. [PMID: 36576419 PMCID: PMC9835832 DOI: 10.1021/acs.jctc.2c01031] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Several approximations are introduced and tested to reduce the computational expenses of the explicitly correlated coupled-cluster singles and doubles with perturbative triples [CCSD(T)] method for both closed and open-shell species. First, the well-established frozen natural orbital (FNO) technique is adapted to explicitly correlated CC approaches. Second, our natural auxiliary function (NAF) scheme is employed to reduce the size of the auxiliary basis required for the density fitting approximation regularly used in explicitly correlated calculations. Third, a new approach, termed the natural auxiliary basis (NAB) approximation, is proposed to decrease the size of the auxiliary basis needed for the expansion of the explicitly correlated geminals. The performance of the above approximations and that of the combined FNO-NAF-NAB approach are tested for atomization and reaction energies. Our results show that overall speedups of 7-, 5-, and 3-times can be achieved with double-, triple-, and quadruple-ζ basis sets, respectively, without any loss in accuracy. The new method can provide, e.g., reaction energies and barrier heights well within chemical accuracy for molecules with more than 40 atoms within a few days using a few dozen processor cores, and calculations with 50+ atoms are still feasible. These routinely affordable computations considerably extend the reach of explicitly correlated CCSD(T).
Collapse
Affiliation(s)
- Mihály Kállay
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary,ELKH-BME
Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary,MTA-BME
Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary,
| | - Réka A. Horváth
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary,ELKH-BME
Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary,MTA-BME
Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| | - László Gyevi-Nagy
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary,ELKH-BME
Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary,MTA-BME
Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| | - Péter R. Nagy
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary,ELKH-BME
Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary,MTA-BME
Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| |
Collapse
|
8
|
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
|
9
|
Nagy PR, Gyevi-Nagy L, Lőrincz BD, Kállay M. Pursuing the basis set limit of CCSD(T) non-covalent interaction energies for medium-sized complexes: case study on the S66 compilation. Mol Phys 2022. [DOI: 10.1080/00268976.2022.2109526] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/15/2022]
Affiliation(s)
- Péter R. Nagy
- Faculty of Chemical Technology and Biotechnology, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest, Hungary
- ELKH-BME Quantum Chemistry Research Group, Budapest, Hungary
| | - László Gyevi-Nagy
- Faculty of Chemical Technology and Biotechnology, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest, Hungary
- ELKH-BME Quantum Chemistry Research Group, Budapest, Hungary
| | - Balázs D. Lőrincz
- Faculty of Chemical Technology and Biotechnology, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest, Hungary
- ELKH-BME Quantum Chemistry Research Group, Budapest, Hungary
| | - Mihály Kállay
- Faculty of Chemical Technology and Biotechnology, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest, Hungary
- ELKH-BME Quantum Chemistry Research Group, Budapest, Hungary
| |
Collapse
|
10
|
Schleich P, Kottmann JS, Aspuru-Guzik A. Improving the accuracy of the variational quantum eigensolver for molecular systems by the explicitly-correlated perturbative [2] R12-correction. Phys Chem Chem Phys 2022; 24:13550-13564. [PMID: 35638435 DOI: 10.1039/d2cp00247g] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
We provide an integration of the universal, perturbative explicitly correlated [2]R12-correction in the context of the Variational Quantum Eigensolver (VQE). This approach is able to increase the accuracy of the underlying reference method significantly while requiring no additional quantum resources. The proposed approach only requires knowledge of the one- and two-particle reduced density matrices (RDMs) of the reference wavefunction; these can be measured after having reached convergence in the VQE. This computation comes at a cost that scales as the sixth power of the number of electrons. We explore the performance of the VQE + [2]R12 approach using both conventional Gaussian basis sets and our recently proposed directly determined pair-natural orbitals obtained by multiresolution analysis (MRA-PNOs). Both Gaussian orbital and PNOs are investigated as a potential set of complementary basis functions in the computation of [2]R12. In particular the combination of MRA-PNOs with [2]R12 has turned out to be very promising - persistently throughout our data, this allowed very accurate simulations at a quantum cost of a minimal basis set. Additionally, we found that the deployment of PNOs as complementary basis can greatly reduce the number of complementary basis functions that enter the computation of the correction at a complexity.
Collapse
Affiliation(s)
- Philipp Schleich
- Department of Computer Science, University of Toronto, Canada. .,Applied and Computational Mathematics, Department of Mathematics, RWTH Aachen University, Aachen, Germany.,Vector Institute for Artificial Intelligence, Toronto, Canada
| | - Jakob S Kottmann
- Department of Computer Science, University of Toronto, Canada. .,Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Canada
| | - Alán Aspuru-Guzik
- Department of Computer Science, University of Toronto, Canada. .,Vector Institute for Artificial Intelligence, Toronto, Canada.,Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Canada.,Canadian Institute for Advanced Research (CIFAR) Lebovic Fellow, Toronto, Canada
| |
Collapse
|
11
|
Nakatsuji H, Nakashima H, Kurokawa YI. Accurate scaling functions of the scaled Schrödinger equation. J Chem Phys 2022; 156:014113. [PMID: 34998320 DOI: 10.1063/5.0077495] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The scaling function g of the scaled Schrödinger equation (SSE) is generalized to obtain accurate solutions of the Schrödinger equation (SE) with the free complement (FC) theory. The electron-nuclear and electron-electron scaling functions, giA and gij, respectively, are generalized. From the relations between SE and SSE at the inter-particle distances being zero and infinity, the scaling function must satisfy the collisional (or coalescent) condition and the asymptotic condition, respectively. Based on these conditions, general scaling functions are classified into "correct" (satisfying both conditions), "reasonable" (satisfying only collisional condition), and "approximate but still useful" (not satisfying collisional condition) classes. Several analytical scaling functions are listed for each class. Popular functions riA and rij belong to the reasonable class. The qualities of many electron-electron scaling functions are examined variationally for the helium atom using the FC theory. Although the complement functions of FC theory are produced generally from both the potential and kinetic operators in the Hamiltonian, those produced from the kinetic operator were shown to be less important than those produced from the potential operator. Hence, we used only the complement functions produced from the potential operator and showed that the correct-class gij functions gave most accurate results and the reasonable-class functions were less accurate. Among the examined correct and reasonable functions, the conventional function rij was worst in accuracy, although it was still very accurate. Thus, we have many potentially accurate "correct" scaling functions for use in FC theory to solve the SEs of atoms and molecules.
Collapse
Affiliation(s)
- Hiroshi Nakatsuji
- Quantum Chemistry Research Institute, Kyoto Technoscience Center 16, 14 Yoshida Kawara-machi, Sakyo-ku, Kyoto 606-8305, Japan
| | - Hiroyuki Nakashima
- Quantum Chemistry Research Institute, Kyoto Technoscience Center 16, 14 Yoshida Kawara-machi, Sakyo-ku, Kyoto 606-8305, Japan
| | - Yusaku I Kurokawa
- Quantum Chemistry Research Institute, Kyoto Technoscience Center 16, 14 Yoshida Kawara-machi, Sakyo-ku, Kyoto 606-8305, Japan
| |
Collapse
|
12
|
Nagy PR, Gyevi-Nagy L, Kállay M. Basis set truncation corrections for improved frozen natural orbital CCSD(T) energies. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1963495] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- Péter R. Nagy
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest, Hungary
| | - László Gyevi-Nagy
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest, Hungary
| | - Mihály Kállay
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest, Hungary
| |
Collapse
|
13
|
Chen J, Kato J, Harper JB, Shao Y, Ho J. On the Accuracy of QM/MM Models: A Systematic Study of Intramolecular Proton Transfer Reactions of Amino Acids in Water. J Phys Chem B 2021; 125:9304-9316. [PMID: 34355564 DOI: 10.1021/acs.jpcb.1c04876] [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/29/2022]
Abstract
This work presents a systematic assessment of QM/QM' and QM/MM models with respect to direct QM calculations for the tautomerization (neutral to zwitterion) reactions of amino acids (glycine, alanine, valine, aspartate, and neutral and protonated histidine) solvated in a 160 water cluster. The effect of varying QM region size and choice of embedding potentials, including fixed-charge and polarizable molecular mechanics force fields (TIP3P and EFP) and various semiempirical QM methods (PM7, GFN2-xTB, DFTBA, DFTB3, HF-3c, and PBEh-3c), on the accuracy of the models was examined. A surprising finding was that molecular mechanics force fields outperformed many of the semiempirical methods. Generally, the errors in the QM/QM' and QM/MM models converge slowly with respect to the QM region size, requiring 50 or more waters to be included in the QM region before the error in the model falls below 1 kcal mol-1 of its pure QM result. Different QM region selection schemes were also compared, and it was found that selection based on Natural Population Analysis (NPA) atomic charges significantly reduced the error in the QM/QM' and QM/MM models particularly if a low-quality embedding potential was used. It is envisaged that these results will be useful for the development of future hybrid QM models.
Collapse
Affiliation(s)
- Junbo Chen
- School of Chemistry, The University of New South Wales, Sydney, NSW 2052, Australia
| | - Jin Kato
- School of Chemistry, The University of New South Wales, Sydney, NSW 2052, Australia
| | - Jason B Harper
- School of Chemistry, The University of New South Wales, Sydney, NSW 2052, Australia
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, United States
| | - Junming Ho
- School of Chemistry, The University of New South Wales, Sydney, NSW 2052, Australia
| |
Collapse
|
14
|
Kállay M, Horváth RA, Gyevi-Nagy L, Nagy PR. Size-consistent explicitly correlated triple excitation correction. J Chem Phys 2021; 155:034107. [PMID: 34293884 DOI: 10.1063/5.0057426] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
A new approach is proposed to reduce the basis set incompleteness error of the triple excitation correction in explicitly correlated coupled-cluster singles and doubles with perturbative triples calculations. Our method is similar to the intuitive triples correction approach of Knizia et al. [J. Chem. Phys. 130, 054104 (2009)] but, in contrast to the latter, is size-consistent. The new approximation is easy to implement, and its overhead is negligible with respect to the conventional (T) correction. The performance of the approach is assessed for atomization, reaction, and interaction energies as well as for bond lengths and harmonic vibrational frequencies. The advantages of its size consistency are also demonstrated.
Collapse
Affiliation(s)
- Mihály Kállay
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, Budapest H-1521, Hungary
| | - Réka A Horváth
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, Budapest H-1521, Hungary
| | - László Gyevi-Nagy
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, Budapest H-1521, Hungary
| | - Péter R Nagy
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, Budapest H-1521, Hungary
| |
Collapse
|
15
|
Guther K, Cohen AJ, Luo H, Alavi A. Binding curve of the beryllium dimer using similarity-transformed FCIQMC: Spectroscopic accuracy with triple-zeta basis sets. J Chem Phys 2021; 155:011102. [PMID: 34241381 DOI: 10.1063/5.0055575] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We demonstrate how similarity-transformed full configuration interaction quantum Monte Carlo (FCIQMC) based on the transcorrelated Hamiltonian can be applied to make highly accurate predictions for the binding curve of the beryllium dimer, marking the first case study of a molecular system with this method. In this context, the non-Hermitian transcorrelated Hamiltonian, resulting from a similarity transformation with a Jastrow factor, serves the purpose to effectively address dynamic correlation beyond the used basis set and thus allows for obtaining energies close to the complete basis set limit from FCIQMC already with moderate basis sets and computational effort. Building on results from other explicitly correlated methods, we discuss the role of the Jastrow factor and its functional form, as well as potential sources for size consistency errors, and arrive at Jastrow forms that allow for high accuracy calculations of the vibrational spectrum of the beryllium dimer.
Collapse
Affiliation(s)
- Kai Guther
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Federal Republic of Germany
| | - Aron J Cohen
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Federal Republic of Germany
| | - Hongjun Luo
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Federal Republic of Germany
| | - Ali Alavi
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Federal Republic of Germany
| |
Collapse
|
16
|
Doran AE, Hirata S. Stochastic evaluation of fourth-order many-body perturbation energies. J Chem Phys 2021; 154:134114. [PMID: 33832241 DOI: 10.1063/5.0047798] [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/14/2022] Open
Abstract
A scalable, stochastic algorithm evaluating the fourth-order many-body perturbation (MP4) correction to energy is proposed. Three hundred Goldstone diagrams representing the MP4 correction are computer generated and then converted into algebraic formulas expressed in terms of Green's functions in real space and imaginary time. They are evaluated by the direct (i.e., non-Markov, non-Metropolis) Monte Carlo (MC) integration accelerated by the redundant-walker and control-variate algorithms. The resulting MC-MP4 method is efficiently parallelized and is shown to display O(n5.3) size-dependence of cost, which is nearly two ranks lower than the O(n7) dependence of the deterministic MP4 algorithm. It evaluates the MP4/aug-cc-pVDZ energy for benzene, naphthalene, phenanthrene, and corannulene with the statistical uncertainty of 10 mEh (1.1% of the total basis-set correlation energy), 38 mEh (2.6%), 110 mEh (5.5%), and 280 mEh (9.0%), respectively, after about 109 MC steps.
Collapse
Affiliation(s)
- Alexander E Doran
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - So Hirata
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| |
Collapse
|
17
|
Förster A, Visscher L. Low-Order Scaling G0W0 by Pair Atomic Density Fitting. J Chem Theory Comput 2020; 16:7381-7399. [PMID: 33174743 PMCID: PMC7726916 DOI: 10.1021/acs.jctc.0c00693] [Citation(s) in RCA: 33] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2020] [Indexed: 12/18/2022]
Abstract
We derive a low-scaling G0W0 algorithm for molecules using pair atomic density fitting (PADF) and an imaginary time representation of the Green's function and describe its implementation in the Slater type orbital (STO)-based Amsterdam density functional (ADF) electronic structure code. We demonstrate the scalability of our algorithm on a series of water clusters with up to 432 atoms and 7776 basis functions and observe asymptotic quadratic scaling with realistic threshold qualities controlling distance effects and basis sets of triple-ζ (TZ) plus double polarization quality. Also owing to a very small prefactor, a G0W0 calculation for the largest of these clusters takes only 240 CPU hours with these settings. We assess the accuracy of our algorithm for HOMO and LUMO energies in the GW100 database. With errors of 0.24 eV for HOMO energies on the quadruple-ζ level, our implementation is less accurate than canonical all-electron implementations using the larger def2-QZVP GTO-type basis set. Apart from basis set errors, this is related to the well-known shortcomings of the GW space-time method using analytical continuation techniques as well as to numerical issues of the PADF approach of accurately representing diffuse atomic orbital (AO) products. We speculate that these difficulties might be overcome by using optimized auxiliary fit sets with more diffuse functions of higher angular momenta. Despite these shortcomings, for subsets of medium and large molecules from the GW5000 database, the error of our approach using basis sets of TZ and augmented double-ζ (DZ) quality is decreasing with system size. On the augmented DZ level, we reproduce canonical, complete basis set limit extrapolated reference values with an accuracy of 80 meV on average for a set of 20 large organic molecules. We anticipate our algorithm, in its current form, to be very useful in the study of single-particle properties of large organic systems such as chromophores and acceptor molecules.
Collapse
Affiliation(s)
- Arno Förster
- Theoretical Chemistry, Vrije Universiteit, De Boelelaan 1083, NL-1081 HV Amsterdam, The Netherlands
| | - Lucas Visscher
- Theoretical Chemistry, Vrije Universiteit, De Boelelaan 1083, NL-1081 HV Amsterdam, The Netherlands
| |
Collapse
|
18
|
Xu E, Uejima M, Ten-No SL. Towards Near-Exact Solutions of Molecular Electronic Structure: Full Coupled-Cluster Reduction with a Second-Order Perturbative Correction. J Phys Chem Lett 2020; 11:9775-9780. [PMID: 33146534 DOI: 10.1021/acs.jpclett.0c03084] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
We introduce a new augmented adaptation of the recently developed full coupled-cluster reduction (FCCR) with a second-order perturbative correction, abbreviated as FCCR(2). FCCR is a selected coupled-cluster expansion aimed at optimally reducing the excitation manifold and commutator expansions for high-rank excitations for obtaining accurate solutions of the electronic Schödinger equation in a size-extensive manner. The present FCCR(2) enables estimating the residual correlation of FCCR by the second-order perturbative correction E(2) from the complementary space of the FCCR projection manifold. The linear relationship between E(2) and the energy of FCCR(2) allows accurate estimates of near-exact energies for a wide variety of molecules with strong electron correlation. The potential of the method is demonstrated using challenging cases, the ground-state electronic energy of the benzene molecule in equilibrium and stretched geometries, and the isomerization energy of the transition metal complex [Cu(NH3)]2O22+.
Collapse
Affiliation(s)
- Enhua Xu
- Graduate School of Science, Technology, and Innovation, Kobe University, Kobe 657-8501, Japan
| | - Motoyuki Uejima
- Graduate School of Science, Technology, and Innovation, Kobe University, Kobe 657-8501, Japan
| | - Seiichiro L Ten-No
- Graduate School of Science, Technology, and Innovation, Kobe University, Kobe 657-8501, Japan
- Graduate School of System Informatics, Kobe University, Kobe 657-8501, Japan
| |
Collapse
|
19
|
Abstract
We introduce the transcorrelated Density Matrix Renormalization Group (tcDMRG) theory for the efficient approximation of the energy for strongly correlated systems. tcDMRG encodes the wave function as a product of a fixed Jastrow or Gutzwiller correlator and a matrix product state. The latter is optimized by applying the imaginary-time variant of time-dependent (TD) DMRG to the non-Hermitian transcorrelated Hamiltonian. We demonstrate the efficiency of tcDMRG with the example of the two-dimensional Fermi-Hubbard Hamiltonian, a notoriously difficult target for the DMRG algorithm, for different sizes, occupation numbers, and interaction strengths. We demonstrate fast energy convergence of tcDMRG, which indicates that tcDMRG could increase the efficiency of standard DMRG beyond quasi-monodimensional systems and provides a generally powerful approach toward the dynamic correlation problem of DMRG.
Collapse
Affiliation(s)
- Alberto Baiardi
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Markus Reiher
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
20
|
Bauer B, Bravyi S, Motta M, Chan GKL. Quantum Algorithms for Quantum Chemistry and Quantum Materials Science. Chem Rev 2020; 120:12685-12717. [DOI: 10.1021/acs.chemrev.9b00829] [Citation(s) in RCA: 134] [Impact Index Per Article: 33.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Affiliation(s)
- Bela Bauer
- Microsoft Quantum, Station Q, University of California
, Santa Barbara, California 93106, United States
| | - Sergey Bravyi
- IBM Quantum, IBM T. J. Watson Research Center
, Yorktown Heights, New York 10598, United States
| | - Mario Motta
- IBM Quantum, IBM Research Almaden
, San Jose, California 95120, United States
| | - Garnet Kin-Lic Chan
- Division of Chemistry and Chemical Engineering, California Institute of Technology
, Pasadena, California 91125, United States
| |
Collapse
|
21
|
Doran AE, Hirata S. Convergence acceleration of Monte Carlo many-body perturbation methods by direct sampling. J Chem Phys 2020; 153:104112. [DOI: 10.1063/5.0020583] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Alexander E. Doran
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - So Hirata
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| |
Collapse
|
22
|
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
|
23
|
Doran AE, Hirata S. Convergence acceleration of Monte Carlo many-body perturbation methods by using many control variates. J Chem Phys 2020; 153:094108. [DOI: 10.1063/5.0020584] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Affiliation(s)
- Alexander E. Doran
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - So Hirata
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| |
Collapse
|
24
|
Muolo A, Baiardi A, Feldmann R, Reiher M. Nuclear-electronic all-particle density matrix renormalization group. J Chem Phys 2020; 152:204103. [PMID: 32486651 DOI: 10.1063/5.0007166] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
We introduce the Nuclear-Electronic All-Particle Density Matrix Renormalization Group (NEAP-DMRG) method for solving the time-independent Schrödinger equation simultaneously for electrons and other quantum species. In contrast to the already existing multicomponent approaches, in this work, we construct from the outset a multi-reference trial wave function with stochastically optimized non-orthogonal Gaussian orbitals. By iterative refining of the Gaussians' positions and widths, we obtain a compact multi-reference expansion for the multicomponent wave function. We extend the DMRG algorithm to multicomponent wave functions to take into account inter- and intra-species correlation effects. The efficient parameterization of the total wave function as a matrix product state allows NEAP-DMRG to accurately approximate the full configuration interaction energies of molecular systems with more than three nuclei and 12 particles in total, which is currently a major challenge for other multicomponent approaches. We present the NEAP-DMRG results for two few-body systems, i.e., H2 and H3 +, and one larger system, namely, BH3.
Collapse
Affiliation(s)
- Andrea Muolo
- ETH Zürich, Laboratory of Physical Chemistry, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Alberto Baiardi
- ETH Zürich, Laboratory of Physical Chemistry, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Robin Feldmann
- ETH Zürich, Laboratory of Physical Chemistry, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Markus Reiher
- ETH Zürich, Laboratory of Physical Chemistry, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
25
|
Förster A, Franchini M, van Lenthe E, Visscher L. A Quadratic Pair Atomic Resolution of the Identity Based SOS-AO-MP2 Algorithm Using Slater Type Orbitals. J Chem Theory Comput 2020; 16:875-891. [PMID: 31930915 PMCID: PMC7027358 DOI: 10.1021/acs.jctc.9b00854] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2019] [Indexed: 01/04/2023]
Abstract
We report a production level implementation of pair atomic resolution of the identity (PARI) based second-order Møller-Plesset perturbation theory (MP2) in the Slater type orbital (STO) based Amsterdam Density Functional (ADF) code. As demonstrated by systematic benchmarks, dimerization and isomerization energies obtained with our code using STO basis sets of triple-ζ-quality show mean absolute deviations from Gaussian type orbital, canonical, basis set limit extrapolated, global density fitting (DF)-MP2 results of less than 1 kcal/mol. Furthermore, we introduce a quadratic scaling atomic orbital based spin-opposite-scaled (SOS)-MP2 approach with a very small prefactor. Due to a worst-case scaling of [Formula: see text], our implementation is very fast already for small systems and shows an exceptionally early crossover to canonical SOS-PARI-MP2. We report computational wall time results for linear as well as for realistic three-dimensional molecules and show that triple-ζ quality calculations on molecules of several hundreds of atoms are only a matter of a few hours on a single compute node, the bottleneck of the computations being the SCF rather than the post-SCF energy correction.
Collapse
Affiliation(s)
- Arno Förster
- Theoretical Chemistry, Vrije
Universiteit, De Boelelaan 1083, NL-1081 HV Amsterdam, The
Netherlands
| | - Mirko Franchini
- Theoretical Chemistry, Vrije
Universiteit, De Boelelaan 1083, NL-1081 HV Amsterdam, The
Netherlands
- Scientific Computing & Modelling
NV, De Boelelaan 1083, NL-1081 HV Amsterdam, The
Netherlands
| | - Erik van Lenthe
- Scientific Computing & Modelling
NV, De Boelelaan 1083, NL-1081 HV Amsterdam, The
Netherlands
| | - Lucas Visscher
- Theoretical Chemistry, Vrije
Universiteit, De Boelelaan 1083, NL-1081 HV Amsterdam, The
Netherlands
| |
Collapse
|
26
|
Shaw RA, Hill JG. A Linear-Scaling Method for Noncovalent Interactions: An Efficient Combination of Absolutely Localized Molecular Orbitals and a Local Random Phase Approximation Approach. J Chem Theory Comput 2019; 15:5352-5369. [PMID: 31465215 DOI: 10.1021/acs.jctc.9b00615] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A novel method for the accurate and efficient calculation of interaction energies in weakly bound complexes composed of a large number of molecules is presented. The new ALMO+RPAd method circumvents the prohibitive scaling of coupled cluster singles and doubles while still providing similar accuracy across a diverse range of weakly bound chemical systems. Linear-scaling procedures for the Fock build are given utilizing absolutely localized molecular orbitals (ALMOs), resulting in the a priori exclusion of basis set superposition errors. A bespoke data structure and algorithm using density fitting are described, leading to linear scaling for the storage and computation of the two-electron integrals. Electron correlation is included through a new, linear-scaling pairwise local random phase approximation approach, including exchange interactions, and decomposed into purely dispersive excitations (RPAxd). Collectively, these allow meaningful decomposition of the interaction energy into physically distinct contributions: electrostatic, polarization, charge transfer, and dispersion. Comparison with symmetry-adapted perturbation theory shows good qualitative agreement. Tests on various dimers and the S66 benchmark set demonstrate results within 0.5 kcal mol-1 of coupled cluster singles and doubles results. On a large cluster of water molecules, we achieve calculations involving over 3500 orbital and 12,000 auxiliary basis functions in under 10 min on a single CPU core.
Collapse
Affiliation(s)
- Robert A Shaw
- Department of Chemistry , University of Sheffield , Sheffield S3 7HF , United Kingdom
| | - J Grant Hill
- Department of Chemistry , University of Sheffield , Sheffield S3 7HF , United Kingdom
| |
Collapse
|
27
|
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
|
28
|
Nudejima T, Ikabata Y, Seino J, Yoshikawa T, Nakai H. Machine-learned electron correlation model based on correlation energy density at complete basis set limit. J Chem Phys 2019; 151:024104. [DOI: 10.1063/1.5100165] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Affiliation(s)
- Takuro Nudejima
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| | - Yasuhiro Ikabata
- Waseda Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| | - Junji Seino
- Waseda Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
- PRESTO, Japan Science and Technology Agency, 4-1-8 Honcho, Kawaguchi, Saitama 332-0012, Japan
| | - Takeshi Yoshikawa
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| | - Hiromi Nakai
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
- Waseda Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
- Elements Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University, Katsura, Kyoto 615-8520, Japan
| |
Collapse
|
29
|
Teke NK, Pavošević F, Peng C, Valeev EF. Explicitly correlated renormalized second-order Green’s function for accurate ionization potentials of closed-shell molecules. J Chem Phys 2019; 150:214103. [DOI: 10.1063/1.5090983] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Nakul K. Teke
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Fabijan Pavošević
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Chong Peng
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Edward F. Valeev
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| |
Collapse
|
30
|
Höfener S, Schieschke N, Klopper W, Köhn A. The extended explicitly-correlated second-order approximate coupled-cluster singles and doubles ansatz suitable for response theory. J Chem Phys 2019; 150:184110. [DOI: 10.1063/1.5094434] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Affiliation(s)
- Sebastian Höfener
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), P.O. Box 6980, D-76049 Karlsruhe, Germany
| | - Nils Schieschke
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), P.O. Box 6980, D-76049 Karlsruhe, Germany
| | - Wim Klopper
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), P.O. Box 6980, D-76049 Karlsruhe, Germany
| | - Andreas Köhn
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| |
Collapse
|
31
|
Almásy L, Bende A. Intermolecular Interaction in Methylene Halide ( CH2F2, CH2Cl2, CH2Br2 and CH2I2) Dimers. Molecules 2019; 24:E1810. [PMID: 31083344 PMCID: PMC6539938 DOI: 10.3390/molecules24091810] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2019] [Revised: 04/27/2019] [Accepted: 05/01/2019] [Indexed: 02/06/2023] Open
Abstract
The intermolecular interaction in difluoromethane, dichloromethane, dibromomethane, and diiodomethane dimers has been investigated using high level quantum chemical methods. The potential energy curve of intermolecular interaction along the C⋯C bond distance obtained using the coupled-cluster theory with singles, doubles, and perturbative triples excitations CCSD(T) were compared with values given by the same method, but applying the local (LCCSD(T)) and the explicitly correlated (CCSD(T)-F12) approximations. The accuracy of other theoretical methods-Hartree-Fock (HF), second order Møller-Plesset perturbation (MP2), and dispersion corrected DFT theory-were also presented. In the case of MP2 level, the canonical and the local-correlation cases combined with the density-fitting technique (DF-LMP2)theories were considered, while for the dispersion-corrected DFT, the empirically-corrected BLYP-D and the M06-2Xexchange-correlation functionals were applied. In all cases, the aug-cc-pVTZ basis set was used, and the results were corrected for the basis set superposition error (BSSE) using the counterpoise method. For each molecular system, several dimer geometries were found, and their mutual orientations were compared with the nearest neighbor orientations obtained in recent neutron scattering studies. The nature of the intermolecular interaction energy was discussed.
Collapse
MESH Headings
- Dimerization
- Hydrocarbons, Brominated/analysis
- Hydrocarbons, Brominated/chemistry
- Hydrocarbons, Chlorinated/analysis
- Hydrocarbons, Chlorinated/chemistry
- Hydrocarbons, Fluorinated/analysis
- Hydrocarbons, Fluorinated/chemistry
- Hydrocarbons, Iodinated/analysis
- Hydrocarbons, Iodinated/chemistry
- Inorganic Chemicals/analysis
- Inorganic Chemicals/chemistry
- Models, Chemical
- Models, Molecular
- Quantum Theory
Collapse
Affiliation(s)
- László Almásy
- State Key Laboratory of Environment-friendly Energy Materials, Southwest University of Science and Technology, Mianyang 621010, China.
- Institute for Solid State Physics and Optics, Wigner Research Centre for Physics, Konkoly Thege út 29-33, 1121 Budapest, Hungary.
| | - Attila Bende
- Molecular and Biomolecular Physics Department, National Institute for Research and Development of Isotopic and Molecular Technologies, Donat Street, No. 67-103, Ro-400293 Cluj-Napoca, Romania.
| |
Collapse
|
32
|
Feldt M, Phung QM, Pierloot K, Mata RA, Harvey JN. Limits of Coupled-Cluster Calculations for Non-Heme Iron Complexes. J Chem Theory Comput 2019; 15:922-937. [PMID: 30605326 DOI: 10.1021/acs.jctc.8b00963] [Citation(s) in RCA: 36] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
In a large variety of studies, the coupled-cluster method with singles, doubles, and perturbative triples (CCSD(T)) is used as a reference for benchmarking the performance of density functional theory (DFT) functionals. In the case of open-shell species, this theory can be applied in different forms depending on the restricted or unrestricted treatment of spin. In this study, we show that these different approaches can produce results which deviate by ∼5 kcal/mol for different species on the potential energy surfaces. This was demonstrated for a simple model of the C-H activation carried out by non-heme iron enzymes. Assessing the limits of CCSD(T) prior to its use as a general benchmark tool is warranted. This was done using higher-order coupled-cluster calculations as well as multiconfigurational second-order perturbation theory (CASPT2), since iron-oxo species present some multireference character. Furthermore, we tested two different implementations of the local coupled-cluster method and compared them to the CCSD(T) results, showing that even though these novel approaches are promising, without further developments they appear not to be suitable for describing two-state reactivity of the system investigated in the current study. Additionally, we implemented and assessed the performance of the hotspot approach for the local unrestricted CCSD(T) scheme which aims at reducing the pair error for systems containing transition metals.
Collapse
Affiliation(s)
- Milica Feldt
- Department of Chemistry , KU Leuven , Celestijnenlaan 200f, Box 2404, 3001 Leuven , Belgium
| | - Quan Manh Phung
- Department of Chemistry , KU Leuven , Celestijnenlaan 200f, Box 2404, 3001 Leuven , Belgium
| | - Kristine Pierloot
- Department of Chemistry , KU Leuven , Celestijnenlaan 200f, Box 2404, 3001 Leuven , Belgium
| | - Ricardo A Mata
- Institut für Physikalische Chemie , Universität Göttingen , Tammannstrasse 6 , D-37077 Göttingen , Germany
| | - Jeremy N Harvey
- Department of Chemistry , KU Leuven , Celestijnenlaan 200f, Box 2404, 3001 Leuven , Belgium
| |
Collapse
|
33
|
Al-Hamdani YS, Tkatchenko A. Understanding non-covalent interactions in larger molecular complexes from first principles. J Chem Phys 2019; 150:010901. [PMID: 30621423 PMCID: PMC6910608 DOI: 10.1063/1.5075487] [Citation(s) in RCA: 38] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2018] [Accepted: 12/05/2018] [Indexed: 01/02/2023] Open
Abstract
Non-covalent interactions pervade all matter and play a fundamental role in layered materials, biological systems, and large molecular complexes. Despite this, our accumulated understanding of non-covalent interactions to date has been mainly developed in the tens-of-atoms molecular regime. This falls considerably short of the scales at which we would like to understand energy trends, structural properties, and temperature dependencies in materials where non-covalent interactions have an appreciable role. However, as more reference information is obtained beyond moderately sized molecular systems, our understanding is improving and we stand to gain pertinent insights by tackling more complex systems, such as supramolecular complexes, molecular crystals, and other soft materials. In addition, accurate reference information is needed to provide the drive for extending the predictive power of more efficient workhorse methods, such as density functional approximations that also approximate van der Waals dispersion interactions. In this perspective, we discuss the first-principles approaches that have been used to obtain reference interaction energies for beyond modestly sized molecular complexes. The methods include quantum Monte Carlo, symmetry-adapted perturbation theory, non-canonical coupled cluster theory, and approaches based on the random-phase approximation. By considering the approximations that underpin each method, the most accurate theoretical references for supramolecular complexes and molecular crystals to date are ascertained. With these, we also assess a handful of widely used exchange-correlation functionals in density functional theory. The discussion culminates in a framework for putting into perspective the accuracy of high-level wavefunction-based methods and identifying future challenges.
Collapse
Affiliation(s)
- Yasmine S Al-Hamdani
- Physics and Materials Science Research Unit, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - Alexandre Tkatchenko
- Physics and Materials Science Research Unit, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| |
Collapse
|
34
|
Bischoff FA. Computing accurate molecular properties in real space using multiresolution analysis. ADVANCES IN QUANTUM CHEMISTRY 2019. [DOI: 10.1016/bs.aiq.2019.04.003] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
35
|
Johnson CM, Doran AE, Ten-no SL, Hirata S. Monte Carlo explicitly correlated many-body Green’s function theory. J Chem Phys 2018; 149:174112. [DOI: 10.1063/1.5054610] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Cole M. Johnson
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Alexander E. Doran
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Seiichiro L. Ten-no
- Graduate School of Science, Technology, and Innovation, Kobe University, Nada-ku, Kobe 657-8501, Japan
- Graduate School of System Informatics, Kobe University, Nada-ku, Kobe 657-8501, Japan
| | - So Hirata
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| |
Collapse
|
36
|
Li MW, Zimmerman PM. Stepwise basis set selection. J Comput Chem 2018; 39:2153-2162. [DOI: 10.1002/jcc.25363] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2018] [Accepted: 05/01/2018] [Indexed: 11/08/2022]
Affiliation(s)
- M. W. Li
- Department of ChemistryUniversity of Michigan930 N. University Ave, Ann Arbor MI48109
| | - P. M. Zimmerman
- Department of ChemistryUniversity of Michigan930 N. University Ave, Ann Arbor MI48109
| |
Collapse
|
37
|
Tew DP, Kats D. Relaxing Constrained Amplitudes: Improved F12 Treatments of Orbital Optimization and Core–Valence Correlation Energies. J Chem Theory Comput 2018; 14:5435-5440. [DOI: 10.1021/acs.jctc.8b00792] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- David P. Tew
- Max Planck Institute for Solid State
Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Daniel Kats
- Max Planck Institute for Solid State
Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| |
Collapse
|
38
|
Nakatsuji H, Nakashima H, Kurokawa YI. Solving the Schrödinger equation of atoms and molecules: Chemical-formula theory, free-complement chemical-formula theory, and intermediate variational theory. J Chem Phys 2018; 149:114105. [PMID: 30243277 DOI: 10.1063/1.5040376] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Chemistry is governed by the principle of quantum mechanics as expressed by the Schrödinger equation (SE) and Dirac equation (DE). The exact general theory for solving these fundamental equations is therefore a key for formulating accurately predictive theory in chemical science. The free-complement (FC) theory for solving the SE of atoms and molecules proposed by one of the authors is such a general theory. On the other hand, the working theory most widely used in chemistry is the chemical formula that refers to the molecular structural formula and chemical reaction formula, collectively. There, the central concepts are the local atomic concept, transferability, and from-atoms-to-molecule concept. Since the chemical formula is the most successful working theory in chemistry ever existed, we formulate our FC theory to have the structure reflecting the chemical formula. Our basic postulate is that as far as the SE is the principle of chemistry, its solutions for chemistry should have the structure that can be related to the chemical formulas. So, in this paper, we first formulate a theory that designs the wave function to reflect the structure of the chemical formula. We call this theory chemical formula theory (CFT). In the CFT, we place the valence ground and excited states of each atom at each position of the chemical formula of the molecule and let them interact using their free valences to form the ground and excited states of the molecule. The principle there is the variational principle so that the ground and excited states obtained satisfy the orthogonality and Hamiltonian-orthogonality relations. Then, we formulate the exact FC theory starting from the initial functions produced by the CFT. This FC theory is referred to as free-complement chemical-formula theory (FC-CFT), which is expected to describe efficiently the solution of the SE by the above reason. The FC-CFT wave function is modified from that of CFT. Since this modification is done by the exact SE, its analysis may give some insights to chemists that assist their chemistry. Thus, this theory would be not only exact but also conceptually useful. Furthermore, the intermediate theory between CFT and FC-CFT would also be useful. There, we use only integratable functions and apply the variational principle so that we refer to this theory as FC-CFT-variational (FC-CFT-V). It is an advanced theory of CFT. Since the variational method is straightforward and powerful, we can do extensive chemical studies in a reasonable accuracy. After finishing such studies, if we still need an exact level of solutions, we add the remaining functions of the FC-CFT and perform the exact calculations. Furthermore, when we deal with large and even giant molecules, the inter-exchange (iExg) theory for the antisymmetry rule introduced previously leads to a large simplification. There, the inter-exchanges between distant electron pairs fade away so that only Coulombic interactions survive. Further in giant systems, even an electrostatic description becomes possible. Then, the FC-CFT for exactly solving the SE would behave essentially to order N for large and giant molecular systems, though the pre-factor should be very large and must be minimized.
Collapse
Affiliation(s)
- Hiroshi Nakatsuji
- Quantum Chemistry Research Institute, Kyoto Technoscience Center 16, 14 Yoshida Kawaramachi, Sakyo-ku, Kyoto 606-8305, Japan
| | - Hiroyuki Nakashima
- Quantum Chemistry Research Institute, Kyoto Technoscience Center 16, 14 Yoshida Kawaramachi, Sakyo-ku, Kyoto 606-8305, Japan
| | - Yusaku I Kurokawa
- Quantum Chemistry Research Institute, Kyoto Technoscience Center 16, 14 Yoshida Kawaramachi, Sakyo-ku, Kyoto 606-8305, Japan
| |
Collapse
|
39
|
Wiles TC, Manby FR. Wavefunction-like Correlation Model for Use in Hybrid Density Functionals. J Chem Theory Comput 2018; 14:4590-4599. [PMID: 30080967 DOI: 10.1021/acs.jctc.8b00337] [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/30/2022]
Abstract
We present Unsöld-W12 (UW12), an approximation to the correlation energy of molecules that is an explicit functional of the single-particle reduced-density matrix. The approximation resembles one part of modern explicitly correlated second-order Møller-Plesset (MP2) theory and is intended as an alternative to MP2 in double-hybrid exchange-correlation functionals. Orbital optimization with UW12 is straightforward, and the UW12 energy is evaluated without a double summation over unoccupied orbitals, leading to a faster basis-set convergence than is seen in double-hybrid functionals. We suggest a one-parameter hybrid exchange-correlation functional XCH-BLYP-UW12. XCH-BLYP-UW12 is similar to double-hybrid functionals, but contains UW12 correlation instead of MP2 correlation. We find that XCH-BLYP-UW12 is more accurate than the existing double-hybrid functional B2-PLYP for small-molecule main-group reaction barrier heights and has roughly the same accuracy as the existing hybrid functional B3LYP for atomization energies.
Collapse
Affiliation(s)
- Timothy C Wiles
- Centre for Computational Chemistry, School of Chemistry , University of Bristol , Bristol BS8 1TS , United Kingdom
| | - Frederick R Manby
- Centre for Computational Chemistry, School of Chemistry , University of Bristol , Bristol BS8 1TS , United Kingdom
| |
Collapse
|
40
|
Bayne MG, Scher JA, Ellis BH, Chakraborty A. Linked-Cluster Formulation of Electron–Hole Interaction Kernel in Real-Space Representation without Using Unoccupied States. J Chem Theory Comput 2018; 14:3656-3666. [DOI: 10.1021/acs.jctc.8b00123] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Michael G. Bayne
- Department of Chemistry, Syracuse University, Syracuse, New York 13244 United States
| | - Jeremy A. Scher
- Department of Chemistry, Syracuse University, Syracuse, New York 13244 United States
| | - Benjamin H. Ellis
- Department of Chemistry, Syracuse University, Syracuse, New York 13244 United States
| | - Arindam Chakraborty
- Department of Chemistry, Syracuse University, Syracuse, New York 13244 United States
| |
Collapse
|
41
|
Affiliation(s)
- Robert A. Shaw
- Department of Chemistry, University of Sheffield, Sheffield, UK
| | - J. Grant Hill
- Department of Chemistry, University of Sheffield, Sheffield, UK
| |
Collapse
|
42
|
Luo H, Alavi A. Combining the Transcorrelated Method with Full Configuration Interaction Quantum Monte Carlo: Application to the Homogeneous Electron Gas. J Chem Theory Comput 2018; 14:1403-1411. [DOI: 10.1021/acs.jctc.7b01257] [Citation(s) in RCA: 36] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Hongjun Luo
- Max-Planck-Institut for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Ali Alavi
- Max-Planck-Institut for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
43
|
Temelso B, Morrison EF, Speer DL, Cao BC, Appiah-Padi N, Kim G, Shields GC. Effect of Mixing Ammonia and Alkylamines on Sulfate Aerosol Formation. J Phys Chem A 2018; 122:1612-1622. [DOI: 10.1021/acs.jpca.7b11236] [Citation(s) in RCA: 42] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Berhane Temelso
- Provost’s
Office and Department of Chemistry, Furman University, Greenville, South Carolina 29613, United States
- Dean’s
Office, College of Arts and Sciences, and Department of Chemistry, Bucknell University, Lewisburg, Pennsylvania 17837, United States
| | - Elizabeth F. Morrison
- Dean’s
Office, College of Arts and Sciences, and Department of Chemistry, Bucknell University, Lewisburg, Pennsylvania 17837, United States
| | - David L. Speer
- Dean’s
Office, College of Arts and Sciences, and Department of Chemistry, Bucknell University, Lewisburg, Pennsylvania 17837, United States
| | - Bobby C. Cao
- Dean’s
Office, College of Arts and Sciences, and Department of Chemistry, Bucknell University, Lewisburg, Pennsylvania 17837, United States
| | - Nana Appiah-Padi
- Dean’s
Office, College of Arts and Sciences, and Department of Chemistry, Bucknell University, Lewisburg, Pennsylvania 17837, United States
| | - Grace Kim
- Dean’s
Office, College of Arts and Sciences, and Department of Chemistry, Bucknell University, Lewisburg, Pennsylvania 17837, United States
| | - George C. Shields
- Provost’s
Office and Department of Chemistry, Furman University, Greenville, South Carolina 29613, United States
- Dean’s
Office, College of Arts and Sciences, and Department of Chemistry, Bucknell University, Lewisburg, Pennsylvania 17837, United States
| |
Collapse
|
44
|
Temelso B, Klein KL, Mabey JW, Pérez C, Pate BH, Kisiel Z, Shields GC. Exploring the Rich Potential Energy Surface of (H2O)11 and Its Physical Implications. J Chem Theory Comput 2018; 14:1141-1153. [DOI: 10.1021/acs.jctc.7b00938] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Affiliation(s)
- Berhane Temelso
- Provost’s
Office and Department of Chemistry, Furman University, Greenville, South Carolina 29613, United States
- Dean’s
Office, College of Arts and Sciences, and Department of Chemistry, Bucknell University, Lewisburg, Pennsylvania 17837, United States
| | - Katurah L. Klein
- Dean’s
Office, College of Arts and Sciences, and Department of Chemistry, Bucknell University, Lewisburg, Pennsylvania 17837, United States
| | - Joel W. Mabey
- Dean’s
Office, College of Arts and Sciences, and Department of Chemistry, Bucknell University, Lewisburg, Pennsylvania 17837, United States
| | - Cristóbal Pérez
- Department
of Chemistry, University of Virginia, McCormick Road, Charlottesville, Virginia 22904-4319, United States
- Max Planck Institute for the Structure and Dynamics of Matter, Luruper Chausse 149, D-22761 Hamburg, Germany
| | - Brooks H. Pate
- Department
of Chemistry, University of Virginia, McCormick Road, Charlottesville, Virginia 22904-4319, United States
| | - Zbigniew Kisiel
- Institute
of Physics, Polish Academy of Sciences, Al. Lotników 32/46, 02-668 Warszawa, Poland
| | - George C. Shields
- Provost’s
Office and Department of Chemistry, Furman University, Greenville, South Carolina 29613, United States
- Dean’s
Office, College of Arts and Sciences, and Department of Chemistry, Bucknell University, Lewisburg, Pennsylvania 17837, United States
| |
Collapse
|
45
|
Segarra-Martí J, Garavelli M, Aquilante F. Converging many-body correlation energies by means of sequence extrapolation. J Chem Phys 2018; 148:034107. [DOI: 10.1063/1.5000783] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Affiliation(s)
- J. Segarra-Martí
- Laboratoire de Chimie UMR 5182, ENS de Lyon, 46 Allée d’Italie, 69364 Lyon Cedex 07, France
| | - M. Garavelli
- Dipartimento di Chimica Industriale “Toso Montanari,” Università di Bologna, Viale del Risorgimento 4, 40136 Bologna, Italy
| | - F. Aquilante
- Dipartimento di Chimica “Giacomo Ciamician,” Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy
| |
Collapse
|
46
|
Harris FE. Atomic Electronic Structure Computations With Hylleraas-CI Wave Functions. ADVANCES IN QUANTUM CHEMISTRY 2018. [DOI: 10.1016/bs.aiq.2017.04.001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
47
|
Barca GM, Loos PF. Recurrence Relations for Four-Electron Integrals Over Gaussian Basis Functions. ADVANCES IN QUANTUM CHEMISTRY 2018. [DOI: 10.1016/bs.aiq.2017.03.004] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
48
|
Sylvetsky N, Kesharwani MK, Martin JML. The aug-cc-pVnZ-F12 basis set family: Correlation consistent basis sets for explicitly correlated benchmark calculations on anions and noncovalent complexes. J Chem Phys 2017; 147:134106. [DOI: 10.1063/1.4998332] [Citation(s) in RCA: 41] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Nitai Sylvetsky
- Department of Organic Chemistry, Weizmann Institute of Science, 76100 Reḥovot, Israel
| | - Manoj K. Kesharwani
- Department of Organic Chemistry, Weizmann Institute of Science, 76100 Reḥovot, Israel
| | - Jan M. L. Martin
- Department of Organic Chemistry, Weizmann Institute of Science, 76100 Reḥovot, Israel
| |
Collapse
|
49
|
Harris FE. Integrals arising in atomic Hylleraas methods. Mol Phys 2017. [DOI: 10.1080/00268976.2016.1277593] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- Frank E. Harris
- Department of Physics, University of Utah, Salt Lake City, UT, USA
- Quantum Theory Project, University of Florida, Gainesville, FL, USA
| |
Collapse
|
50
|
|