Time-dependent optimized coupled-cluster method with doubles and perturbative triples for first principles simulation of multielectron dynamics

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(N⁷) 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.

1 Introduction

Recent years witnessed unprecedented progress in laser technologies, which made it possible to observe the motions of electrons at the attosecond time scale (Itatani et al. (2004); Corkum and Krausz (2007); Krausz and Ivanov (2009); Baker et al. (2006)). On the other hand, various theoretical and numerical methods have been developed for interpreting, understanding, and predicting the experiments.

The multi-configuration time-dependent Hartree-Fock (MCTDHF) method (Caillat et al. (2005); Kato and Kono (2004); Nest et al. (2005); Haxton et al. (2011); Hochstuhl and Bonitz (2011)), and the time-dependent complete-active-space self-consistent-field (TD-CASSCF) method (Sato and Ishikawa (2013); Sato et al. (2016); Sato et al. (2018a)) are the most rigorous approaches to solve time-dependent Schrödinger equation (TDSE) of many-electron systems, where the wavefunction is given by the full configuration interaction (FCI) expansion,

$\mathrm{\Psi }\left(t\right)=\sum _{\mathit{I}}{C}_I\left(t\right){\mathrm{\Phi }}_{\mathit{I}}\left(t\right),(1)$

with both CI coefficients {C_I(t)} and orbital functions {ψ_p(t)} constituting Slater determinants {Φ_I(t)} are propagated in time according to the time-dependent variational principle (TDVP). The TD-CASSCF method broadens the applicability of the MCTDHF method by flexibly classifying the orbital subspace into frozen-core, dynamical-core, and active. Unfortunately, the factorial computational scaling impedes large-scale applications. There are reports of various affordable size-inextensive methods (Miyagi and Madsen (2013, 2014); Haxton and McCurdy (2015); Sato and Ishikawa (2015)) developed by limiting the CI expansion of the wavefunction. Alternatively, the size-extensive coupled-cluster method, which relies on an exponential wavefunction, is a superior choice to address these problems with a polynomial cost-scaling (Kümmel (2003); Shavitt and Bartlett (2009)). We have developed an explicitly time-dependent coupled-cluster method considering optimized orthonormal orbitals within the flexibly chosen active space, called the time-dependent optimized coupled-cluster (TD-OCC) method, (Sato et al. (2018b)) including double (TD-OCCD) and double and triple excitation amplitudes (TD-OCCDT). Our method is a time-dependent formulation of the stationary optimized coupled-cluster method (Scuseria and Schaefer (1987); Sherrill et al. (1998); Krylov et al. (1998)). Kvaal (Kvaal (2012)) also developed an orbital adaptive time-dependent coupled-cluster (OATDCC) method using biorthogonal orbitals. We take note of a few reports on the time-dependent coupled-cluster methods (Huber and Klamroth (2011); Pigg et al. (2012); Nascimento and DePrince (2016)), using time-independent orbitals, and their interpretation (Pedersen and Kvaal (2019); Pedersen et al. (2021)), including the very initial attempts (Schonhammer (1978); Hoodbhoy and Negele (1978, 1979)).

The TD-OCCDT scales as O(N⁸) (N= the number of active orbitals), not ideally suited for applications to larger chemical systems. Therefore, we have developed a few lower cost methods in the TD-OCC framework (Pathak et al. (2020b,c,a, 2021)). We find triple excitations are necessary, including perfect optimization of the orbitals. Therefore, we are interested in developing affordable TD-OCC methods retaining a part of the triples. The most popular coupled-cluster method that treats the triple excitation amplitudes approximately is called CCSD(T) (Raghavachari et al. (1989); Watts et al. (1993)). Bozkaya et al, (Bozkaya and Schaefer (2012)) included various symmetric and asymmetric triple excitation corrections to their optimized double (OD) method.

In this communication, we report the formulation and implementation of the CCSD(T) method in the time-dependent optimized coupled-cluster framework, TD-OCCD(T). Following our previous works (Sato et al. (2018b); Pathak et al. (2020b,c, 2021)), we exclude single excitation amplitudes but optimize the orbitals according to time-dependent variational principle (TDVP). As the first application of this method, we study electron dynamics in Kr using intense near-infrared laser fields.

2 Methods

The second quantization representation of the Hamiltonian, including the laser field, is as follows,

$\widehat{H}=h_{\nu }^{\mu }\left(t\right){\widehat{c}}_{\mu }^{†}{\widehat{c}}_{\nu }+\frac{1}{2}{u}_{\nu \lambda }^{\mu \gamma }{\widehat{c}}_{\mu }^{†}{\widehat{c}}_{\gamma }^{†}{\widehat{c}}_{\lambda }{\widehat{c}}_{\nu }(2)$

