## Abstract

We study steady states of semiconductor nanowires subjected to strong resonant time-periodic drives. The steady states arise from the balance between electron-phonon scattering, electron-hole recombination via photoemission, and Auger scattering processes. We show that tuning the strength of the driving field drives a transition between an electron-hole metal (EHM) phase and a Floquet insulator (FI) phase. We study the critical point controlling this transition. The EHM-to-FI transition can be observed by monitoring the presence of peaks in the density-density response function, which are associated with the Fermi momentum of the EHM phase and are absent in the FI phase. Our results may help guide future studies toward inducing exotic nonequilibrium phases of matter by periodic driving.

## INTRODUCTION

Coherent time-periodic driving provides a versatile tool for inducing novel properties in solid state and atomic systems (*1*–*12*). Prominent applications include Floquet engineering of band topology, light-induced superconductivity, and ultrafast spintronics (*13*–*29*). In many contexts, interesting phenomena may be observed in the short-time dynamics of driven systems (*30*–*32*). Under appropriate conditions, at long times, Floquet systems may also exhibit nontrivial steady-state characteristics such as topological responses or time-crystalline behavior (*33*–*47*).

Here, we investigate transitions between distinct phases realized in the steady states of a periodically driven semiconductor nanowire. We study the case where the drive frequency is larger than the band gap of the material. For weak (incoherent) driving, the steady state is established by a competition between excitation and relaxation rates that lead to a population inversion (i.e., the steady state hosts a nonequilibrium density of excited electrons at the bottom of the conduction band and holes at the top of the valence band, as indicated in Fig. 1A) (*48*, *49*).

When the drive is sufficiently strong, the conduction and valence bands are coherently hybridized for momenta near the “resonance momenta” ±*k*_{R} (see Fig. 1B). In this situation, the Floquet bands provide a good basis for identifying the physical properties and response characteristics of the many-body steady state. Here, we study such strongly driven systems, in the regime where the steady state is well described in terms of electronic populations of the system’s Floquet-Bloch states, with a nearly insulator-like filling in the Floquet basis. Such a system serves as a prototype for a Floquet topological insulator in which a resonant drive is used to induce an effective band inversion (*19*, *46*, *50*, *51*).

An ideal Floquet insulator (FI) is characterized by having a set of Floquet bands that are fully filled, while the remaining Floquet bands are empty. Away from the resonance momenta, the band inversion featured by the ideal FI implies that it hosts a nonequilibrium density of electrons and holes in the conduction and valence bands of the system’s equilibrium band structure. In this respect, the steady state bears similarity to that of a weakly driven system as shown in Fig. 1A. The natural relaxation of these electrons and holes through radiative recombination is manifested in the Floquet picture as interband transitions that create excitations away from the ideal FI state (making holes in the nominally filled Floquet band and putting electrons into the nominally empty Floquet band). Spontaneous emission therefore contributes a source of “quantum heating” in the Floquet basis (see Fig. 1B, wiggly arrow) (*52*). Similarly, inter-Floquet band transitions arising from electron-electron interactions may also create excitations away from the ideal FI state. At the same time, spontaneous electron-phonon scattering may lead to interband transitions that reduce the number of excitations, helping to “cool” the system toward the FI state (Fig. 1B, straight arrows). The steady state is determined by the competition between these scattering processes.

The steady states of nonequilibrium systems are generally sensitive to the rates of scattering processes, which, in turn, depend on a large number of microscopic parameters of the system and the baths to which it is coupled (in contrast to equilibrium where statistical mechanics provide a simple description of steady states in terms of a small number of thermodynamic parameters such as chemical potential and temperature). Nevertheless, we show that the electronic steady states of resonantly driven semiconductor nanowires exhibit two phases, distinguished by their properties for momenta near the resonance: (i) an electron-hole metal (EHM) phase, which features sharp Fermi surfaces for electron and hole excitations in the nominally empty and filled Floquet bands, respectively, and (ii) an FI phase, in which the electron and hole excitations are distributed as a nondegenerate Fermi gas in the Floquet basis. Which of these phases is obtained is determined by the dependence of the inter-Floquet band relaxation rates on the excitation density. A bottleneck for relaxation arises in the situation where these relaxation rates are suppressed at low densities, leading to a pileup of excitations and the formation of an EHM.

