Source author record

Naoto Tsuji

Naoto Tsuji 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

25works
7topics
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

25 published item(s)

preprint2022arXiv

Optical response of the Leggett mode in multiband superconductors in the linear response regime

We study optical responses of Leggett modes in multiband superconductors in the linear response regime. The Leggett mode is a collective mode unique to multiband superconductors that arises from relative phase fluctuations of superconducting orders for different bands. We use the Ginzburg-Landau (GL) description to study the collective modes in multiband systems. We find that multiband superconductors generally allow a linear coupling between the Leggett mode and external electric fields due to the presence of a cross term between different components of superconducting orders in the GL theory. The presence of a linear coupling for the Leggett mode is in sharp contrast with the absence of that for Higgs (amplitude) modes in single-band superconductors that only support nonlinear optical responses such as third harmonic generation (THG). We further confirm such a linear coupling in multiband superconductors by a more microscopic description based on a diagrammatic approach. We study the collective modes within the random phase approximation (RPA) and compute their contribution to the linear optical conductivity. These findings suggest a new route to observe the Leggett mode by optical absorption.

preprint2021arXiv

Tachyonic and Plasma Instabilities of $η$-Pairing States Coupled to Electromagnetic Fields

Cooper pairs featuring a nonzero center-of-mass crystal momentum ${\boldsymbol Q}=(π,π, \dots)$ and an off-diagonal long-range order ($η$-pairing states) constitute exact eigenstates of a Hubbard model [C. N. Yang, Phys. Rev. Lett. 63, 2144 (1989)]. Here we show that the $η$-pairing states are rendered dynamically unstable via coupling to dynamical electromagnetic fields. The instability is caused by "tachyonic" electromagnetic fields and unstable plasma modes for attractive and repulsive interactions, respectively. The typical time scale of the growth of the instability is of the order of femtoseconds for electron systems in solids, which places a strict bound on the lifetime of the $η$-pairing states. The decay of the $η$-pairing states leads to enhanced light emission with frequencies shifted from the Hubbard interaction strength for the repulsive case and unattenuated electromagnetic penetration for the attractive case.

preprint2020arXiv

Dynamical sign reversal of magnetic correlations in dissipative Hubbard models

In quantum magnetism, the virtual exchange of particles mediates an interaction between spins. Here, we show that an inelastic Hubbard interaction fundamentally alters the magnetism of the Hubbard model due to dissipation in spin-exchange processes, leading to sign reversal of magnetic correlations in dissipative quantum dynamics. This mechanism is applicable to both fermionic and bosonic Mott insulators, and can naturally be realized with ultracold atoms undergoing two-body inelastic collisions. The dynamical reversal of magnetic correlations can be detected by using a double-well optical lattice or quantum-gas microscopy, the latter of which facilitates the detection of the magnetic correlations in one-dimensional systems because of spin-charge separation. Our results open a new avenue toward controlling quantum magnetism by dissipation.

preprint2019arXiv

Higgs Mode in Superconductors

When a continuous symmetry of a physical system is spontaneously broken, two types of collective modes typically emerge: the amplitude and phase modes of the order-parameter fluctuation. For superconductors, the amplitude mode is recently referred to as the ''Higgs mode'' as it is a condensed-matter analogue of a Higgs boson in particle physics. Higgs mode is a scalar excitation of the order parameter, distinct from charge or spin fluctuations, and thus does not couple to electromagnetic fields linearly. This is why the Higgs mode in superconductors has evaded experimental observations over a half century after the initial theoretical prediction, except for a charge-density-wave coexisting system. With the advance of nonlinear and time-resolved terahertz spectroscopy techniques, however, it has become possible to study the Higgs mode through the nonlinear light-Higgs coupling. In this review, we overview recent progresses on the study of the Higgs mode in superconductors.

preprint2016arXiv

Brillouin-Wigner theory for high-frequency expansion in periodically driven systems: Application to Floquet topological insulators