where ${\widehat{c}}_{\mu }^{†}$ $({\widehat{c}}_{\mu })$ represents a creation (annihilation) operator in a complete, orthonormal set of 2n_bas time-dependent spin-orbitals {ψ_μ(t)}. n_bas is the number of basis functions used for expanding the spatial part of ψ_μ, which, in the present real-space implementation, corresponds to the number of grid points, and

$h_{\nu }^{\mu }\left(t\right)=\int d{x}_1{\psi }_{\mu }^{\ast }\left(x_1\right)\left[h_0+V_{\mathit{ext}}\right]{\psi }_{\nu }\left(x_1\right),(3)$

$u_{\nu \lambda }^{\mu \gamma }=\iint d{x}_{1}d{x}_2\frac{{\psi }_{\mu }^{\ast }\left(x_1\right){\psi }_{\gamma }^{\ast }\left(x_2\right){\psi }_{\nu }\left(x_1\right){\psi }_{\lambda }\left(x_2\right)}{|{\mathit{r}}_1-{\mathit{r}}_2|},(4)$

where x_i = (r_i, σ_i) represents a composite spatial-spin coordinate. h₀ is the field free one-electronic Hamiltonian and V_ext = A(t)p_z in the velocity gauge, A(t) = −∫^tE(t′)dt′ is the vector potential, with E(t) being the laser electric field linearly polarized along the z axis.

The complete set of 2n_bas spin-orbitals (labeled with μ, ν, γ, λ) is divided into n_occ occupied (o, p, q, r, s) and 2n_bas − n_occ virtual spin-orbitals. The coupled-cluster (or CI) wavefunction is constructed only with occupied spin-orbitals, which are time-dependent in general, and virtual spin-orbitals form the orthogonal complement of the occupied spin-orbital space. The occupied spin-orbitals are classified into n_core core spin-orbitals, which are occupied in the reference Φ and kept uncorrelated, and N = n_occ − n_core active spin-orbitals (t, u, v, w) among which the active electrons are correlated. The active spin-orbitals are further split into those in the hole space (i, j, k, l) and the particle space (a, b, c, d), which are defined as those occupied and unoccupied, respectively, in the reference Φ. The core spin-orbitals can further be split into frozen-core space (i′′, j′′), fixed in time and the dynamical-core space (i′, j′), propagated in time (Sato and Ishikawa (2013)) (See. Figure 1 in Sato et al. (2018b) for a pictorial illustration).

FIGURE 1

FIGURE 1. Time evolution of dipole moment of Kr irradiated by a laser pulse with a wavelength of 800 nm and a peak intensity of 2 × 10¹⁴ W/cm² calculated with TDHF, TD-OCCD, TD-OCCD(T), and TD-CASSCF methods.

The real action formulation of the TDVP with orthonormal orbitals is our guiding principle, (Sato et al. (2018b))

$S=\mathrm{R}\mathrm{e}{\int }_{t_0}^{t_1}Ldt=\frac{1}{2}{\int }_{t_0}^{t_1}\left(L+L^{\ast }\right)dt,(5)$

$L=⟨\mathrm{\Phi }|\left(1+\widehat{\mathrm{\Lambda }}\right)e^{-\widehat{T}}\left(\widehat{H}-i\frac{\partial }{\partial t}\right)e^{\widehat{T}}|\mathrm{\Phi }⟩,(6)$

$\widehat{T}={\widehat{T}}_2+{\widehat{T}}_3\cdots ={\tau }_{ij}^{ab}{\widehat{E}}_{ij}^{ab}+{\tau }_{\mathit{ijk}}^{\mathit{abc}}{\widehat{E}}_{\mathit{ijk}}^{\mathit{abc}}\cdots ,(7)$

$\widehat{\mathrm{\Lambda }}={\widehat{\mathrm{\Lambda }}}_2+{\widehat{\mathrm{\Lambda }}}_3\cdots ={\lambda }_{ab}^{ij}{\widehat{E}}_{ab}^{ij}+{\lambda }_{\mathit{abc}}^{\mathit{ijk}}{\widehat{E}}_{\mathit{abc}}^{\mathit{ijk}}\cdots ,(8)$

where ${\tau }_{ij\cdots }^{ab\cdots }$ $({\lambda }_{ab\cdots }^{ij\cdots })$ are (de-)excitation amplitudes, and ${\widehat{E}}_{ij\cdots }^{ab\cdots }={\widehat{c}}_a^{†}{\widehat{b}}_b^{†}\cdots {\widehat{c}}_j{\widehat{c}}_i$. The stationary conditions, δS = 0, with respect to the variation of the parameters of the wavefunction ($\delta {\tau }_{ij\cdots }^{ab\cdots }$, $\delta {\lambda }_{ab\cdots }^{ij\cdots }$, and δψ_μ) gives us the corresponding equations of motions (EOMs), δψ_μ is orthonormality-conserving orbital variation.

