Reply to Comment on Discrete Time Crystals: Rigidity Criticality and Realizations

This is a reply to the comment from Khemani, Moessner and Sondhi (KMS) [arXiv:2109.00551] on our manuscript [Phys. Rev. Lett. 118, 030401 (2017)]. The main new claim in KMS is that the short-ranged model does not support an MBL DTC phase. We show that, even for the parameter values they consider and the system sizes they study, the claim is an artifact of an unusual choice of range for the crucial plots. Conducting a standard finite-size scaling analysis on the same data strongly suggests that the system is in fact a many-body localized (MBL) discrete time crystal (DTC). Furthermore, we have carried out additional simulations at larger scales, and provide an analytic argument, which fully support the conclusions of our original paper. We also show that the effect of boundary conditions, described as essential by KMS, is exactly what one would expect, with boundary effects decreasing with increasing system size. The other points in KMS are either a rehashing of points already in the literature (for the long-ranged model) or are refuted by a proper finite-size scaling analysis.

The Variational Power of Quantum Circuit Tensor Networks

We characterize the variational power of quantum circuit tensor networks in the representation of physical many-body ground-states. Such tensor networks are formed by replacing the dense block unitaries and isometries in standard tensor networks by local quantum circuits. We explore both quantum circuit matrix product states and the quantum circuit multi-scale entanglement renormalization ansatz, and introduce an adaptive method to optimize the resulting circuits to high fidelity with more than $10^4$ parameters. We benchmark their expressiveness against standard tensor networks, as well as other common circuit architectures, for both the energy and correlation functions of the 1D Heisenberg and Fermi-Hubbard models in the gapless regime. We find quantum circuit tensor networks to be substantially more expressive than other quantum circuits for these problems, and that they can even be more compact than standard tensor networks. Extrapolating to circuit depths which can no longer be emulated classically, this suggests a region of quantum advantage with respect to expressiveness in the representation of physical ground-states.

Superuniversality from disorder at two-dimensional topological phase transitions

