Search

Quantum many-body scarred systems host nonthermal excited eigenstates immersed in a sea of thermal ones. In cases where exact expressions for these special eigenstates are not known, it is computationally demanding to distinguish them from their exponentially many thermal neighbors. We propose a matrix-product-state (MPS) algorithm, dubbed DMRG-S, to extract such states at system sizes far beyond the scope of exact diagonalization. Using this technique, we obtain scarred eigenstates in Rydberg-blockaded chains of up to 80 sites and perform a finite-size scaling study to address the lingering question of the stability for the Néel state revivals in the thermodynamic limit. Our method also provides a systematic way to obtain exact MPS representations for scarred eigenstates near the target energy without a priori knowledge. In particular, we find several new scarred eigenstates with exact MPS representations in kinetically constrained spin and clock models. The combination of numerical and analytical investigations in our work provides a new methodology for future studies of quantum many-body scars.

In this Letter, we study quantum many-body teleportation, where a single qubit is teleported through a strongly interacting quantum system, as a result of a scrambling unitary and local measurements on a few qubits. Usual many-body teleportation protocols require a double copy of the system, and backward time evolution, we demonstrate that teleportation is possible in the 2D spin-1/2 XY model, without these constraints. The necessary long-range entanglement for teleportation is generated from the model hosting special eigenstates known as rainbow scars. We outline a specific protocol for preparing this highly entangled state by evolving a product state and performing iterative measurements on only two qubits with feedback control.

We argue that chaotic power-law interacting systems have emergent limits on information propagation, analogous to relativistic light cones, which depend on the spatial dimension d and the exponent α governing the decay of interactions. Using the dephasing nature of quantum chaos, we map the problem to a stochastic model with a known phase diagram. A linear light cone results for α≥d+1/2. We also provide a Lévy flight (long-range random walk) interpretation of the results and show consistent numerical data for 1D long-range spin models with 200 sites.

We study quantum information scrambling, specifically the growth of Heisenberg operators, in large disordered spin chains using matrix product operator dynamics to scan across the thermalization-localization quantum phase transition. We observe ballistic operator growth for weak disorder, and a sharp transition to a phase with subballistic operator spreading. The critical disorder strength for the ballistic to subballistic transition is well below the many body localization phase transition, as determined from finite size scaling of energy eigenstate entanglement entropy in small chains. In contrast, we find that the transition from subballistic to logarithmic behavior at the actual eigenstate localization transition is not resolved in our finite numerics. These data are discussed in the context of a universal form for the growing operator shape and substantiated with a simple phenomenological model of rare regions.

The interaction effects in ultracold Fermi gases with SU(N) symmetry are studied nonperturbatively in half filled one-dimensional lattices by employing quantum Monte Carlo simulations. We find that, as N increases, weak and strong interacting systems are driven to a crossover region, but from opposite directions as a convergence of itinerancy and Mottness. In the weak interaction region, particles are nearly itinerant, and interparticle collisions are enhanced by N, resulting in the amplification of interaction effects. In contrast, in the strong coupling region, increasing N softens the Mott-insulating background through the enhanced virtual hopping processes. The crossover region exhibits nearly N-independent physical quantities, including the relative bandwidth, Fermi distribution, and the spin structure factor. The difference between even-N and odd-N systems is most prominent at small N's with strong interactions, since the odd case allows local real hopping with an energy scale much larger than the virtual one. The above effects can be experimentally tested in ultracold atom experiments with alkaline-earth(-like) fermions such as ^{87}Sr (^{173}Yb).

Crystal structures and the Bloch theorem play a fundamental role in condensed matter physics. We extend the static crystal to the dynamic "space-time" crystal characterized by the general intertwined space-time periodicities in D+1 dimensions, which include both the static crystal and the Floquet crystal as special cases. A new group structure dubbed a "space-time" group is constructed to describe the discrete symmetries of a space-time crystal. Compared to space and magnetic groups, the space-time group is augmented by "time-screw" rotations and "time-glide" reflections involving fractional translations along the time direction. A complete classification of the 13 space-time groups in one-plus-one dimensions (1+1D) is performed. The Kramers-type degeneracy can arise from the glide time-reversal symmetry without the half-integer spinor structure, which constrains the winding number patterns of spectral dispersions. In 2+1D, nonsymmorphic space-time symmetries enforce spectral degeneracies, leading to protected Floquet semimetal states. We provide a general framework for further studying topological properties of the (D+1)-dimensional space-time crystal.