We show that a simple phase diagram emerges that is controlled by just two effective parameters, which capture the competition between the processes described above. The structure of the phase diagram is governed by a quantum critical point, with a critical region separating the two phases (see Fig. 1C). The transition between the EHM and the FI across the critical region is reminiscent of a finite-temperature crossover. The properties of the crossover are determined by the effective temperature of the electron and hole excitations in the Floquet bands.

As explained above, our main goal is to identify the important parameters that control the phase diagram and to characterize the critical point that separates the EHM and FI phases. To this end, we develop a simple analytical model based on a rate equation approach that captures the key mechanisms that determine the steady state of the driven system. For concreteness, we begin by providing a fully specified model for the driven system and the coupling to the environment. After obtaining the analytical results, we use numerical simulations of this microscopic model to verify the phase diagram and scaling behavior predicted by the rate equation approach. Last, we show that the phase diagram can be explored by tuning the driving frequency and strength. We further show that the EHM phase can be experimentally identified by observing peaks in the density-density response associated with the Fermi momentum of the excited electrons. This response gives rise to Friedel oscillations in the density induced by local inhomogeneities or an external potential.

## RESULTS

### Model of a periodically driven nanowire

To study the EHM and FI phases, we use a tight-binding model describing a uniformly, periodically driven nanowire with the Hamiltonian *k* along the wire. We write the single-particle Bloch Hamiltonian *H _{k}*(

*t*) in the form (

*53*–

*57*)

^{x}, σ

^{y}, and σ

^{z}are Pauli matrices in the orbital basis, with σ

^{z}∣ A〉 = ∣A〉 and σ

^{z}∣ℬ〉 = −∣ℬ〉, and

*A*,

*B*, and

*M*are constants. The periodic drive induces a local coupling between the orbitals, with strength

*V*. Throughout this work, we consider a half-filled system, which is a band insulator in the absence of the drive. The Floquet eigenstates of the time-dependent problem satisfy

*(*

_{kν}*t*)〉 =

*e*

^{−iεkνt/ħ}∣ϕ

_{k}_{ν}(

*t*)〉. Here,

_{k}_{ν}is the quasi-energy. Throughout, we use the convention −ħΩ/2 ≤ ε

_{k}_{ν}< ħΩ/2.

We study the regime where 2ħΩ is less than the total bandwidth (2∣*M*∣ + 2∣*B*∣). In this regime, the drive only resonantly couples states in the two bands via single-photon resonances; these resonances occur at crystal momentum values *k* = ±*k*_{R}, where ħΩ matches the splitting between valence and conduction bands of the nondriven system. The resulting Floquet spectrum exhibits a gap, Δ_{F}, separating the upper (ν = +) and lower (ν = −) Floquet bands. The form of the drive that appears in Eq. 1 is chosen such that the magnitude of the gap is equal to the driving field strength, up to corrections of order *V*^{2}/ħΩ. Other forms of the drive, such as a linearly polarized electric field, would have a different relation between the drive amplitude and Δ_{F} and would yield qualitatively similar results. A plot of a generic quasi-energy band structure for the Hamiltonian in Eq. 1 is shown in Fig. 1B.

### System-bath coupling

In addition to the coherent effects of the drive, captured in Eq. 1, we also describe the key dissipative processes that govern the steady states of the system. To this end, we incorporate in the model couplings between the electrons of the nanowire and acoustic phonons of the *d*-dimensional substrate upon which it sits (with *d* ≥ 2), as well as coupling of the electronic system to its (three-dimensional) electromagnetic environment. For simplicity, we consider a single-polarization mode for each bosonic bath (phonons or photons). The electromagnetic coupling allows radiative recombination of electron-hole pairs via spontaneous photon emission, which provides the primary source of heating in the Floquet basis. The effect of electron-electron interactions on the steady state is discussed in section S1.

