Probing site-resolved correlations in a spin system of ultracold molecules.

Synthetic quantum systems with interacting constituents play an important role in quantum information processing and in explaining fundamental phenomena in many-body physics. Following impressive advances in cooling and trapping techniques, ensembles of ultracold polar molecules have emerged as a promising platform that combines several advantageous properties1-11. These include a large set of internal states with long coherence times12-17 and long-range, anisotropic interactions. These features could enable the exploration of intriguing phases of correlated quantum matter, such as topological superfluids18, quantum spin liquids19, fractional Chern insulators20 and quantum magnets21,22. Probing correlations in these phases is crucial to understanding their properties, necessitating the development of new experimental techniques. Here we use quantum gas microscopy23 to measure the site-resolved dynamics of quantum correlations of polar 23Na87Rb molecules confined in a two-dimensional optical lattice. By using two rotational states of the molecules, we realize a spin-1/2 system with dipolar interactions between particles, producing a quantum spin-exchange model21,22,24,25. We study the evolution of correlations during the thermalization process of an out-of-equilibrium spin system for both spatially isotropic and anisotropic interactions. Furthermore, we examine the correlation dynamics of a spin-anisotropic Heisenberg model engineered from the native spin-exchange model by using periodic microwave pulses26-28. These experiments push the frontier of probing and controlling interacting systems of ultracold molecules, with prospects for exploring new regimes of quantum matter and characterizing entangled states that are useful for quantum computation29,30 and metrology31.

Hole statistics of equilibrium 2D and 3D hard-sphere crystals.

The probability of finding a spherical "hole" of a given radius r contains crucial structural information about many-body systems. Such hole statistics, including the void conditional nearest-neighbor probability functions GV(r), have been well studied for hard-sphere fluids in d-dimensional Euclidean space Rd. However, little is known about these functions for hard-sphere crystals for values of r beyond the hard-sphere diameter, as large holes are extremely rare in crystal phases. To overcome these computational challenges, we introduce a biased-sampling scheme that accurately determines hole statistics for equilibrium hard spheres on ranges of r that far extend those that could be previously explored. We discover that GV(r) in crystal and hexatic states exhibits oscillations whose amplitudes increase rapidly with the packing fraction, which stands in contrast to GV(r) that monotonically increases with r for fluid states. The oscillations in GV(r) for 2D crystals are strongly correlated with the local orientational order metric in the vicinity of the holes, and variations in GV(r) for 3D states indicate a transition between tetrahedral and octahedral holes, demonstrating the power of GV(r) as a probe of local coordination geometry. To further study the statistics of interparticle spacing in hard-sphere systems, we compute the local packing fraction distribution f(ϕl) of Delaunay cells and find that, for d ≤ 3, the excess kurtosis of f(ϕl) switches sign at a certain transitional global packing fraction. Our accurate methods to access hole statistics in hard-sphere crystals at the challenging intermediate length scales reported here can be applied to understand the important problem of solvation and hydrophobicity in water at such length scales.

Many-Body Resonances in the Avalanche Instability of Many-Body Localization.

Many-body localized (MBL) systems fail to reach thermal equilibrium under their own dynamics, even though they are interacting, nonintegrable, and in an extensively excited state. One instability toward thermalization of MBL systems is the so-called "avalanche," where a locally thermalizing rare region is able to spread thermalization through the full system. The spreading of the avalanche may be modeled and numerically studied in finite one-dimensional MBL systems by weakly coupling an infinite-temperature bath to one end of the system. We find that the avalanche spreads primarily via strong many-body resonances between rare near-resonant eigenstates of the closed system. Thus we find and explore a detailed connection between many-body resonances and avalanches in MBL systems.

Many-Body Quantum Chaos and Emergence of Ginibre Ensemble.

