Properties Of Quantum Many-body Systems Research Articles

Exponentially improved efficient machine learning for quantum many-body states with provable guarantees

Solving the ground state and the ground-state properties of quantum many-body systems is generically a hard task for classical algorithms. For a family of Hamiltonians defined on an m-dimensional space of physical parameters, the ground state and its properties at an arbitrary parameter configuration can be predicted via a machine learning protocol up to a prescribed prediction error ɛ, provided that a sample set (of size N) of the states can be efficiently prepared and measured. In a recent work [Huang , ], a rigorous guarantee for such a generalization was proved. Unfortunately, an exponential scaling for the provable sample complexity, N=mO(1ɛ), was found to be universal for generic gapped Hamiltonians. This result applies to the situation where the dimension of the parameter space is large, while the scaling with the accuracy is not an urgent factor, not entering the realm of more precise learning and prediction. In this work, we consider an alternative relevant scenario, where the effective dimension m is a finite, not necessarily large constant, while the scaling with the prediction error becomes the central concern. By jointly preserving the fundamental properties of density matrices in the learning protocol and utilizing the continuity of quantum states in the parameter range of interest, we rigorously obtain a polynomial sample complexity for predicting quantum many-body states and their properties, with respect to the prediction error ɛ and the number of qubits, n, with N=poly(ɛ−1,n,ln1δ), where poly denotes a polynomial function, and (1−δ) is the probability of success. Moreover, if restricted to learning local quantum-state properties, the number of samples can be further reduced to N=poly(ɛ−1,lnnδ). Numerical demonstrations confirm our findings, and an alternative approach utilizing statistical learning theory with reproducing kernel Hilbert space achieves consistent results. The mere continuity assumption indicates that our results are not restricted to gapped Hamiltonian systems and properties within the same phase. Published by the American Physical Society 2024

Deconfinement Dynamics of Fractons in Tilted Bose-Hubbard Chains.

Fractonic constraints can lead to exotic properties of quantum many-body systems. Here, we investigate the dynamics of fracton excitations on top of the ground states of a one-dimensional, dipole-conserving Bose-Hubbard model. We show that nearby fractons undergo a collective motion mediated by exchanging virtual dipole excitations, which provides a powerful dynamical tool to characterize the underlying ground-state phases. We find that, in the gapped Mott insulating phase, fractons are confined to each other as motion requires the exchange of massive dipoles. When crossing the phase transition into a gapless Luttinger liquid of dipoles, fractons deconfine. Their transient deconfinement dynamics scales diffusively and exhibits strong but subleading contributions described by a quantum Lifshitz model. We examine prospects for the experimental realization in tilted Bose-Hubbard chains by numerically simulating the adiabatic state preparation and subsequent time evolution and find clear signatures of the low-energy fracton dynamics.

Robust Extraction of Thermal Observables from State Sampling and Real-Time Dynamics on Quantum Computers

Simulating properties of quantum materials is one of the most promising applications of quantum computation, both near- and long-term. While real-time dynamics can be straightforwardly implemented, the finite temperature ensemble involves non-unitary operators that render an implementation on a near-term quantum computer extremely challenging. Recently, Lu, Bañuls and Cirac \cite{Lu2021} suggested a "time-series quantum Monte Carlo method" which circumvents this problem by extracting finite temperature properties from real-time simulations via Wick's rotation and Monte Carlo sampling of easily preparable states. In this paper, we address the challenges associated with the practical applications of this method, using the two-dimensional transverse field Ising model as a testbed. We demonstrate that estimating Boltzmann weights via Wick's rotation is very sensitive to time-domain truncation and statistical shot noise. To alleviate this problem, we introduce a technique that imposes constraints on the density of states, most notably its non-negativity, and show that this way, we can reliably extract Boltzmann weights from noisy time series. In addition, we show how to reduce the statistical errors of Monte Carlo sampling via a reweighted version of the Wolff cluster algorithm. Our work enables the implementation of the time-series algorithm on present-day quantum computers to study finite temperature properties of many-body quantum systems.