The electron-boson coupling Hamiltonian *s* for “sound”) reads**q** = (*q*, **q**_{⊥}), and frequency ω = *v*_{λ} ∣**q**∣, where *v*_{λ} is the speed of light or sound. The Hamiltonian describing the heat baths reads **q**, denoted *q*, is the crystal momentum component parallel to the wire, and *q*_{⊥} represents its orthogonal component(s). (For photons, *q*_{⊥} has two components, while for phonons, *q*_{⊥} has one or two components, depending on whether *d* = 2 or *d* = 3.) The microscopic details of the electron-photon and electron-phonon couplings are captured by the functions *U*_{λ}(*q*, ω).

Respecting the particle-hole symmetry of the system, for small *q*, we take the coupling between electrons and acoustic phonons polarized along the wire to be (*57*, *58*)*g*_{s} is a coupling parameter, and *U*_{s} is proportional to the unity matrix in the 𝒜, B orbital space. For the electron-photon coupling, we take the simple (**q** independent) form *U*_{𝓁} = *g*_{𝓁}σ^{x}, where *g*_{𝓁} is a coupling parameter and σ^{x} is an orbital-space Pauli matrix (see text below Eq. 1). Our results can easily be extended to describe systems with other forms of electron-phonon coupling.

### Floquet kinetic equations

We work in the regime of weak system-bath coupling, where close to the steady state the electronic density matrix is well described in terms of populations of the Floquet eigenstates (*14*, *46*, *50*). Specifically, we will consider the regime ħ*W*_{scat} ≪ Δ_{F}, where *W*_{scat} is a typical scattering rate for an electron in the steady state (*50*). These populations are given by *e*^{−iεkν}* ^{t}*∣ϕ

_{k}_{ν}(

*t*)〉.

The system-bath coupling induces transitions between Floquet states. As a result, the populations {*f _{kν}*} evolve according to the kinetic equation

*k*in Floquet band ν to state

*k*

^{′}in Floquet band ν

^{′}. For a zero temperature bath, the intrinsic rate

*= ε*

_{l}

_{k}_{ν}− ε

_{k}_{′ν′}−

*l*ħΩ (see Eq. 16).

