Probing entanglement in a 2D hard-core Bose–Hubbard lattice
Amir H. Karamlou, Ilan T. Rosen, Sarah E. Muschinske, Cora N. Barrett, Agustín Di Paolo, Leon Ding, P. M. Harrington, Max Hays, Rabindra Nath Das, David K. Kim, Bethany M. Niedzielski, Meghan Schuldt, Kyle Serniak, Mollie E. Schwartz, Jonilyn Yoder, Simon Gustavsson, Yariv Yanay, Jeffrey A. Grover, William D. Oliver
Abstract Entanglement and its propagation are central to understanding many physical properties of quantum systems 1–3 . Notably, within closed quantum many-body systems, entanglement is believed to yield emergent thermodynamic behaviour 4–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 systems 8–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.
Abstract Entanglement and its propagation are central to understanding many physical properties of quantum systems 1–3 . Notably, within closed quantum many-body systems, entanglement is believed to yield emergent thermodynamic behaviour 4–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 systems 8–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.
Abstract The pursuit of superconducting-based quantum computers has advanced the fabrication of and experimentation with custom lattices of qubits and resonators. Here, we describe a roadmap to use present experimental capabilities to simulate an interacting many-body system of bosons and measure quantities that are exponentially difficult to calculate numerically. We focus on the two-dimensional hard-core Bose–Hubbard model implemented as an array of floating transmon qubits. We describe a control scheme for such a lattice that can perform individual qubit readout and show how the scheme enables the preparation of a highly excited many-body state, in contrast with atomic implementations restricted to the ground state or thermal equilibrium. We discuss what observables could be accessed and how they could be used to better understand the properties of many-body systems, including the observation of the transition of eigenstate entanglement entropy scaling from area-law behavior to volume-law behavior.
We describe an algorithm for studying the entanglement entropy and spectrum of two-dimensional (2D) systems, as a coupled array of N one-dimensional chains in their continuum limit. Using the algorithm to study the quantum Ising model in 2D (both in its disordered phase and near criticality), we confirm the existence of an area law for the entanglement entropy and show that near criticality there is an additive piece scaling as cefflog(N)/6 with ${c}_{\mathrm{eff}}\ensuremath{\approx}1$. Studying the entanglement spectrum, we show that entanglement gap scaling can be used to detect the critical point of the 2D model. When short-range (area law) entanglement dominates we find (numerically and perturbatively) that this spectrum reflects the energy spectrum of a single quantum Ising chain.
The area law relates the bipartite entanglement entropy of a quantum many-body ground state to the size of the boundary between the subsystems, but the geometry of this boundary is rarely discussed. We inspect this in the 1D Bose-Hubbard model at fixed density by comparing four spatial bipartitions of the periodic lattice: first half, second half, even sites, and odd sites; sharing the same number of sites but differing in how the boundary is arranged. We find analytical limits with perturbation theory: in the Mott insulator the contiguous cut obeys the area lay while the alternating cut obeys a volume law S∝Ns, in this sense an anomaly, and for the superfluid both cuts colapse to the binomial saturation due to delocalization of the state. We formulate these limits as a statement about the many-body problem using a generalized slave-boson approach based on mean-field with quantum fluctuations while verifying with Exact Diagonalization (ED) for small lattice sizes and Densitiy Matrix Renormalization Group (DMRG) simulations for Ns≫1. The slave-boson Gaussian ground state allows to compute the entanglement entropy from a reduced correlation matrix for any desired bipartition consistent with ED and DMRG results. Using slave bosons the computational cost is set by the local cutoff nmax rather than the Hilbert space dimension, so we can reach lattice sizes far beyond ED. Our method is capable of establishing the partition-dependent scaling laws as a many-body feature, not only a finite-size effect, in great agreement with the ED for Ns∈[4,10] and DMRG for larger lattice sizes.
In recent years, neural quantum states have emerged as a powerful variational approach, achieving state-of-the-art accuracy when representing the ground-state wave function of a great variety of quantum many-body systems, including spin lattices, interacting fermions or continuous-variable systems. However, accurate neural representations of the ground state of interacting bosons on a lattice have remained elusive. We introduce a neural backflow Jastrow Ansatz, in which occupation factors are dressed with translationally equivariant many-body features generated by a deep neural network. We show that this neural quantum state is able to faithfully represent the ground state of the 2D Bose-Hubbard Hamiltonian across all values of the interaction strength. We scale our simulations to lattices of dimension up to 20 × 20 while achieving the best variational energies reported for this model. This enables us to investigate the scaling of the entanglement entropy across the superfluid-to-Mott quantum phase transition, a quantity hard to extract with non-variational approaches.
While the eigenstate entanglement entropy has been extensively studied for fermionic systems, much less is known about bosonic systems. Here, we study the entanglement entropy of mid-spectrum eigenstates of Bose-Hubbard models, focusing on weakly disordered models with and without particle-number conservation, and contrasting them with the translationally-invariant model. We analyze the volume-law and O(1) contributions to the entanglement entropy via the averages over mid-spectrum eigenstates and the corresponding distributions. We derive the volume-law coefficient of the entanglement entropy by generalizing the mean-field approach from [Phys. Rev. Lett. 119, 220603 (2017)] to many-body systems with a tunable local bosonic cutoff, which agrees with previous analytical and numerical results from [Phys. Rev. B 110, 235154 (2024)]. We show that the volume-law contribution to the entanglement entropy does not change upon breaking translational invariance via on-site disorder. We then numerically study the role of the subleading O(1) contribution to the entanglement entropy. We find that, in the particle-number conserving case, it exhibits a nontrivial dependence on the particle-number density and the local bosonic cutoff, while without particle-number conservation, results suggest the emergence of a universal O(1) contribution beyond the random pure state predictions.
Interacting many-body quantum systems and their dynamics, while fundamental to modern science and technology, are formidable to simulate and understand. However, by discovering their symmetries, conservation laws, and integrability, one can unravel their intricacies. Here, using up to 124 qubits of a fully programmable quantum computer, we uncover local conservation laws and integrability in one- and two-dimensional periodically-driven spin lattices in a regime previously inaccessible to such detailed analysis. We focus on the paradigmatic example of disorder-induced ergodicity breaking, where we first benchmark the system crossover into a localized regime through anomalies in the one-particle-density-matrix spectrum and other hallmark signatures. We then demonstrate that this regime stems from hidden local integrals of motion by faithfully reconstructing their quantum operators, thus providing a more detailed portrait of the system’s integrable dynamics. Our results demonstrate a versatile strategy for extracting the hidden dynamical structure from noisy experiments on large-scale quantum computers. Local integrals of motion are useful for understanding emergent integrability in many-body localized systems. Here the authors use a large-scale superconducting quantum processor with up to 124 qubits to simulate many-body dynamics in 1D and 2D systems and demonstrate extraction of local integrals of motion.