Nearly-frustration-free ground state preparation

Solving for quantum ground states is important for understanding the properties of quantum many-body systems, and quantum computers are potentially well-suited for solving for quantum ground states. Recent work [1] has presented a nearly optimal scheme that prepares ground states on a quantum computer for completely generic Hamiltonians, whose query complexity scales as δ−1, i.e. inversely with their normalized gap. Here we consider instead the ground state preparation problem restricted to a special subset of Hamiltonians, which includes those which we term "nearly-frustration-free": the class of Hamiltonians for which the ground state energy of their block-encoded and hence normalized Hamiltonian α−1H is within δy of -1, where δ is the spectral gap of α−1H and 0≤y≤1. For this subclass, we describe an algorithm whose dependence on the gap is asymptotically better, scaling as δy/2−1, and show that this new dependence is optimal up to factors of log⁡δ. In addition, we give examples of physically motivated Hamiltonians which live in this subclass. Finally, we describe an extension of this method which allows the preparation of excited states both for generic Hamiltonians as well as, at a similar speedup as the ground state case, for those which are nearly frustration-free.

A fast time domain solver for the equilibrium Dyson equation

We consider the numerical solution of the real-time equilibrium Dyson equation, which is used in calculations of the dynamical properties of quantum many-body systems. We show that this equation can be written as a system of coupled, nonlinear, convolutional Volterra integro-differential equations, for which the kernel depends self-consistently on the solution. As is typical in the numerical solution of Volterra-type equations, the computational bottleneck is the quadratic-scaling cost of history integration. However, the structure of the nonlinear Volterra integral operator precludes the use of standard fast algorithms. We propose a quasilinear-scaling FFT-based algorithm which respects the structure of the nonlinear integral operator. The resulting method can reach large propagation times and is thus well-suited to explore quantum many-body phenomena at low energy scales. We demonstrate the solver with two standard model systems: the Bethe graph and the Sachdev-Ye-Kitaev model.

Dynamical properties of quantum many-body systems with long-range interactions

Dynamical properties of quantum many-body systems with long-range interactions

Scrambling and quantum chaos indicators from long-time properties of operator distributions

Scrambling is a key concept in the analysis of nonequilibrium properties of quantum many-body systems. Most studies focus on its characterization via out-of-time-ordered correlation functions (OTOCs), particularly through the early-time decay of the OTOC. However, scrambling is a complex process which involves operator spreading and operator entanglement, and a full characterization requires one to access more refined information on the operator dynamics at several timescales. In this work we analyze operator scrambling by expanding the target operator in a complete basis and studying the structure of the expansion coefficients treated as a coarse-grained probability distribution in the space of operators. We study different features of this distribution, such as its mean, variance, and participation ratio, for the Ising model with longitudinal and transverse fields, kicked collective spin models, and random circuit models. We show that the long-time properties of the operator distribution display common features across these cases, and discuss how these properties can be used as a proxy for the onset of quantum chaos. Finally, we discuss the connection with OTOCs and analyze the cost of probing the operator distribution experimentally using these correlation functions.

Overcoming the entanglement barrier in quantum many-body dynamics via space-time duality

Describing nonequilibrium properties of quantum many-body systems is challenging due to high entanglement in the wave function. We describe the evolution of local observables via the influence matrix (IM), which encodes the effects of a many-body system as an environment for local subsystems. Recent works found that in many dynamical regimes the IM of an infinite system has low temporal entanglement and can be efficiently represented as a matrix-product state (MPS). Yet, direct iterative constructions of the IM encounter highly entangled intermediate states---a temporal entanglement barrier (TEB). We argue that TEB is ubiquitous, and elucidate its physical origin via a semiclassical quasiparticle picture that exactly captures the behavior of integrable spin chains. Further, we show that a TEB also arises in chaotic spin chains, which lack well-defined quasiparticles. Based on these insights, we formulate an alternative light-cone growth algorithm, which provably avoids TEB, thus providing an efficient construction of the thermodynamic-limit IM as a MPS. This work uncovers the origin of the efficiency of the IM approach for studying thermalization and transport.

