Slow Quantum Thermalization and Many-Body Revivals from Mixed Phase Space

Quantum many-body scar states are special eigenstates of nonintegrable models with distinctive entanglement features that give rise to infinitely long-lived coherent dynamics under quantum quenches from certain initial states. We elaborate on a construction of quantum many-body scar states in which they emerge from Einstein-Podolsky-Rosen states in systems with two layers, wherein the two layers are maximally entangled. We apply this construction to spin systems as well as systems of itinerant fermions and bosons and demonstrate how symmetries can be harnessed to enhance its versatility. We show that several well-known examples of quantum many-body scars, including the tower of states in the spin-1 $XY$ model and the $\ensuremath{\eta}$-pairing states in the Fermi-Hubbard model, can be understood within this formalism. We also demonstrate how an infinite tower of many-body scar states can emerge in bilayer Bose-Hubbard models with charge conservation.

We review the fictitious integrable system approach which predicts dynamical tunneling rates from regular states to the chaotic region in systems with a mixed phase space. It is based on the introduction of a fictitious integrable system that resembles the regular dynamics within the regular island. We focus on the direct regular-to-chaotic tunneling process which dominates if nonlinear resonances within the regular island are not relevant. For quantum maps, billiard systems, and optical microcavities, we find excellent agreement with numerical rates for all regular states.

Dynamics of a system in general depends on its initial state and how the system is driven, but in many-body systems the memory is usually averaged out during evolution. Here, interacting quantum systems without external relaxations are shown to retain long-time memory effects in steady states. To identify memory effects, we first show quasi-steady state currents form in finite, isolated Bose and Fermi Hubbard models driven by interaction imbalance and they become steady-state currents in the thermodynamic limit. By comparing the steady state currents from different initial states or ramping rates of the imbalance, long-time memory effects can be quantified. While the memory effects of initial states are more ubiquitous, the memory effects of switching protocols are mostly visible in interaction-induced transport in lattices. Our simulations suggest the systems enter a regime governed by a generalized Fick's law and memory effects lead to initial-state dependent diffusion coefficients. We also identify conditions for enhancing memory effects and discuss possible experimental implications.

Whereas much work in the mathematical literature on quantum chaos has focused on phenomena such as quantum ergodicity and scarring, relatively little is known at the rigorous level about the existence of eigenfunctions whose morphology is more complex. Quantum systems whose dynamics is intermediate between certain regimes - for example, at the transition between Anderson localized and delocalized eigenfunctions, or in systems whose classical dynamics is intermediate between integrability and chaos - have been conjectured in the physics literature to have eigenfunctions exhibiting multifractal, self-similar structure. To-date, no rigorous mathematical results have been obtained about systems of this kind in the context of quantum chaos. We give here the first rigorous proof of the existence of multifractal eigenfunctions for a widely studied class of intermediate quantum systems. Specifically, we derive an analytical formula for the Renyi entropy associated with the eigenfunctions of arithmetic Seba billiards, in the semiclassical limit, as the associated eigenvalues tend to infinity. We also prove multifractality of the ground state for more general, non-arithmetic billiards and show that the fractal exponent in this regime satisfies a symmetry relation, similar to the one predicted in the physics literature, by establishing a connection with the functional equation for Epstein's zeta function.

The quantum dynamics of initial coherent states is studied in the Dicke model and correlated with the dynamics, regular or chaotic, of their classical limit. Analytical expressions for the survival probability, i.e. the probability of finding the system in its initial state at time t, are provided in the regular regions of the model. The results for regular regimes are compared with those of the chaotic ones. It is found that initial coherent states in regular regions have a much longer equilibration time than those located in chaotic regions. The properties of the distributions for the initial coherent states in the Hamiltonian eigenbasis are also studied. It is found that for regular states the components with no negligible contribution are organized in sequences of energy levels distributed according to Gaussian functions. In the case of chaotic coherent states, the energy components do not have a simple structure and the number of participating energy levels is larger than in the regular cases.

