Quantum thermalization through entanglement in an isolated many-body system

Probing spin hydrodynamics on a superconducting quantum simulator

Characterizing the nature of hydrodynamical transport properties in quantum dynamics provides valuable insights into the fundamental understanding of exotic non-equilibrium phases of matter. Experimentally simulating infinite-temperature transport on large-scale complex quantum systems is of considerable interest. Here, using a controllable and coherent superconducting quantum simulator, we experimentally realize the analog quantum circuit, which can efficiently prepare the Haar-random states, and probe spin transport at infinite temperature. We observe diffusive spin transport during the unitary evolution of the ladder-type quantum simulator with ergodic dynamics. Moreover, we explore the transport properties of the systems subjected to strong disorder or a tilted potential, revealing signatures of anomalous subdiffusion in accompany with the breakdown of thermalization. Our work demonstrates a scalable method of probing infinite-temperature spin transport on analog quantum simulators, which paves the way to study other intriguing out-of-equilibrium phenomena from the perspective of transport.

Localized Virtual Purification

Analog and digital quantum simulators can efficiently simulate quantum many-body systems that appear in natural phenomena. However, experimental limitations of near-term devices still make it challenging to perform the entire process of quantum simulation. The purification-based quantum simulation methods can alleviate the limitations in experiments such as the cooling temperature and noise from the environment, while this method has the drawback that it requires global entangled measurement with a prohibitively large number of measurements that scales exponentially with the system size. In this Letter, we propose that we can overcome these problems by restricting the entangled measurements to the vicinity of the local observables to be measured, when the locality of the system can be exploited. We provide theoretical guarantees that the global purification operation can be replaced with local operations under some conditions, in particular for the task of cooling and error mitigation. We furthermore give a numerical verification that the localized purification is valid even when conditions are not satisfied. Our method bridges the fundamental concept of locality with quantum simulators, and therefore is expected to open a path to unexplored quantum many-body phenomena.

Local Readout and Control of Current and Kinetic Energy Operators in Optical Lattices

Quantum gas microscopes have revolutionized quantum simulations with ultracold atoms, allowing one to measure local observables and snapshots of quantum states. However, measurements so far were mostly carried out in the occupation basis. Here, we demonstrate how all kinetic operators, such as kinetic energy or current operators, can be measured and manipulated with single-bond resolution. Beyond simple expectation values of these observables, the single-shot measurements allow one to access full counting statistics and complex correlation functions. Our work paves the way for the implementation of efficient quantum state tomography and hybrid quantum computing protocols for itinerant particles on a lattice. In addition, we demonstrate how site-resolved programmable potentials enable a spatially selective, parallel readout in different bases as well as the engineering of arbitrary initial states.

Random matrix theory approach to quantum Fisher information in quantum ergodic systems

We theoretically investigate quantum parameter estimation in quantum chaotic systems. Our analysis is based on an effective description of quantum ergodic systems in terms of a random matrix Hamiltonian. Based on this approach, we derive an analytical expression for the time evolution of the quantum Fisher information (QFI), which we find to have three distinct timescales. Initially, the QFI increase is quadratic in time, characterizing the timescale over which initial information is extractable from local measurements only. This quickly passes into linear increase with slope determined by the decay rate of the measured spin observable. When the information is fully spread among all degrees of freedom, a second quadratic timescale determines the long-time behavior of the QFI. We test our random matrix theory prediction with the exact diagonalization of a nonintegrable spin system, focusing on the estimation of a local magnetic field by measurements of the many-body state. Our numerical calculations agree with the effective random matrix theory approach and show that the information on the local Hamiltonian parameter is distributed throughout the quantum system during the quantum thermalization process.

Identifying Sources of Entanglement Loss in Photodriven Molecular Electron Spin Teleportation

We report on an electron donor-electron acceptor-stable radical (D-A-R•) molecule in which an electron spin state first prepared on R• is followed by photogeneration of an entangled singlet 1[D•+-A•-] spin pair to produce D•+-A•--R•. Since the A•- and R• spins within D•+-A•--R• are uncorrelated, spin teleportation from R• to D•+ occurs with a maximal 25% efficiency only for the singlet pair 1(A•--R•) by spin-allowed electron transfer from A•- to R•. However, since 1[D•+-A•-] is sufficiently long-lived, coherent spin mixing involving the unreactive 3(A•--R•) population affects entanglement and teleportation within D•+-A•--R•. Pulse electron paramagnetic resonance experiments show a direct correlation between electron spin flip-flops and entanglement loss, providing information for designing molecular materials to serve as nanoscale quantum device interconnects.

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.

Observing the Quantum Mpemba Effect in Quantum Simulations

The nonequilibrium physics of many-body quantum systems harbors various unconventional phenomena. In this Letter, we experimentally investigate one of the most puzzling of these phenomena-the quantum Mpemba effect, where a tilted ferromagnet restores its symmetry more rapidly when it is farther from the symmetric state compared to when it is closer. We present the first experimental evidence of the occurrence of this effect in a trapped-ion quantum simulator. The symmetry breaking and restoration are monitored through entanglement asymmetry, probed via randomized measurements, and postprocessed using the classical shadows technique. Our findings are further substantiated by measuring the Frobenius distance between the experimental state and the stationary thermal symmetric theoretical state, offering direct evidence of subsystem thermalization.