We construct a systematic high-frequency expansion for periodically driven quantum systems based on the Brillouin-Wigner (BW) perturbation theory, which generates an effective Hamiltonian on the projected zero-photon subspace in the Floquet theory, reproducing the quasienergies and eigenstates of the original Floquet Hamiltonian up to desired order in $1/ω$, with $ω$ being the frequency of the drive. The advantage of the BW method is that it is not only efficient in deriving higher-order terms, but even enables us to write down the whole infinite series expansion, as compared to the van Vleck degenerate perturbation theory. The expansion is also free from a spurious dependence on the driving phase, which has been an obstacle in the Floquet-Magnus expansion. We apply the BW expansion to various models of noninteracting electrons driven by circularly polarized light. As the amplitude of the light is increased, the system undergoes a series of Floquet topological-to-topological phase transitions, whose phase boundary in the high-frequency regime is well explained by the BW expansion. As the frequency is lowered, the high-frequency expansion breaks down at some point due to band touching with nonzero-photon sectors, where we find numerically even more intricate and richer Floquet topological phases spring out. We have then analyzed, with the Floquet dynamical mean-field theory, the effects of electron-electron interaction and energy dissipation. We have specifically revealed that phase transitions from Floquet-topological to Mott insulators emerge, where the phase boundaries can again be captured with the high-frequency expansion.

preprint2016arXiv

First-Principles Design of a Half-Filled Flat Band of the Kagome Lattice in Two-Dimensional Metal-Organic Frameworks

We design from first principles a new type of two-dimensional metal-organic frameworks (MOFs) using phenalenyl-based ligands to exhibit a half-filled flat band of the kagome lattice, which is one of the lattice family that shows Lieb-Mielke-Tasaki's flat-band ferromagnetism. Among various MOFs, we find that $\textit{trans}$-Au-THTAP(trihydroxytriaminophenalenyl) has such an ideal band structure, where the Fermi energy is adjusted right at the flat band due to unpaired electrons of radical phenalenyl. The spin-orbit coupling opens a band gap giving a non-zero Chern number to the nearly flat band, as confirmed by the presence of the edge states in first-principles calculations and by fitting to the tight-binding model. This is a novel and realistic example of a system in which a nearly flat band is both ferromagnetic $\textit{and}$ topologically non-trivial.

preprint2016arXiv

Nonequilibrium Dynamical Cluster Approximation study of the Falicov-Kimball Model

We use a nonequilibrium implementation of the dynamical cluster approximation (DCA) to study the effect of short-range correlations on the dynamics of the two-dimensional Falicov-Kimball model after an interaction quench. As in the case of single-site dynamical mean field theory, thermalization is absent in DCA simulations, and for quenches across the metal-insulator boundary, nearest-neighbor charge correlations in the nonthermal steady state are found to be larger than in the thermal state with identical energy. We investigate to what extent it is possible to define an effective temperature of the trapped state after a quench. Based on the ratio between the lesser and retarded Green's function we conclude that a roughly thermal distribution is reached within the energy intervals corresponding to the momentum-patch dependent subbands of the spectral function. The effectively different chemical potentials of these distributions however lead to a very hot, or even negative, effective temperature in the energy intervals between these subbands.

preprint2016arXiv

Nonlinear light-Higgs coupling in superconductors beyond BCS: Effects of the retarded phonon-mediated interaction

We study the contribution of the Higgs amplitude mode on the nonlinear optical response of superconductors beyond the BCS approximation by taking into account the retardation effect in the phonon-mediated attractive interaction. To evaluate the vertex correction in nonlinear optical susceptibilities that contains the effect of collective modes, we propose an efficient scheme which we call the "dotted DMFT" based on the nonequilibrium dynamical mean-field theory (nonequilibrium DMFT) to go around the difficulty of solving the Bethe-Salpeter equation and analytical continuation. The vertex correction is represented by the derivative of the self-energy with respect to the external driving field, which is self-consistently determined by the differentiated ("dotted") DMFT equations. We apply the method to the Holstein model, a prototypical electron-phonon-coupled system, to calculate the susceptibility for the third-harmonic generation including the vertex correction. The results show that, in sharp contrast to the BCS theory, the Higgs mode can contribute to the third-harmonic generation for general polarization of the laser field with an order of magnitude comparable to the contribution from the pair breaking or charge density fluctuations. The physical origin is traced back to the nonlinear resonant light-Higgs coupling, which has been absent in the BCS approximation.

preprint2015arXiv

Electronic structure of helicoidal graphene: massless Dirac particles on a curved surface with a screw symmetry

