Source author record

Takeshi Sato

Takeshi Sato appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

21works
13topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

21 published item(s)

preprint2026arXiv

Time-dependent Hole States in Multiconfigurational Time-Dependent Hartree-Fock Approaches: A Time-Domain Generalization of Extended Koopmans' Theorem

We introduce a framework for resolving electron-hole dynamics within wavefunction-based multiconfigurational time-dependent Hartree-Fock (MCTDHF) theory. Central to this framework is a time-domain generalization of the extended Koopmans' theorem, which rigorously defines time-dependent hole states through single-electron removal. From this foundation, we prove the existence of exact equations of motion for time-dependent Dyson orbitals, enabling instantaneous construction of photofragments' reduced density matrices. The formalism further yields a systematic procedure to extract hole-resolved observables, such as channel-resolved photoelectron momentum distributions, directly from time-dependent \textit{ab initio} wavefunctions. As a demonstration, we employ an attosecond $ω-2ω$ laser strategy to control hole dynamics, thereby resolving a long-standing challenge in MCTDHF simulations. This advance opens a pathway for exploring correlated multielectron dynamics in atoms and molecules under ultrafast laser fields.

preprint2025arXiv

Time-dependent Hole States in Multiconfigurational Time-Dependent Hartree-Fock Approaches: Applications in Photoionization of Water Molecule

By simulating the real-time multielectron wavefunction with the multi-configurational time-dependent Hartree-Fock (MCTDHF) approach, we conduct an \textit{ab initio} study of the single-photon ionization process of a body-fixed water molecule ($\mathrm{H_2O}$) driven by attosecond pulses. To this end, we present a full-dimensional implementation of the MCTDHF method based on one-center expansions, allowing for the simulation of arbitrarily polarized lasers and multi-center polyatomic potentials. With a rigorous definition of the time-dependent hole state (TDHS) using the time-domain generalization of extended Koopmans' theorem (TD-EKT), we derive the reduced ion density matrix within the MCTDHF framework, which inherently encodes the total and channel-resolved photoionization cross sections of $\mathrm{H_2O}$. The cross sections obtained are benchmarked against existing experimental and theoretical results, validating the TDHS formalism. Furthermore, by adjusting the phase delay and intensity ratio of a pair of orthogonally polarized attosecond pulses, we explore the ultrafast control of attosecond coherence between electronic states of $\mathrm{H_2O^+}$.

preprint2023arXiv

Efficient simulation of multielectron dynamics in molecules under intense laser pulses: Implementation of the multiconfiguration time-dependent Hartree-Fock method based on the adaptive finite element method

We present an implementation of the multiconfiguration time-dependent Hartree-Fock method based on the adaptive finite element method for molecules under intense laser pulses. For efficient simulations, orbital functions are propagated by a stable propagator using the short iterative Arnoldi scheme and our implementation is parallelized for distributed memory computing. This is demonstrated by simulating high-harmonic generation from a water molecule and achieves a simulation of multielectron dynamics with overwhelmingly less computational time, compared to our previous work.

preprint2022arXiv

Time-dependent optimized coupled-cluster method with doubles and perturbative triples [TD-OCCD(T)] for first principles simulation of multielectron dynamics

We report the formulation of a new, cost-effective approximation method in the time-dependent optimized coupled-cluster (TD-OCC) framework [T. Sato et al., J. Chem. Phys. 148, 051101 (2018)] for first-principles simulations of multielectron dynamics in an intense laser field. The method, designated as TD-OCCD(T), is a time-dependent, orbital-optimized extension of the "gold-standard" CCSD(T) method in the ground-state electronic structure theory. The equations of motion for the orbital functions and the coupled-cluster amplitudes are derived based on the real-valued time-dependent variational principle using the fourth-order Lagrangian. The TD-OCCD(T) is size extensive and gauge invariant, and scales as O(N^7) with respect to the number of active orbitals N. The pilot application of the TD-OCCD(T) method to the strong-field ionization and high-order harmonic generation from a Kr atom is reported in comparison with the results of the previously developed methods, such as the time-dependent complete-active-space self-consistent field (TD-CASSCF), TD-OCC with double and triple excitations (TD-OCCDT), TD-OCC with double excitations (TD-OCCD), and the time-dependent Hartree-Fock (TDHF) methods.

