Search

Nonintegrable many-body quantum systems typically thermalize at long times through the mechanism of quantum chaos. However, some exceptional systems, such as those harboring quantum scars, break thermalization, serving as testbeds for foundational problems of quantum statistical physics. Here, we investigate a class of nonintegrable bond-staggered models that is endowed with a large number of zero-energy eigenstates and possesses a non-Abelian internal symmetry. We use character theory to give a lower bound on the zero-energy degeneracy, which matches exact diagonalization results, and is found to grow exponentially with the system size. We also show that few-magnon zero-energy states have an exact analytical description, allowing us to build a basis of low-entangled fixed-separation states, which is stable to most perturbations found in experiments. This remarkable dynamical stability of special states elucidates our understanding of nonequilibrium processes in non-Abelian chaotic quantum models., Comment: 10 pages, 7 figures

Periodically driven quantum systems host exotic phenomena which often do not have any analog in undriven systems. Floquet prethermalization and dynamical freezing of certain observables, via the emergence of conservation laws, are realized by controlling the drive frequency. These dynamical regimes can be leveraged to construct quantum memories and have potential applications in quantum information processing. Solid state and cold atom experimental architectures have opened avenues for implementing local Floquet engineering which can achieve spatially modulated quantum control of states. Here, we uncover the novel memory effects of local periodic driving in a nonintegrable spin-half staggered Heisenberg chain. For a boundary-driven protocol at the dynamical freezing frequency, we show the formation of an approximate strong zero mode, a prethermal quasi-local operator, due to the emergence of a discrete global $\mathbb{Z}_2$ symmetry. This is captured by constructing an accurate effective Floquet Hamiltonian using a higher-order partially resummed Floquet-Magnus expansion. The lifetime of the boundary spin can be exponentially enhanced by enlarging the set of suitably chosen driven sites. We demonstrate that in the asymptotic limit, achieved by increasing the number of driven sites, a strong zero mode emerges, where the lifetime of the boundary spin grows exponentially with system size. The non-local processes in the Floquet Hamiltonian play a pivotal role in the total freezing of the boundary spin in the thermodynamic limit. The novel dynamics of the boundary spin is accompanied by a rich structure of entanglement in the Floquet eigenstates where specific bipartitions yield an area-law scaling while the entanglement for random bipartitions scales as a volume-law., Comment: 22 pages, 15 figures, 1 table

Scar eigenstates in a many-body system refers to a small subset of non-thermal finite energy density eigenstates embedded into an otherwise thermal spectrum. This novel non-thermal behaviour has been seen in recent experiments simulating a one-dimensional PXP model with a kinetically-constrained local Hilbert space realized by a chain of Rydberg atoms. We probe these small sets of special eigenstates starting from particular initial states by computing the spread complexity associated to time evolution of the PXP hamiltonian. Since the scar subspace in this model is embedded only loosely, the scar states form a weakly broken representation of the Lie Algebra. We demonstrate why a careful usage of the Forward Scattering Approximation (or similar strategies thereof) is required to extract an appropriate set of Lanczos coefficients in this case as the consequence of this approximate symmetry. This leads to a well defined notion of a closed Krylov subspace and consequently, that of spread complexity. We show how the spread complexity shows approximate revivals starting from both $|\mathbb{Z}_2\rangle$ and $|\mathbb{Z}_3\rangle$ states and how these revivals can be made more accurate by adding optimal perturbations to the bare Hamiltonian. We also investigate the case of the vacuum as the initial state, where revivals can be stabilized using an iterative process of adding few-body terms., Comment: 19 pages, 11 figures

Thermalization and its breakdown in isolated systems has led to a deeper understanding of non-equilibrium quantum states and their dependence on initial conditions. This is prominently highlighted by the existence of quantum scars, special athermal states with an underlying effective superspin structure, embedded in an otherwise chaotic many-body spectrum. Spin $XXZ$ models and their variants in one and higher dimension have been shown to host exact quantum scars, exhibiting perfect revivals of spin helix states that are realizable in synthetic and condensed matter systems. Motivated by these advances, we propose experimentally accessible local, time-dependent protocols to explore the spatial thermalization profile and highlight how different parts of the system thermalize and affect the fate of the superspin. We identify distinct parametric regimes for the ferromagnetic ($X$-polarized) initial state based on the interplay between the driven spin and the rest, including local athermal behavior where the driven spin effectively decouples, acting like a "cold" spot while being instrumental in heating up the other spins. We develop a real and Floquet space picture that explains our numerical observations, and make predictions that can be tested in various experimental setups., Comment: 4+ pages, 4 figures. 2 pages of Supplemental Material