Structure of the Hamiltonian of mean force

The Hamiltonian of mean force is an effective Hamiltonian that allows a quantum system, nonweakly coupled to an environment, to be written in an effective Gibbs state. We present results on the structure of the Hamiltonian of mean force in extended quantum systems with local interactions. We show that its spatial structure exhibits a "skin effect"-its difference from the system Hamiltonian dies off exponentially with distance from the system-environment boundary. For spin systems, we identify the terms that can appear in the Hamiltonian of mean force at different orders in the inverse temperature.

Bath-induced interactions and transient dynamics in open quantum systems at strong coupling: Effective Hamiltonian approach

Understanding the dynamics of dissipative quantum systems, particularly beyond the weak coupling approximation, is central to various quantum applications. While numerically exact methods provide accurate solutions, they often lack the analytical insight provided by theoretical approaches. In this study, we employ the recently developed method dubbed the effective Hamiltonian theory to understand the dynamics of system-bath configurations without resorting to a perturbative description of the system-bath coupling energy. Through a combination of mapping steps and truncation, the effective Hamiltonian theory offers both analytical insights into signatures of strong couplings in open quantum systems and a straightforward path for numerical simulations. To validate the accuracy of the method, we apply it to two canonical models: a single spin immersed in a bosonic bath and two noninteracting spins in a common bath. In both cases, we study the transient regime and the steady state limit at nonzero temperature and spanning system-bath interactions from the weak to the strong regime. By comparing the results of the effective Hamiltonian theory with numerically exact simulations, we show that although the former overlooks non-Markovian features in the transient equilibration dynamics, it correctly captures non-perturbative bath-generated couplings between otherwise non-interacting spins, as observed in their synchronization dynamics and correlations. Altogether, the effective Hamiltonian theory offers a powerful approach for understanding strong coupling dynamics and thermodynamics, capturing the signatures of such interactions in both relaxation dynamics and in the steady state limit.

Chaotic route to classical thermalization: A real-space analysis

Most of the previous studies on classical thermalization focus on the wave-vector space, encountering limitations when extended beyond quasi-integrable regions. In this study, we propose a scheme to study the thermalization of the classical Hamiltonian chain of interacting oscillators in real space by developing a thermalization indicator proposed by Parisi [Europhys. Lett. 40, 357 (1997)0295-507510.1209/epl/i1997-00471-9], which approaches zero in the thermal state. Upon reaching the steady state characterized by the generalized Gibbs ensemble for a harmonic chain, a quench protocol is implemented to change the Hamiltonian to a nonintegrable form instantaneously, thereby preparing nonequilibrium initial states. This approach enables investigations of thermalization in real space, particularly valuable for exploring regions beyond quasi-integrability. For the FPUT-β lattice, we observe that the thermalization time as a function of the nonintegrable strength follows a -2 scaling law in the quasi-integrable region and -1/4 in the strongly integrable region. Moreover, numerical results reveal the thermalization time is proportional to the Lyapunov time, which bridges microscopic chaotic dynamics and the macroscopic thermalization process.

Selective Detection of Dynamics-Complete Set of Correlations via Quantum Channels

Correlations of fluctuations are essential to understanding many-body systems and key information for advancing quantum technologies. To fully describe the dynamics of a physical system, all time-ordered correlations (TOCs), i.e., the dynamics-complete set of correlations are needed. The current measurement techniques can only access a limited set of TOCs, and there has been no systematic and feasible solution for extracting the dynamic-complete set of correlations hitherto. Here we propose a platform-universal protocol to selectively detect arbitrary types of TOCs via quantum channels. In our method, the quantum channels are synthesized with various controls, and engineer the evolution of a sensor-target system along a specific path that corresponds to a desired correlation. Using nuclear magnetic resonance, we experimentally demonstrate this protocol by detecting a specific type of fourth-order TOC that has never been accessed previously. We also show that the knowledge of the TOCs can be used to significantly improve the precision of quantum optimal control. Our method provides a new toolbox for characterizing the quantum many-body states and quantum noise, and hence for advancing the fields of quantum sensing and quantum computing.

Subsystem eigenstate thermalization hypothesis for translation invariant systems

The eigenstate thermalization hypothesis for translation invariant quantum spin systems has been proved recently by using random matrices. In this paper, we study the subsystem version of the eigenstate thermalization hypothesis for translation invariant quantum systems without referring to random matrices. We first find a relation between the quantum variance and the Belavkin-Staszewski relative entropy. Then, by showing the small upper bounds on the quantum variance and the Belavkin-Staszewski relative entropy, we prove the subsystem eigenstate thermalization hypothesis for translation invariant quantum systems with an algebraic speed of convergence in an elementary way. The proof holds for most of the translation invariant quantum lattice models with exponential or algebraic decays of correlations.

An atomic boson sampler

