1
|
Mehmood A, Silfies MC, Durden AS, Allison TK, Levine BG. Simulating ultrafast transient absorption spectra from first principles using a time-dependent configuration interaction probe. J Chem Phys 2024; 161:044107. [PMID: 39041880 DOI: 10.1063/5.0215890] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2024] [Accepted: 06/30/2024] [Indexed: 07/24/2024] Open
Abstract
Transient absorption spectroscopy (TAS) is among the most common ultrafast photochemical experiments, but its interpretation remains challenging. In this work, we present an efficient and robust method for simulating TAS signals from first principles. Excited-state absorption and stimulated emission (SE) signals are computed using time-dependent complete active space configuration interaction (TD-CASCI) simulations, leveraging the robustness of time-domain simulation to minimize electronic structure failure. We demonstrate our approach by simulating the TAS signal of 1'-hydroxy-2'-acetonapthone (HAN) from ab initio multiple spawning nonadiabatic molecular dynamics simulations. Our results are compared to gas-phase TAS data recorded from both jet-cooled (T ∼ 40 K) and hot (∼403 K) molecules via cavity-enhanced TAS (CE-TAS). Decomposition of the computed spectrum allows us to assign a rise in the SE signal to excited-state proton transfer and the ultimate decay of the signal to relaxation through a twisted conical intersection. The total cost of computing the observable signal (∼1700 graphics processing unit hours for ∼4 ns of electron dynamics) was markedly less than that of performing the ab initio multiple spawning calculations used to compute the underlying nonadiabatic dynamics.
Collapse
Affiliation(s)
- Arshad Mehmood
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
- Institute for Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794, USA
| | - Myles C Silfies
- Department of Physics, Stony Brook University, Stony Brook, New York 11794, USA
| | - Andrew S Durden
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
- Institute for Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794, USA
| | - Thomas K Allison
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
- Department of Physics, Stony Brook University, Stony Brook, New York 11794, USA
| | - Benjamin G Levine
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
- Institute for Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794, USA
| |
Collapse
|
2
|
Bartlett RJ. Perspective on Coupled-cluster Theory. The evolution toward simplicity in quantum chemistry. Phys Chem Chem Phys 2024; 26:8013-8037. [PMID: 38390989 DOI: 10.1039/d3cp03853j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/24/2024]
Abstract
Coupled-cluster theory has revolutionized quantum chemistry. It has provided the framework to effectively solve the problem of electron correlation, the main focus of the field for over 60 years. This has enabled ab initio quantum chemistry to provide predictive quality results for most quantities of interest that are obtainable from first-principle calculations. The best that one can do in a basis is the 'full CI,' the exact solution of the non-relativistic Schrödinger equation or, if need be, the relativistic Dirac equation. With due regard to converging the basis set and adequate consideration of higher clusters and relativity in a calculation, virtually predictive results can be obtained. But in addition to its numerical performance, coupled-cluster theory also offers a conceptually new, many-body foundation for the theory that should be appreciated by all practitioners. The latter is emphasized in this perspective, leading to the 'evolution toward simplicity' in the title. The ultimate theory will benefit from the several features that are uniquely exact in coupled-cluster theory and its equation-of-motion (EOM-CC) extensions.
Collapse
Affiliation(s)
- Rodney J Bartlett
- Quantum Theory Project, Department of Chemistry, University of Florida, P. O. Box 117200, Gainesville, Florida, USA.
| |
Collapse
|
3
|
Crisci L, Coppola F, Petrone A, Rega N. Tuning ultrafast time-evolution of photo-induced charge-transfer states: A real-time electronic dynamics study in substituted indenotetracene derivatives. J Comput Chem 2024; 45:210-221. [PMID: 37706600 DOI: 10.1002/jcc.27231] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2023] [Revised: 08/31/2023] [Accepted: 09/05/2023] [Indexed: 09/15/2023]
Abstract
Photo-induced charge transfer (CT) states are pivotal in many technological and biological processes. A deeper knowledge of such states is mandatory for modeling the charge migration dynamics. Real-time time-dependent density functional theory (RT-TD-DFT) electronic dynamics simulations are employed to explicitly observe the electronic density time-evolution upon photo-excitation. Asymmetrically substituted indenotetracene molecules, given their potential application as n-type semiconductors in organic photovoltaic materials, are here investigated. Effects of substituents with different electron-donating characters are analyzed in terms of the overall electronic energy spacing and resulting ultrafast CT dynamics through linear response (LR-)TD-DFT and RT-TD-DFT based approaches. The combination of the computational techniques here employed provided direct access to the electronic density reorganization in time and to its spatial and rational representation in terms of molecular orbital occupation time evolution. Such results can be exploited to design peculiar directional charge dynamics, crucial when photoactive materials are used for light-harvesting applications.
Collapse
Affiliation(s)
- Luigi Crisci
- Department of Chemical Sciences, University of Napoli Federico II, Complesso Universitario di M.S. Angelo, Naples, Italy
- Scuola Normale Superiore di Pisa, Pisa, Italy
| | | | - Alessio Petrone
- Department of Chemical Sciences, University of Napoli Federico II, Complesso Universitario di M.S. Angelo, Naples, Italy
- Scuola Superiore Meridionale, Naples, Italy
- Istituto Nazionale Di Fisica Nucleare, Sezione di Napoli, Complesso Universitario di M.S. Angelo ed. 6, Naples, Italy
| | - Nadia Rega
- Department of Chemical Sciences, University of Napoli Federico II, Complesso Universitario di M.S. Angelo, Naples, Italy
- Scuola Superiore Meridionale, Naples, Italy
- Istituto Nazionale Di Fisica Nucleare, Sezione di Napoli, Complesso Universitario di M.S. Angelo ed. 6, Naples, Italy
| |
Collapse
|
4
|
Zhang C, Zheng X, Liu J, Asthana A, Cheng L. Analytic gradients for relativistic exact-two-component equation-of-motion coupled-cluster singles and doubles method. J Chem Phys 2023; 159:244113. [PMID: 38153147 DOI: 10.1063/5.0175041] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2023] [Accepted: 11/20/2023] [Indexed: 12/29/2023] Open
Abstract
A first implementation of analytic gradients for spinor-based relativistic equation-of-motion coupled-cluster singles and doubles method using an exact two-component Hamiltonian augmented with atomic mean-field spin-orbit integrals is reported. To demonstrate its applicability, we present calculations of equilibrium structures and harmonic vibrational frequencies for the electronic ground and excited states of the radium mono-amide molecule (RaNH2) and the radium mono-methoxide molecule (RaOCH3). Spin-orbit coupling is shown to quench Jahn-Teller effects in the first excited state of RaOCH3, resulting in a C3v equilibrium structure. The calculations also show that the radium atoms in these molecules serve as efficient optical cycling centers.
Collapse
Affiliation(s)
- Chaoqun Zhang
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Xuechen Zheng
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Junzi Liu
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Ayush Asthana
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Lan Cheng
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, USA
| |
Collapse
|
5
|
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
|
6
|
Buttarazzi E, Perrella F, Rega N, Petrone A. Watching the Interplay between Photoinduced Ultrafast Charge Dynamics and Nuclear Vibrations. J Chem Theory Comput 2023; 19:8751-8766. [PMID: 37991892 PMCID: PMC10720350 DOI: 10.1021/acs.jctc.3c00855] [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/04/2023] [Revised: 11/02/2023] [Accepted: 11/03/2023] [Indexed: 11/24/2023]
Abstract
Here is presented the ultrafast hole-electron dynamics of photoinduced metal to ligand charge-transfer (MLCT) states in a Ru(II) complex, [Ru(dcbpy)2(NCS)2]4- (dcbpy = 4,4'-dicarboxy-2,2'-bipyridine), a photoactive molecule employed in dye sensitized solar cells. Via cutting-edge computational techniques, a tailored computational protocol is here presented and developed to provide a detailed analysis of the electronic manifold coupled with nuclear vibrations to better understand the nonradiative pathways and the resulting overall dye performances in light-harvesting processes (electron injection). Thus, the effects of different vibrational modes were investigated on both the electronic levels and charge transfer dynamics through a theoretical-computational approach. First, the linear response time-dependent density functional (LR-TDDFT) formalism was employed to characterize excitation energies and spacing among electronic levels (the electronic layouts). Then, to understand the ultrafast (femtosecond) charge dynamics on the molecular scale, we relied on the nonperturbative mean-field quantum electronic dynamics via real-time (RT-) TDDFT. Three vibrational modes were selected, representative for collective nuclear movements that can have a significant influence on the electronic structure: two involving NCS- ligands and one involving dcbpy ligands. As main results, we observed that such MLCT states, under vibrational distortions, are strongly affected and a faster interligand electron transfer mechanism is observed along with an increasing MLCT character of the adiabatic electronic states approaching closer in energy due to the vibrations. Such findings can help both in providing a molecular picture of multidimensional vibro-electronic spectroscopic techniques, used to characterize ultrafast coherent and noncoherent dynamics of complex systems, and to improve dye performances with particular attention to the study of energy or charge transport processes and vibronic couplings.
Collapse
Affiliation(s)
- Edoardo Buttarazzi
- Scuola
Superiore Meridionale, Largo San Marcellino 10, I-80138 Napoli, Italy
- Department
of Chemical Sciences, University of Napoli
Federico II, Complesso Universitario di Monte S. Angelo, Via Cintia 21, I-80126 Napoli, Italy
| | - Fulvio Perrella
- Scuola
Superiore Meridionale, Largo San Marcellino 10, I-80138 Napoli, Italy
| | - Nadia Rega
- Scuola
Superiore Meridionale, Largo San Marcellino 10, I-80138 Napoli, Italy
- Department
of Chemical Sciences, University of Napoli
Federico II, Complesso Universitario di Monte S. Angelo, Via Cintia 21, I-80126 Napoli, Italy
- Istituto
Nazionale Di Fisica Nucleare, sezione di Napoli, Complesso Universitario
di Monte S. Angelo ed. 6, Via Cintia, I-80126 Napoli, Italy
| | - Alessio Petrone
- Scuola
Superiore Meridionale, Largo San Marcellino 10, I-80138 Napoli, Italy
- Department
of Chemical Sciences, University of Napoli
Federico II, Complesso Universitario di Monte S. Angelo, Via Cintia 21, I-80126 Napoli, Italy
- Istituto
Nazionale Di Fisica Nucleare, sezione di Napoli, Complesso Universitario
di Monte S. Angelo ed. 6, Via Cintia, I-80126 Napoli, Italy
| |
Collapse
|
7
|
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
|
8
|
Jensen AB, Højlund MG, Zoccante A, Madsen NK, Christiansen O. Efficient time-dependent vibrational coupled cluster computations with time-dependent basis sets at the two-mode coupling level: Full and hybrid TDMVCC[2]. J Chem Phys 2023; 159:204106. [PMID: 38010335 DOI: 10.1063/5.0175506] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2023] [Accepted: 11/05/2023] [Indexed: 11/29/2023] Open
Abstract
The computation of the nuclear quantum dynamics of molecules is challenging, requiring both accuracy and efficiency to be applicable to systems of interest. Recently, theories have been developed for employing time-dependent basis functions (denoted modals) with vibrational coupled cluster theory (TDMVCC). The TDMVCC method was introduced along with a pilot implementation, which illustrated good accuracy in benchmark computations. In this paper, we report an efficient implementation of TDMVCC, covering the case where the wave function and Hamiltonian contain up to two-mode couplings. After a careful regrouping of terms, the wave function can be propagated with a cubic computational scaling with respect to the number of degrees of freedom. We discuss the use of a restricted set of active one-mode basis functions for each mode, as well as two interesting limits: (i) the use of a full active basis where the variational modal determination amounts essentially to the variational determination of a time-dependent reference state for the cluster expansion; and (ii) the use of a single function as an active basis for some degrees of freedom. The latter case defines a hybrid TDMVCC/TDH (time-dependent Hartree) approach that can obtain even lower computational scaling. The resulting computational scaling for hybrid and full TDMVCC[2] is illustrated for polyaromatic hydrocarbons with up to 264 modes. Finally, computations on the internal vibrational redistribution of benzoic acid (39 modes) are used to show the faster convergence of TDMVCC/TDH hybrid computations towards TDMVCC compared to simple neglect of some degrees of freedom.
Collapse
Affiliation(s)
| | - Mads Greisen Højlund
- Department of Chemistry, Aarhus University, Langelandsgade 140, DK-8000 Aarhus C, Denmark
| | - Alberto Zoccante
- Dipartimento di Scienze e Innovazione Tecnologica, Università del Piemonte Orientale (UPO), Via T. Michel 11, 15100 Alessandria, Italy
| | - Niels Kristian Madsen
- Department of Chemistry, Aarhus University, Langelandsgade 140, DK-8000 Aarhus C, Denmark
| | - Ove Christiansen
- Department of Chemistry, Aarhus University, Langelandsgade 140, DK-8000 Aarhus C, Denmark
| |
Collapse
|
9
|
Balbi A, Skeidsvoll AS, Koch H. Coupled Cluster Simulation of Impulsive Stimulated X-ray Raman Scattering. J Phys Chem A 2023; 127:8676-8684. [PMID: 37812082 PMCID: PMC10591507 DOI: 10.1021/acs.jpca.3c03678] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2023] [Revised: 08/22/2023] [Indexed: 10/10/2023]
Abstract
Time-dependent equation-of-motion coupled cluster (TD-EOM-CC) is used to simulate impulsive stimulated X-ray Raman scattering (ISXRS) of ultrashort laser pulses by neon, carbon monoxide, pyrrole, and p-aminophenol. The TD-EOM-CC equations are expressed in the basis of field-free EOM-CC states, where the calculation of the core-excited states is simplified through the use of the core-valence separation (CVS) approximation. The transfer of electronic population from the ground state to the core- and valence-excited states is calculated for different numbers of included core- and valence-excited states, as well as for electric field pulses with different polarizations and carrier frequencies. The results indicate that Gaussian pulses can transfer significant electronic populations to the valence states through the Raman process. The sensitivity of this population transfer to the model parameters is analyzed. The time-dependent electronic density for p-aminophenol is also showcased, supporting the interpretation that ISXRS involves localized core excitations and can be used to rapidly generate valence wavepackets.
Collapse
Affiliation(s)
- Alice Balbi
- Scuola
Normale Superiore, Piazza dei Cavalieri, 7, I-56126 Pisa, Italy
| | - Andreas S. Skeidsvoll
- Department
of Chemistry, Norwegian University of Science
and Technology, 7491 Trondheim, Norway
| | - Henrik Koch
- Scuola
Normale Superiore, Piazza dei Cavalieri, 7, I-56126 Pisa, Italy
- Department
of Chemistry, Norwegian University of Science
and Technology, 7491 Trondheim, Norway
| |
Collapse
|
10
|
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
|
11
|
Yuwono SH, Cooper BC, Zhang T, Li X, DePrince AE. Time-dependent equation-of-motion coupled-cluster simulations with a defective Hamiltonian. J Chem Phys 2023; 159:044113. [PMID: 37497820 DOI: 10.1063/5.0157852] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2023] [Accepted: 07/11/2023] [Indexed: 07/28/2023] Open
Abstract
Simulations of laser-induced electron dynamics in a molecular system are performed using time-dependent (TD) equation-of-motion (EOM) coupled-cluster (CC) theory. The target system has been chosen to highlight potential shortcomings of truncated TD-EOM-CC methods [represented in this work by TD-EOM-CC with single and double excitations (TD-EOM-CCSD)], where unphysical spectroscopic features can emerge. Specifically, we explore driven resonant electronic excitations in magnesium fluoride in the proximity of an avoided crossing. Near the avoided crossing, the CCSD similarity-transformed Hamiltonian is defective, meaning that it has complex eigenvalues, and oscillator strengths may take on negative values. When an external field is applied to drive transitions to states exhibiting these traits, unphysical dynamics are observed. For example, the stationary states that make up the time-dependent state acquire populations that can be negative, exceed one, or even complex-valued.
Collapse
Affiliation(s)
- Stephen H Yuwono
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, USA
| | - Brandon C Cooper
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, USA
| | - Tianyuan Zhang
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - A Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, USA
| |
Collapse
|
12
|
Pathak H, Panyala A, Peng B, Bauman NP, Mutlu E, Rehr JJ, Vila FD, Kowalski K. Real-Time Equation-of-Motion Coupled-Cluster Cumulant Green's Function Method: Heterogeneous Parallel Implementation Based on the Tensor Algebra for Many-Body Methods Infrastructure. J Chem Theory Comput 2023; 19:2248-2257. [PMID: 37096369 DOI: 10.1021/acs.jctc.3c00045] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/26/2023]
Abstract
We report the implementation of the real-time equation-of-motion coupled-cluster (RT-EOM-CC) cumulant Green's function method [ J. Chem. Phys. 2020, 152, 174113] within the Tensor Algebra for Many-body Methods (TAMM) infrastructure. TAMM is a massively parallel heterogeneous tensor library designed for utilizing forthcoming exascale computing resources. The two-body electron repulsion matrix elements are Cholesky-decomposed, and we imposed spin-explicit forms of the various operators when evaluating the tensor contractions. Unlike our previous real algebra Tensor Contraction Engine (TCE) implementation, the TAMM implementation supports fully complex algebra. The RT-EOM-CC singles (S) and doubles (D) time-dependent amplitudes are propagated using a first-order Adams-Moulton method. This new implementation shows excellent scalability tested up to 500 GPUs using the Zn-porphyrin molecule with 655 basis functions, with parallel efficiencies above 90% up to 400 GPUs. The TAMM RT-EOM-CCSD was used to study core photoemission spectra in the formaldehyde and ethyl trifluoroacetate (ESCA) molecules. Simulations of the latter involve as many as 71 occupied and 649 virtual orbitals. The relative quasiparticle ionization energies and overall spectral functions agree well with available experimental results.
Collapse
Affiliation(s)
- Himadri Pathak
- Advanced Computing, Mathematics, and Data Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Ajay Panyala
- Advanced Computing, Mathematics, and Data Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Bo Peng
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Nicholas P Bauman
- Physical Sciences Division, 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
| | - John J Rehr
- Department of Physics, University of Washington, Seattle, Washington 98195, United States
| | - Fernando D Vila
- Department of Physics, University of Washington, Seattle, Washington 98195, United States
| | - Karol Kowalski
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| |
Collapse
|
13
|
Ofstad BS, Kristiansen HE, Aurbakken E, Schøyen ØS, Kvaal S, Pedersen TB. Adiabatic extraction of nonlinear optical properties from real-time time-dependent electronic-structure theory. J Chem Phys 2023; 158:2882246. [PMID: 37093994 DOI: 10.1063/5.0145521] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2023] [Accepted: 03/27/2023] [Indexed: 04/26/2023] Open
Abstract
Real-time simulations of laser-driven electron dynamics contain information about molecular optical properties through all orders in response theory. These properties can be extracted by assuming convergence of the power series expansion of induced electric and magnetic multipole moments. However, the accuracy relative to analytical results from response theory quickly deteriorates for higher-order responses due to the presence of high-frequency oscillations in the induced multipole moment in the time domain. This problem has been ascribed to missing higher-order corrections. We here demonstrate that the deviations are caused by nonadiabatic effects arising from the finite-time ramping from zero to full strength of the external laser field. Three different approaches, two using a ramped wave and one using a pulsed wave, for extracting electrical properties from real-time time-dependent electronic-structure simulations are investigated. The standard linear ramp is compared to a quadratic ramp, which is found to yield highly accurate results for polarizabilities, and first and second hyperpolarizabilities, at roughly half the computational cost. Results for the third hyperpolarizability are presented along with a simple, computable measure of reliability.
Collapse
Affiliation(s)
- Benedicte Sverdrup Ofstad
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, Oslo, Norway
| | - Håkon Emil Kristiansen
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, Oslo, Norway
| | - Einar Aurbakken
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, Oslo, Norway
| | | | - Simen Kvaal
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, Oslo, Norway
| | - Thomas Bondo Pedersen
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, Oslo, Norway
| |
Collapse
|
14
|
Motta M, Jones GO, Rice JE, Gujarati TP, Sakuma R, Liepuoniute I, Garcia JM, Ohnishi YY. Quantum chemistry simulation of ground- and excited-state properties of the sulfonium cation on a superconducting quantum processor. Chem Sci 2023; 14:2915-2927. [PMID: 36937596 PMCID: PMC10016331 DOI: 10.1039/d2sc06019a] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2022] [Accepted: 02/14/2023] [Indexed: 02/17/2023] Open
Abstract
The computational description of correlated electronic structure, and particularly of excited states of many-electron systems, is an anticipated application for quantum devices. An important ramification is to determine the dominant molecular fragmentation pathways in photo-dissociation experiments of light-sensitive compounds, like sulfonium-based photo-acid generators used in photolithography. Here we simulate the static and dynamical electronic structure of the H3S+ molecule, taken as a minimal model of a triply-bonded sulfur cation, on a superconducting quantum processor of the IBM Falcon architecture. To this end, we generalize a qubit reduction technique termed entanglement forging or EF [A. Eddins et al., Phys. Rev. X Quantum, 2022, 3, 010309], currently restricted to the evaluation of ground-state energies, to the treatment of molecular properties. While in a conventional quantum simulation a qubit represents a spin-orbital, within EF a qubit represents a spatial orbital, reducing the number of required qubits by half. We combine the generalized EF with quantum subspace expansion [W. Colless et al., Phys. Rev. X, 2018, 8, 011021], a technique used to project the time-independent Schrodinger equation for ground- and excited-states in a subspace. To enable experimental demonstration of this algorithmic workflow, we deploy a sequence of error-mitigation techniques. We compute dipole structure factors and partial atomic charges along ground- and excited-state potential energy curves, revealing the occurrence of homo- and heterolytic fragmentation. This study is an important step towards the computational description of photo-dissociation on near-term quantum devices, as it can be generalized to other photodissociation processes and naturally extended in different ways to achieve more realistic simulations.
Collapse
Affiliation(s)
- Mario Motta
- IBM Quantum, IBM Research - Almaden 650 Harry Road San Jose 95120 CA USA
| | - Gavin O Jones
- IBM Quantum, IBM Research - Almaden 650 Harry Road San Jose 95120 CA USA
| | - Julia E Rice
- IBM Quantum, IBM Research - Almaden 650 Harry Road San Jose 95120 CA USA
| | - Tanvi P Gujarati
- IBM Quantum, IBM Research - Almaden 650 Harry Road San Jose 95120 CA USA
| | - Rei Sakuma
- Materials Informatics Initiative, RD Technology & Digital Transformation Center, JSR Corporation 3-103-9, Tonomachi, Kawasaki-ku Kawasaki 210-0821 Kanagawa Japan
| | - Ieva Liepuoniute
- IBM Quantum, IBM Research - Almaden 650 Harry Road San Jose 95120 CA USA
| | - Jeannette M Garcia
- IBM Quantum, IBM Research - Almaden 650 Harry Road San Jose 95120 CA USA
| | - Yu-Ya Ohnishi
- Materials Informatics Initiative, RD Technology & Digital Transformation Center, JSR Corporation 3-103-9, Tonomachi, Kawasaki-ku Kawasaki 210-0821 Kanagawa Japan
| |
Collapse
|
15
|
Kaye J, Barnett A, Greengard L, De Giovannini U, Rubio A. Eliminating Artificial Boundary Conditions in Time-Dependent Density Functional Theory Using Fourier Contour Deformation. J Chem Theory Comput 2023; 19:1409-1420. [PMID: 36786824 DOI: 10.1021/acs.jctc.2c01013] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/15/2023]
Abstract
We present an efficient method for propagating the time-dependent Kohn-Sham equations in free space, based on the recently introduced Fourier contour deformation (FCD) approach. For potentials which are constant outside a bounded domain, FCD yields a high-order accurate numerical solution of the time-dependent Schrödinger equation directly in free space, without the need for artificial boundary conditions. Of the many existing artificial boundary condition schemes, FCD is most similar to an exact nonlocal transparent boundary condition, but it works directly on Cartesian grids in any dimension, and runs on top of the fast Fourier transform rather than fast algorithms for the application of nonlocal history integral operators. We adapt FCD to time-dependent density functional theory (TDDFT), and describe a simple algorithm to smoothly and automatically truncate long-range Coulomb-like potentials to a time-dependent constant outside of a bounded domain of interest, so that FCD can be used. This approach eliminates errors originating from the use of artificial boundary conditions, leaving only the error of the potential truncation, which is controlled and can be systematically reduced. The method enables accurate simulations of ultrastrong nonlinear electronic processes in molecular complexes in which the interference between bound and continuum states is of paramount importance. We demonstrate results for many-electron TDDFT calculations of absorption and strong field photoelectron spectra for one and two-dimensional models, and observe a significant reduction in the size of the computational domain required to achieve high quality results, as compared with the popular method of complex absorbing potentials.
Collapse
Affiliation(s)
- Jason Kaye
- Center for Computational Mathematics, Flatiron Institute, 162 5th Avenue, New York, New York 10010, United States.,Center for Computational Quantum Physics, Flatiron Institute, 162 5th Avenue, New York, New York 10010, United States
| | - Alex Barnett
- Center for Computational Mathematics, Flatiron Institute, 162 5th Avenue, New York, New York 10010, United States
| | - Leslie Greengard
- Center for Computational Mathematics, Flatiron Institute, 162 5th Avenue, New York, New York 10010, United States.,Courant Institute of Mathematical Sciences, New York University, 251 Mercer Street, New York, New York 10012, United States
| | - Umberto De Giovannini
- Dipartimento di Fisica e Chimica - Emilio Segrè, Università degli Studi di Palermo, Via Archirafi 36, I-90123, Palermo, Italy.,Max Planck Institute for the Structure and Dynamics of Matter and Center for Free Electron Laser Science, 22761, Hamburg, Germany
| | - Angel Rubio
- Center for Computational Quantum Physics, Flatiron Institute, 162 5th Avenue, New York, New York 10010, United States.,Max Planck Institute for the Structure and Dynamics of Matter and Center for Free Electron Laser Science, 22761, Hamburg, Germany
| |
Collapse
|
16
|
Perrella F, Petrone A, Rega N. Understanding Charge Dynamics in Dense Electronic Manifolds in Complex Environments. J Chem Theory Comput 2023; 19:626-639. [PMID: 36602443 PMCID: PMC9878732 DOI: 10.1021/acs.jctc.2c00794] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2022] [Indexed: 01/06/2023]
Abstract
Photoinduced charge transfer (CT) excited states and their relaxation mechanisms can be highly interdependent on the environment effects and the consequent changes in the electronic density. Providing a molecular interpretation of the ultrafast (subpicosecond) interplay between initial photoexcited states in such dense electronic manifolds in condensed phase is crucial for improving and understanding such phenomena. Real-time time-dependent density functional theory is here the method of choice to observe the charge density, explicitly propagated in an ultrafast time domain, along with all time-dependent properties that can be easily extracted from it. A designed protocol of analysis for real-time electronic dynamics to be applied to time evolving electronic density related properties to characterize both in time and in space CT dynamics in complex systems is here introduced and validated, proposing easy to be read cross-correlation maps. As case studies to test such tools, we present the photoinduced charge-transfer electronic dynamics of 5-benzyluracil, a mimic of nucleic acid/protein interactions, and the metal-to-ligand charge-transfer electronic dynamics in water solution of [Ru(dcbpy)2(NCS)2]4-, dcbpy = (4,4'-dicarboxy-2,2'-bipyridine), or "N34-", a dye sensitizer for solar cells. Electrostatic and explicit ab initio treatment of solvent molecules have been compared in the latter case, revealing the importance of the accurate modeling of mutual solute-solvent polarization on CT kinetics. We observed that explicit quantum mechanical treatment of solvent slowed down the charge carriers mobilities with respect to the gas-phase. When all water molecules were modeled instead as simpler embedded point charges, the electronic dynamics appeared enhanced, with a reduced hole-electron distance and higher mean velocities due to the close fixed charges and an artificially increased polarization effect. Such analysis tools and the presented case studies can help to unveil the influence of the electronic manifold, as well as of the finite temperature-induced structural distortions and the environment on the ultrafast charge motions.
Collapse
Affiliation(s)
- Fulvio Perrella
- Department
of Chemical Sciences, University of Napoli
Federico II, Complesso Universitario di
M.S. Angelo, via Cintia 21, I-80126, Napoli, Italy
| | - Alessio Petrone
- Department
of Chemical Sciences, University of Napoli
Federico II, Complesso Universitario di
M.S. Angelo, via Cintia 21, I-80126, Napoli, Italy
- Scuola
Superiore Meridionale, Largo San Marcellino 10, I-80138, Napoli, Italy
- Istituto
Nazionale Di Fisica Nucleare, sezione di Napoli, Complesso Universitario di Monte S. Angelo ed. 6, via Cintia, 80126, Napoli, Italia
| | - Nadia Rega
- Department
of Chemical Sciences, University of Napoli
Federico II, Complesso Universitario di
M.S. Angelo, via Cintia 21, I-80126, Napoli, Italy
- Scuola
Superiore Meridionale, Largo San Marcellino 10, I-80138, Napoli, Italy
- Istituto
Nazionale Di Fisica Nucleare, sezione di Napoli, Complesso Universitario di Monte S. Angelo ed. 6, via Cintia, 80126, Napoli, Italia
- CRIB,
Centro Interdipartimentale di Ricerca sui Biomateriali, Piazzale Tecchio 80, I-80125, Napoli, Italy
| |
Collapse
|
17
|
Pathak H, Sato T, Ishikawa KL. Time-dependent optimized coupled-cluster method with doubles and perturbative triples for first principles simulation of multielectron dynamics. Front Chem 2022; 10:982120. [PMID: 36176891 PMCID: PMC9513851 DOI: 10.3389/fchem.2022.982120] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2022] [Accepted: 08/09/2022] [Indexed: 12/05/2022] Open
Abstract
We report the formulation of a new, cost-effective approximation method in the time-dependent optimized coupled-cluster (TD-OCC) framework [T. Sato et al., J. Chem. Phys. 148, 051101 (2018)] for first-principles simulations of multielectron dynamics in an intense laser field. The method, designated as TD-OCCD(T), is a time-dependent, orbital-optimized extension of the “gold-standard” CCSD(T) method in the ground-state electronic structure theory. The equations of motion for the orbital functions and the coupled-cluster amplitudes are derived based on the real-valued time-dependent variational principle using the fourth-order Lagrangian. The TD-OCCD(T) is size extensive and gauge invariant, and scales as O(N7) with respect to the number of active orbitals N. The pilot application of the TD-OCCD(T) method to the strong-field ionization and high-order harmonic generation from a Kr atom is reported in comparison with the results of the previously developed methods, such as the time-dependent complete-active-space self-consistent field (TD-CASSCF), TD-OCC with double and triple excitations (TD-OCCDT), TD-OCC with double excitations (TD-OCCD), and the time-dependent Hartree-Fock (TDHF) methods.
Collapse
Affiliation(s)
- Himadri Pathak
- Department of Nuclear Engineering and Management, School of Engineering, The University of Tokyo, Tokyo, Japan
| | - Takeshi Sato
- Department of Nuclear Engineering and Management, School of Engineering, The University of Tokyo, Tokyo, Japan
- Photon Science Center, School of Engineering, The University of Tokyo, Tokyo, Japan
- Research Institute for Photon Science and Laser Technology, The University of Tokyo, Tokyo, Japan
- *Correspondence: Takeshi Sato,
| | - Kenichi L. Ishikawa
- Department of Nuclear Engineering and Management, School of Engineering, The University of Tokyo, Tokyo, Japan
- Photon Science Center, School of Engineering, The University of Tokyo, Tokyo, Japan
- Research Institute for Photon Science and Laser Technology, The University of Tokyo, Tokyo, Japan
| |
Collapse
|
18
|
Park YC, Perera A, Bartlett RJ. Density functionals for core excitations. J Chem Phys 2022; 157:094107. [PMID: 36075732 DOI: 10.1063/5.0111095] [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
The core excitation energies and related principal ionization energies are obtained for selected molecules using several density functionals and compared with benchmark equation-of-motion coupled cluster (EOM-CC) results. Both time-dependent and time-independent formulations of excitation spectra in the time-dependent density functional theory and the EOM-CC are employed to obtain excited states that are not always easily accessible with the time-independent method. Among those functionals, we find that the QTP(00) functional, which is only parameterized to reproduce the five IPs of water, provides excellent core IPs and core excitation energies, consistently yielding better excitation and ionization energies. We show that orbital eigenvalues of KS density functional theory play an important role in determining the accuracy of the excitation and photoelectron spectra.
Collapse
Affiliation(s)
- Young Choon Park
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611-8435, USA
| | - Ajith Perera
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611-8435, USA
| | - Rodney J Bartlett
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611-8435, USA
| |
Collapse
|
19
|
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
|
20
|
Vila FD, Rehr JJ, Pathak H, Peng B, Panyala A, Mutlu E, Bauman NP, Kowalski K. Real-time equation-of-motion CC cumulant and CC Green's function simulations of photoemission spectra of water and water dimer. J Chem Phys 2022; 157:044101. [DOI: 10.1063/5.0099192] [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
Newly developed coupled-cluster (CC) methods enable simulations of ionization potentials and spectral functions of molecular systems in a wide range of energy scales ranging from core-binding to valence. This paper discusses results obtained with the real-time equation-of-motion CC cumulant approach (RT-EOM-CC), and CC Green's function (CCGF) approaches in applications to the water and water dimer molecules. We compare the ionization potentials obtained with these methods for the valence region with the results obtained with the CCSD(T) formulation as a difference of energies for N and N-1 electron systems. All methods show good agreement with each other. They also agree well with experiment, with errors usually below 0.1 eV for the ionization potentials.We also analyze unique features of the spectral functions, associated with the position of satellite peaks, obtained with the RT-EOM-CC and CCGF methods employing single and double excitations, as a function of the monomer OH bond length and the proton transfer coordinate in the dimer. Finally, we analyze the impact of the basis set effects on the quality of calculated ionization potentials and find that the basis set effects are less pronounced for the augmented-type sets.
Collapse
Affiliation(s)
| | - John J. Rehr
- Department of Physics, University of Washington College of Arts and Sciences, United States of America
| | - Himadri Pathak
- Pacific Northwest National Laboratory, Pacific Northwest National Laboratory, United States of America
| | - Bo Peng
- Pacific Northwest National Laboratory, United States of America
| | - Ajay Panyala
- Pacific Northwest National Laboratory, United States of America
| | - Erdal Mutlu
- Pacific Northwest National Laboratory, United States of America
| | | | - Karol Kowalski
- Environmental Molecular Sciences Laboratory, Pacific Northwest National Laboratory, United States of America
| |
Collapse
|
21
|
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
|
22
|
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]
|
23
|
Vila FD, Kowalski K, Peng B, Kas JJ, Rehr JJ. Real-Time Equation-of-Motion CCSD Cumulant Green's Function. J Chem Theory Comput 2022; 18:1799-1807. [PMID: 35157796 DOI: 10.1021/acs.jctc.1c01179] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Many-body excitations in X-ray photoemission spectra have been difficult to simulate from first principles. We have recently developed a cumulant-based one-electron Green's function method using the real-time coupled-cluster-singles equation-of-motion approach (RT-EOM-CCS) that provides a general framework for treating these problems. Here we extend this approach to include double excitations in the ground-state energy and in the coupled cluster amplitudes, which have been implemented using subroutines generated by the Tensor Contraction Engine (TCE). As in the case of the singles approximation, RT-EOM-CCSD yields a nonperturbative cumulant form of the Green's function in terms of the time-dependent cluster amplitudes, adding nonlinear corrections to the traditional cumulant forms. The extended approach is applied to the core-hole spectral function for small molecular systems. We find that, when core-optimized basis sets are used, the doubles contributions reduce the mean absolute errors in the core binding energies of the 10e systems from 0.8 to 0.3 eV. They also significantly improve the quasiparticle-satellite gap by reducing its overestimation from about 3-5 to about 0-1 eV in CH4, NH3, and H2O, and also improving the overall shape of the satellite features. Finally, we demonstrate the application of the new implementation to the larger, classical XPS ESCA series of molecules and show that the singles approximation can be paired with a modest basis set to study carbon speciation.
Collapse
Affiliation(s)
- F D Vila
- Department of Physics, University of Washington, Seattle, Washington 98195, United States
| | - K Kowalski
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - B Peng
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - J J Kas
- Department of Physics, University of Washington, Seattle, Washington 98195, United States
| | - J J Rehr
- Department of Physics, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
24
|
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
|
25
|
Park YC, Perera A, Bartlett RJ. Equation of motion coupled-cluster study of core excitation spectra II: Beyond the dipole approximation. J Chem Phys 2021; 155:094103. [PMID: 34496593 DOI: 10.1063/5.0059276] [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 present the time-independent (TI) and time-dependent (TD) equation of motion coupled-cluster (EOM-CC) oscillator strengths not limited to those obtained by the dipole approximation. For the conventional TI-EOM-CC, we implement all the terms in the multipole expansion through second order that contributes to the oscillator strength. These include contributions such as magnetic dipole, electric quadrupole, electric octupole, and magnetic quadrupole. In TD-EOM-CC, we only include the quadrupole moment contributions. This augments our previous work [Y. C. Park, A. Perera, and R. J. Bartlett, J. Chem. Phys. 151, 164117 (2019)]. The inclusion of the quadrupole contributions (and all the other contributions through second order in the case of TI-EOM-CCSD) enables us to obtain the intensities for the pre-edge transitions in the metal K-edge spectra, which are dipole inactive. The TI-EOM-CCSD and TD-EOM-CCSD spectra of Ti4+ atoms are used to showcase the implementation of the second-order oscillator strengths. The origin of 1s → e and 1s → t2 in core spectra from iron tetrachloride and titanium tetrachloride is discussed and compared with the experiment.
Collapse
Affiliation(s)
- Young Choon Park
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
| | - Ajith Perera
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
| | - Rodney J Bartlett
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
| |
Collapse
|
26
|
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
|
27
|
Pathak H, Sato T, Ishikawa KL. Time-dependent optimized coupled-cluster method for multielectron dynamics. IV. Approximate consideration of the triple excitation amplitudes. J Chem Phys 2021; 154:234104. [PMID: 34241273 DOI: 10.1063/5.0054743] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
We present a cost-effective treatment of the triple excitation amplitudes in the time-dependent optimized coupled-cluster (TD-OCC) framework called TD-OCCDT(4) for studying intense laser-driven multielectron dynamics. It considers triple excitation amplitudes correct up to the fourth-order in many-body perturbation theory and achieves a computational scaling of O(N7), with N being the number of active orbital functions. This method is applied to the electron dynamics in Ne and Ar atoms exposed to an intense near-infrared laser pulse with various intensities. We benchmark our results against the TD complete-active-space self-consistent field (TD-CASSCF), TD-OCC with double and triple excitations (TD-OCCDT), TD-OCC with double excitations (TD-OCCD), and TD Hartree-Fock (TDHF) methods to understand how this approximate scheme performs in describing nonperturbatively nonlinear phenomena, such as field-induced ionization and high-harmonic generation. We find that the TD-OCCDT(4) method performs equally well as the TD-OCCDT method, almost perfectly reproducing the results of the fully correlated TD-CASSCF with a more favorable computational scaling.
Collapse
Affiliation(s)
- 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
| | - Takeshi Sato
- 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
- Department of Nuclear Engineering and Management, School of Engineering, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8656, Japan
| |
Collapse
|
28
|
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
|
29
|
Moreno Carrascosa A, Yang M, Yong H, Ma L, Kirrander A, Weber PM, Lopata K. Mapping static core-holes and ring-currents with X-ray scattering. Faraday Discuss 2021; 228:60-81. [PMID: 33605956 DOI: 10.1039/d0fd00124d] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Measuring the attosecond movement of electrons in molecules is challenging due to the high temporal and spatial resolutions required. X-ray scattering-based methods are promising, but many questions remain concerning the sensitivity of the scattering signals to changes in density, as well as the means of reconstructing the dynamics from these signals. In this paper, we present simulations of stationary core-holes and electron dynamics following inner-shell ionization of the oxazole molecule. Using a combination of time-dependent density functional theory simulations along with X-ray scattering theory, we demonstrate that the sudden core-hole ionization produces a significant change in the X-ray scattering response and how the electron currents across the molecule should manifest as measurable modulations to the time dependent X-ray scattering signal. This suggests that X-ray scattering is a viable probe for measuring electronic processes at time scales faster than nuclear motion.
Collapse
Affiliation(s)
| | - Mengqi Yang
- Department of Chemistry, 232 Choppin Hall, Baton Rouge, Louisiana 70803, USA
| | - Haiwang Yong
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, USA
| | - Lingyu Ma
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, USA
| | - Adam Kirrander
- EaStCHEM, School of Chemistry, University of Edinburgh, David Brewster Road, EH9 3FJ Edinburgh, UK
| | - Peter M Weber
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, USA
| | - Kenneth Lopata
- Department of Chemistry, 232 Choppin Hall, Baton Rouge, Louisiana 70803, USA and Center for Computation and Technology, Louisiana State University, Baton Roug, Louisiana 70803, USA.
| |
Collapse
|
30
|
Church MS, Rubenstein BM. Real-time dynamics of strongly correlated fermions using auxiliary field quantum Monte Carlo. J Chem Phys 2021; 154:184103. [PMID: 34241020 DOI: 10.1063/5.0049116] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Spurred by recent technological advances, there is a growing demand for computational methods that can accurately predict the dynamics of correlated electrons. Such methods can provide much-needed theoretical insights into the electron dynamics probed via time-resolved spectroscopy experiments and observed in non-equilibrium ultracold atom experiments. In this article, we develop and benchmark a numerically exact Auxiliary Field Quantum Monte Carlo (AFQMC) method for modeling the dynamics of correlated electrons in real time. AFQMC has become a powerful method for predicting the ground state and finite temperature properties of strongly correlated systems mostly by employing constraints to control the sign problem. Our initial goal in this work is to determine how well AFQMC generalizes to real-time electron dynamics problems without constraints. By modeling the repulsive Hubbard model on different lattices and with differing initial electronic configurations, we show that real-time AFQMC is capable of accurately capturing long-lived electronic coherences beyond the reach of mean field techniques. While the times to which we can meaningfully model decrease with increasing correlation strength and system size as a result of the exponential growth of the dynamical phase problem, we show that our technique can model the short-time behavior of strongly correlated systems to very high accuracy. Crucially, we find that importance sampling, combined with a novel adaptive active space sampling technique, can substantially lengthen the times to which we can simulate. These results establish real-time AFQMC as a viable technique for modeling the dynamics of correlated electron systems and serve as a basis for future sampling advances that will further mitigate the dynamical phase problem.
Collapse
Affiliation(s)
- Matthew S Church
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, USA
| | - Brenda M Rubenstein
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, USA
| |
Collapse
|
31
|
Liu J, Cheng L. Relativistic coupled‐cluster and equation‐of‐motion coupled‐cluster methods. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2021. [DOI: 10.1002/wcms.1536] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- Junzi Liu
- Department of Chemistry The Johns Hopkins University Baltimore Maryland USA
| | - Lan Cheng
- Department of Chemistry The Johns Hopkins University Baltimore Maryland USA
| |
Collapse
|
32
|
Pedersen TB, Kristiansen HE, Bodenstein T, Kvaal S, Schøyen ØS. Interpretation of Coupled-Cluster Many-Electron Dynamics in Terms of Stationary States. J Chem Theory Comput 2021; 17:388-404. [PMID: 33337895 PMCID: PMC7808707 DOI: 10.1021/acs.jctc.0c00977] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2020] [Indexed: 01/06/2023]
Abstract
We demonstrate theoretically and numerically that laser-driven many-electron dynamics, as described by bivariational time-dependent coupled-cluster (CC) theory, may be analyzed in terms of stationary-state populations. Projectors heuristically defined from linear response theory and equation-of-motion CC theory are proposed for the calculation of stationary-state populations during interaction with laser pulses or other external forces, and conservation laws of the populations are discussed. Numerical tests of the proposed projectors, involving both linear and nonlinear optical processes for He and Be atoms and for LiH, CH+, and LiF molecules show that the laser-driven evolution of the stationary-state populations at the coupled-cluster singles-and-doubles (CCSD) level is very close to that obtained by full configuration interaction (FCI) theory, provided that all stationary states actively participating in the dynamics are sufficiently well approximated. When double-excited states are important for the dynamics, the quality of the CCSD results deteriorates. Observing that populations computed from the linear response projector may show spurious small-amplitude, high-frequency oscillations, the equation-of-motion projector emerges as the most promising approach to stationary-state populations.
Collapse
Affiliation(s)
- Thomas Bondo Pedersen
- Hylleraas
Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, N-0315 Oslo, Norway
| | - Håkon Emil Kristiansen
- Hylleraas
Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, N-0315 Oslo, Norway
| | - Tilmann Bodenstein
- Hylleraas
Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, N-0315 Oslo, Norway
| | - Simen Kvaal
- Hylleraas
Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, N-0315 Oslo, Norway
| | | |
Collapse
|
33
|
Vila FD, Rehr JJ, Kas JJ, Kowalski K, Peng B. Real-Time Coupled-Cluster Approach for the Cumulant Green's Function. J Chem Theory Comput 2020; 16:6983-6992. [PMID: 33108872 DOI: 10.1021/acs.jctc.0c00639] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
Green's function methods within many-body perturbation theory provide a general framework for treating electronic correlations in excited states and spectra. Here, we develop the cumulant form of the one-electron Green's function using a real-time coupled-cluster equation-of-motion approach, in an extension of our previous study (Rehr J.; et al. J. Chem. Phys. 2020, 152, 174113). The approach yields a nonperturbative expression for the cumulant in terms of the solution to a set of coupled first-order, nonlinear differential equations. The method thereby adds nonlinear corrections to traditional cumulant methods, which are linear in the self-energy. The approach is applied to the core-hole Green's function and is illustrated for a number of small molecular systems. For these systems, we find that the nonlinear contributions yield significant improvements, both for quasiparticle properties such as core-level binding energies and for inelastic losses that correspond to satellites observed in photoemission spectra.
Collapse
Affiliation(s)
- F D Vila
- Department of Physics, University of Washington, Seattle, Washington 98195, United States
| | - J J Rehr
- Department of Physics, University of Washington, Seattle, Washington 98195, United States
| | - J J Kas
- Department of Physics, University of Washington, Seattle, Washington 98195, United States
| | - K Kowalski
- William R. Wiley Environmental Molecular Sciences Laboratory, Battelle, Pacific Northwest National Laboratory, K8-91, P.O. Box 999, Richland, Washington 99352, United States
| | - B Peng
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| |
Collapse
|
34
|
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
|
35
|
Darapaneni P, Meyer AM, Sereda M, Bruner A, Dorman JA, Lopata K. Simulated field-modulated x-ray absorption in titania. J Chem Phys 2020; 153:054110. [PMID: 32770877 DOI: 10.1063/5.0009677] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this paper, we present a method to compute the x-ray absorption near-edge structure (XANES) spectra of solid-state transition metal oxides using real-time time-dependent density functional theory, including spin-orbit coupling effects. This was performed on bulk-mimicking anatase titania (TiO2) clusters, which allows for the use of hybrid functionals and atom-centered all electron basis sets. Furthermore, this method was employed to calculate the shifts in the XANES spectra of the Ti L-edge in the presence of applied electric fields to understand how external fields can modify the electronic structure, and how this can be probed using x-ray absorption spectroscopy. Specifically, the onset of t2g peaks in the Ti L-edge was observed to red shift and the eg peaks were observed to blue shift with increasing fields, attributed to changes in the hybridization of the conduction band (3d) orbitals.
Collapse
Affiliation(s)
- Pragathi Darapaneni
- Cain Department of Chemical Engineering, Louisiana State University, Baton Rouge, Louisiana 70803, USA
| | - Alexander M Meyer
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, USA
| | - Mykola Sereda
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, USA
| | - Adam Bruner
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, USA
| | - James A Dorman
- Cain Department of Chemical Engineering, Louisiana State University, Baton Rouge, Louisiana 70803, USA
| | - Kenneth Lopata
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, USA
| |
Collapse
|
36
|
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
|
37
|
Hansen MB, Madsen NK, Christiansen O. Extended vibrational coupled cluster: Stationary states and dynamics. J Chem Phys 2020; 153:044133. [DOI: 10.1063/5.0015413] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Mads Bøttger Hansen
- Department of Chemistry, Aarhus University, Langelandsgade 140, DK-8000 Aarhus C, Denmark
| | - Niels Kristian Madsen
- Department of Chemistry, Aarhus University, Langelandsgade 140, DK-8000 Aarhus C, Denmark
| | - Ove Christiansen
- Department of Chemistry, Aarhus University, Langelandsgade 140, DK-8000 Aarhus C, Denmark
| |
Collapse
|
38
|
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
|
39
|
Rehr JJ, Vila FD, Kas JJ, Hirshberg NY, Kowalski K, Peng B. Equation of motion coupled-cluster cumulant approach for intrinsic losses in x-ray spectra. J Chem Phys 2020; 152:174113. [PMID: 32384843 DOI: 10.1063/5.0004865] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a combined equation of motion coupled-cluster cumulant Green's function approach for calculating and understanding intrinsic inelastic losses in core level x-ray absorption spectra (XAS) and x-ray photoemission spectra. The method is based on a factorization of the transition amplitude in the time domain, which leads to a convolution of an effective one-body absorption spectrum and the core-hole spectral function. The spectral function characterizes intrinsic losses in terms of shake-up excitations and satellites using a cumulant representation of the core-hole Green's function that simplifies the interpretation. The one-body spectrum also includes orthogonality corrections that enhance the XAS at the edge.
Collapse
Affiliation(s)
- J J Rehr
- Department of Physics, University of Washington Seattle, Seattle, Washington 98195, USA
| | - F D Vila
- Department of Physics, University of Washington Seattle, Seattle, Washington 98195, USA
| | - J J Kas
- Department of Physics, University of Washington Seattle, Seattle, Washington 98195, USA
| | - N Y Hirshberg
- Department of Physics, University of Washington Seattle, Seattle, Washington 98195, USA
| | - K Kowalski
- Physical Sciences Division, Battelle, Pacific Northwest National Laboratory, K8-91, PO Box 999, Richland, Washington 99352, USA
| | - B Peng
- Physical Sciences Division, Battelle, Pacific Northwest National Laboratory, K8-91, PO Box 999, Richland, Washington 99352, USA
| |
Collapse
|
40
|
Pathak H, Sato T, Ishikawa KL. Time-dependent optimized coupled-cluster method for multielectron dynamics. II. A coupled electron-pair approximation. J Chem Phys 2020; 152:124115. [DOI: 10.1063/1.5143747] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- 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
| | - Takeshi Sato
- Department of Nuclear Engineering and Management, School of Engineering, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8656, Japan
- Photon Science Center, School of Engineering, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8656, Japan
- Research Institute for Photon Science and Laser Technology, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan
| | - Kenichi L. Ishikawa
- Department of Nuclear Engineering and Management, School of Engineering, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8656, Japan
- Photon Science Center, School of Engineering, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8656, Japan
- Research Institute for Photon Science and Laser Technology, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan
| |
Collapse
|
41
|
Kristiansen HE, Schøyen ØS, Kvaal S, Pedersen TB. Numerical stability of time-dependent coupled-cluster methods for many-electron dynamics in intense laser pulses. J Chem Phys 2020; 152:071102. [DOI: 10.1063/1.5142276] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Håkon Emil Kristiansen
- 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
| | - Thomas Bondo Pedersen
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway
| |
Collapse
|
42
|
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
|
43
|
Lingerfelt DB, Ganesh P, Jakowski J, Sumpter BG. Understanding Beam-Induced Electronic Excitations in Materials. J Chem Theory Comput 2020; 16:1200-1214. [DOI: 10.1021/acs.jctc.9b00792] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- David B. Lingerfelt
- Nanomaterials Theory Institute, Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| | - Panchapakesan Ganesh
- Nanomaterials Theory Institute, Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| | - Jacek Jakowski
- Nanomaterials Theory Institute, Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
- Computational Sciences and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| | - Bobby G. Sumpter
- Nanomaterials Theory Institute, Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
- Computational Sciences and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| |
Collapse
|
44
|
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
|
45
|
Schriber JB, Evangelista FA. Time dependent adaptive configuration interaction applied to attosecond charge migration. J Chem Phys 2019; 151:171102. [DOI: 10.1063/1.5126945] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jeffrey B. Schriber
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30318, USA
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Francesco A. Evangelista
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| |
Collapse
|
46
|
Park YC, Perera A, Bartlett RJ. Equation of motion coupled-cluster for core excitation spectra: Two complementary approaches. J Chem Phys 2019; 151:164117. [PMID: 31675901 DOI: 10.1063/1.5117841] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
This paper presents core excitation spectra from coupled-cluster (CC) theory obtained from both a time-independent and a new time-dependent formalism. The conventional time-independent CC formulation for excited states is the equation-of-motion (EOM-CC) method whose eigenvalues and eigenvectors describe the core excited states. An alternative computational procedure is offered by a time-dependent CC description. In that case, the dipole transition operator is expressed in the CC effective Hamiltonian form and propagated with respect to time. The absorption spectrum is obtained from the CC dipole autocorrelation function via a Fourier transformation. Comparisons are made among the time-dependent results obtained from second-order perturbation theory, to coupled cluster doubles and their linearized forms (CCD and LCCD), to CC singles and doubles (CCSD) and the linearized form (LCCSD). In the time-independent case, considerations of triples (EOM-CCSDT) and quadruples (EOM-CCSDTQ) are used to approach sub-electron volt accuracy. A particular target is the allyl radical, as an example of an open-shell molecule. As the results have to ultimately be the same, the two procedures offer a complementary approach toward analyzing experimental results.
Collapse
Affiliation(s)
- Young Choon Park
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
| | - Ajith Perera
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
| | - Rodney J Bartlett
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
| |
Collapse
|
47
|
Hansen MB, Madsen NK, Zoccante A, Christiansen O. Time-dependent vibrational coupled cluster theory: Theory and implementation at the two-mode coupling level. J Chem Phys 2019; 151:154116. [DOI: 10.1063/1.5117207] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Mads Bøttger Hansen
- Department of Chemistry, Aarhus University, Langelandsgade 140, DK–8000 Aarhus C, Denmark
| | - Niels Kristian Madsen
- Department of Chemistry, Aarhus University, Langelandsgade 140, DK–8000 Aarhus C, Denmark
| | - Alberto Zoccante
- Department of Chemistry, Aarhus University, Langelandsgade 140, DK–8000 Aarhus C, Denmark
| | - Ove Christiansen
- Department of Chemistry, Aarhus University, Langelandsgade 140, DK–8000 Aarhus C, Denmark
| |
Collapse
|
48
|
Koulias LN, Williams-Young DB, Nascimento DR, DePrince AE, Li X. Relativistic Real-Time Time-Dependent Equation-of-Motion Coupled-Cluster. J Chem Theory Comput 2019; 15:6617-6624. [DOI: 10.1021/acs.jctc.9b00729] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Lauren N. Koulias
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - David B. Williams-Young
- Computational Research Division, Lawrence Berkeley National Laboratory, 1 Cyclotron Road MS 50A-3111, Berkeley, California 94720, United States
| | - Daniel R. Nascimento
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306, United States
| | - A. Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306, United States
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
49
|
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
|
50
|
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
|