Quantum Many-body Systems Research Articles

Quantum Fluctuation Dynamics of Dispersive Superradiant Pulses in a Hybrid Light-Matter System.

We consider theoretically a driven-dissipative quantum many-body system consisting of an atomic ensemble in a single-mode optical cavity as described by the open Tavis-Cummings model. In this hybrid light-matter system, the interplay between coherent and dissipative processes leads to superradiant pulses with a buildup of strong correlations, even for systems comprising hundreds to thousands of particles. A central feature of the mean-field dynamics is a self-reversal of two spin degrees of freedom due to an underlying time-reversal symmetry, which is broken by quantum fluctuations. We demonstrate a quench protocol that can maintain highly non-Gaussian states over long timescales. This general mechanism offers interesting possibilities for the generation and control of complex fluctuation patterns, as suggested for the improvement of quantum sensing protocols for dissipative spin amplification.

Ballistic macroscopic fluctuation theory

We introduce a new universal framework describing fluctuations and correlations in quantum and classical many-body systems, at the Euler hydrodynamic scale of space and time. The framework adapts the ideas of the conventional macroscopic fluctuation theory (MFT) to systems that support ballistic transport. The resulting “ballistic MFT” (BMFT) is solely based on the Euler hydrodynamics data of the many-body system. Within this framework, mesoscopic observables are classical random variables depending only on the fluctuating conserved densities, and Euler-scale fluctuations are obtained by deterministically transporting thermodynamic fluctuations via the Euler hydrodynamics. Using the BMFT, we show that long-range correlations in space generically develop over time from long-wavelength inhomogeneous initial states in interacting models. This result, which we verify by numerical calculations, challenges the long-held paradigm that at the Euler scale, fluid cells may be considered uncorrelated. We also show that the Gallavotti-Cohen fluctuation theorem for non-equilibrium ballistic transport follows purely from time-reversal invariance of the Euler hydrodynamics. We check the validity of the BMFT by applying it to integrable systems, and in particular the hard-rod gas, with extensive simulations that confirm our analytical results.

Extensive long-range entanglement in a nonequilibrium steady state

Entanglement measures constitute powerful tools in the quantitative description of quantum many-body systems out of equilibrium. We study entanglement in the current-carrying steady state of a paradigmatic one-dimensional model of noninteracting fermions at zero temperature in the presence of a scatterer. We show that disjoint intervals located on opposite sides of the scatterer, and within similar distances from it, maintain volume-law entanglement regardless of their separation, as measured by their fermionic negativity and coherent information. The mutual information of the intervals, which quantifies the total correlations between them, follows a similar scaling. Interestingly, this scaling entails in particular that if the position of one of the intervals is kept fixed, then the correlation measures depend non-monotonically on the distance between the intervals. By deriving exact expressions for the extensive terms of these quantities, we prove their simple functional dependence on the scattering probabilities, and demonstrate that the strong long-range entanglement is generated by the coherence between the transmitted and reflected parts of propagating particles within the bias-voltage window. The generality and simplicity of the model suggest that this behavior should characterize a large class of nonequilibrium steady states.

Complexity of frustration: A new source of non-local non-stabilizerness

We advance the characterization of complexity in quantum many-body systems by examining WW-states embedded in a spin chain. Such states show an amount of non-stabilizerness or “magic”, measured as the Stabilizer Rényi Entropy, that grows logarithmically with the number of qubits/spins. We focus on systems whose Hamiltonian admits a classical point with extensive degeneracy. Near these points, a Clifford circuit can convert the ground state into a WW-state, while in the rest of the phase to which the classical point belongs, it is dressed with local quantum correlations. Topological frustrated quantum spin-chains host phases with the desired phenomenology, and we show that their ground state’s Stabilizer Rényi Entropy is the sum of that of the WW-states plus an extensive local contribution. Our work reveals that WW-states/frustrated ground states display a non-local degree of complexity that can be harvested as a quantum resource and has no counterpart in GHZ states/non-frustrated systems.

Residual Tensor Train: A Quantum-Inspired Approach for Learning Multiple Multilinear Correlations

States of quantum many-body systems are defined in a high-dimensional Hilbert space, where rich and complex interactions among subsystems can be modelled. In machine learning, complex multiple multilinear correlations may also exist within input features. In this paper, we present a quantum-inspired multilinear model, named Residual Tensor Train (ResTT), to capture the multiple multilinear correlations of features, from low to high orders, within a single model. ResTT is able to build a robust decision boundary in a high-dimensional space for solving fitting and classification tasks. In particular, we prove that the fully-connected layer and the Volterra series can be taken as special cases of ResTT. Furthermore, we derive the rule for weight initialization that stabilizes the training of ResTT based on a mean-field analysis. We prove that such a rule is much more relaxed than that of TT, which means ResTT can easily address the vanishing and exploding gradient problem that exists in the existing TT models. Numerical experiments demonstrate that ResTT outperforms the state-of-the-art tensor network and benchmark deep learning models on MNIST and Fashion-MNIST datasets. Moreover, ResTT achieves better performance than other statistical methods on two practical examples with limited data which are known to have complex feature interactions.