Error-resilient Monte Carlo quantum simulation of imaginary time

Computing the ground-state properties of quantum many-body systems is a promising application of near-term quantum hardware with a potential impact in many fields. The conventional algorithm quantum phase estimation uses deep circuits and requires fault-tolerant technologies. Many quantum simulation algorithms developed recently work in an inexact and variational manner to exploit shallow circuits. In this work, we combine quantum Monte Carlo with quantum computing and propose an algorithm for simulating the imaginary-time evolution and solving the ground-state problem. By sampling the real-time evolution operator with a random evolution time according to a modified Cauchy-Lorentz distribution, we can compute the expected value of an observable in imaginary-time evolution. Our algorithm approaches the exact solution given a circuit depth increasing polylogarithmically with the desired accuracy. Compared with quantum phase estimation, the Trotter step number, i.e. the circuit depth, can be thousands of times smaller to achieve the same accuracy in the ground-state energy. We verify the resilience to Trotterisation errors caused by the finite circuit depth in the numerical simulation of various models. The results show that Monte Carlo quantum simulation is promising even without a fully fault-tolerant quantum computer.

Deep Reinforcement Learning for Preparation of Thermal and Prethermal Quantum States

We propose a method based on deep reinforcement learning that efficiently prepares a quantum many-body pure state in thermal or prethermal equilibrium. The main physical intuition underlying the method is that the information on the equilibrium states can be efficiently encoded/extracted by focusing on only a few local observables, relying on the typicality of equilibrium states. Instead of resorting to the expensive preparation protocol that adopts global features such as the quantum state fidelity, we show that the equilibrium states can be efficiently prepared only by learning the expectation values of local observables. We demonstrate our method by preparing two illustrative examples: Gibbs ensembles in non-integrable systems and generalized Gibbs ensembles in integrable systems. Pure states prepared solely from local observables are numerically shown to successfully encode the macroscopic properties of the equilibrium states. Furthermore, we find that the preparation errors, with respect to the system size, decay exponentially for Gibbs ensembles and polynomially for generalized Gibbs ensembles, which are in agreement with the finite-size fluctuation within thermodynamic ensembles. Our method paves a path toward studying the thermodynamic and statistical properties of quantum many-body systems in quantum hardware.

Variational quantum simulation of thermal statistical states on a superconducting quantum processer

Quantum computers promise to solve finite-temperature properties of quantum many-body systems, which is generally challenging for classical computers due to high computational complexities. Here, we report experimental preparations of Gibbs states and excited states of Heisenberg XX and XXZ models by using a 5-qubit programmable superconducting processor. In the experiments, we apply a hybrid quantum–classical algorithm to generate finite temperature states with classical probability models and variational quantum circuits. We reveal that the Hamiltonians can be fully diagonalized with optimized quantum circuits, which enable us to prepare excited states at arbitrary energy density. We demonstrate that the approach has a self-verifying feature and can estimate fundamental thermal observables with a small statistical error. Based on numerical results, we further show that the time complexity of our approach scales polynomially in the number of qubits, revealing its potential in solving large-scale problems.

Life and death of the Bose polaron

Spectroscopic and interferometric measurements complement each other in extracting the fundamental properties of quantum many-body systems. While spectroscopy provides precise measurements of equilibrated energies, interferometry can elucidate the dynamical evolution of the system. For an impurity immersed in a bosonic medium, both are equally important for understanding the quasiparticle physics of the Bose polaron. Here, we compare the interferometric and spectroscopic timescales to the underlying dynamical regimes of the impurity dynamics and the polaron lifetime, highlighting the capability of the interferometric approach to clearly resolve polaron dynamics. In particular, interferometric measurements of the coherence amplitude at strong interactions reveal faster quantum dynamics at large repulsive interaction strengths than at unitarity. These observations are in excellent agreement with a short-time theoretical prediction including both the continuum and the attractive polaron branch. For longer times, qualitative agreement with a many-body theoretical prediction which includes both branches is obtained. Moreover, the polaron energy is extracted from interferometric measurements of the observed phase velocity in agreement with previous spectroscopic results from weak to strong attractive interactions. Finally, the phase evolution allows for the measurement of an energetic equilibration timescale, describing the initial approach of the phase velocity to the polaron energy. Theoretically, this is shown to lie within the regime of universal dynamics revealing a fast initial evolution towards the formation of polarons. Our results give a comprehensive picture of the many-body physics governing the Bose polaron and thus validates the quasiparticle framework for further studies.

