1
|
Lechner MH, Papadopoulos A, Sivalingam K, Auer AA, Koslowski A, Becker U, Wennmohs F, Neese F. Code generation in ORCA: progress, efficiency and tight integration. Phys Chem Chem Phys 2024; 26:15205-15220. [PMID: 38767596 DOI: 10.1039/d4cp00444b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/22/2024]
Abstract
An improved version of ORCA's automated generator environment (ORCA-AGE II) is presented. The algorithmic improvements and the move to C++ as the programming language lead to a performance gain of up to two orders of magnitude compared to the previously developed PYTHON toolchain. Additionally, the restructured modular design allows for far more complex code engines to be implemented readily. Importantly, we have realised an extremely tight integration with the ORCA host program. This allows for a workflow in which only the wavefunction Ansatz is part of the source code repository while all actual high-level code is generated automatically, inserted at the appropriate place in the host program before it is compiled and linked together with the hand written code parts. This construction ensures longevity and uniform code quality. Furthermore the new developments allow ORCA-AGE II to generate parallelised production-level code for highly complex theories, such as fully internally contracted multireference coupled-cluster theory (fic-MRCC) with an enormous number of contributing tensor contractions. We also discuss the automated implementation of nuclear gradients for arbitrary theories. All these improvements enable the implementation of theories that are too complex for the human mind and also reduce development times by orders of magnitude. We hope that this work enables researchers to concentrate on the intellectual content of the theories they develop rather than be concerned with technical details of the implementation.
Collapse
Affiliation(s)
- Marvin H Lechner
- Department of Molecular Theory and Spectroscopy, Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany.
| | - Anastasios Papadopoulos
- Department of Molecular Theory and Spectroscopy, Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany.
| | - Kantharuban Sivalingam
- Department of Molecular Theory and Spectroscopy, Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany.
| | - Alexander A Auer
- Department of Molecular Theory and Spectroscopy, Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany.
| | - Axel Koslowski
- Department of Molecular Theory and Spectroscopy, Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany.
| | - Ute Becker
- Department of Molecular Theory and Spectroscopy, Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany.
| | - Frank Wennmohs
- Department of Molecular Theory and Spectroscopy, Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany.
| | - Frank Neese
- Department of Molecular Theory and Spectroscopy, Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany.
| |
Collapse
|
2
|
Evangelista FA. Automatic derivation of many-body theories based on general Fermi vacua. J Chem Phys 2022; 157:064111. [PMID: 35963725 DOI: 10.1063/5.0097858] [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
This paper describes Wick&d, an implementation of the algebra of second-quantized operators normal ordered with respect to general correlated references and the corresponding Wick theorem [D. Mukherjee, Chem. Phys. Lett. 274, 561 (1997) and W. Kutzelnigg and D. Mukherjee, J. Chem. Phys. 107, 432 (1997)]. Wick&d employs a compact representation of operators and a backtracking algorithm to efficiently evaluate Wick contractions. Since Wick&d can handle both fully and partially contracted terms, it can be applied to both projective and Fock-space many-body formalisms. To demonstrate the usefulness of Wick&d, we use it to evaluate the single-reference coupled cluster equations up to octuple excitations and report an automated derivation and implementation of the second-order driven similarity renormalization group multireference perturbation theory.
Collapse
Affiliation(s)
- Francesco A Evangelista
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| |
Collapse
|
3
|
Christlmaier EM, Kats D, Alavi A, Usvyat D. Full Configuration Interaction Quantum Monte Carlo treatment of fragments embedded in a periodic mean field. J Chem Phys 2022; 156:154107. [DOI: 10.1063/5.0084040] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present an embedded fragment approach for high-level quantum chemical calculations on local features in periodic systems. The fragment is defined as a set of localized orbitals (occupied and virtual) corresponding to a converged periodic Hartree-Fock solution. These orbitals serve as the basis for the in-fragment post-Hartree Fock treatment. The embedding field for the fragment, consisting of the Coulomb and exchange potential from the rest of the crystal, is included in the fragment's one-electron Hamiltonian. As an application of the embedded fragment approach we investigate the performanceof full configuration interaction quantum Monte Carlo (FCIQMC) with the adaptive shift. As the orbital choice we use the natural orbitals from the distinguishable cluster method with singles and doubles. FCIQMC is a stochastic approximation to the full CI method and can be routinely applied to much larger active spaces than the latter. This makes this method especially attractive in the context of open shell defects in crystals, where fragments of adequate size can be ratherlarge. As a test case we consider dissociation of a fluorine atom from a fluorographane surface. This process poses a challenge for high-level electronic structure models as both the static and dynamic correlations are essential here. Furthermore the active space for an adequate fragment (32 electrons in 173 orbitals) is already quite large even for FCIQMC. Despite this, FCIQMC delivers accurate dissociation and total energies.
Collapse
Affiliation(s)
| | - Daniel Kats
- Max-Planck-Institute for Solid State Research, Germany
| | - Ali Alavi
- Max-Planck-Institute for Solid State Research, Germany
| | - Denis Usvyat
- Institute of Chemistry, Humboldt University of Berlin, Germany
| |
Collapse
|
4
|
Hansen AS, Aurbakken E, Pedersen TB. Smooth potential-energy surfaces in fragmentation-based local correlation methods for periodic systems. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1896046] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
Affiliation(s)
- A. S. Hansen
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, Oslo, Norway
| | - E. Aurbakken
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, Oslo, Norway
| | - T. B. Pedersen
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, Oslo, Norway
| |
Collapse
|
5
|
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
|
6
|
Calvin JA, Peng C, Rishi V, Kumar A, Valeev EF. Many-Body Quantum Chemistry on Massively Parallel Computers. Chem Rev 2020; 121:1203-1231. [DOI: 10.1021/acs.chemrev.0c00006] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Affiliation(s)
- Justus A. Calvin
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Chong Peng
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Varun Rishi
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Ashutosh Kumar
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Edward F. Valeev
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| |
Collapse
|
7
|
Scott CJC, Di Remigio R, Crawford TD, Thom AJW. Theory and implementation of a novel stochastic approach to coupled cluster. J Chem Phys 2020; 153:144117. [DOI: 10.1063/5.0026513] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Charles J. C. Scott
- Department of Physics, King’s College London, Strand, London WC2R 2LS, United Kingdom
| | - Roberto Di Remigio
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - T. Daniel Crawford
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
- Molecular Sciences Software Institute, Blacksburg, Virginia 24060, USA
| | - Alex J. W. Thom
- Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| |
Collapse
|
8
|
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
|
9
|
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: 502] [Impact Index Per Article: 125.5] [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
|
10
|
Park JW, Al-Saadon R, MacLeod MK, Shiozaki T, Vlaisavljevich B. Multireference Electron Correlation Methods: Journeys along Potential Energy Surfaces. Chem Rev 2020; 120:5878-5909. [PMID: 32239929 DOI: 10.1021/acs.chemrev.9b00496] [Citation(s) in RCA: 77] [Impact Index Per Article: 19.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Multireference electron correlation methods describe static and dynamical electron correlation in a balanced way and, therefore, can yield accurate and predictive results even when single-reference methods or multiconfigurational self-consistent field theory fails. One of their most prominent applications in quantum chemistry is the exploration of potential energy surfaces. This includes the optimization of molecular geometries, such as equilibrium geometries and conical intersections and on-the-fly photodynamics simulations, both of which depend heavily on the ability of the method to properly explore the potential energy surface. Because such applications require nuclear gradients and derivative couplings, the availability of analytical nuclear gradients greatly enhances the scope of quantum chemical methods. This review focuses on the developments and advances made in the past two decades. A detailed account of the analytical nuclear gradient and derivative coupling theories is presented. Emphasis is given to the software infrastructure that allows one to make use of these methods. Notable applications of multireference electron correlation methods to chemistry, including geometry optimizations and on-the-fly dynamics, are summarized at the end followed by a discussion of future prospects.
Collapse
Affiliation(s)
- Jae Woo Park
- Department of Chemistry, Chungbuk National University, Chungdae-ro 1, Cheongju 28644, Korea
| | - Rachael Al-Saadon
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States
| | - Matthew K MacLeod
- Workday, 4900 Pearl Circle East, Suite 100, Boulder, Colorado 80301, United States
| | - Toru Shiozaki
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States.,Quantum Simulation Technologies, Inc., 625 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Bess Vlaisavljevich
- Department of Chemistry, University of South Dakota, 414 East Clark Street, Vermillion, South Dakota 57069, United States
| |
Collapse
|
11
|
Kats D, Werner HJ. Multi-state local complete active space second-order perturbation theory using pair natural orbitals (PNO-MS-CASPT2). J Chem Phys 2019; 150:214107. [DOI: 10.1063/1.5097644] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Daniel Kats
- Max-Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Hans-Joachim Werner
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| |
Collapse
|
12
|
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
|
13
|
Kats D, Köhn A. On the distinguishable cluster approximation for triple excitations. J Chem Phys 2019; 150:151101. [DOI: 10.1063/1.5096343] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- 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, D-70569 Stuttgart, Germany
| |
Collapse
|
14
|
Krause C, Werner HJ. Scalable Electron Correlation Methods. 6. Local Spin-Restricted Open-Shell Second-Order Møller-Plesset Perturbation Theory Using Pair Natural Orbitals: PNO-RMP2. J Chem Theory Comput 2019; 15:987-1005. [PMID: 30571916 DOI: 10.1021/acs.jctc.8b01012] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present a (near) linear scaling implementation of high-spin open-shell Møller-Plesset perturbation theory using pair natural orbitals (PNO-RMP2). The theory is based on a new variant of open-shell MP2 which is fully spin-adapted and uses a single set of spin-free amplitudes, as in closed-shell MP2. This method, denoted SROMP2, is invariant to unitary orbital transformations within the closed, open, and virtual orbital subspaces. Accordingly, only a single set of PNOs per spatial orbital pair is needed, and the efficiency is similar to closed-shell calculations. The PNOs are obtained using a semicanonical approximation with large domains of projected atomic orbitals (PAOs). Linear scaling is achieved provided that the open-shell orbitals are local, and distant pairs are treated by multipole approximations. The method is efficiently parallelized. The convergence of ionization and reaction energies as a function of the PAO and PNO domain sizes is demonstrated and found to be very similar as for closed-shell calculations. The suitability of the PNOs for explicitly correlated PNO-RCCSD-F12 calculations is also tested. So far, this method is only simulated using a conventional program with appropriate projections to the PAO and PNO subspaces. It is demonstrated for radical stabilization energies as well as ionization potentials that the errors caused by the local domain approximations with our default thresholds are negligible.
Collapse
Affiliation(s)
- Christine Krause
- 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
|
15
|
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
|
16
|
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
|
17
|
Ma Q, Werner H. Explicitly correlated local coupled‐cluster methods using pair natural orbitals. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2018. [DOI: 10.1002/wcms.1371] [Citation(s) in RCA: 100] [Impact Index Per Article: 16.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Affiliation(s)
- Qianli Ma
- Institute for Theoretical ChemistryUniversity of StuttgartStuttgartGermany
| | | |
Collapse
|
18
|
Rolik Z, Kállay M. Novel strategy to implement active-space coupled-cluster methods. J Chem Phys 2018; 148:124108. [PMID: 29604813 DOI: 10.1063/1.5004971] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A new approach is presented for the efficient implementation of coupled-cluster (CC) methods including higher excitations based on a molecular orbital space partitioned into active and inactive orbitals. In the new framework, the string representation of amplitudes and intermediates is used as long as it is beneficial, but the contractions are evaluated as matrix products. Using a new diagrammatic technique, the CC equations are represented in a compact form due to the string notations we introduced. As an application of these ideas, a new automated implementation of the single-reference-based multi-reference CC equations is presented for arbitrary excitation levels. The new program can be considered as an improvement over the previous implementations in many respects; e.g., diagram contributions are evaluated by efficient vectorized subroutines. Timings for test calculations for various complete active-space problems are presented. As an application of the new code, the weak interactions in the Be dimer were studied.
Collapse
Affiliation(s)
- Zoltán Rolik
- MTA-BME "Lendület" Quantum Chemistry Research Group, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, H-1521 Budapest, Hungary
| | - Mihály Kállay
- MTA-BME "Lendület" Quantum Chemistry Research Group, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, H-1521 Budapest, Hungary
| |
Collapse
|
19
|
Usvyat D, Maschio L, Schütz M. Periodic and fragment models based on the local correlation approach. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2018. [DOI: 10.1002/wcms.1357] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Denis Usvyat
- Institut für ChemieHumboldt‐Universität zu BerlinBerlinGermany
| | - Lorenzo Maschio
- Dipartimento di Chimica and NIS (Nanostructured Interfaces and Surfaces) CentreUniversità di TorinoTorinoItaly
| | - Martin Schütz
- Institut für ChemieHumboldt‐Universität zu BerlinBerlinGermany
| |
Collapse
|
20
|
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
|
21
|
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
|
22
|
Peng B, Kowalski K. Highly Efficient and Scalable Compound Decomposition of Two-Electron Integral Tensor and Its Application in Coupled Cluster Calculations. J Chem Theory Comput 2017; 13:4179-4192. [DOI: 10.1021/acs.jctc.7b00605] [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)
- Bo Peng
- William R. Wiley Environmental
Molecular Sciences Laboratory, Battelle, Pacific Northwest National Laboratory, K8-91, P. O. Box 999, Richland, Washington 99352, United States
| | - Karol Kowalski
- William R. Wiley Environmental
Molecular Sciences Laboratory, Battelle, Pacific Northwest National Laboratory, K8-91, P. O. Box 999, Richland, Washington 99352, United States
| |
Collapse
|
23
|
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
|
24
|
Schütz M, Maschio L, Karttunen AJ, Usvyat D. Exfoliation Energy of Black Phosphorus Revisited: A Coupled Cluster Benchmark. J Phys Chem Lett 2017; 8:1290-1294. [PMID: 28248525 DOI: 10.1021/acs.jpclett.7b00253] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Black phosphorus (black-P) consists of phosphorene sheets, stacked by van der Waals dispersion. In a recent study based on periodic local second-order Møller-Plesset perturbation theory (LMP2) with higher-order corrections evaluated on finite clusters, we obtained a value of -151 meV/atom for the exfoliation energy. This is almost twice as large as another recent theoretical result (around -80 meV/atom) obtained with quantum Monte Carlo (QMC). Here, we revisit this system on the basis of the recently implemented, periodically embedded ring coupled cluster (rCCD) model instead of LMP2. Higher-order coupled cluster corrections on top of rCCD are obtained from finite clusters by utilizing our new "unit-cell-in-cluster" scheme. Our new value of -92 meV/atom is noticeably lower than that based on LMP2 and in reasonably close agreement with the QMC result. However, in contrast to QMC, no strong effect from the second-neighbor and farther layers in black-P are observed in our calculations.
Collapse
Affiliation(s)
- Martin Schütz
- Institut für Chemie, Humboldt-Universität zu Berlin , Brook-Taylor-Strasse 2, D-12489 Berlin, Germany
| | - Lorenzo Maschio
- Dipartimento di Chimica and NIS (Nanostructured Interfaces and Surfaces) Centre, Università di Torino , via Giuria 5, I-10125 Torino, Italy
| | - Antti J Karttunen
- Department of Chemistry and Materials Science, Aalto University , Kemistintie 1, 02150 Espoo, Finland
| | - Denis Usvyat
- Institut für Chemie, Humboldt-Universität zu Berlin , Brook-Taylor-Strasse 2, D-12489 Berlin, Germany
| |
Collapse
|
25
|
Manzer S, Epifanovsky E, Krylov AI, Head-Gordon M. A General Sparse Tensor Framework for Electronic Structure Theory. J Chem Theory Comput 2017; 13:1108-1116. [DOI: 10.1021/acs.jctc.6b00853] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Samuel Manzer
- Kenneth
S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California − Berkeley, Berkeley, California 94720, United States
- Chemical
Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Evgeny Epifanovsky
- Q-Chem, Inc., 6601 Owens
Drive, Suite 105, Pleasanton, California 94588, United States
| | - Anna I. Krylov
- Department
of Chemistry, University of Southern California, Los Angeles, California 90089, United States
| | - Martin Head-Gordon
- Kenneth
S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California − Berkeley, Berkeley, California 94720, United States
- Chemical
Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
26
|
Menezes F, Kats D, Werner HJ. Local complete active space second-order perturbation theory using pair natural orbitals (PNO-CASPT2). J Chem Phys 2016; 145:124115. [DOI: 10.1063/1.4963019] [Citation(s) in RCA: 68] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
27
|
Masur O, Schütz M, Maschio L, Usvyat D. Fragment-Based Direct-Local-Ring-Coupled-Cluster Doubles Treatment Embedded in the Periodic Hartree–Fock Solution. J Chem Theory Comput 2016; 12:5145-5156. [DOI: 10.1021/acs.jctc.6b00651] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Oliver Masur
- Institut
für Physikalische und Theoretische Chemie, Universität Regensburg, 93049 Regensburg, Germany
| | - Martin Schütz
- Institut
für Physikalische und Theoretische Chemie, Universität Regensburg, 93049 Regensburg, Germany
| | - Lorenzo Maschio
- Dipartimento
di Chimica, and NIS (Nanostructured Interfaces and Surfaces) Centre, Università di Torino, via Giuria 5, I-10125 Torino, Italy
| | - Denis Usvyat
- Institut
für Physikalische und Theoretische Chemie, Universität Regensburg, 93049 Regensburg, Germany
| |
Collapse
|
28
|
Kats D. Speeding up local correlation methods: System-inherent domains. J Chem Phys 2016; 145:014103. [DOI: 10.1063/1.4954963] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Daniel Kats
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| |
Collapse
|
29
|
Kats D. The distinguishable cluster approach from a screened Coulomb formalism. J Chem Phys 2016; 144:044102. [PMID: 26827197 DOI: 10.1063/1.4940398] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The distinguishable cluster doubles equations have been derived starting from an effective screened Coulomb formalism and a particle-hole symmetric formulation of the Fock matrix. A perturbative triples correction to the distinguishable cluster with singles and doubles (DCSD) has been introduced employing the screened integrals. It is shown that the resulting DCSD(T) method is more accurate than DCSD for reaction energies and is less sensitive to the static correlation than coupled cluster with singles and doubles with a perturbative triples correction.
Collapse
Affiliation(s)
- Daniel Kats
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| |
Collapse
|
30
|
Usvyat D, Maschio L, Schütz M. Periodic local MP2 method employing orbital specific virtuals. J Chem Phys 2015; 143:102805. [DOI: 10.1063/1.4921301] [Citation(s) in RCA: 41] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Denis Usvyat
- Institute for Physical and Theoretical Chemistry, Universität Regensburg, Universitätsstraße 31, D-93040 Regensburg, Germany
| | - Lorenzo Maschio
- Dipartimento di Chimica, and Centre of Excellence NIS (Nanostructured Interfaces and Surfaces), Università di Torino, via Giuria 5, I-10125 Torino, Italy
| | - Martin Schütz
- Institute for Physical and Theoretical Chemistry, Universität Regensburg, Universitätsstraße 31, D-93040 Regensburg, Germany
| |
Collapse
|
31
|
Tanaka S, Mochizuki Y, Komeiji Y, Okiyama Y, Fukuzawa K. Electron-correlated fragment-molecular-orbital calculations for biomolecular and nano systems. Phys Chem Chem Phys 2015; 16:10310-44. [PMID: 24740821 DOI: 10.1039/c4cp00316k] [Citation(s) in RCA: 189] [Impact Index Per Article: 21.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Recent developments in the fragment molecular orbital (FMO) method for theoretical formulation, implementation, and application to nano and biomolecular systems are reviewed. The FMO method has enabled ab initio quantum-mechanical calculations for large molecular systems such as protein-ligand complexes at a reasonable computational cost in a parallelized way. There have been a wealth of application outcomes from the FMO method in the fields of biochemistry, medicinal chemistry and nanotechnology, in which the electron correlation effects play vital roles. With the aid of the advances in high-performance computing, the FMO method promises larger, faster, and more accurate simulations of biomolecular and related systems, including the descriptions of dynamical behaviors in solvent environments. The current status and future prospects of the FMO scheme are addressed in these contexts.
Collapse
Affiliation(s)
- Shigenori Tanaka
- Graduate School of System Informatics, Kobe University, 1-1 Rokkodai, Nada-ku, Kobe 657-8501, Japan.
| | | | | | | | | |
Collapse
|
32
|
Vereecken L, Glowacki DR, Pilling MJ. Theoretical Chemical Kinetics in Tropospheric Chemistry: Methodologies and Applications. Chem Rev 2015; 115:4063-114. [DOI: 10.1021/cr500488p] [Citation(s) in RCA: 146] [Impact Index Per Article: 16.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Luc Vereecken
- Max Planck Institute for Chemistry, 55128 Mainz, Germany
| | - David R. Glowacki
- PULSE
Institute and Department of Chemistry, Stanford University, Stanford, California 94305, United States
- SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
- School
of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
- Department
of Computer Science, University of Bristol, Bristol BS8 1UB, United Kingdom
| | | |
Collapse
|
33
|
Kats D, Kreplin D, Werner HJ, Manby FR. Accurate thermochemistry from explicitly correlated distinguishable cluster approximation. J Chem Phys 2015; 142:064111. [DOI: 10.1063/1.4907591] [Citation(s) in RCA: 44] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Affiliation(s)
- Daniel Kats
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - David Kreplin
- 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
| | - Frederick R. Manby
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS,United Kingdom
| |
Collapse
|
34
|
MacLeod MK, Shiozaki T. Communication: Automatic code generation enables nuclear gradient computations for fully internally contracted multireference theory. J Chem Phys 2015; 142:051103. [DOI: 10.1063/1.4907717] [Citation(s) in RCA: 78] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023] Open
Affiliation(s)
- Matthew K. MacLeod
- Department of Chemistry, Northwestern University, 2145 Sheridan Rd., Evanston, Illinois 60208, USA
| | - Toru Shiozaki
- Department of Chemistry, Northwestern University, 2145 Sheridan Rd., Evanston, Illinois 60208, USA
| |
Collapse
|
35
|
Affiliation(s)
- Daniel Kats
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| |
Collapse
|
36
|
Kats D. Communication: The distinguishable cluster approximation. II. The role of orbital relaxation. J Chem Phys 2014; 141:061101. [DOI: 10.1063/1.4892792] [Citation(s) in RCA: 44] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Affiliation(s)
- Daniel Kats
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| |
Collapse
|
37
|
Schütz M, Masur O, Usvyat D. Efficient and accurate treatment of weak pairs in local CCSD(T) calculations. II. Beyond the ring approximation. J Chem Phys 2014; 140:244107. [DOI: 10.1063/1.4884156] [Citation(s) in RCA: 39] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Martin Schütz
- Institut für Physikalische und Theoretische Chemie, Universität Regensburg, Regensburg D-93040, Germany
| | - Oliver Masur
- Institut für Physikalische und Theoretische Chemie, Universität Regensburg, Regensburg D-93040, Germany
| | - Denis Usvyat
- Institut für Physikalische und Theoretische Chemie, Universität Regensburg, Regensburg D-93040, Germany
| |
Collapse
|
38
|
Rolik Z, Szegedy L, Ladjánszki I, Ladóczki B, Kállay M. An efficient linear-scaling CCSD(T) method based on local natural orbitals. J Chem Phys 2014; 139:094105. [PMID: 24028100 DOI: 10.1063/1.4819401] [Citation(s) in RCA: 288] [Impact Index Per Article: 28.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
An improved version of our general-order local coupled-cluster (CC) approach [Z. Rolik and M. Kállay, J. Chem. Phys. 135, 104111 (2011)] and its efficient implementation at the CC singles and doubles with perturbative triples [CCSD(T)] level is presented. The method combines the cluster-in-molecule approach of Li and co-workers [J. Chem. Phys. 131, 114109 (2009)] with frozen natural orbital (NO) techniques. To break down the unfavorable fifth-power scaling of our original approach a two-level domain construction algorithm has been developed. First, an extended domain of localized molecular orbitals (LMOs) is assembled based on the spatial distance of the orbitals. The necessary integrals are evaluated and transformed in these domains invoking the density fitting approximation. In the second step, for each occupied LMO of the extended domain a local subspace of occupied and virtual orbitals is constructed including approximate second-order Mo̸ller-Plesset NOs. The CC equations are solved and the perturbative corrections are calculated in the local subspace for each occupied LMO using a highly-efficient CCSD(T) code, which was optimized for the typical sizes of the local subspaces. The total correlation energy is evaluated as the sum of the individual contributions. The computation time of our approach scales linearly with the system size, while its memory and disk space requirements are independent thereof. Test calculations demonstrate that currently our method is one of the most efficient local CCSD(T) approaches and can be routinely applied to molecules of up to 100 atoms with reasonable basis sets.
Collapse
Affiliation(s)
- Zoltán Rolik
- MTA-BME Lendület Quantum Chemistry Research Group, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, H-1521 Budapest, P.O. Box 91, Hungary
| | | | | | | | | |
Collapse
|
39
|
Masur O, Usvyat D, Schütz M. Efficient and accurate treatment of weak pairs in local CCSD(T) calculations. J Chem Phys 2013; 139:164116. [DOI: 10.1063/1.4826534] [Citation(s) in RCA: 37] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|