A boson sampler implements a restricted model of quantum computing. It is defined by the ability to sample from the distribution resulting from the interference of identical bosons propagating according to programmable, non-interacting dynamics1. An efficient exact classical simulation of boson sampling is not believed to exist, which has motivated ground-breaking boson sampling experiments in photonics with increasingly many photons2-12. However, it is difficult to generate and reliably evolve specific numbers of photons with low loss, and thus probabilistic techniques for postselection7 or marked changes to standard boson sampling10-12 are generally used. Here, we address the above challenges by implementing boson sampling using ultracold atoms13,14 in a two-dimensional, tunnel-coupled optical lattice. This demonstration is enabled by a previously unrealized combination of tools involving high-fidelity optical cooling and imaging of atoms in a lattice, as well as programmable control of those atoms using optical tweezers. When extended to interacting systems, our work demonstrates the core abilities required to directly assemble ground and excited states in simulations of various Hubbard models15,16.

Probing entanglement in a 2D hard-core Bose-Hubbard lattice

Entanglement and its propagation are central to understanding many physical properties of quantum systems1-3. Notably, within closed quantum many-body systems, entanglement is believed to yield emergent thermodynamic behaviour4-7. However, a universal understanding remains challenging owing to the non-integrability and computational intractability of most large-scale quantum systems. Quantum hardware platforms provide a means to study the formation and scaling of entanglement in interacting many-body systems8-14. Here we use a controllable 4 × 4 array of superconducting qubits to emulate a 2D hard-core Bose-Hubbard (HCBH) lattice. We generate superposition states by simultaneously driving all lattice sites and extract correlation lengths and entanglement entropy across its many-body energy spectrum. We observe volume-law entanglement scaling for states at the centre of the spectrum and a crossover to the onset of area-law scaling near its edges.

Irreversible entropy transport enhanced by fermionic superfluidity

The nature of particle and entropy flow between two superfluids is often understood in terms of reversible flow carried by an entropy-free, macroscopic wavefunction. While this wavefunction is responsible for many intriguing properties of superfluids and superconductors, its interplay with excitations in non-equilibrium situations is less understood. Here we observe large concurrent flows of both particles and entropy through a ballistic channel connecting two strongly interacting fermionic superfluids. Both currents respond nonlinearly to chemical potential and temperature biases. We find that the entropy transported per particle is much larger than the prediction of superfluid hydrodynamics in the linear regime and largely independent of changes in the channel's geometry. By contrast, the timescales of advective and diffusive entropy transport vary significantly with the channel geometry. In our setting, superfluidity counterintuitively increases the speed of entropy transport. Moreover, we develop a phenomenological model describing the nonlinear dynamics within the framework of generalized gradient dynamics. Our approach for measuring entropy currents may help elucidate mechanisms of heat transfer in superfluids and superconducting devices.

Benchmarking highly entangled states on a 60-atom analogue quantum simulator

Quantum systems have entered a competitive regime in which classical computers must make approximations to represent highly entangled quantum states1,2. However, in this beyond-classically-exact regime, fidelity comparisons between quantum and classical systems have so far been limited to digital quantum devices2-5, and it remains unsolved how to estimate the actual entanglement content of experiments6. Here, we perform fidelity benchmarking and mixed-state entanglement estimation with a 60-atom analogue Rydberg quantum simulator, reaching a high-entanglement entropy regime in which exact classical simulation becomes impractical. Our benchmarking protocol involves extrapolation from comparisons against an approximate classical algorithm, introduced here, with varying entanglement limits. We then develop and demonstrate an estimator of the experimental mixed-state entanglement6, finding our experiment is competitive with state-of-the-art digital quantum devices performing random circuit evolution2-5. Finally, we compare the experimental fidelity against that achieved by various approximate classical algorithms, and find that only the algorithm we introduce is able to keep pace with the experiment on the classical hardware we use. Our results enable a new model for evaluating the ability of both analogue and digital quantum devices to generate entanglement in the beyond-classically-exact regime, and highlight the evolving divide between quantum and classical systems.

Estimating entanglement entropy via variational quantum circuits with classical neural networks

Entropy plays a crucial role in both physics and information science, encompassing classical and quantum domains. In this paper, we present the quantum neural entropy estimator (QNEE), an approach that combines classical neural network (NN) with variational quantum circuits to estimate the von Neumann and Rényi entropies of a quantum state. QNEE provides accurate estimates of entropy while also yielding the eigenvalues and eigenstates of the input density matrix. Leveraging the capabilities of classical NN, QNEE can classify different phases of quantum systems that accompany the changes of entanglement entropy. Our numerical simulation demonstrates the effectiveness of QNEE by applying it to the 1D XXZ Heisenberg model. In particular, QNEE exhibits high sensitivity in estimating entanglement entropy near the phase transition point. We expect that QNEE will serve as a valuable tool for quantum entropy estimation and phase classification.

Stable quantum-correlated many-body states through engineered dissipation

