1
|
Chen T, Zhang C, Cheng L, Ng KB, Malbrunot-Ettenauer S, Flambaum VV, Lasner Z, Doyle JM, Yu P, Conn CJ, Zhang C, Hutzler NR, Jayich AM, Augenbraun B, DeMille D. Relativistic Exact Two-Component Coupled-Cluster Study of Molecular Sensitivity Factors for Nuclear Schiff Moments. J Phys Chem A 2024. [PMID: 39047199 DOI: 10.1021/acs.jpca.4c02640] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/27/2024]
Abstract
Relativistic exact two-component coupled-cluster calculations of molecular sensitivity factors for nuclear Schiff moments (NSMs) are reported. We focus on molecules containing heavy nuclei, especially octupole-deformed nuclei. Analytic relativistic coupled-cluster gradient techniques are used and serve as useful tools for identifying candidate molecules that sensitively probe for physics beyond the Standard Model in the hadronic sector. Notably, these tools enable straightforward "black-box" calculations. Two competing chemical mechanisms that contribute to the NSM are analyzed, illuminating the physics of ligand effects on NSM sensitivity factors.
Collapse
Affiliation(s)
- Tianxiang Chen
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Chaoqun Zhang
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Lan Cheng
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Kia Boon Ng
- TRIUMF, 4004 Wesbrook Mall, Vancouver, BC V6T 2A3, Canada
| | - Stephan Malbrunot-Ettenauer
- TRIUMF, 4004 Wesbrook Mall, Vancouver, BC V6T 2A3, Canada
- Department of Physics, University of Toronto, Toronto M5S 1A7, Canada
| | - Victor V Flambaum
- School of Physics, University of New South Wales, Sydney 2052, Australia
| | - Zack Lasner
- Department of Physics, Harvard University, Cambridge, Massachusetts 02138, United States
- Harvard-MIT Center for Ultracold Atoms, Cambridge, Massachusetts 02138, United States
| | - John M Doyle
- Department of Physics, Harvard University, Cambridge, Massachusetts 02138, United States
- Harvard-MIT Center for Ultracold Atoms, Cambridge, Massachusetts 02138, United States
| | - Phelan Yu
- Division of Physics, Mathematics, and Astronomy, California Institute of Technology, Pasadena, California 91125, United States
| | - Chandler J Conn
- Division of Physics, Mathematics, and Astronomy, California Institute of Technology, Pasadena, California 91125, United States
| | - Chi Zhang
- Division of Physics, Mathematics, and Astronomy, California Institute of Technology, Pasadena, California 91125, United States
| | - Nicholas R Hutzler
- Division of Physics, Mathematics, and Astronomy, California Institute of Technology, Pasadena, California 91125, United States
| | - Andrew M Jayich
- Department of Physics, University of California, Santa Barbara, California 93106, United States
| | - Benjamin Augenbraun
- Department of Chemistry, Williams College, 47 Lab Campus Drive, Williamstown, Massachusetts 01267, United States
| | - David DeMille
- Department of Physics, University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
2
|
Sun Y. Accelerating Density Matrix Embedding with Stochastic Density Fitting Theory: An Application to Hydrogen Bonded Clusters. J Chem Theory Comput 2024. [PMID: 39029091 DOI: 10.1021/acs.jctc.4c00489] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/21/2024]
Abstract
In this work, we demonstrate how using semistochastic density fitting (ss-DF) can accelerate self-consistent density matrix embedding theory (DMET) calculations by reducing the number of auxiliary orbitals in the three-indexed DF integrals. This reduction results in significant time savings when building the Hartree-Fock (HF) Coulomb and Exchange Matrices and in transforming integrals from the atomic orbital (AO) basis to the embedding orbital (EO) basis. We apply ss-DF to a range of hydrogen-bonded clusters to showcase its effectiveness. First, we examine how the amount of deterministic space impacts the quality of the calculation in a (H2O)10 cluster. Next, we test the computational efficiency of ss-DF compared to deterministic DF (d-DF) in water clusters containing 6-30 water molecules using a triple-ζ basis set. Finally, we perform numerical structural optimizations on water and hydrogen fluoride clusters, revealing that DMET can recover weak interactions using a back-transformed energy formula. This work demonstrates the potential of using stochastic resolution of identity in quantum embedding theories and highlights its capability to recover weak interactions effectively.
Collapse
Affiliation(s)
- Yi Sun
- Department of Chemistry, Chicago Center for Theoretical Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
3
|
Grazioli L, Schleicher LT, Stopkowicz S, Gauss J. Theoretical prediction of closed-shell paramagnetism for scandium and yttrium hydride. J Comput Chem 2024; 45:1215-1223. [PMID: 38334014 DOI: 10.1002/jcc.27305] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2023] [Revised: 12/19/2023] [Accepted: 12/23/2023] [Indexed: 02/10/2024]
Abstract
Following chemical intuition, one would expect that all closed-shell molecules are diamagnetic. However, it is known that this is not the case for some second-row hydrides with low-lying unoccupied π orbitals due to an unquenching of the total angular momentum in the presence of an external magnetic field. In this article, the transition-metal hydrides ScH and YH are investigated, assuming a similar unquenching effect involving low-lying unoccupied π and δ orbitals formed from the metal d orbitals rather than the p orbitals. We are comparing results obtained with various quantum-chemical methods (HF, CCSD, CCSD(T), CCSDT) and basis sets. The obtained positive values for the magnetizabilities clearly indicate paramagnetic behavior. Vibrational effects on the magnetizability tensor are also considered, but these effects are small and do not change the overall conclusion that both ScH and YH are further examples for closed-shell paramagnetism.
Collapse
Affiliation(s)
- Laura Grazioli
- Department Chemie, Johannes Gutenberg-Universität Mainz, Mainz, Germany
| | - Luca T Schleicher
- Department Chemie, Johannes Gutenberg-Universität Mainz, Mainz, Germany
| | - Stella Stopkowicz
- Fachrichtung Chemie, Universität des Saarlandes, Saarbrücken, Germany
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, Oslo, Norway
| | - Jürgen Gauss
- Department Chemie, Johannes Gutenberg-Universität Mainz, Mainz, Germany
| |
Collapse
|
4
|
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
|
5
|
Rutskoy B, Ozerov G, Bezrukov D. The Role of Bond Functions in Describing Intermolecular Electron Correlation for Van der Waals Dimers: A Study of (CH 4) 2 and Ne 2. Int J Mol Sci 2024; 25:1472. [PMID: 38338750 PMCID: PMC10855067 DOI: 10.3390/ijms25031472] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2023] [Revised: 01/19/2024] [Accepted: 01/22/2024] [Indexed: 02/12/2024] Open
Abstract
We present a study of the intermolecular interactions in van der Waals complexes of methane and neon dimers within the framework of the CCSD method. This approach was implemented and applied to calculate and examine the behavior of the contracted two-particle reduced density matrix (2-RDM). It was demonstrated that the region near the minimum of the two-particle density matrix correlation part, corresponding to the primary bulk of the Coulomb hole contribution, exerts a significant influence on the dispersion interaction energetics of the studied systems. As a result, the bond functions approach was applied to improve the convergence performance for the intermolecular correlation energy results with respect to the size of the atomic basis. For this, substantial acceleration was achieved by introducing an auxiliary basis of bond functions centered on the minima of the 2-RDM. For both methane and neon dimers, this general conclusion was confirmed with a series of CCSD calculations for the 2-RDM and the correlation energies.
Collapse
Affiliation(s)
- Bogdan Rutskoy
- National Research Centre “Kurchatov Institute”, Moscow 123182, Russia;
- Institute of Nuclear Physics and Technology, National Research Nuclear University “MEPhI” (Moscow Engineering Physics Institute), Moscow 115409, Russia
- Chemistry Department, M.V. Lomonosov Moscow State University, Moscow 119991, Russia;
| | - Georgiy Ozerov
- Chemistry Department, M.V. Lomonosov Moscow State University, Moscow 119991, Russia;
| | - Dmitry Bezrukov
- Chemistry Department, M.V. Lomonosov Moscow State University, Moscow 119991, Russia;
| |
Collapse
|
6
|
Damour Y, Quintero-Monsebaiz R, Caffarel M, Jacquemin D, Kossoski F, Scemama A, Loos PF. Ground- and Excited-State Dipole Moments and Oscillator Strengths of Full Configuration Interaction Quality. J Chem Theory Comput 2023; 19:221-234. [PMID: 36548519 DOI: 10.1021/acs.jctc.2c01111] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
We report ground- and excited-state dipole moments and oscillator strengths (computed in different "gauges" or representations) of full configuration interaction (FCI) quality using the selected configuration interaction method known as Configuration Interaction using a Perturbative Selection made Iteratively (CIPSI). Thanks to a set encompassing 35 ground- and excited-state properties computed in 11 small molecules, the present near-FCI estimates allow us to assess the accuracy of high-order coupled-cluster (CC) calculations including up to quadruple excitations. In particular, we show that incrementing the excitation degree of the CC expansion (from CC with singles and doubles (CCSD) to CC with singles, doubles, and triples (CCSDT) or from CCSDT to CC with singles, doubles, triples, and quadruples (CCSDTQ)) reduces the average error with respect to the near-FCI reference values by approximately 1 order of magnitude.
Collapse
Affiliation(s)
- Yann Damour
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Raúl Quintero-Monsebaiz
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Michel Caffarel
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Denis Jacquemin
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France.,Institut Universitaire de France (IUF), F-75005 Paris, France
| | - Fábris Kossoski
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Anthony Scemama
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| |
Collapse
|
7
|
Lesiuk M. When Gold Is Not Enough: Platinum Standard of Quantum Chemistry with N7 Cost. J Chem Theory Comput 2022; 18:6537-6556. [PMID: 36314739 PMCID: PMC9648181 DOI: 10.1021/acs.jctc.2c00460] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
![]()
In this paper, we extend the rank-reduced coupled-cluster
formalism
to the calculation of non-iterative energy corrections due to quadruple
excitations. There are two major components of the proposed formalism.
The first is an approximate compression of the quadruple excitation
amplitudes using the Tucker format. The second is a modified functional
used for the evaluation of the corrections which gives exactly the
same results for the exact amplitudes, but is less susceptible to
errors resulting from the aforementioned compression. We show, both
theoretically and numerically, that the computational cost of the
proposed method scales as the seventh power of the system size. Using
reference results for a set of small molecules, the method is calibrated
to deliver relative accuracy of a few percent in energy corrections.
To illustrate the potential of the theory, we calculate the isomerization
energy of ortho/meta benzyne (C6H4) and the barrier height for the Cope rearrangement
in bullvalene (C10H10). The method retains a
near-black-box nature of the conventional coupled-cluster formalism
and depends on only one additional parameter that controls the accuracy.
Collapse
Affiliation(s)
- Michał Lesiuk
- Faculty of Chemistry, University of Warsaw, Pasteura 1, Warsaw 02-093, Poland
| |
Collapse
|
8
|
Zheng X, Zhang C, Liu J, Cheng L. Geometry Optimizations with Spinor-Based Relativistic Coupled-Cluster Theory. J Chem Phys 2022; 156:151101. [DOI: 10.1063/5.0086281] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Development of analytic gradients for relativistic coupled-cluster singles and doubles augmented with a non-iterative triples [CCSD(T)] method using an all-electron exact two-component Hamiltonian with atomic mean-field spin-orbit integrals (X2CAMF) is reported. This enables efficient CC geometry optimizations with spin-orbit coupling included in orbitals. The applicability of the implementation is demonstrated using benchmark X2CAMF-CCSD(T) calculations of equilibrium structures and harmonic vibrational frequencies for methyl halides, CH3X, X=Br, I, At, as well as calculations of rotational constants and infrared spectrum for RaSH+, a radioactive molecular ion of interest to spectroscopic study.
Collapse
Affiliation(s)
- Xuechen Zheng
- Johns Hopkins University Department of Chemistry, United States of America
| | - Chaoqun Zhang
- Johns Hopkins University Department of Chemistry, United States of America
| | - Junzi Liu
- Chemistry, Johns Hopkins University Department of Chemistry, United States of America
| | - Lan Cheng
- Chemistry, Johns Hopkins University Department of Chemistry, United States of America
| |
Collapse
|
9
|
Alagöz Y, Ünal A, Bozkaya U. Efficient implementations of the symmetric and asymmetric triple excitation corrections for the orbital-optimized coupled-cluster doubles method with the density-fitting approximation. J Chem Phys 2021; 155:114104. [PMID: 34551547 DOI: 10.1063/5.0061351] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Efficient implementations of the symmetric and asymmetric triple excitation corrections for the orbital-optimized coupled-cluster doubles (OCCD) method with the density-fitting approach, denoted by DF-OCCD(T) and DF-OCCD(T)Λ, are presented. The computational cost of the DF-OCCD(T) method is compared with that of the conventional OCCD(T). In the conventional OCCD(T) and OCCD(T)Λ methods, one needs to perform four-index integral transformations at each coupled-cluster doubles iterations, which limits its applications to large chemical systems. Our results demonstrate that DF-OCCD(T) provides dramatically lower computational costs compared to OCCD(T), and there are more than 68-fold reductions in the computational time for the C5H12 molecule with the cc-pVTZ basis set. Our results show that the DF-OCCD(T) and DF-OCCD(T)Λ methods are very helpful for the study of single bond-breaking problems. Performances of the DF-OCCD(T) and DF-OCCD(T)Λ methods are noticeably better than that of the coupled-cluster singles and doubles with perturbative triples [CCSD(T)] method for the potential energy surfaces of the molecules considered. Specifically, the DF-OCCD(T)Λ method provides dramatic improvements upon CCSD(T), and there are 8-14-fold reductions in nonparallelity errors. Overall, we conclude that the DF-OCCD(T)Λ method is very promising for the study of challenging chemical systems, where the CCSD(T) fails.
Collapse
Affiliation(s)
- Yavuz Alagöz
- Department of Chemistry, Hacettepe University, Ankara 06800, Turkey
| | - Aslı Ünal
- Department of Chemistry, Hacettepe University, Ankara 06800, Turkey
| | - Uğur Bozkaya
- Department of Chemistry, Hacettepe University, Ankara 06800, Turkey
| |
Collapse
|
10
|
Liu J, Zheng X, Asthana A, Zhang C, Cheng L. Analytic evaluation of energy first derivatives for spin-orbit coupled-cluster singles and doubles augmented with noniterative triples method: General formulation and an implementation for first-order properties. J Chem Phys 2021; 154:064110. [PMID: 33588557 DOI: 10.1063/5.0038779] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A formulation of analytic energy first derivatives for the coupled-cluster singles and doubles augmented with noniterative triples [CCSD(T)] method with spin-orbit coupling included at the orbital level and an implementation for evaluation of first-order properties are reported. The standard density-matrix formulation for analytic CC gradient theory adapted to complex algebra has been used. The orbital-relaxation contributions from frozen core, occupied, virtual, and frozen virtual orbitals to analytic spin-orbit CCSD(T) gradients are fully taken into account and treated efficiently, which is of importance to calculations of heavy elements. Benchmark calculations of first-order properties including dipole moments and electric-field gradients using the corresponding exact two-component property integrals are presented for heavy-element containing molecules to demonstrate the applicability and usefulness of the present analytic scheme.
Collapse
Affiliation(s)
- Junzi Liu
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Xuechen Zheng
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Ayush Asthana
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Chaoqun Zhang
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Lan Cheng
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, USA
| |
Collapse
|
11
|
Bozkaya U, Ünal A, Alagöz Y. Energy and analytic gradients for the orbital-optimized coupled-cluster doubles method with the density-fitting approximation: An efficient implementation. J Chem Phys 2020; 153:244115. [PMID: 33380091 DOI: 10.1063/5.0035811] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Efficient implementations of the orbital-optimized coupled-cluster doubles (or simply "optimized CCD," OCCD, for short) method and its analytic energy gradients with the density-fitting (DF) approach, denoted by DF-OCCD, are presented. In addition to the DF approach, the Cholesky-decomposed variant (CD-OCCD) is also implemented for energy computations. The computational cost of the DF-OCCD method (available in a plugin version of the DFOCC module of PSI4) is compared with that of the conventional OCCD (from the Q-CHEM package). The OCCD computations were performed with the Q-CHEM package in which OCCD are denoted by OD. In the conventional OCCD method, one needs to perform four-index integral transformations at each of the CCD iterations, which limits its applications to large chemical systems. Our results demonstrate that DF-OCCD provides dramatically lower computational costs compared to OCCD, and there are almost eightfold reductions in the computational time for the C6H14 molecule with the cc-pVTZ basis set. For open-shell geometries, interaction energies, and hydrogen transfer reactions, DF-OCCD provides significant improvements upon DF-CCD. Furthermore, the performance of the DF-OCCD method is substantially better for harmonic vibrational frequencies in the case of symmetry-breaking problems. Moreover, several factors make DF-OCCD more attractive compared to CCSD: (1) for DF-OCCD, there is no need for orbital relaxation contributions in analytic gradient computations; (2) active spaces can readily be incorporated into DF-OCCD; (3) DF-OCCD provides accurate vibrational frequencies when symmetry-breaking problems are observed; (4) in its response function, DF-OCCD avoids artificial poles; hence, excited-state molecular properties can be computed via linear response theory; and (5) symmetric and asymmetric triples corrections based on DF-OCCD [DF-OCCD(T)] have a significantly better performance in near degeneracy regions.
Collapse
Affiliation(s)
- Uğur Bozkaya
- Department of Chemistry, Hacettepe University, Ankara 06800, Turkey
| | - Aslı Ünal
- Department of Chemistry, Hacettepe University, Ankara 06800, Turkey
| | - Yavuz Alagöz
- Department of Chemistry, Hacettepe University, Ankara 06800, Turkey
| |
Collapse
|
12
|
Eriksen JJ, Stopkowicz S, Jagau TC, Helgaker T. Foreword: Prof. Gauss Festschrift. Mol Phys 2020. [DOI: 10.1080/00268976.2020.1817247] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Affiliation(s)
| | - Stella Stopkowicz
- Department Chemie, Johannes Gutenberg-Universität Mainz, Mainz, Germany
| | - Thomas-C. Jagau
- Departement Chemie, Katholieke Universiteit Leuven, Leuven, Belgium
| | - Trygve Helgaker
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, Oslo, Norway
| |
Collapse
|
13
|
Matthews DA. Analytic Gradients of Approximate Coupled Cluster Methods with Quadruple Excitations. J Chem Theory Comput 2020; 16:6195-6206. [DOI: 10.1021/acs.jctc.0c00522] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
14
|
Ye HZ, Tran HK, Van Voorhis T. Bootstrap Embedding For Large Molecular Systems. J Chem Theory Comput 2020; 16:5035-5046. [PMID: 32589842 DOI: 10.1021/acs.jctc.0c00438] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
Recent developments in quantum embedding theories have provided attractive approaches to correlated calculations for large systems. In this work, we extend our previous work [J. Chem. Theory Comput. 2019, 15, 4497-4506; J. Phys. Chem. Lett. 2019, 10, 6368-6374] on bootstrap embedding (BE) to enable correlated ab initio calculations at the coupled cluster with singles and doubles (CCSD) level for large molecules. We introduce several new algorithmic developments that significantly reduce the computational cost of BE, while maintaining its accuracy. The resulting implementation scales as O(N3) for the integral transform and O(N) for the CCSD calculation. Numerical results on a series of conjugated molecules suggest that BE with reasonably sized fragments can recover more than 99.5% of the total correlation energy of a full CCSD calculation, while the required computational resources (time and storage) compare favorably to one popular local correlation scheme: domain localized pair natural orbital (DLPNO). The largest BE calculation in this work involves ∼2900 basis functions and can be performed on a single node with 16 CPU cores and 64 GB of memory in a few days. We anticipate that these developments represent an important step toward the application of BE to solve practical problems.
Collapse
Affiliation(s)
- Hong-Zhou Ye
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Henry K Tran
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Troy Van Voorhis
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
15
|
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
|
16
|
Warden CE, Smith DGA, Burns LA, Bozkaya U, Sherrill CD. Efficient and automated computation of accurate molecular geometries using focal-point approximations to large-basis coupled-cluster theory. J Chem Phys 2020; 152:124109. [DOI: 10.1063/5.0004863] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Constance E. Warden
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Daniel G. A. Smith
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Lori A. Burns
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Uğur Bozkaya
- Department of Chemistry, Hacettepe University, Ankara 06800, Turkey
| | - C. David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| |
Collapse
|
17
|
Hampe F, Gross N, Stopkowicz S. Full triples contribution in coupled-cluster and equation-of-motion coupled-cluster methods for atoms and molecules in strong magnetic fields. Phys Chem Chem Phys 2020; 22:23522-23529. [DOI: 10.1039/d0cp04169f] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Finite-field EOM-CCSDT: a highly accurate method for the theoretical prediction of excitation energies and electronic spectra in strong magnetic fields.
Collapse
Affiliation(s)
- Florian Hampe
- Department für Chemie, Johannes Gutenberg-Universität Mainz
- D-55099 Mainz
- Germany
- Department of Chemistry, University of Southern California
- Los Angeles
| | - Niklas Gross
- Department für Chemie, Johannes Gutenberg-Universität Mainz
- D-55099 Mainz
- Germany
| | - Stella Stopkowicz
- Department für Chemie, Johannes Gutenberg-Universität Mainz
- D-55099 Mainz
- Germany
| |
Collapse
|
18
|
Wang S, Li C, Evangelista FA. Analytic gradients for the single-reference driven similarity renormalization group second-order perturbation theory. J Chem Phys 2019; 151:044118. [PMID: 31370522 DOI: 10.1063/1.5100175] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Abstract
We derive and implement analytic energy gradients for the single-reference driven similarity renormalization group second-order perturbation theory (DSRG-PT2). The resulting equations possess an asymptotic scaling that is identical to that of the second-order Møller-Plesset perturbation theory (MP2), indicating that the exponential regularizer in the DSRG equations does not introduce formal difficulties in the gradient theory. We apply the DSRG-PT2 method to optimizing the geometries of 15 small molecules. The equilibrium bond lengths computed with DSRG-PT2 are found similar to those of MP2, yielding a mean absolute error of 0.0033 Å and a standard deviation of 0.0045 Å when compared with coupled cluster with singles, doubles, and perturbative triples.
Collapse
Affiliation(s)
- Shuhe Wang
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - 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
|
19
|
Feng X, Epifanovsky E, Gauss J, Krylov AI. Implementation of analytic gradients for CCSD and EOM-CCSD using Cholesky decomposition of the electron-repulsion integrals and their derivatives: Theory and benchmarks. J Chem Phys 2019; 151:014110. [DOI: 10.1063/1.5100022] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Affiliation(s)
- Xintian Feng
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, USA
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Evgeny Epifanovsky
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Jürgen Gauss
- Institut für Physikalische Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, D-55128 Mainz, Germany
| | - Anna I. Krylov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, USA
- Institut für Physikalische Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, D-55128 Mainz, Germany
| |
Collapse
|
20
|
Breidung J, Thiel W. Equilibrium Structures of the Phosphorus Trihalides PF 3 and PCl 3, and the Phosphoranes PH 3F 2, PF 5, PCl 3F 2, and PCl 5. J Phys Chem A 2019; 123:5600-5612. [PMID: 31181157 PMCID: PMC6750832 DOI: 10.1021/acs.jpca.9b04406] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2019] [Revised: 06/07/2019] [Indexed: 11/30/2022]
Abstract
Among the title species, a reliable and accurate equilibrium geometry ( re structure) is available only for PF3, which has been determined experimentally more than 20 years ago. Here, we report accurate re structures for all title molecules, which were obtained using a composite computational approach based on explicitly correlated coupled-cluster theory (CCSD(T)-F12b) in conjunction with a large correlation-consistent basis set (cc-pCVQZ-F12) to take core-valence electron correlation into account. Additional terms were included to correct for the effects of iterative triple excitations (CCSDT), noniterative quadruple excitations (CCSDT(Q)), and scalar relativistic contributions (DKH2-CCSD(T)). The performance of this computational procedure was established through test calculations on selected small molecules (PH, PF, PCl, PH2, PF2, and PH3). For PF3, PCl3, PH3F2, and PF5 sufficiently accurate experimental ground-state rotational constants from the literature were used to determine semiexperimental re structures, which were found to be in excellent agreement with the corresponding best estimates from the current composite approach. The recommended equilibrium structural parameters are for PCl3, re(PCl) = 203.94 pm and θe(ClPCl) = 100.18°; for PH3F2, re(PHeq) = 138.38 pm and re(PFax) = 164.15 pm; for PF5, re(PFeq) = 153.10 pm and re(PFax) = 157.14 pm; for PCl3F2, re(PCleq) = 200.21 pm and re(PFax) = 159.37 pm; and for PCl5, re(PCleq) = 201.29 pm and re(PClax) = 211.83 pm. The associated uncertainties are estimated to be ±0.10 pm and ±0.10°, respectively.
Collapse
Affiliation(s)
- Jürgen Breidung
- Max-Planck-Institut für
Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Walter Thiel
- Max-Planck-Institut für
Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
21
|
Zhou X, Cao Z, Wang F. Analytical energy gradients for ionized states using equation-of-motion coupled-cluster theory with spin-orbit coupling. J Chem Phys 2019; 150:154114. [PMID: 31005096 DOI: 10.1063/1.5088639] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Spin-orbit coupling (SOC) may have a significant effect on the structure and harmonic frequencies of particularly heavy p-block element compounds. However, reports on analytical energy gradients with SOC are scarce, especially for excited states. In this work, we implemented analytical energy gradients for ionized states using the equation-of-motion coupled-cluster (CC) theory at the CC singles and doubles level (EOM-IP-CCSD) with SOC. Effects of SOC on structure and harmonic frequencies as well as properties for both the ground and some excited states of open-shell compounds with one unpaired electron can be investigated efficiently with the present implementation. A closed-shell reference is required in the calculations, and SOC is included in post-Hartree-Fock treatment. Relativistic effective core potentials are employed in dealing with both scalar relativistic effects and SOC, and we treat perturbations that are even under time reversal in this work. Both time-reversal symmetry and double point group symmetry for D2h * and its subgroups are exploited in the implementation. The method is applicable to states which can be reached by removing one electron from a closed-shell reference state. The results of some open-shell cations indicate the importance of SOC on structures and harmonic frequencies of heavy element compounds.
Collapse
Affiliation(s)
- Xiaojun Zhou
- Institute of Atomic and Molecular Physics, Key Laboratory of High Energy Density Physics and Technology, Ministry of Education, Sichuan University, Chengdu, People's Republic of China
| | - Zhanli Cao
- School of Science, Xi'an University of Posts and Telecommunications, Xi'an 710121, People's Republic of China
| | - Fan Wang
- Institute of Atomic and Molecular Physics, Key Laboratory of High Energy Density Physics and Technology, Ministry of Education, Sichuan University, Chengdu, People's Republic of China
| |
Collapse
|
22
|
Luzon I, Livshits E, Gope K, Baer R, Strasser D. Making Sense of Coulomb Explosion Imaging. J Phys Chem Lett 2019; 10:1361-1367. [PMID: 30840457 DOI: 10.1021/acs.jpclett.9b00576] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
A multifaceted agreement between ab initio theoretical predictions and experimental measurements, including branching ratios, channel-specific kinetic energy release, and three-body momentum correlation spectra, leads to the identification of new mechanisms in Coulomb-explosion (CE) induced two- and three-body breakup processes in methanol. These identified mechanisms include direct nonadiabatic Coulomb explosion responsible for CO bond-breaking, a long-range " inverse harpooning" dominating the production of H2+ + HCOH+, a transient proton migration leading to surprising energy partitioning in three-body fragmentation and other complex dynamics forming products such as H2O+ and H3+. These mechanisms provide general concepts that should be useful for analyzing future time-resolved Coulomb explosion imaging of methanol as well as other molecular systems. These advances are enabled by a combination of recently developed experimental and computational techniques, using weak ultrafast EUV pulses to initiate the CE and a high-level quantum chemistry approach to follow the resulting field-free nonadiabatic molecular dynamics.
Collapse
Affiliation(s)
- Itamar Luzon
- Institute of Chemistry , The Hebrew University of Jerusalem , Jerusalem 91904 , Israel
| | - Ester Livshits
- Fritz Haber Center for Molecular Dynamics and the Institute of Chemistry , The Hebrew University of Jerusalem , Jerusalem 91904 , Israel
| | - Krishnendu Gope
- Institute of Chemistry , The Hebrew University of Jerusalem , Jerusalem 91904 , Israel
| | - Roi Baer
- Fritz Haber Center for Molecular Dynamics and the Institute of Chemistry , The Hebrew University of Jerusalem , Jerusalem 91904 , Israel
| | - Daniel Strasser
- Institute of Chemistry , The Hebrew University of Jerusalem , Jerusalem 91904 , Israel
| |
Collapse
|
23
|
Datta D, Gauss J. Accurate Prediction of Hyperfine Coupling Tensors for Main Group Elements Using a Unitary Group Based Rigorously Spin-Adapted Coupled-Cluster Theory. J Chem Theory Comput 2019; 15:1572-1592. [PMID: 30698956 DOI: 10.1021/acs.jctc.8b01048] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present the development of a perturbative triples correction scheme for the previously reported unitary group based spin-adapted combinatoric open-shell coupled-cluster (CC) singles and doubles (COS-CCSD) approach and report on the applications of the newly developed method, termed "COS-CCSD(T)", to the calculation of hyperfine coupling (HFC) tensors for radicals consisting of hydrogen, second- and third-row elements. The COS-CCSD(T) method involves a single noniterative step with [Formula: see text] scaling of the computational cost for the calculation of triples corrections to the energy. The key feature of this development is the use of spatial semicanonical orbitals generated from standard restricted open-shell Hartree-Fock (ROHF) orbitals, which allows the unperturbed Hamiltonian operator to be defined in terms of a diagonal spin-free Fock operator. The HFC tensors are computed as a first-order property via implementation of an analytic derivative scheme. The required one-particle spin density matrix is computed by using one- and two-particle spin-free density matrices that are obtained from the analytic derivative implementation, in this way avoiding the use of any spin-dependent operator and maintaining spin adaptation of the CC wavefunction. Benchmark calculations of HFC tensors for a set of 21 radicals indicate reasonably good agreement of the COS-CCSD(T) results with experiment and a consistent improvement over the COS-CCSD method. We demonstrate that the accuracies of the isotropic hyperfine coupling constants obtained in unrestricted HF (UHF) reference based spin-orbital CCSD(T) calculations deteriorate when spin contamination in the UHF wavefunction is large, and the results may even become qualitatively incorrect when spin polarization is the driving mechanism. Within a similar noniterative perturbative treatment of triple excitations, the spin-adapted COS-CCSD(T) approach produces accurate results, thus ensuring cost-effectiveness together with reliability.
Collapse
Affiliation(s)
- Dipayan Datta
- Institut für Physikalische Chemie , Johannes Gutenberg-Universität Mainz , Duesbergweg 10-14 , 55128 Mainz , Germany
| | - Jürgen Gauss
- Institut für Physikalische Chemie , Johannes Gutenberg-Universität Mainz , Duesbergweg 10-14 , 55128 Mainz , Germany
| |
Collapse
|
24
|
Crawford TD, Kumar A, Bazanté AP, Di Remigio R. Reduced‐scaling coupled cluster response theory: Challenges and opportunities. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2019. [DOI: 10.1002/wcms.1406] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- T. Daniel Crawford
- Department of Chemistry Virginia Tech, Blacksburg Virginia
- The Molecular Sciences Software Institute Blacksburg Virginia
| | - Ashutosh Kumar
- Department of Chemistry Virginia Tech, Blacksburg Virginia
| | | | - Roberto Di Remigio
- Department of Chemistry Virginia Tech, Blacksburg Virginia
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry University of Tromsø ‐ The Arctic University of Norway Tromsø Norway
| |
Collapse
|
25
|
Mitxelena I, Piris M, Ugalde JM. Advances in approximate natural orbital functional theory. ADVANCES IN QUANTUM CHEMISTRY 2019. [DOI: 10.1016/bs.aiq.2019.04.001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
26
|
Schlimgen AW, Mazziotti DA. Analytical gradients of variational reduced-density-matrix and wavefunction-based methods from an overlap-reweighted semidefinite program. J Chem Phys 2018; 149:164111. [DOI: 10.1063/1.5043104] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Anthony W. Schlimgen
- Department of Chemistry and The James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| | - David A. Mazziotti
- Department of Chemistry and The James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
27
|
Bozkaya U, Sherrill CD. Analytic energy gradients for the coupled-cluster singles and doubles with perturbative triples method with the density-fitting approximation. J Chem Phys 2018; 147:044104. [PMID: 28764345 DOI: 10.1063/1.4994918] [Citation(s) in RCA: 30] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
An efficient implementation of analytic gradients for the coupled-cluster singles and doubles with perturbative triples [CCSD(T)] method with the density-fitting (DF) approximation, denoted as DF-CCSD(T), is reported. For the molecules considered, the DF approach substantially accelerates conventional CCSD(T) analytic gradients due to the reduced input/output time and the acceleration of the so-called "gradient terms": formation of particle density matrices (PDMs), computation of the generalized Fock-matrix (GFM), solution of the Z-vector equation, formation of the effective PDMs and GFM, back-transformation of the PDMs and GFM, from the molecular orbital to the atomic orbital (AO) basis, and computation of gradients in the AO basis. For the largest member of the molecular test set considered (C6H14), the computational times for analytic gradients (with the correlation-consistent polarized valence triple-ζ basis set in serial) are 106.2 [CCSD(T)] and 49.8 [DF-CCSD(T)] h, a speedup of more than 2-fold. In the evaluation of gradient terms, the DF approach completely avoids the use of four-index two-electron integrals. Similar to our previous studies on DF-second-order Møller-Plesset perturbation theory and DF-CCSD gradients, our formalism employs 2- and 3-index two-particle density matrices (TPDMs) instead of 4-index TPDMs. Errors introduced by the DF approximation are negligible for equilibrium geometries and harmonic vibrational frequencies.
Collapse
Affiliation(s)
- Uğur Bozkaya
- Department of Chemistry, Hacettepe University, Ankara 06800, Turkey
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| |
Collapse
|
28
|
Bozkaya U, Sherrill CD. Analytic energy gradients for the coupled-cluster singles and doubles method with the density-fitting approximation. J Chem Phys 2017; 144:174103. [PMID: 27155621 DOI: 10.1063/1.4948318] [Citation(s) in RCA: 30] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
An efficient implementation is presented for analytic gradients of the coupled-cluster singles and doubles (CCSD) method with the density-fitting approximation, denoted DF-CCSD. Frozen core terms are also included. When applied to a set of alkanes, the DF-CCSD analytic gradients are significantly accelerated compared to conventional CCSD for larger molecules. The efficiency of our DF-CCSD algorithm arises from the acceleration of several different terms, which are designated as the "gradient terms": computation of particle density matrices (PDMs), generalized Fock-matrix (GFM), solution of the Z-vector equation, formation of the relaxed PDMs and GFM, back-transformation of PDMs and GFM to the atomic orbital (AO) basis, and evaluation of gradients in the AO basis. For the largest member of the alkane set (C10H22), the computational times for the gradient terms (with the cc-pVTZ basis set) are 2582.6 (CCSD) and 310.7 (DF-CCSD) min, respectively, a speed up of more than 8-folds. For gradient related terms, the DF approach avoids the usage of four-index electron repulsion integrals. Based on our previous study [U. Bozkaya, J. Chem. Phys. 141, 124108 (2014)], our formalism completely avoids construction or storage of the 4-index two-particle density matrix (TPDM), using instead 2- and 3-index TPDMs. The DF approach introduces negligible errors for equilibrium bond lengths and harmonic vibrational frequencies.
Collapse
Affiliation(s)
- Uğur Bozkaya
- Department of Chemistry, Hacettepe University, Ankara 06800, Turkey
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| |
Collapse
|
29
|
Mitxelena I, Piris M. Analytic gradients for natural orbital functional theory. J Chem Phys 2017; 146:014102. [PMID: 28063441 DOI: 10.1063/1.4973271] [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/14/2022] Open
Abstract
The analytic energy gradients with respect to nuclear motion are derived for the natural orbital functional (NOF) theory. The resulting equations do not require resorting to linear-response theory, so the computation of NOF energy gradients is analogous to gradient calculations at the Hartree-Fock level of theory. The structures of 15 spin-compensated systems, composed of first- and second-row atoms, are optimized employing the conjugate gradient algorithm. As functionals, two orbital-pairing approaches were used, namely, the fifth and sixth Piris NOFs (PNOF5 and PNOF6). For the latter, the obtained equilibrium geometries are compared with coupled cluster singles and doubles calculations and accurate empirical data.
Collapse
Affiliation(s)
- Ion Mitxelena
- Kimika Fakultatea, Euskal Herriko Unibertsitatea (UPV/EHU), P.K. 1072, 20080 Donostia, Spain
| | - Mario Piris
- Kimika Fakultatea, Euskal Herriko Unibertsitatea (UPV/EHU), P.K. 1072, 20080 Donostia, Spain
| |
Collapse
|
30
|
Datta D, Kossmann S, Neese F. Analytic energy derivatives for the calculation of the first-order molecular properties using the domain-based local pair-natural orbital coupled-cluster theory. J Chem Phys 2016. [DOI: 10.1063/1.4962369] [Citation(s) in RCA: 51] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Dipayan Datta
- Max-Planck-Institut für Chemische Energiekonversion, Stiftstrasse 34-36, D-45470 Mülheim an der Ruhr, Germany
| | - Simone Kossmann
- Max-Planck-Institut für Chemische Energiekonversion, Stiftstrasse 34-36, D-45470 Mülheim an der Ruhr, Germany
| | - Frank Neese
- Max-Planck-Institut für Chemische Energiekonversion, Stiftstrasse 34-36, D-45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
31
|
Bozkaya U. A noniterative asymmetric triple excitation correction for the density-fitted coupled-cluster singles and doubles method: Preliminary applications. J Chem Phys 2016; 144:144108. [DOI: 10.1063/1.4945706] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Affiliation(s)
- Uğur Bozkaya
- Department of Chemistry, Hacettepe University, Ankara 06800, Turkey
| |
Collapse
|
32
|
McCarthy MC, Martinez O, McGuire BA, Crabtree KN, Martin-Drumel MA, Stanton JF. Isotopic studies of trans- and cis-HOCO using rotational spectroscopy: Formation, chemical bonding, and molecular structures. J Chem Phys 2016; 144:124304. [DOI: 10.1063/1.4944070] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Affiliation(s)
- Michael C. McCarthy
- Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, Massachusetts 02138, USA and School of Engineering and Applied Sciences, Harvard University, 29 Oxford Street, Cambridge, Massachusetts 02138, USA
| | - Oscar Martinez
- Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, Massachusetts 02138, USA and School of Engineering and Applied Sciences, Harvard University, 29 Oxford Street, Cambridge, Massachusetts 02138, USA
| | - Brett A. McGuire
- Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, Massachusetts 02138, USA and School of Engineering and Applied Sciences, Harvard University, 29 Oxford Street, Cambridge, Massachusetts 02138, USA
- National Radio Astronomy Observatory, Charlottesville, Virginia 22901, USA
| | - Kyle N. Crabtree
- Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, Massachusetts 02138, USA and School of Engineering and Applied Sciences, Harvard University, 29 Oxford Street, Cambridge, Massachusetts 02138, USA
| | - Marie-Aline Martin-Drumel
- Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, Massachusetts 02138, USA and School of Engineering and Applied Sciences, Harvard University, 29 Oxford Street, Cambridge, Massachusetts 02138, USA
| | - John F. Stanton
- Department of Chemistry and Biochemistry, The University of Texas at Austin, 1 University Station A5300, Austin, Texas 78712-0165, USA
| |
Collapse
|
33
|
Eriksen JJ, Jørgensen P, Gauss J. On the convergence of perturbative coupled cluster triples expansions: Error cancellations in the CCSD(T) model and the importance of amplitude relaxation. J Chem Phys 2015; 142:014102. [DOI: 10.1063/1.4904754] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Janus J. Eriksen
- qLEAP Center for Theoretical Chemistry, Department of Chemistry, Aarhus University, Langelandsgade 140, DK–8000 Aarhus C, Denmark
| | - Poul Jørgensen
- qLEAP Center for Theoretical Chemistry, Department of Chemistry, Aarhus University, Langelandsgade 140, DK–8000 Aarhus C, Denmark
| | - Jürgen Gauss
- Institut für Physikalische Chemie, Johannes Gutenberg–Universität Mainz, Duesbergweg 10–14, D–55128 Mainz, Germany
| |
Collapse
|
34
|
Datta D, Gauss J. Analytic first derivatives for a spin-adapted open-shell coupled cluster theory: Evaluation of first-order electrical properties. J Chem Phys 2014; 141:104102. [DOI: 10.1063/1.4894773] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022] Open
Affiliation(s)
- Dipayan Datta
- Institut für Physikalische Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, D-55128 Mainz, Germany
| | - Jürgen Gauss
- Institut für Physikalische Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, D-55128 Mainz, Germany
| |
Collapse
|
35
|
Taylor DE, Rice BM. Quantum-Informed Multiscale M&S for Energetic Materials. ADVANCES IN QUANTUM CHEMISTRY 2014. [DOI: 10.1016/b978-0-12-800345-9.00005-2] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
36
|
Stopkowicz S, Cheng L, Harding ME, Puzzarini C, Gauss J. The bromine nuclear quadrupole moment revisited. Mol Phys 2013. [DOI: 10.1080/00268976.2013.796072] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- Stella Stopkowicz
- a Institut für Physikalische Chemie , Universität Mainz , Mainz , D-55099 , Germany
- b Centre for Theoretical and Computational Chemistry, Department of Chemistry , University of Oslo , Oslo , N-0315 , Norway
| | - Lan Cheng
- a Institut für Physikalische Chemie , Universität Mainz , Mainz , D-55099 , Germany
- c Institute for Theoretical Chemistry, Department of Chemistry and Biochemistry , The University of Texas at Austin , Austin , TX , 78712 , USA
| | - Michael E. Harding
- a Institut für Physikalische Chemie , Universität Mainz , Mainz , D-55099 , Germany
- d Institut für Nanotechnologie, Karlsruher Institut für Technologie (KIT) , Karlsruhe , D-76021 , Germany
| | - Cristina Puzzarini
- e Dipartimento di Chimica ‘Giacomo Ciamician’ , Università di Bologna , Via Selmi 2, Bologna , I-40126 , Italy
| | - Jürgen Gauss
- a Institut für Physikalische Chemie , Universität Mainz , Mainz , D-55099 , Germany
| |
Collapse
|
37
|
Cheng L, Stopkowicz S, Stanton JF, Gauss J. The route to high accuracy inab initiocalculations of Cu quadrupole-coupling constants. J Chem Phys 2012; 137:224302. [DOI: 10.1063/1.4767767] [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
|
38
|
Bozkaya U, Schaefer HF. Symmetric and asymmetric triple excitation corrections for the orbital-optimized coupled-cluster doubles method: improving upon CCSD(T) and CCSD(T)(Λ): preliminary application. J Chem Phys 2012; 136:204114. [PMID: 22667547 DOI: 10.1063/1.4720382] [Citation(s) in RCA: 49] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Symmetric and asymmetric triple excitation corrections for the orbital-optimized coupled-cluster doubles (OO-CCD or simply "OD" for short) method are investigated. The conventional symmetric and asymmetric perturbative triples corrections [(T) and (T)(Λ)] are implemented, the latter one for the first time. Additionally, two new triples corrections, denoted as OD(Λ) and OD(Λ)(T), are introduced. We applied the new methods to potential energy surfaces of the BH, HF, C(2), N(2), and CH(4) molecules, and compare the errors in total energies, with respect to full configuration interaction, with those from the standard coupled-cluster singles and doubles (CCSD), with perturbative triples [CCSD(T)], and asymmetric triples correction (CCSD(T)(Λ)) methods. The CCSD(T) method fails badly at stretched geometries, the corresponding nonparallelity error is 7-281 kcal mol(-1), although it gives reliable results near equilibrium geometries. The new symmetric triples correction, CCSD(Λ), noticeably improves upon CCSD(T) (by 4-14 kcal mol(-1)) for BH, HF, and CH(4); however, its performance is worse than CCSD(T) (by 1.6-4.2 kcal mol(-1)) for C(2) and N(2). The asymmetric triples corrections, CCSD(T)(Λ) and CCSD(Λ)(T), perform remarkably better than CCSD(T) (by 5-18 kcal mol(-1)) for the BH, HF, and CH(4) molecules, while for C(2) and N(2) their results are similar to those of CCSD(T). Although the performance of CCSD and OD is similar, the situation is significantly different in the case of triples corrections, especially at stretched geometries. The OD(T) method improves upon CCSD(T) by 1-279 kcal mol(-1). The new symmetric triples correction, OD(Λ), enhances the OD(T) results (by 0.01-2.0 kcal mol(-1)) for BH, HF, and CH(4); however, its performance is worse than OD(T) (by 1.9-2.3 kcal mol(-1)) for C(2) and N(2). The asymmetric triples corrections, OD(T)(Λ) and OD(Λ)(T), perform better than OD(T) (by 2.0-6.2 kcal mol(-1)). The latter method is slightly better for the BH, HF, and CH(4) molecules. However, for C(2) and N(2) the new results are similar to those of OD(T). For the BH, HF, and CH(4) molecules, OD(Λ)(T) provides the best potential energy curves among the considered methods, while for C(2) and N(2) the OD(T) method prevails. Hence, for single-bond breaking the OD(Λ)(T) method appears to be superior, whereas for multiple-bond breaking the OD(T) method is better.
Collapse
Affiliation(s)
- Uğur Bozkaya
- Department of Chemistry, Atatürk University, Erzurum 25240, Turkey.
| | | |
Collapse
|
39
|
Cazzoli G, Puzzarini C, Gauss J. Rotational spectrum of silyl chloride: hyperfine structure and equilibrium geometry. Mol Phys 2012. [DOI: 10.1080/00268976.2012.680518] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Affiliation(s)
- Gabriele Cazzoli
- a Dipartimento di Chimica ‘G. Ciamician’ , Università degli Studi di Bologna , Via Selmi 2, I-40126 Bologna , Italy
| | - Cristina Puzzarini
- a Dipartimento di Chimica ‘G. Ciamician’ , Università degli Studi di Bologna , Via Selmi 2, I-40126 Bologna , Italy
| | - Jürgen Gauss
- b Institut für Physikalische Chemie , Universität Mainz , D-55099 Mainz , Germany
| |
Collapse
|
40
|
RAVICHANDRAN LALITHA, BHATTACHARYA DEBARATI, VAVAL NAYANA, PAL SOURAV. Fock-space multi-reference coupled-cluster response with the effect of triples on dipole moment of ClO and SF radicals#. J CHEM SCI 2012. [DOI: 10.1007/s12039-012-0220-x] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
|
41
|
Ascik PN, Rugango R, Simmonett AC, Compaan KR, Schaefer HF. The Beryllium Pentamer: Trailing an Uneven Sequence of Dissociation Energies. Chemphyschem 2012; 13:1255-60. [DOI: 10.1002/cphc.201101008] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2011] [Indexed: 11/09/2022]
|
42
|
Helgaker T, Coriani S, Jørgensen P, Kristensen K, Olsen J, Ruud K. Recent Advances in Wave Function-Based Methods of Molecular-Property Calculations. Chem Rev 2012; 112:543-631. [DOI: 10.1021/cr2002239] [Citation(s) in RCA: 463] [Impact Index Per Article: 38.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Trygve Helgaker
- Centre for Theoretical and Computational Chemistry, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway
| | - Sonia Coriani
- Dipartimento di Scienze Chimiche e Farmaceutiche, Università degli Studi di Trieste, Via Giorgieri 1, I-34127 Trieste, Italy
| | - Poul Jørgensen
- Lundbeck Center for Theoretical Chemistry, Department of Chemistry, Aarhus University, 8000 Aarhus C, Denmark
| | - Kasper Kristensen
- Lundbeck Center for Theoretical Chemistry, Department of Chemistry, Aarhus University, 8000 Aarhus C, Denmark
| | - Jeppe Olsen
- Lundbeck Center for Theoretical Chemistry, Department of Chemistry, Aarhus University, 8000 Aarhus C, Denmark
| | - Kenneth Ruud
- Centre for Theoretical and Computational Chemistry, Department of Chemistry, University of Tromsø, N-9037 Tromsø, Norway
| |
Collapse
|
43
|
Kahn K, Kirtman B, Hagen A, Noga J. Communication: Convergence of anharmonic infrared intensities of hydrogen fluoride in traditional and explicitly correlated coupled cluster calculations. J Chem Phys 2011; 135:131103. [DOI: 10.1063/1.3647566] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023] Open
Affiliation(s)
- Kalju Kahn
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, Santa Barbara, California 93106, USA
| | - Bernard Kirtman
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, Santa Barbara, California 93106, USA
| | - Alexander Hagen
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, Santa Barbara, California 93106, USA
| | - Jozef Noga
- Faculty of Natural Sciences, Department of Inorganic Chemistry, Comenius University, Mlynská dolina CH2, SK-84215 Bratislava, Slovakia
- Institute of Inorganic Chemistry, Slovak Academy of Sciences, SK-84536 Bratislava, Slovakia
| |
Collapse
|
44
|
Bozkaya U, Turney JM, Yamaguchi Y, Schaefer HF, Sherrill CD. Quadratically convergent algorithm for orbital optimization in the orbital-optimized coupled-cluster doubles method and in orbital-optimized second-order Møller-Plesset perturbation theory. J Chem Phys 2011; 135:104103. [DOI: 10.1063/1.3631129] [Citation(s) in RCA: 97] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
45
|
Cheng L, Gauss J. Analytical evaluation of first-order electrical properties based on the spin-free Dirac-Coulomb Hamiltonian. J Chem Phys 2011; 134:244112. [DOI: 10.1063/1.3601056] [Citation(s) in RCA: 37] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023] Open
|
46
|
Ravichandran L, Vaval N, Pal S. Effect of Triples to Dipole Moments in Fock-Space Multireference Coupled Cluster Method. J Chem Theory Comput 2011; 7:876-83. [PMID: 26606338 DOI: 10.1021/ct1007263] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
In this paper, we present the new implementation of partial triples for the dipole moment of doublet radicals in Lagrangian formulation of Fock-space multireference coupled cluster (Λ-FSMRCC) response method. We have implemented a specific scheme of noniterative triples, in addition to singles and doubles schemes, which accounts for the effects appearing at least at the third order in dipole moments. The method is applied to the ground states of OH, OOH, HCOO, CN, CH, and PO radicals.
Collapse
Affiliation(s)
| | - Nayana Vaval
- Physical Chemistry Division, National Chemical Laboratory , Pune-411008, India
| | - Sourav Pal
- Physical Chemistry Division, National Chemical Laboratory , Pune-411008, India
| |
Collapse
|
47
|
Ascik PN, Wilke JJ, Simmonett AC, Yamaguchi Y, Schaefer HF. The Beryllium tetramer: Profiling an elusive molecule. J Chem Phys 2011; 134:074110. [DOI: 10.1063/1.3553366] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
|
48
|
Korona T. The effect of local approximations on first-order properties from expectation-value coupled cluster theory. Theor Chem Acc 2010. [DOI: 10.1007/s00214-010-0872-x] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
49
|
Kahn K, Kirtman B, Noga J, Ten-no S. Anharmonic vibrational analysis of water with traditional and explicitly correlated coupled cluster methods. J Chem Phys 2010; 133:074106. [DOI: 10.1063/1.3464837] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
|
50
|
Prochnow E, Harding ME, Gauss J. Parallel Calculation of CCSDT and Mk-MRCCSDT Energies. J Chem Theory Comput 2010; 6:2339-47. [DOI: 10.1021/ct1002016] [Citation(s) in RCA: 37] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Eric Prochnow
- Institut für Physikalische Chemie, Universität Mainz, D-55099 Mainz, Germany
| | - Michael E. Harding
- Institut für Physikalische Chemie, Universität Mainz, D-55099 Mainz, Germany
| | - Jürgen Gauss
- Institut für Physikalische Chemie, Universität Mainz, D-55099 Mainz, Germany
| |
Collapse
|