QUANTUM GASES. Observation of many-body localization of interacting fermions in a quasirandom optical lattice

Quantum Walks and Correlated Dynamics in an Interacting Synthetic Rydberg Lattice

Coherent dynamics of interacting quantum particles plays a central role in the study of strongly correlated quantum matter and the pursuit of quantum information processors. Here, we present the state space of interacting Rydberg atoms as a synthetic landscape on which to control and observe coherent and correlated dynamics. With full control of the coupling strengths and energy offsets between the pairs of sites in a nine-site synthetic lattice, we realize quantum walks, Bloch oscillations, and dynamics in an Escher-type "continuous staircase." In the interacting regime, we observe correlated quantum walks, Bloch oscillations, and confinement of particle pairs. Additionally, we simultaneously tilt our lattice both up and down to achieve coherent pair oscillations. When combined with a few straightforward upgrades, this work establishes synthetic Rydberg lattices of interacting atom arrays as a promising platform for programmable quantum many-body dynamics with access to features that are difficult to realize in real-space lattices.

Ultraslow Growth of Number Entropy in an ℓ-Bit Model of Many-Body Localization

We demonstrate that slow growth of the number entropy following a quench from a local product state is consistent with many-body localization. To do this, we construct a novel random circuit ℓ-bit model with exponentially localized ℓ-bits and exponentially decaying interactions between them. We observe an ultraslow growth of the number entropy starting from a Néel state, saturating at a value that grows with system size. This suggests that the observation of such growth in microscopic models is not sufficient to rule out many-body localization.

Interaction-Driven Instabilities in the Random-Field XXZ Chain

Despite enormous efforts devoted to the study of the many-body localization (MBL) phenomenon, the nature of the high-energy behavior of the Heisenberg spin chain in a strong random magnetic field is lacking consensus. Here, we take a step back by exploring the weak interaction limit starting from the Anderson localized (AL) insulator. Through shift-invert diagonalization, we find that, below a certain disorder threshold h^{*}, weak interactions necessarily lead to an ergodic instability, whereas at strong disorder the AL insulator directly turns into MBL, in agreement with a simple interpretation of the avalanche theory for restoration of ergodicity. We further map the phase diagram for the generic XXZ model in the disorder h-interaction Δ plane. Taking advantage of the total magnetization conservation, our results unveil the remarkable behavior of the spin-spin correlation functions: in the regime indicated as MBL by standard observables, their exponential decay undergoes an inversion of orientation ξ_{z}>ξ_{x}. We find that the longitudinal length ξ_{z} is a key quantity for capturing ergodic instabilities, as it increases with system size near the thermal phase, in sharp contrast to its transverse counterpart.

Interaction Effects on the Dynamical Anderson Metal-Insulator Transition Using Kicked Quantum Gases

Understanding the interplay of interaction and disorder in quantum transport poses long-standing scientific challenges for theory and experiment. While highly controlled ultracold atomic platforms combining atomic interactions with spatially disordered lattices have led to remarkable advances, the extension of such controlled studies to phenomena in high-dimensional disordered systems, such as the three-dimensional Anderson metal-insulator transition has been limited. Kicked quantum gases provide an alternate experimental platform that captures the Anderson model in momentum space and features dynamical localization as the analog of Anderson localization. Here, we utilize a momentum space lattice platform using quasiperiodically kicked ultracold atomic gases to experimentally investigate interaction effects on the three-dimensional dynamical Anderson metal-insulator transition. We observe interaction-driven subdiffusion and a divergence of delocalization onset time on approaching the phase boundary. Mean-field numerical simulations show qualitative agreement with experimental observations, but with significant quantitative deviations.

Survival Probability, Particle Imbalance, and Their Relationship in Quadratic Models

We argue that the dynamics of particle imbalance in quadratic fermionic models is, for the majority of initial many-body product states in the site occupation basis, virtually indistinguishable from the dynamics of survival probabilities of single-particle states. We then generalize our statement to a similar relationship between the non-equal time and space density correlation functions in many-body states, and the transition probabilities of single-particle states at nonzero distances. Finally, we study the equal-time connected density-density correlation functions in many-body states, which exhibit certain qualitative analogies with the survival and transition probabilities of single-particle states. Our results are numerically tested for two paradigmatic models of single-particle localization: the 3D Anderson model and the 1D Aubry-André model. This work gives an affirmative answer to the question of whether it is possible to measure features of single-particle survival and transition probabilities by the dynamics of observables in many-body states.

Superconducting Quantum Simulation for Many-Body Physics beyond Equilibrium

Quantum computing is an exciting field that uses quantum principles, such as quantum superposition and entanglement, to tackle complex computational problems. Superconducting quantum circuits, based on Josephson junctions, is one of the most promising physical realizations to achieve the long-term goal of building fault-tolerant quantum computers. The past decade has witnessed the rapid development of this field, where many intermediate-scale multi-qubit experiments emerged to simulate nonequilibrium quantum many-body dynamics that are challenging for classical computers. Here, we review the basic concepts of superconducting quantum simulation and their recent experimental progress in exploring exotic nonequilibrium quantum phenomena emerging in strongly interacting many-body systems, e.g., many-body localization, quantum many-body scars, and discrete time crystals. We further discuss the prospects of quantum simulation experiments to truly solve open problems in nonequilibrium many-body systems.

Unveiling the Stabilizer Group of a Matrix Product State