Engineered dissipative reservoirs have the potential to steer many-body quantum systems toward correlated steady states useful for quantum simulation of high-temperature superconductivity or quantum magnetism. Using up to 49 superconducting qubits, we prepared low-energy states of the transverse-field Ising model through coupling to dissipative auxiliary qubits. In one dimension, we observed long-range quantum correlations and a ground-state fidelity of 0.86 for 18 qubits at the critical point. In two dimensions, we found mutual information that extends beyond nearest neighbors. Lastly, by coupling the system to auxiliaries emulating reservoirs with different chemical potentials, we explored transport in the quantum Heisenberg model. Our results establish engineered dissipation as a scalable alternative to unitary evolution for preparing entangled many-body states on noisy quantum processors.

Many-body quantum chaos in mixtures of multiple species

We study spectral correlations in many-body quantum mixtures of fermions, bosons, and qubits with periodically kicked spreading and mixing of species. We take two types of mixing, namely, Jaynes-Cummings and Rabi, respectively, satisfying and breaking the conservation of a total number of species. We analytically derive the generating Hamiltonians whose spectral properties determine the spectral form factor in the leading order. We further analyze the system-size (L) scaling of Thouless time t^{*}, beyond which the spectral form factor follows the prediction of random matrix theory. The L dependence of t^{*} crosses over from lnL to L^{2} with an increasing Jaynes-Cummings mixing between qubits and fermions or bosons in a finite-size chain, and it finally settles to t^{*}∝O(L^{2}) in the thermodynamic limit for any mixing strength. The Rabi mixing between qubits and fermions leads to t^{*}∝O(lnL), previously predicted for single species of qubits or fermions without total-number conservation.

Mesoscopic fluctuations in entanglement dynamics

Understanding fluctuation phenomena plays a dominant role in the development of many-body physics. The time evolution of entanglement is essential to a broad range of subjects in many-body physics, ranging from exotic quantum matter to quantum thermalization. Stemming from various dynamical processes of information, fluctuations in entanglement evolution differ conceptually from out-of-equilibrium fluctuations of traditional physical quantities. Their studies remain elusive. Here we uncover an emergent random structure in the evolution of the many-body wavefunction in two classes of integrable-either interacting or noninteracting-lattice models. It gives rise to out-of-equilibrium entanglement fluctuations which fall into the paradigm of mesoscopic fluctuations of wave interference origin. Specifically, the entanglement entropy variance obeys a universal scaling law in each class, and the full distribution displays a sub-Gaussian upper and a sub-Gamma lower tail. These statistics are independent of both the system's microscopic details and the choice of entanglement probes, and broaden the class of mesoscopic universalities. They have practical implications for controlling entanglement in mesoscopic devices.

Logical quantum processor based on reconfigurable atom arrays

Suppressing errors is the central challenge for useful quantum computing1, requiring quantum error correction (QEC)2-6 for large-scale processing. However, the overhead in the realization of error-corrected 'logical' qubits, in which information is encoded across many physical qubits for redundancy2-4, poses substantial challenges to large-scale logical quantum computing. Here we report the realization of a programmable quantum processor based on encoded logical qubits operating with up to 280 physical qubits. Using logical-level control and a zoned architecture in reconfigurable neutral-atom arrays7, our system combines high two-qubit gate fidelities8, arbitrary connectivity7,9, as well as fully programmable single-qubit rotations and mid-circuit readout10-15. Operating this logical processor with various types of encoding, we demonstrate improvement of a two-qubit logic gate by scaling surface-code6 distance from d = 3 to d = 7, preparation of colour-code qubits with break-even fidelities5, fault-tolerant creation of logical Greenberger-Horne-Zeilinger (GHZ) states and feedforward entanglement teleportation, as well as operation of 40 colour-code qubits. Finally, using 3D [[8,3,2]] code blocks16,17, we realize computationally complex sampling circuits18 with up to 48 logical qubits entangled with hypercube connectivity19 with 228 logical two-qubit gates and 48 logical CCZ gates20. We find that this logical encoding substantially improves algorithmic performance with error detection, outperforming physical-qubit fidelities at both cross-entropy benchmarking and quantum simulations of fast scrambling21,22. These results herald the advent of early error-corrected quantum computation and chart a path towards large-scale logical processors.

Quantifying spatiotemporal patterns in classical and quantum systems out of equilibrium

A rich variety of nonequilibrium dynamical phenomena and processes unambiguously calls for the development of general numerical techniques to probe and estimate a complex interplay between spatial and temporal degrees of freedom in many-body systems of completely different nature. In this work we provide a solution to this problem by adopting a structural complexity measure to quantify spatiotemporal patterns in the time-dependent digital representation of a system. On the basis of very limited amount of data our approach allows us to distinguish different dynamical regimes and define critical parameters in both classical and quantum systems. By the example of the discrete time crystal realized in nonequilibrium quantum systems we provide a complete low-level characterization of this nontrivial dynamical phase with only processing bitstrings, which can be considered as a valuable alternative to previous studies based on the calculations of qubit correlation functions.

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.

Finite-Time Teleportation Phase Transition in Random Quantum Circuits

