We test the eigenstate thermalization hypothesis (ETH) in 1+1-dimensional SU(2) lattice gauge theory (LGT) with one flavor of dynamical fermions. Using the loop-string-hadron framework of the LGT with a bosonic cut-off, we exactly diagonalize the Hamiltonian for finite size systems and calculate matrix elements (MEs) in the eigenbasis for both local and non-local operators. We analyze different indicators to identify the parameter space for quantum chaos at finite lattice sizes and investigate how the ETH behavior emerges in both the diagonal and off-diagonal MEs. Our investigations allow us to study various time scales of thermalization and the emergence of random matrix behavior, and highlight the interplays of the several diagnostics with each other. Furthermore, from the off-diagonal MEs, we extract a smooth function that is closely related to the spectral function for both local and non-local operators. We find numerical evidence of the spectral gap and the memory peak in the non-local operator case. Finally, we investigate aspects of subsystem ETH in the lattice gauge theory and identify certain features in the subsystem reduced density matrix that are unique to gauge theories.

Accurately determining ground-state properties of quantum many-body systems remains one of the major challenges of quantum simulation. In this work, we present a protocol for estimating the ground-state energy using only global time evolution under a target Hamiltonian. This avoids the need for controlled operations that are typically required in conventional quantum phase estimation and extends the algorithm applicability to analog simulators. Our method extracts energy differences from measurements of the Loschmidt echo over an initial ground-state approximation, combines them with direct energy measurements, and solves a set of equations to infer the individual eigenenergies. We benchmark this protocol on free-fermion systems, showing orders-of-magnitude precision gains over direct energy measurements on the initial state, with accuracy improving rapidly with initial-state fidelity and persisting for hundreds of modes. We further demonstrate applicability to the 2D Ising and Fermi-Hubbard models and show that the approach extends naturally to other observables such as order parameters. Finally, we analyze the effect of experimental imperfections and propose error-mitigation strategies. These results establish a practical route to compute physically relevant quantities with high precision using globally controlled quantum simulators.

Matrix-product unitaries (MPUs) are many-body unitary operators that, as a consequence of their tensor-network structure, preserve the entanglement area law in 1D systems. However, it is unknown how to implement an MPU as a quantum circuit since the individual tensors describing the MPU are not unitary. In this paper, we show that a large class of MPUs can be implemented with a polynomial-depth quantum circuit. For an $N$-site MPU built from a repeated bulk tensor with open boundary, we explicitly construct a quantum circuit of polynomial depth $T = O(N^{\alpha})$ realizing the MPU, where the constant $\alpha$ depends only on the bulk and boundary tensor and not the system size $N$. We show that this class includes nontrivial unitaries that generate long-range entanglement and, in particular, contains a large class of unitaries constructed from representations of $C^*$-weak Hopf algebras. Furthermore, we also adapt our construction to nonuniform translationally-varying MPUs and show that they can be implemented by a circuit of depth $O(N^{\beta} \, \mathrm{poly}\, D)$ where $\beta \le 1 + \log_2 \sqrt{D}/ s_{\min}$, with $D$ being the bond dimension and $s_{\min}$ is the smallest nonzero Schmidt value of the normalized Choi state corresponding to the MPU.

We propose applying the adiabatic algorithm to prepare high-energy eigenstates of integrable models on a quantum computer. We first review the standard adiabatic algorithm to prepare ground states in each magnetization sector of the prototypical XXZ Heisenberg chain. Based on the thermodynamic Bethe ansatz, we show that the algorithm circuit depth is polynomial in the number of qubits $N$, outperforming previous methods explicitly relying on integrability. Next, we propose a protocol to prepare arbitrary eigenstates of integrable models that satisfy certain conditions. For a given target eigenstate, we construct a suitable parent Hamiltonian written in terms of a complete set of local conserved quantities. We propose using such Hamiltonian as an input for an adiabatic algorithm. After benchmarking this construction in the case of the non-interacting XY spin chain, where we can rigorously prove its efficiency, we apply it to prepare arbitrary eigenstates of the Richardson-Gaudin models. In this case, we provide numerical evidence that the circuit depth of our algorithm is polynomial in $N$ for all eigenstates, despite the models being interacting.