Because of the large speed of light, the momentum transfer *q* = *k* − *k*^{′} associated with spontaneous photon emission is small: ∣*q*∣ ≪ ω/*v*_{𝓁}. We thus take the density of states of the photon reservoir, integrated over perpendicular momentum components, as *v*_{𝓁}∣*q*∣ and ρ_{𝓁}(*q*, ω) = 0 otherwise. Here, *q* and ω.

Similarly, allowing substrates of different dimensionality, the phonon density of states reads*v*_{s}∣*q*∣< ω < ω_{D} and ρ_{s}(*q*, ω) = 0 otherwise. For acoustic phonons in a *d*-dimensional homogeneous substrate, η = *d* − 2. The value of η also depends on the dispersion of phonons. Here,

Note that ρ_{s}(*q*, ω) is cut off at the Debye frequency, ω_{D}. Throughout this work, take ω_{D} to be in the range Δ_{F} < ħω_{D} < Δ_{E}, where Δ_{E} is the gap at the quasi-energy zone edge ε = ħΩ/2 (see Fig. 1B). The condition ħω_{D} < Δ_{E} ensures the absence of phonon-mediated Floquet-Umklapp processes, corresponding to transitions with *l* > 0 (see text below Eq. 4) (*50*).

### Analysis of the steady state

We now analyze the steady state that emerges from the interplay between “excitation” processes associated with spontaneous emission of photons and “relaxation” processes mediated by the emission of acoustic phonons. From this analysis, we will extract the two key parameters that control the phase diagram of the system and characterize the critical point that separates the EHM and FI phases.

In the steady state, most of the excited electrons pile up in the two “valleys” near the minima of the upper Floquet band [see Fig. 1B and (*50*)]. The pileup results from a bottleneck in the excitation-relaxation cycle. This bottleneck is due to small intrinsic rates of inter-Floquet band relaxation for electrons close to the band minima compared with the rates of phonon-mediated intra-Floquet band relaxation of highly excited electrons within the upper Floquet band (see detailed discussion below). Because of particle-hole symmetry in our model, the holes form a mirror image population in the lower Floquet band.

Within the regime of interest, the distributions of electrons near the bottom of the upper Floquet band (*f*_{k+}) and holes near the top of the lower Floquet band (1 − *f*_{k−}) can, to a good degree, be approximately described by separate Fermi-Dirac distributions related by *f*_{k+} = 1 − *f*_{k−}, where*T* are the effective chemical potential counted from the bottom of the band and temperature of the electrons in the upper Floquet band (in energy units, *k*_{B} = 1, where *k*_{B} is Boltzmann’s constant). We note that although most of the excitations occupy the band minima according to Eq. 6, a small but finite density of excitations occupying high quasi-energy levels is necessary to maintain this distribution.

The EHM and FI phases are distinguished by whether μ lies above the bottom of the upper Floquet band, giving rise to a sharp Fermi surface for μ ≫ *T* (EHM phase), or below the band (i.e., in the Floquet gap), giving the FI phase. We thus seek to characterize the dependence of μ and *T* on the parameters of the system and the heat baths. To this end, we first develop a phenomenological model that captures the phase structure of the system and allows us to extract these dependencies. We will then corroborate these predictions with numerical simulations of the full kinetic equation (Eq. 4). In the main text, our analysis is done for zero bath temperature. The effects of finite bath temperature are discussed in section S2.

*Balance equation for the excitation density*

To build the phenomenological model, we seek two balance equations that must be satisfied by the populations of the Floquet bands in the steady state. The first balance equation fixes the value of the total excitation density of electrons in the upper Floquet band, *f _{k}* ≡

*f*

_{k+}. The equation for

*n*describes the balance between interband excitation and relaxation processes (

*48*,

*49*). A given value of

*n*corresponds to many different combinations of μ and

*T*. A second equation, following from the balance of intraband and interband relaxation processes, provides the relation between μ and

*T*.

We now analyze the processes leading to the steady state in more detail. Electron-hole excitations are predominantly generated by photon-mediated Floquet-Umklapp processes. In the low-excitation regime, *na* ≪ 1, which we consider throughout, photon-mediated processes excite electrons from an almost full to an almost empty Floquet band. Therefore, the excitation rate *k* = ± *k*_{R} and accumulate there. A mirror image process applies for the holes.

Annihilation of Floquet electron-hole pairs, which is primarily mediated by phonons, occurs on longer time scales than the relaxation to the band minima. The rate of such processes is proportional to the density of electron excitations in the upper Floquet band and hole excitations in the lower Floquet band, i.e., *50*) (here, we used that the electron and hole densities are equal for half filling). The rate equation for the density of excited electrons due to inter- and intra-Floquet band processes thus reads*59*).

*Balance equations for* μ *and T*

Next, we split the excitations in the upper band into two subpopulations of high and low quasi-energy excitations, which will allow us to characterize the form of the distribution (i.e., to determine μ and *T*). We partition the Brillouin zone into a momentum sector of width 2Δ*k* around the minimum of the upper band, *k* to be such that *k* with ε_{k}_{R + Δk, +} − ε_{k}_{R, +} ≪ *T* such that the distribution at *f*_{kR + Δk} can be approximated by *f*_{kR}. The partition of the momentum domain yields the partition of the density in the upper band, *n*, into two contributions: *n*_{R} and δ*n*, corresponding to the total density of electrons within *n* = *n* − *n*_{R} (see Fig. 2). The factor 2 in the definition of *n*_{R} accounts for the contributions of the two valleys.

The dominant source rate for *n*_{R} arises from the intraband scattering of electrons occupying quasi-energy states outside the region _{intra} captures the average intrinsic rate of collisions across the boundary of

To estimate Λ_{intra}, we need to estimate the energy and momentum transferred between the electrons and the phonons in this scattering process. We first note that as *f _{k}* decays over a quasi-energy window of order

*T*, we expect the densities δ

*n*and

*h*

_{R}to be mostly concentrated in this window. Therefore, we estimate the energy of the phonons emitted in the processes contributing to

*T*. The density of states for such phonons is nonvanishing for momentum transfers ∣

*q*∣ <

*T*/ħ

*v*

_{s}. Therefore, large momentum-transfer intraband processes connecting populations in the valleys near

*k*

_{R}and −

*k*