We introduce a disorder-free model of $S=1/2$ spins on the square lattice in a constrained Hilbert space where two up-spins are not allowed simultaneously on any two neighboring sites of the lattice. The interactions are given by ring-exchange terms on elementary plaquettes that conserve both the total magnetization as well as dipole moment. We show that this model provides a tractable example of strong Hilbert space fragmentation in two dimensions with typical initial states evading thermalization with respect to the full Hilbert space. Given any product state, the system can be decomposed into disjoint spatial regions made of edge and/or vertex sharing plaquettes that we dub as ``quantum drums''. These quantum drums come in many shapes and sizes and specifying the plaquettes that belong to a drum fixes its spectrum. The spectra of some small drums is calculated analytically. We study two bigger quasi-one-dimensional drums numerically, dubbed ``wire'' and a ``junction of two wires'' respectively. We find that these possess a chaotic spectrum but also support distinct families of quantum many-body scars that cause periodic revivals from different initial states. The wire is shown to be equivalent to the one-dimensional PXP chain with open boundaries, a paradigmatic model for quantum many-body scarring; while the junction of two wires represents a distinct constrained model., Comment: v3: 57 pages, 30 figs; resubmission to SciPost Physics

Periodic driving has been established as a powerful technique for engineering novel phases of matter and intrinsically out-of-equilibrium phenomena such as time crystals. Recent work by Bluvstein et al. [Science 371, 1355 (2021)] has demonstrated that periodic driving can also lead to a significant enhancement of quantum many-body scarring, whereby certain non-integrable systems can display persistent quantum revivals from special initial states. Nevertheless, the mechanisms behind driving-induced scar enhancement remain poorly understood. Here we report a detailed study of the effect of periodic driving on the PXP model describing Rydberg atoms in the presence of a strong Rydberg blockade - the canonical static model of quantum many-body scarring. We show that periodic modulation of the chemical potential gives rise to a rich phase diagram, with at least two distinct types of scarring regimes that we distinguish by examining their Floquet spectra. We formulate a toy model, based on a sequence of square pulses, that accurately captures the details of the scarred dynamics and allows for analytical treatment in the large-amplitude and high-frequency driving regimes. Finally, we point out that driving with a spatially inhomogeneous chemical potential allows to stabilize quantum revivals from arbitrary initial states in the PXP model, via a mechanism similar to prethermalization., Comment: 21 pages, 17 figures

We study the stroboscopic dynamics of a periodically driven finite Rydberg chain with staggered ($\Delta$) and time-dependent uniform ($\lambda(t)$) detuning terms using exact diagonalization (ED). We show that at intermediate drive frequencies ($\omega_D$), the presence of a finite $\Delta$ results in violation of the eigenstate thermalization hypothesis (ETH) via clustering of Floquet eigenstates. Such clustering is lost at special commensurate drive frequencies for which $\hbar \omega_d=n \Delta$ ($n \in Z$) leading to restoration of ergodicity. The violation of ETH in these driven finite-sized chains is also evident from the dynamical freezing displayed by the density-density correlation function at specific $\omega_D$. Such a correlator exhibits stable oscillations with perfect revivals when driven close to the freezing frequencies for initial all spin-down ($|0\rangle$) or Neel ($|{\mathbb Z}_2\rangle$, with up-spins on even sites) states. The amplitudes of these oscillations vanish at the freezing frequencies and reduces upon increasing $\Delta$; their frequencies, however, remains pinned to $\Delta/\hbar$ in the large $\Delta$ limit. In contrast, for the $|{\bar {\mathbb Z}_2}\rangle$ (time-reversed partner of $|{\mathbb Z}_2\rangle$) initial state, we find complete absence of such oscillations leading to freezing for a range of $\omega_D$; this range increases with $\Delta$. We also study the properties of quantum many-body scars in the Floquet spectrum of the model as a function of $\Delta$ and show the existence of novel mid-spectrum scars at large $\Delta$. We supplement our numerical results with those from an analytic Floquet Hamiltonian computed using Floquet perturbation theory (FPT) and also provide a semi-analytic computation of the quantum scar states within a forward scattering approximation (FSA)., Comment: v1; 19 pages, 19 figs

