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 tha
t, 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.
We show that a class of random all-to-all spin models, realizable in systems of atoms coupled to an optical cavity, gives rise to a rich dynamical phase diagram due to the pairwise separable nature of the couplings. By controlling the experimental pa
rameters, one can tune between integrable and chaotic dynamics on the one hand, and between classical and quantum regimes on the other hand. For two special values of a spin-anisotropy parameter, the model exhibits rational-Gaudin type integrability and it is characterized by an extensive set of spin-bilinear integrals of motion, independent of the spin size. More generically, we find a novel integrable structure with conserved charges that are not purely bilinear. Instead, they develop `dressing tails of higher-body terms, reminiscent of the dressed local integrals of motion found in Many-Body Localized phases. Surprisingly, this new type of integrable dynamics found in finite-size spin-1/2 systems disappears in the large-$S$ limit, giving way to classical chaos. We identify parameter regimes for characterizing these different dynamical behaviors in realistic experiments, in light of the limitations set by cavity dissipation.
Recent work by De Roeck et al. [Phys. Rev. B 95, 155129 (2017)] has argued that many-body localization (MBL) is unstable in two and higher dimensions due to a thermalization avalanche triggered by rare regions of weak disorder. To examine these argum
ents, we construct several models of a finite ergodic bubble coupled to an Anderson insulator of non-interacting fermions. We first describe the ergodic region using a GOE random matrix and perform an exact diagonalization study of small systems. The results are in excellent agreement with a refined theory of the thermalization avalanche that includes transient finite-size effects, lending strong support to the avalanche scenario. We then explore the limit of large system sizes by modeling the ergodic region via a Hubbard model with all-to-all random hopping: the combined system, consisting of the bubble and the insulator, can be reduced to an effective Anderson impurity problem. We find that the spectral function of a local operator in the ergodic region changes dramatically when coupling to a large number of localized fermionic states---this occurs even when the localized sites are weakly coupled to the bubble. In principle, for a given size of the ergodic region, this may arrest the avalanche. However, this back-action effect is suppressed and the avalanche can be recovered if the ergodic bubble is large enough. Thus, the main effect of the back-action is to renormalize the critical bubble size.
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 wh
ich 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.
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.
Classical satisfiability (SAT) and quantum satisfiability (QSAT) are complete problems for the complexity classes NP and QMA which are believed to be intractable for classical and quantum computers, respectively. Statistical ensembles of instances of
these problems have been studied previously in an attempt to elucidate their typical, as opposed to worst case, behavior. In this paper we introduce a new statistical ensemble that interpolates between classical and quantum. For the simplest 2-SAT/2-QSAT ensemble we find the exact boundary that separates SAT and UNSAT instances. We do so by establishing coincident lower and upper bounds, in the limit of large instances, on the extent of the UNSAT and SAT regions, respectively.
We numerically study the interplay between superconductivity and disorder on the graphene honeycomb lattice with on-site Hubbard attractive interactions U using a spatially inhomogeneous self-consistent Bogoliubov-de Gennes (BdG) approach. In the abs
ence of disorder there are two phases at charge neutrality. Below a critical value Uc for attractive interactions there is a Dirac semimetal phase and above it there is a superconducting phase. We add scalar potential disorder to the system, while remaining at charge neutrality on average. Numerical solution of the BdG equations suggests that while in the strong attraction regime (U > Uc) disorder has the usual effect of suppressing superconductivity, in the weak attraction regime (U < Uc) weak disorder enhances superconductivity. In the weak attraction regime, disorder that is too strong eventually suppresses superconductivity, i.e., there is an optimal disorder strength that maximizes the critical temperature Tc. Our numerical results also suggest that in the weakly disordered regime, mesoscopic inhomogeneities enhance superconductivity significantly more than what is predicted by a spatially uniform mean-field theory a` la Abrikosov-Gorkov. In this regime, superconductivity consists of rare phase-coherent superconducting islands. We also study the enhancement of the superconducting proximity effect by disorder and mesoscopic inhomogeneities, and obtain typical spatial plots of the tunneling density of states and the superfluid susceptibility that can be directly compared to scanning tunneling miscroscopy (STM) experiments on proximity-induced superconductivity in graphene.