We present a novel classical algorithm designed to learn the stabilizer group-namely, the group of Pauli strings for which a state is a ±1 eigenvector-of a given matrix product state (MPS). The algorithm is based on a clever and theoretically grounded biased sampling in the Pauli (or Bell) basis. Its output is a set of independent stabilizer generators whose total number is directly associated with the stabilizer nullity, notably a well-established nonstabilizer monotone. We benchmark our method on T-doped states randomly scrambled via Clifford unitary dynamics, demonstrating very accurate estimates up to highly entangled MPS with bond dimension χ∼10^{3}. Our method, thanks to a very favorable scaling O(χ^{3}), represents the first effective approach to obtain a genuine magic monotone for MPS, enabling systematic investigations of quantum many-body physics out of equilibrium.

Probing the design principles of photosynthetic systems through fluorescence noise measurement

Elucidating the energetic processes which govern photosynthesis, the engine of life on earth, are an essential goal both for fundamental research and for cutting-edge biotechnological applications. Fluorescent signal of photosynthetic markers has long been utilised in this endeavour. In this research we demonstrate the use of fluorescent noise analysis to reveal further layers of intricacy in photosynthetic energy transfer. While noise is a common tool analysing dynamics in physics and engineering, its application in biology has thus far been limited. Here, a distinct behaviour in photosynthetic pigments across various chemical and biological environments is measured. These changes seem to elucidate quantum effects governing the generation of oxidative radicals. Although our method offers insights, it is important to note that the interpretation should be further validated expertly to support as conclusive theory. This innovative method is simple, non-invasive, and immediate, making it a promising tool to uncover further, more complex energetic events in photosynthesis, with potential uses in environmental monitoring, agriculture, and food-tech.

Scaling and universality at ramped quench dynamical quantum phase transitions

The nonequilibrium dynamics of a periodically driven extended XY model, in the presence of linear time dependent magnetic field, is investigated using the notion of dynamical quantum phase transitions (DQPTs). Along the similar lines to the equilibrium phase transition, the main purpose of this work is to search fundamental concepts such as scaling and universality at the ramped quench DQPTs. We have shown that the critical points of the model, where the gap closing occurs, can be moved by tuning the driven frequency and consequently the presence of or absence of DQPTs can be flexibly controlled by adjusting the driven frequency. We have uncovered that, for a ramp across the single quantum critical point, the critical mode at which DQPTs occur is classified into three regions: the Kibble-Zurek (KZ) region, where the critical mode scales linearly with the square root of the sweep velocity, the pre-saturated (PS) region, and the saturated (S) region where the critical mode makes a plateau versus the sweep velocity. While for a ramp that crosses two critical points, the critical modes disclose just the KZ and PS regions. On the basis of numerical simulations, we find that the dynamical free energy scales linearly with time, as approaches to DQPT time, with the exponentν=1±0.01for all sweep velocities and driven frequencies.

Quantum criticality of generalized Aubry-André models with exact mobility edges using fidelity susceptibility

In this study, we explore the quantum critical phenomena in generalized Aubry-André models, with a particular focus on the scaling behavior at various filling states. Our approach involves using quantum fidelity susceptibility to precisely identify the mobility edges in these systems. Through a finite-size scaling analysis of the fidelity susceptibility, we are able to determine both the correlation-length critical exponent and the dynamical critical exponent at the critical point of the generalized Aubry-André model. Based on the Diophantine equation conjecture, we can determines the number of subsequences of the Fibonacci sequence and the corresponding scaling functions for a specific filling fraction, as well as the universality class. Our findings demonstrate the effectiveness of employing the generalized fidelity susceptibility for the analysis of unconventional quantum criticality and the associated universal information of quasiperiodic systems in cutting-edge quantum simulation experiments.

Stark Many-Body Localization in Interacting Infinite Dimensional Systems

We study bulk particle transport in a Fermi-Hubbard model on an infinite-dimensional Bethe lattice, driven by a constant electric field. Previous numerical studies showed that one dimensional analogs of this system exhibit a breakdown of diffusion due to Stark many-body localization at least up to time that scales exponentially with the system size. Here, we consider systems initially in a spin density wave state using a combination of numerically exact and approximate techniques. We show that for sufficiently weak electric fields, the wave's momentum component decays exponentially with time in a way consistent with normal diffusion. By studying different wavelengths, we extract the dynamical exponent and the generalized diffusion coefficient at each field strength. Interestingly, we find a nonmonotonic dependence of the dynamical exponent on the electric field. As the field increases toward a critical value proportional to the Hubbard interaction strength, transport slows down, becoming subdiffusive. At large interaction strengths, however, transport speeds up again with increasing field, exhibiting superdiffusive characteristics when the electric field is comparable to the interaction strength. Eventually, at the large field limit, localization occurs and the current through the system is suppressed.

Dimerized Hofstadter model in two-leg ladder quasi-crystals

We theoretically study topological features, band structure, and localization properties of a dimerized two-leg ladder with an oscillating on-site potential. The periodicity of the on-site potential can take either rational or irrational values. We consider two types of dimerized configurations; symmetric and asymmetric models. For rational values of the periodicity as long as inversion symmetry is preserved both symmetric and asymmetric ladders can host topological phases. Additionally, the energy spectrum of the models exhibits a fractal structure known as the Hofstadter butterfly spectrum, dependent on the dimerization of the hopping and the strength of the on-site potential. In the case of irrational values for the periodicity, a metal-insulator phase transition occurs with small values of the critical strength of the on-site potential in the dimerized cases. Our models incorporate the effects of lattice configuration and quasi-periodicity, paving the way for establishing platforms that host both topological and non-topological phase transitions.

Observation of Subdiffusive Dynamic Scaling in a Driven and Disordered Bose Gas