Quantum simulation of bosons with the contracted quantum eigensolver

Quantum computers are promising tools for simulating many-body quantum systems due to their potential scaling advantage over classical computers. While significant effort has been expended on many-fermion systems, here we simulate a model entangled many-boson system with the contracted quantum eigensolver (CQE). We generalize the CQE to many-boson systems by encoding the bosonic wavefunction on qubits. The CQE provides a compact ansatz for the bosonic wave function whose gradient is proportional to the residual of a contracted Schrödinger equation. We apply the CQE to a bosonic system, where N quantum harmonic oscillators are coupled through a pairwise quadratic repulsion. The model is relevant to the study of coupled vibrations in molecular systems on quantum devices. Results demonstrate the potential efficiency of the CQE in simulating bosonic processes such as molecular vibrations with good accuracy and convergence even in the presence of noise.

Late-time universal distribution functions of observables in one-dimensional many-body quantum systems.

We study the probability distribution function of the long-time values of observables being time-evolved by Hamiltonians modeling clean and disordered one-dimensional chains of many spin-1/2 particles. In particular, we analyze the return probability and its version for a completely extended initial state, the so-called spectral form factor. We complement our analysis with the spin autocorrelation and connected spin-spin correlation functions, both of interest in experiments with quantum simulators. We show that the distribution function has a universal shape provided the central limit theorem holds. Explicitly, the shape is exponential for the return probability and spectral form factor, meanwhile it is Gaussian for the few-body observables. We also discuss implications over the so-called many-body localization. Remarkably, our approach requires only a single sample of the dynamics and small system sizes, which could be quite advantageous when dealing specially with disordered systems.

Many-body quantum sign structures as non-glassy Ising models

The non-trivial phase structure of the eigenstates of many-body quantum systems severely limits the applicability of quantum Monte Carlo, variational, and machine learning methods. Here, we study real-valued signful ground-state wave functions of frustrated quantum spin systems and, assuming that the tasks of finding wave function amplitudes and signs can be separated, show that the signs can be easily bootstrapped from the amplitudes. We map the problem of finding the sign structure to an auxiliary classical Ising model defined on a subset of the Hilbert space basis. We show that the Ising model does not exhibit significant frustrations even for highly frustrated parental quantum systems, and is solvable with a fully deterministic O(Klog K)-time combinatorial algorithm (where K is the Ising model size). Given the ground state amplitudes, we reconstruct the signs of the ground states of several frustrated quantum models, thereby revealing the hidden simplicity of many-body sign structures.

String-net formulation of Hamiltonian lattice Yang-Mills theories and quantum many-body scars in a nonabelian gauge theory

We study the Hamiltonian lattice Yang-Mills theory based on spin networks that provide a useful basis to represent the physical states satisfying the Gauss law constraints. We focus on SU(2) Yang-Mills theory in (2 + 1) dimensions. Following the string-net model, we introduce a regularization of the Kogut-Susskind Hamiltonian of lattice Yang-Mills theory based on the q deformation, which respects the (discretized) SU(2) gauge symmetry as quantum group, i.e., SU(2)k, and enables implementation of the lattice Yang-Mills theory both in classical and quantum algorithms by referring to those of the string-net model. Using the regularized Hamiltonian, we study quantum scars in a nonabelian gauge theory. Quantum scars are nonthermal energy eigenstates arising in the constrained quantum many-body systems. We find that quantum scars from zero modes, which have been found in abelian gauge theories arise even in a nonabelian gauge theory. We also show the spectrum of a single-plaquette model for SU(2)k and SU(3)k with naive cutoff and that based on the q-deformation to discuss cutoff dependence of the formulation.

Unveiling order from chaos by approximate 2-localization of random matrices

Quantum many-body systems are typically endowed with a tensor product structure. A structure they inherited from probability theory, where the probability of two independent events is the product of the probabilities. The tensor product structure of a Hamiltonian thus gives a natural decomposition of the system into independent smaller subsystems. It is interesting to understand whether a given Hamiltonian is compatible with some particular tensor product structure. In particular, we ask, is there a basis in which an arbitrary Hamiltonian has a 2-local form, i.e., it contains only pairwise interactions? Here we show, using analytical and numerical calculations, that a generic Hamiltonian (e.g., a large random matrix) can be approximately written as a linear combination of two-body interaction terms with high precision; that is, the Hamiltonian is 2-local in a carefully chosen basis. Moreover, we show that these Hamiltonians are not fine-tuned, meaning that the spectrum is robust against perturbations of the coupling constants. Finally, by analyzing the adjacency structure of the couplings [Formula: see text], we suggest a possible mechanism for the emergence of geometric locality from quantum chaos.