How long does it take to entangle two distant qubits in a quantum circuit evolved by generic unitary dynamics? We show that if the time evolution is followed by measurements of all but two infinitely separated test qubits, then the entanglement between them can undergo a phase transition and become nonzero at a finite critical time t_{c}. The fidelity of teleporting a quantum state from an input qubit to an infinitely distant output qubit shows the same critical onset. Specifically, these finite-time transitions occur in short-range interacting two-dimensional random unitary circuits and in sufficiently long-range interacting one-dimensional circuits. The phase transition is understood by mapping the random continuous-time evolution to a finite-temperature thermal state of an effective spin Hamiltonian, where the inverse temperature equals the evolution time in the circuit. In this framework, the entanglement between two distant qubits at times t>t_{c} corresponds to the emergence of long-range ferromagnetic spin correlations below the critical temperature. We verify these predictions using numerical simulation of Clifford circuits and propose potential realizations in existing platforms for quantum simulation.

Commensurate and incommensurate 1D interacting quantum systems

Single-atom imaging resolution of many-body quantum systems in optical lattices is routinely achieved with quantum-gas microscopes. Key to their great versatility as quantum simulators is the ability to use engineered light potentials at the microscopic level. Here, we employ dynamically varying microscopic light potentials in a quantum-gas microscope to study commensurate and incommensurate 1D systems of interacting bosonic Rb atoms. Such incommensurate systems are analogous to doped insulating states that exhibit atom transport and compressibility. Initially, a commensurate system with unit filling and fixed atom number is prepared between two potential barriers. We deterministically create an incommensurate system by dynamically changing the position of the barriers such that the number of available lattice sites is reduced while retaining the atom number. Our systems are characterised by measuring the distribution of particles and holes as a function of the lattice filling, and interaction strength, and we probe the particle mobility by applying a bias potential. Our work provides the foundation for preparation of low-entropy states with controlled filling in optical-lattice experiments.

Particle-hole thermalization in a composite superconducting and normal-conducting nanowire

The mechanisms by which isolated condensed matter systems thermalize is a topic of growing interest. Thermalization is known to be linked to the emergence of chaos in the dynamics of a system. We show that a solid state scattering system, containing superconducting elements, can thermalize scattered states without affecting the degree of entanglement of the scattered states. We consider a composite NSNSNSNSN nanowire composed of Bi_{2}Sr_{2}CaCu_{2}O_{8+x} superconducting segments (S) and normal conducting segments (N). We consider parameter regimes where all current flow is due to tunneling currents that are facilitated by quasibound state resonances inside the SNSNSNS structure. At certain energies, scattered pure states approach ergodicity, even though they remain pure.

Unifying Emergent Hydrodynamics and Lindbladian Low-Energy Spectra across Symmetries, Constraints, and Long-Range Interactions

We identify emergent hydrodynamics governing charge transport in Brownian random time evolution with various symmetries, constraints, and ranges of interactions. This is accomplished via a mapping between the averaged dynamics and the low-energy spectrum of a Lindblad operator, which acts as an effective Hamiltonian in a doubled Hilbert space. By explicitly constructing dispersive excited states of this effective Hamiltonian using a single-mode approximation, we provide a comprehensive understanding of diffusive, subdiffusive, and superdiffusive relaxation in many-body systems with conserved multipole moments and variable interaction ranges. Our approach further allows us to identify exotic Krylov-space-resolved diffusive relaxation despite the presence of dipole conservation, which we verify numerically. Therefore, we provide a general and versatile framework to qualitatively understand the dynamics of conserved operators under random unitary time evolution.

Sunburst quantum Ising model under interaction quench: Entanglement and role of initial state coherence

We study the nonequilibrium dynamics of an isolated bipartite quantum system, the sunburst quantum Ising model, under interaction quench. The prequench limit of this model is two noninteracting integrable systems, namely a transverse Ising chain and finite number of isolated qubits. As a function of interaction strength, the spectral fluctuation property goes from Poisson to Wigner-Dyson statistics. We chose entanglement entropy as a probe to study the approach to thermalization or lack of it in postquench dynamics. In the near-integrable limit, as expected, the linear entropy displays oscillatory behavior, while in the chaotic limit it saturates. Along with the chaotic nature of the time evolution generator, we show the importance of the role played by the coherence of the initial state in deciding the nature of thermalization. We further show that these findings are general by replacing the Ising ring with a disordered XXZ model with disorder strength putting it in the many-body localized phase.

Entangling Dynamics from Effective Rotor-Spin-Wave Separation in U(1)-Symmetric Quantum Spin Models

The nonequilibrium dynamics of quantum spin models is a most challenging topic, due to the exponentiality of Hilbert space, and it is central to the understanding of the many-body entangled states that can be generated by state-of-the-art quantum simulators. A particularly important class of evolutions is the one governed by U(1)-symmetric Hamiltonians, initialized in a state that breaks the U(1) symmetry-the paradigmatic example being the evolution of the so-called one-axis-twisting (OAT) model, featuring infinite-range interactions between spins. In this Letter, we show that the dynamics of the OAT model can be closely reproduced by systems with power-law-decaying interactions, thanks to an effective separation between the zero-momentum degrees of freedom, associated with the so-called Anderson tower of states, and reconstructing an OAT model, as well as finite-momentum ones, associated with spin-wave excitations. This mechanism explains quantitatively the recent numerical observation of spin squeezing and Schrödinger cat-state generation in the dynamics of dipolar Hamiltonians, and it paves the way for the extension of this observation to a much larger class of models of immediate relevance for quantum simulations.