We investigate the effects of quenched randomness on topological quantum phase transitions in strongly interacting two-dimensional systems. We focus first on transitions driven by the condensation of a subset of fractionalized quasiparticles (`anyons) identified with `electric charge excitations of a phase with intrinsic topological order. All other anyons have nontrivial mutual statistics with the condensed subset and hence become confined at the anyon condensation transition. Using a combination of microscopically exact duality transformations and asymptotically exact real-space renormalization group techniques applied to these two-dimensional disordered gauge theories, we argue that the resulting critical scaling behavior is `superuniversal across a wide range of such condensation transitions, and is controlled by the same infinite-randomness fixed point as that of the 2D random transverse-field Ising model. We validate this claim using large-scale quantum Monte Carlo simulations that allow us to extract zero-temperature critical exponents and correlation functions in (2+1)D disordered interacting systems. We discuss generalizations of these results to a large class of ground-state and excited-state topological transitions in systems with intrinsic topological order as well as those where topological order is either protected or enriched by global symmetries. When the underlying topological order and the symmetry group are Abelian, our results provide prototypes for topological phase transitions between distinct many-body localized phases.

Quantum Hall network models as Floquet topological insulators

Network models for equilibrium integer quantum Hall (IQH) transitions are described by unitary scattering matrices, that can also be viewed as representing non-equilibrium Floquet systems. The resulting Floquet bands have zero Chern number, and are instead characterized by a chiral Floquet (CF) winding number. This begs the question: How can a model without Chern number describe IQH systems? We resolve this apparent paradox by showing that non-zero Chern number is recovered from the network model via the energy dependence of network model scattering parameters. This relationship shows that, despite their topologically distinct origins, IQH and CF topology-changing transitions share identical universal scaling properties.

Tracking the quantized information transfer at the edge of a chiral Floquet phase

Two-dimensional arrays of periodically driven qubits can host inherently dynamical topological phases with anomalous chiral edge dynamics. These chiral Floquet phases are formally characterized by a dynamical topological invariant, the chiral unitary index. Introducing a quantity called the chiral mutual information, we show that this invariant can be precisely interpreted in terms of a quantized chiral transfer of quantum information along the edge of the system, and devise a physical setup to measure it.

String order parameters for 1d Floquet Symmetry Protected Topological Phases

Floquet symmetry protected topological (FSPT) phases are non-equilibrium topological phases enabled by time-periodic driving. FSPT phases of 1d chains of bosons, spins, or qubits host dynamically protected edge states that can store quantum information without decoherence, making them promising for use as quantum memories. While FSPT order cannot be detected by any local measurement, here we construct non-local string order parameters that directly measure general 1d FSPT order. We propose a superconducting-qubit array based realization of the simplest Ising-FSPT, which can be implemented with existing quantum computing hardware. We devise an interferometric scheme to directly measure the non-local string order using only simple one- and two- qubit operations and single-qubit measurements.

Logarithmically Slow Relaxation in Quasi-Periodically Driven Random Spin Chains

We simulate the dynamics of a disordered interacting spin chain subject to a quasi-periodic time-dependent drive, corresponding to a stroboscopic Fibonacci sequence of two distinct Hamiltonians. Exploiting the recursive drive structure, we can efficiently simulate exponentially long times. After an initial transient, the system exhibits a long-lived glassy regime characterized by a logarithmically slow growth of entanglement and decay of correlations analogous to the dynamics at the many-body delocalization transition. Ultimately, at long time-scales, which diverge exponentially for weak or rapid drives, the system thermalizes to infinite temperature. The slow relaxation enables metastable dynamical phases, exemplified by a time quasi-crystal in which spins exhibit persistent oscillations with a distinct quasi-periodic pattern from that of the drive. We show that in contrast with Floquet systems, a high-frequency expansion strictly breaks down above fourth order, and fails to produce an effective static Hamiltonian that would capture the pre-thermal glassy relaxation.

An infinite family of 3d Floquet topological paramagnets

We uncover an infinite family of time-reversal symmetric 3d interacting topological insulators of bosons or spins, in time-periodically driven systems, which we term Floquet topological paramagnets (FTPMs). These FTPM phases exhibit intrinsically dynamical properties that could not occur in thermal equilibrium, and are governed by an infinite set of $Z_2$-valued topological invariants, one for each prime number. The topological invariants are physically characterized by surface magnetic domain walls that act as unidirectional quantum channels, transferring quantized packets of information during each driving period. We construct exactly solvable models realizing each of these phases, and discuss the anomalous dynamics of their topologically protected surface states. Unlike previous encountered examples of Floquet SPT phases, these 3d FTPMs are not captured by group cohomology methods, and cannot be obtained from equilibrium classifications simply by treating the discrete time-translation as an ordinary symmetry. The simplest such FTPM phase can feature anomalous $Z_2$ (toric code) surface topological order, in which the gauge electric and magnetic excitations are exchanged in each Floquet period, which cannot occur in a pure 2d system without breaking time reversal symmetry.

Floquet symmetry-protected topological phases in cold atomic systems

We propose and analyze two distinct routes toward realizing interacting symmetry-protected topological (SPT) phases via periodic driving. First, we demonstrate that a driven transverse-field Ising model can be used to engineer complex interactions which enable the emulation of an equilibrium SPT phase. This phase remains stable only within a parametric time scale controlled by the driving frequency, beyond which its topological features break down. To overcome this issue, we consider an alternate route based upon realizing an intrinsically Floquet SPT phase that does not have any equilibrium analog. In both cases, we show that disorder, leading to many-body localization, prevents runaway heating and enables the observation of coherent quantum dynamics at high energy densities. Furthermore, we clarify the distinction between the equilibrium and Floquet SPT phases by identifying a unique micromotion-based entanglement spectrum signature of the latter. Finally, we propose a unifying implementation in a one-dimensional chain of Rydberg-dressed atoms and show that protected edge modes are observable on realistic experimental time scales.

Discrete time crystals: rigidity, criticality, and realizations

Despite being forbidden in equilibrium, spontaneous breaking of time translation symmetry can occur in periodically driven, Floquet systems with discrete time-translation symmetry. The period of the resulting discrete time crystal is quantized to an integer multiple of the drive period, arising from a combination of collective synchronization and many body localization. Here, we consider a simple model for a one dimensional discrete time crystal which explicitly reveals the rigidity of the emergent oscillations as the drive is varied. We numerically map out its phase diagram and compute the properties of the dynamical phase transition where the time crystal melts into a trivial Floquet insulator. Moreover, we demonstrate that the model can be realized with current experimental technologies and propose a blueprint based upon a one dimensional chain of trapped ions. Using experimental parameters (featuring long-range interactions), we identify the phase boundaries of the ion-time-crystal and propose a measurable signature of the symmetry breaking phase transition.