Massless Dirac particles on the helicoid are theoretically investigated. With its possible application being helical graphene, we explore how the peculiarities of Dirac particles appear on the curved, screw-symmetric surface. The zweibein is used to derive the massless Dirac equation on the helicoid, as well as general curved surfaces. Bound states of massless Dirac electrons are shown to be absent on the helicoid, and then the scattering probabilities and the phase shifts on the surface are obtained from numerically calculated wave functions. We find the local density of states and the phase shifts idiosyncratic especially around the axis of the helicoid. Bound states of massive Dirac electrons on the surface are also shown to be absent as an extension of the above result on massless Dirac electrons.

preprint2015arXiv

FLEX+DMFT approach to the $d$-wave superconducting phase diagram of the two-dimensional Hubbard model

The dynamical mean-field theory (DMFT) combined with the fluctuation exchange (FLEX) method, namely FLEX+DMFT, is an approach for correlated electron systems to incorporate both local and non-local long-range correlations in a self-consistent manner. We formulate FLEX+DMFT in a systematic way starting from a Luttinger-Ward functional, and apply it to study the $d$-wave superconductivity in the two-dimensional repulsive Hubbard model. The critical temperature ($T_c$) curve obtained in the FLEX+DMFT exhibits a dome structure as a function of the filling, which has not been clearly observed in the FLEX approach alone. We trace back the origin of the dome to the local vertex correction from DMFT that renders a filling dependence in the FLEX self-energy. We compare the results with those of GW+DMFT, where the $T_c$-dome structure is qualitatively reproduced due to the same vertex correction effect, but a crucial difference from FLEX+DMFT is that $T_c$ is always estimated below the Néel temperature in GW+DMFT. The single-particle spectral function obtained with FLEX+DMFT exhibits a double-peak structure as a precursor of the Hubbard bands at temperature above $T_c$.

preprint2015arXiv

Interaction-Driven Topological Insulator in Fermionic Cold Atoms on an Optical Lattice: A Design with a Density Functional Formalism

We design an interaction-driven topological insulator for fermionic cold atoms in an optical lattice, that is, we pose the question of whether we can realize in a continuous space a spontaneous symmetry breaking induced by the inter-atom interaction into a topological Chern insulator. Such a state, sometimes called a "topological Mott insulator", has yet to be realized in solid-state systems, since this requires, in the tight-binding model, large offsite interactions on top of a small onsite interaction. Here we overcome the difficulty by introducing a spin-dependent potential, where a spin-selective occupation of fermions in $A$ and $B$ sublattices makes the onsite interaction Pauli-forbidden, while a sizeable inter-site interaction is achieved by a shallow optical potential with a large overlap between neighboring Wannier orbitals. This puts the system away from the tight-binding model, so that we adopt the density functional theory for cold-atoms, here extended to accommodate non-collinear spin structures emerging in the topological regime, to quantitatively demonstrate the phase transition to the topological Mott insulator.

preprint2015arXiv

Nonequilibrium dynamical cluster theory

We study the effect of spatially nonlocal correlations on the nonequilibrium dynamics of interacting fermions by constructing the nonequilibrium dynamical cluster theory, a cluster generalization of the nonequilibrium dynamical mean-field theory (DMFT). The formalism is applied to interaction quenches in the Hubbard model in one and two dimensions, and the results are compared with data from single-site DMFT, the time-dependent density matrix renormalization group, and lattice perturbation theory. Both in one and two dimensions the double occupancy quickly thermalizes, while the momentum distribution relaxes only on much longer time scales. For the two-dimensional square lattice we find a strongly momentum-dependent evolution of the momentum distribution around the Fermi energy, with a much faster relaxation near the momenta $(0,π)$ and $(π,0)$ than near $(π/2,π/2)$. This result is interpreted as reflecting the momentum-anisotropic quasiparticle lifetime of the marginal Fermi liquid. The method is further applied to the two-dimensional Hubbard model driven by a dc electric field, where the damping of the Bloch oscillation of the current is found to be less effective than predicted by DMFT and lattice perturbation theory.

preprint2015arXiv

Theory of Anderson pseudospin resonance with Higgs mode in superconductors

