Efficient Tensor Network Ansatz for High-Dimensional Quantum Many-Body Problems

Anomalous suppression of large-scale density fluctuations in classical and quantum spin liquids

Classical spin liquids (CSLs) are intriguing states of matter that do not exhibit long-range magnetic order and are characterized by an extensive ground-state degeneracy. Adding quantum fluctuations, which induce dynamics between these different classical ground states, can give rise to quantum spin liquids (QSLs). QSLs are highly entangled quantum phases of matter characterized by fascinating emergent properties, such as fractionalized excitations and topological order. One such exotic quantum liquid is the [Formula: see text] QSL, which can be regarded as a resonating valence bond (RVB) state formed from superpositions of dimer coverings of an underlying lattice. In this work, we unveil a hidden large-scale structural property of archetypal CSLs and QSLs known as hyperuniformity, i.e., normalized infinite-wavelength density fluctuations are completely suppressed in these systems. In particular, we first demonstrate that classical ensembles of close-packed dimers and their corresponding quantum RVB states are perfectly hyperuniform in general. Subsequently, we focus on a ruby-lattice spin liquid that was recently realized in a Rydberg-atom quantum simulator, and show that the QSL remains effectively hyperuniform even in the presence of a finite density of spinon and vison excitations, as long as the dimer constraint is still largely preserved. Moreover, we demonstrate that metrics based on the framework of hyperuniformity can be used to distinguish the QSL from other proximate quantum phases. These metrics can help identify potential QSL candidates, which can then be further analyzed using more advanced, computationally intensive quantum numerics to confirm their status as true QSLs.

A multi-layer multi-configurational time-dependent Hartree approach to lattice models beyond one dimension

The multi-layer multi-configurational time-dependent Hartree (MCTDH) approach is an efficient method to study quantum dynamics in real and imaginary time. The present work explores its potential to describe quantum fluids. The multi-layer MCTDH approach in second quantization representation is used to study lattice models beyond one dimension at finite temperatures. A scheme to map the lattice sites onto the MCTDH tree representation for multi-dimensional lattice models is proposed. A statistical sampling scheme previously used in MCTDH calculations is adapted to facilitate an efficient description of the thermal ensemble. As example, a two-dimensional hard-core Bose-Hubbard model is studied considering up to 64 × 64 lattice sites. The single particle function basis set size required to obtain converged results is found to not increase with the lattice size. The numerical results properly simulate the finite temperature Berezinskii-Kosterlitz-Thouless phase transition.

Tensor renormalization group for fermions

We review the basic ideas of the tensor renormalization group method and show how they can be applied for lattice field theory models involving relativistic fermions and Grassmann variables in arbitrary dimensions. We discuss recent progress for entanglement filtering, loop optimization, bond-weighting techniques and matrix product decompositions for Grassmann tensor networks. The new methods are tested with two-dimensional Wilson-Majorana fermions and multi-flavor Gross-Neveu models. We show that the methods can also be applied to the fermionic Hubbard model in 1+1 and 2+1 dimensions.

Entanglement in the quantum phases of an unfrustrated Rydberg atom array

Recent experimental advances have stimulated interest in the use of large, two-dimensional arrays of Rydberg atoms as a platform for quantum information processing and to study exotic many-body quantum states. However, the native long-range interactions between the atoms complicate experimental analysis and precise theoretical understanding of these systems. Here we use new tensor network algorithms capable of including all long-range interactions to study the ground state phase diagram of Rydberg atoms in a geometrically unfrustrated square lattice array. We find a greatly altered phase diagram from earlier numerical and experimental studies, revealed by studying the phases on the bulk lattice and their analogs in experiment-sized finite arrays. We further describe a previously unknown region with a nematic phase stabilized by short-range entanglement and an order from disorder mechanism. Broadly our results yield a conceptual guide for future experiments, while our techniques provide a blueprint for converging numerical studies in other lattices.

Efficient simulations of epidemic models with tensor networks: Application to the one-dimensional susceptible-infected-susceptible model

