
We present the timedependent completeactivespace selfconsistentfield
(TDCASSCF) 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 frozencore (tightly bound
electrons with no response to the field), dynamicalcore (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.
Infiniterange exterior complex scaling in addition to maskfunction boundary
is adopted as an efficient absorbing boundary. We show numerical examples and
illustrate how to extract relevant physical quantities such as ionization
yield, highharmonic spectrum, and photoelectron spectrum from our
fulldimensional implementation for atoms. The TDCASSCF method will open a way
to the ab initio simulation study of ultrafast intense laser science in
realistic atoms and molecules.

We propose a gaugeinvariant formulation of the channel orbitalbased
timedependent 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 gaugetransformed orbitals, not fixed orbitals, that are
equivalent to the conventional EOMs in the length gauge using fixed orbitals.
The new velocitygauge EOMs avoid the use of the lengthgauge dipole operator,
which diverges at large distance, and allows to exploit computational
advantages of the velocitygauge treatment over the lengthgauge one, e.g, a
faster convergence in simulations with intense and longwavelength lasers, and
the feasibility of exterior complex scaling as an absorbing boundary. The
reformulated TDCIS method is applied to an exactly solvable model of
onedimensional 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
highharmonic generation.

Timedependent coupledcluster method with timevarying orbital functions,
called timedependent optimized coupled cluster (TDOCC) 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 (TDOCCD) and double and triple excitations (TDOCCDT)
within the optimized active orbitals. The present method is size extensive and
gauge invariant, a polynomial costscaling alternative to the timedependent
multiconfiguration selfconsistentfield method. The first application of the
TDOCC method to intenselaser driven correlated electron dynamics in Ar atom
is reported.

We present a numerical implementation of the infiniterange exterior complex
scaling (irECS) [Phys. Rev. A 81, 053845 (2010)] as an efficient absorbing
boundary to the timedependent completeactivespace selfconsistent field
(TDCASSCF) method [Phys. Rev. A 94, 023405 (2016)] for multielectron atoms
subject to an intense laser pulse. We introduce GaussLaguerreRadau quadrature
points to construct discrete variable representation basis functions in the
last radial finite element extending to infinity. This implementation is
applied to strongfield ionization and highharmonic 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 highharmonic generation from Ne,
for example, 80% cost reduction is achieved, compared to a maskfunction
absorption boundary.

We investigate the multielectron effects on highharmonic generation from
solidstate materials using the timedependent HartreeFock theory. We find
qualitative change in harmonic spectra, in particular, multipleplateau
formation at significantly lower laser intensities than within the
independentelectron approximation. We reveal its origin in terms of interband
polarization, i.e, electronhole 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
freeelectron 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 timedependent multiconfiguration
selfconsistentfield 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 timedependent spinorbitals prepared for each particle kind.
We derive equations of motion followed by configurationinteraction (CI)
coefficients and spinorbitals for general, not restricted to fullCI,
configuration spaces. The present method provides a flexible framework for the
firstprinciples 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 highharmonic generation
(HHG), using allelectron firstprinciples simulations for a onedimensional
(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
threedimensional 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 singleelectron nature in most cases.

We theoretically study highharmonic generation (HHG) from solids driven by
intense laser pulses using a onedimensional model periodic crystal. By
numerically solving the timedependent Schr\"{o}dinger equation directly on a
realspace grid, we successfully reproduce experimentally observed unique
features of solidstate HHG such as the linear cutoffenergy scaling and the
sudden transition from a single to multipleplateau structure. Based on the
simulation results, we propose a simple model that incorporates
vectorpotentialinduced intraband displacement, interband tunneling, and
recombination with the valenceband hole. One key parameter is the
valleytopeak amplitude of the pulse vector potential, which determines the
crystal momentum displacement during the half cycle. When the maximum
peaktovalley 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 multipleplateau 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
timefrequency structure of HHG without elaborate quantummechanical
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 gasphase case. Our model can be viewed as a solidstate and
momentumspace counterpart of the familiar threestep model, highly successful
for gasphase HHG, and provide a unified basis to understand HHG from
solidstate materials and gaseous media.

The accurate description of the nonlinear response of manyelectron systems
to stronglaser 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 threedimensional
implementation of the timedependent twoparticle reduced density matrix
(TD2RDM) method for manyelectron atoms. We benchmark this approach by a
comparison with multiconfigurational timedependent HartreeFock (MCTDHF)
results for the harmonic spectra of beryllium and neon. We show that the
TD2RDM is very wellsuited to describe the nonlinear atomic response and to
reveal the influence of electroncorrelation effects.

We present the numerical implementation of the timedependent
completeactivespace selfconsistentfield (TDCASSCF) 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 gaugeinvariant frozencore approximation, an efficient evaluation
of the Coulomb mean field scaling linearly with the number of basis functions,
and a splitoperator method specifically designed for stable propagation of
stiff spatial derivative operators. We apply this method to highharmonic
generation in helium, beryllium, and neon and explore the role of electron
correlations.

We report a threedimensional numerical implementation of multiconfiguration
timedependent HartreeFock (MCTDHF) based on a multiresolution Cartesian
grid, with no need to assume any symmetry of molecular structure. We
successfully compute highharmonic generation (HHG) of H2 and H2O. The present
implementation will open a way to the firstprinciple theoretical study of
intensefield and attosecondpulse induced ultrafast phenomena in general
molecules.

In parallel with the evolution of femtosecond and attosecond laser as well as
freeelectron 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
wavefunctionbased numerical approaches to simulate multielectron dynamics in
atoms and molecules driven by intense longwavelength and/or ultrashort
shortwavelength laser pulses. Direct solution of the timedependent
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 electronelectron correlation.
Multiconfiguration selfconsistentfield (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. Timedependent
extension of the $R$matrix theory, originally develop for electronatom
collision, can realistically and accurately describe laserdriven complex
multielectron atoms.

The timedependent multiconfiguration selfconsistentfield method based on
the occupationrestricted multiple active space model is proposed (TDORMAS)
for multielectron dynamics in intense laser fields. Extending the previously
proposed timedependent completeactivespace selfconsistentfield method
[TDCASSCF; Phys. Rev. A, {\bf 88}, 023402 (2013)], which divides the occupied
orbitals into core and active orbitals, the TDORMAS 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
largeactivespace 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 timedependent variational principle, and an
efficient algorithm is proposed to solve for the orbital time derivatives.
Indepth descriptions of the computational implementation are given in a
readily programmable manner. The numerical application to the onedimensional
lithium hydride cluster models demonstrates that the high flexibility of the
TDORMAS framework allows for the costeffective simulations of multielectron
dynamics, by exploiting systematic series of approximations to the TDCASSCF
method.

Describing timedependent manybody systems where correlation effects play an
important role remains a major theoretical challenge. In this paper we develop
a timedependent manybody theory that is based on the twoparticle reduced
density matrix (2RDM). We develop a closed equation of motion for the 2RDM
employing a novel reconstruction functional for the threeparticle reduced
density matrix (3RDM) 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 longrange Coulomb interactions we employ
the onedimensional model for lithium hydride (LiH) in strong infrared laser
fields. We probe both oneparticle observables such as the timedependent
dipole moment and twoparticle observables such as the pair density and mean
electronelectron interaction energy. Our results are in very good agreement
with numerically exact solutions for the $N$electron wavefunction obtained
from the multiconfigurational timedependent HartreeFock method.

We theoretically investigate the mechanism of enhanced ionization in
twoelectron molecules by analyzing Bohmian trajectories for a onedimensional
H2 in an intense laser field. We identify both types of ionizing trajectories
corresponding to the ejection from the upfield and downfield 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
twocolor twophoton nearthreshold ionization of hydrogen and noble gas (He,
Ne, and Ar) atoms by a combined action of femtosecond extreme ultraviolet (EUV)
and nearinfrared (IR) laser pulses. Using the secondorder timedependent
perturbation theory, we clarify how the twophoton ionization process depends
on EUVIR pulse delay and how it is connected to the interplay between resonant
and nonresonant ionization paths. Furthermore, by solving the timedependent
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 pulsedelay effect is more pronounced for $p$shell ionization than for
$s$shell ionization.

The timedependent completeactivespace selfconsistentfield (TDCASSCF)
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 frozencore (to model tightly bound electrons with no
response to the field), dynamicalcore (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 manyelectron systems. The classification into the
subspaces can be done flexibly, according to simulated physical situations and
desired accuracy, and the multiconfiguration timedependent HartreeFock
(MCTDHF) approach is included as a special case. To assess its performance, we
apply the TDCASSCF method to the ionization dynamics of onedimensional
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 TDCASSCF method
will open a way to the firstprinciple theoretical study of intensefield
induced ultrafast phenomena in realistic atoms and molecules.

We theoretically study the pulsewidth dependence of the photoelectron
angular distribution (PAD) from the resonanceenhanced twophoton single
ionization of He by femtosecond ($\lesssim 20$ fs) extremeultraviolet pulses,
based on the timedependent perturbation theory and simulations with the full
timedependent 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 abovethreshold 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 timedependent 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
highorder 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.