We explore the dynamics of a tuneable box-trapped Bose gas under strong periodic forcing in the presence of weak disorder. In absence of interparticle interactions, the interplay of the drive and disorder results in an isotropic nonthermal momentum distribution that shows subdiffusive dynamic scaling, with sublinear energy growth and the universal scaling function captured well by a compressed exponential. We explain that this subdiffusion in momentum space can naturally be understood as a random walk in energy space. We also experimentally show that for increasing interaction strength, the gas behavior smoothly crosses over to wave turbulence characterized by a power-law momentum distribution, which opens new possibilities for systematic studies of the interplay of disorder and interactions in driven quantum systems.

Classical route to ergodicity and scarring in collective quantum systems

Ergodicity, a fundamental concept in statistical mechanics, is not yet a fully understood phenomena for closed quantum systems, particularly its connection with the underlying chaos. In this review, we consider a few examples of collective quantum systems to unveil the intricate relationship of ergodicity as well as its deviation due to quantum scarring phenomena with their classical counterpart. A comprehensive overview of classical and quantum chaos is provided, along with the tools essential for their detection. Furthermore, we survey recent theoretical and experimental advancements in the domain of ergodicity and its violations. This review aims to illuminate the classical perspective of quantum scarring phenomena in interacting quantum systems.

Quantum Scars and Regular Eigenstates in a Chaotic Spinor Condensate

Quantum many-body scars consist of a few low-entropy eigenstates in an otherwise chaotic many-body spectrum, and can weakly break ergodicity resulting in robust oscillatory dynamics. The notion of quantum many-body scars follows the original single-particle scars introduced within the context of quantum billiards, where scarring manifests in the form of a quantum eigenstate concentrating around an underlying classical unstable periodic orbit. A direct connection between these notions remains an outstanding problem. Here, we study a many-body spinor condensate that, owing to its collective interactions, is amenable to the diagnostics of scars. We characterize the system's rich dynamics, spectrum, and phase space, consisting of both regular and chaotic states. The former are low in entropy, violate the eigenstate thermalization hypothesis, and can be traced back to integrable effective Hamiltonians, whereas most of the latter are scarred by the underlying semiclassical unstable periodic orbits, while satisfying the eigenstate thermalization hypothesis. We outline an experimental proposal to probe our theory in trapped spin-1 Bose-Einstein condensates.

Coulomb Blockade in a Nonthermalized Quantum Dot

We investigate nonequilibrium transport properties of a quantum dot in the Coulomb blockade regime under the condition of negligible inelastic scattering during the dwelling time of the electrons in the dot. Using the quantum kinetic equation we show that the absence of thermalization leads to a double step in the distribution function of electrons on the dot, provided that it is symmetrically coupled to the leads. This drastically changes nonlinear transport through the dot resulting in an additional (compared to the thermalized case) jump in the conductance at voltages close to the charging energy, which could serve as an experimental manifestation of the absence of thermalization.

Observation of an exotic insulator to insulator transition upon electron doping the Mott insulator CeMnAsO

A promising route to discover exotic electronic states in correlated electron systems is to vary the hole or electron doping away from a Mott insulating state. Important examples include quantum criticality and high-temperature superconductivity in cuprates. Here, we report the surprising discovery of a quantum insulating state upon electron doping the Mott insulator CeMnAsO, which emerges below a distinct critical transition temperature, TII. The insulator-insulator transition is accompanied by a significant reduction in electron mobility as well as a colossal Seebeck effect and slow dynamics due to decoupling of the electrons from the lattice phonons. The origin of the transition is tentatively interpreted in terms of many-body localization, which has not been observed previously in a solid-state material.

Critical Phase Dualities in 1D Exactly Solvable Quasiperiodic Models

We propose a solvable class of 1D quasiperiodic tight-binding models encompassing extended, localized, and critical phases, separated by nontrivial mobility edges. Limiting cases include the Aubry-André model and the models of Sriram Ganeshan, J. H. Pixley, and S. Das Sarma [Phys. Rev. Lett. 114, 146601 (2015)PRLTAO0031-900710.1103/PhysRevLett.114.146601] and J. Biddle and S. Das Sarma [Phys. Rev. Lett. 104, 070601 (2010)PRLTAO0031-900710.1103/PhysRevLett.104.070601]. The analytical treatment follows from recognizing these models as a novel type of fixed points of the renormalization group procedure recently proposed in Phys. Rev. B 108, L100201 (2023)10.1103/PhysRevB.108.L100201 for characterizing phases of quasiperiodic structures. Beyond known limits, the proposed class of models extends previously encountered localized-delocalized duality transformations to points within multifractal critical phases. Besides an experimental confirmation of multifractal duality, realizing the proposed class of models in optical lattices allows stabilizing multifractal critical phases and nontrivial mobility edges in an undriven system without the need for the unbounded potentials required by previous proposals.

Quasicrystalline Bose Glass in the Absence of Disorder and Quasidisorder

We study the low-temperature phases of interacting bosons on a two-dimensional quasicrystalline lattice. By means of numerically exact path integral Monte Carlo simulations, we show that for sufficiently weak interactions the system is a homogeneous Bose-Einstein condensate that develops density modulations for increasing filling factor. The simultaneous occurrence of sizeable condensate fraction and density modulation can be interpreted as the analogous, in a quasicrystalline lattice, of supersolid phases occurring in conventional periodic lattices. For sufficiently large interaction strength and particle density, global condensation is lost and quantum exchanges are restricted to specific spatial regions. The emerging quantum phase is therefore a Bose glass, which here is stabilized in the absence of any source of disorder or quasidisorder, purely as a result of the interplay between quantum effects, particle interactions and quasicrystalline substrate. This finding clearly indicates that (quasi)disorder is not essential to observe Bose glass physics. Our results are of interest for ongoing experiments on (quasi)disorder-free quasicrystalline lattices.

