Short Iterative Lanczos Integration in Time-Dependent Equation-of-Motion Coupled-Cluster Theory

Approximate Exponential Integrators for Time-Dependent Equation-of-Motion Coupled Cluster Theory

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.

Cost-Efficient High-Resolution Linear Absorption Spectra through Extrapolating the Dipole Moment from Real-Time Time-Dependent Electronic-Structure Theory

We present a novel function fitting method for approximating the propagation of the time-dependent electric dipole moment from real-time electronic structure calculations. Real-time calculations of the electronic absorption spectrum require discrete Fourier transforms of the electric dipole moment. The spectral resolution is determined by the total propagation time, i.e., the trajectory length of the dipole moment, causing a high computational cost. Our developed method uses function fitting on shorter trajectories of the dipole moment, achieving arbitrary spectral resolution through extrapolation. Numerical testing shows that the fitting method can reproduce high-resolution spectra by using short dipole trajectories. The method converges with as little as 100 a.u. dipole trajectories for some systems, though the difficulty converging increases with the spectral density. We also introduce an error estimate of the fit, reliably assessing its convergence and hence the quality of the approximated spectrum.

Time-dependent equation-of-motion coupled-cluster simulations with a defective Hamiltonian

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.

Adiabatic extraction of nonlinear optical properties from real-time time-dependent electronic-structure theory

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.

Accelerating Real-Time Coupled Cluster Methods with Single-Precision Arithmetic and Adaptive Numerical Integration

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.

Real-time equation-of-motion CC cumulant and CC Green's function simulations of photoemission spectra of water and water dimer

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.

Linear and Nonlinear Optical Properties from TDOMP2 Theory

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.