We show that non-Hermitian Ginibre random matrix behaviors emerge in spatially extended many-body quantum chaotic systems in the space direction, just as Hermitian random matrix behaviors emerge in chaotic systems in the time direction. Starting with translational invariant models, which can be associated with dual transfer matrices with complex-valued spectra, we show that the linear ramp of the spectral form factor necessitates that the dual spectra have nontrivial correlations, which in fact fall under the universality class of the Ginibre ensemble, demonstrated by computing the level spacing distribution and the dissipative spectral form factor. As a result of this connection, the exact spectral form factor for the Ginibre ensemble can be used to universally describe the spectral form factor for translational invariant many-body quantum chaotic systems in the scaling limit where t and L are large, while the ratio between L and L_{Th}, the many-body Thouless length is fixed. With appropriate variations of Ginibre models, we analytically demonstrate that our claim generalizes to models without translational invariance as well. The emergence of the Ginibre ensemble is a genuine consequence of the strongly interacting and spatially extended nature of the quantum chaotic systems we consider, unlike the traditional emergence of Hermitian random matrix ensembles.

Emergence of Fermi's Golden Rule.

Fermi's golden rule applies in the limit where an initial quantum state is weakly coupled to a continuum of other final states overlapping its energy. Here we investigate what happens away from this limit, where the set of final states is discrete, with a nonzero mean level spacing; this question arises in a number of recently investigated many-body systems. For different symmetry classes, we analytically and/or numerically calculate the universal crossovers in the average decay of the initial state as the level spacing is varied, with the golden rule emerging in the limit of a continuum. Among the corrections to the exponential decay of the initial state given by Fermi's golden rule is the appearance of the spectral form factor in the longtime regime for small but nonzero level spacing.

Field Theory of Charge Sharpening in Symmetric Monitored Quantum Circuits.

Monitored quantum circuits (MRCs) exhibit a measurement-induced phase transition between area-law and volume-law entanglement scaling. MRCs with a conserved charge additionally exhibit two distinct volume-law entangled phases that cannot be characterized by equilibrium notions of symmetry-breaking or topological order, but rather by the nonequilibrium dynamics and steady-state distribution of charge fluctuations. These include a charge-fuzzy phase in which charge information is rapidly scrambled leading to slowly decaying spatial fluctuations of charge in the steady state, and a charge-sharp phase in which measurements collapse quantum fluctuations of charge without destroying the volume-law entanglement of neutral degrees of freedom. By taking a continuous-time, weak-measurement limit, we construct a controlled replica field theory description of these phases and their intervening charge-sharpening transition in one spatial dimension. We find that the charge fuzzy phase is a critical phase with continuously evolving critical exponents that terminates in a modified Kosterlitz-Thouless transition to the short-range correlated charge-sharp phase. We numerically corroborate these scaling predictions also hold for discrete-time projective-measurement circuit models using large-scale matrix-product state simulations, and discuss generalizations to higher dimensions.

Operator Scaling Dimensions and Multifractality at Measurement-Induced Transitions.

Repeated local measurements of quantum many-body systems can induce a phase transition in their entanglement structure. These measurement-induced phase transitions (MIPTs) have been studied for various types of dynamics, yet most cases yield quantitatively similar critical exponents, making it unclear how many distinct universality classes are present. Here, we probe the properties of the conformal field theories governing these MIPTs using a numerical transfer-matrix method, which allows us to extract the effective central charge, as well as the first few low-lying scaling dimensions of operators at these critical points for (1+1)-dimensional systems. Our results provide convincing evidence that the generic and Clifford MIPTs for qubits lie in different universality classes and that both are distinct from the percolation transition for qudits in the limit of large on-site Hilbert space dimension. For the generic case, we find strong evidence of multifractal scaling of correlation functions at the critical point, reflected in a continuous spectrum of scaling dimensions.

Quantum Ergodicity in the Many-Body Localization Problem.

We generalize Page's result on the entanglement entropy of random pure states to the many-body eigenstates of realistic disordered many-body systems subject to long-range interactions. This extension leads to two principal conclusions: first, for increasing disorder the "shells" of constant energy supporting a system's eigenstates fill only a fraction of its full Fock space and are subject to intrinsic correlations absent in synthetic high-dimensional random lattice systems. Second, in all regimes preceding the many-body localization transition individual eigenstates are thermally distributed over these shells. These results, corroborated by comparison to exact diagonalization for an SYK model, are at variance with the concept of "nonergodic extended states" in many-body systems discussed in the recent literature.