Exact New Mobility Edges between Critical and Localized States

The disorder systems host three types of fundamental quantum states, known as the extended, localized, and critical states, of which the critical states remain being much less explored. Here we propose a class of exactly solvable models which host a novel type of exact mobility edges (MEs) separating localized states from robust critical states, and propose experimental realization. Here the robustness refers to the stability against both single-particle perturbation and interactions in the few-body regime. The exactly solvable one-dimensional models are featured by a quasiperiodic mosaic type of both hopping terms and on-site potentials. The analytic results enable us to unambiguously obtain the critical states which otherwise require arduous numerical verification including the careful finite size scalings. The critical states and new MEs are shown to be robust, illustrating a generic mechanism unveiled here that the critical states are protected by zeros of quasiperiodic hopping terms in the thermodynamic limit. Further, we propose a novel experimental scheme to realize the exactly solvable model and the new MEs in an incommensurate Rydberg Raman superarray. This Letter may pave a way to precisely explore the critical states and new ME physics with experimental feasibility.

Shadow Tomography from Emergent State Designs in Analog Quantum Simulators

We introduce a method that allows one to infer many properties of a quantum state-including nonlinear functions such as Rényi entropies-using only global control over the constituent degrees of freedom. In this protocol, the state of interest is first entangled with a set of ancillas under a fixed global unitary, before projective measurements are made. We show that when the unitary is sufficiently entangling, a universal relationship between the statistics of the measurement outcomes and properties of the state emerges, which can be connected to the recently discovered phenomeonon of emergent quantum state designs in chaotic systems. Thanks to this relationship, arbitrary observables can be reconstructed using the same number of experimental repetitions that would be required in classical shadow tomography [Huang et al., Nat. Phys. 16, 1050 (2020)NPAHAX1745-247310.1038/s41567-020-0932-7]. Unlike previous approaches to shadow tomography, our protocol can be implemented using only global Hamiltonian evolution, as opposed to qubit-selective logic gates, which makes it particularly well suited to analog quantum simulators, including ultracold atoms in optical lattices and arrays of Rydberg atoms.

Late-time universal distribution functions of observables in one-dimensional many-body quantum systems

We study the probability distribution function of the long-time values of observables being time-evolved by Hamiltonians modeling clean and disordered one-dimensional chains of many spin-1/2 particles. In particular, we analyze the return probability and its version for a completely extended initial state, the so-called spectral form factor. We complement our analysis with the spin autocorrelation and connected spin-spin correlation functions, both of interest in experiments with quantum simulators. We show that the distribution function has a universal shape provided the central limit theorem holds. Explicitly, the shape is exponential for the return probability and spectral form factor, meanwhile it is Gaussian for the few-body observables. We also discuss implications over the so-called many-body localization. Remarkably, our approach requires only a single sample of the dynamics and small system sizes, which could be quite advantageous when dealing specially with disordered systems.

Spontaneous symmetry breaking and localization in nonequilibrium steady states of interactive quantum systems

The time evolution of a physical system is generally described by a differential equation, which can be solved numerically by adopting a difference scheme with space-time discretization. This discretization, as a numerical artifact, results in accumulated errors during evolution and thus usually plays a negative role in simulations. In a quantum circuit, however, the "evolution time" is represented by the depth of the circuit layer, and thus is intrinsically discrete. Hence, the discretization-induced error therein is not a numerical artifact, but a physical observable effect responsible for remarkable nonequilibrium phenomena absent in conventional quantum dynamics. In this paper, we show that the combination of measurement feedback and temporal discretization can give rise to a new type of quantum dynamics. As physical consequences of this interactive quantum dynamics, a nonequilibrium steady state with spontaneous symmetry breaking is revealed in a zero-dimensional (single-qubit) system. A localization mechanism distinct from that in the well-established Anderson localization has also been proposed in a one-dimensional interactive quantum system.

One-dimensional Lévy quasicrystal

Space-fractional quantum mechanics (SFQM) is a generalization of the standard quantum mechanics when the Brownian trajectories in Feynman path integrals are replaced by Lévy flights. We introduce Lévy quasicrystal by discretizing the space-fractional Schrödinger equation using the Grünwald-Letnikov derivatives and adding on-site quasiperiodic potential. The discretized version of the usual Schrödinger equation maps to the Aubry-André (AA) Hamiltonian, which supports localization-delocalization transition even in one dimension. We find the similarities between Lévy quasicrystal and the AA model with power-law hopping, and show that the Lévy quasicrystal supports a delocalization-localization transition as one tunes the quasiperiodic potential strength and shows the coexistence of localized and delocalized states separated by mobility edge. Hence, a possible realization of SFQM in optical experiments should be a new experimental platform to test the predictions of AA models in the presence of power-law hopping.

Reviving product states in the disordered Heisenberg chain

When a generic quantum system is prepared in a simple initial condition, it typically equilibrates toward a state that can be described by a thermal ensemble. A known exception is localized systems that are non-ergodic and do not thermalize; however, local observables are still believed to become stationary. Here we demonstrate that this general picture is incomplete by constructing product states that feature periodic high-fidelity revivals of the full wavefunction and local observables that oscillate indefinitely. The system neither equilibrates nor thermalizes. This is analogous to the phenomenon of weak ergodicity breaking due to many-body scars and challenges aspects of the current phenomenology of many-body localization, such as the logarithmic growth of the entanglement entropy. To support our claim, we combine analytic arguments with large-scale tensor network numerics for the disordered Heisenberg chain. Our results hold for arbitrarily long times in chains of 160 sites up to machine precision.

