Kinetically constrained models were originally introduced to capture slow relaxation in glassy systems, where dynamics are hindered by local constraints instead of energy barriers. Their quantum counterparts have recently drawn attention for exhibiting highly degenerate eigenstates at zero energy-known as zero modes-stemming from chiral symmetry. Yet, the structure and implications of these zero modes remain poorly understood. In this work, we focus on the properties of the zero mode subspace in quantum kinetically constrained models with a U(1) particle-conservation symmetry. We use the U(1) East, which lacks inversion symmetry, and the inversion-symmetric U(1) East-West models to illustrate our two main results. First, we observe that the simultaneous presence of constraints and chiral symmetry generally leads to a parametric increase in the number of zero modes due to the fragmentation of the many-body Hilbert space into disconnected sectors. Second, we generalize the concept of compact localized states from single-particle physics and introduce the notion of collective bound states, a special kind of nonergodic eigenstates that are robust to enlarging the system size. We formulate sufficient criteria for their existence, arguing that the degenerate zero mode subspace plays a central role, and demonstrate bound states in both example models and in a two-dimensional model, the U(1) North-East, and in the pairflip model, a system without particle conservation. Our results motivate a systematic study of bound states and their relation to ergodicity breaking, transport, and other properties of quantum kinetically constrained models.

We present and test a protocol to learn the matrix-product operator (MPO) representation of an experimentally prepared quantum state. The protocol takes as input classical shadows corresponding to local randomized measurements, and outputs the tensors of an MPO maximizing a suitably defined fidelity with the experimental state. The tensor optimization is carried out sequentially, similarly to the well-known density matrix renormalization group algorithm. Our approach is provably efficient under certain technical conditions expected to be met in short-range correlated states and in typical noisy experimental settings. Under the same conditions, we also provide an efficient scheme to estimate fidelities between the learned and the experimental states. We experimentally demonstrate our protocol by learning entangled quantum states of up to N = 96 qubits in a superconducting quantum processor. Our method upgrades classical shadows to large-scale quantum computation and simulation experiments.

Describing general quantum many-body dynamics is a challenging task due to the exponential growth of the Hilbert space with system size. The time-dependent variational principle (TDVP) provides a powerful tool to tackle this task by projecting quantum evolution onto a classical dynamical system within a variational manifold. In classical systems, periodic orbits play a crucial role in understanding the structure of the phase space and the long-term behavior of the system. However, finding periodic orbits is generally difficult, and their existence and properties in generic TDVP dynamics over matrix product states have remained largely unexplored. In this work, we develop an algorithm to systematically identify and characterize periodic orbits in TDVP dynamics. Applying our method to the periodically kicked Ising model, we uncover both stable and unstable periodic orbits. We characterize the Kolmogorov-ArnoldMoser tori in the vicinity of stable periodic orbits and track the change of the periodic orbits as we modify the Hamiltonian parameters. We observe that periodic orbits exist at any value of the coupling constant of the kicked Ising model between prethermal and fully thermalizing regimes, but their relevance to quantum dynamics and imprint on quantum eigenstates diminishes as the system leaves the prethermal regime. Our results demonstrate that periodic orbits provide valuable insights into the TDVP approximation of quantum many-body evolution and establish a closer connection between quantum and classical chaos.

Non-Hermitian many-body localization (NH MBL) has emerged as a possible scenario for stable localization in open systems, as suggested by spectral indicators identifying a putative transition for finite system sizes. In this work, we shift the focus to dynamical probes, specifically the steady-state spin current, to investigate transport properties in a disordered, non-Hermitian XXZ spin chain. Through exact diagonalization for small systems and tensor-network methods for larger chains, we demonstrate that the steady-state current remains finite and decays exponentially with disorder strength, showing no evidence of a transition up to disorder values far beyond the previously claimed critical point. Our results reveal a stark discrepancy between spectral indicators, which suggest localization, and transport behavior, which indicates delocalization. This highlights the importance of dynamical observables in characterizing NH MBL and suggests that traditional spectral measures may not fully capture the physics of non-Hermitian systems. Additionally, we observe a noncommutativity of limits in system size and time, further complicating the interpretation of finite-size studies. These findings challenge the existence of NH MBL in the studied model and underscore the need for alternative approaches to understanding localization in non-Hermitian settings.

Eigenstates of quantum many-body systems are often used to define phases of matter in and out of equilibrium,. however, experimentally accessing highly excited eigenstates is a challenging task, calling for alternative strategies to dynamically probe nonequilibrium phases. In this work, we characterize the dynamical properties of a disordered spin chain, focusing on the spin-glass regime. Using tensor-network simulations, we observe oscillatory behavior of local expectation values and bipartite entanglement entropy. We explain these oscillations deep in the many-body localized spin-glass regime via a simple theoretical model. From perturbation theory, we predict the timescales up to which our analytical description is valid and confirm it with numerical simulations. Finally, we study the correlation length dynamics, which, after a long-time plateau, resume growing in line with renormalization group (RG) expectations. Our work suggests that RG predictions can be quantitatively tested against numerical simulations and experiments, potentially enabling microscopic descriptions of dynamical phases in large systems.

Persistent revivals recently observed in Rydberg atom simulators have challenged our understanding of thermalization and attracted much interest to the concept of quantum many-body scars (QMBSs). QMBSs are non-thermal highly excited eigenstates that coexist with typical eigenstates in the spectrum of manybody Hamiltonians, and have since been reported in multiple theoretical models, including the so-called PXP model, approximately realized by Rydberg simulators. At the same time, questions of how common QMBSs are and in what models they are physically realized remain open. In this Letter, we demonstrate that QMBSs exist in a broader family of models that includes and generalizes PXP to longer-range constraints and states with different periodicity. We show that in each model, multiple QMBS families can be found. Each of them relies on a different approximate u(2) algebra, leading to oscillatory dynamics in all cases. However, in contrast to the PXP model, their observation requires launching dynamics from weakly entangled initial states rather than from a product state. QMBSs reported here may be experimentally probed using Rydberg atom simulator in the regime of longer-range Rydberg blockades.

Estimating global properties of many-body quantum systems such as entropy or bipartite entanglement is a notoriously difficult task, typically requiring a number of measurements or classical postprocessing resources growing exponentially in the system size. In this work, we address the problem of estimating global entropies and mixed-state entanglement via partial-transposed (PT) moments and show that efficient estimation strategies exist under the assumption that all the spatial correlation lengths are finite. Focusing on one-dimensional systems, we identify a set of approximate factorization conditions (AFCs) on the system density matrix, which allow us to reconstruct entropies and PT moments from information on local subsystems. This identification yields a simple and efficient strategy for entropy and entanglement estimation. Our method could be implemented in different ways, depending on how information on local subsystems is extracted. Focusing on randomized measurements providing a practical and common measurement scheme, we prove that our protocol requires only polynomially many measurements and postprocessing operations, assuming that the state to be measured satisfies the AFCs. We prove that the AFCs hold for finite-depth quantum-circuit states and translation-invariant matrix-product density operators and provide numerical evidence that they are satisfied in more general, physically interesting cases, including thermal states of local Hamiltonians. We argue that our method could be practically useful to detect bipartite mixed-state entanglement for large numbers of qubits available in today's quantum platforms.