preprint2022arXiv

Variational Quantum Simulation for Periodic Materials

We present a quantum-classical hybrid algorithm that simulates electronic structures of periodic systems such as ground states and quasiparticle band structures. By extending the unitary coupled cluster (UCC) theory to describe crystals in arbitrary dimensions, for a hydrogen chain, we numerically demonstrate that the UCC ansatz implemented on a quantum circuit can be successfully optimized with a small deviation from the exact diagonalization over the entire range of the potential energy curves. Furthermore, by using the quantum subspace expansion method, in which we truncate the Hilbert space within the linear response regime from the ground state, the quasiparticle band structure is computed as charged excited states. Our work establishes a powerful interface between the rapidly developing quantum technology and modern material science.

preprint2020arXiv

Learning the constitutive relation of polymeric flows with memory

We develop a learning strategy to infer the constitutive relation for the stress of polymeric flows with memory. We make no assumptions regarding the functional form of the constitutive relations, except that they should be expressible in differential form as a function of the local stress- and strain-rate tensors. In particular, we use a Gaussian Process regression to infer the constitutive relations from stress trajectories generated from small-scale (fixed strain-rate) microscopic polymer simulations. For simplicity, a Hookean dumbbell representation is used as a microscopic model, but the method itself can be generalized to incorporate more realistic descriptions. The learned constitutive relation is then used to perform macroscopic flow simulations, allowing us to update the stress distribution in the fluid in a manner that accounts for the microscopic polymer dynamics. The results using the learned constitutive relation are in excellent agreement with full Multi-Scale Simulations, which directly couple micro/macro degrees of freedom, as well as the exact analytical solution given by the Maxwell constitutive relation. We are able to fully capture the history dependence of the flow, as well as the elastic effects in the fluid. We expect the proposed learning/simulation approach to be used not only to study the dynamics of entangled polymer flows, but also for the complex dynamics of other Soft Matter systems, which possess a similar hierarchy of length- and time-scales.

preprint2020arXiv

Multiscale simulation of a polymer melt flow between two coaxial cylinders under nonisothermal conditions

We successfully extend a multiscale simulation (MSS) method to nonisothermal well-entangled polymer melt flows between two coaxial cylinders. In the multiscale simulation, the macroscopic flow system is connected to a number of microscopic systems through the velocity gradient tensor, stress tensor and temperature. At the macroscopic level, in addition to the momentum balance equation, we consider the energy balance equation, where heat generation plays an important role not only in the temperature distribution but also in the flow profile. At the microscopic level, a dual slip-link model is employed for well-entangled polymers. To incorporate the temperature effect into the microscopic systems, we used the time-temperature superposition rule for the slip-link model, in which the temperature dependence of the parameters is not known; on the other hand, the way to take into account the temperature effect in the macroscopic equations has been well established. We find that the extended multiscale simulation method is quite effective in revealing the relation between nonisothermal polymeric flows for both steady and transient cases and the microscopic states of polymer chains expressed by primitive paths and slip-links. It is also found that the temperature-dependent reptation-time-based Weissenberg number is a suitable measure for understanding the extent of the polymer chain deformation in the range of the shear rate used in this study.

preprint2020arXiv

Strong-Field Electron Dynamics in Solids