Phenomenology of the Prethermal Many-Body Localized Regime

The dynamical phase diagram of interacting disordered systems has seen substantial revision over the past few years. Theory must now account for a large prethermal many-body localized regime in which thermalization is extremely slow, but not completely arrested. We derive a quantitative description of these dynamics in short-ranged one-dimensional systems using a model of successive many-body resonances. The model explains the decay timescale of mean autocorrelators, the functional form of the decay-a stretched exponential-and relates the value of the stretch exponent to the broad distribution of resonance timescales. The Jacobi method of matrix diagonalization provides numerical access to this distribution, as well as a conceptual framework for our analysis. The resonance model correctly predicts the stretch exponents for several models in the literature. Successive resonances may also underlie slow thermalization in strongly disordered systems in higher dimensions, or with long-range interactions.

Proposal for Observing Yang-Lee Criticality in Rydberg Atomic Arrays

Yang-Lee edge singularities (YLES) are the edges of the partition function zeros of an interacting spin model in the space of complex control parameters. They play an important role in understanding non-Hermitian phase transitions in many-body physics, as well as characterizing the corresponding nonunitary criticality. Even though such partition function zeroes have been measured in dynamical experiments where time acts as the imaginary control field, experimentally demonstrating such YLES criticality with a physical imaginary field has remained elusive due to the difficulty of physically realizing non-Hermitian many-body models. We provide a protocol for observing the YLES by detecting kinked dynamical magnetization responses due to broken PT symmetry, thus enabling the physical probing of nonunitary phase transitions in nonequilibrium settings. In particular, scaling analyses based on our nonunitary time evolution circuit with matrix product states accurately recover the exponents uniquely associated with the corresponding nonunitary CFT. We provide an explicit proposal for observing YLES criticality in Floquet quenched Rydberg atomic arrays with laser-induced loss, which paves the way towards a universal platform for simulating non-Hermitian many-body dynamical phenomena.

Ergodicity breaking in rapidly rotating C60 fullerenes

Ergodicity, the central tenet of statistical mechanics, requires an isolated system to explore all available phase space constrained by energy and symmetry. Mechanisms for violating ergodicity are of interest for probing nonequilibrium matter and protecting quantum coherence in complex systems. Polyatomic molecules have long served as a platform for probing ergodicity breaking in vibrational energy transport. Here, we report the observation of rotational ergodicity breaking in an unprecedentedly large molecule, 12C60, determined from its icosahedral rovibrational fine structure. The ergodicity breaking occurs well below the vibrational ergodicity threshold and exhibits multiple transitions between ergodic and nonergodic regimes with increasing angular momentum. These peculiar dynamics result from the molecule's distinctive combination of symmetry, size, and rigidity, highlighting its relevance to emergent phenomena in mesoscopic quantum systems.

Absence of localization in interacting spin chains with a discrete symmetry

Novel paradigms of strong ergodicity breaking have recently attracted significant attention in condensed matter physics. Understanding the exact conditions required for their emergence or breakdown not only sheds more light on thermalization and its absence in closed quantum many-body systems, but it also has potential benefits for applications in quantum information technology. A case of particular interest is many-body localization whose conditions are not yet fully settled. Here, we prove that spin chains symmetric under a combination of mirror and spin-flip symmetries and with a non-degenerate spectrum show finite spin transport at zero total magnetization and infinite temperature. We demonstrate this numerically using two prominent examples: the Stark many-body localization system (Stark-MBL) and the symmetrized many-body localization system (symmetrized-MBL). We provide evidence of delocalization at all energy densities and show that delocalization persists when the symmetry is broken. We use our results to construct two localized systems which, when coupled, delocalize each other. Our work demonstrates the dramatic effect symmetries can have on disordered systems, proves that the existence of exact resonances is not a sufficient condition for delocalization, and opens the door to generalization to higher spatial dimensions and different conservation laws.

Many-body correlations in one-dimensional optical lattices with alkaline-earth(-like) atoms

We explore the rich nature of correlations in the ground state of ultracold atoms trapped in state-dependent optical lattices. In particular, we consider interacting fermionic ytterbium or strontium atoms, realizing a two-orbital Hubbard model with two spin components. We analyze the model in one-dimensional setting with the experimentally relevant hierarchy of tunneling and interaction amplitudes by means of exact diagonalization and matrix product states approaches, and study the correlation functions in density, spin, and orbital sectors as functions of variable densities of atoms in the ground and metastable excited states. We show that in certain ranges of densities these atomic systems demonstrate strong density-wave, ferro- and antiferromagnetic, as well as antiferroorbital correlations.

Thermodynamic Phase Diagram of Two-Dimensional Bosons in a Quasicrystal Potential

Quantum simulation of quasicrystals in synthetic bosonic matter now paves the way for the exploration of these intriguing systems in wide parameter ranges. Yet thermal fluctuations in such systems compete with quantum coherence and significantly affect the zero-temperature quantum phases. Here we determine the thermodynamic phase diagram of interacting bosons in a two-dimensional, homogeneous quasicrystal potential. We find our results using quantum Monte Carlo simulations. Finite-size effects are carefully taken into account and the quantum phases are systematically distinguished from thermal phases. In particular, we demonstrate stabilization of a genuine Bose glass phase against the normal fluid in sizable parameter ranges. We interpret our results for strong interactions using a fermionization picture and discuss experimental relevance.

Controlling Local Thermalization Dynamics in a Floquet-Engineered Dipolar Ensemble