Observation of antiferromagnetic correlations in the Hubbard model with ultracold atoms.

Ultracold atoms in optical lattices have great potential to contribute to a better understanding of some of the most important issues in many-body physics, such as high-temperature superconductivity. The Hubbard model--a simplified representation of fermions moving on a periodic lattice--is thought to describe the essential details of copper oxide superconductivity. This model describes many of the features shared by the copper oxides, including an interaction-driven Mott insulating state and an antiferromagnetic (AFM) state. Optical lattices filled with a two-spin-component Fermi gas of ultracold atoms can faithfully realize the Hubbard model with readily tunable parameters, and thus provide a platform for the systematic exploration of its phase diagram. Realization of strongly correlated phases, however, has been hindered by the need to cool the atoms to temperatures as low as the magnetic exchange energy, and also by the lack of reliable thermometry. Here we demonstrate spin-sensitive Bragg scattering of light to measure AFM spin correlations in a realization of the three-dimensional Hubbard model at temperatures down to 1.4 times that of the AFM phase transition. This temperature regime is beyond the range of validity of a simple high-temperature series expansion, which brings our experiment close to the limit of the capabilities of current numerical techniques, particularly at metallic densities. We reach these low temperatures using a compensated optical lattice technique, in which the confinement of each lattice beam is compensated by a blue-detuned laser beam. The temperature of the atoms in the lattice is deduced by comparing the light scattering to determinant quantum Monte Carlo simulations and numerical linked-cluster expansion calculations. Further refinement of the compensated lattice may produce even lower temperatures which, along with light scattering thermometry, would open avenues for producing and characterizing other novel quantum states of matter, such as the pseudogap regime and correlated metallic states of the two-dimensional Hubbard model.

Scalable Probes of Measurement-Induced Criticality.

We uncover a local order parameter for measurement-induced phase transitions: the average entropy of a single reference qubit initially entangled with the system. Using this order parameter, we identify scalable probes of measurement-induced criticality that are immediately applicable to advanced quantum computing platforms. We test our proposal on a 1+1 dimensional stabilizer circuit model that can be classically simulated in polynomial time. We introduce the concept of a "decoding light cone" to establish the local and efficiently measurable nature of this probe. We also estimate bulk and surface critical exponents for the transition. Developing scalable probes of measurement-induced criticality in more general models may be a useful application of noisy intermediate scale quantum devices, as well as point to more efficient realizations of fault-tolerant quantum computation.

Localization as an Entanglement Phase Transition in Boundary-Driven Anderson Models.

The Anderson localization transition is one of the most well studied examples of a zero temperature quantum phase transition. On the other hand, many open questions remain about the phenomenology of disordered systems driven far out of equilibrium. Here we study the localization transition in the prototypical three-dimensional, noninteracting Anderson model when the system is driven at its boundaries to induce a current carrying nonequilibrium steady state. Recently we showed that the diffusive phase of this model exhibits extensive mutual information of its nonequilibrium steady-state density matrix. We show that this extensive scaling persists in the entanglement and at the localization critical point, before crossing over to a short-range (area-law) scaling in the localized phase. We introduce an entanglement witness for fermionic states that we name the mutual coherence, which, for fermionic Gaussian states, is also a lower bound on the mutual information. Through a combination of analytical arguments and numerics, we determine the finite-size scaling of the mutual coherence across the transition. These results further develop the notion of entanglement phase transitions in open systems, with direct implications for driven many-body localized systems, as well as experimental studies of driven-disordered systems.

Weyl Semimetal to Metal Phase Transitions Driven by Quasiperiodic Potentials.

