Recent advances in quantum simulators allow direct experimental access to ensembles of pure states generated by measuring part of an isolated quantum many-body system. These projected ensembles encode fine-grained information beyond thermal expectation values and provide a new window into quantum thermalization. In chaotic dynamics, projected ensembles exhibit universal statistics governed by maximum-entropy principles, known as deep thermalization. At infinite temperature, this universality is characterized by Haar-random ensembles. More generally, physical constraints such as finite temperature or conservation laws lead to Scrooge ensembles, which are maximally entropic distributions of pure states consistent with these constraints. Here, we introduce Scrooge $k$-designs, which approximate Scrooge ensembles, and we use this framework to sharpen the conditions under which Scrooge-like behavior emerges. We first show that global Scrooge designs arise from long-time chaotic unitary dynamics alone, without measurements. Second, we show that measuring a complementary subsystem of a scrambled global state drawn from a global Scrooge $2k$-design induces a local Scrooge $k$-design. Third, we show that a local Scrooge $k$-design arises from an arbitrary entangled state when the complementary system is measured in a scrambled basis induced by a unitary drawn from a Haar $2k$-design. These results show that the resources required to generate approximate Scrooge ensembles scale only with the desired degree of approximation, enabling efficient implementations. Complementing our analytical results, numerical simulations identify coherence, entanglement, nonstabilizerness, and information scrambling as essential ingredients for the emergence of local Scrooge-like behavior. Together, our findings advance theoretical explanations for maximally entropic, information-stingy randomness in quantum many-body systems.

For closed quantum systems, topological orders are understood through the equivalence classes of ground states of gapped local Hamiltonians. The generalization of this conceptual paradigm to open quantum systems, however, remains elusive, often relying on operational definitions without fundamental principles. Here, we fill this gap by proposing an approach based on three axioms: (i) local recoverability, (ii) absence of long-range correlations, and (iii) spatial uniformity. States that satisfy these axioms are fixed points; requiring the axioms only after coarse-graining promotes each fixed point to an equivalence class, i.e., a phase, presenting the first step toward the axiomatic classification of mixed-state phases of matter: the mixed-state bootstrap program. From these axioms, a rich set of topological data naturally emerges; importantly, these data are robust under relaxation of axioms. For example, each topological mixed state supports locally indistinguishable classical and/or quantum logical memories with distinct responses to topological operations. These data label distinct mixed-state phases, allowing one to distinguish them. We further uncover a hierarchy of secret-sharing constraints: In non-Abelian phases, reliable recovery—even of information that looks purely classical—demands a specific coordination among spatial subregions, a requirement different across non-Abelian classes. This originates from non-Abelian fusion rules that can stay robust under decoherence. Finally, we performed large-scale numerical simulations to corroborate stability: Weakly decohered fixed points respect the axioms once coarse-grained. Surprisingly, the axiom violations at different length scales exhibit an exact crossing at the critical point, enabling its identification without any rescaling or data-collapse ambiguities. These results lay the foundation for a systematic classification of topological states in open quantum systems.

We introduce a variational method for simulating the dynamics of interacting open quantum bosonic systems deep in the quantum regime. The method is based on a multidimensional Wigner phase-space representation and employs a variational multi-Gaussian Ansatz, whose accuracy is systematically controlled by the number of Gaussian components. The variational equations of motion are derived from the Dirac-Frenkel principle and evaluated efficiently by combining the analytical structure of Gaussian functions with automatic differentiation. As a key physical application, we study a driven-dissipative two-dimensional Bose-Hubbard lattice with two-boson coherent driving and two-body losses. Using our dynamical approach, we compute the finite-size scaling of the Liouvillian spectral gap—extracted from the relaxation dynamics—which vanishes in the thermodynamic limit. Our results reveal critical slowing down with dynamical exponents of the 2D quantum Ising universality class, demonstrating the power of our method to capture complex quantum dynamics in large open systems.

