A tremendous amount of recent attention has focused on characterizing the dynamical properties of periodically driven many-body systems. Here, we use a novel numerical tool termed "density matrix truncation" (DMT) to investigate the late-time dynamics of large-scale Floquet systems. We find that DMT accurately captures two essential pieces of Floquet physics, namely, prethermalization and late-time heating to infinite temperature. Moreover, by implementing a spatially inhomogeneous drive, we demonstrate that an interplay between Floquet heating and diffusive transport is crucial to understanding the system's dynamics. Finally, we show that DMT also provides a powerful method for quantitatively capturing the emergence of hydrodynamics in static (undriven) Hamiltonians; in particular, by simulating the dynamics of generic, large-scale quantum spin chains (up to L=100), we are able to directly extract the energy diffusion coefficient.
A tremendous amount of recent attention has focused on characterizing the dynamical properties of periodically driven many-body systems. Here, we use a novel numerical tool termed "density matrix truncation" (DMT) to investigate the late-time dynamics of large-scale Floquet systems. We find that DMT accurately captures two essential pieces of Floquet physics, namely, prethermalization and late-time heating to infinite temperature. Moreover, by implementing a spatially inhomogeneous drive, we demonstrate that an interplay between Floquet heating and diffusive transport is crucial to understanding the system's dynamics. Finally, we show that DMT also provides a powerful method for quantitatively capturing the emergence of hydrodynamics in static (undriven) Hamiltonians; in particular, by simulating the dynamics of generic, large-scale quantum spin chains (up to L=100), we are able to directly extract the energy diffusion coefficient.
Extracting late-time hydrodynamic transport properties directly from microscopic quantum models, especially in dimensions beyond one, is believed to be classically intractable. This owes, in part, to the extensive growth of entanglement in higher dimensions and often a mismatch between geometric and tensor network locality. Here, we extend the Density Matrix Truncation (DMT) algorithm to accurately and quantitatively predict high-temperature transport coefficients (e.g. diffusion constants) in up to three dimensions. Additionally, we demonstrate the ability to capture the slow crossover dynamics characteristic of weakly broken integrability. The intuition underlying our approach is to preserve the geometrically local information crucial for hydrodynamics throughout the time evolution, despite using a one-dimensional matrix product state ansatz for higher-dimensional calculations. Our work pushes the boundary of classical computability and opens the door to benchmarking the high-temperature, late-time dynamics of generic quantum simulators.
In ergodic quantum spin chains, locally conserved quantities such as energy or particle number generically evolve according to hydrodynamic equations as they relax to equilibrium. We investigate the complexity of simulating hydrodynamics at infinite temperature with multiple methods: time evolving block decimation (TEBD), TEBD with density matrix truncation (DMT), the recursion method with a universal operator growth hypothesis, and operator-size truncated (OST) dynamics. Density matrix truncation and the OST dynamics give consistent dynamical correlations to t=60/J and diffusion coefficients agreeing within 1. TEBD only converges for $t\ensuremath{\lesssim}20$, but still produces diffusion coefficients accurate within 1%. The universal operator growth hypothesis fails to converge and only matches other methods on short times. We see no evidence of long-time tails in either DMT or OST calculations of the current-current correlator, although we cannot rule out that they appear at longer times. We do observe power-law corrections to the energy density correlator. At finite wavelength, we observe a crossover from purely diffusive, overdamped decay of the energy density to underdamped oscillatory behavior similar to that of cold-atom experiments. We dub this behavior ``hot band second sound'' and offer a microscopically motivated toy model.
Prethermalization, by introducing emergent quasiconserved observables, plays a crucial role in protecting periodically driven (Floquet) many-body phases over an exponentially long time, while the ultimate fate of such quasiconserved operators can signal thermalization to infinite temperature. To elucidate the properties of prethermal quasiconservation in many-body Floquet systems, here we systematically analyze infinite-temperature correlations between observables. We numerically show that the late-time behavior of the autocorrelations unambiguously distinguishes quasiconserved observables from nonconserved ones, allowing one to single out a set of linearly independent quasiconserved observables. By investigating two Floquet spin models, we identify two different mechanisms underlying the quasiconservation law. First, we numerically verify energy quasiconservation when the driving frequency is large, so that the system dynamics is approximately described by a static prethermal Hamiltonian. More interestingly, under moderate driving frequency, another quasiconserved observable can still persist if the Floquet driving contains a large global rotation. We show theoretically how to calculate this conserved observable and provide numerical verification. Having systematically identified all quasiconserved observables, we can finally investigate their behavior in the infinite-time limit and thermodynamic limit, using autocorrelations obtained from both numerical simulation and experiments in solid-state nuclear magnetic resonance systems.
Periodically driven Floquet quantum many-body systems have revealed new insights into the rich interplay of thermalization and growth of entanglement. The phenomenology of , whereby a translationally invariant many-body system exhibits emergent conservation laws and a slow growth of entanglement entropy at certain fixed ratios of a drive amplitude and frequency, presents a novel paradigm for retaining memory of an initial state up to late times. Previous studies of dynamical freezing have largely been restricted to a high-frequency Floquet-Magnus expansion and numerical exact diagonalization. Both techniques are unable to capture the slow approach to thermalization, or lack thereof, in a systematic fashion. By employing Floquet , where the time-dependent part of the Hamiltonian is gradually decoupled from the effective Hamiltonian using a sequence of unitary transformations, we unveil the universal approach to dynamical freezing and beyond, at asymptotically late times. We analyze the behavior associated with the flow renormalization at and near freezing using both exact-diagonalization and tensor-network-based methods and contrast the results with the conventional prethermal phenomenon. For a generic nonintegrable spin Hamiltonian with a periodic cosine wave drive, the flow approaches an unstable fixed point with an approximate emergent symmetry. We observe that at freezing the thermalization timescales are delayed compared to away from freezing, and the flow trajectory undergoes a series of events. Our numerical results are supported by analytical solutions to the flow equations.
We investigate a class of periodically driven many-body systems that allows us to extend the phenomenon of prethermalization to the vicinity of isolated intermediate-to-low drive frequencies away from the high-frequency limit. We provide numerical evidence for the formation of a parametrically long-lived prethermal plateau, captured by an effective Floquet Hamiltonian computed using the replica inverse-frequency expansion, and demonstrate its stability w.r.t.~random perturbations in the drive period. Considering exclusively nonintegrable Floquet Hamiltonians, we find that heating rates are non-universal: we observe Fermi's Golden Rule scaling, power-law scaling inconsistent with the Golden Rule, and non-power law scaling, depending on the drive. Despite the asymptotic character of the inverse-frequency expansion, we show that it describes the thermostatic properties of the state all along the evolution up to infinite temperature, with higher-order terms improving the accuracy. Our results suggest a dynamical mechanism to gradually increase temperature in isolated quantum simulators, such as ultracold atoms, and open up an alternative possibility to investigate thermal phase transitions and the interplay between thermal and quantum criticality using Floquet drives.
Periodically-driven (Floquet) quantum systems are ubiquitous in science and technology. For example, when a laser illuminates a material or an AC voltage is applied to a device, the system is well-described by a time-periodic Hamiltonian. In recent years, periodic driving has been proposed, not just as a tool to excite and probe devices, but actually as a mechanism of engineering new phases of matter, some of which have no equilibrium analog. However, with this promise comes a serious problem. Intuitively, if energy is injected into and distributed throughout a system, it is no surprise that it tends to heat up indefinitely to infinite temperature. In this thesis, we study the mechanisms of heating, i.e. the process of thermalization, in Floquet systems and propose methods to control them. Specifically, for non-interacting Floquet systems that are coupled to external bosonic and fermionic baths (e.g. laser-driven electrons in a semiconductor that interact with phonons and an external lead), we classify the relevant scattering processes that contribute to cooling/heating in the Floquet bands and suggest methods to suppress heating via bandwidth-restrictions on the baths. We find that is possible, with appropriate dissipative engineering, to stabilize a controlled incompressible nonequilibrium steady-state resembling a ground state - a state we term the "Floquet insulator." We extend this analysis to include short-range interactions that contribute additional heating processes and show, under the same framework, that heating can be controlled with dissipation. In the process, we develop a simple effective model for the Floquet band densities that captures the essence of all the Floquet scattering processes and that is useful for ballparking experimentally-relevant estimates of heating. Next, we turn our attention to strongly-interacting closed Floquet systems and study how heating emerges through a proliferation of resonances. We find a novel integrable point governing the strong-interaction limit of the Floquet system and examine the breakdown of integrability via the proliferation of resonances. We observe two distinct scaling regimes, attributed to non-thermal and thermal behavior, and discover a power-law scaling of the crossover between them as a function of system size. The lingering ergodicity-breaking effects of the conserved quantities in the vicinity (in parameter space) of the integrable point at finite size is a phenomena we term "near-integrability." These results suggest that small quantum systems, which are accessible currently in many platforms (e.g. trapped ions, cold atoms, superconducting devices), intrinsically host non-thermal states that one may be able to utilize to avoid heating. Furthermore, our results suggest a "dual" interpretation, in the thermodynamic limit, that a periodically-driven system exhibits prethermalization as a power-law in interaction strength.