For deriving the TD-OCCD(T) method, we first construct a fourth-order Lagrangian defined in Pathak et al. (2021). We make a further approximation to the Lagrangian and write separating it into two parts,

$L_{\text{CCD}\left(\mathrm{T}\right)}^{\left(4\right)}=L_0+⟨\mathrm{\Phi }|\left(1+{\widehat{\mathrm{\Lambda }}}_2\right){\left[\left(\bar{f}+\widehat{v}\right)e^{{\widehat{T}}_2}\right]}_c|\mathrm{\Phi }⟩-i{\lambda }_{ab}^{ij}{\dot{\tau }}_{ij}^{ab}(\mathrm{9a})$

$+⟨\mathrm{\Phi }|{\widehat{\mathrm{\Lambda }}}_2{\left[\left(\bar{f}+\widehat{v}\right){\widehat{T}}_3\right]}_c|\mathrm{\Phi }⟩+⟨\mathrm{\Phi }|{\widehat{\mathrm{\Lambda }}}_3{\left(\bar{f}{\widehat{T}}_3\right)}_c|\mathrm{\Phi }⟩+⟨\mathrm{\Phi }|{\widehat{\mathrm{\Lambda }}}_3{\left(\widehat{v}{\widehat{T}}_2\right)}_c|\mathrm{\Phi }⟩-i{\lambda }_{\mathit{abc}}^{\mathit{ijk}}{\dot{\tau }}_{\mathit{ijk}}^{\mathit{abc}},(\mathrm{9b})$

where $\bar{f}=\widehat{f}-i\widehat{X}$, $\widehat{f}=(h_q^p+v_{qj}^{pj})\{{\widehat{E}}_q^p\}$, $\widehat{v}=v_{qs}^{pr}\{{\widehat{E}}_{qs}^{pr}\}/4$, and $v_{qs}^{pr}=u_{qs}^{pr}-u_{sq}^{pr}$, $\widehat{X}=X_{\nu }^{\mu }{\widehat{E}}_{\nu }^{\mu }$, and $X_{\nu }^{\mu }=⟨{\psi }_{\mu }|{\dot{\psi }}_{\nu }⟩$ is anti-Hermitian. The double amplitudes are obtained by making $L_{\text{CCD}\left(\text{T}\right)}^{(4)}$ of Eq. 9a stationary with respect to $\delta S/\delta {\lambda }_{ab}^{ij}(t)=0$, $\delta S/\delta {\tau }_{ij}^{ab}(t)=0$, the triples by making Eq. 9b stationary with respect to $\delta S/\delta {\lambda }_{\mathit{abc}}^{\mathit{ijk}}(t)=0$, and $\delta S/\delta {\tau }_{\mathit{ijk}}^{\mathit{abc}}(t)=0$,

$\begin{array}{ccc}\hfill i{\dot{\tau }}_{ij}^{ab}& =\hfill & v_{ij}^{ab}-p\left(ij\right){\bar{f}}_j^k{\tau }_{ik}^{ab}+p\left(ab\right){\bar{f}}_c^a{\tau }_{ij}^{cb}\hfill \\ \hfill & +\hfill & \frac{1}{2}{v}_{cd}^{ab}{\tau }_{ij}^{cd}+\frac{1}{2}{v}_{ij}^{kl}{\tau }_{kl}^{ab}+p\left(ij\right)p\left(ab\right)v_{ic}^{ak}{\tau }_{kj}^{cb}\hfill \\ \hfill & -\hfill & \frac{1}{2}p\left(ij\right){\tau }_{ik}^{ab}{\tau }_{jl}^{cd}{v}_{cd}^{kl}+\frac{1}{2}p\left(ab\right){\tau }_{ij}^{bc}{\tau }_{kl}^{ad}{v}_{cd}^{kl}\hfill \\ \hfill & +\hfill & \frac{1}{4}{\tau }_{kl}^{ab}{\tau }_{ij}^{cd}{v}_{cd}^{kl}+\frac{1}{2}p\left(ij\right)p\left(ab\right){\tau }_{il}^{bc}{\tau }_{jk}^{ad}{v}_{cd}^{kl}\hfill \end{array}(10)$