Motivated by previous works on a Floquet version of the PXP model [Mukherjee {\it et al.} Phys. Rev. B 102, 075123 (2020), Mukherjee {\it et al.} Phys. Rev. B 101, 245107 (2020)], we study a one-dimensional spin-$1/2$ lattice model with three-spin interactions in the same constrained Hilbert space (where all configurations with two adjacent $S^z=\uparrow$ spins are excluded). We show that this model possesses an extensive fragmentation of the Hilbert space which leads to a breakdown of thermalization upon unitary evolution starting from a large class of simple initial states. Despite the non-integrable nature of the Hamiltonian, many of its high-energy eigenstates admit a quasiparticle description. A class of these, which we dub as "bubble eigenstates", have integer eigenvalues (including mid-spectrum zero modes) and strictly localized quasiparticles while another class contains mobile quasiparticles leading to a dispersion in momentum space. Other anomalous eigenstates that arise due to a {\it secondary} fragmentation mechanism, including those that lead to flat bands in momentum space due to destructive quantum interference, are also discussed. The consequences of adding a (non-commuting) staggered magnetic field and a PXP term respectively to this model, where the former preserves the Hilbert space fragmentation while the latter destroys it, are discussed. A Floquet version with time-dependent staggered field also evades thermalization with additional features like freezing of exponentially many states at special drive frequencies. Finally, we map the model to a $U(1)$ lattice gauge theory coupled to dynamical fermions and discuss the interpretation of some of these anomalous states in this language. A class of gauge-invariant states show reduced mobility of the elementary charged excitations with only certain charge-neutral objects being mobile suggesting a connection to fractons., Comment: v3; 22 pages, 11 figures, corrected some typos, close to published version

Eigenstate Thermalization Hypothesis(ETH) has played a pivotal role in understanding ergodicity and its breaking in isolated quantum many-body systems. Recent experiment on 51-atom Rydberg quantum simulator and subsequent theoretical analysis have shown that hardcore kinetic constraint can lead to weak ergodicity breaking. In this work, we demonstrate, using 1d spin-1 PXP chains, that miscellaneous type of ergodicity can be realized by adjusting the hardcore constraints between different components of nearest neighbor spins. This includes ETH violation due to emergent shattering of Hilbert space into exponentially many subsectors of various sizes, a novel form of non-integrability with an extensive number of local conserved quantities and strong ergodicity. We analyze these different forms of ergodicity and study their impact on the non-equilibrium dynamics of a Z2 initial state. We use forward scattering approximation (FSA) to understand the amount of Z2-oscillation present in these models. Our work shows that not only ergodicity breaking but an appropriate choice of constraints can lead to restoration of ergodicity as well.

In this study, an order by disorder mechanism has been proposed in a two-dimensional PXP model, where the extensive degeneracy of the classical ground-state manifold is due to strict occupation constraints instead of geometrical frustrations. By performing an unbias large-scale quantum monte carlos simulation, we find that local quantum fluctuations, which usually work against long-range ordering, lift the macroscopic classical degeneracy and give rise to a compressible ground state with charge-density-wave long-range order. A simple trial wavefunction has been proposed to capture the essence of the ground-state of the two-dimensional PXP model. The finite temperature properties of this model have also been studied, and we find a thermal phase transition with an universality class of two-dimensional Ising model., Comment: 4 pages

We compute the Floquet Hamiltonian $H_F$ for weakly interacting fermions subjected to a continuous periodic drive using a Floquet perturbation theory (FPT) with the interaction amplitude being the perturbation parameter. This allows us to address the dynamics of the system at intermediate drive frequencies $\hbar \omega_D \ge V_0 \ll {\mathcal J}_0$, where ${\mathcal J}_0$ is the amplitude of the kinetic term, $\omega_D$ is the drive frequency, and $V_0$ is the typical interaction strength between the fermions. We compute, for random initial states, the fidelity $F$ between wavefunctions after a drive cycle obtained using $H_F$ and that obtained using exact diagonalization (ED). We find that FPT yields a substantially larger value of $F$ compared to its Magnus counterpart for $V_0\le \hbar \omega_D$ and $V_0\ll {\mathcal J}_0$. We use the $H_F$ obtained to study the nature of the steady state of an weakly interacting fermion chain; we find a wide range of $\omega_D$ which leads to subthermal or superthermal steady states for finite chains. The driven fermionic chain displays perfect dynamical localization for $V_0=0$; we address the fate of this dynamical localization in the steady state of a finite interacting chain and show that there is a crossover between localized and delocalized steady states. We discuss the implication of our results for thermodynamically large chains and chart out experiments which can test our theory., Comment: 11 pages,8 figures