Fractional quantum Hall states are promising platforms for topological quantum computation due to their capacity to encode quantum information in topologically degenerate ground states and in the fusion space of non-Abelian anyons. We investigate how the information encoded in two paradigmatic states, the Laughlin and Moore-Read states, is affected by density decoherence—coupling of local charge density to nonthermal noise. We identify a critical filling factor ${\nu }_c$, above which the quantum information remains fully recoverable for arbitrarily strong decoherence. The $\nu =1/3$ Laughlin state and $\nu =1/2$ Moore-Read state both lie within this range. Below ${\nu }_c$ both classes of states undergo a decoherence-induced Berezinskii-Kosterlitz-Thouless (BKT) transition into a critical decohered phase. For Laughlin states, information encoded in the topological ground state manifold degrades continuously with decoherence strength inside this critical phase, vanishing only in the limit of infinite decoherence strength. On the other hand, quantum information encoded in the fusion space of non-Abelian anyons of the Moore-Read states remains fully recoverable under arbitrarily strong decoherence even beyond the BKT transition. These results lend further support to the promise of non-Abelian FQH states as platforms for topological quantum computation. They reveal how the interplay between chiral topological order and global symmetry can give rise to novel mixed-state phases under decoherence, distinguished by their capacity for storing quantum information.

By combining the assumptions of Bell locality with those of generalized noncontextuality, we define classes of noncontextuality inequalities for correlations arising in a bipartite Bell circuit. These classes are distinguished by which subsets of the full set of operational identities are taken as input to the principle of noncontextuality; certain natural subsets form a hierarchy that provides a way of understanding and classifying different forms of quantum correlations, including entanglement, steering, and nonlocality. Each level of this hierarchy gives rise to a corresponding class of noncontextuality inequalities whose violation witnesses one of these forms of bipartite quantum resourcefulness, thereby yielding different sufficient conditions for entanglement. The resulting entanglement certification paradigm requires no prior characterization of the measurements, delivers a verdict independent of tomographic gauge freedom, and can certify any entangled state without auxiliary entangled sources. To illustrate the power of this paradigm, we show that noncontextuality inequalities can certify entanglement for families of two-qubit isotropic states for which certification by Bell or steering inequalities is known to be impossible. We also show that, compared with the Bell test, this approach certifies a much larger fraction of entangled states, while the associated membership problem is computationally more tractable. On the experimental side, we describe techniques to ensure nontrivial operational identities in the presence of noisy and imperfect implementations. We also identify the sufficient condition under which these techniques are valid, namely, a particular notion of tomographic completeness of the implemented Bell circuit, which ensures that the operational identities are independent of tomographic gauge. Finally, we provide an experimental demonstration of the superior performance of this entanglement certification technique using polarization-entangled photons.

Quantum symmetrization is the task of transforming a nonstrictly increasing list of $n$ integers into an equal superposition of all permutations of the list (or more generally, performing this operation coherently on a superposition of such lists). This task plays a key role in initial state preparation for first-quantized simulations. Motivated by an application to fermionic systems, various algorithms have been proposed to solve a weaker version of symmetrization in which the input list is strictly increasing, but the general symmetrization problem with repetitions in the input list has not been well studied. We present the first efficient quantum algorithms for the general symmetrization problem. If $m$ is the greatest possible value of the input list, our first algorithm symmetrizes any single classical input list using $\tilde{O}(\log n)$ depth and $O(n\log n+\log m)$ ancilla qubits, and our second algorithm symmetrizes an arbitrary superposition of input lists using $\tilde{O}({\log }^{3}n)$ depth and $O(n\log n)$ ancilla qubits. Our algorithms enable efficient simulation of bosonic quantum systems in first quantization (in particular, we describe a simulation of the Bose-Hubbard model with complexity $\mathrm{poly}(n,\log m)$) and can prepare (superpositions of) Dicke states of any Hamming weight in $\tilde{O}(\log n)$ depth (respectively, $\tilde{O}({\log }^{3}n)$ depth) using $O(n\log n)$ ancilla qubits. We also propose an $\tilde{O}({\log }^{3}n)$-depth quantum algorithm to transform second-quantized states to first-quantized states. Using this algorithm, QFT-based quantum telescope arrays can image brighter photon sources, extending quantum interferometric imaging systems to a new regime.

Quantum routers (QRouters) are essential components of bucket-brigade quantum random access memory (QRAM), enabling quantum applications such as Grover’s search and quantum machine learning. Despite significant theoretical advances, achieving scalable and coherent QRouters experimentally remains challenging. Here, we demonstrate coherent quantum routers using a superconducting quantum processor, laying a practical foundation for scalable QRAM systems. The quantum router at the core of our implementation utilizes the transition composite gate scheme, wherein auxiliary energy levels temporarily mediate conditional interactions. By leveraging two-qubit primitives in its construction, the quantum router preserves straightforward initialization and calibration while significantly reducing circuit depth and duration relative to traditional gate-based decompositions. Moreover, by encoding routing addresses in the nonadjacent qutrit states $|0⟩$ and $|2⟩$, our design inherently enables erasure-detection capability, providing efficient postselection to mitigate routing errors. Experimentally, we achieve an average fidelity of 94.8% across three individual QRouters, and we validate scalability through a two-layer quantum routing network achieving a fidelity of 82.4%. Our results represent a significant advancement in quantum routing technology, providing enhanced fidelity, built-in error resilience, and practical scalability crucial for the development of future QRAM and large-scale quantum computing architectures.

