Floquet Replica Space

• Floquet Replica Space is an enlarged Hilbert space construction that combines physical states with periodic time harmonics to convert time-periodic dynamics into a static eigenvalue problem.
• It enables analysis of quasienergy folding, symmetry actions, and inter-replica coupling, providing insights into the structure of driven quantum systems.
• Its applications span driven quantum circuits, space-time crystals, and dissipative systems, offering a framework for exploring topological phases and operator growth.

Floquet replica space is the enlarged Hilbert-space construction used to represent a time-periodic problem as a static eigenvalue problem. In its standard form, it is the tensor product of the physical Hilbert space with the space of periodic time harmonics, so that each integer harmonic labels a “replica” or Floquet sector shifted by a drive quantum. Across the literature, closely related constructions are called Sambe space, extended Hilbert space, Floquet-Hilbert space, or, in specialized settings, Floquet phase space and symplectic-Floquet space. In this language, quasienergies are defined modulo the drive frequency, periodic driving appears as inter-replica coupling, and questions of symmetry, topology, transport, chaos, dissipation, and operator growth are reformulated as structure in an extended space (Jin et al., 2021, Qin et al., 2017, Braver et al., 2024).

1. Formal definition and basic structure

For a Hamiltonian with period $T$, $H(t+T)=H(t)$, the one-period evolution operator is

$$U_F=\mathcal{T}\exp\!\left[-i\int_0^T H(t)\,dt\right],$$

and quasienergies are defined modulo $\Omega=2\pi/T$. In the standard extended-space construction one introduces a basis $|\alpha,m\rangle$, or equivalently $|j,m\rangle=|j\rangle e^{im\Omega t}$, where $|\alpha\rangle$ or $|j\rangle$ is a physical basis state and $m\in\mathbb Z$ is the Floquet harmonic index. The Floquet operator in this space is

$\mathcal{H}_F = H(t)-i\partial_t,$

with matrix elements

$H(t+T)=H(t)$0

Each replica is therefore a copy of the physical Hilbert space shifted by $H(t+T)=H(t)$1, and the time-dependent problem becomes a static eigenvalue problem in an infinite matrix indexed by $H(t+T)=H(t)$2 (Jin et al., 2021, Qin et al., 2017).

The same structure appears in Green’s-function formulations. A two-time Green’s function can be rewritten as a Floquet matrix $H(t+T)=H(t)$3, with $H(t+T)=H(t)$4 restricted to the first Floquet Brillouin zone and the replica indices carrying the harmonic content. In this representation, higher Floquet bands are simply components with larger $H(t+T)=H(t)$5, and photon-assisted processes appear as off-diagonal matrix elements between distinct replicas (Qin et al., 2017).

A closely related notation is the Floquet-Hilbert-space basis $H(t+T)=H(t)$6, in which the quasienergy operator $H(t+T)=H(t)$7 has diagonal blocks $H(t+T)=H(t)$8 and off-diagonal blocks given by the Fourier components $H(t+T)=H(t)$9. This makes explicit that quasienergy spectra are ladders $$U_F=\mathcal{T}\exp\!\left[-i\int_0^T H(t)\,dt\right],$$0, all physically equivalent modulo $$U_F=\mathcal{T}\exp\!\left[-i\int_0^T H(t)\,dt\right],$$1 (Braver et al., 2024).

2. Replica space beyond the textbook Sambe construction

In discrete-time systems the same idea survives, but the “replicas” need not originate from Fourier harmonics alone. A notable example is the class of same-gate quantum circuits built from translationally invariant nearest-neighbor two-site gates. For periodic boundary conditions, every simple circuit is spectrally equivalent to a canonical propagator $$U_F=\mathcal{T}\exp\!\left[-i\int_0^T H(t)\,dt\right],$$2, and in a Floquet-theoretic interpretation one writes

$$U_F=\mathcal{T}\exp\!\left[-i\int_0^T H(t)\,dt\right],$$3

Here the physical Floquet period is built from $$U_F=\mathcal{T}\exp\!\left[-i\int_0^T H(t)\,dt\right],$$4 applications of a root operator $$U_F=\mathcal{T}\exp\!\left[-i\int_0^T H(t)\,dt\right],$$5, followed by a translation. Within each eigenspace of $$U_F=\mathcal{T}\exp\!\left[-i\int_0^T H(t)\,dt\right],$$6, the spectrum of $$U_F=\mathcal{T}\exp\!\left[-i\int_0^T H(t)\,dt\right],$$7 is a folded $$U_F=\mathcal{T}\exp\!\left[-i\int_0^T H(t)\,dt\right],$$8-fold image of the spectrum of $$U_F=\mathcal{T}\exp\!\left[-i\int_0^T H(t)\,dt\right],$$9. This suggests a replica picture in which the layer index $\Omega=2\pi/T$0 acts as a discrete Floquet-harmonic label and the physical spectrum is a superposition of $\Omega=2\pi/T$1 symmetry-related quasi-energy sectors (Duh et al., 2024).