$\begin{array}{ccc}\hfill -i{\dot{\lambda }}_{ab}^{ij}& =\hfill & v_{ab}^{ij}-p\left(ij\right){\bar{f}}_k^i{\lambda }_{ab}^{kj}+p\left(ab\right){\bar{f}}_a^c{\lambda }_{cb}^{ij}\hfill \\ \hfill & +\hfill & \frac{1}{2}{v}_{ab}^{cd}{\lambda }_{cd}^{ij}+\frac{1}{2}{v}_{kl}^{ij}{\lambda }_{ab}^{kl}+p\left(ij\right)p\left(ab\right)v_{kb}^{cj}{\lambda }_{ac}^{ik}\hfill \\ \hfill & -\hfill & \frac{1}{2}p\left(ij\right){\lambda }_{cd}^{ik}{\tau }_{kl}^{cd}{v}_{ab}^{jl}+\frac{1}{2}p\left(ab\right){\lambda }_{bc}^{kl}{\tau }_{kl}^{cd}{v}_{ad}^{ij}\hfill \\ \hfill & +\hfill & \frac{1}{4}{\lambda }_{ab}^{kl}{\tau }_{kl}^{cd}{v}_{cd}^{ij}+\frac{1}{2}p\left(ij\right)p\left(ab\right){\lambda }_{ac}^{jk}{\tau }_{kl}^{cd}{v}_{bd}^{il}\hfill \\ \hfill & -\hfill & \frac{1}{2}p\left(ij\right){\lambda }_{ab}^{ik}{\tau }_{kl}^{cd}{v}_{cd}^{jl}\hfill \\ \hfill & +\hfill & \frac{1}{2}p\left(ab\right){\lambda }_{bc}^{ij}{\tau }_{kl}^{cd}{v}_{ad}^{kl}+\frac{1}{4}{\lambda }_{cd}^{ij}{\tau }_{kl}^{cd}{v}_{ab}^{kl}\hfill \end{array}(11)$

$\begin{array}{ccc}\hfill i{\dot{\tau }}_{\mathit{ijk}}^{\mathit{abc}}& =\hfill & p\left(k/ij\right)p\left(a/bc\right)v_{dk}^{bc}{t}_{ij}^{ad}-p\left(i/jk\right)p\left(c/ab\right)v_{jk}^{lc}{t}_{il}^{ab}\hfill \\ \hfill & -\hfill & p\left(k/ij\right){\bar{f}}_k^l{\tau }_{\mathit{ijl}}^{\mathit{abc}}+p\left(c/ab\right){\bar{f}}_d^c{\tau }_{\mathit{ijk}}^{\mathit{abd}},\hfill \end{array}(12)$

$\begin{array}{ccc}\hfill -i{\dot{\lambda }}_{\mathit{abc}}^{\mathit{ijk}}& =\hfill & p\left(k/ij\right)p\left(a/bc\right)v_{bc}^{dk}{\lambda }_{ad}^{ij}-p\left(c/ab\right)p\left(i/jk\right)v_{lc}^{jk}{\lambda }_{ab}^{ij}\hfill \\ \hfill & +\hfill & p\left(c/ab\right){\bar{f}}_c^d{\lambda }_{\mathit{abd}}^{\mathit{ijk}}-p\left(k/ij\right){\bar{f}}_l^k{\lambda }_{\mathit{abc}}^{\mathit{ijl}}\hfill \\ \hfill & +\hfill & p\left(i/jk\right)p\left(a/bc\right){\bar{f}}_a^i{\lambda }_{bc}^{jk},\hfill \end{array}(13)$

where p(μν) and p(μ|νγ) are the permutation operators; p(μν)A_μν = A_μν − A_νμ, and p(μ/νγ) = 1 − p(μν) − p(μγ).

The EOM for the orbitals can be written down in the following form Sato et al. (2016),

$i|\dot{{\psi }_p}\rangle =\left(\widehat{1}-\widehat{P}\right)\widehat{F}|{\psi }_p\rangle +i|{\psi }_q\rangle X_p^q,(14)$

where $\widehat{1}={\sum }_{\mu }|{\psi }_{\mu }\rangle \langle {\psi }_{\mu }|$ is the identity operator within the space spanned by the given basis, $\widehat{P}={\sum }_q|{\psi }_q\rangle \langle {\psi }_q|$ is the projector onto the occupied spin-orbital space, and

$\widehat{F}|{\psi }_p\rangle =\widehat{h}|{\psi }_p\rangle +{\widehat{W}}_s^r|{\psi }_q\rangle P_{or}^{qs}{\left(D^{-1}\right)}_p^o,(15)$

where D and P are Hermitialized one- (1RDM) and two- (2RDM) particle reduced density matrices defined in Sato et al. (2018b), and $W_s^r$ is the mean-field operator (Sato and Ishikawa (2013)). The matrix element $X_p^q$ includes orbital rotations among various subspaces. Non-redundant orbital rotations are determined by $i\left({\delta }_b^a{D}_i^j-D_b^a{\delta }_i^j\right)X_j^b=F_p^a{D}_i^p-D_p^a{F}_p^{i\ast }-\frac{i}{8}{\dot{\tau }}_{\mathit{ijk}}^{\mathit{abc}}{\lambda }_{bc}^{jk}-\frac{i}{8}{\tau }_{\mathit{ijk}}^{\mathit{abc}}{\dot{\lambda }}_{bc}^{jk}$. Redundant orbital rotations $\{X_{j^{\prime }}^{i^{\prime }}\}$, $\{X_j^i\}$, and $\{X_b^a\}$ can be arbitrary antiHermitian matrix elements. The general expressions for the RDMs are the same as in the TD-OCCDT(4) method (Pathak et al. (2021)).

