This paper investigates the transverse Ising model on a discretization of two-dimensional anti--de Sitter space. We use classical and quantum algorithms to simulate real-time evolution and measure out-of-time-ordered correlators (OTOC). The latter can probe thermalization and scrambling of quantum information under time evolution. We compared tensor network-based methods both with simulation on gate-based superconducting quantum devices and analog quantum simulation using Rydberg arrays. While studying this system's thermalization properties, we observed different regimes depending on the radius of curvature of the space. In particular, we find a region of parameter space where the thermalization time depends only logarithmically on the number of degrees of freedom.
This paper investigates the transverse Ising model on a discretization of two-dimensional anti--de Sitter space. We use classical and quantum algorithms to simulate real-time evolution and measure out-of-time-ordered correlators (OTOC). The latter can probe thermalization and scrambling of quantum information under time evolution. We compared tensor network-based methods both with simulation on gate-based superconducting quantum devices and analog quantum simulation using Rydberg arrays. While studying this system's thermalization properties, we observed different regimes depending on the radius of curvature of the space. In particular, we find a region of parameter space where the thermalization time depends only logarithmically on the number of degrees of freedom.
Quantum Ising Model on (2+1)−Dimensional Anti−de Sitter Space using Tensor Networks
We study the quantum Ising model on (2+1)-dimensional anti-de Sitter space using Matrix Product States (MPS) and Matrix Product Operators (MPOs). We explore the bulk phase diagram of the theory on regular tessellations of hyperbolic space with coordination number seven and find disordered and ordered phases separated by a phase transition. We find that the boundary-boundary spin correlation function exhibits power law scaling deep in the disordered phase of the Ising model consistent with holography. At the critical point, we find the boundary entanglement entropy scales logarithmically with subsystem size but away from this, we see a linear scaling. In comparison, the full system exhibits a volume law scaling, which is expected in chaotic and/or highly connected systems. We also measure Out of time Ordered Correlators (OTOCs) to explore the scrambling behavior of the theory.
The quantum dynamics of many-qubit systems is an outstanding problem that has recently driven significant advances in both numerical methods and programmable quantum processing units. In this work, we employ a comprehensive toolbox of state-of-the-art numerical approaches to classically simulate the dynamics of the two-dimensional transverse field Ising model. Our methods include three different tensor network techniques -- matrix product states, tree-tensor networks, and two-dimensional tensor-networks under the belief propagation approximation -- as well as time-dependent variational Monte Carlo with Neural Quantum States. We focus on two paradigmatic dynamical protocols: (i) quantum annealing through a critical point and (ii) post-quench dynamics. Our extensive results show the quantitative predictions of various state-of-the-art numerical methods providing a benchmark for future numerical investigations and experimental studies with the aim to push the limitations on classical and QPUs. In particular, our work connects classical simulability to different regimes associated with quantum dynamics in Rydberg arrays - namely, quasi-adiabatic dynamics, the Kibble-Zurek mechanism, and quantum quenches.
Out-of-time-ordered correlators (OTOCs) are a key observable in a wide range of interconnected fields including many-body physics, quantum information science, and quantum gravity. Measuring OTOCs using near-term quantum simulators will extend our ability to explore fundamental aspects of these fields and the subtle connections between them. Here, we demonstrate an experimental method to measure OTOCs at finite temperatures and use the method to study their temperature dependence. These measurements are performed on a digital quantum computer running a simulation of the transverse field Ising model. Our flexible method, based on the creation of a thermofield double state, can be extended to other models and enables us to probe the OTOC's temperature-dependent decay rate. Measuring this decay rate opens up the possibility of testing the fundamental temperature-dependent bounds on quantum information scrambling.
The efficient numerical simulation of nonequilibrium real-time evolution in isolated quantum matter constitutes a key challenge for current computational methods. This holds in particular in the regime of two spatial dimensions, whose experimental exploration is currently pursued with strong efforts in quantum simulators. In this work we present a versatile and efficient machine learning inspired approach based on a recently introduced artificial neural network encoding of quantum many-body wave functions. We identify and resolve key challenges for the simulation of time evolution, which previously imposed significant limitations on the accurate description of large systems and long-time dynamics. As a concrete example, we study the dynamics of the paradigmatic two-dimensional transverse-field Ising model, as recently also realized experimentally in systems of Rydberg atoms. Calculating the nonequilibrium real-time evolution across a broad range of parameters, we, for instance, observe collapse and revival oscillations of ferromagnetic order and demonstrate that the reached timescales are comparable to or exceed the capabilities of state-of-the-art tensor network methods.
A closed quantum system never forgets its initial state, but the encoded information can get scrambled and become inaccessible without measuring a large fraction of all the system degrees of freedom. This scrambling can be diagnosed by studying the spatial spreading of initially local operators under the Heisenberg time evolution, and the decay of the out-of-time-ordered correlators (OTOC). What insights can OTOCs provide to understand the dynamics of quantum many-body systems? What are the characteristic behaviors of OTOCs during the time evolution? How is information scrambling affected by the dissipation in open quantum many-body systems? We first study slow scrambling in many-body localized systems via calculating various correlators, two-point retarded correlators and OTOCs. Comparing with retarded correlators, OTOCs provide more information about the dynamics. We find that disorder slows and partially halts the onset of information scrambling. Instead of ballistic spreading, propagation of information forms a logarithmic light cone. Next, we study the finite-size scaling of OTOCs at late times in generic thermalizing quantum many-body systems. When energy is conserved, the late-time saturation value of the OTOC of generic traceless local operators scales as an inverse polynomial in the system size. This is in contrast to the inverse exponential scaling expected for chaotic dynamics without energy conservation. We also study information scrambling in open quantum many-body systems. We define a dissipative version of OTOC and study its behaviors in a prototypical chaotic quantum chain with dissipation. We find that dissipation leads to not only the overall decay of the scrambled information due to leaking, but also structural changes so that the information light cone can only reach a finite distance even when the effect of overall decay is removed. Finally, we construct a family of local Hamiltonians for understanding the asymmetric information scrambling. Our models live on a one-dimensional lattice and exhibit asymmetric butterfly light cone between the left and right spatial directions.
Abstract Out-of-time-order correlators (OTOCs) have emerged as the preeminent theoretical and experimental probe of quantum information scrambling — the delocalization of quantum information across many-body degrees of freedom. Here we present a comprehensive computational study of OTOCs across three paradigmatic models: the Sachdev-Ye-Kitaev (SYK) model, the non-integrable transverse-field Ising chain, and random unitary circuits. Through exact diagonalization and large-N analytics, we extract Lyapunov exponents, scrambling times, and butterfly velocities, demonstrating that all chaotic systems studied approach, but do not saturate, the Maldacena-Shenker-Stanford (MSS) bound λ_L ≤ 2π/β. Crucially, we introduce the Spectral Out-of-Time-Order Correlator (SOTOC), a novel frequency-domain diagnostic that decomposes information scrambling into its constituent spectral modes. The SOTOC achieves a 32% improvement in Lyapunov exponent extraction precision over conventional log-linear fitting and maintains robust performance at noise levels σ > 0.30, where traditional methods fail. Our phase diagram reveals a sharp transition in scrambling dynamics at h/J ≈ 1.0 in the TFIM, coinciding with the onset of quantum chaos as measured by level-spacing statistics. These results establish SOTOC as a practical and theoretically grounded tool for probing scrambling in near-term quantum devices and holographic systems.