A broader generalization appears in space-time crystals. There the fundamental reciprocal lattice is not the rectangular $\Omega=2\pi/T$2 lattice of conventional Floquet-Bloch theory, but an oblique space-time lattice generated by mixed reciprocal vectors. The corresponding space-time Floquet operator evolves the system over a fraction of the physical period,

$\Omega=2\pi/T$3

Its eigenvalues define a space-time band structure that unfolds ordinary Floquet bands. In this formulation, conventional Floquet replicas are a folded projection of a more fundamental mixed momentum-frequency replica lattice (Melkani et al., 18 Oct 2025).

A different, but mathematically parallel, reinterpretation appears in “Floquet phase space.” For a periodically driven particle near an $\Omega=2\pi/T$4 resonance, a secular approximation yields an effective Hamiltonian

$\Omega=2\pi/T$5

with $\Omega=2\pi/T$6 a slow angle variable. This suggests a synthetic lattice in the angular coordinate, and the resulting Bloch-like bands in $\Omega=2\pi/T$7 are the phase-space analogue of a replica ladder (Zhang et al., 2021).

3. Symmetry actions in replica space

Replica space is often the natural arena in which space-time symmetries become algebraically transparent. In driven nonlinear photonic crystals with a space-time screw symmetry $\Omega=2\pi/T$8, the action on a Floquet basis state is

$\Omega=2\pi/T$9

where $|\alpha,m\rangle$0 is the static $|\alpha,m\rangle$1 eigenvalue. The factor $|\alpha,m\rangle$2 comes entirely from the time translation by $|\alpha,m\rangle$3. Replica index $|\alpha,m\rangle$4 therefore directly changes the symmetry eigenvalue, and only replica states with the same $|\alpha,m\rangle$5 index can hybridize. In that setting, symmetry indicators and nested Wannier bands are defined on Floquet bands in the enlarged space, not on the undriven bands (Jin et al., 2021).

For time-glide-symmetric Floquet systems, the enlarged frequency-domain Hamiltonian

$|\alpha,m\rangle$6

inherits an ordinary reflection symmetry in replica space. In a two-replica truncation near quasienergy $|\alpha,m\rangle$7, the effective reflection operator is

$|\alpha,m\rangle$8

with $|\alpha,m\rangle$9 acting on the replica pair. This mapping allows a time-glide Floquet problem to be analyzed as a static reflection-symmetric higher-order topological insulator in enlarged space (Peng et al., 2018).

A related prethermal construction shows that dynamical space-time symmetries can map onto the projective static symmetry group of a prethermal Hamiltonian. In that framework the micromotion operator and the effective Hamiltonian transform covariantly under an extended group that includes both a finite-order Floquet unitary and order-two dynamical symmetries, so symmetry constraints become statements about the structure of the interaction-picture replica space rather than only about stroboscopic evolution (Na et al., 2024).

4. Spectral folding, roots, and effective Hamiltonians

One recurring use of replica-space reasoning is to distinguish a physically fundamental spectrum from its folded Floquet image. In same-gate quantum circuits, if

$|j,m\rangle=|j\rangle e^{im\Omega t}$0

then in a momentum sector of $|j,m\rangle=|j\rangle e^{im\Omega t}$1 the eigenphases of the full propagator satisfy

$|j,m\rangle=|j\rangle e^{im\Omega t}$2

The map $|j,m\rangle=|j\rangle e^{im\Omega t}$3 folds the root spectrum onto the unit circle. As $|j,m\rangle=|j\rangle e^{im\Omega t}$4 increases, this folding destroys level correlations, so the raw Floquet spectrum can look increasingly Poisson-like even when the root operator has circular-unitary-ensemble statistics. In this setting, the appropriate desymmetrization is not a conventional symmetry-sector decomposition of $|j,m\rangle=|j\rangle e^{im\Omega t}$5, but analysis of the root $|j,m\rangle=|j\rangle e^{im\Omega t}$6 itself (Duh et al., 2024).

A different but related spectral problem arises in resonantly driven interacting systems, where standard high-frequency expansions fail because resonant processes make denominators small. Extended degenerate perturbation theory addresses this directly in Floquet-Hilbert space by enlarging the “degenerate” subspace to an entire Floquet zone. The reordered quasienergy matrix $|j,m\rangle=|j\rangle e^{im\Omega t}$7 is built so that each $|j,m\rangle=|j\rangle e^{im\Omega t}$8-th diagonal block contains energies reduced to the $|j,m\rangle=|j\rangle e^{im\Omega t}$9-th Floquet zone, and the perturbation theory then resembles a van Vleck expansion while retaining the exact intra-zone structure. The resulting effective Hamiltonian is more accurate than conventional DPT while remaining less costly than exact Floquet diagonalization (Braver et al., 2024).