The contact process is an emblematic model of a nonequilibrium system, containing a phase transition between inactive and active dynamical regimes. In the epidemiological context, the model is known as the susceptible-infected-susceptible model, and it is widely used to describe contagious spreading. In this work, we demonstrate how accurate and efficient representations of the full probability distribution over all configurations of the contact process on a one-dimensional chain can be obtained by means of matrix product states (MPSs). We modify and adapt MPS methods from many-body quantum systems to study the classical distributions of the driven contact process at late times. We give accurate and efficient results for the distribution of large gaps, and illustrate the advantage of our methods over Monte Carlo simulations. Furthermore, we study the large deviation statistics of the dynamical activity, defined as the total number of configuration changes along a trajectory, and investigate quantum-inspired entropic measures, based on the second Rényi entropy.

Highly Resolved Spectral Functions of Two-Dimensional Systems with Neural Quantum States

Spectral functions are central to link experimental probes to theoretical models in condensed matter physics. However, performing exact numerical calculations for interacting quantum matter has remained a key challenge especially beyond one spatial dimension. In this work, we develop a versatile approach using neural quantum states to obtain spectral properties based on simulations of the dynamics of excitations initially localized in real or momentum space. We apply this approach to compute the dynamical structure factor in the vicinity of quantum critical points (QCPs) of different two-dimensional quantum Ising models, including one that describes the complex density wave orders of Rydberg atom arrays. When combined with deep network architectures we find that our method reliably describes dynamical structure factors of arrays with up to 24×24 spins, including the diverging timescales at critical points. Our approach is broadly applicable to interacting quantum lattice models in two dimensions and consequently opens up a route to compute spectral properties of correlated quantum matter in yet inaccessible regimes.

Emergent Glassy Behavior in a Kagome Rydberg Atom Array

We present large-scale quantum Monte Carlo simulation results on a realistic Hamiltonian of kagome-lattice Rydberg atom arrays. Although the system has no intrinsic disorder, intriguingly, our analyses of static and dynamic properties on large system sizes reveal emergent glassy behavior in a region of parameter space located between two valence bond solid phases. The extent of this glassy region is demarcated using the Edwards-Anderson order parameter, and its phase transitions to the two proximate valence bond solids-as well as the crossover towards a trivial paramagnetic phase-are identified. We demonstrate the intrinsically slow (imaginary) time dynamics deep inside the glassy phase and discuss experimental considerations for detecting such a quantum disordered phase with numerous nearly degenerate local minima. Our proposal paves a new route to the study of real-time glassy phenomena and highlights the potential for quantum simulation of a distinct phase of quantum matter beyond solids and liquids in current-generation Rydberg platforms.

Quantum neuronal sensing of quantum many-body states on a 61-qubit programmable superconducting processor

Classifying many-body quantum states with distinct properties and phases of matter is one of the most fundamental tasks in quantum many-body physics. However, due to the exponential complexity that emerges from the enormous numbers of interacting particles, classifying large-scale quantum states has been extremely challenging for classical approaches. Here, we propose a new approach called quantum neuronal sensing. Utilizing a 61-qubit superconducting quantum processor, we show that our scheme can efficiently classify two different types of many-body phenomena: namely the ergodic and localized phases of matter. Our quantum neuronal sensing process allows us to extract the necessary information coming from the statistical characteristics of the eigenspectrum to distinguish these phases of matter by measuring only one qubit and offers better phase resolution than conventional methods, such as measuring the imbalance. Our work demonstrates the feasibility and scalability of quantum neuronal sensing for near-term quantum processors and opens new avenues for exploring quantum many-body phenomena in larger-scale systems.

Emergent Z_{2} Gauge Theories and Topological Excitations in Rydberg Atom Arrays

Strongly interacting arrays of Rydberg atoms provide versatile platforms for exploring exotic many-body phases and dynamics of correlated quantum systems. Motivated by recent experimental advances, we show that the combination of Rydberg interactions and appropriate lattice geometries naturally leads to emergent Z_{2} gauge theories endowed with matter fields. Based on this mapping, we describe how Rydberg platforms could realize two distinct classes of topological Z_{2} quantum spin liquids, which differ in their patterns of translational symmetry fractionalization. We also discuss the natures of the fractionalized excitations of these Z_{2} spin liquid states using both fermionic and bosonic parton theories and illustrate their rich interplay with proximate solid phases.

Computing time-periodic steady-state currents via the time evolution of tensor network states