A superconductor illuminated by an ac electric field with frequency $Ω$ is theoretically found to generate a collective precession of Anderson's pseudospins, and hence a coherent amplitude oscillation of the order parameter, with a doubled frequency $2Ω$ through a nonlinear light-matter coupling. We provide a fundamental theory, based on the mean-field formalism, to show that the induced pseudospin precession resonates with the Higgs amplitude mode of the superconductor at $2Ω=2Δ$ with $2Δ$ being the superconducting gap. The resonant precession is accompanied by a divergent enhancement of the third-harmonic generation (THG). By decomposing the THG susceptibility into the bare one and vertex correction, we find that the enhancement of the THG cannot be explained by individual quasiparticle excitations (pair breaking), so that the THG serves as a smoking gun for an identification of the collective Higgs mode. We further explore the effect of electron-electron scattering on the pseudospin resonance by applying the nonequilibrium dynamical mean-field theory to the attractive Hubbard model driven by ac electric fields. The result indicates that the pseudospin resonance is robust against electron correlations, although the resonance width is broadened due to electron scattering, which determines the lifetime of the Higgs mode.

preprint2014arXiv

Interaction quench in the Holstein model: Thermalization crossover from electron- to phonon-dominated relaxation

We study the relaxation of the Holstein model after a sudden switch-on of the interaction by means of the nonequilibrium dynamical mean field theory, with the self-consistent Migdal approximation as an impurity solver. We show that there exists a qualitative change in the thermalization dynamics as the interaction is varied in the weak-coupling regime. On the weaker interaction side of this crossover, the phonon oscillations are damped more rapidly than the electron thermalization timescale, as determined from the relaxation of the electron momentum distribution function. On the stronger interaction side, the relaxation of the electrons becomes faster than the phonon damping. In this regime, despite long-lived phonon oscillations, a thermalized momentum distribution is realized temporarily. The origin of the thermalization crossover found here is traced back to different behaviors of the electron and phonon self-energies as a function of the electron-phonon coupling. In addition, the importance of the phonon dynamics is demonstrated by comparing the self-consistent Migdal results with those obtained with a simpler Hatree-Fock impurity solver that neglects the phonon self-energy. The latter scheme does not properly describe the evolution and thermalization of isolated electron-phonon systems.

preprint2014arXiv

Nonequilibrium dynamical mean-field theory and its applications

The study of nonequilibrium phenomena in correlated lattice systems has developed into an active and exciting branch of condensed matter physics. This research field provides rich new insights that could not be obtained from the study of equilibrium situations, and the theoretical understanding of the physics often requires the development of new concepts and methods. On the experimental side, ultra-fast pump-probe spectroscopies enable studies of excitation and relaxation phenomena in correlated electron systems, while ultra-cold atoms in optical lattices provide a new way to control and measure the time-evolution of interacting lattice systems with a vastly different characteristic timescale compared to electron systems. A theoretical description of these phenomena is challenging because, firstly, we have to compute the quantum-mechanical time-evolution of many-body systems out of equilibrium, and secondly, deal with strong-correlation effects which can be of nonperturbative nature. In this review, we discuss the nonequilibrium extension of the dynamical mean field theory (DMFT), which treats quantum fluctuations in the time domain and works directly in the thermodynamic limit. The method reduces the complexity of the calculation via a mapping to a self-consistent impurity problem. Particular emphasis is placed on a detailed derivation of the formalism, and on a discussion of numerical techniques, which enable solutions of the effective nonequilibrium DMFT impurity problem. We summarize the insights gained into the properties of the infinite-dimensional Hubbard model under strong non-equilibrium conditions. These examples illustrate the current ability of the theoretical framework to reproduce and understand fundamental nonequilibrium phenomena, such as the dielectric breakdown of Mott insulators, photo-doping, and collapse-and-revival oscillations in quenched systems.

preprint2014arXiv

Nonequilibrium dynamical mean-field theory based on weak-coupling perturbation expansions: Application to dynamical symmetry breaking in the Hubbard model