A common misconception is that the one-period propagator is always the correct object for spectral diagnosis. The circuit classification results show that hidden space-time symmetry can make the full propagator a folded power of a more primitive root, so physically relevant level statistics may only emerge after resolving the appropriate replica or root structure (Duh et al., 2024).

5. Topology, transport, and experimental signatures

Replica space is also where Floquet topology is most naturally defined. In driven quadrupole photonic crystals, the pair of Floquet bands below the drive-induced gap carries a quadrupole moment determined by the space-time screw indices at high-symmetry momenta,

$|\alpha\rangle$0

and the same phase is confirmed by nested Wannier bands. The resulting Floquet quadrupole phase has $|\alpha\rangle$1 and supports corner states inside the Floquet gap (Jin et al., 2021).

In a 1D tight-binding chain under a uniform electric field, the extended structure can be reorganized as a momentum-frequency Brillouin zone. Electric translation operators shift momentum and frequency in a projective fashion, and the corresponding electric Floquet-Bloch wavefunctions are labeled by $|\alpha\rangle$2 rather than by a momentum plus an unbounded harmonic index. In this formulation, ordinary Floquet replicas become bands on a compact momentum-frequency torus, and the Zak phase winds by $|\alpha\rangle$3 around the frequency cycle (Ke et al., 2024).

For space-time crystals, the topology of the unfolded space-time bands is encoded by winding numbers $|\alpha\rangle$4 and $|\alpha\rangle$5. These determine quantized transport over one fundamental space-time step,

$|\alpha\rangle$6

while the transport over one full Floquet period is the integer

$|\alpha\rangle$7

This is the paper’s fractional version of adiabatic charge transport: the full-period response is integer, but the natural space-time step can carry a rational topological charge (Melkani et al., 18 Oct 2025).

Experimental access to replica structure is particularly clear in graphene. Time- and angle-resolved photoemission has directly observed replicas of the Dirac cone in monolayer epitaxial graphene, and the polarization dependence identifies their origin through scattering between Floquet-Bloch states and Volkov states (Choi et al., 2024). A complementary driven-dissipative analysis shows that in a low-frequency, high-amplitude regime the Dirac-point gap

$|\alpha\rangle$8

can exceed the Floquet-zone width, so the dominant occupied Floquet replicas are separated by more than $|\alpha\rangle$9. This allows the electron distribution across replicas to distinguish Floquet replicas from laser-assisted photoemission replicas in realistic trARPES conditions (Broers et al., 2021).

6. Interacting, dissipative, and operator-space extensions

The replica construction remains useful when the physical Hilbert space fragments or when the dynamics are open. In a periodically kicked XXZ spin chain, exact conservation of the absolute magnetization

$|j\rangle$0

and approximate conservation of the domain-wall number $|j\rangle$1 produce a fine block structure of the Floquet operator. In replica-space language, each Floquet replica splits into many fragments $|j\rangle$2. Strong fragmentation localizes Floquet eigenstates within these blocks, and $|j\rangle$3-pairs inside small fragments stabilize discrete time-crystalline order without disorder (Tang et al., 16 Dec 2025).

For interacting driven fermions, real-space Floquet DMFT promotes Green’s functions, self-energies, and Dyson equations to matrices in replica space. Keeping only the central replica reproduces a high-frequency effective Hamiltonian, whereas keeping several replicas captures higher Floquet bands, sidebands, and non-equilibrium steady-state occupations. This clarifies the difference between a driven steady state and an equilibrium state of an effective static model (Qin et al., 2017).

Open-system extensions may require non-Hermitian generators. In the dynamical Casimir problem, the relevant enlarged space is a symplectic-Floquet space $|j\rangle$4, where the physical degrees of freedom are canonical mode operators rather than wavefunctions. The resulting Floquet-Liouvillian is non-Hermitian, its eigenmodes are multimode Bogoliubov modes, and the competition between parametric amplification and dissipation is encoded in complex quasifrequencies in the replica-resolved spectrum (Tanaka et al., 2020).

Finally, there are operator-space analogues of Floquet replica space. The Floquet operator Krylov construction treats the adjoint action $|j\rangle$5 as motion on an emergent one-dimensional chain indexed by Krylov depth rather than by harmonic number. This maps stroboscopic operator dynamics onto a Floquet inhomogeneous transverse-field Ising model whose couplings are Krylov angles; the construction is a different extended-space representation, but it plays a role parallel to replica space by converting a driven many-body problem into structured motion on a synthetic dimension (Yeh et al., 2024, Yeh et al., 2023).

Across these settings, Floquet replica space is not merely a bookkeeping device for Fourier harmonics. It is the framework in which hidden roots, mixed space-time symmetries, topology, fragmentation, dissipation, and even operator growth become manifest.