Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2402.01558 · Physical Review Letters

Confinement in the Transverse Field Ising Model on the Heavy Hex Lattice

Joseph Tindall, Dries Sels

Published Oct 29, 2024• 5 citations
arXivHTMLPDF

Abstract

Inspired by a recent quantum computing experiment [Y. Kim et al., Nature (London), 618, 500-5 (2023)NATUAS0028-083610.1038/s41586-023-06096-3], we study the emergence of confinement in the transverse field Ising model on a decorated hexagonal lattice. Using an infinite tensor network state optimized with belief propagation we show how a quench from a broken symmetry state leads to striking nonthermal behavior underpinned by persistent oscillations and saturation of the entanglement entropy. We explain this phenomenon by constructing a minimal model based on the confinement of elementary excitations. Our model is in excellent agreement with our numerical results. For quenches to larger values of the transverse field and/or from nonsymmetry broken states, our numerical results display the expected signatures of thermalization: a linear growth of entanglement entropy in time, propagation of correlations, and the saturation of observables to their thermal averages. These results provide a physical explanation for the unexpected classical simulability of the quantum dynamics.

DOI
10.1103/physrevlett.133.180402

Related papers

Entanglement dynamics in confining spin chains

Stefano Scopa, Pasquale Calabrese, Alvise Bastianello

The confinement of elementary excitations induces distinctive features in the non-equilibrium quench dynamics. One of the most remarkable is the suppression of entanglement entropy, which in several instances turns out to oscillate rather than grow indefinitely. While the qualitative physical origin of this behavior is clear, till now no quantitative understanding away from the field theory limit was available. Here we investigate this problem in the weak quench limit, when mesons are excited at rest, hindering entropy growth and exhibiting persistent oscillations. We provide analytical predictions of the entire entanglement dynamics based on a Gaussian approximation of the many-body state, which captures numerical data with great accuracy and is further simplified to a semiclassical quasiparticle picture in the regime of weak confinement. Our methods are valid in general and we apply explicitly to two prototypical models: the Ising chain in a tilted field and the experimentally relevant long-range Ising model.

arXiv:2111.11483•Physical review. B./Physical review. B•Mar 16, 2022•29 citations

Simulating dynamics of the two-dimensional transverse-field Ising model: a comparative study of large-scale classical numerics

Joseph Vovrosh, Sergi Julià-Farré, Wladislaw Krinitsin +11 more

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.

quant-pharXiv:2511.19340•Physical Review Research•Nov 25, 2025•1 citation

Constrained dynamics and confinement in the two-dimensional quantum Ising model

Luka Pavešič, Daniel Jaschke, Simone Montangero

We investigate the dynamics of the quantum Ising model on two-dimensional square lattices up to 16 $\ifmmode\times\else\texttimes\fi{}$ 16 spins. In the ordered phase, the model exhibits dynamically constrained dynamics, leading to confinement of elementary excitations and slow thermalization. After demonstrating the signatures of confinement, we probe the dynamics of interfaces in the constrained regime through sudden quenches of product states with domains of opposite magnetization. We find that the nature of excitations can be captured by perturbation theory throughout the confining regime and identify the crossover to the deconfining regime. We systematically explore the effect of the transverse field on the modes propagating along flat interfaces and investigate the crossover from resonant to diffusive melting of a square of flipped spins embedded in a larger lattice.

arXiv:2406.11979•Physical review. B./Physical review. B•Apr 10, 2025•9 citations

Confinement and entanglement dynamics on a digital quantum computer

Joseph Vovrosh, Johannes Knolle