Bounding the correlations predicted by quantum theory is an important challenge in quantum information science. Today’s leading approach is semidefinite programming relaxations, but existing methods still cannot account for many relevant types of constraints. Here, we propose a general semidefinite relaxation methodology that can incorporate a breadth of constraints needed in various quantum correlation problems, thereby generalising the seminal Navascués-Pironio-Acín hierarchy. It yields useful results at reasonable computational cost. We showcase the methodology and its features by using it to address five different quantum information problems. These are (i) entanglement witnessing from imperfect measurement devices, (ii) certifying measurements from fidelity-constrained sources, (iii) computing dimensionality in genuine multiparticle entangled states, (iv) benchmarking dimensionality for state preparation devices, and (v) finding uncertainty relations for nearly anticommuting observables. These applications reflect both the usefulness and versatility of the methodology, as well as its potential for broader relevance in the field.

The generation of non-Gaussian quantum states is a key requirement for universal continuous-variable quantum information processing. We report the experimental generation of large-amplitude squeezed coherent-state superpositions (squeezed cat states) on free-space optical pulses, reaching an amplitude of $\alpha =2.47$, which, to our knowledge, exceeds all previously reported values. Our protocol relies on the controlled mixing of the Fock states $|1⟩$ and $|2⟩$ through a tunable beam splitter, followed by heralding via homodyne detection. The resulting state displays three well-resolved negative regions in its Wigner function and achieves a fidelity of 0.53 with the target state $\propto \widehat{S}(z)(|\alpha ⟩-|-\alpha ⟩)$, with $\alpha =2.47$ and squeezing parameter $z=0.56$. These results constitute a significant milestone for temporal breeding protocols and for the iterative generation of optical Gottesman-Kitaev-Preskill states, opening perspectives for scalable and fault-tolerant photonic quantum architectures.

Analog quantum simulators offer a route to exploring strongly correlated many-body dynamics beyond classical computation, but their predictive power remains limited by the absence of quantitative error estimation. Establishing rigorous uncertainty bounds is essential for elevating such devices from qualitative demonstrations to quantitative scientific tools. Here we introduce a general framework for bounded-error quantum simulation, which provides predictions for many-body observables with experimentally quantifiable uncertainties. The approach combines Hamiltonian and Lindbladian Learning—a statistically rigorous inference of the coherent and dissipative generators governing the dynamics—with the propagation of their uncertainties into the simulated observables, yielding confidence bounds directly derived from experimental data. We demonstrate this framework on trapped-ion quantum simulators implementing long-range Ising interactions with up to 51 ions. We analyze error bounds on two levels. First, we learn an open-system model from experimental data collected in an initial time window of quench dynamics, simulate the corresponding master equation, and quantitatively verify consistency between theoretical predictions and measured dynamics at long times. Second, we discuss short-time error bounds directly from experimental measurements alone, without relying on classical simulation—crucial for entering regimes of quantum advantage. In both cases, the learned models reproduce the experimental evolution within the predicted bounds, demonstrating quantitative reliability and internal consistency. By integrating statistical learning, open-system modeling, and precise experimental control, bounded-error quantum simulation provides a scalable foundation for trusted analog quantum computation, bridging the gap between experimental quantum platforms and predictive many-body physics. The techniques presented here directly extend to digital quantum simulation.

We construct a family of static, geometrically local Hamiltonians that inherit eigenstate properties of periodically driven (Floquet) systems. Our construction is a variation of the Feynman-Kitaev clock—a well-known mapping between quantum circuits and local Hamiltonians—where the clock register is given periodic boundary conditions. Assuming the eigenstate thermalization hypothesis (ETH) holds for the input circuit, our construction yields Hamiltonians whose eigenstates have properties characteristic of infinite temperature, like volume-law entanglement entropy, across the whole spectrum—including the ground state. We then construct a family of exactly solvable Floquet quantum circuits whose eigenstates are shown to obey the ETH at infinite temperature. Combining the two constructions yields a new family of local Hamiltonians with provably volume-law entangled ground states, and the first such construction where the volume law holds for all contiguous subsystems.

