1
|
Dang DK, Wilson LW, Zimmerman PM. The numerical evaluation of Slater integrals on graphics processing units. J Comput Chem 2022; 43:1680-1689. [PMID: 35861566 DOI: 10.1002/jcc.26968] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2021] [Revised: 05/13/2022] [Accepted: 06/27/2022] [Indexed: 11/11/2022]
Abstract
This article presents SlaterGPU, a graphics processing unit (GPU) accelerated library that uses OpenACC to numerically compute Slater-type orbital (STO) integrals. The electron repulsion integrals (ERI) are computed under the RI approximation using the Coulomb potential of the Slater basis function. To fully realize the performance capabilities of modern GPUs, the Slater integrals are evaluated in mixed-precision, resulting in speedups for the ERIs of over 80×. Parallelization on multiple GPUs allows for integral throughput of over 3 million integrals per second. This places STO integral throughput within reach of single-threaded, conventional Gaussian integration schemes. To test the quality of the integrals, the fluorine exchange reaction barrier in fluoromethane was computed using heat-bath configuration interaction (HBCI). In addition, the singlet-triplet gap of cyclobutadiene was examined using HBCI in a triple- ζ $$ \zeta $$ , polarized basis set. These benchmarks demonstrate the library's ability to generate the full set of integrals necessary for configuration interaction with up to 6 h $$ 6h $$ functions in the auxiliary basis.
Collapse
Affiliation(s)
- Duy-Khoi Dang
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan, USA
| | - Leighton W Wilson
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan, USA.,Department of Mathematics, University of Michigan, Ann Arbor, Michigan, USA
| | - Paul M Zimmerman
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan, USA
| |
Collapse
|
2
|
Wang S, Li C, Evangelista FA. Analytic Energy Gradients for the Driven Similarity Renormalization Group Multireference Second-Order Perturbation Theory. J Chem Theory Comput 2021; 17:7666-7681. [PMID: 34839660 DOI: 10.1021/acs.jctc.1c00980] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We derive analytic energy gradients of the driven similarity renormalization group (DSRG) multireference second-order perturbation theory (MRPT2) using the method of Lagrange multipliers. In the Lagrangian, we impose constraints for a complete-active-space self-consistent-field reference wave function and the semicanonical orthonormal molecular orbitals. Solving for the associated Lagrange multipliers is found to share the same asymptotic scaling of a single DSRG-MRPT2 energy computation. A pilot implementation of the DSRG-MRPT2 analytic gradients is used to optimize the geometry of the singlet and triplet states of p-benzyne. The equilibrium bond lengths and angles are similar to those computed via other MRPT2s and Mukherjee's multireference coupled cluster theory. An approximate DSRG-MRPT2 method that neglects the contributions of the three-body density cumulant is found to introduce negligible errors in the geometry of p-benzyne, lending itself to a promising low-cost approach for molecular geometry optimizations using large active spaces.
Collapse
Affiliation(s)
- Shuhe Wang
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Chenyang Li
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Francesco A Evangelista
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
3
|
Relativistic Fock Space Coupled Cluster Method for Many-Electron Systems: Non-Perturbative Account for Connected Triple Excitations. Symmetry (Basel) 2020. [DOI: 10.3390/sym12071101] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Abstract
The Fock space relativistic coupled cluster method (FS-RCC) is one of the most promising tools of electronic structure modeling for atomic and molecular systems containing heavy nuclei. Until recently, capabilities of the FS-RCC method were severely restricted by the fact that only single and double excitations in the exponential parametrization of the wave operator were considered. We report the design and the first computer implementation of FS-RCC schemes with full and simplified non-perturbative account for triple excitations in the cluster operator. Numerical stability of the new computational scheme and thus its applicability to a wide variety of molecular electronic states is ensured using the dynamic shift technique combined with the extrapolation to zero-shift limit. Pilot applications to atomic (Tl, Pb) and molecular (TlH) systems reported in the paper indicate that the breakthrough in accuracy and predictive power of the electronic structure calculations for heavy-element compounds can be achieved. Moreover, the described approach can provide a firm basis for high-precision modeling of heavy molecular systems with several open shells, including actinide compounds.
Collapse
|
4
|
Li C, Evangelista FA. Connected three-body terms in single-reference unitary many-body theories: Iterative and perturbative approximations. J Chem Phys 2020; 152:234116. [PMID: 32571040 DOI: 10.1063/5.0008353] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
This work introduces various approaches to include connected three-body terms in unitary many-body theories, focusing on the driven similarity renormalization group (DSRG). Starting from the least approximate method-the linearized DSRG truncated to one-, two-, and three-body operators [LDSRG(3)]-we develop several approximate LDSRG(3) models with reduced computational cost. Through a perturbative analysis, we motivate a family of iterative LDSRG(3)-n and -n' (n = 1, 2, 3, 4) methods that contain a subset of the LDSRG(3) diagrams. Among these variants, the LDSRG(3)-2 scheme has the same computational complexity of coupled cluster theory with singles, doubles, and triples (CCSDT), but it outperforms CCSDT in the accuracy of the predicted correlation energies. We also propose and implement two perturbative triples corrections based on the linearized DSRG truncated to one- and two-body operators augmented with recursive semi-quadratic commutators [qDSRG(2)]. The resulting qDSRG(2)+(T) approach matches the accuracy of the "gold-standard" coupled cluster theory with singles, doubles, and perturbative triples model on the energetics of twenty-eight closed-shell atoms and small molecules.
Collapse
Affiliation(s)
- Chenyang Li
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Francesco A Evangelista
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| |
Collapse
|
5
|
Matthews DA, Cheng L, Harding ME, Lipparini F, Stopkowicz S, Jagau TC, Szalay PG, Gauss J, Stanton JF. Coupled-cluster techniques for computational chemistry: The CFOUR program package. J Chem Phys 2020; 152:214108. [DOI: 10.1063/5.0004837] [Citation(s) in RCA: 214] [Impact Index Per Article: 53.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Affiliation(s)
- Devin A. Matthews
- Department of Chemistry, Southern Methodist University, Dallas, Texas 75275, USA
| | - Lan Cheng
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Michael E. Harding
- Institut für Physikalische Chemie, Karlsruher Institut für Technologie (KIT), Kaiserstr. 12, D-76131 Karlsruhe, Germany
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, I-56124 Pisa, Italy
| | - Stella Stopkowicz
- Department Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, D-55128 Mainz, Germany
| | - Thomas-C. Jagau
- Department of Chemistry, University of Munich (LMU), Butenandtstr. 5-13, D-81377 Munich, Germany
| | - Péter G. Szalay
- ELTE Eötvös Loránd University, Institute of Chemistry, Pázmány Péter sétány 1/A, H-1117 Budapest, Hungary
| | - Jürgen Gauss
- Department Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, D-55128 Mainz, Germany
| | - John F. Stanton
- Quantum Theory Project, Departments of Chemistry and Physics, University of Florida, Gainesville, Florida 32611, USA
| |
Collapse
|
6
|
Smith DGA, Burns LA, Simmonett AC, Parrish RM, Schieber MC, Galvelis R, Kraus P, Kruse H, Di Remigio R, Alenaizan A, James AM, Lehtola S, Misiewicz JP, Scheurer M, Shaw RA, Schriber JB, Xie Y, Glick ZL, Sirianni DA, O’Brien JS, Waldrop JM, Kumar A, Hohenstein EG, Pritchard BP, Brooks BR, Schaefer HF, Sokolov AY, Patkowski K, DePrince AE, Bozkaya U, King RA, Evangelista FA, Turney JM, Crawford TD, Sherrill CD. Psi4 1.4: Open-source software for high-throughput quantum chemistry. J Chem Phys 2020; 152:184108. [PMID: 32414239 PMCID: PMC7228781 DOI: 10.1063/5.0006002] [Citation(s) in RCA: 333] [Impact Index Per Article: 83.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2020] [Accepted: 04/12/2020] [Indexed: 12/13/2022] Open
Abstract
PSI4 is a free and open-source ab initio electronic structure program providing implementations of Hartree-Fock, density functional theory, many-body perturbation theory, configuration interaction, density cumulant theory, symmetry-adapted perturbation theory, and coupled-cluster theory. Most of the methods are quite efficient, thanks to density fitting and multi-core parallelism. The program is a hybrid of C++ and Python, and calculations may be run with very simple text files or using the Python API, facilitating post-processing and complex workflows; method developers also have access to most of PSI4's core functionalities via Python. Job specification may be passed using The Molecular Sciences Software Institute (MolSSI) QCSCHEMA data format, facilitating interoperability. A rewrite of our top-level computation driver, and concomitant adoption of the MolSSI QCARCHIVE INFRASTRUCTURE project, makes the latest version of PSI4 well suited to distributed computation of large numbers of independent tasks. The project has fostered the development of independent software components that may be reused in other quantum chemistry programs.
Collapse
Affiliation(s)
| | - Lori A. Burns
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| | - Andrew C. Simmonett
- National Institutes of Health – National Heart,
Lung and Blood Institute, Laboratory of Computational Biology, Bethesda,
Maryland 20892, USA
| | - Robert M. Parrish
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| | - Matthew C. Schieber
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| | | | - Peter Kraus
- School of Molecular and Life Sciences, Curtin
University, Kent St., Bentley, Perth, Western Australia 6102,
Australia
| | - Holger Kruse
- Institute of Biophysics of the Czech Academy of
Sciences, Královopolská 135, 612 65 Brno, Czech
Republic
| | - Roberto Di Remigio
- Department of Chemistry, Centre for Theoretical
and Computational Chemistry, UiT, The Arctic University of Norway, N-9037
Tromsø, Norway
| | - Asem Alenaizan
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| | - Andrew M. James
- Department of Chemistry, Virginia
Tech, Blacksburg, Virginia 24061, USA
| | - Susi Lehtola
- Department of Chemistry, University of
Helsinki, P.O. Box 55 (A. I. Virtasen aukio 1), FI-00014 Helsinki,
Finland
| | - Jonathon P. Misiewicz
- Center for Computational Quantum Chemistry,
University of Georgia, Athens, Georgia 30602, USA
| | - Maximilian Scheurer
- Interdisciplinary Center for Scientific
Computing, Heidelberg University, D-69120 Heidelberg,
Germany
| | - Robert A. Shaw
- ARC Centre of Excellence in Exciton Science,
School of Science, RMIT University, Melbourne, VIC 3000,
Australia
| | - Jeffrey B. Schriber
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| | - Yi Xie
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| | - Zachary L. Glick
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| | - Dominic A. Sirianni
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| | - Joseph Senan O’Brien
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| | - Jonathan M. Waldrop
- Department of Chemistry and Biochemistry, Auburn
University, Auburn, Alabama 36849, USA
| | - Ashutosh Kumar
- Department of Chemistry, Virginia
Tech, Blacksburg, Virginia 24061, USA
| | - Edward G. Hohenstein
- SLAC National Accelerator Laboratory, Stanford
PULSE Institute, Menlo Park, California 94025,
USA
| | | | - Bernard R. Brooks
- National Institutes of Health – National Heart,
Lung and Blood Institute, Laboratory of Computational Biology, Bethesda,
Maryland 20892, USA
| | - Henry F. Schaefer
- Center for Computational Quantum Chemistry,
University of Georgia, Athens, Georgia 30602, USA
| | - Alexander Yu. Sokolov
- Department of Chemistry and Biochemistry, The
Ohio State University, Columbus, Ohio 43210, USA
| | - Konrad Patkowski
- Department of Chemistry and Biochemistry, Auburn
University, Auburn, Alabama 36849, USA
| | - A. Eugene DePrince
- Department of Chemistry and Biochemistry,
Florida State University, Tallahassee, Florida 32306-4390,
USA
| | - Uğur Bozkaya
- Department of Chemistry, Hacettepe
University, Ankara 06800, Turkey
| | - Rollin A. King
- Department of Chemistry, Bethel
University, St. Paul, Minnesota 55112, USA
| | | | - Justin M. Turney
- Center for Computational Quantum Chemistry,
University of Georgia, Athens, Georgia 30602, USA
| | | | - C. David Sherrill
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| |
Collapse
|
7
|
Wu CH, Magers DB, Harding LB, Klippenstein SJ, Allen WD. Reaction Profiles and Kinetics for Radical-Radical Hydrogen Abstraction via Multireference Coupled Cluster Theory. J Chem Theory Comput 2020; 16:1511-1525. [PMID: 32073856 DOI: 10.1021/acs.jctc.9b00966] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Radical-radical abstractions in hydrocarbon oxidation chemistry are disproportionation reactions that are generally exothermic with little or no barrier yet are underappreciated and poorly studied. Such challenging multireference electronic structure problems are tackled here using the recently developed state-specific multireference coupled cluster methods Mk-MRCCSD and Mk-MRCCSD(T), as well as the companion perturbation theory Mk-MRPT2 and the established MRCISD, MRCISD+Q, and CASPT2 approaches. Reaction paths are investigated for five prototypes involving radical-radical hydrogen abstraction: H + BeH → H2+ Be, H + NH2 → H2 + NH, CH3 + C2H5 → CH4 + C2H4, H + C2H5 → H2 + C2H4, and H + HCO → H2 + CO. Full configuration interaction (FCI) benchmark computations for the H + BeH, H + NH2, and H + HCO reactions prove that Mk-MRCCSD(T) provides superior accuracy for the interaction energies in the entrance channel, with mean absolute errors less than 0.3 kcal mol-1 and percentage deviations less than 10% over the fragment separations of relevance to kinetics. To facilitate combustion studies, energetics for the CH3 + C2H5, H + C2H5, and H + HCO reactions were computed at each level of theory with correlation-consistent basis sets (cc-pVXZ, X = T, Q, 5) and extrapolated to the complete basis set (CBS) limit. These CBS energies were coupled with CASPT2 projected vibrational frequencies along a minimum energy path to obtain rate constants for these three reactions. The rigorous Mk-MRCCSD(T)/CBS results demonstrate unequivocally that these three reactions proceed with no barrier in the entrance channel, contrary to some earlier predictions. Mk-MRCCSD(T) also reveals that the economical CASPT2 method performs well for large interfragment separations but may deteriorate substantially at shorter distances.
Collapse
Affiliation(s)
- Chia-Hua Wu
- Center for Computational Quantum Chemistry and Department of Chemistry, University of Georgia, Athens, Georgia 30602, United States
| | - D Brandon Magers
- Center for Computational Quantum Chemistry and Department of Chemistry, University of Georgia, Athens, Georgia 30602, United States.,Department of Chemistry and Physics, Belhaven University, Jackson, Mississippi 39202, United States
| | - Lawrence B Harding
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Stephen J Klippenstein
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Wesley D Allen
- Center for Computational Quantum Chemistry and Department of Chemistry, University of Georgia, Athens, Georgia 30602, United States
| |
Collapse
|
8
|
Musiał M, Meissner L, Cembrzynska J. The intermediate Hamiltonian Fock-space coupled-cluster method with approximate evaluation of the three-body effects. J Chem Phys 2019; 151:184102. [DOI: 10.1063/1.5124806] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Monika Musiał
- Institute of Chemistry, University of Silesia in Katowice, Szkolna 9, 40-006 Katowice, Poland
| | - Leszek Meissner
- Institute of Physics, Nicholaus Copernicus University, Grudziadzka 5/7, 87-100 Toruń, Poland
| | - Justyna Cembrzynska
- Institute of Physics, Nicholaus Copernicus University, Grudziadzka 5/7, 87-100 Toruń, Poland
| |
Collapse
|
9
|
Manna S, Ray SS, Chattopadhyay S, Chaudhuri RK. A simplified account of the correlation effects to bond breaking processes: The Brillouin-Wigner perturbation theory using a multireference formulation. J Chem Phys 2019. [DOI: 10.1063/1.5097657] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Shovan Manna
- Department of Chemistry, Indian Institute of Engineering Science and Technology, Shibpur, Howrah 711103, India
| | - Suvonil Sinha Ray
- Department of Chemistry, University of Calcutta, Kolkata 700009, India
| | - Sudip Chattopadhyay
- Department of Chemistry, Indian Institute of Engineering Science and Technology, Shibpur, Howrah 711103, India
| | | |
Collapse
|
10
|
Zhang T, Li C, Evangelista FA. Improving the Efficiency of the Multireference Driven Similarity Renormalization Group via Sequential Transformation, Density Fitting, and the Noninteracting Virtual Orbital Approximation. J Chem Theory Comput 2019; 15:4399-4414. [PMID: 31268704 DOI: 10.1021/acs.jctc.9b00353] [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
This study examines several techniques to improve the efficiency of the linearized multireference driven similarity renormalization group truncated to one- and two-body operators [MR-LDSRG(2)]. We propose a sequential MR-LDSRG(2) [sq-MR-LDSRG(2)] scheme, in which one-body substitutions are folded exactly into the Hamiltonian. This new approach is combined with density fitting (DF) to reduce the storage cost of two-electron integrals. To further avoid storage of large four-index intermediates, we propose a noninteracting virtual orbital (NIVO) approximation of the Baker-Campbell-Hausdorff series that neglects commutators terms with three and four virtual indices. The NIVO approximation reduces the computational prefactor of the MR-LDSRG(2), bringing it closer to that of coupled cluster with singles and doubles (CCSD). We test the effect of the DF and NIVO approximations on the MR-LDSRG(2) and sq-MR-LDSRG(2) methods by computing properties of eight diatomic molecules. The diatomic constants obtained by DF-sq-MR-LDSRG(2)+NIVO are found to be as accurate as those from the original MR-LDSRG(2) and coupled cluster theory with singles, doubles, and perturbative triples. Finally, we demonstrate that the DF-sq-MR-LDSRG(2)+NIVO scheme can be applied to chemical systems with more than 550 basis functions by computing the automerization energy of cyclobutadiene with a quintuple-ζ basis set. The predicted automerization energy is found to be similar to the value computed with Mukherjee's state-specific multireference coupled cluster theory with singles and doubles.
Collapse
Affiliation(s)
- Tianyuan Zhang
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation , Emory University , Atlanta , Georgia 30322 , United States
| | - Chenyang Li
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation , Emory University , Atlanta , Georgia 30322 , United States
| | - Francesco A Evangelista
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation , Emory University , Atlanta , Georgia 30322 , United States
| |
Collapse
|
11
|
Lang J, Brabec J, Saitow M, Pittner J, Neese F, Demel O. Perturbative triples correction to domain-based local pair natural orbital variants of Mukherjee's state specific coupled cluster method. Phys Chem Chem Phys 2019; 21:5022-5038. [PMID: 30762044 DOI: 10.1039/c8cp03577f] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In this article we report an implementation of the perturbative triples correction to Mukherjee's state-specific multireference coupled cluster method based on the domain-based pair natural orbital approach (DLPNO-MkCC). We tested the performance of DLPNO-MkCCSD(T) in calculations involving tetramethyleneethane and isomers of naphthynes. These tests show that more than 97% of triples energy was recovered with respect to the canonical MkCCSD(T) method, which together with the DLPNO-MkCCSD part accounts for about 99.70-99.85% of the total correlation energy. The applicability of the method was demonstrated on calculations of singlet-triplet gaps for several large systems: triangulene, dynemicin A, and a beryllium complex.
Collapse
Affiliation(s)
- Jakub Lang
- J. Heyrovský Institute of Physical Chemistry, v.v.i., Academy of Sciences of the Czech Republic, Dolejškova 3, 18223 Prague 8, Czech Republic.
| | | | | | | | | | | |
Collapse
|
12
|
Tada K, Maruyama T, Koga H, Okumura M, Tanaka S. Extent of Spin Contamination Errors in DFT/Plane-wave Calculation of Surfaces: A Case of Au Atom Aggregation on a MgO Surface. Molecules 2019; 24:molecules24030505. [PMID: 30704148 PMCID: PMC6385026 DOI: 10.3390/molecules24030505] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/28/2018] [Revised: 01/28/2019] [Accepted: 01/30/2019] [Indexed: 11/16/2022] Open
Abstract
The aggregation of Au atoms onto a Au dimer (Au₂) on a MgO (001) surface was calculated by restricted (spin-un-polarized) and unrestricted (spin-polarized) density functional theory calculations with a plane-wave basis and the approximate spin projection (AP) method. The unrestricted calculations included spin contamination errors of 0.0⁻0.1 eV, and the errors were removed using the AP method. The potential energy curves for the aggregation reaction estimated by the restricted and unrestricted calculations were different owing to the estimation of the open-shell structure by the unrestricted calculations. These results show the importance of the open-shell structure and correction of the spin contamination error for the calculation of small-cluster-aggregations and molecule dimerization on surfaces.
Collapse
Affiliation(s)
- Kohei Tada
- Research Institute of Electrochemical Energy, National Institute of Advanced Industrial Science and Technology (AIST), 1-8-31, Midorigaoka, Ikeda, Osaka 563-8577, Japan.
| | - Tomohiro Maruyama
- Department of Chemistry, Graduate School of Science, Osaka University, 1-1, Machikaneyama, Toyonaka, Osaka 560-0043, Japan.
| | - Hiroaki Koga
- Elements Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University, 1-30 Goryo Ohara, Nishikyo, Kyoto 615-8245, Japan.
| | - Mitsutaka Okumura
- Department of Chemistry, Graduate School of Science, Osaka University, 1-1, Machikaneyama, Toyonaka, Osaka 560-0043, Japan.
- Elements Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University, 1-30 Goryo Ohara, Nishikyo, Kyoto 615-8245, Japan.
| | - Shingo Tanaka
- Research Institute of Electrochemical Energy, National Institute of Advanced Industrial Science and Technology (AIST), 1-8-31, Midorigaoka, Ikeda, Osaka 563-8577, Japan.
| |
Collapse
|
13
|
Tsuchimochi T, Ten-No SL. Extending spin-symmetry projected coupled-cluster to large model spaces using an iterative null-space projection technique. J Comput Chem 2018; 40:265-278. [PMID: 30520115 DOI: 10.1002/jcc.25587] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2018] [Revised: 07/31/2018] [Accepted: 08/01/2018] [Indexed: 11/06/2022]
Abstract
Recently, we introduced an orbital-invariant approximate coupled-cluster (CC) method in the spin-projection manifold. The multi-determinantal property of spin-projection means that the parametrization in the spin-extended CC (ECC) ansatz is nonorthogonal and overcomplete. Therefore, the linear dependencies must be removed by an orthogonalization procedure to obtain meaningful solutions. Multi-reference methods often achieve this by diagonalizing a metric of the equation system, but this is not feasible with ECC because of the enormous size of the metric, a consequence of the incomplete active space of the spin-projected Hartree-Fock reference. As a result, the applicability of ECC has been limited to small benchmark systems, for which the ansatz was shown to be superior to the configuration interaction and linearized approximations. In this article, we provide a solution to this problem that completely avoids the metric diagonalization by iteratively projecting out its null-space from the working equations. As the additional computational cost required for this iterative projection is only marginal, it greatly expands the application range of ECC. We demonstrate the potential of approximate ECC by studying the complete basis set limit of F2 and transition metal complexes such as NiO, Mn2 , and [Cu2 O2 ]2+ , which have all been hindered by the prohibitively large metric size. We also identify the potential inadequacy of the molecular orbitals given by spin-projected Hartree-Fock in some cases, and propose possible solutions. © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Takashi Tsuchimochi
- Graduate School of Science, Technology, and Innovation, Kobe University, Kobe, Hyogo, 657-0025, Japan
| | - Seiichiro L Ten-No
- Graduate School of Science, Technology, and Innovation, Kobe University, Kobe, Hyogo, 657-0025, Japan
| |
Collapse
|
14
|
Tsuchimochi T, Ten-No SL. Orbital-invariant spin-extended approximate coupled-cluster for multi-reference systems. J Chem Phys 2018; 149:044109. [PMID: 30068163 DOI: 10.1063/1.5036542] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present an approximate treatment of spin-extended coupled-cluster (ECC) based on the spin-projection of the broken-symmetry coupled-cluster (CC) ansatz. ECC completely eliminates the spin-contamination of unrestricted CC and is therefore expected to provide better descriptions of dynamical and static correlation effects, but introduces two distinct problems. The first issue is the emergence of non-terminating amplitude equations, which are caused by the de-excitation effects inherent in symmetry projection operators. In this study, we take a minimalist approach and truncate the Taylor series of the exponential ansatz at a certain order such that the approximation safely recovers the traditional CC without spin-projection. The second issue is that the nonlinear equations of ECC become underdetermined, although consistent, yielding an infinitude of solutions. This problem arises because of the redundancies in the excitation manifold, as is common in other multi-reference approaches. We remove the linear dependencies in ECC by employing an orthogonal projection manifold. We also propose an efficient solver for our method, in which the components are usually sparse but not diagonal-dominant. It is shown that our approach is rigorously orbital-invariant and provides more accurate results than its configuration interaction and linearized CC analogues for chemical systems.
Collapse
Affiliation(s)
- Takashi Tsuchimochi
- Graduate School of Science, Technology, and Innovation, Kobe University, Kobe, Hyogo 657-0025, Japan
| | - Seiichiro L Ten-No
- Graduate School of Science, Technology, and Innovation, Kobe University, Kobe, Hyogo 657-0025, Japan
| |
Collapse
|
15
|
Evangelista FA. Perspective: Multireference coupled cluster theories of dynamical electron correlation. J Chem Phys 2018; 149:030901. [DOI: 10.1063/1.5039496] [Citation(s) in RCA: 69] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Affiliation(s)
- Francesco A. Evangelista
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| |
Collapse
|
16
|
Pastorczak E, Pernal K. Correlation Energy from the Adiabatic Connection Formalism for Complete Active Space Wave Functions. J Chem Theory Comput 2018; 14:3493-3503. [DOI: 10.1021/acs.jctc.8b00213] [Citation(s) in RCA: 26] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Ewa Pastorczak
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland
| | - Katarzyna Pernal
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland
| |
Collapse
|
17
|
Veis L, Antalík A, Legeza Ö, Alavi A, Pittner J. The Intricate Case of Tetramethyleneethane: A Full Configuration Interaction Quantum Monte Carlo Benchmark and Multireference Coupled Cluster Studies. J Chem Theory Comput 2018; 14:2439-2445. [DOI: 10.1021/acs.jctc.8b00022] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Libor Veis
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Andrej Antalík
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
- Faculty of Mathematics and Physics, Charles University, 11636 Prague, Czech Republic
| | - Örs Legeza
- Strongly Correlated Systems “Lendület” Research group, Wigner Research Centre for Physics, H-1525 Budapest, Hungary
| | - Ali Alavi
- Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, United Kingdom
- Max Planck Institüt für Festkörperforschung, 70569 Stuttgart, Germany
| | - Jiří Pittner
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| |
Collapse
|
18
|
Combined complete active space configuration interaction and perturbation theory applied to conformational energy prototypes: Rotation and inversion barriers. COMPUT THEOR CHEM 2017. [DOI: 10.1016/j.comptc.2017.10.003] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
19
|
Klippenstein SJ, Harding LB, Ruscic B. Ab Initio Computations and Active Thermochemical Tables Hand in Hand: Heats of Formation of Core Combustion Species. J Phys Chem A 2017; 121:6580-6602. [DOI: 10.1021/acs.jpca.7b05945] [Citation(s) in RCA: 106] [Impact Index Per Article: 15.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Stephen J. Klippenstein
- Chemical Sciences and Engineering
Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - Lawrence B. Harding
- Chemical Sciences and Engineering
Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - Branko Ruscic
- Chemical Sciences and Engineering
Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| |
Collapse
|
20
|
Lang J, Švaňa M, Demel O, Brabec J, Kedžuch S, Noga J, Kowalski K, Pittner J. A MRCC study of the isomerisation of cyclopropane. Mol Phys 2017. [DOI: 10.1080/00268976.2017.1317112] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Affiliation(s)
- Jakub Lang
- J. Heyrovsky Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, Prague 8, Czech Republic
- Faculty of Science, Charles University in Prague, Prague 2, Czech Republic
| | - Matej Švaňa
- J. Heyrovsky Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, Prague 8, Czech Republic
| | - Ondřej Demel
- J. Heyrovsky Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, Prague 8, Czech Republic
| | - Jiri Brabec
- J. Heyrovsky Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, Prague 8, Czech Republic
| | - Stanislav Kedžuch
- Institute of Inorganic Chemistry, Slovak Academy of Sciences, Bratislava, Slovakia
| | - Jozef Noga
- Institute of Inorganic Chemistry, Slovak Academy of Sciences, Bratislava, Slovakia
- Department of Inorganic Chemistry, Faculty of Natural Sciences, Comenius University, Bratislava, Slovakia
| | - Karol Kowalski
- William R. Wiley Environmental Molecular Sciences Laboratory, Battelle, Pacific Northwest National Laboratory, Richland, WA 99352, USA
| | - Jiří Pittner
- J. Heyrovsky Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, Prague 8, Czech Republic
| |
Collapse
|
21
|
Garniron Y, Giner E, Malrieu JP, Scemama A. Alternative definition of excitation amplitudes in multi-reference state-specific coupled cluster. J Chem Phys 2017; 146:154107. [DOI: 10.1063/1.4980034] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
22
|
Li C, Evangelista FA. Driven similarity renormalization group: Third-order multireference perturbation theory. J Chem Phys 2017; 146:124132. [DOI: 10.1063/1.4979016] [Citation(s) in RCA: 31] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Chenyang Li
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Francesco A. Evangelista
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| |
Collapse
|
23
|
Sinha Ray S, Ghosh A, Shit A, Chaudhuri RK, Chattopadhyay S. A simplified ab initio treatment of diradicaloid structures produced from stretching and breaking chemical bonds. Phys Chem Chem Phys 2017; 19:22282-22301. [DOI: 10.1039/c7cp03564k] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
With a proper choice of active spaces, the single root perturbation theory employing improved virtual orbitals can flawlessly describe the ground, excited, ionized, and dissociated states having varying degrees of degeneracy at the expense of low computational cost.
Collapse
Affiliation(s)
- Suvonil Sinha Ray
- Department of Chemistry
- Indian Institute of Engineering Science and Technology
- Howrah 711103
- India
| | - Anirban Ghosh
- Department of Chemistry
- Indian Institute of Engineering Science and Technology
- Howrah 711103
- India
| | | | | | - Sudip Chattopadhyay
- Department of Chemistry
- Indian Institute of Engineering Science and Technology
- Howrah 711103
- India
| |
Collapse
|
24
|
Lefrancois D, Rehn DR, Dreuw A. Accurate adiabatic singlet-triplet gaps in atoms and molecules employing the third-order spin-flip algebraic diagrammatic construction scheme for the polarization propagator. J Chem Phys 2016; 145:084102. [DOI: 10.1063/1.4961298] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Daniel Lefrancois
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Dirk R. Rehn
- Departments of Physics, Chemistry and Biology, Linköping University, SE-581 83 Linköping, Sweden
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| |
Collapse
|
25
|
Sinha Ray S, Ghosh A, Chattopadhyay S, Chaudhuri RK. Taming the Electronic Structure of Diradicals through the Window of Computationally Cost Effective Multireference Perturbation Theory. J Phys Chem A 2016; 120:5897-916. [DOI: 10.1021/acs.jpca.6b03211] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Suvonil Sinha Ray
- Department of Chemistry, Indian Institute of Engineering Science and Technology, Shibpur, Howrah 711103, India
| | - Anirban Ghosh
- Department of Chemistry, Indian Institute of Engineering Science and Technology, Shibpur, Howrah 711103, India
| | - Sudip Chattopadhyay
- Department of Chemistry, Indian Institute of Engineering Science and Technology, Shibpur, Howrah 711103, India
| | - Rajat K. Chaudhuri
- Theoretical Physics, Indian Institute of Astrophysics, Bangalore 560034, India
| |
Collapse
|
26
|
Mullinax JW, Sokolov AY, Schaefer HF. Can density cumulant functional theory describe static correlation effects? J Chem Theory Comput 2016; 11:2487-95. [PMID: 26575548 DOI: 10.1021/acs.jctc.5b00346] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We evaluate the performance of density cumulant functional theory (DCT) for capturing static correlation effects. In particular, we examine systems with significant multideterminant character of the electronic wave function, such as the beryllium dimer, diatomic carbon, m-benzyne, 2,6-pyridyne, twisted ethylene, as well as the barrier for double-bond migration in cyclobutadiene. We compute molecular properties of these systems using the ODC-12 and DC-12 variants of DCT and compare these results to multireference configuration interaction and multireference coupled-cluster theories, as well as single-reference coupled-cluster theory with single, double (CCSD), and perturbative triple excitations [CCSD(T)]. For all systems the DCT methods show intermediate performance between that of CCSD and CCSD(T), with significant improvement over the former method. In particular, for the beryllium dimer, m-benzyne, and 2,6-pyridyne, the ODC-12 method along with CCSD(T) correctly predict the global minimum structures, while CCSD predictions fail qualitatively, underestimating the multireference effects. Our results suggest that the DC-12 and ODC-12 methods are capable of describing emerging static correlation effects but should be used cautiously when highly accurate results are required. Conveniently, the appearance of multireference effects in DCT can be diagnosed by analyzing the DCT natural orbital occupations, which are readily available at the end of the energy computation.
Collapse
Affiliation(s)
- J Wayne Mullinax
- Center for Computational Quantum Chemistry, University of Georgia , Athens, Georgia 30602, United States
| | - Alexander Yu Sokolov
- Center for Computational Quantum Chemistry, University of Georgia , Athens, Georgia 30602, United States.,Department of Chemistry, Princeton University , Princeton, New Jersey 08544, United States
| | - Henry F Schaefer
- Center for Computational Quantum Chemistry, University of Georgia , Athens, Georgia 30602, United States
| |
Collapse
|
27
|
Giner E, David G, Scemama A, Malrieu JP. A simple approach to the state-specific MR-CC using the intermediate Hamiltonian formalism. J Chem Phys 2016; 144:064101. [DOI: 10.1063/1.4940781] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
|
28
|
Aprà E, Kowalski K. Implementation of High-Order Multireference Coupled-Cluster Methods on Intel Many Integrated Core Architecture. J Chem Theory Comput 2016; 12:1129-38. [DOI: 10.1021/acs.jctc.5b00957] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- E. Aprà
- William R. Wiley Environmental
Molecular Sciences Laboratory, Battelle, Pacific Northwest National Laboratory, K8-91, P.O. Box 999, Richland, Washington 99352, United States
| | - K. 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
|
29
|
Qiu Y, Wu CH, Schaefer III HF, Allen WD, Agarwal J. σ Bond activation through tunneling: formation of the boron hydride cations BHn+ (n = 2, 4, 6). Phys Chem Chem Phys 2016; 18:4063-70. [DOI: 10.1039/c5cp05505a] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The network of H2 additions to B+ and subsequent insertion reactions serve as a tractable model for hydrogen storage in elementary boron-containing compounds.
Collapse
Affiliation(s)
- Yudong Qiu
- Center for Computational Quantum Chemistry
- University of Georgia
- Athens
- Georgia
| | - Chia-Hua Wu
- Center for Computational Quantum Chemistry
- University of Georgia
- Athens
- Georgia
| | | | - Wesley D. Allen
- Center for Computational Quantum Chemistry
- University of Georgia
- Athens
- Georgia
| | - Jay Agarwal
- Center for Computational Quantum Chemistry
- University of Georgia
- Athens
- Georgia
| |
Collapse
|
30
|
Xu E, Zhao D, Li S. Multireference Second Order Perturbation Theory with a Simplified Treatment of Dynamical Correlation. J Chem Theory Comput 2015; 11:4634-43. [DOI: 10.1021/acs.jctc.5b00495] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Enhua Xu
- School of Chemistry and Chemical
Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of
Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing, 210093, P. R. China
| | - Dongbo Zhao
- School of Chemistry and Chemical
Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of
Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing, 210093, P. R. China
| | - Shuhua Li
- School of Chemistry and Chemical
Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of
Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing, 210093, P. R. China
| |
Collapse
|
31
|
Bytautas L, Scuseria GE, Ruedenberg K. Seniority number description of potential energy surfaces: Symmetric dissociation of water, N2, C2, and Be2. J Chem Phys 2015; 143:094105. [DOI: 10.1063/1.4929904] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Laimutis Bytautas
- Department of Chemistry, Galveston College, 4015 Ave. Q, Galveston, Texas 77550, USA
| | - Gustavo E. Scuseria
- Department of Chemistry, Rice University, Houston, Texas 77005, USA
- Chemistry Department, Faculty of Science, King Abdulaziz University, Jeddah 21589, Saudi Arabia
| | - Klaus Ruedenberg
- Department of Chemistry, Iowa State University, Ames, Iowa 50011, USA
| |
Collapse
|
32
|
Demel O, Pittner J, Neese F. A Local Pair Natural Orbital-Based Multireference Mukherjee’s Coupled Cluster Method. J Chem Theory Comput 2015; 11:3104-14. [DOI: 10.1021/acs.jctc.5b00334] [Citation(s) in RCA: 50] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Ondřej Demel
- J.
Heyrovský Institute of Physical Chemistry, v.v.i., Academy of Sciences of the Czech Republic, Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Jiří Pittner
- J.
Heyrovský Institute of Physical Chemistry, v.v.i., Academy of Sciences of the Czech Republic, Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Frank Neese
- Max Planck Institute
of Chemical Energy Conversion, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
33
|
Banik S, Ravichandran L, Brabec J, Hubač I, Kowalski K, Pittner J. Iterative universal state selective correction for the Brillouin-Wigner multireference coupled-cluster theory. J Chem Phys 2015; 142:114106. [PMID: 25796230 DOI: 10.1063/1.4914311] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
As a further development of the previously introduced a posteriori Universal State-Selective (USS) corrections [K. Kowalski, J. Chem. Phys. 134, 194107 (2011); J. Brabec et al., ibid. 136, 124102 (2012)], we suggest an iterative form of the USS correction by means of correcting effective Hamiltonian matrix elements. We also formulate USS corrections via the left Bloch equations. The convergence of the USS corrections with excitation level towards the full configuration interaction (FCI) limit is also investigated. Various forms of the USS and simplified diagonal USS corrections at the singles and doubles and perturbative triple levels are numerically assessed on several model systems and on the ozone and tetramethyleneethane molecules. It is shown that the iterative USS correction can successfully replace the previously developed a posteriori Brillouin-Wigner coupled cluster size-extensivity correction, while it is not sensitive to intruder states and performs well also in other cases when the a posteriori one fails, like, e.g., for the asymmetric vibration mode of ozone.
Collapse
Affiliation(s)
- Subrata Banik
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, CZ-18223 Prague 8, Czech Republic
| | - Lalitha Ravichandran
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, CZ-18223 Prague 8, Czech Republic
| | - Jiří Brabec
- Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Ivan Hubač
- Department of Chemical Physics, Faculty of Mathematics and Physics, Comenius University, 84215 Bratislava, Slovakia
| | - Karol Kowalski
- William R. Wiley Environmental Molecular Sciences Laboratory, Battelle, Pacific Northwest National Laboratory, K8-91, P.O. Box 999, Richland, Washington 99352, USA
| | - Jiří Pittner
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, CZ-18223 Prague 8, Czech Republic
| |
Collapse
|
34
|
Chattopadhyay S, Chaudhuri RK, Mahapatra US. State-specific multireference perturbation theory with improved virtual orbitals: Taming the ground state of F2, Be2,and N2. J Comput Chem 2015; 36:907-25. [DOI: 10.1002/jcc.23873] [Citation(s) in RCA: 32] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2014] [Revised: 01/09/2015] [Accepted: 01/11/2015] [Indexed: 11/07/2022]
Affiliation(s)
- Sudip Chattopadhyay
- Department of Chemistry; Indian Institute of Engineering Science and Technology; Shibpur Howrah 711103 India
| | | | - Uttam Sinha Mahapatra
- Department of Physics; Maulana Azad College; 8 Rafi Ahmed Kidwai Road Kolkata 700013 India
| |
Collapse
|
35
|
Abstract
A method is suggested to build simple multiconfigurational wave functions specified uniquely by an energy cutoff Λ. These are constructed from a model space containing determinants with energy relative to that of the most stable determinant no greater than Λ. The resulting Λ-CI wave function is adaptive, being able to represent both single-reference and multireference electronic states. We also consider a more compact wave function parameterization (Λ+SD-CI), which is based on a small Λ-CI reference and adds a selection of all the singly and doubly excited determinants generated from it. We report two heuristic algorithms to build Λ-CI wave functions. The first is based on an approximate prescreening of the full configuration interaction space, while the second performs a breadth-first search coupled with pruning. The Λ-CI and Λ+SD-CI approaches are used to compute the dissociation curve of N2 and the potential energy curves for the first three singlet states of C2. Special attention is paid to the issue of energy discontinuities caused by changes in the size of the Λ-CI wave function along the potential energy curve. This problem is shown to be solvable by smoothing the matrix elements of the Hamiltonian. Our last example, involving the Cu2O2(2+) core, illustrates an alternative use of the Λ-CI method: as a tool to both estimate the multireference character of a wave function and to create a compact model space to be used in subsequent high-level multireference coupled cluster computations.
Collapse
Affiliation(s)
- Francesco A Evangelista
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| |
Collapse
|
36
|
Boguslawski K, Tecmer P, Bultinck P, De Baerdemacker S, Van Neck D, Ayers PW. Nonvariational Orbital Optimization Techniques for the AP1roG Wave Function. J Chem Theory Comput 2014; 10:4873-82. [DOI: 10.1021/ct500759q] [Citation(s) in RCA: 45] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Affiliation(s)
- Katharina Boguslawski
- Department
of Chemistry and Chemical Biology, McMaster University, 1280 Main
Street West, Hamilton, L8S 4M1 Ontario, Canada
| | - Paweł Tecmer
- Department
of Chemistry and Chemical Biology, McMaster University, 1280 Main
Street West, Hamilton, L8S 4M1 Ontario, Canada
| | - Patrick Bultinck
- Department
of Inorganic and Physical Chemistry, Ghent University, Krijgslaan
281 (S3), 9000 Gent, East Flanders, Belgium
| | - Stijn De Baerdemacker
- Center
for Molecular Modelling, Ghent University, Technologiepark 903, 9052 Gent, East
Flanders, Belgium
| | - Dimitri Van Neck
- Center
for Molecular Modelling, Ghent University, Technologiepark 903, 9052 Gent, East
Flanders, Belgium
| | - Paul W. Ayers
- Department
of Chemistry and Chemical Biology, McMaster University, 1280 Main
Street West, Hamilton, L8S 4M1 Ontario, Canada
| |
Collapse
|
37
|
Aspects of size extensivity in unitary group adapted multi-reference coupled cluster theories: the role of cumulant decomposition of spin-free reduced density matrices. Theor Chem Acc 2014. [DOI: 10.1007/s00214-014-1522-5] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
38
|
Demel O, Kedžuch S, Noga J, Pittner J. Perturbative triples correction for explicitly correlated Mukherjee’s state-specific coupled cluster method. Mol Phys 2013. [DOI: 10.1080/00268976.2013.809488] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
39
|
Bhaskaran-Nair K, Kowalski K. Bridging single and multireference coupled cluster theories with universal state selective formalism. J Chem Phys 2013; 138:204114. [DOI: 10.1063/1.4806768] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
40
|
|
41
|
Bhaskaran-Nair K, Ma W, Krishnamoorthy S, Villa O, van Dam HJJ, Aprà E, Kowalski K. Noniterative Multireference Coupled Cluster Methods on Heterogeneous CPU–GPU Systems. J Chem Theory Comput 2013; 9:1949-57. [DOI: 10.1021/ct301130u] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Kiran Bhaskaran-Nair
- William R. Wiley Environmental Molecular
Sciences Laboratory,
Battelle, Pacific Northwest National Laboratory, K8-91, P.O. Box 999,
Richland, Washington 99352, United States
| | - Wenjing Ma
- William R. Wiley Environmental Molecular
Sciences Laboratory,
Battelle, Pacific Northwest National Laboratory, K8-91, P.O. Box 999,
Richland, Washington 99352, United States
| | - Sriram Krishnamoorthy
- William R. Wiley Environmental Molecular
Sciences Laboratory,
Battelle, Pacific Northwest National Laboratory, K8-91, P.O. Box 999,
Richland, Washington 99352, United States
| | - Oreste Villa
- William R. Wiley Environmental Molecular
Sciences Laboratory,
Battelle, Pacific Northwest National Laboratory, K8-91, P.O. Box 999,
Richland, Washington 99352, United States
| | - Hubertus J. J. van Dam
- William R. Wiley Environmental Molecular
Sciences Laboratory,
Battelle, Pacific Northwest National Laboratory, K8-91, P.O. Box 999,
Richland, Washington 99352, United States
| | - Edoardo Aprà
- 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
|
42
|
Yang KR, Jalan A, Green WH, Truhlar DG. Which Ab Initio Wave Function Methods Are Adequate for Quantitative Calculations of the Energies of Biradicals? The Performance of Coupled-Cluster and Multi-Reference Methods Along a Single-Bond Dissociation Coordinate. J Chem Theory Comput 2012; 9:418-31. [DOI: 10.1021/ct3009528] [Citation(s) in RCA: 33] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Ke R. Yang
- Department of Chemistry and
Supercomputing Institute, University of Minnesota, 207 Pleasant Street S. E., Minneapolis, Minnesota 55455-0431, United
States
| | - Amrit Jalan
- Department of Chemical
Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts
02139, United States
| | - William H. Green
- Department of Chemical
Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts
02139, United States
| | - Donald G. Truhlar
- Department of Chemistry and
Supercomputing Institute, University of Minnesota, 207 Pleasant Street S. E., Minneapolis, Minnesota 55455-0431, United
States
| |
Collapse
|
43
|
Melnichuk A, Bartlett RJ. Relaxed active space: Fixing tailored-CC with high order coupled cluster. I. J Chem Phys 2012; 137:214103. [DOI: 10.1063/1.4767900] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
44
|
Shen J, Piecuch P. Merging Active-Space and Renormalized Coupled-Cluster Methods via the CC(P;Q) Formalism, with Benchmark Calculations for Singlet–Triplet Gaps in Biradical Systems. J Chem Theory Comput 2012; 8:4968-88. [DOI: 10.1021/ct300762m] [Citation(s) in RCA: 46] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jun Shen
- Department of Chemistry, Michigan State University,
East Lansing, Michigan 48824, United States
| | - Piotr Piecuch
- Department of Chemistry, Michigan State University,
East Lansing, Michigan 48824, United States
| |
Collapse
|
45
|
Köhn A, Hanauer M, Mück LA, Jagau TC, Gauss J. State-specific multireference coupled-cluster theory. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2012. [DOI: 10.1002/wcms.1120] [Citation(s) in RCA: 83] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
|
46
|
Bhaskaran-Nair K, Brabec J, Aprà E, van Dam HJJ, Pittner J, Kowalski K. Implementation of the multireference Brillouin-Wigner and Mukherjee's coupled cluster methods with non-iterative triple excitations utilizing reference-level parallelism. J Chem Phys 2012; 137:094112. [DOI: 10.1063/1.4747698] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
47
|
Shen J, Piecuch P. Combining active-space coupled-cluster methods with moment energy corrections via the CC(P;Q) methodology, with benchmark calculations for biradical transition states. J Chem Phys 2012; 136:144104. [PMID: 22502498 DOI: 10.1063/1.3700802] [Citation(s) in RCA: 52] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We have recently suggested the CC(P;Q) methodology that can correct energies obtained in the active-space coupled-cluster (CC) or equation-of-motion (EOM) CC calculations, which recover much of the nondynamical and some dynamical electron correlation effects, for the higher-order, mostly dynamical, correlations missing in the active-space CC/EOMCC considerations. It is shown that one can greatly improve the description of biradical transition states, both in terms of the resulting energy barriers and total energies, by combining the CC approach with singles, doubles, and active-space triples, termed CCSDt, with the CC(P;Q)-style correction due to missing triple excitations defining the CC(t;3) approximation.
Collapse
Affiliation(s)
- Jun Shen
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | | |
Collapse
|
48
|
Jagau TC, Gauss J. Linear-response theory for Mukherjee's multireference coupled-cluster method: Excitation energies. J Chem Phys 2012; 137:044116. [DOI: 10.1063/1.4734309] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
49
|
Jagau TC, Gauss J. Linear-response theory for Mukherjee's multireference coupled-cluster method: Static and dynamic polarizabilities. J Chem Phys 2012; 137:044115. [DOI: 10.1063/1.4734308] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
50
|
Maitra R, Sinha D, Mukherjee D. Unitary group adapted state-specific multi-reference coupled cluster theory: Formulation and pilot numerical applications. J Chem Phys 2012; 137:024105. [DOI: 10.1063/1.4731341] [Citation(s) in RCA: 33] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|