The localization of eigenfunctions around classical periodic orbits is studied numerically for the H-atom in a strong magnetic field by calculating their Husimi distribution in phase space. In contrast to the configuration space representation, the phase space distributions are simply structured: about 90% of eigenstates may be unambigously related to fixed points and invariant manifolds of periodic orbits, indicating that scars are the rule rather than the exception. In order to measure the influence of one particular orbit, we calculate the integrals of the energetically lowest 500 Husimi distributions along the orbit. Their incoherent superposition defines the scar strength distribution for the particular periodic orbit which is analyzed by Fourier transformation. The Husimi distribution at (q, p) in phase space may be represented as a scalar product of the wavefunction with a coherent state of the unperturbed system, i.e., a radial Gaussian wave packet located at (q, p) in the (regularized) Coulomb system. This simplifies the actual calculation of the Husimi distribution and allows to treat their incoherent superposition within Gutzwillers theory extended to matrix elements of an operator A, if we choose A to be the projector on a coherent state.

Understanding stickiness and power-law behavior of Poincaré recurrence statistics is an open problem for higher-dimensional systems, in contrast to the well-understood case of systems with two degrees of freedom. We study such intermittent behavior of chaotic orbits in three-dimensional volume-preserving maps using the example of the Arnold-Beltrami-Childress map. The map has a mixed phase space with a cylindrical regular region surrounded by a chaotic sea for the considered parameters. We observe a characteristic overall power-law decay of the cumulative Poincaré recurrence statistics with significant oscillations superimposed. This slow decay is caused by orbits which spend long times close to the surface of the regular region. Representing such long-trapped orbits in frequency space shows clear signatures of partial barriers and reveals that coupled resonances play an essential role. Using a small number of the most relevant resonances allows for classifying long-trapped orbits. From this the Poincaré recurrence statistics can be divided into different exponentially decaying contributions, which very accurately explains the overall power-law behavior including the oscillations.

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.

Understanding the behavior of quantum many-body systems under decoherence is essential for developing robust quantum technologies. Here, we examine the fate of weak ergodicity breaking in systems hosting quantum many-body scars when subject to local pure dephasing-an experimentally relevant form of environmental noise. Focusing on a large class of models with an approximate su(2)-structured scar subspace, we show that scarred eigenmodes of the Liouvillian exhibit a transition reminiscent of spontaneous PT-symmetry breaking as the dephasing strength increases. Unlike previously studied non-Hermitian mechanisms, this transition arises from a distinct quantum jump effect. Remarkably, in platforms such as the XY spin ladder and PXP model of Rydberg atom arrays, the critical dephasing rate shows only weak dependence on the system size, revealing an unexpected robustness of scarred dynamics in noisy environments.

Motivated by the existence of exact many-body quantum scars in the Affleck-Kennedy-Lieb-Tasaki (AKLT) chain, in this work we explore the connection between matrix product state (MPS) wave functions and many-body quantum scarred Hamiltonians. We provide a method to systematically search for and construct parent Hamiltonians with towers of exact eigenstates composed of quasiparticles on top of an MPS wave function. These exact eigenstates have low entanglement in spite of being in the middle of the spectrum, thus violating the strong eigenstate thermalization hypothesis. Using our approach, we recover the AKLT chain starting from the MPS of its ground state, and we derive the most general nearest-neighbor Hamiltonian that shares the AKLT quasiparticle tower of exact eigenstates. We further apply this formalism to other simple MPS wave functions, and derive families of Hamiltonians that exhibit AKLT-like quantum scars. As a consequence, we also construct a scar-preserving deformation that connects the AKLT chain to the integrable spin-1 pure biquadratic model. Finally, we also derive other families of Hamiltonians that exhibit types of exact quantum scars, including a $\\textit{U}$(1)-invariant perturbed Potts model.