Quantifying physical concepts in terms of the ultimate performance of a given task has been central to theoretical progress, as illustrated by thermodynamic entropy and entanglement entropy, which quantify irreversibility and quantum correlations, respectively. Symmetry breaking is equally universal yet lacks such an operational quantification. While an operational characterization of symmetry breaking through asymptotic state-conversion efficiency is a central goal of the resource theory of asymmetry (RTA), such a characterization has so far been completed only for the U(1) group among continuous symmetries. Here, we identify the complete measure of symmetry breaking for a general continuous symmetry described by any compact Lie group. Specifically, we show that the asymptotic conversion rate between many copies of pure states in RTA is determined by the quantum geometric tensor, thereby establishing it as the complete measure of symmetry breaking. As an immediate consequence of our conversion rate formula, we also resolve the Marvian-Spekkens conjecture on conditions for reversible conversion in RTA, which has remained unproven for over a decade. Leveraging the connection between symmetry breaking and the theory of quantum reference frames, we also systematically introduce a standardized reference state for frameness based on our asymptotic conversion theory. In addition, by applying our analysis to a standard quantum-thermodynamic scenario, we show that asymptotic state conversion in contact with heat baths generally requires macroscopic coherence in the thermodynamic limit.

We introduce configuration space path integrals for quantum fields interacting with classical fields. We show that this can be done consistently by proving that the dynamics are completely positive directly, without resorting to master equation methods. These path integrals allow one to readily impose space-time symmetries, including Lorentz invariance or diffeomorphism invariance. They generalize and combine the Feynman-Vernon path integral of open quantum systems and the stochastic path integral of classical stochastic dynamics while respecting symmetry principles. We introduce a path integral formulation of general relativity where the space-time metric is treated classically. The theory is a candidate for a fundamental theory that reconciles general relativity with quantum mechanics. The theory is manifestly covariant, and may be inequivalent to the theory derived using master equation methods. We prove that entanglement cannot be created via the classical field, reinforcing proposals to test the quantum nature of gravity via entanglement generation.

The distribution of entanglement across distant qubits is a central challenge for the operation of scalable quantum computers and large-scale quantum networks. Existing approaches rely on deterministic state transfer, or probabilistic protocols that require active control or measurements and postselection. Here, we demonstrate a fundamentally different, fully autonomous process, where two remote qubits are entangled through their coupling to a quantum-correlated photonic reservoir. In our experiment, a Josephson parametric converter produces a Gaussian, continuous-variable entangled state of propagating microwave fields that drives two spatially separated superconducting transmon qubits into a stationary, discrete-variable entangled state. We also show how qubit tomography unlocks a direct and sensitive verification of two-mode squeezing in the microwave domain. These results establish networks of qubits interfaced with distributed continuous-variable entangled states as a powerful platform for foundational studies and quantum-technology applications.

Remote entanglement between widely separated qubits is a fundamental quantum phenomenon and a critical resource for quantum information applications. Generating entanglement between independent qubits separated by arbitrary, potentially large distances requires propagating quantum states and is typically achieved using pulsed protocols combining distinct steps of local entanglement generation followed by distribution. This necessity raises an intriguing question: Can remote entanglement be stabilized indefinitely, instead of periodically regenerated and redistributed only after decay? Here, we demonstrate that this is indeed possible, reporting autonomous stabilization of entanglement between two separate superconducting-qubit devices. Combining nonreciprocal waveguide coupling and local driving, we experimentally realize a symmetry-based coherent quantum-absorber scheme in a cascaded network. We quantify the degree of entanglement through quantum state tomography, finding that the protocol’s entangling power is severely limited by imperfections that break the required symmetry. We show, however, that a modified protocol based on an alternate symmetry is far more robust, enabling us to achieve a concurrence approaching 0.5, a limit set by only local loss in the network. Our results enable on-demand delivery of high-fidelity entanglement in modular quantum processors and networks and pave the way for autonomously protecting distributed quantum information.

The quest to identify quantum advantages lies at the heart of quantum technology. While quantum devices promise extraordinary capabilities, from exponential computational speedups to unprecedented measurement precision, distinguishing genuine advantages from mere illusions remains a formidable challenge. In this endeavor, quantum theorists are like prophets attempting to foretell the future, yet the boundary between visionary insight and unfounded fantasy is perilously thin. In this Perspective, we examine our mathematical tools for navigating the vast world of quantum advantages across computation, learning, sensing, and communication. We explore five keystone properties, predictability, typicality, robustness, verifiability, and usefulness that define an ideal quantum advantage, and envision what new quantum advantages could arise in a future with ubiquitous quantum technology. We prove that some quantum advantages are inherently unpredictable using classical resources alone, suggesting a landscape far richer than what we can currently foresee. While mathematical rigor remains our indispensable guide, the ultimate power of quantum technologies may emerge from advantages we cannot yet conceive.

