We investigate the numerical stability of time-dependent coupled-cluster theory for many-electron dynamics in intense laser pulses, comparing two coupled-cluster formulations with full configuration interaction theory. Our numerical experiments show that orbital-adaptive time-dependent coupled-cluster doubles (OAT-DCCD) theory offers significantly improved stability compared with the conventional Hartree-Fock-based time-dependent coupled-cluster singles-and-doubles (TDCCSD) formulation. The improved stability stems from greatly reduced oscillations in the doubles amplitudes, which, in turn, can be traced to the dynamic biorthonormal reference determinants of OATDCCD theory. As long as these are good approximations to the Brueckner determinant, OATDCCD theory is numerically stable. We propose the reference weight as a diagnostic quantity to identify situations where the TDCCSD and OATDCCD theories become unstable.This document gives additional results to the main article. Included are demonstrations of non-collapsing OATDCCD and TDCCSD simulations of Be in the cc-pVDZ basis, discussion of controllable error in the twoparticle case for TDFCI and OATDCCD, Be in the cc-pVTZ basis, and systematic breaking of the TDCCSD and OATDCCD methods for LiH in the aug-cc-pVDZ basis by increasing the field strength.
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.
Recent years have witnessed an increasing interest in time‐dependent coupled‐cluster (TDCC) theory for simulating laser‐driven electronic dynamics in atoms and molecules, and for simulating molecular vibrational dynamics. Starting from the time‐dependent bivariational principle, we review different flavors of single‐reference TDCC theory with either orthonormal static, orthonormal time‐dependent, or biorthonormal time‐dependent spin orbitals. The time‐dependent extension of equation‐of‐motion coupled‐cluster theory is also discussed, along with the applications of TDCC methods to the calculation of linear absorption spectra, linear and low‐order nonlinear response functions, highly nonlinear high harmonic generation spectra and ionization dynamics. In addition, the role of TDCC theory in finite‐temperature many‐body quantum mechanics is briefly described along with a few other application areas.This article is categorized under: Electronic Structure Theory > Ab Initio Electronic Structure Methods Theoretical and Physical Chemistry > Spectroscopy Software > Simulation Methods
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.
We have developed a McMurchie-Davidson-like recursion formula for efficient evaluation of the Coulomb attraction and interaction matrix elements between two-dimensional primitive Cartesian Gaussian type orbitals. We also present recurrence schemes for combined position and differential operator integrals, and three-center Gaussian integrals. The Cartesian Gaussian orbitals are isotropic in the exponent, but with arbitrary centers and angular momentum.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.