General Class of Continuous Variable Entanglement Criteria

We present a general class of entanglement criteria for continuous variable systems. Our criteria are based on the Husimi Q distribution and allow for optimization over the set of all concave functions rendering them extremely general and versatile. We show that several entropic criteria and second moment criteria are obtained as special cases. Our criteria reveal entanglement of families of states undetected by any commonly used criteria and provide clear advantages under typical experimental constraints such as finite detector resolution and measurement statistics.

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.

Benchmarking Quantum Simulators Using Ergodic Quantum Dynamics

We propose and analyze a sample-efficient protocol to estimate the fidelity between an experimentally prepared state and an ideal target state, applicable to a wide class of analog quantum simulators without advanced spatiotemporal control. Our protocol relies on universal fluctuations emerging from generic Hamiltonian dynamics, which we discover in the present work. It does not require fine-tuned control over state preparation, quantum evolution, or readout capability, while achieving near optimal sample complexity: a percent-level precision is obtained with ∼10^{3} measurements, independent of system size. Furthermore, the accuracy of our fidelity estimation improves exponentially with increasing system size. We numerically demonstrate our protocol in a variety of quantum simulator platforms, including quantum gas microscopes, trapped ions, and Rydberg atom arrays. We discuss applications of our method for tasks such as multiparameter estimation of quantum states and processes.

Realizing spin squeezing with Rydberg interactions in an optical clock

Neutral-atom arrays trapped in optical potentials are a powerful platform for studying quantum physics, combining precise single-particle control and detection with a range of tunable entangling interactions1-3. For example, these capabilities have been leveraged for state-of-the-art frequency metrology4,5 as well as microscopic studies of entangled many-particle states6-11. Here we combine these applications to realize spin squeezing-a widely studied operation for producing metrologically useful entanglement-in an optical atomic clock based on a programmable array of interacting optical qubits. In this demonstration of Rydberg-mediated squeezing with a neutral-atom optical clock, we generate states that have almost four decibels of metrological gain. In addition, we perform a synchronous frequency comparison between independent squeezed states and observe a fractional-frequency stability of 1.087(1) × 10-15 at one-second averaging time, which is 1.94(1) decibels below the standard quantum limit and reaches a fractional precision at the 10-17 level during a half-hour measurement. We further leverage the programmable control afforded by optical tweezer arrays to apply local phase shifts to explore spin squeezing in measurements that operate beyond the relative coherence time with the optical local oscillator. The realization of this spin-squeezing protocol in a programmable atom-array clock will enable a wide range of quantum-information-inspired techniques for optimal phase estimation and Heisenberg-limited optical atomic clocks12-16.

Real-time broadening of bath-induced density profiles from closed-system correlation functions

The Lindblad master equation is one of the main approaches to open quantum systems. While it has been widely applied in the context of condensed matter systems to study properties of steady states in the limit of long times, the actual route to such steady states has attracted less attention yet. Here, we investigate the nonequilibrium dynamics of spin chains with a local coupling to a single Lindblad bath and analyze the transport properties of the induced magnetization. Combining typicality and equilibration arguments with stochastic unraveling, we unveil for the case of weak driving that the dynamics in the open system can be constructed on the basis of correlation functions in the closed system, which establishes a connection between the Lindblad approach and linear response theory at finite times. In this way, we provide a particular example where closed and open approaches to quantum transport agree strictly. We demonstrate this fact numerically for the spin-1/2 XXZ chain at the isotropic point and in the easy-axis regime, where superdiffusive and diffusive scaling is observed, respectively.

Deep learning of quantum entanglement from incomplete measurements

The quantification of the entanglement present in a physical system is of paramount importance for fundamental research and many cutting-edge applications. Now, achieving this goal requires either a priori knowledge on the system or very demanding experimental procedures such as full state tomography or collective measurements. Here, we demonstrate that, by using neural networks, we can quantify the degree of entanglement without the need to know the full description of the quantum state. Our method allows for direct quantification of the quantum correlations using an incomplete set of local measurements. Despite using undersampled measurements, we achieve a quantification error of up to an order of magnitude lower than the state-of-the-art quantum tomography. Furthermore, we achieve this result using networks trained using exclusively simulated data. Last, we derive a method based on a convolutional network input that can accept data from various measurement scenarios and perform, to some extent, independently of the measurement device.

Quantum simulation of thermodynamics in an integrated quantum photonic processor

One of the core questions of quantum physics is how to reconcile the unitary evolution of quantum states, which is information-preserving and time-reversible, with evolution following the second law of thermodynamics, which, in general, is neither. The resolution to this paradox is to recognize that global unitary evolution of a multi-partite quantum state causes the state of local subsystems to evolve towards maximum-entropy states. In this work, we experimentally demonstrate this effect in linear quantum optics by simultaneously showing the convergence of local quantum states to a generalized Gibbs ensemble constituting a maximum-entropy state under precisely controlled conditions, while introducing an efficient certification method to demonstrate that the state retains global purity. Our quantum states are manipulated by a programmable integrated quantum photonic processor, which simulates arbitrary non-interacting Hamiltonians, demonstrating the universality of this phenomenon. Our results show the potential of photonic devices for quantum simulations involving non-Gaussian states.