We explore the stability of three-dimensional Weyl and Dirac semimetals subject to quasiperiodic potentials. We present numerical evidence that the semimetal is stable for weak quasiperiodic potentials, despite being unstable for weak random potentials. As the quasiperiodic potential strength increases, the semimetal transitions to a metal, then to an "inverted" semimetal, and then finally to a metal again. The semimetal and metal are distinguished by the density of states at the Weyl point, as well as by level statistics, transport, and the momentum-space structure of eigenstates near the Weyl point. The critical properties of the transitions in quasiperiodic systems differ from those in random systems: we do not find a clear critical scaling regime in energy; instead, at the quasiperiodic transitions, the density of states appears to jump abruptly (and discontinuously to within our resolution).

Light-Cone Spreading of Perturbations and the Butterfly Effect in a Classical Spin Chain.

We find that the effects of a localized perturbation in a chaotic classical many-body system-the classical Heisenberg chain at infinite temperature-spread ballistically with a finite speed even when the local spin dynamics is diffusive. We study two complementary aspects of this butterfly effect: the rapid growth of the perturbation, and its simultaneous ballistic (light-cone) spread, as characterized by the Lyapunov exponents and the butterfly speed, respectively. We connect this to recent studies of the out-of-time-ordered commutators (OTOC), which have been proposed as an indicator of chaos in a quantum system. We provide a straightforward identification of the OTOC with a natural correlator in our system and demonstrate that many of its interesting qualitative features are present in the classical system. Finally, by analyzing the scaling forms, we relate the growth, spread, and propagation of the perturbation with the growth of one-dimensional interfaces described by the Kardar-Parisi-Zhang equation.

Two Universality Classes for the Many-Body Localization Transition.

We provide a systematic comparison of the many-body localization (MBL) transition in spin chains with nonrandom quasiperiodic versus random fields. We find evidence suggesting that these belong to two separate universality classes: the first dominated by "intrinsic" intrasample randomness, and the second dominated by external intersample quenched randomness. We show that the effects of intersample quenched randomness are strongly growing, but not yet dominant, at the system sizes probed by exact-diagonalization studies on random models. Thus, the observed finite-size critical scaling collapses in such studies appear to be in a preasymptotic regime near the nonrandom universality class, but showing signs of the initial crossover towards the external-randomness-dominated universality class. Our results provide an explanation for why exact-diagonalization studies on random models see an apparent scaling near the transition while also obtaining finite-size scaling exponents that strongly violate Harris-Chayes bounds that apply to disorder-driven transitions. We also show that the MBL phase is more stable for the quasiperiodic model as compared to the random one, and the transition in the quasiperiodic model suffers less from certain finite-size effects.

Thermal inclusions: how one spin can destroy a many-body localized phase.

Many-body localized (MBL) systems lie outside the framework of statistical mechanics, as they fail to equilibrate under their own quantum dynamics. Even basic features of MBL systems, such as their stability to thermal inclusions and the nature of the dynamical transition to thermalizing behaviour, remain poorly understood. We study a simple central spin model to address these questions: a two-level system interacting with strength J with N≫1 localized bits subject to random fields. On increasing J, the system transitions from an MBL to a delocalized phase on the vanishing scale Jc(N)∼1/N, up to logarithmic corrections. In the transition region, the single-site eigenstate entanglement entropies exhibit bimodal distributions, so that localized bits are either 'on' (strongly entangled) or 'off' (weakly entangled) in eigenstates. The clusters of 'on' bits vary significantly between eigenstates of the same sample, which provides evidence for a heterogeneous discontinuous transition out of the localized phase in single-site observables. We obtain these results by perturbative mapping to bond percolation on the hypercube at small J and by numerical exact diagonalization of the full many-body system. Our results support the arguments that the MBL phase is unstable in systems with short-range interactions and quenched randomness in dimensions d that are high but finite.This article is part of the themed issue 'Breakdown of ergodicity in quantum systems: from solids to synthetic matter'.

Obtaining highly excited eigenstates of the localized XX chain via DMRG-X.

