Time-dependent coupled-cluster method with time-varying orbital functions, called time-dependent optimized coupled-cluster (TD-OCC) method, is formulated for multielectron dynamics in an intense laser field. We have successfully derived the equations of motion for CC amplitudes and orthonormal orbital functions based on the real action functional, and implemented the method including double excitations (TD-OCCD) and double and triple excitations (TD-OCCDT) within the optimized active orbitals. The present method is size extensive and gauge invariant, a polynomial cost-scaling alternative to the time-dependent multiconfiguration self-consistent-field method. The first application of the TD-OCC method of intense-laser driven correlated electron dynamics in Ar atom is reported.

Photoinduced phase transitions have been intensively studied owing to the potential capability to control a material of interest in the ultrafast manner, which can induce exotic phases unable to be attained at equilibrium. However, the key mechanisms are still under debate, and it is currently a central issue as to how the couplings between the electron, lattice, and spin degrees of freedom are evolving during photoinduced phase transitions. Here, we use a recently developed analysis method, which we call frequency-domain angle-resolved photoemission spectroscopy (FDARPES), and reveal mode- and band-selective electron-phonon couplings during the photoinduced insulator-to-metal transition for ${\mathrm{Ta}}_2\mathrm{Ni}{\mathrm{Se}}_5$. We find that the lattice modulation corresponding to the 2 THz phonon mode, where the Ta lattice is sheared along the $a$ axis, is the most relevant for the emergence of photoinduced semimetallic state. Furthermore, we find that the semimetallic and semiconducting bands coexist in the transient state, and demonstrate that FDARPES spectra can selectively detect the phonon-specific couplings to the two coexistent band structures during the photoinduced phase transition by resolving them in the frequency domain.

We theoretically study high-harmonic generation (HHG) from transition-metal elements Mn and ${\mathrm{Mn}}^{+}$ using full-dimensional, all-electron, first-principles simulations. The HHG spectra calculated with the time-dependent complete-active-space self-consistent-field (TD-CASSCF) and occupation-restricted multiple-active-space (TD-ORMAS) methods exhibit a prominent peak at $\sim 50$ eV, successfully reproducing resonant enhancement observed in previous experiments [ Opt. Express 20, 25239 (2012)]. Artificially freezing $3p$ orbitals in simulations results in its disappearance, which shows the essential role played by $3p$ electrons in the resonant harmonics (RH). Further transition-resolved analysis unambiguously identifies constructively interfering $3p\text{−}3d$ ($m=0,\pm 1$) giant resonance transitions as the origin of the RH, as also implied by its position in the spectra. Time-frequency analysis indicates that the recolliding electron combines with the parent ion to form the upper state of the transitions. In addition, this study shows that the TD-CASSCF and TD-ORMAS methods can be applied to open-shell atoms with many unpaired inner electrons.

Quantum mechanically, photoionization can be fully described by the complex photoionization amplitudes that describe the transition between the ground state and the continuum state. Knowledge of the value of the phase of these amplitudes has been a central interest in photoionization studies and newly developing attosecond science, since the phase can reveal important information about phenomena such as electron correlation. We present a new attosecond-precision interferometric method of angle-resolved measurement for the phase of the photoionization amplitudes, using two phase-locked extreme ultraviolet pulses of frequency $\omega$ and $2\omega$, from a free-electron laser. Phase differences $\mathrm{\Delta }\tilde{\eta }$ between one- and two-photon ionization channels, averaged over multiple wave packets, are extracted for neon $2p$ electrons as a function of the emission angle at photoelectron energies 7.9, 10.2, and 16.6 eV. $\mathrm{\Delta }\tilde{\eta }$ is nearly constant for emission parallel to the electric vector but increases at 10.2 eV for emission perpendicular to the electric vector. We model our observations with both perturbation and ab initio theory and find excellent agreement. In the existing method for attosecond measurement, reconstruction of attosecond beating by interference of two-photon transitions (RABBITT), a phase difference between two-photon pathways involving absorption and emission of an infrared photon is extracted. Our method can be used for extraction of a phase difference between single-photon and two-photon pathways and provides a new tool for attosecond science, which is complementary to RABBITT.

We investigate the multielectron effects on high-order harmonic generation from solid-state materials using the time-dependent Hartree-Fock theory. We find qualitative change in harmonic spectra, in particular, multiple-plateau formation at significantly lower laser intensities than within the independent-electron approximation. We reveal its origin in terms of interband polarization, i.e., electron-hole polarization, enabling interband excitation at remote crystal momenta via Coulomb potential.

We theoretically study multielectron effects in high-harmonic generation (HHG), using all-electron first-principles simulations for a one-dimensional model atom. In addition to the usual plateau and cutoff (from a cation in the present case, since the neutral is immediately ionized), we find a prominent resonance peak far above the plateau and a second plateau extended beyond the first cutoff. These features originate from the dication response enhanced by orders of magnitude due to the action of the Coulomb force from the rescattering electron, and, hence, are a clear manifestation of electron correlation. Although the present simulations are done in 1D, the physical mechanism underlying the dramatic enhancement is expected to hold also for three-dimensional real systems. This will provide new possibilities to explore dynamical electron correlation in intense laser fields using HHG, which is usually considered to be of single-electron nature in most cases.