Understanding the microscopic mechanisms of thermalization in closed quantum systems is among the key challenges in modern quantum many-body physics. We demonstrate a method to probe local thermalization in a large-scale many-body system by exploiting its inherent disorder and use this to uncover the thermalization mechanisms in a three-dimensional, dipolar-interacting spin system with tunable interactions. Utilizing advanced Hamiltonian engineering techniques to explore a range of spin Hamiltonians, we observe a striking change in the characteristic shape and timescale of local correlation decay as we vary the engineered exchange anisotropy. We show that these observations originate from the system's intrinsic many-body dynamics and reveal the signatures of conservation laws within localized clusters of spins, which do not readily manifest using global probes. Our method provides an exquisite lens into the tunable nature of local thermalization dynamics and enables detailed studies of scrambling, thermalization, and hydrodynamics in strongly interacting quantum systems.

Work Extraction from Unknown Quantum Sources

Energy extraction is a central task in thermodynamics. In quantum physics, ergotropy measures the amount of work extractable under cyclic Hamiltonian control. As its full extraction requires perfect knowledge of the initial state, however, it does not characterize the work value of unknown or untrusted quantum sources. Fully characterizing such sources would require quantum tomography, which is prohibitively costly in experiments due to the exponential growth of required measurements and operational limitations. Here, we therefore derive a new notion of ergotropy applicable when nothing is known about the quantum states produced by the source, apart from what can be learned by performing only a single type of coarse-grained measurement. We find that in this case the extracted work is defined by the Boltzmann and observational entropy in cases where the measurement outcomes are, or are not, used in the work extraction, respectively. This notion of ergotropy represents a realistic measure of extractable work, which can be used as the relevant figure of merit to characterize a quantum battery.

Engineering random spin models with atoms in a high-finesse cavity

All-to-all interacting, disordered quantum many-body models have a wide range of applications across disciplines, from spin glasses in condensed-matter physics over holographic duality in high-energy physics to annealing algorithms in quantum computing. Typically, these models are abstractions that do not find unambiguous physical realizations in nature. Here we realize an all-to-all interacting, disordered spin system by subjecting an atomic cloud in a cavity to a controllable light shift. Adjusting the detuning between atom resonance and cavity mode, we can tune between disordered versions of a central-mode model and a Lipkin-Meshkov-Glick model. By spectroscopically probing the low-energy excitations of the system, we explore the competition of interactions with disorder across a broad parameter range. We show how disorder in the central-mode model breaks the strong collective coupling, making the dark-state manifold cross over to a random distribution of weakly mixed light-matter, 'grey', states. In the Lipkin-Meshkov-Glick model, the ferromagnetic finite-sized ground state evolves towards a paramagnet as disorder is increased. In that regime, semi-localized eigenstates emerge, as we observe by extracting bounds on the participation ratio. These results present substantial steps towards freely programmable cavity-mediated interactions for the design of arbitrary spin Hamiltonians.

Dynamical Gauge Fields with Bosonic Codes

The quantum simulation of dynamical gauge field theories offers the opportunity to study complex high-energy physics with controllable low-energy devices. For quantum computation, bosonic codes promise robust error correction that exploits multiparticle redundancy in bosons. Here, we demonstrate how bosonic codes can be used to simulate dynamical gauge fields. We encode both matter and dynamical gauge fields in a network of resonators that are coupled via three-wave mixing. The mapping to a Z_{2} dynamical lattice gauge theory is established when the gauge resonators operate as Schrödinger cat states. We explore the optimal conditions under which the system preserves the required gauge symmetries. Our findings promote realizing high-energy models using bosonic codes.

Exploring the Regime of Fragmentation in Strongly Tilted Fermi-Hubbard Chains

Intriguingly, quantum many-body systems may defy thermalization even without disorder. One example is so-called fragmented models, where the many-body Hilbert space fragments into dynamically disconnected subspaces that are not determined by the global symmetries of the model. In this Letter we demonstrate that the tilted one-dimensional Fermi-Hubbard model naturally realizes distinct effective Hamiltonians that are expected to support nonergodic behavior due to fragmentation, even at resonances between the tilt energy and the Hubbard on site interaction. We find that the effective description captures the observed dynamics in experimentally accessible parameter ranges of moderate tilt values. Specifically, we observe a pronounced dependence of the relaxation dynamics on the initial doublon fraction, which directly reveals the microscopic processes of the fragmented model. Our results pave the way for future studies of nonergodic behavior in higher dimensions.

Down-conversion of a single photon as a probe of many-body localization

Decay of a particle into more particles is a ubiquitous phenomenon to interacting quantum systems, taking place in colliders, nuclear reactors or solids. In a nonlinear medium, even a single photon would decay by down-converting (splitting) into lower-frequency photons with the same total energy¹, at a rate given by Fermi's golden rule. However, the energy-conservation condition cannot be matched precisely if the medium is finite and only supports quantized modes. In this case, the fate of the photon becomes the long-standing question of many-body localization, originally formulated as a gedanken experiment for the lifetime of a single Fermi-liquid quasiparticle confined to a quantum dot². Here we implement such an experiment using a superconducting multimode cavity, the nonlinearity of which was tailored to strongly violate the photon-number conservation. The resulting interaction attempts to convert a single photon excitation into a shower of low-energy photons but fails owing to the many-body localization mechanism, which manifests as a striking spectral fine structure of multiparticle resonances at the standing-wave-mode frequencies of the cavity. Each resonance was identified as a many-body state of radiation composed of photons from a broad frequency range and not obeying Fermi's golden rule theory. Our result introduces a new platform to explore the fundamentals of many-body localization without having to control many atoms or qubits^3-9.

Restoring Ergodicity in a Strongly Disordered Interacting Chain