Entanglement spread area law in gapped ground states

In this work, we make a connection between two seemingly different problems. The first problem involves characterizing the properties of entanglement in the ground state of gapped local Hamiltonians, which is a central topic in quantum many-body physics. The second problem is on the quantum communication complexity of testing bipartite states with EPR assistance, a well-known question in quantum information theory. We construct a communication protocol for testing (or measuring) the ground state and use its communication complexity to reveal a new structural property for the ground state entanglement. This property, known as the entanglement spread, roughly measures the ratio between the largest and the smallest Schmidt coefficients across a cut in the ground state. Our main result shows that gapped ground states possess limited entanglement spread across any cut, exhibiting an "area law" behavior. Our result quite generally applies to any interaction graph with an improved bound for the special case of lattices. This entanglement spread area law includes interaction graphs constructed in [Aharonov et al., FOCS'14] that violate a generalized area law for the entanglement entropy. Our construction also provides evidence for a conjecture in physics by Li and Haldane on the entanglement spectrum of lattice Hamiltonians [Li and Haldane, PRL'08]. On the technical side, we use recent advances in Hamiltonian simulation algorithms along with quantum phase estimation to give a new construction for an approximate ground space projector (AGSP) over arbitrary interaction graphs.

Probing precise interatomic potentials by nonadiabatic nonlinear phonons

Accurately probing the interatomic potential (IP) in crystals is essential for understanding the dynamic mechanisms in phase transitions and chemical reactions. Here, by interrogating laser-induced coherent phonon dynamics, we propose a new kind of optical microscopy for determining effective interatomic potentials beyond the harmonic and adiabatic approximation. This technique is tested against available experimental and first-principles data, by which the interatomic potential of graphene and a few other materials are successfully reconstructed in a large configuration space and beyond the ground state. A significant nonadiabatic effect emerges in nonlinear phonon dynamics and the IP reconstruction, which corrects the anharmonic part of IP much stronger (up to 50%) than the harmonic part. The nonadiabaticity-modulated potentials can lead to ∼10% correction on thermal expansion coefficients and ∼0.3 eV on phase transition barriers. This work offers a new strategy for probing complex interactions in quantum materials and exemplifies the universality of nonadiabaticity, which deserves a full consideration for precisely determining physical properties of quantum many-body systems.

Classical algorithms for many-body quantum systems at finite energies

Properties of quantum many-body systems at finite energy densities, such as microcanonical averages, are in general difficult to compute with classical algorithms. Quantum algorithms that combine quantum dynamics with filtering and sampling techniques can however access some of these quantities. Here, the authors show that simulating the time evolution part with tensor networks provides quantum-inspired classical algorithms that reach much larger systems than other numerical methods.

Entanglement of Skeletal Regions.

The entanglement entropy (EE) encodes key properties of quantum many-body systems. It is usually calculated for subregions of finite volume (or area in 2D). Here, we study the EE of skeletal regions that have no volume, such as a line in 2D. We show that skeletal entanglement displays new behavior compared with its bulk counterpart, and leads to distinct universal quantities. We provide nonperturbative bounds for the skeletal area-law coefficient of a large family of quantum states. We then explore skeletal scaling for the toric code, conformal bosons and Dirac fermions, Lifshitz critical points, and Fermi liquids. We discover signatures including skeletal topological EE, novel corner terms, and strict area-law scaling for metals. These findings suggest that skeletal entropy serves as a measure for the range of entanglement. Finally, we outline open questions relating to other systems and measures such as the logarithmic negativity.