We show how to perform scalable fault-tolerant non-Clifford gates in two dimensions by introducing domain walls between the surface code and a non-Abelian topological code whose code space is stabilized by Clifford operators. We formulate a path-integral framework which provides both a macroscopic picture for different logical gates as well as a way to derive the associated microscopic circuits. We also show an equivalence between our approach and prior proposals where a 2D array of qubits reproduces the action of a transversal gate in a 3D stabilizer code over time, thus establishing a connection between 3D codes and 2D non-Abelian topological phases. We prove a threshold theorem for our protocols under local stochastic circuit noise using a just-in-time decoder to correct the non-Abelian code.

BosonSampling is the leading candidate for demonstrating quantum computational advantage in photonic systems. While we have recently seen many impressive experimental demonstrations, there is still a formidable distance between the complexity-theoretic hardness arguments and current experiments. One of the largest gaps involves the ratio of particles to modes—all current hardness evidence assumes a “dilute” regime in which the number of linear optical modes scales at least quadratically in the number of particles. By contrast, current experiments operate in a “saturated” regime with a linear number of modes. In this paper, we bridge this gap, bringing the hardness evidence for experiments in the saturated regime to the same level as had been previously established for the dilute regime. This involves proving a worst-to-average-case reduction for computing the permanent, which is robust both to large numbers of row repetitions and also to distributions over matrices with correlated entries. We also apply similar arguments to give evidence for hardness of Gaussian BosonSampling in the saturated regime.

Nonequilibrium dynamics of quantum many-body systems is challenging for classical computing, providing opportunities for demonstrating practical quantum computational advantage with analog quantum simulators. Owing to the intimate connection with a random matrix ensemble, it is proposed to be classically intractable to sample the driven thermalized many-body states of a Bose-Hubbard system and further extract multipoint correlations from the output strings for characterizing quantum systems. Here, leveraging dedicated precise manipulations and atom-number-resolved detection through a quantum gas microscope with bichromatic superlattices, we perform sampling of the driven Hubbard chains and two-leg ladders in the thermalized phase involving up to 64 sites with 20 atoms, yielding a Hilbert space dimension of ${10}^{19}$ and outpacing the most powerful supercomputer in terms of sampling rate by 3 orders of magnitude. The volume law scaling of the Rényi entanglement entropy in the thermalized phase is observed, which hinders efficient classical simulation for large systems. We employ the Bayesian tests to verify that our prepared systems operate in the driven thermalized phase. Multipoint correlations of up to 14th-order extracted from the experimental samples offer clear distinctions between the thermalized and many-body-localized phases, where classical computations such as tensor network fail to give accurate and faithful predictions within a reasonable time cost. Our work demonstrates the sampling of an interacting chaotic system performed on a quantum processor of ultracold atoms and opens the door of utilizable quantum computational advantage in simulating Floquet dynamics of many-body systems.

There is a received wisdom about where to draw the boundary between classical and nonclassical for various types of quantum processes. For multipartite states, it is the divide between separable and entangled; for channels, the divide between entanglement breaking and not; for sets of measurements, the divide between compatible and incompatible; for assemblages, the divide between unsteerable and steerable. However, these choices have not been motivated by any unified notion of what it means to be classically explainable. One well-motivated notion of classical explainability is the one based on generalized noncontextuality: A set of circuits is classically explainable if a generalized-noncontextual ontological model can realize the statistics they generate. In this work, we show that this notion can be leveraged to define a classical-nonclassical divide for individual quantum processes of arbitrary type. A set of measurements, for example, is judged to be classical if and only if a particular set of circuits—the one obtained by contracting these measurements with every possible quantum state—is classically explainable in the sense just articulated. We begin the task of characterizing where the classical-nonclassical divide lies according to this proposal for a variety of different types of processes. In particular, we show that all of the following are judged to be nonclassical: every entangled state, every set of incompatible measurements, every non-entanglement-breaking channel, and every steerable assemblage. Our proposal differs from the received wisdom, however, insofar as it also judges certain subsets of the complementary classes to be nonclassical, including certain separable states, compatible sets of measurements, entanglement-breaking channels, and unsteerable assemblages. Finally, we prove structure theorems characterizing the classical-nonclassical divide based on whether a process admits of a specific type of frame representation.