_{R}are forbidden for low-temperature steady states, which satisfy

*T*< 2

*k*

_{R}ħ

*v*

_{s}.

Because of the momentum dependence of *U*_{s}(*q*, ω) (Eqs. 2 and 3), the main contributions to Λ_{intra} come from the largest momentum transfers allowed within each valley, which we label *q*_{*}. Because at low densities *q*_{*} is small on the characteristic momentum scale of the Floquet dispersion, the overlap of the states before and after the transition is almost 1. Thus, we evaluate the electron-phonon coupling matrix element between these states by *q*_{*}*u*_{s}(ω) (see Eq. 3). Using Floquet-Fermi’s golden rule (see Eq. 16), we estimate

The source rate for *n*_{R} is balanced by annihilation of electrons comprising the population *n*_{R} with holes in the lower Floquet band, accompanied by phonon emission. These inter-Floquet band relaxation processes act as a sink for the density of excitations; the total rate of these depletion processes is proportional to the population and the total density of holes in the lower Floquet band, which is given by *n*. We thus estimate

The net interband scattering rate has two main contributions. The first arises from “vertical” (intravalley) processes (shown by a black arrow in Fig. 1B) and the second from “diagonal” (intervalley) processes (shown by an orange arrow). The former (intravalley) processes involve small momentum transfers, *q* ≲ *q*_{*}, while the latter processes require large momentum transfers, *q* ≈ 2*k*_{R}, where 2*k*_{R} is the separation in momentum of the two valleys. In both cases, the energy of the emitted phonon is approximately equal to the Floquet gap, ħω ≈ Δ_{F}. The intervalley processes are therefore kinematically constrained due to the large momentum transfers involved and only become active when the Floquet gap is above a critical value Δ_{F,c} ≡ 2ħ*v*_{s}*k*_{R}.

As a first step, we consider the case where the drive strength is low enough that the Floquet gap is below the critical value, Δ_{F} < Δ_{F,c}, where only the intravalley phonon-mediated inter-Floquet band scattering processes are active. To evaluate Λ_{inter}, we must estimate the electron-phonon matrix element governing these processes. We note that such processes connect nearly orthogonal electronic states (in different bands, at nearby momenta). Given the matrix form of the electron-phonon coupling in Eq. 3, for small *q*, the matrix element of the intravalley process is proportional to *q*^{2}*u*_{s}(Δ_{F}). In this estimate, we have used that the energy transfer is approximately equal to the Floquet gap, Δ_{F}.

Because of the fact that the matrix element above is suppressed for small *q*, the intravalley interband transition rate will be dominated by processes with the maximal momentum transfer, *q*_{*}. Thus, from Fermi’s golden rule, we obtain_{s}(Δ_{F}) is the phonon density of states (Eq. 5) evaluated at the energy transfer Δ_{F}. Here, we have defined α to highlight the total power of _{inter} in Eq. 10. Our model, Eq. 1, gives α = 2. Below, we will discuss the relation between *q*_{*} and the excitation density in each phase. Through this relation, we will see that α characterizes how the interband transition rate is suppressed at small excitation density. We thus refer to α as the “bottleneck parameter,” as it controls the ratio Λ_{inter}/Λ_{intra}.

Combining all of the source and sink rates for *n*_{R} found above, we arrive at the rate equation

Before turning to the analysis of the two phases, we note that for the specific model and regime discussed so far, the bottleneck parameter takes the value α = 2. However, the rate equation analysis applies more broadly, when α takes other values. In particular, we note that the initial and final electronic states of the intervalley interband scattering processes typically do not suffer from the near-orthogonality discussed for the intravalley processes. Moreover, for these processes, the electron-phonon coupling *U*_{s}(*q*, ω) is evaluated at a momentum transfer *q* ∼ 2*k*_{R}. Hence, the rates of intervalley interband scattering processes are governed by *q*_{*} appearing in

Below, we treat α and η as free parameters and map out a phase diagram (Fig. 1C) that captures the behavior across different parameter regimes (with a transition controlled by α − η). In particular, we will demonstrate the existence of a phase transition between the EHM and FI phases by obtaining critical exponents that characterize the phase transition in terms of α and η.

*The EHM phase*