We study the dynamics of the periodically driven Rydberg chain starting from the state with zero Rydberg excitations (vacuum state denoted by $|0\rangle$) using a square pulse protocol in the high drive amplitude limit. We show, using exact diagonalization for finite system sizes ($L\le 26$), that the Floquet Hamiltonian of the system, within a range of drive frequencies which we chart out, hosts a set of quantum scars which have large overlap with the $|0\rangle$ state. These scars are distinct from their counterparts having high overlap with the maximal Rydberg excitation state ($|\mathbb{Z}_2\rangle$); they coexist with the latter class of scars and lead to persistent coherent oscillations of the density-density correlator starting from the $|0\rangle$ state. We also identify special drive frequencies at which the system undergoes perfect dynamic freezing and provide an analytic explanation for this phenomenon. Finally, we demonstrate that for a wide range of drive frequencies, the system reaches a steady state with sub-thermal values of the density-density correlator. The presence of such sub-thermal steady states, which are absent for dynamics starting from the $|\mathbb{Z}_2\rangle$ state, imply a weak violation of the eigenstate thermalization hypothesis in finite sized Rydberg chains distinct from that due to the scar-induced persistent oscillations reported earlier. We conjecture that in the thermodynamic limit such states may exist as pre-thermal steady states that show anomalously slow relaxation. We supplement our numerical results by deriving an analytic expression for the Floquet Hamiltonian using a Floquet perturbation theory in the high amplitude limit which provides an analytic, albeit qualitative, understanding of these phenomena at arbitrary drive frequencies. We discuss experiments which can test our theory., Comment: 18 pages and 17 figures including appendices

We study the unitary dynamics of randomly or quasi-periodically driven tilted Bose-Hubbard (tBH) model in one dimension deep inside its Mott phase starting from a $\mathbb{Z}_2$ symmetry-broken state. The randomness is implemented via a telegraph noise protocol in the drive period while the quasi-periodic drive is chosen to correspond to a Thue-Morse sequence. The periodically driven tBH model (with a square pulse protocol characterized by a time period $T$) is known to exhibit transitions from dynamical regimes with long-time coherent oscillations to those with rapid thermalization. Here we show that starting from a regime where the periodic drive leads to rapid thermalization, a random drive, which consists of a random sequence of square pulses with period $T+\alpha dT$, where $\alpha=\pm 1$ is a random number and $dT$ is the amplitude of the noise, restores long-time coherent oscillations for special values of $dT$. A similar phenomenon can be seen for a quasi-periodic drive following a Thue-Morse sequence where such coherent behavior is shown to occur for a larger number of points in the $(T, dT)$ plane due to the additional structure of the drive protocol. We chart out the dynamics of the system in the presence of such aperiodic drives, provide a qualitative analytical understanding of this phenomenon, point out the role of quantum scars behind it, and discuss experiments which can test our theory., Comment: v2; 10 pages, 11 Figs; references corrected

The presence of quantum scars, athermal eigenstates of a many-body Hamiltonian with finite energy density, leads to absence of ergodicity and long-time coherent dynamics in closed quantum systems starting from simple initial states. Such non-ergodic coherent dynamics, where the system does not explore its entire phase space, has been experimentally observed in a chain of ultracold Rydberg atoms. We show, via study of a periodically driven Rydberg chain, that the drive frequency acts as a tuning parameter for several reentrant transitions between ergodic and non-ergodic regimes. The former regime shows rapid thermalization of correlation functions and absence of scars in the spectrum of the system's Floquet Hamiltonian. The latter regime, in contrast, has scars in its Floquet spectrum which control the long-time coherent dynamics of correlation functions. Our results open a new possibility of drive frequency-induced tuning between ergodic and non-ergodic dynamics in experimentally realizable disorder-free quantum many-body systems., Comment: v2; 14 pages and 8 figures including supplementary material. Expanded version with additional results and discussions; a co-author added

We study Floquet topological transition in irradiated graphene when the polarization of incident light changes randomly with time. We numerically confirm that the noise averaged time evolution operator approaches a steady value in the limit of exact Trotter decomposition of the whole period where incident light has different polarization at each interval of the decomposition. This steady limit is found to coincide with time-evolution operator calculated from the noise-averaged Hamiltonian. We observe that at the six corners (Dirac($K$) point) of the hexagonal Brillouin zone of graphene random Gaussian noise strongly modifies the phaseband structure induced by circularly polarized light whereas in zone-center ($\Gamma$ point) even a strong noise isn't able to do the same. This can be understood by analyzing the deterministic noise averaged Hamiltonian which has a different Fourier structure as well as lesser no of symmetries compared to the noise-free one. In 1D systems noise is found to renormalize the drive amplitude only., Comment: 8 Pages, 8 figures