Solid-state materials have recently emerged as a new stage of strong-field physics and attosecond science. The mechanism of the electron dynamics driven by an ultrashort intense laser pulse is under intensive discussion. Here we theoretically discuss momentum-space strong-field electron dynamics in graphene and crystalline dielectrics and semiconductors. First, within massless Dirac fermion and tight-binding models for graphene, we rigorously derive intraband displacement and interband transition, which form the basis for understanding solid-state strong-field physics including high-harmonic generation (HHG). Then, based on the time-dependent Schrödinger equation for a one-dimensional model crystal, we introduce a simple, multiband, momentum-space three-step model that incorporates intraband displacement, interband tunneling, and recombination with a valence band hole. We also analyze how the model is modified by electron-hole interaction. Finally, actual three-dimensional materials are investigated. We present a time-dependent density-matrix method whose results for HHG are compared with experimental measurement results. Moreover, we describe the dynamical Franz-Keldysh effect in femtosecond time resolution, i.e., the time-dependent modulation of a dielectric function under an intense laser field, using a real-time time-dependent density functional theory.

preprint2020arXiv

Study of laser-driven multielectron dynamics of Ne atom using time-dependent optimized second-order many-body perturbation theory

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 ($5\times10^{14}$, $8\times10^{14}$ , and $1\times10^{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 $5\times10^{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.

preprint2020arXiv

Time-dependent optimized coupled-cluster method for multielectron dynamics II. A coupled electron-pair approximation

We report the implementation of a cost-effective approximation method within the framework of time-dependent optimized coupled-cluster (TD-OCC) method [J. Chem. Phys. 148, 051101 (2018)] for real-time simulations of intense laser-driven multielectron dynamics. The method, designated as TD-OCEPA0, is a time-dependent extension of the simplest version of the coupled-electron pair approximation with optimized orbitals [J. Chem. Phys. 139, 054104 (2013)]. It is size extensive, gauge invariant, and computationally much more efficient than the TD-OCC with double excitations (TD-OCCD). We employed this method to simulate the electron dynamics in Ne and Ar atoms exposed to intense near infrared laser pulses with various intensities. The computed results, including high-harmonic generation spectra and ionization yields, are compared with those of various other methods ranging from uncorrelated time-dependent Hartree-Fock (TDHF) to fully-correlated (within the active orbital space) time-dependent complete-active-space self-consistent-field (TD-CASSCF). The TD-OCEPA0 results show a good agreement with TD-CASSCF ones for moderate laser intensities. For higher intensities, however, TD-OCEPA0 tends to overestimate the correlation effect, as occasionally observed for CEPA0 in the ground-state correlation energy calculations.

preprint2020arXiv

Time-dependent optimized coupled-cluster method for multielectron dynamics III: A second-order many-body perturbation approximation

We report successful implementation of the time-dependent second-order many-body perturbation theory using optimized orthonormal orbital functions called time-dependent optimized second-order many-body perturbation theory [TD-OMP2] to reach out to relatively larger chemical systems for the study of intense-laser-driven multielectron dynamics. We apply this method to strong-field ionization and high-order harmonic generation (HHG) of Ar. The calculation results are benchmarked against ab initio time-dependent complete-active-space self-consistent field (TD-CASSCF), time-dependent optimized coupled-cluster double (TD-OCCD), and time-dependent Hartree-Fock (TDHF) methods, as well as a single active electron (SAE) model to explore the role of electron correlation.

preprint2016arXiv

Implementation of the multiconfiguration time-dependent Hatree-Fock method for general molecules on a multi-resolution Cartesian grid

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.

preprint2016arXiv

Multiscale simulation of polymer melt spinning by using the dumbbell model

We investigated the spinning process of a polymeric material by using a multiscale simulation method which connects the macroscopic and microscopic states through the stress and strain-rate tensor fields, by using Lagrangian particles (filled with polymer chains) along the spinning line. We introduce a large number of Lagrangian fluid particles into the fluid, each containing Np-Hookean-dumbbells to mimic the polymer chains (Np=$10^4$), which is equivalent to the upper convected Maxwell fluid in the limit that Np$\rightarrow \infty$. Depending on the Reynolds number Re, we studied the dynamical behaviors of fibers for the (a) Re =0 and (b) finite Re cases, for different draw ratios Dr, ranging from 10 to 30, and two typical Deborah numbers De=$10^{-3}$ and De=$10^{-2}$. In the limit Re$\rightarrow$ 0 (a), as the Deborah number De increases, the elastic effect makes the system stable. At finite Re (b), we found that inertial effects play an important role in determining the dynamical behavior of the spinning process, and for Dr=$10^{-2}$ the system is quite stable, at least up to a draw ratio of Dr=30. We also found that the fiber velocity and cross section area are determined solely by the draw ratio. By comparing the velocity and cross section area profiles with the end-point distribution for the dumbbell connective vectors, for dumbbells located in Lagrangian particles along typical places along the spinning line, we show that our multiscale simulation method successfully bridges the microscopic state of the system with its simultaneous macroscopic flow behavior. It is also confirmed that the present schemes gives good agreements with the results obtained by the Maxwell constitutive equation.

preprint2016arXiv

Time-dependent complete-active-space self-consistent-field method for atoms: Application to high-harmonic generation

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.

preprint2015arXiv

A Review on Ab Initio Approaches for Multielectron Dynamics

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ödinger 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.

preprint2015arXiv

Time-dependent multiconfiguration self-consistent-field method based on occupation restricted multiple active space model for multielectron dynamics in intense laser fields

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.

preprint2014arXiv

Analysis of strong-field enhanced ionization of molecules using Bohmian trajectories

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.

preprint2014arXiv

The structure of approximate two electron wavefunctions in intense laser driven ionization dynamics

The structure of approximate two electron wavefunction is deeply investigated, both theoretically and numerically, in the strong-field driven ionization dynamics. Theoretical analyses clarify that for two electron singlet systems, the previously proposed time-dependent extended Hartree-Fock (TD-EHF) method [Phys. Rev. A 51, 3999 (1995)] is equivalent to the multiconfiguration time-dependent Hartree-Fock method with two occupied orbitals. The latter wavefunction is further transformed into the natural expansion form, enabling the direct propagation of the natural orbitals (NOs). These methods, as well as the conventional time-dependent Hartree-Fock (TDHF) method, are numerically assessed for the description of ionization dynamics of one-dimensional helium atom model. This numerical analysis (i) explains the reason behind the well-known failure of TDHF method to describe tunneling ionization, (ii) demonstrates the interpretive power of the TD-EHF wavefunction both in the original nonorthogonal and the NO-based formulations, and (iii) highlights different manifestations of the electron correlation (effect beyond the single determinant description), in tunneling ionization, high harmonic generation, and nonsequential double ionization. Possible extensions of the NO basis approach to multielectron systems are briefly discussed.

preprint2013arXiv

Time-dependent complete-active-space self-consistent field method for multielectron dynamics in intense laser fields

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.

preprint2012arXiv

Recent progress in silica aerogel Cherenkov radiator

In this paper, we present recent progress in the development of hydrophobic silica aerogel as a Cherenkov radiator. In addition to the conventional method, the recently developed pin-drying method for producing high-refractive-index aerogels with high transparency was studied in detail. Optical qualities and large tile handling for crack-free aerogels were investigated. Sufficient photons were detected from high-performance aerogels in a beam test.

preprint2003arXiv

D-brane Dynamics and Creations of Open and Closed Strings after Recommbination

A quantum-mechanical technique is used within the framework of U(2) super-Yang-Mills theory to investigate what happens in the process after recombination of two D-p-branes at one angle. Two types of initial conditions are considered, one of which with $p=4$ is a candidate of inflation mechanism. It is observed that the branes' shapes come to have three extremes due to localization of tachyon condensation. Furthermore, open string pairs connecting the decaying D-branes are shown to be created; most part of the released energy is used to create them. It also strongly suggests that creation of closed strings happens afterward. Closed strings as gravitational radiation from the D-branes are also shown to be created. A few speculations are also given on implications of the above phenomena for an inflation model.