We begin by considering the situation when interband relaxation is the rate-limiting step in the relaxation cascade of excited Floquet electron-hole pairs. This regime realizes the EHM phase, where the excitations on top of the FI state exhibit sharp Fermi surfaces, μ ≫ *T*. We will use Eqs. 8 to 11 to obtain the conditions on α and η where this regime is realized.

We now use Eq. 11 to obtain a relation between μ and *T* in the steady state (when _{inter} and Λ_{intra} by their estimates as functions of *q*_{*} and *T*, respectively (Eqs. 9 and 10). Under the assumption μ ≫ *T*, we have that *f _{k}* equals

*n*

_{R}, this yields Δ

*k*= π

*n*/2. We evaluate the densities

*n*≈

*n*

_{R}≈ 4μ

*D*(μ) and δ

*n*≈

*h*

_{R}≈ 2

*TD*(μ), where

*D*(ε) is the electronic density of states near the parabolic minimum of the upper Floquet band. Using this in Eq. 11 and solving for the steady state gives (2μ/

*T*)

^{2}= Λ

_{intra}/Λ

_{inter}. Next, we evaluate the maximal momentum transfer as

*q*

_{*}≈ π

*n*. Using

*D*

_{0}is a constant, and ρ

_{s}~ ω

^{η}(see Eq. 5), we obtain

*c*is a numerical constant of the order of unity.

By definition, the existence of an EHM phase with a sharp Fermi surface and a small density of excitations requires μ/*T* ≫ 1 and *T*/Δ_{F} ≪ 1. In order for the relation in Eq. 12 to be satisfied along with these conditions in the limit of small effective temperatures, the exponent (η − α)/(α + 1) must be negative; therefore, the EHM phase requires α < η. Again, using *n* ≈ 4μ*D*(μ) and expressing *n* in terms of κ as in Eq. 8, Eq. 12 yields relations for *T* and μ in terms of_{0} = (Δ_{F}*D*_{0})^{2}. Explicit expressions for ν_{T} and ν_{μ} appear in the table in Fig. 3A.

Another characteristic feature of the EHM phase is the dependence of the excitation density on the recombination rate Γ_{rec} (with all other parameters fixed). Taking into account the implicit dependence of Λ_{inter} on the excitation density through Eq. 10, from Eq. 8 we find

*The FI phase*

As explained below Eq. 12, for α < η, the EHM is not a consistent steady state of the rate Eqs. 7 and 11. We will now analyze the rate Eqs. 7 and 11 assuming that the distribution is a nondegenerate Fermi gas with μ < 0 (i.e., the chemical potential lies within the Floquet gap, below the band), which is characteristic of the FI phase. We will show that this distribution is a consistent solution of the rate Eqs. 7 and 11 for α < η.

We approximate the distribution as *f*_{k} ≈ *ze*^{–εk+/T}, where *z* = *e*^{-|μ|/T} is the fugacity. Our target here is to characterize the FI phase by solving for *z*, which gives the occupation at the bottom of the upper Floquet band and is related to the excitation density via *z* ≪ 1, we fix the value of Δ*k* in the definition of *k* ≪ *D*(*T*)*T* (see Fig. 2B). We note that other choices of Δ*k* will lead to the same power-law dependence. For this choice of Δ*k*, *f _{k}* is approximately constant within the region

*n*

_{R}≈ 4

*z*Δ

*k*,

*h*

_{R}≈ 4Δ

*k*, and δ

*n*=

*n*−

*n*

_{R}≈

*n*. Substituting these relations into Eq. 11 gives

*z*= Λ

_{intra}/Λ

_{inter}.

The spread of the density δ*n* over the momentum domain determines the typical momentum transfer, *q*_{*} ≈ *TD*(*T*). Given this value of *q*_{*} in the definitions of Λ_{intra} and Λ_{inter} (Eqs. 9 and 10), we solve for the steady state of Eq. 11, *c*′ is a numerical coefficient of order unity, generically different from *c* in Eq. 12.

Combining Eq. 14 with the expressions *n* ≈ κ^{1/2} yields a dependence of *z* and *T* on κ for the FI phase in the form of power laws with exponents ν_{T} and ν_{z}

Another characteristic feature of the EHM phase is the dependence of the excitation density on the recombination rate Γ_{rec} (with all other parameters fixed). Taking into account the implicit dependence of Λ_{inter} on the excitation density through Eq. 10, from Eq. 8 we find

Note the difference between the exponents ν_{T} in the FI and EHM phases. This difference, together with the power laws for *z* in the FI phase and μ in the EHM phase, indicate discontinuous dependencies of these quantities on κ (at small κ). This nonanalytic behavior implies the existence of a quantum critical point at α = η when κ → 0. Increasing κ increases the effective temperature of the steady state, *T*, giving rise to a finite effective temperature crossover above the critical point in the κ − α plane. The phase diagram (in the κ − α plane) characterizing the steady state of the system is depicted in Fig. 1C.

### Numerical results

To support the analysis above, we numerically solve for the steady state of the full Floquet kinetic equation (Eq. 4) on a lattice of 5000 *k* points at half filling. Because the effective value of α depends in a complicated way on the system and drive parameters, in the numerical simulations, we verified the critical behavior near the EHM-to-FI phase transition by fixing α and varying η. We fixed α = 2 by setting Δ_{F} < Δ_{F,c} and swept the value of η (see Eq. 5) from η = 0 to η = 4 across the critical value at η = α = 2. For each fixed set of parameters, we fit the steady-state distribution of electrons in the upper Floquet band to a Fermi-Dirac distribution (Eq. 6) to extract *T*, μ, and *z*. We repeat this procedure for several values of *g*_{𝓁} and a fixed *g*_{s}, yielding values of κ over several decades, and extract the power-law scalings of these three quantities as functions of κ. The resulting numerically extracted exponents ν_{T}, ν_{μ}, and ν_{z} are plotted as functions of η in Fig. 3B. We find good agreement between our analytical and numerical results.

### Inducing and observing the phase transition

In this section, we discuss how to induce and observe the transition between the EHM and FI phase in a system with fixed band structure and bath parameters. By tuning the driving amplitude, which controls the Floquet gap Δ_{F}, the kinematic constraints that control the relative contributions of intravalley and intervalley interband scattering processes can be modified in situ. So far, we discussed the regime Δ_{F} < Δ_{F,c}, in which only the intravalley processes, corresponding to α = 2, are active (see the discussion above Eq. 10). For Δ_{F} > Δ_{F,c}, however, the kinematic constraint on the phonons that forbids the intervalley processes is lifted. In this case, the interband relaxation rate consists of a sum of the intravalley and intervalley processes (see Fig. 1B).

As discussed above, the intervalley interband relaxation processes are generically characterized by α = 0. In the presence of both processes, we expect an intermediate behavior between that described by the power laws in Eqs. 13 and 15 with α = 2 and α = 0. However, for small values of κ, processes with the lowest value of α dominate. Therefore, for η = 1, corresponding to phonons in three dimensions, the system should undergo a sharp transition from the EHM to the FI phase as Δ_{F} is varied across Δ_{F,c} by changing the drive amplitude. This sharp transition is depicted schematically by the dark red arrow in Fig. 1C (note that increasing Δ_{F,c} also decreases the value of κ, as it increases the density of states for phonon-assisted interband relaxation).

The physical difference between the EHM and FI phases is manifested in the response of the system to density perturbations. For demonstration, we numerically compute the density response to potential charge fluctuations, described by χ(*q*, 𝜔) (see Eq. 18). We examine the behavior of χ(*q*, 𝜔) as a function of driving amplitude in the vicinity of the EHM-to-FI transition. For details of the numerical simulations, see section S3. The inset to Fig. 4A shows the real part of the response function χ′(*q*) = Re[χ(*q*, 0)]. In the FI phase, χ(*q*) exhibits large peaks at zero wave number and at wave numbers connecting the two valleys at *k _{R}* and –

*k*. Because of the sharp Fermi surface in the EHM phase, each peak splits into two peaks separated by 2

_{R}*k*

_{F}≡ π

*n*. In particular, χ′(

*q*) exhibits two peaks at

*q*= ± 2

*k*

_{F}. The splitting of these peaks is absent in the FI phase (see Fig. 4B and C). Experimental signatures of the split peaks are Friedel oscillations in the screening potential, which can be observed, for example, using x-ray scattering (

*60*) and scanning tunneling microscopy (

*61*–

*63*) or through a minimally invasive electrostatic probe as in (

*64*). In Fig. 4A, we plot the “strength” of the 2

*k*

_{F}peaks, defined as

*V*(i.e., along the trajectory depicted by the dark red arrow in Fig. 1C). The peaks disappear abruptly when the driving field imposes Δ

_{F}> Δ

_{F,c}, for which the steady state is in the FI phase.

## DISCUSSION

To appreciate the physical scales, we estimate the value of κ in periodically driven semiconductors in the FI and EHM phases. To this end, we use a typical value for the radiative recombination lifetime, τ_{rr} ∼ 1 to 10 ns, and use the hot-electron lifetime as a typical time scale for phonon-mediated relaxation, τ_{he} ∼ 10 to 100 fs, respectively (*65*–*69*). We further take *k*_{R} ∼ 0.01π/*a*. In the FI phase, this yields the estimate _{0} ≡ (Δ_{F}*D*_{0})^{2} . We estimate *v*_{R} is the band velocity of the semiconductor at the resonance momentum, *k*_{R}. In the EHM phase, the highest possible κ_{0} is obtained for Δ_{F} ≈ Δ_{F,c} ≡ 2ħ*k*_{R}*v*_{s} (see discussion in the second paragraph below Eq. 9). We thus estimate κ_{0}*a*^{2} ≈ 2(*m*_{*}*v*_{s}*a*/πħ)^{2}, where *m*_{*} ≡ *ħk*_{R}/*v*_{R} is the effective mass of the semiconductor. Therefore, large sound velocities and effective masses are favorable for realization of the EHM phase.

In this work, we uncovered a transition between EHM and FI phases in a driven one-dimensional electronic system. Our results can be generalized to other one- and two-dimensional resonantly driven Floquet-Bloch systems. In two-dimensional systems, we expect that the system still supports EHM and FI phases that can be accessed via the drive strength. Studying the features of the transition in two-dimensional systems and possibilities for controlling their behavior by reservoir engineering is an important direction for future study.

## MATERIALS AND METHODS

### Evaluation of the intrinsic scattering rates

We have evaluated the intrinsic scattering rates, * _{l}* ≡ ε

_{kν}− ε

_{k′ν′}−

*l*ħΩ. For the full derivation, see section S4.

### Floquet-Lindhard function

The density response to potential charge fluctuations, χ(*x*, *t*), is evaluated by the Floquet-Lindhard function. Following the Kubo formula, we define χ(*x*, *t*) through the density-density response function averaged over the driving period, 𝒯 ≡ 2π/Ω*t*) = 1 for *t* > 0 and Θ(*t*) = 0; otherwise, *70*)*71*).

## SUPPLEMENTARY MATERIALS

Supplementary material for this article is available at http://advances.sciencemag.org/cgi/content/full/6/35/eaay4922/DC1

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is **not** for commercial advantage and provided the original work is properly cited.

## REFERENCES AND NOTES

**Acknowledgments:**We would like to thank T. Gulden, G. K. Gupta, V. Kalnitsky, G. Matsman, and A. Turner for discussions and D. Cohen and Y. Katz for technical support.

**Funding:**N.H.L. acknowledges support from the European Research Council (ERC) under the European Union Horizon 2020 Research and Innovation Programme (grant agreement no. 639172) and from the Israeli Center of Research Excellence (I-CORE) “Circle of Light.” M.S.R. acknowledges the support of the European Research Council (ERC) under the European Union Horizon 2020 Research and Innovation Programme (grant agreement no. 678862) and the Villum Foundation. I.E. acknowledges support from the Ministry of Science and Technology, Israel.

**Author contributions:**I.E. performed the numerical and analytic calculations. All authors contributed to the data analysis and wrote the manuscript.

**Competing interests:**The authors declare that they have no competing interests.

**Data and materials availability:**All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors.

- Copyright © 2020 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution NonCommercial License 4.0 (CC BY-NC).