Witnessing light-driven entanglement using time-resolved resonant inelastic X-ray scattering

Characterizing and controlling entanglement in quantum materials is crucial for the development of next-generation quantum technologies. However, defining a quantifiable figure of merit for entanglement in macroscopic solids is theoretically and experimentally challenging. At equilibrium the presence of entanglement can be diagnosed by extracting entanglement witnesses from spectroscopic observables and a nonequilibrium extension of this method could lead to the discovery of novel dynamical phenomena. Here, we propose a systematic approach to quantify the time-dependent quantum Fisher information and entanglement depth of transient states of quantum materials with time-resolved resonant inelastic x-ray scattering. Using a quarter-filled extended Hubbard model as an example, we benchmark the efficiency of this approach and predict a light-enhanced many-body entanglement due to the proximity to a phase boundary. Our work sets the stage for experimentally witnessing and controlling entanglement in light-driven quantum materials via ultrafast spectroscopic measurements.

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.

A Differential-Geometric Approach to Quantum Ignorance Consistent with Entropic Properties of Statistical Mechanics

In this paper, we construct the metric tensor and volume for the manifold of purifications associated with an arbitrary reduced density operator ρS. We also define a quantum coarse-graining (CG) to study the volume where macrostates are the manifolds of purifications, which we call surfaces of ignorance (SOI), and microstates are the purifications of ρS. In this context, the volume functions as a multiplicity of the macrostates that quantifies the amount of information missing from ρS. Using examples where the SOI are generated using representations of SU(2), SO(3), and SO(N), we show two features of the CG: (1) A system beginning in an atypical macrostate of smaller volume evolves to macrostates of greater volume until it reaches the equilibrium macrostate in a process in which the system and environment become strictly more entangled, and (2) the equilibrium macrostate takes up the vast majority of the coarse-grained space especially as the dimension of the total system becomes large. Here, the equilibrium macrostate corresponds to a maximum entanglement between the system and the environment. To demonstrate feature (1) for the examples considered, we show that the volume behaves like the von Neumann entropy in that it is zero for pure states, maximal for maximally mixed states, and is a concave function with respect to the purity of ρS. These two features are essential to typicality arguments regarding thermalization and Boltzmann's original CG.

Quench-Probe Setup as Analyzer of Fractionalized Entanglement Spreading

We propose a novel spatially inhomogeneous setup for revealing quench-induced fractionalized excitations in entanglement dynamics. In this quench-probe setting, the region undergoing a quantum quench is tunnel coupled to a static region, the probe. Subsequently, the time-dependent entanglement signatures of a tunable subset of excitations propagating to the probe are monitored by energy selectivity. We exemplify the power of this generic approach by identifying a unique dynamical signature associated with the presence of an isolated Majorana zero mode in the postquench Hamiltonian. In this case excitations emitted from the topological part of the system give rise to a fractionalized jump of log(2)/2 in the entanglement entropy of the probe. This dynamical effect is highly sensitive to the localized nature of the Majorana zero mode, but does not require the preparation of a topological initial state.

Observation of entanglement transition of pseudo-random mixed states

Random quantum states serve as a powerful tool in various scientific fields, including quantum supremacy and black hole physics. It has been theoretically predicted that entanglement transitions may happen for different partitions of multipartite random quantum states; however, the experimental observation of these transitions is still absent. Here, we experimentally demonstrate the entanglement transitions witnessed by negativity on a fully connected superconducting processor. We apply parallel entangling operations, that significantly decrease the depth of the pseudo-random circuits, to generate pseudo-random pure states of up to 15 qubits. By quantum state tomography of the reduced density matrix of six qubits, we measure the negativity spectra. Then, by changing the sizes of the environment and subsystems, we observe the entanglement transitions that are directly identified by logarithmic entanglement negativities based on the negativity spectra. In addition, we characterize the randomness of our circuits by measuring the distance between the distribution of output bit-string probabilities and the Porter-Thomas distribution. Our results show that superconducting processors with all-to-all connectivity constitute a promising platform for generating random states and understanding the entanglement structure of multipartite quantum systems.

Entropy and temperature in finite isolated quantum systems

We investigate how the temperature calculated from the microcanonical entropy compares with the canonical temperature for finite isolated quantum systems. We concentrate on systems with sizes that make them accessible to numerical exact diagonalization. We thus characterize the deviations from ensemble equivalence at finite sizes. We describe multiple ways to compute the microcanonical entropy and present numerical results for the entropy and temperature computed in these various ways. We show that using an energy window whose width has a particular energy dependence results in a temperature with minimal deviations from the canonical temperature.

Many-Body Interference at the Onset of Chaos

We unveil the signature of many-body interference across dynamical regimes of the Bose-Hubbard model. Increasing the particles' indistinguishability enhances the temporal fluctuations of few-body observables, with a dramatic amplification at the onset of quantum chaos. By resolving the exchange symmetries of partially distinguishable particles, we explain this amplification as the fingerprint of the initial state's coherences in the eigenbasis.