We study the Floquet phase diagram of two-dimensional Dirac materials such as graphene and the one-dimensional (1D) spin-1/2 $XY$ model in a transverse field in the presence of periodic time-varying terms in their Hamiltonians in the low drive frequency ($\omega$) regime where standard $1/\omega$ perturbative expansions fail. For graphene, such periodic time dependent terms are generated via the application of external radiation of amplitude $A_0$ and time period $T = 2\pi/\omega$, while for the 1D $XY$ model, they result from a two-rate drive protocol with time-dependent magnetic field and nearest-neighbor couplings between the spins. Using the adiabatic-impulse method, we provide several semi-analytic criteria for the occurrence of changes in the topology of the phase bands of such systems. For irradiated graphene, we point out the role of the symmetries of $H(t)$ and $U$ behind such topology changes. Our analysis reveals that at low frequencies, phase band topology changes may also happen at $t= T/3, 2T/3$ (apart from $t=T$). We chart out the phase diagrams at $t=T/3, 2T/3,\, {\rm and }\, T$ as a function of $A_0$ and $T$ using exact numerics, and compare them with the prediction of the adiabatic-impulse method. We show that several characteristics of these phase diagrams can be analytically understood from results obtained using the adiabatic-impulse method and point out the crucial contribution of the high-symmetry points in the graphene Brillouin zone to these diagrams. Finally we study the 1D $XY$ model with a two-rate driving protocol using the adiabatic-impulse method and exact numerics revealing a phase band crossing at $t=T/2$ and $k=\pi/2$. We also study the anomalous end modes generated by such a drive. We suggest experiments to test our theory., Comment: v1; 26 pages, 19 Figs

We chart out the phase diagram of ultracold `spin-half' bosons in a one-dimensional optical lattice in the presence of Aubry-Andr\'e (AA) potential and with spin-orbit (SO) and Raman couplings investigating the transition from superfluid (SF) to localized phases and the existence of density wave phase for nearest-neighbor interaction (NNI). We show that the presence of SO coupling and AA potential leads to a novel spin-split momentum distribution of the bosons in the localized phase near the boundary with the SF phase, which can act as a signature of such a transition. We also obtain the level statistics of the bosons in the superfluid phase with finite NNI and demonstrate its change from Gaussian Unitary Ensemble (GUE) to Gaussian Orthogonal Ensemble (GOE) as a function of the Raman coupling. We discuss experiments which can test our theory., Comment: 4 pages, 4 Figs, Supplementary material; v1

We present generic conditions for phase band crossings for a class of periodically driven integrable systems represented by free fermionic models subjected to arbitrary periodic drive protocols characterized by a frequency $\omega_D$. These models provide a representation for the Ising and $XY$ models in $d=1$, the Kitaev model in $d=2$, several kinds of superconductors, and Dirac fermions in graphene and atop topological insulator surfaces. Our results demonstrate that the presence of a critical point/region in the system Hamiltonian (which is traversed at a finite rate during the dynamics) may change the conditions for phase band crossings that occur at the critical modes. We also show that for $d>1$, phase band crossings leave their imprint on the equal-time off-diagonal fermionic correlation functions of these models; the Fourier transforms of such correlation functions, $F_{\vec k_0}( \omega_0)$, have maxima and minima at specific frequencies which can be directly related to $\omega_D$ and the time at which the phase bands cross at $\vec k = \vec k_0$. We discuss the significance of our results in the contexts of generic Hamiltonians with $N>2$ phase bands and the underlying symmetry of the driven Hamiltonian., Comment: v1; 7 figs, 12 pages

We study periodically driven closed quantum systems where two parameters of the system Hamiltonian are driven with frequencies $\omega_1$ and $\omega_2=r \omega_1$. We show that such drives may be used to tune towards dynamics induced freezing where the wavefunction of the state of the system after a drive cycle at time $T= 2\pi/\omega_1$ has almost perfect overlap with the initial state. We locate regions in the $(\omega_1 ,r)$ plane where the freezing is near exact for a class of integrable and a specific non-integrable model. The integrable models that we study encompass Ising and XY models in $d=1$, Kitaev model in $d=2$, and Dirac fermions in graphene and atop a topological insulator surface whereas the non-integrable model studied involves the experimentally realized one-dimensional (1D) tilted Bose-Hubbard model in an optical lattice. In addition, we compute the relevant correlation functions of such driven systems and describe their characteristics in the region of $(\omega_1,r)$ plane where the freezing is near-exact. We supplement our numerical analysis with semi-analytic results for integrable driven systems within adiabatic-impulse approximation and discuss experiments which may test our theory., Comment: v1: 9 pages 9 figs