We consider a chain of interacting fermions with random disorder that was intensively studied in the context of many-body localization. We show that only a small fraction of the two-body interaction represents a true local perturbation to the Anderson insulator. While this true perturbation is nonzero at any finite disorder strength W, it decreases with increasing W. This establishes a view that the strongly disordered system should be viewed as a weakly perturbed integrable model, i.e., a weakly perturbed Anderson insulator. As a consequence, the latter can hardly be distinguished from a strictly integrable system in finite-size calculations at large W. We then introduce a rescaled model in which the true perturbation is of the same order of magnitude as the other terms of the Hamiltonian, and show that the system remains ergodic at arbitrary large disorder.

Canonical Density Matrices from Eigenstates of Mixed Systems

One key issue of the foundation of statistical mechanics is the emergence of equilibrium ensembles in isolated and closed quantum systems. Recently, it was predicted that in the thermodynamic (N→∞) limit of large quantum many-body systems, canonical density matrices emerge for small subsystems from almost all pure states. This notion of canonical typicality is assumed to originate from the entanglement between subsystem and environment and the resulting intrinsic quantum complexity of the many-body state. For individual eigenstates, it has been shown that local observables show thermal properties provided the eigenstate thermalization hypothesis holds, which requires the system to be quantum-chaotic. In the present paper, we study the emergence of thermal states in the regime of a quantum analog of a mixed phase space. Specifically, we study the emergence of the canonical density matrix of an impurity upon reduction from isolated energy eigenstates of a large but finite quantum system the impurity is embedded in. Our system can be tuned by means of a single parameter from quantum integrability to quantum chaos and corresponds in between to a system with mixed quantum phase space. We show that the probability for finding a canonical density matrix when reducing the ensemble of energy eigenstates of the finite many-body system can be quantitatively controlled and tuned by the degree of quantum chaos present. For the transition from quantum integrability to quantum chaos, we find a continuous and universal (i.e., size-independent) relation between the fraction of canonical eigenstates and the degree of chaoticity as measured by the Brody parameter or the Shannon entropy.

Dynamics in Systems with Modulated Symmetries

We extend the notions of multipole and subsystem symmetries to more general spatially modulated symmetries. We uncover two instances with exponential and (quasi)periodic modulations and provide simple microscopic models in one, two, and three dimensions. Seeking to understand their effect on the long-time dynamics, we numerically study a stochastic cellular automaton evolution that obeys such symmetries. We prove that, in one dimension, the periodically modulated symmetries lead to a diffusive scaling of correlations modulated by a finite microscopic momentum. In higher dimensions, these symmetries take the form of lines and surfaces of conserved momenta. These give rise to exotic forms of subdiffusive behavior with a rich spatial structure influenced by lattice-scale features. Exponential modulation, on the other hand, can lead to correlations that are infinitely long-lived at the boundary while decaying exponentially in the bulk.

Bath Engineering Enhanced Quantum Critical Engines

Driving a quantum system across quantum critical points leads to non-adiabatic excitations in the system. This in turn may adversely affect the functioning of a quantum machine which uses a quantum critical substance as its working medium. Here we propose a bath-engineered quantum engine (BEQE), in which we use the Kibble-Zurek mechanism and critical scaling laws to formulate a protocol for enhancing the performance of finite-time quantum engines operating close to quantum phase transitions. In the case of free fermionic systems, BEQE enables finite-time engines to outperform engines operating in the presence of shortcuts to adiabaticity, and even infinite-time engines under suitable conditions, thus showing the remarkable advantages offered by this technique. Open questions remain regarding the use of BEQE based on non-integrable models.

Topological Phase Transitions and Mobility Edges in Non-Hermitian Quasicrystals

Non-Hermiticity significantly enriches the properties of topological models, leading to exotic features such as the non-Hermitian skin effects and non-Bloch bulk-boundary correspondence that have no counterparts in Hermitian settings. Its impact is particularly illustrating in non-Hermitian quasicrystals where the interplay between non-Hermiticity and quasiperiodicity results in the concurrence of the delocalization-localization transition, the parity-time (PT)-symmetry breaking, and the onset of the non-Hermitian skin effects. Here, we experimentally simulate non-Hermitian quasicrystals using photonic quantum walks. Using dynamic observables, we demonstrate that the system can transit from a delocalized, PT-symmetry broken phase that features non-Hermitian skin effects, to a localized, PT-symmetry unbroken phase with no non-Hermitian skin effects. The measured critical point is consistent with the theoretical prediction through a spectral winding number, confirming the topological origin of the phase transition. More interestingly, we also provide the first experimental evidence for mobility edges which are induced by non-Hermiticity. Our Letter opens the avenue of investigating the interplay of non-Hermiticity, quasiperiodicity, and spectral topology in open quantum systems.

Emergence of Hilbert Space Fragmentation in Ising Models with a Weak Transverse Field

The transverse-field Ising model is one of the fundamental models in quantum many-body systems, yet a full understanding of its dynamics remains elusive in higher than one dimension. Here, we show for the first time the breakdown of ergodicity in d-dimensional Ising models with a weak transverse field in a prethermal regime. We demonstrate that novel Hilbert-space fragmentation occurs in the effective nonintegrable model with d≥2 as a consequence of only one emergent global conservation law of the domain wall number. Our results indicate nontrivial initial-state dependence for nonequilibrium dynamics of the Ising models with a weak transverse field.

Localization Detection Based on Quantum Dynamics

