1
|
Garner SM, Upadhyay S, Li X, Hammes-Schiffer S. Nuclear-Electronic Orbital Time-Dependent Configuration Interaction Method. J Phys Chem Lett 2024; 15:6017-6023. [PMID: 38815051 DOI: 10.1021/acs.jpclett.4c00805] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/01/2024]
Abstract
Combining real-time electronic structure with the nuclear-electronic orbital (NEO) method has enabled the simulation of complex nonadiabatic chemical processes. However, accurate descriptions of hydrogen tunneling and double excitations require multiconfigurational treatments. Herein, we develop and implement the real-time NEO time-dependent configuration interaction (NEO-TDCI) approach. Comparison to NEO-full CI calculations of absorption spectra for a molecular system shows that the NEO-TDCI approach can accurately capture the tunneling splitting associated with the electronic ground state as well as vibronic progressions corresponding to double electron-proton excitations associated with excited electronic states. Both of these features are absent from spectra obtained with single reference real-time NEO methods. Our simulations of hydrogen tunneling dynamics illustrate the oscillation of the proton density from one side to the other via a delocalized, bilobal proton wave function. These results indicate that the NEO-TDCI approach is highly suitable for studying hydrogen tunneling and other inherently multiconfigurational systems.
Collapse
Affiliation(s)
- Scott M Garner
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States
| | - Shiv Upadhyay
- 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, Princeton University, Princeton, New Jersey 08544, United States
| |
Collapse
|
2
|
Dong X, Thompson LM. Time propagation of electronic wavefunctions using nonorthogonal determinant expansions. J Chem Phys 2024; 160:024106. [PMID: 38189613 DOI: 10.1063/5.0179601] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2023] [Accepted: 12/15/2023] [Indexed: 01/09/2024] Open
Abstract
The use of truncated configuration interaction in real-time time-dependent simulations of electron dynamics provides a balance of computational cost and accuracy, while avoiding some of the failures associated with real-time time-dependent density functional theory. However, low-order truncated configuration interaction also has limitations, such as overestimation of polarizability in configuration interaction singles, even when perturbative doubles are included. Increasing the size of the determinant expansion may not be computationally feasible, and so, in this work, we investigate the use of nonorthogonality in the determinant expansion to establish the extent to which higher-order substitutions can be recovered, providing an improved description of electron dynamics. Model systems are investigated to quantify the extent to which different methods accurately reproduce the (hyper)polarizability, including the high-harmonic generation spectrum of H2, water, and butadiene.
Collapse
Affiliation(s)
- Xinju Dong
- Department of Chemistry, University of Louisville, Louisville, Kentucky 40205, USA
| | - Lee M Thompson
- Department of Chemistry, University of Louisville, Louisville, Kentucky 40205, USA
| |
Collapse
|
3
|
Williams-Young DB, Yuwono SH, DePrince III AE, Yang C. Approximate Exponential Integrators for Time-Dependent Equation-of-Motion Coupled Cluster Theory. J Chem Theory Comput 2023; 19:9177-9186. [PMID: 38086060 PMCID: PMC10753770 DOI: 10.1021/acs.jctc.3c00911] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2023] [Revised: 10/09/2023] [Accepted: 10/16/2023] [Indexed: 12/27/2023]
Abstract
With a growing demand for time-domain simulations of correlated many-body systems, the development of efficient and stable integration schemes for the time-dependent Schrödinger equation is of keen interest in modern electronic structure theory. In this work, we present two approaches for the formation of the quantum propagator for time-dependent equation-of-motion coupled cluster theory based on the Chebyshev and Arnoldi expansions of the complex, nonhermitian matrix exponential, respectively. The proposed algorithms are compared with the short-iterative Lanczos method of Cooper et al. [J. Phys. Chem. A 2021 125, 5438-5447], the fourth-order Runge-Kutta method, and exact dynamics for a set of small but challenging test problems. For each of the cases studied, both of the proposed integration schemes demonstrate superior accuracy and efficiency relative to the reference simulations.
Collapse
Affiliation(s)
- David B. Williams-Young
- Applied
Mathematics and Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Stephen H. Yuwono
- Department
of Chemistry and Biochemistry, Florida State
University, Tallahassee, Florida 32306, United States
| | - A. Eugene DePrince III
- Department
of Chemistry and Biochemistry, Florida State
University, Tallahassee, Florida 32306, United States
| | - Chao Yang
- Applied
Mathematics and Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
4
|
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
|
5
|
Langkabel F, Albrecht PA, Bande A, Krause P. Making Optical Excitations Visible - an Exciton Wavefunction Extension to the Time-dependent Configuration Interaction Method. Chem Phys 2022. [DOI: 10.1016/j.chemphys.2022.111502] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
6
|
Durden AS, Levine BG. Floquet Time-Dependent Configuration Interaction for Modeling Ultrafast Electron Dynamics. J Chem Theory Comput 2022; 18:795-806. [DOI: 10.1021/acs.jctc.1c01009] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Affiliation(s)
- Andrew S. Durden
- Department of Chemistry and Institute for Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794, United States
| | - Benjamin G. Levine
- Department of Chemistry and Institute for Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794, United States
| |
Collapse
|
7
|
Cooper BC, Koulias LN, Nascimento DR, Li X, DePrince AE. Short Iterative Lanczos Integration in Time-Dependent Equation-of-Motion Coupled-Cluster Theory. J Phys Chem A 2021; 125:5438-5447. [PMID: 34121405 DOI: 10.1021/acs.jpca.1c01102] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A time-dependent (TD) formulation of equation-of-motion coupled-cluster (EOM-CC) theory can provide excited-state information over an arbitrarily wide energy window with a reduced memory footprint relative to conventional, frequency-domain EOM-CC theory. However, the floating-point costs of the time-integration required by TD-EOM-CC are generally far larger than those of the frequency-domain form of the approach. This work considers the potential of the short iterative Lanczos (SIL) integration scheme [J. Chem. Phys. 1986, 85, 5870-5876] to reduce the floating-point costs of TD-EOM-CC simulations. Low-energy and K-edge absorption features for small molecules are evaluated using TD-EOM-CC with single and double excitations, with the time-integrations carried out via SIL and fourth-order Runge-Kutta (RK4) schemes. Spectra derived from SIL- and RK4-driven simulations are nearly indistinguishable, and with an appropriately chosen subspace dimension, the SIL requires far fewer floating-point operations than are required by RK4. For K-edge spectra, SIL is the more efficient scheme by an average factor of 7.2.
Collapse
Affiliation(s)
- Brandon C Cooper
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, United States
| | - Lauren N Koulias
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Daniel R Nascimento
- Department of Chemistry, The University of Memphis, Memphis, Tennessee 38152, United States
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - A Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, United States
| |
Collapse
|
8
|
Baiardi A. Electron Dynamics with the Time-Dependent Density Matrix Renormalization Group. J Chem Theory Comput 2021; 17:3320-3334. [PMID: 34043347 DOI: 10.1021/acs.jctc.0c01048] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
In this work, we simulate the electron dynamics in molecular systems with the time-dependent density matrix renormalization group (TD-DMRG) algorithm. We leverage the generality of the so-called tangent-space TD-DMRG formulation and design a computational framework in which the dynamics is driven by the exact nonrelativistic electronic Hamiltonian. We show that by parametrizing the wave function as a matrix product state, we can accurately simulate the dynamics of systems including up to 20 electrons and 32 orbitals. We apply the TD-DMRG algorithm to three problems that are hardly targeted by time-independent methods: the calculation of molecular (hyper)polarizabilities, the simulation of electronic absorption spectra, and the study of ultrafast ionization dynamics.
Collapse
Affiliation(s)
- Alberto Baiardi
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, Zürich 8093, Switzerland
| |
Collapse
|
9
|
Woźniak AP, Lesiuk M, Przybytek M, Efimov DK, Prauzner-Bechcicki JS, Mandrysz M, Ciappina M, Pisanty E, Zakrzewski J, Lewenstein M, Moszyński R. A systematic construction of Gaussian basis sets for the description of laser field ionization and high-harmonic generation. J Chem Phys 2021; 154:094111. [PMID: 33685145 DOI: 10.1063/5.0040879] [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 precise understanding of mechanisms governing the dynamics of electrons in atoms and molecules subjected to intense laser fields has a key importance for the description of attosecond processes such as the high-harmonic generation and ionization. From the theoretical point of view, this is still a challenging task, as new approaches to solve the time-dependent Schrödinger equation with both good accuracy and efficiency are still emerging. Until recently, the purely numerical methods of real-time propagation of the wavefunction using finite grids have been frequently and successfully used to capture the electron dynamics in small one- or two-electron systems. However, as the main focus of attoscience shifts toward many-electron systems, such techniques are no longer effective and need to be replaced by more approximate but computationally efficient ones. In this paper, we explore the increasingly popular method of expanding the wavefunction of the examined system into a linear combination of atomic orbitals and present a novel systematic scheme for constructing an optimal Gaussian basis set suitable for the description of excited and continuum atomic or molecular states. We analyze the performance of the proposed basis sets by carrying out a series of time-dependent configuration interaction calculations for the hydrogen atom in fields of intensity varying from 5 × 1013 W/cm2 to 5 × 1014 W/cm2. We also compare the results with the data obtained using Gaussian basis sets proposed previously by other authors.
Collapse
Affiliation(s)
| | - Michał Lesiuk
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| | - Michał Przybytek
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| | - Dmitry K Efimov
- Institute of Theoretical Physics, Jagiellonian University in Krakow, Łojasiewicza 11, 30-348 Kraków, Poland
| | - Jakub S Prauzner-Bechcicki
- Marian Smoluchowski Institute of Physics, Jagiellonian University in Krakow, Łojasiewicza 11, 30-348 Kraków, Poland
| | - Michał Mandrysz
- Institute of Theoretical Physics, Jagiellonian University in Krakow, Łojasiewicza 11, 30-348 Kraków, Poland
| | - Marcelo Ciappina
- ICFO-Institut de Ciencies Fotoniques, The Barcelona Institute of Science and Technology, Av. Carl Friedrich Gauss 3, 08860, Castelldefels, Barcelona, Spain
| | - Emilio Pisanty
- ICFO-Institut de Ciencies Fotoniques, The Barcelona Institute of Science and Technology, Av. Carl Friedrich Gauss 3, 08860, Castelldefels, Barcelona, Spain
| | - Jakub Zakrzewski
- Institute of Theoretical Physics, Jagiellonian University in Krakow, Łojasiewicza 11, 30-348 Kraków, Poland
| | - Maciej Lewenstein
- ICFO-Institut de Ciencies Fotoniques, The Barcelona Institute of Science and Technology, Av. Carl Friedrich Gauss 3, 08860, Castelldefels, Barcelona, Spain
| | - Robert Moszyński
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| |
Collapse
|
10
|
Li X, Govind N, Isborn C, DePrince AE, Lopata K. Real-Time Time-Dependent Electronic Structure Theory. Chem Rev 2020; 120:9951-9993. [DOI: 10.1021/acs.chemrev.0c00223] [Citation(s) in RCA: 64] [Impact Index Per Article: 16.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195-1700, United States
| | - Niranjan Govind
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Christine Isborn
- Department of Chemistry and Chemical Biology, University of California, Merced, California 95343, United States
| | - A. Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, United States
| | - Kenneth Lopata
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, United States
| |
Collapse
|
11
|
Kowalski K, Bauman NP. Sub-system quantum dynamics using coupled cluster downfolding techniques. J Chem Phys 2020; 152:244127. [DOI: 10.1063/5.0008436] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Karol Kowalski
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| | - Nicholas P. Bauman
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| |
Collapse
|
12
|
Baiardi A, Reiher M. The density matrix renormalization group in chemistry and molecular physics: Recent developments and new challenges. J Chem Phys 2020; 152:040903. [DOI: 10.1063/1.5129672] [Citation(s) in RCA: 85] [Impact Index Per Article: 21.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023] Open
Affiliation(s)
- Alberto Baiardi
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Markus Reiher
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
13
|
Nascimento DR, DePrince AE. A general time-domain formulation of equation-of-motion coupled-cluster theory for linear spectroscopy. J Chem Phys 2019; 151:204107. [DOI: 10.1063/1.5125494] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Affiliation(s)
- Daniel R. Nascimento
- 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
|
14
|
Baiardi A, Reiher M. Large-Scale Quantum Dynamics with Matrix Product States. J Chem Theory Comput 2019; 15:3481-3498. [DOI: 10.1021/acs.jctc.9b00301] [Citation(s) in RCA: 48] [Impact Index Per Article: 9.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Affiliation(s)
- Alberto Baiardi
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Markus Reiher
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
15
|
Liu H, Jenkins AJ, Wildman A, Frisch MJ, Lipparini F, Mennucci B, Li X. Time-Dependent Complete Active Space Embedded in a Polarizable Force Field. J Chem Theory Comput 2019; 15:1633-1641. [DOI: 10.1021/acs.jctc.8b01152] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Hongbin Liu
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Andrew J. Jenkins
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Andrew Wildman
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Michael J. Frisch
- Gaussian Inc., 340 Quinnipiac Street, Building 40, Wallingford, Connecticut 06492, United States
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via Risorgimento 35, 56126 Pisa, Italy
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via Risorgimento 35, 56126 Pisa, Italy
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
16
|
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
|
17
|
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
|
18
|
Ulusoy IS, Stewart Z, Wilson AK. The role of the CI expansion length in time-dependent studies. J Chem Phys 2018; 148:014107. [DOI: 10.1063/1.5004412] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Inga S. Ulusoy
- Department of Chemistry, Michigan State University, 578 S Shaw Lane, East Lansing, Michigan 48824-1322, USA
| | - Zachary Stewart
- Department of Chemistry, Michigan State University, 578 S Shaw Lane, East Lansing, Michigan 48824-1322, USA
| | - Angela K. Wilson
- Department of Chemistry, Michigan State University, 578 S Shaw Lane, East Lansing, Michigan 48824-1322, USA
| |
Collapse
|