Confinement describes the phenomenon when the attraction between two particles grows with their distance, most prominently found in quantum chromodynamics (QCD) between quarks. In condensed matter physics, confinement can appear in quantum spin chains, for example, in the one dimensional transverse field Ising model (TFIM) with an additional longitudinal field, famously observed in the quantum material cobalt niobate or in optical lattices. Here, we establish that state-of-the-art quantum computers have reached capabilities to simulate confinement physics in spin chains. We report quantitative confinement signatures of the TFIM on an IBM quantum computer observed via two distinct velocities for information propagation from domain walls and their mesonic bound states. We also find the confinement induced slow down of entanglement spreading by implementing randomized measurement protocols for the second order Rényi entanglement entropy. Our results are a crucial step for probing non-perturbative interacting quantum phenomena on digital quantum computers beyond the capabilities of classical hardware.

arXiv:2001.03044•Scientific Reports•Jun 2, 2021•74 citations

Confinement and Entanglement Dynamics on a Digital Quantum Computer

Joseph Vovrosh, Johannes Knolle

Abstract Confinement describes the phenomenon when the attraction between two particles grows with their distance, most prominently found in quantum chromodynamics (QCD) between quarks. In condensed matter physics, confinement can appear in quantum spin chains, for example, in the one dimensional transverse field Ising model (TFIM) with an additional longitudinal field, famously observed in the quantum material cobalt niobate or in optical lattices. Here, we establish that state-of-the-art quantum computers have reached quantum simulation capabilities to explore confinement physics in spin chains. We report quantitative confinement signatures of the TFIM on an IBM quantum computer observed via two distinct velocities for information propagation from domain walls and their mesonic bound states. We also find the confinement induced slow down of entanglement spreading by implementing randomized measurement protocols for the second order Renyi entanglement entropy. Our results are a crucial step for probing non-perturbative interacting quantum phenomena on digital quantum computers beyond the capabilities of classical hardware.

•Sep 8, 2020•10 citations

Quantum quench dynamics in the transverse field Ising model of fully connected spins

Ingo Homrighausen

The transverse field Ising model of N fully connected spins provides a mathematical accessible, yet non-trivial, system to investigate the non-thermal dynamics after a sudden quantum quench. The focus in this thesis is three-fold. First, in a 1/N expansion the dynamics of the mean magnetization within its variance is obtained analytically in the large N limit. The variance constitutes a leading order correction to the mean field limit at infinite N. While mean field theories are one of the most accessible approximations to deal with the quantum complexity, its validity in time out of equilibrium has to be checked on a case by case basis. By studying the dynamics of the variance, we see that the mean field approximation can break down at surprisingly early times scaling with the square root of N. The underlying mechanism in the fully connected Ising model is identified as a dephasing effect between nearby effective orbits. The dynamics of the magnetization variance discriminates four qualitatively different regions in the dynamical phase diagram that are indistinguishable by just looking at the mean magnetization. Second, bipartite entanglement between disjoint groups of the spins in the fully connected Ising model is studied. In the large N limit the entanglement Hamiltonian is shown to be a harmonic oscillator. A quantitative relation between the angular frequency of the oscillator and a squeezing measure of the spin state is derived. Explicit expressions for the Renyi entanglement entropies are given. In contrast to the variance, which vanishes in the thermodynamic limit, the entanglement entropy saturate to an N-independent positive value. The dynamics of the entanglement Hamiltonian and the entanglement entropies is discussed. Linear growth as a function of time is linked to the unstable fixed point and the homoclinic orbit of a pitchfork bifurcation, while logarithmic increase is explained due to a dephasing mechanism of nearby effective orbits. A generalization to other models that permit a semiclassical effective description (on a possibly higher dimensional phase space) is given. As a byproduct the connection between entangled Gaussian states and the symplectic eigenvalues of its reduced covariance is derived. An upper bound on the von Neumann entanglement entropy of any bipartite pure state is given in terms of the symplectic properties of its reduced covariance. Third, the existence of dynamical phase transitions (DPTs) in the fully connected Ising model is argued by two different approaches. DPTs are defined as non-analyticities in the time-dependence of the return probability rate function in the thermodynamic limit. Like entanglement entropy, the behavior of the return rate function is an observable-independent quantity used to characterize non-equilibrium states. A numerically based connection between the DPT and the analytic property of the infinite time averaged order parameter as a function of quench strength is discussed. The concept of Feynman history states is used to devise the notion of history entanglement. An information theoretic interpretation of the history entanglement entropy as the precision of a quantum clock follows from Holevo’s bound. The dependence of this entanglement entropy on the observation time constitutes an operator-independent characterization of the quantum dynamics in closed systems. For observation times comparable to the Heisenberg time, the history entanglement becomes sensitive to the statistics of the energy gaps, and thus distinguishes integrable and non-integrable dynamics. This is demonstrated for two examples, a random matrix ensemble, and a one-dimensional hardcore boson model with integrability breaking next nearest neighbor interaction.