3 Numerical results and discussion

Our numerical implementation has an interface with the Gaussian09 program (Frisch et al. (2009)) for checking ground state energy with the standard Gaussian basis results. We study BH molecule with double-ζ plus polarization (DZP). We have reported ground state energy computed by propagating in the imaginary time for OCCD and OCCD(T) methods in Table 1 and compared those with the optimized double and asymmetric triple excitation corrections for the orbital-optimized doubles method of Bozkaya et al., Bozkaya and Schaefer (2012). We also compare our OCCD ground state energy result with Krylov et al.,Krylov et al. (1998) within the chosen active space of six electrons correlated among the six optimized active orbitals. We obtained a perfect agreement for all available values.

TABLE 1

TABLE 1. Comparison of the ground state energy of BH (r_e=2.4 bohr) molecule in DZP basis^a.

We have used a spherical-finite-element-discrete-variable representation (FEDVR) basis for representing orbital functions, Sato et al. (2016); Orimo et al. (2018) ${\chi }_{\mathit{klm}}(r,\theta ,\psi )=\frac{1}{r}{f}_k(r)Y_{lm}(\theta ,\varphi )$ where Y_lm and f_k(r) are spherical harmonics and the normalized radial-FEDVR basis function, respectively. The expansion of the spherical harmonics continued up to the maximum angular momentum L_max, and the radial FEDVR basis supports the range of radial coordinate 0 ≤ r ≤ R_max, with cos^1/4 mask function used as an absorbing boundary for avoiding unphysical reflection from the wall of the simulation box. We have used l_max = 72, and the FEDVR basis supporting the radial coordinate 0 < r < 300 using 78 finite elements each containing 25 DVR functions. The absorbing boundary is switched on at r = 180 in all our simulations. The Fourth-order exponential Runge-Kutta method (Hochbruck and Ostermann (2010)) is used to propagate the EOMs with 20000 time steps for each optical cycle. We run the simulations for a further 6,000 time steps after the end of the pulse. In all correlation calculations, eight electrons of 4s4p orbitals are considered as active and correlated among thirteen active orbitals. We report simulation results computed using a three-cycle laser pulse with a central wavelength of 800 nm having intensity 2 × 10¹⁴ W/cm² and a period of T = 2π/ω₀ ∼ 2.67 fs.

We report the time evolution of dipole moment of Kr in Figure 1 and in Figure 2 single electron ionization probability. Time-dependent dipole moment is evaluated as a trace $⟨{\psi }_p|\widehat{z}|{\psi }_q⟩D_p^q$ using 1RDMs. For the single electron ionization probability, we computed the probability of finding an electron outside a sphere of a radius of 20 a.u. using RDMs defined in Refs. 19; 20; 37. We compare the results of TD-CASSCF, TD-OCCD(T), TD-OCCD, and TDHF methods.

FIGURE 2

FIGURE 2. Time evolution of single ionization probability of Kr irradiated by a laser pulse with a wavelength of 800 nm and a peak intensity of 2 × 10¹⁴ W/cm² calculated with TDHF, TD-OCCD, TD-OCCD(T), and TD-CASSCF methods.

We observe a substantial underestimation (both in Figure 1, and Figure 2) by the TDHF method due to the lack of correlation treatment. All correlation methods perform according to their ability to treat electron correlation. We also computed results using the TD-OCCDT method but not reported here since those results are not identifiable from the TD-CASSCF results within the graphical resolution.