Assigning temperatures to eigenstates

In the study of thermalization in finite isolated quantum systems, an inescapable issue is the definition of temperature. We examine and compare different possible ways of assigning temperatures to energies or equivalently to eigenstates in such systems. A commonly used assignment of temperature in the context of thermalization is based on the canonical energy-temperature relationship, which depends only on energy eigenvalues and not on the structure of eigenstates. For eigenstates, we consider defining temperature by minimizing the distance between (full or reduced) eigenstate density matrices and canonical density matrices. We show that for full eigenstates, the minimizing temperature depends on the distance measure chosen and matches the canonical temperature for the trace distance; however, the two matrices are not close. With reduced density matrices, the minimizing temperature has fluctuations that scale with subsystem and system size but appears to be independent of distance measure. In particular limits, the two matrices become equivalent while the temperature tends to the canonical temperature.

Quantum coherence controls the nature of equilibration and thermalization in coupled chaotic systems

A bipartite system whose subsystems are fully quantum chaotic and coupled by a perturbative interaction with a tunable strength is a paradigmatic model for investigating how isolated quantum systems relax toward an equilibrium. It is found that quantum coherence of the initial product states in the energy eigenbasis of the subsystems-quantified by the off-diagonal elements of the subsystem density matrices-can be viewed as a resource for equilibration and thermalization as manifested by the entanglement generated. Results are given for four distinct perturbation strength regimes, the ultraweak, weak, intermediate, and strong regimes. For each, three types of tensor product states are considered for the initial state: uniform superpositions, random superpositions, and individual subsystem eigenstates. A universal timescale is identified involving the interaction strength parameter. In particular, maximally coherent initial product states (a form of uniform superpositions) thermalize under time evolution for any perturbation strength in spite of the fact that in the ultraweak perturbative regime the underlying eigenstates of the system have a tensor product structure and are not at all thermal-like; though the time taken to thermalize tends to infinity as the interaction vanishes. Moreover, it is shown that in the ultraweak regime the initial entanglement growth of the system whose initial states are maximally coherent is quadratic-in-time, in contrast to the widely observed linear behavior.

Exploring Integrability-Chaos Transition with a Sequence of Independent Perturbations

A gas of interacting particles is a paradigmatic example of chaotic systems. It is shown here that, even if all but one particle are fixed in generic positions, the excited states of the moving particle are chaotic. They are characterized by the number of principal components (NPC)-the number of integrable system eigenstates involved into the nonintegrable one, which increases linearly with the number of strong scatterers. This rule is a particular case of the general effect of an additional perturbation on the system chaotic properties. The perturbation independence criteria supposing the system chaoticity increase are derived here as well. The effect can be observed in experiments with photons or cold atoms as the decay of observable fluctuation variance, which is inversely proportional to NPC and, therefore, to the number of scatterers. This decay indicates that the eigenstate thermalization is approached. The results are confirmed by numerical calculations for a harmonic waveguide with zero-range scatterers along its axis.

Eigenstate thermalization hypothesis in two-dimensional XXZ model with or without SU(2) symmetry

We investigate the eigenstate thermalization properties of the spin-1/2 XXZ model in two-dimensional rectangular lattices of size L_{1}×L_{2} under periodic boundary conditions. Exploiting the symmetry property, we can perform an exact diagonalization study of the energy eigenvalues up to system size 4×7 and of the energy eigenstates up to 4×6. Numerical analysis of the Hamiltonian eigenvalue spectrum and matrix elements of an observable in the Hamiltonian eigenstate basis supports that the two-dimensional XXZ model follows the eigenstate thermalization hypothesis. When the spin interaction is isotropic, the XXZ model Hamiltonian conserves the total spin and has SU(2) symmetry. We show that the eigenstate thermalization hypothesis is still valid within each subspace where the total spin is a good quantum number.

Preparing random states and benchmarking with many-body quantum chaos

Producing quantum states at random has become increasingly important in modern quantum science, with applications being both theoretical and practical. In particular, ensembles of such randomly distributed, but pure, quantum states underlie our understanding of complexity in quantum circuits¹ and black holes², and have been used for benchmarking quantum devices^3,4 in tests of quantum advantage^5,6. However, creating random ensembles has necessitated a high degree of spatio-temporal control^7-12 placing such studies out of reach for a wide class of quantum systems. Here we solve this problem by predicting and experimentally observing the emergence of random state ensembles naturally under time-independent Hamiltonian dynamics, which we use to implement an efficient, widely applicable benchmarking protocol. The observed random ensembles emerge from projective measurements and are intimately linked to universal correlations built up between subsystems of a larger quantum system, offering new insights into quantum thermalization¹³. Predicated on this discovery, we develop a fidelity estimation scheme, which we demonstrate for a Rydberg quantum simulator with up to 25 atoms using fewer than 10⁴ experimental samples. This method has broad applicability, as we demonstrate for Hamiltonian parameter estimation, target-state generation benchmarking, and comparison of analogue and digital quantum devices. Our work has implications for understanding randomness in quantum dynamics¹⁴ and enables applications of this concept in a much wider context^{4,5,9,10,15-20}.