We benchmark a variant of the recently introduced density matrix renormalization group (DMRG)-X algorithm against exact results for the localized random field XX chain. We find that the eigenstates obtained via DMRG-X exhibit a highly accurate l-bit description for system sizes much bigger than the direct, many-body, exact diagonalization in the spin variables is able to access. We take advantage of the underlying free fermion description of the XX model to accurately test the strengths and limitations of this algorithm for large system sizes. We discuss the theoretical constraints on the performance of the algorithm from the entanglement properties of the eigenstates, and its actual performance at different values of disorder. A small but significant improvement to the algorithm is also presented, which helps significantly with convergence. We find that, at high entanglement, DMRG-X shows a bias towards eigenstates with low entanglement, but can be improved with increased bond dimension. This result suggests that one must be careful when applying the algorithm for interacting many-body localized spin models near a transition.This article is part of the themed issue 'Breakdown of ergodicity in quantum systems: from solids to synthetic matter'.

Thermal Equilibrium of a Macroscopic Quantum System in a Pure State.

We consider the notion of thermal equilibrium for an individual closed macroscopic quantum system in a pure state, i.e., described by a wave function. The macroscopic properties in thermal equilibrium of such a system, determined by its wave function, must be the same as those obtained from thermodynamics, e.g., spatial uniformity of temperature and chemical potential. When this is true we say that the system is in macroscopic thermal equilibrium (MATE). Such a system may, however, not be in microscopic thermal equilibrium (MITE). The latter requires that the reduced density matrices of small subsystems be close to those obtained from the microcanonical, equivalently the canonical, ensemble for the whole system. The distinction between MITE and MATE is particularly relevant for systems with many-body localization for which the energy eigenfuctions fail to be in MITE while necessarily most of them, but not all, are in MATE. We note, however, that for generic macroscopic systems, including those with MBL, most wave functions in an energy shell are in both MATE and MITE. For a classical macroscopic system, MATE holds for most phase points on the energy surface, but MITE fails to hold for any phase point.

Speed, dissipation, and error in kinetic proofreading.

Proofreading mechanisms increase specificity in biochemical reactions by allowing for the dissociation of intermediate complexes. These mechanisms disrupt and reset the reaction to undo errors at the cost of increased time of reaction and free energy expenditure. Here, we draw an analogy between proofreading and microtubule growth which share some of the features described above. Our analogy relates the statistics of growth and shrinkage of microtubules in physical space to the cycling of intermediate complexes in the space of chemical states in proofreading mechanisms. Using this analogy, we find a new kinetic regime of proofreading in which an exponential speed-up of the process can be achieved at the cost of a somewhat larger error rate. This regime is analogous to the transition region between two known growth regimes of microtubules (bounded and unbounded) and is sharply defined in the limit of large proofreading networks. We find that this advantageous regime of speed-error tradeoff might be present in proofreading schemes studied earlier in the charging of tRNA by tRNA synthetases, in RecA filament assembly on ssDNA, and in protein synthesis by ribosomes.

Ballistic spreading of entanglement in a diffusive nonintegrable system.

We study the time evolution of the entanglement entropy of a one-dimensional nonintegrable spin chain, starting from random nonentangled initial pure states. We use exact diagonalization of a nonintegrable quantum Ising chain with transverse and longitudinal fields to obtain the exact quantum dynamics. We show that the entanglement entropy increases linearly with time before finite-size saturation begins, demonstrating a ballistic spreading of the entanglement, while the energy transport in the same system is diffusive. Thus, we explicitly demonstrate that the spreading of entanglement is much faster than the energy diffusion in this nonintegrable system.

Weak plaquette valence bond order in the S = 1/2 honeycomb J1 - J2 Heisenberg model.

Using the density matrix renormalization group, we investigate the S = 1/2 Heisenberg model on the honeycomb lattice with first (J(1)) and second (J(2)) neighbor interactions. We are able to study long open cylinders with widths up to 12 lattice spacings. For J(2)/J(1) near 0.3, we find an apparently paramagnetic phase, bordered by an antiferromagnetic phase for J(2) ≲ 0.26 and by a valence bond crystal for J(2) ≳ 0.36. The longest correlation length that we find in this intermediate phase is for plaquette valence bond order. This correlation length grows strongly with cylinder circumference, indicating either quantum criticality or weak plaquette valence bond order.