Adaptive variational quantum minimally entangled typical thermal states for finite temperature simulations

Scalable quantum algorithms for the simulation of quantum many-body systems in thermal equilibrium are important for predicting properties of quantum matter at finite temperatures. Here we describe and benchmark a quantum computing version of the minimally entangled typical thermal states (METTS) algorithm for which we adopt an adaptive variational approach to perform the required quantum imaginary time evolution. The algorithm, which we name AVQMETTS, dynamically generates compact and problem-specific quantum circuits, which are suitable for noisy intermediate-scale quantum (NISQ) hardware. We benchmark AVQMETTS on statevector simulators and perform thermal energy calculations of integrable and nonintegrable quantum spin models in one and two dimensions and demonstrate an approximately linear system-size scaling of the circuit complexity. We further map out the finite-temperature phase transition line of the two-dimensional transverse field Ising model. Finally, we study the impact of noise on AVQMETTS calculations using a phenomenological noise model.

Full counting statistics as probe of measurement-induced transitions in the quantum Ising chain

Non-equilibrium dynamics of many-body quantum systems under the effect of measurement protocols is attracting an increasing amount of attention. It has been recently revealed that measurements may induce different non-equilibrium regimes and an abrupt change in the scaling-law of the bipartite entanglement entropy. However, our understanding of how these regimes appear{, how they affect the statistics of local quantities and}, finally whether they survive in the thermodynamic limit, is much less established. Here we investigate measurement-induced phase transitions in the Quantum Ising chain coupled to a monitoring environment. In particular we show that local projective measurements induce a quantitative modification of the out-of-equilibrium probability distribution function of the local magnetization. Starting from a GHZ state, the relaxation of the paramagnetic and the ferromagnetic order is analysed. In particular we describe how the probability distributions associated to them show different behaviour depending on the measurement rate.

Quantum information geometry of driven CFTs

Driven quantum systems exhibit a large variety of interesting and sometimes exotic phenomena. Of particular interest are driven conformal field theories (CFTs) which describe quantum many-body systems at criticality. In this paper, we develop both a spacetime and a quantum information geometry perspective on driven 2d CFTs. We show that for a large class of driving protocols the theories admit an alternative but equivalent formulation in terms of a CFT defined on a spacetime with a time-dependent metric. We prove this equivalence both in the operator formulation as well as in the path integral description of the theory. A complementary quantum information geometric perspective for driven 2d CFTs employs the so-called Bogoliubov-Kubo-Mori (BKM) metric, which is the counterpart of the Fisher metric of classical information theory, and which is obtained from a perturbative expansion of relative entropy. We compute the BKM metric for the universal sector of Virasoro excitations of a thermal state, which captures a large class of driving protocols, and find it to be a useful tool to classify and characterize different types of driving. For Möbius driving by the SL(2, ℝ) subgroup, the BKM metric becomes the hyperbolic metric on the disk. We show how the non-trivial dynamics of Floquet driven CFTs is encoded in the BKM geometry via Möbius transformations. This allows us to identify ergodic and non-ergodic regimes in the driving. We also explain how holographic driven CFTs are dual to driven BTZ black holes with evolving horizons. The deformation of the black hole horizon towards and away from the asymptotic boundary provides a holographic understanding of heating and cooling in Floquet CFTs.

Exact Dirac-Bogoliubov-de Gennes Dynamics for Inhomogeneous Quantum Liquids.

We study inhomogeneous 1+1-dimensional quantum many-body systems described by Tomonaga-Luttinger-liquid theory with general propagation velocity and Luttinger parameter varying smoothly in space, equivalent to an inhomogeneous compactification radius for free boson conformal field theory. This model appears prominently in low-energy descriptions, including for trapped ultracold atoms, while here we present an application to quantum Hall edges with inhomogeneous interactions. The dynamics is shown to be governed by a pair of coupled continuity equations identical to inhomogeneous Dirac-Bogoliubov-de Gennes equations with a local gap and solved by analytical means. We obtain their exact Green's functions and scattering matrix using a Magnus expansion, which generalize previous results for conformal interfaces and quantum wires coupled to leads. Our results explicitly describe the late-time evolution following quantum quenches, including inhomogeneous interaction quenches, and Andreev reflections between coupled quantum Hall edges, revealing remarkably universal dependence on details at stationarity or at late times out of equilibrium.

Learning ground states of gapped quantum Hamiltonians with Kernel Methods