Robust analytic continuation combining the advantages of the sparse modeling approach and the Padé approximation

Analytic continuation (AC) from the imaginary-time Green's function to the spectral function is a crucial process for numerical studies of the dynamical properties of quantum many-body systems. This process, however, is an ill-posed problem; that is, the obtained spectrum is unstable against the noise of the Green's function. Though several numerical methods have been developed, each of them has its own advantages and disadvantages. The sparse modeling (SpM) AC method, for example, is robust against the noise of the Green's function but suffers from unphysical oscillations in the low-energy region. We propose a new method that combines the SpM AC with the Pad\'{e} approximation. This combination, called SpM-Pad\'{e}, inherits robustness against noise from SpM and low-energy accuracy from Pad\'{e}, compensating for the disadvantages of each. We demonstrate that the SpM- Pad\'{e} method yields low-variance and low-biased results with almost the same computational cost as that of the SpM method.

Subdiffusion and Many-Body Quantum Chaos with Kinetic Constraints.

We investigate the spectral and transport properties of many-body quantum systems with conserved charges and kinetic constraints. Using random unitary circuits, we compute ensemble-averaged spectral form factors and linear-response correlation functions, and find that their characteristic timescales are given by the inverse gap of an effective Hamiltonian-or equivalently, a transfer matrix describing a classical Markov process. Our approach allows us to connect directly the Thouless time, t_{Th}, determined by the spectral form factor, to transport properties and linear-response correlators. Using tensor network methods, we determine the dynamical exponent z for a number of constrained, conserving models. We find universality classes with diffusive, subdiffusive, quasilocalized, and localized dynamics, depending on the severity of the constraints. In particular, we show that quantum systems with "Fredkin" constraints exhibit anomalous transport with dynamical exponent z≃8/3.

On the Accuracy of Random Phase Approximation for Dynamical Structure Factors in Cold Atomic Gases

Many-body physics poses one of the greatest challenges to science in the 21st century. Still more daunting is the problem of accurately calculating the properties of quantum many-body systems in the strongly correlated regime. Cold atomic gases provide an excellent test ground, for both experimentalists and theorists, to study the exotic and sometimes counterintuitive behavior of quantum many-body problems. Of particular interest is the appearance of collective excitations in these systems, such as the famous Goldstone mode and the elusive Higgs mode. It is particularly important to assess the robustness of theoretical and computational techniques to study such excitations. We build on the unprecedented opportunity provided by the fact that, in some cases, exact numerical predictions can be obtained through quantum Monte Carlo. We use these predictions to assess the accuracy of the Random Phase Approximation, which is widely considered to be a method of choice for the study of the collective excitations in a cold atomic Fermi gas modeled with a Fermi–Hubbard Hamiltonian. We found good agreement between the two methodologies for the dynamic properties, particularly for the position of the Goldstone mode. We also explored the possibility of using a renormalized, effective potential in place of the physical potential. We determined that using a renormalized potential is likely too simplistic a method for improving the accuracy of generalized Random Phase Approximation and that a more sophisticated approach is needed.

Local quench spectroscopy of many-body quantum systems

Quench spectroscopy is a relatively new method which enables the investigation of spectral properties of many-body quantum systems by monitoring the out-of-equilibrium dynamics of real-space observables after a quench. So far the approach has been devised for global quenches or using local engineering of momentum-resolved excitations. Here, we extend the quench spectroscopy method to local quenches. We show that it allows us to extract quantitative information about global properties of the system, and in particular the elementary excitation spectrum. Using state-of-the-art numerical methods, we simulate the out-of-equilibrium dynamics of a variety of quantum systems following various local quench protocols and demonstrate a general scheme for designing an appropriate local quench protocol for any chosen model. We provide detailed examples of how the local quench protocol can be realised in realistic current generation experiments, including ultracold atomic gases and trapped ion systems.