In this manuscript, we explore the intersection of QML and TN in the context of the one-dimensional ANNNI model with a transverse field. The study aims to concretely connect QML and TN by combining them in various stages of algorithm construction, focusing on phase diagram reconstruction for the ANNNI model, with supervised and unsupervised techniques. The model's significance lies in its representation of quantum fluctuations and frustrated exchange interactions, making it a paradigm for studying magnetic ordering, frustration, and the presence of a floating phase. It concludes with discussions of the results, including insights from increased system sizes and considerations for future work, such as addressing limitations in QCNN and exploring more realistic implementations of QC.

Rydberg-cavity systems are emerging as promising platforms for quantum simulation and quantum information processing. These hybrid architectures combine two complementary interaction mechanisms: cavity photons mediate collective long-range couplings, while Rydberg excitations generate strong short-range interactions. Together, they offer a setting for engineering many-body phases characterized by a hierarchy of interactions across widely different length scales. In this work, we introduce a minimal and scalable model for such systems. Focusing on the strong Rydberg blockade regime, we restrict the Hilbert space to the subspace enforced by the blockade, yielding a kinetically constrained long-range model in one spatial dimension. This approach both captures the physics of Rydberg-cavity experiments in the regime of strong Rydberg interactions and provides a conceptually transparent framework for studying the interplay of long-range and short-range interactions. At equilibrium, in addition to paramagnetic and Néel-ordered phases, the system supports a blockaded ferromagnetic/superradiant phase, distinct from the conventional superradiant phase. Out of equilibrium, we identify long-range quantum many-body scars, which are atypical nonthermal eigenstates that evade the eigenstate thermalization hypothesis, and giving rise to slow entanglement growth. In contrast to the linear-in-time entanglement growth characteristic of short-range scarred models, these long-range scars exhibit logarithmic entanglement dynamics. Our results establish a minimal yet versatile framework for Rydberg-cavity systems, and provide a stepping stone for future theoretical and experimental studies of this frontier platform in quantum many-body physics.

We show how to optimally reduce the local Hilbert basis of lattice quantum many-body (QMB) Hamiltonians. The basis truncation exploits the most relevant eigenvalues of the estimated single-site reduced density matrix (RDM). It is accurate and numerically stable across different model phases, even close to quantum phase transitions. We apply this procedure to different models, such as the Sine-Gordon model, the $\varphi^{4}$ theory, and lattice gauge theories, namely Abelian $\mathrm{U}(1)$ and non-Abelian $\mathrm{SU}(2)$, in one and two spatial dimensions. Our results reduce state-of-the-art estimates of computational resources for classical and quantum simulations.

Communication in quantum networks suffers notoriously from photon loss. Resulting errors can be mitigated with a suitable measurement herald at the receiving node. However, waiting for a herald and communicating the measurement result back to the sender in a repeat-until-success strategy makes the protocol slow and prone to errors from false heralds such as detector dark counts. Here we implement an entanglement herald at the sending node by employing a cascaded two-photon emission of a single atom into two optical fiber cavities: The polarization of one photon is entangled with the spin of the atom, and the second photon heralds entanglement generation. We show that heralding improves the atom-photon entanglement in-fiber efficiency and fidelity to 68(3)% and 87(2)%, respectively. We highlight the potential of our source for noise-limited long-distance quantum communication by extending the range for constant fidelity or, alternatively, increasing the fidelity for a given distance.

We investigate the time dependence of anti-flatness in the entanglement spectrum, a measure for non-stabilizerness and lower bound for non-local quantum magic, on a subsystem of a linear SU(2) plaquette chain during thermalization. Tracing the time evolution of a large number of initial states, we find that the anti-flatness exhibits a barrier-like maximum during the time period when the entanglement entropy of the subsystem grows rapidly from the initial value to the microcanonical entropy. The location of the peak is strongly correlated with the time when the entanglement exhibits the strongest growth. This behavior is found for generic highly excited initial computational basis states and persists for coupling constants across the ergodic regime, revealing a universal structure of the entanglement spectrum during thermalization. We conclude that quantitative simulations of thermalization for nonabelian gauge theories require quantum computing. We speculate that this property generalizes to other quantum chaotic systems.

