1
|
Majee K, Chakraborty S, Mukhopadhyay T, Nayak MK, Dutta AK. A reduced cost four-component relativistic unitary coupled cluster method for atoms and molecules. J Chem Phys 2024; 161:034101. [PMID: 39007370 DOI: 10.1063/5.0207091] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2024] [Accepted: 06/12/2024] [Indexed: 07/16/2024] Open
Abstract
We present a four-component relativistic unitary coupled cluster method for atoms and molecules. We have used commutator-based non-perturbative approximation using the "Bernoulli expansion" to derive an approximation to the relativistic unitary coupled cluster method. The performance of the full quadratic unitary coupled-cluster singles and doubles method (qUCCSD), as well as a perturbative approximation variant (UCC3), has been reported for both energies and properties. It can be seen that both methods give results comparable to those of the standard relativistic coupled cluster method. The qUCCSD method shows better agreement with experimental results due to the better inclusion of relaxation effects. The relativistic UCC3 and qUCCSD methods can simulate the spin-forbidden transition with easy access to transition properties. A natural spinor-based scheme to reduce the computational cost of relativistic UCC3 and qUCCSD methods has been discussed.
Collapse
Affiliation(s)
- Kamal Majee
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Sudipta Chakraborty
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Tamoghna Mukhopadhyay
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Malaya K Nayak
- Theoretical Chemistry Section, Bhabha Atomic Research Centre, Trombay, Mumbai 400085, India
- Homi Bhabha National Institute, BARC Training School Complex, Anushakti Nagar, Mumbai 400094, India
| | - Achintya Kumar Dutta
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| |
Collapse
|
2
|
Yuan X, Halbert L, Pototschnig JV, Papadopoulos A, Coriani S, Visscher L, Pereira Gomes AS. Formulation and Implementation of Frequency-Dependent Linear Response Properties with Relativistic Coupled Cluster Theory for GPU-Accelerated Computer Architectures. J Chem Theory Comput 2024; 20:677-694. [PMID: 38193434 DOI: 10.1021/acs.jctc.3c00812] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2024]
Abstract
We present the development and implementation of relativistic coupled cluster linear response theory (CC-LR), which allows the determination of molecular properties arising from time-dependent or time-independent electric, magnetic, or mixed electric-magnetic perturbations (within a common gauge origin for the magnetic properties) as well as taking into account the finite lifetime of excited states in the framework of damped response theory. We showcase our implementation, which is capable to offload the computationally intensive tensor contractions characteristic of coupled cluster theory onto graphical processing units, in the calculation of (a) frequency-(in)dependent dipole-dipole polarizabilities of IIB atoms and selected diatomic molecules, with a particular emphasis on the calculation of valence absorption cross sections for the I2 molecule; (b) indirect spin-spin coupling constants for benchmark systems such as the hydrogen halides (HX, X = F-I) as well the H2Se-H2O dimer as a prototypical system containing hydrogen bonds; and (c) optical rotations at the sodium D line for hydrogen peroxide analogues (H2Y2, Y = O, S, Se, Te). Thanks to this implementation, we are able to show the similarities in performance, but often the significant discrepancies, between CC-LR and approximate methods such as density functional theory. Comparing standard CC response theory with the flavor based upon the equation of motion formalism, we find that for valence properties such as polarizabilities, the two frameworks yield very similar results across the periodic table as found elsewhere in the literature; for properties that probe the core region, such as spin-spin couplings, on the other hand, we show a progressive differentiation between the two as relativistic effects become more important. Our results also suggest that as one goes down the periodic table, it may become increasingly difficult to measure pure optical rotation at the sodium D line due to the appearance of absorbing states.
Collapse
Affiliation(s)
- Xiang Yuan
- Univ. Lille, CNRS, UMR 8523─PhLAM─Physique des Lasers Atomes et Molécules, F-59000 Lille, France
- Department of Chemistry and Pharmaceutical Sciences, Faculty of Science, Vrije Universiteit Amsterdam, 1081 HV Amsterdam, The Netherlands
| | - Loïc Halbert
- Univ. Lille, CNRS, UMR 8523─PhLAM─Physique des Lasers Atomes et Molécules, F-59000 Lille, France
| | - Johann Valentin Pototschnig
- Department of Chemistry and Pharmaceutical Sciences, Faculty of Science, Vrije Universiteit Amsterdam, 1081 HV Amsterdam, The Netherlands
| | - Anastasios Papadopoulos
- Department of Chemistry and Pharmaceutical Sciences, Faculty of Science, Vrije Universiteit Amsterdam, 1081 HV Amsterdam, The Netherlands
| | - Sonia Coriani
- DTU Chemistry─Department of Chemistry, Technical University of Denmark, DK-2800 Kongens Lyngby, Denmark
| | - Lucas Visscher
- Department of Chemistry and Pharmaceutical Sciences, Faculty of Science, Vrije Universiteit Amsterdam, 1081 HV Amsterdam, The Netherlands
| | | |
Collapse
|
3
|
Yuan X, Halbert L, Visscher L, Pereira Gomes AS. Frequency-Dependent Quadratic Response Properties and Two-Photon Absorption from Relativistic Equation-of-Motion Coupled Cluster Theory. J Chem Theory Comput 2023; 19:9248-9259. [PMID: 38079602 DOI: 10.1021/acs.jctc.3c01011] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2023]
Abstract
We present the implementation of quadratic response theory based upon the relativistic equation-of-motion coupled cluster method. We showcase our implementation, whose generality allows us to consider both time-dependent and time-independent electric and magnetic perturbations, by considering the static and frequency-dependent hyperpolarizability of hydrogen halides (HX, X = F-At), providing comprehensive insights into their electronic response characteristics. Additionally, we evaluated the Verdet constant for noble gases Xe and Rn and discussed the relative importance of relativistic and electron correlation effects for these magneto-optical properties. Finally, we calculate the two-photon absorption cross sections of transition [ns1S0 → (n + 1)s1S0] of Ga+ and In+, which are suggested as candidates for new ion clocks. As our implementation allows for the use of nonrelativistic Hamiltonians as well, we have compared our EOM-QRCC results to the QR-CC implementation in the DALTON code and show that the differences between CC and EOMCC response are in general smaller than 5% for the properties considered. Collectively, the results underscore the versatility of our implementation and its potential as a benchmark tool for other approximated models, such as density functional theory for higher-order properties.
Collapse
Affiliation(s)
- Xiang Yuan
- Physique des Lasers Atomes et Molecules, Universite de Lille, F-59000 Lille, France
- Department of Chemistry and Pharmaceutical Sciences, Faculty of Science, Vrije Universiteit Amsterdam, 1081 HV Amsterdam, The Netherlands
| | - Loïc Halbert
- Physique des Lasers Atomes et Molecules, Universite de Lille, F-59000 Lille, France
| | - Lucas Visscher
- Department of Chemistry and Pharmaceutical Sciences, Faculty of Science, Vrije Universiteit Amsterdam, 1081 HV Amsterdam, The Netherlands
| | | |
Collapse
|
4
|
Peyton BG, Stewart ZJ, Weidman JD, Wilson AK. Tailoring light-induced charge transfer and intersystem crossing in FeCO using time-dependent spin-orbit configuration interaction. J Chem Phys 2023; 159:204108. [PMID: 38014783 DOI: 10.1063/5.0173529] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2023] [Accepted: 11/02/2023] [Indexed: 11/29/2023] Open
Abstract
Real-time (RT) electronic structure methods provide a natural framework for describing light-matter interactions in arbitrary time-dependent electromagnetic fields (EMF). Optically induced excited state transitions are of particular interest, which require tuned EMF to drive population transfer to and from the specific state(s) of interest. Intersystem crossing, or spin-flip, may be driven through shaped EMF or laser pulses. These transitions can result in long-lived "spin-trapped" excited states, which are especially useful for materials requiring charge separation or protracted excited state lifetimes. Time-dependent configuration interaction (TDCI) is unique among RT methods in that it may be implemented in a basis of eigenstates, allowing for rapid propagation of the time-dependent Schrödinger equation. The recent spin-orbit TDCI (TD-SOCI) enables a real-time description of spin-flip dynamics in an arbitrary EMF and, therefore, provides an ideal framework for rational pulse design. The present study explores the mechanism of multiple spin-flip pathways for a model transition metal complex, FeCO, using shaped pulses designed to drive controlled intersystem crossing and charge transfer. These results show that extremely tunable excited state dynamics can be achieved by considering the dipole transition matrix elements between the states of interest.
Collapse
Affiliation(s)
- Benjamin G Peyton
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Zachary J Stewart
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Jared D Weidman
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Angela K Wilson
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| |
Collapse
|
5
|
Di Felice R, Mayes ML, Richard RM, Williams-Young DB, Chan GKL, de Jong WA, Govind N, Head-Gordon M, Hermes MR, Kowalski K, Li X, Lischka H, Mueller KT, Mutlu E, Niklasson AMN, Pederson MR, Peng B, Shepard R, Valeev EF, van Schilfgaarde M, Vlaisavljevich B, Windus TL, Xantheas SS, Zhang X, Zimmerman PM. A Perspective on Sustainable Computational Chemistry Software Development and Integration. J Chem Theory Comput 2023; 19:7056-7076. [PMID: 37769271 PMCID: PMC10601486 DOI: 10.1021/acs.jctc.3c00419] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/14/2023] [Indexed: 09/30/2023]
Abstract
The power of quantum chemistry to predict the ground and excited state properties of complex chemical systems has driven the development of computational quantum chemistry software, integrating advances in theory, applied mathematics, and computer science. The emergence of new computational paradigms associated with exascale technologies also poses significant challenges that require a flexible forward strategy to take full advantage of existing and forthcoming computational resources. In this context, the sustainability and interoperability of computational chemistry software development are among the most pressing issues. In this perspective, we discuss software infrastructure needs and investments with an eye to fully utilize exascale resources and provide unique computational tools for next-generation science problems and scientific discoveries.
Collapse
Affiliation(s)
- Rosa Di Felice
- Departments
of Physics and Astronomy and Quantitative and Computational Biology, University of Southern California, Los Angeles, California 90089, United States
- CNR-NANO
Modena, Modena 41125, Italy
| | - Maricris L. Mayes
- Department
of Chemistry and Biochemistry, University
of Massachusetts Dartmouth, North Dartmouth, Massachusetts 02747, United States
| | | | | | - Garnet Kin-Lic Chan
- Division
of Chemistry and Chemical Engineering, California
Institute of Technology, Pasadena, California 91125, United States
| | - Wibe A. de Jong
- Lawrence
Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Niranjan Govind
- Physical
Sciences Division, Pacific Northwest National
Laboratory, Richland, Washington 99354, United States
| | - Martin Head-Gordon
- Pitzer Center
for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Matthew R. Hermes
- Department
of Chemistry, Chicago Center for Theoretical Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| | - Karol Kowalski
- Physical
Sciences Division, Pacific Northwest National
Laboratory, Richland, Washington 99354, United States
| | - Xiaosong Li
- Department
of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Hans Lischka
- Department
of Chemistry and Biochemistry, Texas Tech
University, Lubbock, Texas 79409, United States
| | - Karl T. Mueller
- Physical
and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Erdal Mutlu
- Advanced
Computing, Mathematics, and Data Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Anders M. N. Niklasson
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Mark R. Pederson
- Department
of Physics, The University of Texas at El
Paso, El Paso, Texas 79968, United States
| | - Bo Peng
- Physical
Sciences Division, Pacific Northwest National
Laboratory, Richland, Washington 99354, United States
| | - Ron Shepard
- Chemical
Sciences and Engineering Division, Argonne
National Laboratory, Lemont, Illinois 60439, United States
| | - Edward F. Valeev
- Department
of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | | | - Bess Vlaisavljevich
- Department
of Chemistry, University of South Dakota, Vermillion, South Dakota 57069, United States
| | - Theresa L. Windus
- Department
of Chemistry, Iowa State University and
Ames Laboratory, Ames, Iowa 50011, United States
| | - Sotiris S. Xantheas
- Department
of Chemistry, University of Washington, Seattle, Washington 98195, United States
- Advanced
Computing, Mathematics and Data Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Xing Zhang
- Division
of Chemistry and Chemical Engineering, California
Institute of Technology, Pasadena, California 91125, United States
| | - Paul M. Zimmerman
- Department
of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
6
|
Peyton BG, Wang Z, Crawford TD. Reduced Scaling Real-Time Coupled Cluster Theory. J Phys Chem A 2023; 127:8486-8499. [PMID: 37782945 DOI: 10.1021/acs.jpca.3c05151] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/04/2023]
Abstract
Real-time coupled cluster (CC) methods have several advantages over their frequency-domain counterparts, namely, response and equation of motion CC theories. Broadband spectra, strong fields, and pulse manipulation allow for the simulation of complex spectroscopies that are unreachable using frequency-domain approaches. Due to the high-order polynomial scaling, the required numerical time propagation of the CC residual expressions is a computationally demanding process. This scaling may be reduced by local correlation schemes, which aim to reduce the size of the (virtual) orbital space by truncation according to user-defined parameters. We present the first application of local correlation to real-time CC. As in previous studies of locally correlated frequency-domain CC, traditional local correlation schemes are of limited utility for field-dependent properties; however, a perturbation-aware scheme proves promising. A detailed analysis of the amplitude dynamics suggests that the main challenge is a strong time dependence of the wave function sparsity.
Collapse
Affiliation(s)
- Benjamin G Peyton
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Zhe Wang
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - T Daniel Crawford
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| |
Collapse
|
7
|
Barcza B, Szirmai Á, Tajti A, Stanton JF, Szalay PG. Benchmarking Aspects of Ab Initio Fragment Models for Accurate Excimer Potential Energy Surfaces. J Chem Theory Comput 2023; 19:3580-3600. [PMID: 37236166 PMCID: PMC10694823 DOI: 10.1021/acs.jctc.3c00104] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2023] [Indexed: 05/28/2023]
Abstract
While Coupled-Cluster methods have been proven to provide an accurate description of excited electronic states, the scaling of the computational costs with the system size limits the degree for which these methods can be applied. In this work different aspects of fragment-based approaches are studied on noncovalently bound molecular complexes with interacting chromophores of the fragments, such as π-stacked nucleobases. The interaction of the fragments is considered at two distinct steps. First, the states localized on the fragments are described in the presence of the other fragment(s); for this we test two approaches. One method is founded on QM/MM principles, only including the electrostatic interaction between the fragments in the electronic structure calculation with Pauli repulsion and dispersion effects added separately. The other model, a Projection-based Embedding (PbE) using the Huzinaga equation, includes both electrostatic and Pauli repulsion and only needs to be augmented by dispersion interactions. In both schemes the extended Effective Fragment Potential (EFP2) method of Gordon et al. was found to provide an adequate correction for the missing terms. In the second step, the interaction of the localized chromophores is modeled for a proper description of the excitonic coupling. Here the inclusion of purely electrostatic contributions appears to be sufficient: it is found that the Coulomb part of the coupling provides accurate splitting of the energies of interacting chromophores that are separated by more than 4 Å.
Collapse
Affiliation(s)
- Bónis Barcza
- Laboratory
of Theoretical Chemistry, Institute of Chemistry, ELTE Eötvös Loránd University, P.O. Box 32, H-1117 Budapest, Hungary
- György
Hevesy Doctoral School, Institute of Chemistry, ELTE Eötvös Loránd University, H-1117 Budapest, Hungary
| | - Ádám
B. Szirmai
- Laboratory
of Theoretical Chemistry, Institute of Chemistry, ELTE Eötvös Loránd University, P.O. Box 32, H-1117 Budapest, Hungary
- György
Hevesy Doctoral School, Institute of Chemistry, ELTE Eötvös Loránd University, H-1117 Budapest, Hungary
| | - Attila Tajti
- Laboratory
of Theoretical Chemistry, Institute of Chemistry, ELTE Eötvös Loránd University, P.O. Box 32, H-1117 Budapest, Hungary
| | - John F. Stanton
- Quantum
Theory Project, Department of Chemistry, University of Florida, Gainesville, Florida 32611, United States
| | - Péter G. Szalay
- Laboratory
of Theoretical Chemistry, Institute of Chemistry, ELTE Eötvös Loránd University, P.O. Box 32, H-1117 Budapest, Hungary
| |
Collapse
|
8
|
Mehendale SG, Peng B, Govind N, Alexeev Y. Exploring Parameter Redundancy in the Unitary Coupled-Cluster Ansätze for Hybrid Variational Quantum Computing. J Phys Chem A 2023; 127:4526-4537. [PMID: 37193645 DOI: 10.1021/acs.jpca.3c00550] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/18/2023]
Abstract
One of the commonly used chemically inspired approaches in variational quantum computing is the unitary coupled-cluster (UCC) ansätze. Despite being a systematic way of approaching the exact limit, the number of parameters in the standard UCC ansätze exhibits unfavorable scaling with respect to the system size, hindering its practical use on near-term quantum devices. Efforts have been taken to propose some variants of the UCC ansätze with better scaling. In this paper, we explore the parameter redundancy in the preparation of unitary coupled-cluster singles and doubles (UCCSD) ansätze employing spin-adapted formulation, small amplitude filtration, and entropy-based orbital selection approaches. Numerical results of using our approach on some small molecules have exhibited a significant cost reduction in the number of parameters to be optimized and in the time to convergence compared with conventional UCCSD-VQE simulations. We also discuss the potential application of some machine learning techniques in further exploring the parameter redundancy, providing a possible direction for future studies.
Collapse
Affiliation(s)
- Shashank G Mehendale
- Indian Institute of Science Education and Research (IISER), Kolkata, West Bengal 741246, India
| | - Bo Peng
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Niranjan Govind
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Yuri Alexeev
- Computational Science Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| |
Collapse
|
9
|
Sacchetta F, Graf D, Laqua H, Ambroise MA, Kussmann J, Dreuw A, Ochsenfeld C. An effective sub-quadratic scaling atomic-orbital reformulation of the scaled opposite-spin RI-CC2 ground-state model using Cholesky-decomposed densities and an attenuated Coulomb metric. J Chem Phys 2022; 157:104104. [DOI: 10.1063/5.0098719] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
An atomic-orbital reformulation of the Laplace-transformed scaled opposite-spin (SOS) coupled cluster singles and doubles (CC2) model within the resolution of the identity (RI) approximation (SOS-RI-CC2) is presented that extends its applicability to molecules with several hundreds of atoms and triple-zeta basis sets. We exploit sparse linear algebra and an attenuated Coulomb metric to decrease the disk space demands and the computational efforts. In this way, an effective sub-quadratic computational scaling is achieved with our ω-SOS-CDD-RI-CC2 model. Moreover, Cholesky decomposition of the ground-state one-electron density matrix reduces the prefactor, allowing for an early crossover with the molecular orbital formulation. The accuracy and performance of the presented method are investigated for various molecular systems.
Collapse
Affiliation(s)
- F. Sacchetta
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Munich, Germany
| | - D. Graf
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Munich, Germany
| | - H. Laqua
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Munich, Germany
| | - M. A. Ambroise
- Chair of Theoretical and Computational Chemistry, Interdisciplinary Center for Scientific Computing, Heidelberg University, Heidelberg, Germany
| | - J. Kussmann
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Munich, Germany
| | - A. Dreuw
- Chair of Theoretical and Computational Chemistry, Interdisciplinary Center for Scientific Computing, Heidelberg University, Heidelberg, Germany
| | - C. Ochsenfeld
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Munich, Germany
| |
Collapse
|
10
|
D'Cunha R, Crawford TD. Applications of a perturbation-aware local correlation method to coupled cluster linear response properties. Mol Phys 2022. [DOI: 10.1080/00268976.2022.2112627] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/15/2022]
Affiliation(s)
- Ruhee D'Cunha
- Department of Chemistry, Virginia Tech, Blacksburg, VA, USA
| | - T. Daniel Crawford
- Department of Chemistry, Virginia Tech, Blacksburg, VA, USA
- Molecular Sciences Software Institute, Blacksburg, VA, USA
| |
Collapse
|
11
|
Wang Z, Peyton BG, Crawford TD. Accelerating Real-Time Coupled Cluster Methods with Single-Precision Arithmetic and Adaptive Numerical Integration. J Chem Theory Comput 2022; 18:5479-5491. [PMID: 35939815 DOI: 10.1021/acs.jctc.2c00490] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We explore the framework of a real-time coupled cluster method with a focus on improving its computational efficiency. Propagation of the wave function via the time-dependent Schrödinger equation places high demands on computing resources, particularly for high level theories such as coupled cluster with polynomial scaling. Similar to earlier investigations of coupled cluster properties, we demonstrate that the use of single-precision arithmetic reduces both the storage and multiplicative costs of the real-time simulation by approximately a factor of 2 with no significant impact on the resulting UV/vis absorption spectrum computed via the Fourier transform of the time-dependent dipole moment. Additional speedups─of up to a factor of 14 in test simulations of water clusters─are obtained via a straightforward GPU-based implementation as compared to conventional CPU calculations. We also find that further performance optimization is accessible through sagacious selection of numerical integration algorithms, and the adaptive methods, such as the Cash-Karp integrator, provide an effective balance between computing costs and numerical stability. Finally, we demonstrate that a simple mixed-step integrator based on the conventional fourth-order Runge-Kutta approach is capable of stable propagations even for strong external fields, provided the time step is appropriately adapted to the duration of the laser pulse with only minimal computational overhead.
Collapse
Affiliation(s)
- Zhe Wang
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Benjamin G Peyton
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - T Daniel Crawford
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| |
Collapse
|
12
|
Paul AC, Folkestad SD, Myhre RH, Koch H. Oscillator Strengths in the Framework of Equation of Motion Multilevel CC3. J Chem Theory Comput 2022; 18:5246-5258. [PMID: 35921447 PMCID: PMC9476665 DOI: 10.1021/acs.jctc.2c00164] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
We present an efficient implementation of the equation
of motion
oscillator strengths for the closed-shell multilevel coupled cluster
singles and doubles with perturbative triples method (MLCC3) in the
electronic structure program eT. The orbital space is split into an active part treated with
CC3 and an inactive part computed at the coupled cluster singles and
doubles (CCSD) level of theory. Asymptotically, the CC3 contribution
scales as floating-point operations, where nV is the total number of virtual orbitals while nv and no are the
number of active virtual and occupied orbitals, respectively. The
CC3 contribution, thus, only scales linearly with the full system
size and can become negligible compared to the cost of CCSD. We demonstrate
the capabilities of our implementation by calculating the ultraviolet–visible
spectrum of azobenzene and a core excited state of betaine 30 with
more than 1000 molecular orbitals.
Collapse
Affiliation(s)
- Alexander C Paul
- Department of Chemistry, Norwegian University of Science and Technology, NTNU, 7491 Trondheim, Norway
| | - Sarai Dery Folkestad
- Department of Chemistry, Norwegian University of Science and Technology, NTNU, 7491 Trondheim, Norway
| | - Rolf H Myhre
- Department of Chemistry, Norwegian University of Science and Technology, NTNU, 7491 Trondheim, Norway
| | - Henrik Koch
- Department of Chemistry, Norwegian University of Science and Technology, NTNU, 7491 Trondheim, Norway.,Scuola Normale Superiore, Piazza dei Cavaleri 7, 56126 Pisa, Italy
| |
Collapse
|
13
|
Kristiansen HE, Ofstad BS, Hauge E, Aurbakken E, Schøyen ØS, Kvaal S, Pedersen TB. Linear and Nonlinear Optical Properties from TDOMP2 Theory. J Chem Theory Comput 2022; 18:3687-3702. [PMID: 35436120 PMCID: PMC9202312 DOI: 10.1021/acs.jctc.1c01309] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
We present a derivation
of real-time (RT) time-dependent orbital-optimized
Møller–Plesset (TDOMP2) theory and its biorthogonal companion,
time-dependent non-orthogonal OMP2 theory, starting from the time-dependent
bivariational principle and a parametrization based on the exponential
orbital-rotation operator formulation commonly used in the time-independent
molecular electronic structure theory. We apply the TDOMP2 method
to extract absorption spectra and frequency-dependent polarizabilities
and first hyperpolarizabilities from RT simulations, comparing the
results with those obtained from conventional time-dependent coupled-cluster
singles and doubles (TDCCSD) simulations and from its second-order
approximation, TDCC2. We also compare our results with those from
CCSD and CC2 linear and quadratic response theories. Our results indicate
that while TDOMP2 absorption spectra are of the same quality as TDCC2
spectra, including core excitations where optimized orbitals might
be particularly important, frequency-dependent polarizabilities and
hyperpolarizabilities from TDOMP2 simulations are significantly closer
to TDCCSD results than those from TDCC2 simulations.
Collapse
Affiliation(s)
- Håkon Emil Kristiansen
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, Oslo N-0315, Norway
| | - Benedicte Sverdrup Ofstad
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, Oslo N-0315, Norway
| | - Eirill Hauge
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, Oslo N-0315, Norway.,Simula Research Laboratory, Kristian Augusts Gate 23, Oslo 0164, Norway
| | - Einar Aurbakken
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, Oslo N-0315, Norway
| | | | - Simen Kvaal
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, Oslo N-0315, Norway.,Centre for Advanced Study at the Norwegian Academy of Science and Letters, Drammensveien 78, Oslo N-0271, Norway
| | - Thomas Bondo Pedersen
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, Oslo N-0315, Norway.,Centre for Advanced Study at the Norwegian Academy of Science and Letters, Drammensveien 78, Oslo N-0271, Norway
| |
Collapse
|
14
|
Chamoli S, Surjuse K, Jangid B, Nayak MK, Dutta AK. A reduced cost four-component relativistic coupled cluster method based on natural spinors. J Chem Phys 2022; 156:204120. [PMID: 35649878 DOI: 10.1063/5.0085932] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023] Open
Abstract
We present the theory, implementation, and benchmark results for a frozen natural spinors based reduced cost four-component relativistic coupled cluster method. The natural spinors are obtained by diagonalizing the one-body reduced density matrix from a relativistic second-order Møller-Plesset calculation based on a four-component Dirac-Coulomb Hamiltonian. The correlation energy in the coupled cluster method converges more rapidly with respect to the size of the virtual space in the frozen natural spinor basis than that observed in the standard canonical spinors obtained from the Dirac-Hartree-Fock calculation. The convergence of properties is not smooth in the frozen natural spinor basis. However, the inclusion of the perturbative correction smoothens the convergence of the properties with respect to the size of the virtual space in the frozen natural spinor basis and greatly reduces the truncation errors in both energy and property calculations. The accuracy of the frozen natural spinor based coupled cluster methods can be controlled by a single threshold and is a black box to use.
Collapse
Affiliation(s)
- Somesh Chamoli
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Kshitijkumar Surjuse
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Bhavnesh Jangid
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Malaya K Nayak
- Theoretical Chemistry Section, Bhabha Atomic Research Centre, Trombay, Mumbai 400085, India
| | - Achintya Kumar Dutta
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| |
Collapse
|
15
|
Hohenstein EG, Fales BS, Parrish RM, Martínez TJ. Rank-reduced coupled-cluster. III. Tensor hypercontraction of the doubles amplitudes. J Chem Phys 2022; 156:054102. [DOI: 10.1063/5.0077770] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Edward G. Hohenstein
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, USA
| | - B. Scott Fales
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, USA
| | | | - Todd J. Martínez
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, USA
| |
Collapse
|
16
|
Hohenstein EG, Martínez TJ. GPU acceleration of rank-reduced coupled-cluster singles and doubles. J Chem Phys 2021; 155:184110. [PMID: 34773962 DOI: 10.1063/5.0063467] [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
We have developed a graphical processing unit (GPU) accelerated implementation of our recently introduced rank-reduced coupled-cluster singles and doubles (RR-CCSD) method. RR-CCSD introduces a low-rank approximation of the doubles amplitudes. This is combined with a low-rank approximation of the electron repulsion integrals via Cholesky decomposition. The result of these two low-rank approximations is the replacement of the usual fourth-order CCSD tensors with products of second- and third-order tensors. In our implementation, only a single fourth-order tensor must be constructed as an intermediate during the solution of the amplitude equations. Owing in large part to the compression of the doubles amplitudes, the GPU-accelerated implementation shows excellent parallel efficiency (95% on eight GPUs). Our implementation can solve the RR-CCSD equations for up to 400 electrons and 1550 basis functions-roughly 50% larger than the largest canonical CCSD computations that have been performed on any hardware. In addition to increased scalability, the RR-CCSD computations are faster than the corresponding CCSD computations for all but the smallest molecules. We test the accuracy of RR-CCSD for a variety of chemical systems including up to 1000 basis functions and determine that accuracy to better than 0.1% error in the correlation energy can be achieved with roughly 95% compression of the ov space for the largest systems considered. We also demonstrate that conformational energies can be predicted to be within 0.1 kcal mol-1 with efficient compression applied to the wavefunction. Finally, we find that low-rank approximations of the CCSD doubles amplitudes used in the similarity transformation of the Hamiltonian prior to a conventional equation-of-motion CCSD computation will not introduce significant errors (on the order of a few hundredths of an electronvolt) into the resulting excitation energies.
Collapse
Affiliation(s)
- Edward G Hohenstein
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
| | - Todd J Martínez
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
| |
Collapse
|
17
|
D'Cunha R, Crawford TD. PNO++: Perturbed Pair Natural Orbitals for Coupled Cluster Linear Response Theory. J Chem Theory Comput 2021; 17:290-301. [PMID: 33351627 DOI: 10.1021/acs.jctc.0c01086] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Reduced-scaling methods are needed to make accurate and systematically improvable coupled cluster linear response methods for the calculation of molecular properties tractable for large molecules. In this paper, we examine the perturbed pair natural orbital-based PNO++ approach that creates an orbital space optimized for response properties derived from a lower-cost field-perturbed density matrix. We analyze truncation errors in correlation energies, dynamic polarizabilities, and specific rotations from a coupled cluster singles and doubles (CCSD) reference. We find that incorporating a fixed number of orbitals from the pair natural orbital (PNO) space into the PNO++ method-a new method presented here, the "combined PNO++" approach-recovers accuracy in the CCSD correlation energy while preserving the well-behaved convergence behavior of the PNO++ method for linear response properties.
Collapse
Affiliation(s)
- Ruhee D'Cunha
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - T Daniel Crawford
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States.,Molecular Sciences Software Institute, 1880 Pratt Drive, Suite 1100, Blacksburg, Virginia 24060, United States
| |
Collapse
|
18
|
Kumar A, Neese F, Valeev EF. Explicitly correlated coupled cluster method for accurate treatment of open-shell molecules with hundreds of atoms. J Chem Phys 2020; 153:094105. [PMID: 32891102 DOI: 10.1063/5.0012753] [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
We present a near-linear scaling formulation of the explicitly correlated coupled-cluster singles and doubles with the perturbative triples method [CCSD(T)F12¯] for high-spin states of open-shell species. The approach is based on the conventional open-shell CCSD formalism [M. Saitow et al., J. Chem. Phys. 146, 164105 (2017)] utilizing the domain local pair-natural orbitals (DLPNO) framework. The use of spin-independent set of pair-natural orbitals ensures exact agreement with the closed-shell formalism reported previously, with only marginally impact on the cost (e.g., the open-shell formalism is only 1.5 times slower than the closed-shell counterpart for the C160H322 n-alkane, with the measured size complexity of ≈1.2). Evaluation of coupled-cluster energies near the complete-basis-set (CBS) limit for open-shell systems with more than 550 atoms and 5000 basis functions is feasible on a single multi-core computer in less than 3 days. The aug-cc-pVTZ DLPNO-CCSD(T)F12¯ contribution to the heat of formation for the 50 largest molecules among the 348 core combustion species benchmark set [J. Klippenstein et al., J. Phys. Chem. A 121, 6580-6602 (2017)] had root-mean-square deviation (RMSD) from the extrapolated CBS CCSD(T) reference values of 0.3 kcal/mol. For a more challenging set of 50 reactions involving small closed- and open-shell molecules [G. Knizia et al., J. Chem. Phys. 130, 054104 (2009)], the aug-cc-pVQ(+d)Z DLPNO-CCSD(T)F12¯ yielded a RMSD of ∼0.4 kcal/mol with respect to the CBS CCSD(T) estimate.
Collapse
Affiliation(s)
- Ashutosh Kumar
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Frank Neese
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Edward F Valeev
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| |
Collapse
|
19
|
Pavošević F, Tao Z, Culpitt T, Zhao L, Li X, Hammes-Schiffer S. Frequency and Time Domain Nuclear-Electronic Orbital Equation-of-Motion Coupled Cluster Methods: Combination Bands and Electronic-Protonic Double Excitations. J Phys Chem Lett 2020; 11:6435-6442. [PMID: 32658486 DOI: 10.1021/acs.jpclett.0c01891] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
The accurate description of excited vibronic states is important for modeling a wide range of photoinduced processes. The nuclear-electronic orbital (NEO) approach, which treats specified protons on the same level as the electrons, can describe excited electronic-protonic states. Herein the multicomponent equation-of-motion coupled cluster with singles and doubles (NEO-EOM-CCSD) method and its time-domain counterpart, TD-NEO-EOM-CCSD, are developed and implemented. The application of these methods to the HCN molecule highlights their capabilities. These methods predict qualitatively reasonable energies and intensities for a combination band corresponding to simultaneous excitation of two vibrational modes, as well as an overtone. These methods also describe states with double excitation character, such as excited electronic-protonic states corresponding to the simultaneous excitation of an electron and a proton. The ability of the NEO-EOM-CCSD method and its time-dependent counterpart to describe combination bands, overtones, and double excitations will enable a wide range of photochemical applications.
Collapse
Affiliation(s)
- Fabijan Pavošević
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Zhen Tao
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Tanner Culpitt
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Luning Zhao
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Sharon Hammes-Schiffer
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| |
Collapse
|
20
|
Peng C, Lewis CA, Wang X, Clement MC, Pierce K, Rishi V, Pavošević F, Slattery S, Zhang J, Teke N, Kumar A, Masteran C, Asadchev A, Calvin JA, Valeev EF. Massively Parallel Quantum Chemistry: A high-performance research platform for electronic structure. J Chem Phys 2020; 153:044120. [DOI: 10.1063/5.0005889] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Chong Peng
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Cannada A. Lewis
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Xiao Wang
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | | | - Karl Pierce
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Varun Rishi
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Fabijan Pavošević
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Samuel Slattery
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Jinmei Zhang
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Nakul Teke
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Ashutosh Kumar
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Conner Masteran
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Andrey Asadchev
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Justus A. Calvin
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Edward F. Valeev
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| |
Collapse
|
21
|
Peyton BG, Briggs C, D’Cunha R, Margraf JT, Crawford TD. Machine-Learning Coupled Cluster Properties through a Density Tensor Representation. J Phys Chem A 2020; 124:4861-4871. [DOI: 10.1021/acs.jpca.0c02804] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Affiliation(s)
- Benjamin G. Peyton
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Connor Briggs
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Ruhee D’Cunha
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Johannes T. Margraf
- Chair for Theoretical Chemistry, Technische Universität München, Lichtenbergstrasse 4, D-85747 Garching, Germany
| | - T. Daniel Crawford
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| |
Collapse
|
22
|
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: 337] [Impact Index Per Article: 84.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
|
23
|
Pokhilko P, Izmodenov D, Krylov AI. Extension of frozen natural orbital approximation to open-shell references: Theory, implementation, and application to single-molecule magnets. J Chem Phys 2020; 152:034105. [DOI: 10.1063/1.5138643] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Affiliation(s)
- Pavel Pokhilko
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, USA
| | - Daniil Izmodenov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, USA
- Department of Chemistry, Lomonosov Moscow State University, Moscow 119991, Russia
| | - Anna I. Krylov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, USA
| |
Collapse
|
24
|
Mester D, Nagy PR, Kállay M. Reduced-Scaling Correlation Methods for the Excited States of Large Molecules: Implementation and Benchmarks for the Second-Order Algebraic-Diagrammatic Construction Approach. J Chem Theory Comput 2019; 15:6111-6126. [DOI: 10.1021/acs.jctc.9b00735] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Dávid Mester
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, H-1521 Budapest, P.O. Box 91, Hungary
| | - Péter R. Nagy
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, H-1521 Budapest, P.O. Box 91, Hungary
| | - Mihály Kállay
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, H-1521 Budapest, P.O. Box 91, Hungary
| |
Collapse
|
25
|
Pavošević F, Hammes-Schiffer S. Multicomponent coupled cluster singles and doubles and Brueckner doubles methods: Proton densities and energies. J Chem Phys 2019; 151:074104. [PMID: 31438716 DOI: 10.1063/1.5116113] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The nuclear-electronic orbital (NEO) framework enables computationally practical coupled cluster calculations of multicomponent molecular systems, in which all electrons and specified nuclei, typically protons, are treated quantum mechanically. In addition to energies, computing accurate proton densities is essential for the calculation of reliable molecular properties, including vibrationally averaged geometries and vibrational frequencies. Herein, the Lagrangian formalism for the multicomponent coupled cluster with single and double excitations (NEO-CCSD) method is derived and implemented. The multicomponent coupled cluster with double excitations method using optimized Brueckner orbitals, denoted as NEO-BCCD, is also developed. Both of these methods are used to compute the proton densities for two molecular systems. The results illustrate that orbital relaxation effects, which can be included either indirectly with the NEO-CCSD method or directly with the NEO-BCCD method, are critical for computing even qualitatively accurate proton densities. Both methods are also able to provide accurate proton affinities and vibrationally averaged optimized geometries. This Lagrangian formalism will enable the calculation of other properties such as analytical nuclear gradients and Hessians with NEO coupled cluster methods. Moreover, the accuracy of these methods may be improved systematically by the inclusion of higher-order excitations. Thus, this work provides the foundation for a wide range of future methodological developments and applications within the NEO framework.
Collapse
Affiliation(s)
- Fabijan Pavošević
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, USA
| | - Sharon Hammes-Schiffer
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, USA
| |
Collapse
|