A quantum many-body scar system usually contains a special non-thermal subspace (approximately) decoupled from the rest of the Hilbert space. In this work, we propose a general structure called deformed symmetric spaces for the decoupled subspaces hosting quantum many-body scars, which are irreducible sectors of simple Lie groups transformed by matrix-product operators (or projected entangled pair operators), of which the entanglement entropies are proved to obey sub-volume-law scaling and thus violate the eigenstate thermalization hypothesis. A deformed symmetric space, in general, is required to have at least a U(1) sub-Lie-group symmetry to allow coherent periodic dynamics from certain low-entangled initial states. We enumerate several possible deforming transformations based on the sub-group symmetry requirement and recover many existing models whose scar states are not connected by symmetry. In particular, a two-dimensional scar model is proposed, which hosts a periodic dynamical trajectory on which all states are topologically ordered.

The lifetimes of optical modes in whispering-gallery cavities depend crucially on the underlying classical ray dynamics, and they may be spoiled by the presence of classical nonlinear resonances due to resonance-assisted tunneling. Here we present an intuitive semiclassical picture that allows for an accurate prediction of decay rates of optical modes in systems with a mixed phase space. We also extend the perturbative description from near-integrable systems to systems with a mixed phase space, and we find equally good agreement. Both approaches are based on the approximation of the actual ray dynamics by an integrable Hamiltonian, which enables us to perform a semiclassical quantization of the system and to introduce a ray-based description of the decay of optical modes. The coupling between them is determined either perturbatively or semiclassically in terms of complex paths.

The concept of quantum many-body scars has recently been put forward as a route to describe weak ergodicity breaking and violation of the eigenstate thermalization hypothesis. We propose a simple setup to generate quantum many-body scars in a doubly modulated Bose-Hubbard system which can be readily implemented in cold atomic gases. The dynamics are shown to be governed by kinetic constraints which appear via density-assisted tunneling in a high-frequency expansion. We find the optimal driving parameters for the kinetically constrained hopping which leads to small isolated subspaces of scared eigenstates. The experimental signatures and the transition to fully thermalizing behavior as a function of driving frequency are analyzed.

Experimental observation of coherent oscillations in a Rydberg atom chain [H. Bernien et al., Nature (London) 551, 579 (2017)] has led to the discovery of quantum many-body scars (QMBSs) which is a new paradigm for ergodicity breaking. The experimental findings in the Rydberg chain can be well captured by a kinetically constrained model called the ``PXP'' model, which has been shown to host the eigenstate thermalization hypothesis (ETH)-violating scar states in the middle of the spectrum. Much effort has been put into identifying similar kinetically restricted systems that show a violation of ETH. In this paper, we study the QMBSs that can arise in one such model, namely the spin-1 Kitaev chain, where owing to some conserved quantities, the Hilbert space gets fragmented into unequal disconnected subspaces. Recently, You et al. [Phys. Rev. Res. 4, 013103 (2022)] showed that the ground state sector of this chain can be mapped exactly onto the prototypical PXP model and thus hosts QMBSs. Here, we demonstrate that the phenomenon of scarring is also present in other sectors, and in particular, we identify a sector that exhibits substantially more scarring than the ground state one. We propose an initial state and numerically demonstrate that its fidelity revivals are pronounced and longer lived than those in the PXP model.

Dynamical decoupling and Hamiltonian engineering are well-established techniques that have been used to control qubit systems. However, designing the corresponding methods for qudit systems has been challenging due to the lack of a Bloch sphere representation, more complex interactions, and additional control constraints. By identifying several general structures associated with such problems, we develop a formalism for the robust dynamical decoupling and Hamiltonian engineering of strongly interacting qudit systems. Our formalism significantly simplifies qudit pulse-sequence design while naturally incorporating robustness conditions necessary for experimental practicality. We experimentally demonstrate these techniques in a strongly interacting, disordered ensemble of spin-1 nitrogen-vacancy centers, achieving more than an order-of-magnitude improvement in coherence time over existing pulse sequences. We further describe how our techniques enable the engineering of exotic many-body phenomena such as quantum many-body scars, and open up new opportunities for quantum metrology with enhanced sensitivities. These results enable wide-reaching new applications for dynamical decoupling and Hamiltonian engineering in many-body physics and quantum metrology. Published by the American Physical Society 2024