The time-dependent complete-active-space self-consistent-field (TD-CASSCF) method for the description of multielectron dynamics in intense laser fields is presented, and a comprehensive description of the method is given. It introduces the concept of frozen-core (to model tightly bound electrons with no response to the field), dynamical-core (to model electrons tightly bound but responding to the field), and active (fully correlated to describe ionizing electrons) orbital subspaces, allowing compact yet accurate representation of ionization dynamics in many-electron systems. The classification into the subspaces can be done flexibly, according to simulated physical situations and desired accuracy, and the multiconfiguration time-dependent Hartree-Fock (MCTDHF) approach is included as a special case. To assess its performance, we apply the TD-CASSCF method to the ionization dynamics of one-dimensional lithium hydride (LiH) and LiH dimer models, and confirm that the present method closely reproduces rigorous MCTDHF results if active orbital space is chosen large enough to include appreciably ionizing electrons. The TD-CASSCF method will open a way to the first-principles theoretical study of intense-field-induced ultrafast phenomena in realistic atoms and molecules.

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.

The time-dependent multiconfiguration self-consistent-field method based on the occupation-restricted multiple-active-space model is proposed (TD-ORMAS) for multielectron dynamics in intense laser fields. Extending the previously proposed time-dependent complete-active-space self-consistent-field method [TD-CASSCF; Phys. Rev. A 88, 023402 (2013)], which divides the occupied orbitals into core and active orbitals, the TD-ORMAS method further subdivides the active orbitals into an arbitrary number of subgroups and poses the occupation restriction by giving the minimum and maximum number of electrons distributed in each subgroup. This enables highly flexible construction of the configuration-interaction (CI) space, allowing a large-active-space simulation of dynamics, e.g., the core excitation or ionization. The equations of motion for both CI coefficients and spatial orbitals are derived based on the time-dependent variational principle, and an efficient algorithm is proposed to solve for the orbital time derivatives. In-depth descriptions of the computational implementation are given in a readily programmable manner. The numerical application to the one-dimensional lithium hydride cluster models demonstrates that the high flexibility of the TD-ORMAS framework allows for the cost-effective simulations of multielectron dynamics by exploiting systematic series of approximations to the TD-CASSCF method.

We present a numerical implementation of the time-dependent complete-active-space self-consistent-field (TD-CASSCF) method [ Phys. Rev. A 88, 023402 (2013)] for atoms driven by a strong linearly polarized laser pulse. The present implementation treats the problem in its full dimensionality and introduces a gauge-invariant frozen-core approximation, an efficient evaluation of the Coulomb mean field scaling linearly with the number of basis functions, and a split-operator method specifically designed for stable propagation of stiff spatial derivative operators. We apply this method to high-harmonic generation in helium, beryllium, and neon and explore the role of electron correlations.

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.