•Jan 1, 2021•0 citations
Inspired by a recent quantum computing experiment [Y. Kim et al., Nature (London), 618, 500-5 (2023)NATUAS0028-083610.1038/s41586-023-06096-3], we study the emergence of confinement in the transverse field Ising model on a decorated hexagonal lattice. Using an infinite tensor network state optimized with belief propagation we show how a quench from a broken symmetry state leads to striking nonthermal behavior underpinned by persistent oscillations and saturation of the entanglement entropy. We explain this phenomenon by constructing a minimal model based on the confinement of elementary excitations. Our model is in excellent agreement with our numerical results. For quenches to larger values of the transverse field and/or from nonsymmetry broken states, our numerical results display the expected signatures of thermalization: a linear growth of entanglement entropy in time, propagation of correlations, and the saturation of observables to their thermal averages. These results provide a physical explanation for the unexpected classical simulability of the quantum dynamics.
The confinement of elementary excitations induces distinctive features in the non-equilibrium quench dynamics. One of the most remarkable is the suppression of entanglement entropy, which in several instances turns out to oscillate rather than grow indefinitely. While the qualitative physical origin of this behavior is clear, till now no quantitative understanding away from the field theory limit was available. Here we investigate this problem in the weak quench limit, when mesons are excited at rest, hindering entropy growth and exhibiting persistent oscillations. We provide analytical predictions of the entire entanglement dynamics based on a Gaussian approximation of the many-body state, which captures numerical data with great accuracy and is further simplified to a semiclassical quasiparticle picture in the regime of weak confinement. Our methods are valid in general and we apply explicitly to two prototypical models: the Ising chain in a tilted field and the experimentally relevant long-range Ising model.
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.
We investigate the dynamics of the quantum Ising model on two-dimensional square lattices up to 16 $\ifmmode\times\else\texttimes\fi{}$ 16 spins. In the ordered phase, the model exhibits dynamically constrained dynamics, leading to confinement of elementary excitations and slow thermalization. After demonstrating the signatures of confinement, we probe the dynamics of interfaces in the constrained regime through sudden quenches of product states with domains of opposite magnetization. We find that the nature of excitations can be captured by perturbation theory throughout the confining regime and identify the crossover to the deconfining regime. We systematically explore the effect of the transverse field on the modes propagating along flat interfaces and investigate the crossover from resonant to diffusive melting of a square of flipped spins embedded in a larger lattice.
Confinement describes the phenomenon when the attraction between two particles grows with their distance, most prominently found in quantum chromodynamics (QCD) between quarks. In condensed matter physics, confinement can appear in quantum spin chains, for example, in the one dimensional transverse field Ising model (TFIM) with an additional longitudinal field, famously observed in the quantum material cobalt niobate or in optical lattices. Here, we establish that state-of-the-art quantum computers have reached capabilities to simulate confinement physics in spin chains. We report quantitative confinement signatures of the TFIM on an IBM quantum computer observed via two distinct velocities for information propagation from domain walls and their mesonic bound states. We also find the confinement induced slow down of entanglement spreading by implementing randomized measurement protocols for the second order Rényi entanglement entropy. Our results are a crucial step for probing non-perturbative interacting quantum phenomena on digital quantum computers beyond the capabilities of classical hardware.
Abstract Confinement describes the phenomenon when the attraction between two particles grows with their distance, most prominently found in quantum chromodynamics (QCD) between quarks. In condensed matter physics, confinement can appear in quantum spin chains, for example, in the one dimensional transverse field Ising model (TFIM) with an additional longitudinal field, famously observed in the quantum material cobalt niobate or in optical lattices. Here, we establish that state-of-the-art quantum computers have reached quantum simulation capabilities to explore confinement physics in spin chains. We report quantitative confinement signatures of the TFIM on an IBM quantum computer observed via two distinct velocities for information propagation from domain walls and their mesonic bound states. We also find the confinement induced slow down of entanglement spreading by implementing randomized measurement protocols for the second order Renyi entanglement entropy. Our results are a crucial step for probing non-perturbative interacting quantum phenomena on digital quantum computers beyond the capabilities of classical hardware.
The transverse field Ising model of N fully connected spins provides a mathematical accessible, yet non-trivial, system to investigate the non-thermal dynamics after a sudden quantum quench. The focus in this thesis is three-fold. First, in a 1/N expansion the dynamics of the mean magnetization within its variance is obtained analytically in the large N limit. The variance constitutes a leading order correction to the mean field limit at infinite N. While mean field theories are one of the most accessible approximations to deal with the quantum complexity, its validity in time out of equilibrium has to be checked on a case by case basis. By studying the dynamics of the variance, we see that the mean field approximation can break down at surprisingly early times scaling with the square root of N. The underlying mechanism in the fully connected Ising model is identified as a dephasing effect between nearby effective orbits. The dynamics of the magnetization variance discriminates four qualitatively different regions in the dynamical phase diagram that are indistinguishable by just looking at the mean magnetization. Second, bipartite entanglement between disjoint groups of the spins in the fully connected Ising model is studied. In the large N limit the entanglement Hamiltonian is shown to be a harmonic oscillator. A quantitative relation between the angular frequency of the oscillator and a squeezing measure of the spin state is derived. Explicit expressions for the Renyi entanglement entropies are given. In contrast to the variance, which vanishes in the thermodynamic limit, the entanglement entropy saturate to an N-independent positive value. The dynamics of the entanglement Hamiltonian and the entanglement entropies is discussed. Linear growth as a function of time is linked to the unstable fixed point and the homoclinic orbit of a pitchfork bifurcation, while logarithmic increase is explained due to a dephasing mechanism of nearby effective orbits. A generalization to other models that permit a semiclassical effective description (on a possibly higher dimensional phase space) is given. As a byproduct the connection between entangled Gaussian states and the symplectic eigenvalues of its reduced covariance is derived. An upper bound on the von Neumann entanglement entropy of any bipartite pure state is given in terms of the symplectic properties of its reduced covariance. Third, the existence of dynamical phase transitions (DPTs) in the fully connected Ising model is argued by two different approaches. DPTs are defined as non-analyticities in the time-dependence of the return probability rate function in the thermodynamic limit. Like entanglement entropy, the behavior of the return rate function is an observable-independent quantity used to characterize non-equilibrium states. A numerically based connection between the DPT and the analytic property of the infinite time averaged order parameter as a function of quench strength is discussed. The concept of Feynman history states is used to devise the notion of history entanglement. An information theoretic interpretation of the history entanglement entropy as the precision of a quantum clock follows from Holevo’s bound. The dependence of this entanglement entropy on the observation time constitutes an operator-independent characterization of the quantum dynamics in closed systems. For observation times comparable to the Heisenberg time, the history entanglement becomes sensitive to the statistics of the energy gaps, and thus distinguishes integrable and non-integrable dynamics. This is demonstrated for two examples, a random matrix ensemble, and a one-dimensional hardcore boson model with integrability breaking next nearest neighbor interaction.