We present an approach based upon binary tree tensor network (BTTN) states for computing steady-state current statistics for a many-particle 1D ratchet subject to volume exclusion interactions. The ratcheted particles, which move on a lattice with periodic boundary conditions subject to a time-periodic drive, can be stochastically evolved in time to sample representative trajectories via a Gillespie method. In lieu of generating realizations of trajectories, a BTTN state can variationally approximate a distribution over the vast number of many-body configurations. We apply the density matrix renormalization group (DMRG) algorithm to initialize BTTN states, which are then propagated in time via the time-dependent variational principle (TDVP) algorithm to yield the steady-state behavior, including the effects of both typical and rare trajectories. The application of the methods to ratchet currents is highlighted, but the approach extends naturally to other interacting lattice models with time-dependent driving. Though trajectory sampling is conceptually and computationally simpler, we discuss situations for which the BTTN TDVP strategy can be beneficial.

Loop-free tensor networks for high-energy physics

This brief review introduces the reader to tensor network methods, a powerful theoretical and numerical paradigm spawning from condensed matter physics and quantum information science and increasingly exploited in different fields of research, from artificial intelligence to quantum chemistry. Here, we specialize our presentation on the application of loop-free tensor network methods to the study of high-energy physics problems and, in particular, to the study of lattice gauge theories where tensor networks can be applied in regimes where Monte Carlo methods are hindered by the sign problem. This article is part of the theme issue 'Quantum technologies in particle physics'.

Quantum phases of matter on a 256-atom programmable quantum simulator

Motivated by far-reaching applications ranging from quantum simulations of complex processes in physics and chemistry to quantum information processing¹, a broad effort is currently underway to build large-scale programmable quantum systems. Such systems provide insights into strongly correlated quantum matter^2-6, while at the same time enabling new methods for computation^7-10 and metrology¹¹. Here we demonstrate a programmable quantum simulator based on deterministically prepared two-dimensional arrays of neutral atoms, featuring strong interactions controlled by coherent atomic excitation into Rydberg states¹². Using this approach, we realize a quantum spin model with tunable interactions for system sizes ranging from 64 to 256 qubits. We benchmark the system by characterizing high-fidelity antiferromagnetically ordered states and demonstrating quantum critical dynamics consistent with an Ising quantum phase transition in (2 + 1) dimensions¹³. We then create and study several new quantum phases that arise from the interplay between interactions and coherent laser excitation¹⁴, experimentally map the phase diagram and investigate the role of quantum fluctuations. Offering a new lens into the study of complex quantum matter, these observations pave the way for investigations of exotic quantum phases, non-equilibrium entanglement dynamics and hardware-efficient realization of quantum algorithms.

Optimizing radiotherapy plans for cancer treatment with Tensor Networks

We present a novel application of Tensor Network methods in cancer treatment as a potential tool to solve the dose optimization problem in radiotherapy. In particular, the intensity-modulated radiation therapy technique-that allows treating irregular and inhomogeneous tumors while reducing the radiation toxicity on healthy organs-is based on the optimization problem of the beamlets intensities that shall result in a maximal delivery of the therapy dose to cancer while avoiding the organs at risk of being damaged by the radiation. The resulting optimization problem is expressed as a cost function to be optimized. Here, we map the cost function into an Ising-like Hamiltonian, describing a system of long-range interacting qubits. Finally, we solve the dose optimization problem by finding the ground-state of the Hamiltonian using a Tree Tensor Network algorithm. In particular, we present an anatomical scenario exemplifying a prostate cancer treatment. A similar approach can be applied to future hybrid classical-quantum algorithms, paving the way for the use of quantum technologies in future medical treatments.

Lattice quantum electrodynamics in (3+1)-dimensions at finite density with tensor networks

Gauge theories are of paramount importance in our understanding of fundamental constituents of matter and their interactions. However, the complete characterization of their phase diagrams and the full understanding of non-perturbative effects are still debated, especially at finite charge density, mostly due to the sign-problem affecting Monte Carlo numerical simulations. Here, we report the Tensor Network simulation of a three dimensional lattice gauge theory in the Hamiltonian formulation including dynamical matter: Using this sign-problem-free method, we simulate the ground states of a compact Quantum Electrodynamics at zero and finite charge densities, and address fundamental questions such as the characterization of collective phases of the model, the presence of a confining phase at large gauge coupling, and the study of charge-screening effects.