• We present the time-dependent complete-active-space self-consistent-field (TD-CASSCF) method to simulate multielectron dynamics in ultrafast intense laser fields from the first principles. While based on multiconfiguration expansion, it divides the orbital space into frozen-core (tightly bound electrons with no response to the field), dynamical-core (electrons tightly bound but responding to the field), and active (fully correlated to describe highly excited and ejected electrons) orbital subspaces. The subspace decomposition can be done flexibly, conforming to phenomena under investigation and desired accuracy. The method is gauge invariant and size extensive. Infinite-range exterior complex scaling in addition to mask-function boundary is adopted as an efficient absorbing boundary. We show numerical examples and illustrate how to extract relevant physical quantities such as ionization yield, high-harmonic spectrum, and photoelectron spectrum from our full-dimensional implementation for atoms. The TD-CASSCF method will open a way to the ab initio simulation study of ultrafast intense laser science in realistic atoms and molecules.
  • We propose a gauge-invariant formulation of the channel orbital-based time-dependent configuration interaction singles (TDCIS) method [Phys. Rev. A 74, 043420 (2006)], one of the powerful ab initio methods to investigate electron dynamics in atoms and molecules subject to an external laser field. In the present formulation, we derive the equations of motion (EOMs) in the velocity gauge using gauge-transformed orbitals, not fixed orbitals, that are equivalent to the conventional EOMs in the length gauge using fixed orbitals. The new velocity-gauge EOMs avoid the use of the length-gauge dipole operator, which diverges at large distance, and allows to exploit computational advantages of the velocity-gauge treatment over the length-gauge one, e.g, a faster convergence in simulations with intense and long-wavelength lasers, and the feasibility of exterior complex scaling as an absorbing boundary. The reformulated TDCIS method is applied to an exactly solvable model of one-dimensional helium atom in an intense laser field to numerically demonstrate the gauge invariance. We also discuss the consistent method for evaluating the time derivative of an observable, relevant e.g, in simulating high-harmonic generation.
  • 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 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 to intense-laser driven correlated electron dynamics in Ar atom is reported.
  • We present a numerical implementation of the infinite-range exterior complex scaling (irECS) [Phys. Rev. A 81, 053845 (2010)] as an efficient absorbing boundary to the time-dependent complete-active-space self-consistent field (TD-CASSCF) method [Phys. Rev. A 94, 023405 (2016)] for multielectron atoms subject to an intense laser pulse. We introduce Gauss-Laguerre-Radau quadrature points to construct discrete variable representation basis functions in the last radial finite element extending to infinity. This implementation is applied to strong-field ionization and high-harmonic generation in He, Be, and Ne atoms. It efficiently prevents unphysical reflection of photoelectron wave packets at the simulation boundary, enabling accurate simulations with substantially reduced computational cost, even under significant (~ 50%) double ionization. For the case of a simulation of high-harmonic generation from Ne, for example, 80% cost reduction is achieved, compared to a mask-function absorption boundary.
  • We investigate the multielectron effects on high-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.
  • Single photon laser enabled Auger decay (spLEAD) has been redicted theoretically [Phys. Rev. Lett. 111, 083004 (2013)] and here we report its first experimental observation in neon. Using coherent, bichromatic free-electron laser pulses, we have detected the process and coherently controlled the angular distribution of the emitted electrons by varying the phase difference between the two laser fields. Since spLEAD is highly sensitive to electron correlation, this is a promising method for probing both correlation and ultrafast hole migration in more complex systems.
  • We present the fully general time-dependent multiconfiguration self-consistent-field method to describe the dynamics of a system consisting of arbitrary different kinds and numbers of interacting fermions and bosons. The total wave function is expressed as a superposition of different configurations constructed from time-dependent spin-orbitals prepared for each particle kind. We derive equations of motion followed by configuration-interaction (CI) coefficients and spin-orbitals for general, not restricted to full-CI, configuration spaces. The present method provides a flexible framework for the first-principles theoretical study of, e.g., correlated multielectron and multinucleus quantum dynamics in general molecules induced by intense laser fields and attosecond light pulses.
  • We theoretically study multielectron effects in high-harmonic generation (HHG), using all-electron first-principles simulations for a one-dimensional (1D) model atom. In addition to 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.
  • We theoretically study high-harmonic generation (HHG) from solids driven by intense laser pulses using a one-dimensional model periodic crystal. By numerically solving the time-dependent Schr\"{o}dinger equation directly on a real-space grid, we successfully reproduce experimentally observed unique features of solid-state HHG such as the linear cutoff-energy scaling and the sudden transition from a single- to multiple-plateau structure. Based on the simulation results, we propose a simple model that incorporates vector-potential-induced intraband displacement, interband tunneling, and recombination with the valence-band hole. One key parameter is the valley-to-peak amplitude of the pulse vector potential, which determines the crystal momentum displacement during the half cycle. When the maximum peak-to-valley amplitude $A_{\rm peak}$ reaches the half width $\frac{\pi}{a}$ of the Brillouin zone with $a$ being the lattice constant, the HHG spectrum exhibits a transition from a single- to multiple-plateau structure, and even further plateaus appear at $A_{\rm peak} = \frac{2\pi}{a}, \frac{3\pi}{a}, \cdots$. The multiple cutoff positions are given as functions of $A_{\rm peak}$ and the second maximum $A_{\rm peak}^{\prime}$, in terms of the energy difference between different bands. Using our recipe, one can draw electron trajectories in the momentum space, from which one can deduce, for example, the time-frequency structure of HHG without elaborate quantum-mechanical calculations. Finally, we reveal that the cutoff positions depend on not only the intensity and wavelength of the pulse, but also its duration, in marked contrast to the gas-phase case. Our model can be viewed as a solid-state and momentum-space counterpart of the familiar three-step model, highly successful for gas-phase HHG, and provide a unified basis to understand HHG from solid-state materials and gaseous media.
  • The accurate description of the non-linear response of many-electron systems to strong-laser fields remains a major challenge. Methods that bypass the unfavorable exponential scaling with particle number are required to address larger systems. In this paper we present a fully three-dimensional implementation of the time-dependent two-particle reduced density matrix (TD-2RDM) method for many-electron atoms. We benchmark this approach by a comparison with multi-configurational time-dependent Hartree-Fock (MCTDHF) results for the harmonic spectra of beryllium and neon. We show that the TD-2RDM is very well-suited to describe the non-linear atomic response and to reveal the influence of electron-correlation effects.
  • We present the 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 report a three-dimensional numerical implementation of multiconfiguration time-dependent Hartree-Fock (MCTDHF) based on a multi-resolution Cartesian grid, with no need to assume any symmetry of molecular structure. We successfully compute high-harmonic generation (HHG) of H2 and H2O. The present implementation will open a way to the first-principle theoretical study of intense-field and attosecond-pulse induced ultrafast phenomena in general molecules.
  • In parallel with the evolution of femtosecond and attosecond laser as well as free-electron laser technology, a variety of theoretical methods have been developed to describe the behavior of atoms, molecules, clusters, and solids under the action of those laser pulses. Here we review major ab initio wave-function-based numerical approaches to simulate multielectron dynamics in atoms and molecules driven by intense long-wavelength and/or ultrashort short-wavelength laser pulses. Direct solution of the time-dependent Schr\"odinger equation (TDSE), though its applicability is limited to He, ${\rm H}_2$, and Li, can provide an exact description and has been greatly contributing to the understanding of dynamical electron-electron correlation. Multiconfiguration self-consistent-field (MCSCF) approach offers a flexible framework from which a variety of methods can be derived to treat both atoms and molecules, with possibility to systematically control the accuracy. The equations of motion of configuration interaction coefficients and molecular orbitals for general MCSCF ansatz have recently been derived. Time-dependent extension of the $R$-matrix theory, originally develop for electron-atom collision, can realistically and accurately describe laser-driven complex multielectron atoms.
  • 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, {\bf 88}, 023402 (2013)], which divides the occupied orbitals into core and active orbitals, the TD-ORMAS method {\it further} subdivides the active orbitals into an arbitrary number of subgroups, and poses the {\it 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 both for 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.
  • Describing time-dependent many-body systems where correlation effects play an important role remains a major theoretical challenge. In this paper we develop a time-dependent many-body theory that is based on the two-particle reduced density matrix (2-RDM). We develop a closed equation of motion for the 2-RDM employing a novel reconstruction functional for the three-particle reduced density matrix (3-RDM) that preserves norm, energy, and spin symmetries during time propagation. We show that approximately enforcing $N$-representability during time evolution is essential for achieving stable solutions. As a prototypical test case which features long-range Coulomb interactions we employ the one-dimensional model for lithium hydride (LiH) in strong infrared laser fields. We probe both one-particle observables such as the time-dependent dipole moment and two-particle observables such as the pair density and mean electron-electron interaction energy. Our results are in very good agreement with numerically exact solutions for the $N$-electron wavefunction obtained from the multiconfigurational time-dependent Hartree-Fock method.
  • We theoretically investigate the mechanism of enhanced ionization in two-electron molecules by analyzing Bohmian trajectories for a one-dimensional H2 in an intense laser field. We identify both types of ionizing trajectories corresponding to the ejection from the up-field and down-field cores. The trajectories of the two electrons are correlated with each other in the former while correlation is negligible in the latter. The contributions from the two ionization types, though depending on laser intensity and internuclear distance, are comparable to each other.
  • We theoretically study the photoelectron angular distributions (PADs) from two-color two-photon near-threshold ionization of hydrogen and noble gas (He, Ne, and Ar) atoms by a combined action of femtosecond extreme ultraviolet (EUV) and near-infrared (IR) laser pulses. Using the second-order time-dependent perturbation theory, we clarify how the two-photon ionization process depends on EUV-IR pulse delay and how it is connected to the interplay between resonant and nonresonant ionization paths. Furthermore, by solving the time-dependent Schr\"odinger equation, we calculate the anisotropy parameters $\beta_2$ and $\beta_4$ as well as the amplitude ratio and relative phase between partial waves characterizing the PADs. We show that in general these parameters notably depend on the time delay between the EUV and IR pulses, except for He. This dependence is related to the varying relative role of resonant and nonresonant paths of photoionization. Our numerical results for H, He, Ne, and Ar show that the pulse-delay effect is more pronounced for $p$-shell ionization than for $s$-shell ionization.
  • 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-principle theoretical study of intense-field induced ultrafast phenomena in realistic atoms and molecules.
  • We theoretically study the pulse-width dependence of the photoelectron angular distribution (PAD) from the resonance-enhanced two-photon single ionization of He by femtosecond ($\lesssim 20$ fs) extreme-ultraviolet pulses, based on the time-dependent perturbation theory and simulations with the full time-dependent Schr\"odinger equation. In particular, we focus on the competition between resonant and nonresonant ionization paths, which leads to the relative phase $\delta$ between the $S$ and $D$ wave packets distinct from the corresponding scattering phase shift difference. When the spectrally broadened pulse is resonant with an excited level, the competition varies with pulse width, and, therefore, $\delta$ and the PAD also change with it. On the other hand, when the Rydberg manifold is excited, $\delta$ and the PAD do not much vary with the pulse width, except for the very short pulse regime.
  • Within a semiclassical description of above-threshold ionization (ATI) we identify the interplay between intracycle and intercycle interferences. The former is imprinted as a modulation envelope on the discrete multiphoton peaks formed by the latter. This allows to unravel the complex interference pattern observed for the full solution of the time-dependent Schr\"odinger equation (TDSE) in terms of diffraction at a grating in the time domain. These modulations can be clearly seen in the dependence of the ATI spectra on the laser wavelength. Shifts in energy modulation result from the effect of the long Coulomb tail of the atomic potential.
  • We theoretically study the scaling with the driving wavelength of the high-order harmonic generation (HHG) under the simultaneous irradiation of an extreme ultraviolet (XUV) pulse. Surprisingly, when the cutoff energy and ionization yield are fixed, the harmonic yield is nearly independent of the driving wavelength. We identify its origin as the combination of the initial spatial width of the states excited by the XUV pulse, making the wavepacket spreading less prominent, and the shallowing of the ionization potential, which suggests complex nature of the wavelength dependence of HHG.