The Loschmidt echo - the probability of a quantum many-body system to return to its initial state following a dynamical evolution - generally contains key information about a quantum system, relevant across various scientific fields including quantum chaos, quantum many-body physics, or high-energy physics. However, it is typically exponentially small in system size, posing an outstanding challenge for experiments. Here, we experimentally investigate the subsystem Loschmidt echo, a quasi-local observable that captures key features of the Loschmidt echo while being readily accessible experimentally. Utilizing quantum gas microscopy, we study its short- and long-time dynamics. In the short-time regime, we observe a dynamical quantum phase transition arising from genuine higher-order correlations. In the long-time regime, the subsystem Loschmidt echo allows us to quantitatively determine the effective dimension and structure of the accessible Hilbert space in the thermodynamic limit. Performing these measurements in the ergodic regime and in the presence of emergent kinetic constraints, we provide direct experimental evidence for ergodicity breaking due to fragmentation of the Hilbert space. Our results establish the subsystem Loschmidt echo as a novel and powerful tool that allows paradigmatic studies of both non-equilibrium dynamics and equilibrium thermodynamics of quantum many-body systems, applicable to a broad range of quantum simulation and computing platforms.

The high-temperature series expansion for quantum spin models is a well-established tool to compute thermodynamic quantities and equal-time spin correlations, in particular for frustrated interactions. We extend the scope of this expansion to the dynamic Matsubara spin-spin correlator and develop a fully analytic algorithm to compute its expansion coefficients. We focus on Heisenberg models with a single coupling constant J and spin lengths S=1/2,1. The expansion coefficients up to 12th order in J/T are precomputed on all possible ~10^6 graphs embeddable in arbitrary lattices and are provided under this https URL. This enables calculation of static momentum-resolved susceptibilities for arbitrary site-pairs or wavevectors. We test our results for the S=1/2 Heisenberg chain and on the triangular lattice model. Moreover, the analytic frequency dependence in the expansion allows for stable analytic continuation to the real-frequency dynamic structure factor. This important application is discussed in a companion letter.

The interplay between topology and strong interactions gives rise to a variety of exotic quantum phases, including fractional quantum Hall (FQH) states and their lattice analogs - fractional Chern insulators (FCIs). Such topologically ordered states host fractionalized excitations, which for spinful systems are often accompanied by ferromagnetism and skyrmions. Here, we study a Hofstadter-Hubbard model of spinful fermions on a square lattice, extended by nearest-neighbor interactions. Using large-scale density matrix renormalization group (DMRG) simulations, we demonstrate the emergence of a spin-polarized $\frac{1}{3}$-Laughlin-like FCI phase, characterized by a quantized many-body Chern number, a finite charge gap, and hidden off-diagonal long-range order. We further investigate the quantum Hall ferromagnet at $\nu=1$ and its skyrmionic excitations upon doping. In particular, we find that nearest-neighbor repulsion is sufficient to stabilize both particle- and hole-skyrmions in the ground state around $\nu=1$, whereas we do not find such textures around $\nu=\frac{1}{3}$. The diagnostic toolbox presented in this work, based on local densities, correlation functions, and spin-resolved observables, is directly applicable in quantum gas microscopy experiments. Our results open new pathways for experimental exploration of FCIs with spin textures in both ultracold atom and electronic systems.

Parton distribution functions (PDFs) describe universal properties of bound states and allow us to calculate scattering amplitudes in processes with large momentum transfer. Calculating PDFs involves the evaluation of matrix elements with a Wilson line in a light-cone direction. In contrast to Monte Carlo methods in Euclidean spacetime, these matrix elements can be directly calculated in Minkowski-space using the Hamiltonian formalism. The necessary spatial- and time-evolution can be efficiently applied using established tensor network methods. We present PDFs in the Schwinger model calculated with matrix product states.

