51
|
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
|
52
|
Williams‐Young DB, Petrone A, Sun S, Stetina TF, Lestrange P, Hoyer CE, Nascimento DR, Koulias L, Wildman A, Kasper J, Goings JJ, Ding F, DePrince AE, Valeev EF, Li X. The Chronus Quantum software package. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2019. [DOI: 10.1002/wcms.1436] [Citation(s) in RCA: 34] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
- David B. Williams‐Young
- Computational Research Division Lawrence Berkeley National Laboratory Berkeley California
- Department of Chemistry University of Washington Seattle Washington
| | - Alessio Petrone
- Dipartimento di Scienze Chimiche Università di Napoli “Federico II”, Complesso Universitario di M.S. Angelo Naples Italy
| | - Shichao Sun
- Department of Chemistry University of Washington Seattle Washington
| | - Torin F. Stetina
- Department of Chemistry University of Washington Seattle Washington
| | | | - Chad E. Hoyer
- Department of Chemistry University of Washington Seattle Washington
| | - Daniel R. Nascimento
- Department of Chemistry and Biochemistry Florida State University Tallahassee Florida
| | - Lauren Koulias
- Department of Chemistry University of Washington Seattle Washington
| | - Andrew Wildman
- Department of Chemistry University of Washington Seattle Washington
| | - Joseph Kasper
- Department of Chemistry University of Washington Seattle Washington
| | | | - Feizhi Ding
- Division of Chemistry and Chemical Engineering California Institute of Technology Pasadena California
| | - A. Eugene DePrince
- Department of Chemistry and Biochemistry Florida State University Tallahassee Florida
| | | | - Xiaosong Li
- Department of Chemistry University of Washington Seattle Washington
| |
Collapse
|
53
|
Neville SP, Schuurman MS. Efficient calculation of X-ray absorption spectra using Chebyshev-Slepian filter diagonalisation. J Chem Phys 2019; 150:184115. [DOI: 10.1063/1.5092975] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023] Open
Affiliation(s)
- Simon P. Neville
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, 10 Marie Curie, Ottawa, Ontario K1N 6N5, Canada
| | - Michael S. Schuurman
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, 10 Marie Curie, Ottawa, Ontario K1N 6N5, Canada
- National Research Council of Canada, 100 Sussex Drive, Ottawa, Ontario K1A 0R6, Canada
| |
Collapse
|
54
|
Pedersen TB, Kvaal S. Symplectic integration and physical interpretation of time-dependent coupled-cluster theory. J Chem Phys 2019; 150:144106. [DOI: 10.1063/1.5085390] [Citation(s) in RCA: 32] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023] Open
Affiliation(s)
- Thomas Bondo Pedersen
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway
| | - Simen Kvaal
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway
| |
Collapse
|
55
|
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
|
56
|
Maradzike E, DePrince AE. Modeling core-level excitations with variationally optimized reduced-density matrices and the extended random phase approximation. J Chem Phys 2018; 149:234101. [PMID: 30579305 DOI: 10.1063/1.5048924] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The information contained within ground-state one- and two-electron reduced-density matrices (RDMs) can be used to compute wave functions and energies for electronically excited states through the extended random phase approximation (ERPA). The ERPA is an appealing framework for describing excitations out of states obtained via the variational optimization of the two-electron RDM (2-RDM), as the variational 2-RDM (v2RDM) approach itself can only be used to describe the lowest-energy state of a given spin symmetry. The utility of the ERPA for predicting near-edge features relevant to x-ray absorption spectroscopy is assessed for the case that the 2-RDM is obtained from a ground-state v2RDM-driven complete active space self-consistent field (CASSCF) computation. A class of killer conditions for the CASSCF-specific ERPA excitation operator is derived, and it is demonstrated that a reliable description of core-level excitations requires an excitation operator that fulfills these conditions; the core-valence separation (CVS) scheme yields such an operator. Absolute excitation energies evaluated within the CASSCF/CVS-ERPA framework are slightly more accurate than those obtained from the usual random phase approximation (RPA), but the CVS-ERPA is not more accurate than RPA for predicting the relative positions of near-edge features. Nonetheless, CVS-ERPA is established as a reasonable starting point for the treatment of core-level excitations using variationally optimized 2-RDMs.
Collapse
Affiliation(s)
- Elvis Maradzike
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, USA
| | - A Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, USA
| |
Collapse
|
57
|
Konecny L, Kadek M, Komorovsky S, Ruud K, Repisky M. Resolution-of-identity accelerated relativistic two- and four-component electron dynamics approach to chiroptical spectroscopies. J Chem Phys 2018; 149:204104. [PMID: 30501232 DOI: 10.1063/1.5051032] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present an implementation and application of electron dynamics based on real-time time-dependent density functional theory (RT-TDDFT) and relativistic 2-component X2C and 4-component Dirac-Coulomb (4c) Hamiltonians to the calculation of electron circular dichroism and optical rotatory dispersion spectra. In addition, the resolution-of-identity approximation for the Coulomb term (RI-J) is introduced into RT-TDDFT and formulated entirely in terms of complex quaternion algebra. The proposed methodology was assessed on the dimethylchalcogenirane series, C4H8X (X = O, S, Se, Te, Po, Lv), and the spectra obtained by non-relativistic and relativistic methods start to disagree for Se and Te, while dramatic differences are observed for Po and Lv. The X2C approach, even in its simplest one-particle form, reproduces the reference 4c results surprisingly well across the entire series while offering an 8-fold speed-up of the simulations. An overall acceleration of RT-TDDFT by means of X2C and RI-J increases with system size and approaches a factor of almost 25 when compared to the full 4c treatment, without compromising the accuracy of the final spectra. These results suggest that one-particle X2C electron dynamics with RI-J acceleration is an attractive method for the calculation of chiroptical spectra in the valence region.
Collapse
Affiliation(s)
- Lukas Konecny
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Tromsø-The Arctic University of Norway, Tromsø, Norway
| | - Marius Kadek
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Tromsø-The Arctic University of Norway, Tromsø, Norway
| | - Stanislav Komorovsky
- Institute of Inorganic Chemistry, Slovak Academy of Sciences, Bratislava, Slovakia
| | - Kenneth Ruud
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Tromsø-The Arctic University of Norway, Tromsø, Norway
| | - Michal Repisky
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Tromsø-The Arctic University of Norway, Tromsø, Norway
| |
Collapse
|
58
|
Neville SP, Schuurman MS. A general approach for the calculation and characterization of x-ray absorption spectra. J Chem Phys 2018; 149:154111. [DOI: 10.1063/1.5048520] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Affiliation(s)
- Simon P. Neville
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, 10 Marie Curie, Ottawa, Ontario K1N 6N5, Canada
| | - Michael S. Schuurman
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, 10 Marie Curie, Ottawa, Ontario K1N 6N5, Canada
- National Research Council of Canada, 100 Sussex Drive, Ottawa, Ontario K1A 0R6, Canada
| |
Collapse
|
59
|
White AF, Chan GKL. A Time-Dependent Formulation of Coupled-Cluster Theory for Many-Fermion Systems at Finite Temperature. J Chem Theory Comput 2018; 14:5690-5700. [DOI: 10.1021/acs.jctc.8b00773] [Citation(s) in RCA: 39] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Alec F. White
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Garnet Kin-Lic Chan
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| |
Collapse
|
60
|
Peng WT, Fales BS, Levine BG. Simulating Electron Dynamics of Complex Molecules with Time-Dependent Complete Active Space Configuration Interaction. J Chem Theory Comput 2018; 14:4129-4138. [PMID: 29986143 DOI: 10.1021/acs.jctc.8b00381] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Time-dependent electronic structure methods are growing in popularity as tools for modeling ultrafast and/or nonlinear processes, for computing spectra, and as the electronic structure component of mean-field molecular dynamics simulations. Time-dependent configuration interaction (TD-CI) offers several advantages over the widely used real-time time-dependent density functional theory: namely, that it correctly models Rabi oscillations; it offers a spin-pure description of open-shell systems; and a hierarchy of TD-CI methods can be defined that systematically approach the exact solution of the time-dependent Schrodinger equation (TDSE). In this work, we present a novel TD-CI approach that extends TD-CI to large complete active-space configuration expansions. Such extension is enabled by use of a direct configuration interaction approach that eliminates the need to explicitly build, store, or diagonalize the Hamiltonian matrix. Graphics processing unit (GPU) acceleration enables fast solution of the TDSE even for large active spaces-up to 12 electrons in 12 orbitals (853776 determinants) in this work. A symplectic split operator propagator yields long-time norm conservation. We demonstrate the applicability of our approach by computing the response of a large molecule with a strongly correlated ground state, decacene (C42H24), to various pulses (δ-function, transform limited, chirped). Our simulations predict that chirped pulses can be used to induce dipole-forbidden transitions. Simulations of decacene using the 6-31G(d) basis set and a 12 electrons/12 orbitals active space took 20.1 h to propagate for 100 fs with a 1 attosecond time step on a single NVIDIA K40 GPU. Convergence with respect to time step is found to depend on the property being computed and the chosen active space.
Collapse
Affiliation(s)
- Wei-Tao Peng
- Department of Chemistry , Michigan State University , East Lansing , Michigan 48824 , United States
| | - B Scott Fales
- Department of Chemistry and the PULSE Institute , Stanford University , Stanford , California 94305 , United States.,SLAC National Accelerator Laboratory , Menlo Park , California 94025 , United States
| | - Benjamin G Levine
- Department of Chemistry , Michigan State University , East Lansing , Michigan 48824 , United States
| |
Collapse
|
61
|
Abstract
This report presents selected highlights from 2017 final birth data on key demographic, health care utilization, and infant health indicators. General fertility rates (the number of births per 1,000 females aged 15-44 years) and teen birth rates are presented by race and Hispanic origin. The use of Medicaid as the source of payment for the delivery and preterm birth rates are presented by the age of the mother. Data for 2017 are compared with 2016 for each indicator.
Collapse
|
62
|
Norman P, Dreuw A. Simulating X-ray Spectroscopies and Calculating Core-Excited States of Molecules. Chem Rev 2018; 118:7208-7248. [DOI: 10.1021/acs.chemrev.8b00156] [Citation(s) in RCA: 164] [Impact Index Per Article: 27.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Affiliation(s)
- Patrick Norman
- Department of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, SE-106 91 Stockholm, Sweden
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| |
Collapse
|
63
|
Sato T, Pathak H, Orimo Y, Ishikawa KL. Communication: Time-dependent optimized coupled-cluster method for multielectron dynamics. J Chem Phys 2018; 148:051101. [DOI: 10.1063/1.5020633] [Citation(s) in RCA: 47] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023] Open
Affiliation(s)
- Takeshi Sato
- Photon Science Center, School of Engineering, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8656, Japan
- Department of Nuclear Engineering and Management, School of Engineering, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8656, Japan
| | - Himadri Pathak
- Department of Nuclear Engineering and Management, School of Engineering, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8656, Japan
| | - Yuki Orimo
- Department of Nuclear Engineering and Management, School of Engineering, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8656, Japan
| | - Kenichi L. Ishikawa
- Photon Science Center, School of Engineering, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8656, Japan
- Department of Nuclear Engineering and Management, School of Engineering, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8656, Japan
| |
Collapse
|
64
|
Goings JJ, Lestrange PJ, Li X. Real‐time time‐dependent electronic structure theory. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2017. [DOI: 10.1002/wcms.1341] [Citation(s) in RCA: 84] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
| | | | - Xiaosong Li
- Department of ChemistryUniversity of Washington Seattle WA USA
| |
Collapse
|
65
|
Koh KJ, Nguyen-Beck TS, Parkhill J. Accelerating Realtime TDDFT with Block-Orthogonalized Manby–Miller Embedding Theory. J Chem Theory Comput 2017; 13:4173-4178. [DOI: 10.1021/acs.jctc.7b00494] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Kevin J. Koh
- Department of Chemistry and
Biochemistry, The University of Notre Dame, 251 Nieuwland Science Hall, Notre Dame, Indiana 46556, United States
| | - Triet S. Nguyen-Beck
- Department of Chemistry and
Biochemistry, The University of Notre Dame, 251 Nieuwland Science Hall, Notre Dame, Indiana 46556, United States
| | - John Parkhill
- Department of Chemistry and
Biochemistry, The University of Notre Dame, 251 Nieuwland Science Hall, Notre Dame, Indiana 46556, United States
| |
Collapse
|
66
|
Mester D, Nagy PR, Kállay M. Reduced-cost linear-response CC2 method based on natural orbitals and natural auxiliary functions. J Chem Phys 2017; 146:194102. [PMID: 28527453 DOI: 10.1063/1.4983277] [Citation(s) in RCA: 43] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A reduced-cost density fitting (DF) linear-response second-order coupled-cluster (CC2) method has been developed for the evaluation of excitation energies. The method is based on the simultaneous truncation of the molecular orbital (MO) basis and the auxiliary basis set used for the DF approximation. For the reduction of the size of the MO basis, state-specific natural orbitals (NOs) are constructed for each excited state using the average of the second-order Møller-Plesset (MP2) and the corresponding configuration interaction singles with perturbative doubles [CIS(D)] density matrices. After removing the NOs of low occupation number, natural auxiliary functions (NAFs) are constructed [M. Kállay, J. Chem. Phys. 141, 244113 (2014)], and the NAF basis is also truncated. Our results show that, for a triple-zeta basis set, about 60% of the virtual MOs can be dropped, while the size of the fitting basis can be reduced by a factor of five. This results in a dramatic reduction of the computational costs of the solution of the CC2 equations, which are in our approach about as expensive as the evaluation of the MP2 and CIS(D) density matrices. All in all, an average speedup of more than an order of magnitude can be achieved at the expense of a mean absolute error of 0.02 eV in the calculated excitation energies compared to the canonical CC2 results. Our benchmark calculations demonstrate that the new approach enables the efficient computation of CC2 excitation energies for excited states of all types of medium-sized molecules composed of up to 100 atoms with triple-zeta quality basis sets.
Collapse
Affiliation(s)
- Dávid Mester
- MTA-BME Lendület Quantum Chemistry Research Group, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Péter R Nagy
- MTA-BME Lendület Quantum Chemistry Research Group, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Mihály Kállay
- MTA-BME Lendület Quantum Chemistry Research Group, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| |
Collapse
|
67
|
Baudin P, Kjærgaard T, Kristensen K. CC2 oscillator strengths within the local framework for calculating excitation energies (LoFEx). J Chem Phys 2017; 146:144107. [PMID: 28411600 DOI: 10.1063/1.4979713] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
In a recent work [P. Baudin and K. Kristensen, J. Chem. Phys. 144, 224106 (2016)], we introduced a local framework for calculating excitation energies (LoFEx), based on second-order approximated coupled cluster (CC2) linear-response theory. LoFEx is a black-box method in which a reduced excitation orbital space (XOS) is optimized to provide coupled cluster (CC) excitation energies at a reduced computational cost. In this article, we present an extension of the LoFEx algorithm to the calculation of CC2 oscillator strengths. Two different strategies are suggested, in which the size of the XOS is determined based on the excitation energy or the oscillator strength of the targeted transitions. The two strategies are applied to a set of medium-sized organic molecules in order to assess both the accuracy and the computational cost of the methods. The results show that CC2 excitation energies and oscillator strengths can be calculated at a reduced computational cost, provided that the targeted transitions are local compared to the size of the molecule. To illustrate the potential of LoFEx for large molecules, both strategies have been successfully applied to the lowest transition of the bivalirudin molecule (4255 basis functions) and compared with time-dependent density functional theory.
Collapse
Affiliation(s)
- Pablo Baudin
- qLEAP Center for Theoretical Chemistry, Department of Chemistry, Aarhus University, Langelandsgade 140, DK-8000 Aarhus C, Denmark
| | - Thomas Kjærgaard
- qLEAP Center for Theoretical Chemistry, Department of Chemistry, Aarhus University, Langelandsgade 140, DK-8000 Aarhus C, Denmark
| | - Kasper Kristensen
- qLEAP Center for Theoretical Chemistry, Department of Chemistry, Aarhus University, Langelandsgade 140, DK-8000 Aarhus C, Denmark
| |
Collapse
|