Next, we report high-harmonic generation in Figure 3. It is calculated by squaring the modulus I(ω) = |a(ω)|² of the Fourier transform of the expectation value of the dipole acceleration with a modified Ehrenfest expression (Sato et al. (2016)). In panel (c) of Figure 3, we plot the absolute relative deviation (δ(ω), of the spectral amplitude a(ω) from the TD-CASSCF value for each method. All methods qualitatively predict similar HHG spectra with TDHF underestimates the spectral intensity. The relative deviation of results from TD-CASSCF ones follows the general trend TDHF$>$TD-OCCD$>$TD-OCCD(T)$>$TD-OCCDT, the same as what we observe for the time-dependent dipole moment and single ionization probability. We also simulated results with lower and higher intensity. However, the trend remains the same.

FIGURE 3

FIGURE 3. The HHG spectra (A,B) and the relative deviation (C) of the spectral amplitude from the TD-CASSCF spectrum from Kr irradiated by a laser pulse with a wavelength of 800 nm and a peak intensity of 2 × 10¹⁴ W/cm² with various methods.

Finally, we make a tally of computational costs for all the methods considered in this article. All simulations performed using an Intel(R) Xeon(R) Gold 6,230 central processing unit (CPU) with 40 processors with a clock speed of 2.10 GHz, and report total simulations time in Table 2. Further, we report a reduction in the computational cost for various TD-OCC methods relative to the TD-CASSCF. We see a massive 63% cost reduction for the TD-OCCD(T) method, which is larger than for the TD-OCCDT method (58%), and a minimal increase from the TD-OCCD method.

TABLE 2

TABLE 2. Comparison of the total simulation time^a (in min) spent for TD-CASSCF, TD-OCCDT, TDCCD(T), and TD-OCCD methods.

4 Concluding remarks

We have reported the formulation and implementation of the TD-OCCD(T) method. As the first application, we employed this method to study laser-driven dynamics in Kr exposed to an intense near-infrared laser pulse. We observe a 63% cost reduction in comparison to the TD-CASSCF method without losing much accuracy. Therefore, we conclude that TD-OCCD(T) method will certainly be beneficial in exploring highly accurate ab initio simulations of electron dynamics in larger chemical systems.

Data availability statement

The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.

Author contributions

HP and TS formulated the method. HP numerically implemented the method and performed simulations. All the authors analyzed the results and contributed to the submitted version of the manuscript.

Funding

This research was supported in part by a Grant-in-Aid for Scientific Research (Grants No. JP18H03891 and No. JP19H00869) from the Ministry of Education, Culture, Sports, Science and Technology (MEXT) of Japan. This research was also partially supported by JST COI (Grant No. JPMJCE1313), JST CREST (Grant No. JPMJCR15N1), and by MEXT Quantum Leap Flagship Program (MEXT Q-LEAP) Grant Number JPMXS0118067246.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Baker, S., Robinson, J. S., Haworth, C., Teng, H., Smith, R., Chirilă, C., et al. (2006). Probing proton dynamics in molecules on an attosecond time scale. Science 312, 424–427. doi:10.1126/science.1123904

PubMed Abstract | CrossRef Full Text | Google Scholar

Bozkaya, U., and Schaefer, H. F. (2012). Symmetric and asymmetric triple excitation corrections for the orbital-optimized coupled-cluster doubles method: Improving upon ccsd (t) and ccsd (t) λ: Preliminary application. J. Chem. Phys. 136, 204114. doi:10.1063/1.4720382

PubMed Abstract | CrossRef Full Text | Google Scholar

Caillat, J., Zanghellini, J., Kitzler, M., Koch, O., Kreuzer, W., and Scrinzi, A. (2005). Correlated multielectron systems in strong laser fields: A multiconfiguration time-dependent Hartree-Fock approach. Phys. Rev. A . Coll. Park. 71, 012712. doi:10.1103/physreva.71.012712

CrossRef Full Text | Google Scholar

Corkum, P. á., and Krausz, F. (2007). Attosecond science. Nat. Phys. 3, 381–387. doi:10.1038/nphys620

CrossRef Full Text | Google Scholar

[Dataset] Frisch, M., Trucks, G., Schlegel, H., Scuseria, G., Robb, M., Cheeseman, J., et al. (2009). Gaussian 09, revision d. 01.

Google Scholar

Haxton, D. J., Lawler, K. V., and McCurdy, C. W. (2011). Multiconfiguration time-dependent Hartree-Fock treatment of electronic and nuclear dynamics in diatomic molecules. Phys. Rev. A . Coll. Park. 83, 063416. doi:10.1103/physreva.83.063416

CrossRef Full Text | Google Scholar

Haxton, D. J., and McCurdy, C. W. (2015). Two methods for restricted configuration spaces within the multiconfiguration time-dependent Hartree-Fock method. Phys. Rev. A . Coll. Park. 91, 012509. doi:10.1103/physreva.91.012509

CrossRef Full Text | Google Scholar

Hochbruck, M., and Ostermann, A. (2010). Exponential integrators. Acta Numer. 19, 209–286. doi:10.1017/s0962492910000048

CrossRef Full Text | Google Scholar

Hochstuhl, D., and Bonitz, M. (2011). Two-photon ionization of helium studied with the multiconfigurational time-dependent Hartree–Fock method. J. Chem. Phys. 134, 084106. doi:10.1063/1.3553176

PubMed Abstract | CrossRef Full Text | Google Scholar

Hoodbhoy, P., and Negele, J. W. (1978). Time-dependent coupled-cluster approximation to nuclear dynamics. i. application to a solvable model. Phys. Rev. C 18, 2380–2394. doi:10.1103/physrevc.18.2380

CrossRef Full Text | Google Scholar

Hoodbhoy, P., and Negele, J. W. (1979). Time-dependent coupled-cluster approximation to nuclear dynamics. ii. general formulation. Phys. Rev. C 19, 1971–1982. doi:10.1103/physrevc.19.1971

CrossRef Full Text | Google Scholar

Huber, C., and Klamroth, T. (2011). Explicitly time-dependent coupled cluster singles doubles calculations of laser-driven many-electron dynamics. J. Chem. Phys. 134, 054113. doi:10.1063/1.3530807

PubMed Abstract | CrossRef Full Text | Google Scholar

Itatani, J., Levesque, J., Zeidler, D., Niikura, H., Pépin, H., Kieffer, J. C., et al. (2004). Tomographic imaging of molecular orbitals. Nature 432, 867–871. doi:10.1038/nature03183

PubMed Abstract | CrossRef Full Text | Google Scholar

Kato, T., and Kono, H. (2004). Time-dependent multiconfiguration theory for electronic dynamics of molecules in an intense laser field. Chem. Phys. Lett. 392, 533–540. doi:10.1016/j.cplett.2004.05.106

CrossRef Full Text | Google Scholar

Krausz, F., and Ivanov, M. (2009). Attosecond physics. Rev. Mod. Phys. 81, 163–234. doi:10.1103/revmodphys.81.163

CrossRef Full Text | Google Scholar

Krylov, A. I., Sherrill, C. D., Byrd, E. F., and Head-Gordon, M. (1998). Size-consistent wave functions for nondynamical correlation energy: The valence active space optimized orbital coupled–cluster doubles model. J. Chem. Phys. 109, 10669–10678. doi:10.1063/1.477764

CrossRef Full Text | Google Scholar

Kümmel, H. G. (2003). A biography of the coupled cluster method. Int. J. Mod. Phys. B 17, 5311–5325. doi:10.1142/s0217979203020442

CrossRef Full Text | Google Scholar

Kvaal, S. (2012). Ab initio quantum dynamics using coupled-cluster. J. Chem. Phys. 136, 194109. doi:10.1063/1.4718427

PubMed Abstract | CrossRef Full Text | Google Scholar

Miyagi, H., and Madsen, L. B. (2013). Time-dependent restricted-active-space self-consistent-field theory for laser-driven many-electron dynamics. Phys. Rev. A . Coll. Park. 87, 062511. doi:10.1103/physreva.87.062511

CrossRef Full Text | Google Scholar

Miyagi, H., and Madsen, L. B. (2014). Time-dependent restricted-active-space self-consistent-field theory for laser-driven many-electron dynamics. ii. extended formulation and numerical analysis. Phys. Rev. A . Coll. Park. 89, 063416. doi:10.1103/physreva.89.063416

CrossRef Full Text | Google Scholar

Nascimento, D. R., and DePrince, A. E. (2016). Linear absorption spectra from explicitly time-dependent equation-of-motion coupled-cluster theory. J. Chem. Theory Comput. 12, 5834–5840. doi:10.1021/acs.jctc.6b00796

PubMed Abstract | CrossRef Full Text | Google Scholar

Nest, M., Klamroth, T., and Saalfrank, P. (2005). The multiconfiguration time-dependent Hartree–Fock method for quantum chemical calculations. J. Chem. Phys. 122, 124102. doi:10.1063/1.1862243

PubMed Abstract | CrossRef Full Text | Google Scholar

Orimo, Y., Sato, T., Scrinzi, A., and Ishikawa, K. L. (2018). Implementation of the infinite-range exterior complex scaling to the time-dependent complete-active-space self-consistent-field method. Phys. Rev. A. Coll. Park. 97, 023423. doi:10.1103/physreva.97.023423

CrossRef Full Text | Google Scholar

Pathak, H., Sato, T., and Ishikawa, K. L. (2020a). Study of laser-driven multielectron dynamics of ne atom using time-dependent optimised second-order many-body perturbation theory. Mol. Phys. 118, e1813910. doi:10.1080/00268976.2020.1813910

CrossRef Full Text | Google Scholar

Pathak, H., Sato, T., and Ishikawa, K. L. (2021). Time-dependent optimized coupled-cluster method for multielectron dynamics iv: Approximate consideration of the triple excitation amplitudes. J. Chem. Phys. 154, 234104. doi:10.1063/5.0054743

PubMed Abstract | CrossRef Full Text | Google Scholar

Pathak, H., Sato, T., and Ishikawa, K. L. (2020b). Time-dependent optimized coupled-cluster method for multielectron dynamics. ii. a coupled electron-pair approximation. J. Chem. Phys. 152, 124115. doi:10.1063/1.5143747

PubMed Abstract | CrossRef Full Text | Google Scholar

Pathak, H., Sato, T., and Ishikawa, K. L. (2020c). Time-dependent optimized coupled-cluster method for multielectron dynamics. iii. a second-order many-body perturbation approximation. J. Chem. Phys. 153, 034110. doi:10.1063/5.0008789

PubMed Abstract | CrossRef Full Text | Google Scholar

Pedersen, T. B., Kristiansen, H. E., Bodenstein, T., Kvaal, S., and Schøyen, Ø. S. (2021). Interpretation of coupled-cluster many-electron dynamics in terms of stationary states. J. Chem. Theory Comput. 17, 388–404. doi:10.1021/acs.jctc.0c00977

PubMed Abstract | CrossRef Full Text | Google Scholar

Pedersen, T. B., and Kvaal, S. (2019). Symplectic integration and physical interpretation of time-dependent coupled-cluster theory. J. Chem. Phys. 150, 144106. doi:10.1063/1.5085390

PubMed Abstract | CrossRef Full Text | Google Scholar

Pigg, D. A., Hagen, G., Nam, H., and Papenbrock, T. (2012). Time-dependent coupled-cluster method for atomic nuclei. Phys. Rev. C 86, 014308. doi:10.1103/physrevc.86.014308

CrossRef Full Text | Google Scholar

Raghavachari, K., Trucks, G. W., Pople, J. A., and Head-Gordon, M. (1989). A fifth-order perturbation comparison of electron correlation theories. Chem. Phys. Lett. 157, 479–483. doi:10.1016/s0009-2614(89)87395-6

CrossRef Full Text | Google Scholar

Sato, T., Ishikawa, K. L., Březinová, I., Lackner, F., Nagele, S., and Burgdörfer, J. (2016). Time-dependent complete-active-space self-consistent-field method for atoms: Application to high-order harmonic generation. Phys. Rev. A . Coll. Park. 94, 023405. doi:10.1103/physreva.94.023405

CrossRef Full Text | Google Scholar

Sato, T., and Ishikawa, K. L. (2013). Time-dependent complete-active-space self-consistent-field method for multielectron dynamics in intense laser fields. Phys. Rev. A . Coll. Park. 88, 023402. doi:10.1103/physreva.88.023402

CrossRef Full Text | Google Scholar

Sato, T., and Ishikawa, K. L. (2015). Time-dependent multiconfiguration self-consistent-field method based on the occupation-restricted multiple-active-space model for multielectron dynamics in intense laser fields. Phys. Rev. A . Coll. Park. 91, 023417. doi:10.1103/physreva.91.023417

CrossRef Full Text | Google Scholar

Sato, T., Orimo, Y., Teramura, T., Tugs, O., and Ishikawa, K. L. (2018a). “Time-dependent complete-active-space self-consistent-field method for ultrafast intense laser science,” in Progress in ultrafast intense laser science XIV (Berlin, Germany: Springer), 143–171.

CrossRef Full Text | Google Scholar

Sato, T., Pathak, H., Orimo, Y., and Ishikawa, K. L. (2018b). Communication: Time-dependent optimized coupled-cluster method for multielectron dynamics. J. Chem. Phys. 148, 051101. doi:10.1063/1.5020633

PubMed Abstract | CrossRef Full Text | Google Scholar

Schonhammer, K., and Gunnarsson, O. (1978). Time-dependent approach to the calculation of spectral functions. Phys. Rev. B 18, 6606. doi:10.1103/physrevb.18.6606

CrossRef Full Text | Google Scholar

Scuseria, G. E., and Schaefer, H. F. (1987). The optimization of molecular orbitals for coupled cluster wavefunctions. Chem. Phys. Lett. 142, 354–358. doi:10.1016/0009-2614(87)85122-9

CrossRef Full Text | Google Scholar

Shavitt, I., and Bartlett, R. J. (2009). Many-body methods in chemistry and physics: MBPT and coupled-cluster theory. Cambridge: Cambridge University Press.

Google Scholar

Sherrill, C. D., Krylov, A. I., Byrd, E. F., and Head-Gordon, M. (1998). Energies and analytic gradients for a coupled-cluster doubles model using variational brueckner orbitals: Application to symmetry breaking in o 4+. J. Chem. Phys. 109, 4171–4181. doi:10.1063/1.477023

CrossRef Full Text | Google Scholar

Watts, J. D., Gauss, J., and Bartlett, R. J. (1993). Coupled-cluster methods with noniterative triple excitations for restricted open-shell Hartree–Fock and other general single determinant reference functions. energies and analytical gradients. J. Chem. Phys. 98, 8718–8733. doi:10.1063/1.464480

CrossRef Full Text | Google Scholar