We discuss the general formalism and validity of weak-coupling perturbation theory as an impurity solver for nonequilibrium dynamical mean-field theory. The method is implemented and tested in the Hubbard model, using expansions up to fourth order for the paramagnetic phase at half filling and third order for the antiferromagnetic and paramagnetic phase away from half filling. We explore various types of weak-coupling expansions and examine the accuracy and applicability of the methods for equilibrium and nonequilibrium problems. We find that in most cases an expansion of local self-energy diagrams including all the tadpole diagrams with respect to the Weiss Green's function (bare-diagram expansion) gives more accurate results than other schemes such as self-consistent perturbation theory using the fully interacting Green's function (bold-diagram expansion). In the paramagnetic phase at half filling, the fourth-order bare expansion improves the result of the second-order expansion in the weak-coupling regime, while both expansions suddenly fail at some intermediate interaction strength. The higher-order bare perturbation is especially advantageous in the antiferromagnetic phase near half filling. We use the third-order bare perturbation expansion within the nonequilibrium dynamical mean-field theory to study dynamical symmetry breaking from the paramagnetic to the antiferromagnetic phase induced by an interaction ramp in the Hubbard model. The results show that the order parameter, after an initial exponential growth, exhibits an amplitude oscillation around a nonthermal value followed by a slow drift toward the thermal value. The transient dynamics seems to be governed by a nonthermal critical point, associated with a nonthermal universality class, which is distinct from the conventional Ginzburg-Landau theory.

preprint2014arXiv

Supersolid Phase Accompanied by a Quantum Critical Point in the Intermediate Coupling Regime of the Holstein Model

We reveal that electron-phonon systems described by the Holstein model on a bipartite lattice exhibit, away from half filling, a supersolid (SS) phase characterized by coexisting charge order (CO) and superconductivity (SC), and an accompanying quantum critical point (QCP). The SS, demonstrated by the dynamical mean-field theory with a quantum Monte Carlo impurity solver, emerges in the intermediate-coupling regime, where the peak of the Tc dome is located and the metal-insulator crossover occurs. On the other hand, in the weak- and strong-coupling regimes the CO-SC boundary is of first order with no intervening SS phases. The QCP is associated with the continuous transition from SS to SC and characterized by a reentrant behavior of the SS around it. We further show that the SS-SC transition is hallmarked by diverging charge fluctuations and a kink (peak) in the superfluid density.

preprint2013arXiv

Dynamical mean-field analysis of ordered phases in the half-filled Holstein-Hubbard model

We investigate ordered phases (s-wave superconducting (SC), antiferromagnetic (AF) and charge ordered (CO) phases) of the half-filled Holstein-Hubbard model using the dynamical mean-field theory in combination with a continuous time quantum Monte Carlo impurity solver. We particularly focus on the parameter regime where the electron-electron Coulomb interaction and the phonon mediated retarded attractive interaction are both comparable to the bare electron bandwidth to understand the interplay of these interactions and the retardation effects due to finite phonon frequency ω_0. Results for the transition temperatures for SC and AF as a function of the effective interaction U_eff(=U-λwith U the Coulomb repulsion and λthe phonon-mediated attraction) show deviations from the the corresponding results in the anti-adiabatic limit. These are analyzed by comparing with an effective static model in the polaron picture. We also determine the phase diagram around U_eff=0 at low, but nonzero temperatures to show that paramagnetic phases appear between AF and CO when U and λare small or large, while AF and CO directly compete with each other at moderate U and λ.

preprint2013arXiv

Nonthermal antiferromagnetic order and nonequilibrium criticality in the Hubbard model

We study dynamical phase transitions from antiferromagnetic to paramagnetic states driven by an interaction quench in the fermionic Hubbard model using the nonequilibrium dynamical mean-field theory. We identify two dynamical transition points where the relaxation behavior qualitatively changes: one corresponds to the thermal phase transition at which the order parameter decays critically slowly in a power law $\propto t^{-1/2}$, and the other is connected to the existence of nonthermal antiferromagnetic order in systems with effective temperature above the thermal critical temperature. The frequency of the amplitude mode extrapolates to zero as one approaches the nonthermal (quasi)critical point, and thermalization is significantly delayed by the trapping in the nonthermal state. A slow relaxation of the nonthermal order is followed by a faster thermalization process.

preprint2013arXiv

Ordered phases in the Holstein-Hubbard model: Interplay of strong Coulomb interaction and electron-phonon coupling