We study the behavior of magic as a bipartite correlation in the quantum Ising chain across its quantum phase transition, and at finite temperature. In order to quantify the magic of partitions rigorously, we formulate a hybrid scheme that combines stochastic sampling of reduced density matrices via quantum Monte Carlo, with state-of-the-art estimators for the robustness of magic - a {\it bona fide} measure of magic for mixed states. This allows us to compute the mutual robustness of magic for partitions up to 8 sites, embedded into a much larger system. We show how mutual robustness is directly related to critical behaviors: at the critical point, it displays a power law decay as a function of the distance between partitions, whose exponent is related to the partition size. Once finite temperature is included, mutual magic retains its low temperature value up to an effective critical temperature, whose dependence on size is also algebraic. This suggests that magic, differently from entanglement, does not necessarily undergo a sudden death.

Symbolic recovery of differential equations is the ambitious attempt at automating the derivation of governing equations with the use of machine learning techniques. In contrast to classical methods which assume the structure of the equation to be known and focus on the estimation of specific parameters, these algorithms aim to learn the structure and the parameters simultaneously. While the uniqueness and, therefore, the identifiability of parameters of governing equations are a well-addressed problem in the field of parameter estimation, it has not been investigated for symbolic recovery. However, this problem should be even more present in this field since the algorithms aim to cover larger spaces of governing equations. In this paper, we investigate under which conditions a solution of a differential equation does not uniquely determine the equation itself. For various classes of differential equations, we provide both necessary and sufficient conditions for a function to uniquely determine the corresponding differential equation. We then use our results to devise numerical algorithms aiming to determine whether a function solves a differential equation uniquely. Finally, we provide extensive numerical experiments showing that our algorithms can indeed guarantee the uniqueness of the learned governing differential equation, without assuming any knowledge about the analytic form of function, thereby ensuring the reliability of the learned equation.

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 schemes or probabilistic protocols that require active control or measurement and postselection. Here we demonstrate an alternative, 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. Beyond entanglement distribution, we also show that superconducting qubits can be used to directly certify two-mode squeezing, with higher sensitivity and without the need for calibrated noise-subtraction. These results establish networks of qubits interfaced with distributed continuous-variable entangled states as a powerful new platform for both foundational studies and quantum-technology relevant applications.

Out-of-equilibrium phases in many-body systems constitute a new paradigm in quantum matter - they exhibit dynamical properties that may otherwise be forbidden by equilibrium thermodynamics. Among these non-equilibrium phases are periodically driven (Floquet) systems [1-5], which are generically difficult to simulate classically because of their high entanglement. Here we realize a Floquet topologically ordered state theoretically proposed in ref. [6], on an array of superconducting qubits. We image the characteristic dynamics of its chiral edge modes and characterize its emergent anyonic excitations. Devising an interferometric algorithm allows us to introduce and measure a bulk topological invariant to probe the dynamical transmutation of anyons for system sizes up to 58 qubits. Our work demonstrates that quantum processors can provide key insights into the thus-far largely unexplored landscape of highly entangled non-equilibrium phases of matter.

Fermionic quantum processors are a promising platform for quantum simulation of correlated fermionic matter. In this work, we study a hardware-efficient protocol for measuring complex expectation values of the time-evolution operator, commonly referred to as Loschmidt echoes, with fermions in an optical superlattice. We analyze the algorithm for the Fermi-Hubbard model at half-filling as well as at finite doping. The method relies on global quench dynamics and short imaginary time evolution, the latter being realized by architecture-tailored pulse sequences starting from a product state of plaquettes. Our numerical results show that complex Loschmidt echoes can be efficiently obtained for large many-body states over a broad spectral range. This allows one to measure spectral properties of the Fermi-Hubbard model, such as the local density of states, and paves the way for the study of finite-temperature properties in current fermionic quantum simulators.

Quantum spin liquids (QSLs) are long-range entangled phases of frustrated magnets exhibiting fractionalized spin excitations. In two dimensions, there is limited analytical understanding of their excitation spectra beyond parton mean-field theories, which fail to capture many features of the finite frequency dynamical response from recent experimental and numerical works. We use a self-consistent random phase approximation (RPA) for the $J_1$-$J_2$ Heiseneberg model on the triangular lattice to describe the strong spinon-spinon interactions of the U(1) Dirac QSL. We obtain quantitative results for the dynamical spin structure factor and phase diagram compatible with comprehensive numerical efforts. We extend the method to chiral QSLs, and discuss its broad range of applicability to other models and for describing inelastic neutron scattering experiments.