Neural network approaches to approximate the ground state of quantum hamiltonians require the numerical solution of a highly nonlinear optimization problem. We introduce a statistical learning approach that makes the optimization trivial by using kernel methods. Our scheme is an approximate realization of the power method, where supervised learning is used to learn the next step of the power iteration. We show that the ground state properties of arbitrary gapped quantum hamiltonians can be reached with polynomial resources under the assumption that the supervised learning is efficient. Using kernel ridge regression, we provide numerical evidence that the learning assumption is verified by applying our scheme to find the ground states of several prototypical interacting many-body quantum systems, both in one and two dimensions, showing the flexibility of our approach.

Linear growth of circuit complexity from Brownian dynamics

How rapidly can a many-body quantum system generate randomness? Using path integral methods, we demonstrate that Brownian quantum systems have circuit complexity that grows linearly with time. In particular, we study Brownian clusters of N spins or fermions with time-dependent all-to-all interactions, and calculate the Frame Potential to characterize complexity growth in these models. In both cases the problem can be mapped to an effective statistical mechanics problem which we study using path integral methods. Within this framework it is straightforward to show that the kth Frame Potential comes within ϵ of the Haar value after a time of order t ~ kN + k log k + log ϵ−1. Using a bound on the diamond norm, this implies that such circuits are capable of coming very close to a unitary k-design after a time of order t ~ kN. We also consider the same question for systems with a time-independent Hamiltonian and argue that a small amount of time-dependent randomness is sufficient to generate a k-design in linear time provided the underlying Hamiltonian is quantum chaotic. These models provide explicit examples of linear complexity growth that are analytically tractable and are directly applicable to practical applications calling for unitary k-designs.

Quantifying hole-motion-induced frustration in doped antiferromagnets by Hamiltonian reconstruction

Unveiling the microscopic origins of quantum phases dominated by the interplay of spin and motional degrees of freedom constitutes one of the central challenges in strongly correlated many-body physics. When holes move through an antiferromagnetic spin background, they displace the positions of spins, which induces effective frustration in the magnetic environment. However, a concrete characterization of this effect in a quantum many-body system is still an unsolved problem. Here we present a Hamiltonian reconstruction scheme that allows for a precise quantification of hole-motion-induced frustration. We access non-local correlation functions through projective measurements of the many-body state, from which effective spin-Hamiltonians can be recovered after detaching the magnetic background from dominant charge fluctuations. The scheme is applied to systems of mixed dimensionality, where holes are restricted to move in one dimension, but SU(2) superexchange is two-dimensional. We demonstrate that hole motion drives the spin background into a highly frustrated regime, which can quantitatively be described by an effective J1–J2-type spin model. We exemplify the applicability of the reconstruction scheme to ultracold atom experiments by recovering effective spin-Hamiltonians of experimentally obtained 1D Fermi-Hubbard snapshots. Our method can be generalized to fully 2D systems, enabling promising microscopic perspectives on the doped Hubbard model.

Formation of Tesseract Time Crystals on a Quantum Computer

The engineering of new states of matter through Floquet driving has revolutionized the field of condensed matter physics. This technique enables the creation of hybrid topological states and ordered phases that are absent in normal systems. Crystalline structures, exemplifying spatially ordered systems under periodic driving, have been extensively studied. However, recent focus has shifted towards discrete time crystals (DTCs), periodically driven quantum many-body systems that break time translation symmetry under specific conditions. In this paper, the model of discrete time crystals is extended to allow for the formation of time-varying tesseracts, allowing for the investigation of time translational symmetry in pseudo-higher-dimensional lattice 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.

Physical insights from imaginary-time density–density correlation functions

An accurate theoretical description of the dynamic properties of correlated quantum many-body systems, such as the dynamic structure factor S(q, ω), is important in many fields. Unfortunately, highly accurate quantum Monte Carlo methods are usually restricted to the imaginary time domain, and the analytic continuation of the imaginary-time density–density correlation function F(q, τ) to real frequencies is a notoriously hard problem. Here, it is argued that often no such analytic continuation is required because by definition, F(q, τ) contains the same physical information as does S(q, ω), only represented unfamiliarly. Specifically, it is shown how one can directly extract key information such as the temperature or quasi-particle excitation energies from the τ domain, which is highly relevant for equation-of-state measurements of matter under extreme conditions [T. Dornheim et al., Nat. Commun. 13, 7911 (2022)]. As a practical example, ab initio path-integral Monte Carlo results for the uniform electron gas (UEG) are considered, and it is shown that even nontrivial processes such as the roton feature of the UEG at low density [T. Dornheim et al., Commun. Phys. 5, 304 (2022)] are manifested straightforwardly in F(q, τ). A comprehensive overview is given of various useful properties of F(q, τ) and how it relates to the usual dynamic structure factor. In fact, working directly in the τ domain is advantageous for many reasons and opens up multiple avenues for future applications.