No Arabic 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 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 time-dependent complete-active-space self-consistent field (TD-CASSCF), time-dependent optimized coupled-cluster with double and triple excitations (TD-OCCDT), time-dependent optimized coupled-cluster with double excitations (TD-OCCD), and the time-dependent 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 fully-correlated TD-CASSCF with a more favorable computational scaling.
We calculate the high-harmonic generation (HHG) spectra, strong-field ionization, and time-dependent dipole-moment of Ne using explicitly time-dependent optimized second-order many-body perturbation method (TD-OMP2) where both orbitals and amplitudes are time-dependent. We consider near-infrared (800 nm) and mid-infrared (1200 nm) laser pulses with very high intensities ($5times10^{14}$, $8times10^{14}$ , and $1times10^{15}$ W/cm$^2$), required for strong-field experiments with the high-ionization potential (21.6 eV) atom. We compare the result of the TD-OMP2 method with the time-dependent complete-active-space self-consistent field method and the time-dependent Hartree-Fock method. Further, we report the implementation of the TD-CC2 method within the chosen active space, which is also a second-order approximation to the TD-CCSD method, and present results of time-dependent dipole-moment and HHG spectra with an intensity of $5times10^{13}$ W/cm$^2$ at a wavelength of 800 nm. It is found that the TD-CC2 method is not stable in the case with a higher laser intensity, and it does not provide a gauge-invariant description of the physical properties, which makes TD-OMP2 a superior choice to reach out to larger chemical systems, especially for the study of strong-field dynamics. The obtained results indicate that the TD-OMP2 method shows moderate performance, overestimating the response of Ne, while TDHF underestimates it. Nevertheless, it is remarkable that stable computation of such highly nonlinear nonperturbative phenomena is possible within the framework of time-dependent perturbation method, by virtue of the nonperturbative inclusion of the laser-electron interaction and time-dependent optimization of orbitals.
Molecular iodine was photoexcited by a strong 800 nm laser, driving several channels of multiphoton excitation. The motion following photoexcitation was probed using time-resolved X-ray scattering, which produces a scattering map $S(Q,tau)$. Temporal Fourier transform methods were employed to obtain a frequency-resolved X-ray scattering signal $tilde{S}(Q,omega)$. Taken together, $S(Q,tau)$ and $tilde{S}(Q,omega)$ separate different modes of motion, so that mode-specific nuclear oscillatory positions, oscillation amplitudes, directions of motions, and times may be measured accurately. Molecular dissociations likewise have a distinct signature, which may be used to identify both velocities and dissociation time shifts, and also can reveal laser-induced couplings among the molecular potentials.
Greens function methods within many-body perturbation theory provide a general framework for treating electronic correlations in excited states. Here we investigate the cumulant form of the one-electron Greens function based on the coupled-cluster equation of motion approach in an extension of our previous study. The approach yields a non-perturbative expression for the cumulant in terms of the solution to a set of coupled first order, non-linear differential equations. The method thereby adds non-linear corrections to traditional cumulant methods linear in the self energy. The approach is applied to the core-hole Greens function and illustrated for a number of small molecular systems. For these systems we find that the non-linear contributions lead to significant improvements both for quasiparticle properties such as core-level binding energies, as well as the satellites corresponding to inelastic losses observed in photoemission spectra.
Since in coupled-cluster (CC) theory ground-state and excitation energies are eigenvalues of a non-Hermitian matrix, these energies can in principle take on complex values. In this paper we discuss the appearance of complex energy values in CC calculations from a mathematical perspective. We analyze the behaviour of the eigenvalues of Hermitian matrices that are perturbed (in a non-Hermitian manner) by a real parameter. Based on these results we show that for CC calculations with real-valued Hamiltonian matrices the ground-state energy generally takes a real value. Furthermore, we show that in the case of real-valued Hamiltonian matrices complex excitation energies only occur in the context of conical intersections. In such a case, unphysical consequences are encountered such as a wrong dimension of the intersection seam, large numerical deviations from full configuration-interaction (FCI) results, and the square-root-like behaviour of the potential surfaces near the conical intersection. In the case of CC calculations with complex-valued Hamiltonian matrix elements, it turns out that complex energy values are to be expected for ground and excited states when no symmetry is present. We confirm the occurrence of complex energies by sample calculations using a six-state model and by CC calculations for the H2O molecule in a strong magnetic field. We furthermore show that symmetry can prevent the occurrence of complex energy values. Lastly, we demonstrate that in most cases the real part of the complex energy values provides a very good approximation to the FCI energy.
The validity of optimized dynamical decoupling (DD) is extended to analytically time dependent Hamiltonians. As long as an expansion in time is possible the time dependence of the initial Hamiltonian does not affect the efficiency of optimized dynamical decoupling (UDD, Uhrig DD). This extension provides the analytic basis for (i) applying UDD to effective Hamiltonians in time dependent reference frames, for instance in the interaction picture of fast modes and for (ii) its application in hierarchical DD schemes with $pi$ pulses about two perpendicular axes in spin space. to suppress general decoherence, i.e., longitudinal relaxation and dephasing.