We study the Holstein-Hubbard model at half filling to explore ordered phases including su- perconductivity (SC), antiferromagnetism (AF), and charge order (CO) in situations where the electron-electron and electron-phonon interactions are strong (~ electronic bandwidth). The model is solved in the dynamical mean-field approximation using a continuous-time quantum Monte Carlo impurity solver. We determine the superconducting transition temperature (Tc) and the SC order parameter and show that the phonon-induced retardation or the strong Coulomb interaction leads to a significant reduction and shift of the Tc dome if one interprets the system as having an effective static interaction Ueff given by the Hubbard U reduced by the phonon-mediated attraction in the static limit. This behavior is analyzed by comparison to an effective static model in the polaron representation with a reduced bandwidth. We also determine the finite-temperature phase diagram including AF and CO. In the moderate-coupling regime, there is a hysteretic region of AF and CO around Ueff = 0, while the two phases are separated by a paramagnetic metal in the weak-coupling regime and a paramagnetic insulator in the strong coupling regime.

preprint2012arXiv

Nonthermal symmetry broken states in the strongly interacting Hubbard model

We study the time evolution of the antiferromagnetic order parameter after interaction quenches in the Hubbard model. Using the nonequilibrium dynamical mean field formalism, we show that the system, after a quench from intermediate to strong interaction, is trapped in a nonthermal state which is reminiscent of a photo-doped state and protected by the slow decay of doublons. If the effective doping of this state is low enough, it exhibits robust antiferromagnetic order, even if the system is highly excited and the thermal state thus expected to be paramagnetic. We comment on the implication of our findings for the stability of nonthermal superconducting states.

preprint2012arXiv

Repulsion-to-attraction transition in correlated electron systems triggered by a monocycle pulse

We study the time evolution of the Hubbard model driven by a half-cycle or monocycle pulsed electric field F(t) using the nonequilibrium dynamical mean-field theory. We find that for properly chosen pulse shapes the electron-electron interaction can be effectively and permanently switched from repulsive to attractive if there is no energy dissipation. The physics behind the interaction conversion is a nonadiabatic shift $δ$ of the population in momentum space. When $δ\simπ$, the shifted population relaxes to a negative-temperature state, which leads to the interaction switching. Due to electron correlation effects $δ$ deviates from the dynamical phase $ϕ=\int dt F(t)$, which enables the seemingly counterintuitive repulsion-to-attraction transition by a monocycle pulse with $ϕ=0$.

preprint2011arXiv

Dynamical band flipping in fermionic lattice systems: An ac-field-driven change of the interaction from repulsive to attractive

We show theoretically that the sudden application of an appropriate ac field to correlated lattice fermions flips the band structure and effectively switches the interaction from repulsive to attractive. The nonadiabatically driven system is characterized by a negative temperature with a population inversion. We numerically demonstrate the converted interaction in an ac-driven Hubbard model with the nonequilibrium dynamical mean-field theory solved by the continuous-time quantum Monte Carlo method. Based on this, we propose an efficient ramp-up protocol for ac fields that can suppress heating, which leads to an effectively attractive Hubbard model with a temperature below the superconducting transition temperature of the equilibrium system.

preprint2009arXiv

Photoinduced insulator-metal transition and nonlinear optical response of correlated electrons -- a DMFT analysis

We investigate a photoinduced insulator-metal transition in the Falicov-Kimball model with the dynamical mean-field theory (DMFT) extended to nonequilibrium systems in periodic modulations in time. When the photon energy of the pump light is $Ω\approx U$ ($U$: the interaction strength), a Drude-like peak is found to grow in the optical conductivity spectrum, which is an evidence that the system is driven into a metallic state. During the transition, the band gap does not collapse, whereas the distribution function exhibits a non-monotonic behaviour away from the Fermi distribution. This indicates that the transition cannot be accounted for by heating effects, but creation of photo-carriers is responsible.

preprint2009arXiv

Photoinduced insulator-metal transition in correlated electrons -- a Floquet analysis with the dynamical mean-field theory

In order to investigate photoinduced insulator-metal transitions observed in correlated electron systems, we propose a new theoretical method, where we combine a Floquet-matrix method for AC-driven systems with the dynamical mean-field theory. The method can treat nonequilibrium steady states exactly beyond the linear-response regime. We have applied the method to the Falicov-Kimball model coupled to AC electric fields, and numerically obtained the spectral function, the nonequilibrium distribution function and the current-voltage characteristic. The results show that intense AC fields indeed drive Mott-like insulating states into photoinduced metallic states in a nonlinear way.