Detecting many-body localization (MBL) typically requires the calculation of high-energy eigenstates using numerical approaches. This study investigates methods that assume the use of a quantum device to detect disorder-induced localization. Numerical simulations for small systems demonstrate how the magnetization and twist overlap, which can be easily obtained from the measurement of qubits in a quantum device, changing from the thermal phase to the localized phase. The twist overlap evaluated using the wave function at the end of the time evolution behaves similarly to the one evaluated with eigenstates in the middle of the energy spectrum under a specific condition. The twist overlap evaluated using the wave function after time evolution for many disorder realizations is a promising probe for detecting MBL in quantum computing approaches.

Ergodicity Breaking Transition in Zero Dimensions

It is of great current interest to establish toy models of ergodicity breaking transitions in quantum many-body systems. Here, we study a model that is expected to exhibit an ergodic to nonergodic transition in the thermodynamic limit upon tuning the coupling between an ergodic quantum dot and distant particles with spin-1/2. The model is effectively zero dimensional; however, a variant of the model was proposed by De Roeck and Huveneers to describe the avalanche mechanism of ergodicity breaking transition in one-dimensional disordered spin chains. We show that exact numerical results based on the spectral form factor calculation accurately agree with theoretical predictions, and hence unambiguously confirm existence of the ergodicity breaking transition in this model. We benchmark specific properties that represent hallmarks of the ergodicity breaking transition in finite systems.

Digital quantum simulation of Floquet symmetry-protected topological phases

Quantum many-body systems away from equilibrium host a rich variety of exotic phenomena that are forbidden by equilibrium thermodynamics. A prominent example is that of discrete time crystals1-8, in which time-translational symmetry is spontaneously broken in periodically driven systems. Pioneering experiments have observed signatures of time crystalline phases with trapped ions9,10, solid-state spin systems11-15, ultracold atoms16,17 and superconducting qubits18-20. Here we report the observation of a distinct type of non-equilibrium state of matter, Floquet symmetry-protected topological phases, which are implemented through digital quantum simulation with an array of programmable superconducting qubits. We observe robust long-lived temporal correlations and subharmonic temporal response for the edge spins over up to 40 driving cycles using a circuit of depth exceeding 240 and acting on 26 qubits. We demonstrate that the subharmonic response is independent of the initial state, and experimentally map out a phase boundary between the Floquet symmetry-protected topological and thermal phases. Our results establish a versatile digital simulation approach to exploring exotic non-equilibrium phases of matter with current noisy intermediate-scale quantum processors21.

Ultracold Bose Gases in Dynamic Disorder with Tunable Correlation Time

We study experimentally the dissipative dynamics of ultracold bosonic gases in a dynamic disorder potential with tunable correlation time. First, we measure the heating rate of thermal clouds exposed to the dynamic potential and present a model of the heating process, revealing the microscopic origin of dissipation from a thermal, trapped cloud of bosons. Second, for Bose-Einstein condensates, we measure the particle loss rate induced by the dynamic environment. Depending on the correlation time, the losses are either dominated by heating of residual thermal particles or the creation of excitations in the superfluid, a notion we substantiate with a rate model. Our results illuminate the interplay between superfluidity and time-dependent disorder and on more general grounds establish ultracold atoms as a platform for studying spatiotemporal noise and time-dependent disorder.

Observing emergent hydrodynamics in a long-range quantum magnet

Identifying universal properties of nonequilibrium quantum states is a major challenge in modern physics. A fascinating prediction is that classical hydrodynamics emerges universally in the evolution of any interacting quantum system. We experimentally probed the quantum dynamics of 51 individually controlled ions, realizing a long-range interacting spin chain. By measuring space-time-resolved correlation functions in an infinite temperature state, we observed a whole family of hydrodynamic universality classes, ranging from normal diffusion to anomalous superdiffusion, that are described by Lévy flights. We extracted the transport coefficients of the hydrodynamic theory, reflecting the microscopic properties of the system. Our observations demonstrate the potential for engineered quantum systems to provide key insights into universal properties of nonequilibrium states of quantum matter.

Fermionic Many-Body Localization for Random and Quasiperiodic Systems in the Presence of Short- and Long-Range Interactions

We study many-body localization (MBL) for interacting one-dimensional lattice fermions in random (Anderson) and quasiperiodic (Aubry-Andre) models, focusing on the role of interaction range. We obtain the MBL quantum phase diagrams by calculating the experimentally relevant inverse participation ratio (IPR) at half-filling using exact diagonalization methods and extrapolating to the infinite system size. For short-range interactions, our results produce in the phase diagram a qualitative symmetry between weak and strong interaction limits. For long-range interactions, no such symmetry exists as the strongly interacting system is always many-body localized, independent of the effective disorder strength, and the system is analogous to a pinned Wigner crystal. We obtain various scaling exponents for the IPR, suggesting conditions for different MBL regimes arising from interaction effects.

Quantum transitions, ergodicity, and quantum scars in the coupled top model

We consider an interacting collective spin model known as coupled top (CT), exhibiting a rich variety of phenomena related to quantum transitions, ergodicity, and formation of quantum scars, discussed in our previous work [Mondal, Sinha, and Sinha, Phys. Rev. E 102, 020101(R) (2020)2470-004510.1103/PhysRevE.102.020101]. In this work, we present a detailed analysis of the different type of transitions in the CT model, and find their connection with the underlying collective spin dynamics. Apart from the quantum scarring phenomena, we also identify another source of deviation from ergodicity due to the presence of nonergodic multifractal states. The degree of ergodicity of the eigenstates across the energy band is quantified from the relative entanglement entropy as well as multifractal dimensions, which can be probed from nonequilibrium dynamics. Finally, we discuss the detection of nonergodic behavior and different types of quantum scars using "out-of-time-order correlators," which has relevance in the recent experiments.