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:2108.11051 · Communications Physics

Tensor-network study of correlation-spreading dynamics in the two-dimensional Bose-Hubbard model

Ryui Kaneko, Ippei Danshita

Published Mar 21, 2022• 30 citations
arXivHTMLPDF

Abstract

Abstract Recent developments in analog quantum simulators based on cold atoms and trapped ions call for cross-validating the accuracy of quantum-simulation experiments with use of quantitative numerical methods; however, it is particularly challenging for dynamics of systems with more than one spatial dimension. Here we demonstrate that a tensor-network method running on classical computers is useful for this purpose. We specifically analyze real-time dynamics of the two-dimensional Bose-Hubbard model after a sudden quench starting from the Mott insulator by means of the infinite projected entangled pair state algorithm. Calculated single-particle correlation functions are found to be in good agreement with a recent experiment [Y. Takasu et al., Sci. Adv. 6, eaba9255 (2020)]. By estimating the phase and group velocities from the single-particle and density-density correlation functions, we predict how these velocities vary in the moderate interaction region, which serves as a quantitative benchmark for future experiments.

DOI
10.1038/s42005-022-00848-9

Related papers

Tensor network simulation of the quantum Kibble-Zurek quench from the Mott to the superfluid phase in the two-dimensional Bose-Hubbard model

Jacek Dziarmaga, Jakub M. Mazur

Quantum simulations of the Bose-Hubbard model (BHM) at commensurate filling can follow spreading of correlations after a sudden quench for times long enough to estimate their propagation velocities. In this work we perform tensor network simulation of the quantum Kibble-Zurek (KZ) ramp from the Mott towards the superfluid phase in the square lattice BHM and demonstrate that even relatively short ramp/quench times allow one to test the power laws predicted by the KZ mechanism (KZM). They can be verified for the correlation length and the excitation energy but the most reliable test is based on the KZM scaling hypothesis for the single-particle correlation function: scaled correlation functions for different quench times evaluated at the same scaled time collapse to the same scaling function of the scaled distance. The scaling of the space and time variables is done according to the KZ power laws.

arXiv:2302.13347•Physical review. B./Physical review. B•Apr 24, 2023•15 citations

Tensor Network Annealing Algorithm for Two-Dimensional Thermal States

Augustine Kshetrimayum, Matteo Rizzi, Jens Eisert +1 more

Tensor network methods have become a powerful class of tools to capture strongly correlated matter, but methods to capture the experimentally ubiquitous family of models at finite temperature beyond one spatial dimension are largely lacking. We introduce a tensor network algorithm able to simulate thermal states of two-dimensional quantum lattice systems in the thermodynamic limit. The method develops instances of projected entangled pair states and projected entangled pair operators for this purpose. It is the key feature of this algorithm to resemble the cooling down of the system from an infinite temperature state until it reaches the desired finite-temperature regime. As a benchmark, we study the finite-temperature phase transition of the Ising model on an infinite square lattice, for which we obtain remarkable agreement with the exact solution. We then turn to study the finite-temperature Bose-Hubbard model in the limits of two (hard-core) and three bosonic modes per site. Our technique can be used to support the experimental study of actual effectively two-dimensional materials in the laboratory, as well as to benchmark optical lattice quantum simulators with ultracold atoms.

arXiv:1809.08258•Physical Review Letters•Feb 22, 2019•84 citations

Machine learning out-of-equilibrium correlations in the Bose-Hubbard model

Ali Mokhtari-Jazi, Malcolm P. Kennett

Calculating the out-of-equilibrium dynamics of many-body quantum systems theoretically is a challenging problem. Essentially exact results can be obtained for the out-of-equilibrium correlations in the Bose-Hubbard model in one dimension, but higher dimensions require approximate methods. One such method is the two-particle irreducible strong coupling (2PISC) approach [M. R. C. Fitzpatrick and M. P. Kennett, Nucl. Phys. B 930, 1 (2018)]. Calculations of the single-particle correlations using this method yield values of the velocity for correlation spreading that match well with exact methods in one dimension and experiments in one and two dimensions. However, the 2PISC method is less accurate for determining the amplitude of correlations, especially in the regime where interactions are not very strong. Viewing the calculation of the single-particle correlations as an image correction problem, we train a neural network (NN) to take input from the 2PISC approach to reproduce the output of exact diagonalization calculations. We show that the predictions of the NN improve on 2PISC results for parameters outside the training region. Our approach is not specific to the Bose-Hubbard model and may find application to the out-of-equilibrium dynamics of other quantum many-body systems.

arXiv:2409.18225•Physical Review A•Apr 18, 2025•0 citations

Lieb-Robinson-like bounds for the Mott phase of the Bose-Hubbard model in dimensions greater than one

Ali Mokhtari-Jazi, Matthew Fitzpatrick, Malcolm P. Kennett

Lieb-Robinson and related bounds set an upper limit on the rate of spreading of information in non-relativistic quantum systems. Experimentally, they have been observed in the spreading of correlations in the Bose-Hubbard model after a quantum quench. Using a recently developed two particle irreducible (2PI) strong coupling approach to out-of-equilibrium dynamics in the Bose-Hubbard model we calculate both the group and phase velocities for the spreading of single-particle correlations in one, two and three dimensions as a function of interaction strength. Our results are in quantitative agreement with measurements of the velocities for the spreading of single particle correlations in both the one and two dimensional Bose-Hubbard model realized with ultra-cold atoms. They also agree with the claim that the phase velocity rather than the group velocity was observed in recent experiments in two dimensions. We demonstrate that there can be large differences between the phase and group velocities for the spreading of correlations and also explore the variation of the anisotropy in the velocity at which correlations spread across the phase diagram of the Bose-Hubbard model. Our results establish the 2PI strong coupling approach as a powerful tool to study out-of-equilibrium dynamics in the Bose-Hubbard model in dimensions greater than one.

•Jul 31, 2020•0 citations

Space-time correlations in the Bose Hubbard model after a quantum quench

Matthew Fitzpatrick, Malcolm P. Kennett

We study correlations in space and time after a quantum quench in the Bose Hubbard model. We obtain a two particle irreducible effective action within the Schwinger-Keldysh formalism that allows us to study dynamics in both the superfluid and Mott insulating phases. We derive equations of motion for both the superfluid order parameter and two-point correlation functions. We discuss the numerical solutions of these equations and calculate the quasimomentum distribution for quenches in the Mott phase. We demonstrate light-cone like spreading of correlations in the Mott phase in both one and two dimensions and calculate the propagation velocity in each dimension.

•Jun 13, 2016•1 citation

Phase and group velocities for correlation spreading in the Mott phase of the Bose-Hubbard model in dimensions greater than one

Ali Mokhtari-Jazi, Matthew Fitzpatrick, Malcolm P. Kennett

Lieb-Robinson and related bounds set an upper limit on the speed at which information propagates in nonrelativistic quantum systems. Experimentally, light-cone-like spreading has been observed for correlations in the Bose-Hubbard model (BHM) after a quantum quench. Using a two-particle irreducible (2PI) strong-coupling approach to out-of-equilibrium dynamics in the BHM we calculate both the group and phase velocities for the spreading of single-particle correlations in one, two, and three dimensions as a function of interaction strength. Our results are in quantitative agreement with measurements of the speed of spreading of single-particle correlations in both the one- and two-dimensional BHM realized with ultracold atoms. They also are consistent with the claim that the phase velocity rather than the group velocity was observed in recent experiments in two dimensions. We demonstrate that there can be large differences between the phase and group velocities for the spreading of correlations and explore how the anisotropy in the velocity varies across the phase diagram of the BHM. Our results establish the 2PI strong-coupling approach as a powerful tool to study out-of-equilibrium dynamics in the BHM in dimensions greater than one.

arXiv:2008.00015•Physical Review A•Feb 26, 2021•12 citations
Abstract Recent developments in analog quantum simulators based on cold atoms and trapped ions call for cross-validating the accuracy of quantum-simulation experiments with use of quantitative numerical methods; however, it is particularly challenging for dynamics of systems with more than one spatial dimension. Here we demonstrate that a tensor-network method running on classical computers is useful for this purpose. We specifically analyze real-time dynamics of the two-dimensional Bose-Hubbard model after a sudden quench starting from the Mott insulator by means of the infinite projected entangled pair state algorithm. Calculated single-particle correlation functions are found to be in good agreement with a recent experiment [Y. Takasu et al., Sci. Adv. 6, eaba9255 (2020)]. By estimating the phase and group velocities from the single-particle and density-density correlation functions, we predict how these velocities vary in the moderate interaction region, which serves as a quantitative benchmark for future experiments.
Quantum simulations of the Bose-Hubbard model (BHM) at commensurate filling can follow spreading of correlations after a sudden quench for times long enough to estimate their propagation velocities. In this work we perform tensor network simulation of the quantum Kibble-Zurek (KZ) ramp from the Mott towards the superfluid phase in the square lattice BHM and demonstrate that even relatively short ramp/quench times allow one to test the power laws predicted by the KZ mechanism (KZM). They can be verified for the correlation length and the excitation energy but the most reliable test is based on the KZM scaling hypothesis for the single-particle correlation function: scaled correlation functions for different quench times evaluated at the same scaled time collapse to the same scaling function of the scaled distance. The scaling of the space and time variables is done according to the KZ power laws.
Tensor network methods have become a powerful class of tools to capture strongly correlated matter, but methods to capture the experimentally ubiquitous family of models at finite temperature beyond one spatial dimension are largely lacking. We introduce a tensor network algorithm able to simulate thermal states of two-dimensional quantum lattice systems in the thermodynamic limit. The method develops instances of projected entangled pair states and projected entangled pair operators for this purpose. It is the key feature of this algorithm to resemble the cooling down of the system from an infinite temperature state until it reaches the desired finite-temperature regime. As a benchmark, we study the finite-temperature phase transition of the Ising model on an infinite square lattice, for which we obtain remarkable agreement with the exact solution. We then turn to study the finite-temperature Bose-Hubbard model in the limits of two (hard-core) and three bosonic modes per site. Our technique can be used to support the experimental study of actual effectively two-dimensional materials in the laboratory, as well as to benchmark optical lattice quantum simulators with ultracold atoms.
Calculating the out-of-equilibrium dynamics of many-body quantum systems theoretically is a challenging problem. Essentially exact results can be obtained for the out-of-equilibrium correlations in the Bose-Hubbard model in one dimension, but higher dimensions require approximate methods. One such method is the two-particle irreducible strong coupling (2PISC) approach [M. R. C. Fitzpatrick and M. P. Kennett, Nucl. Phys. B 930, 1 (2018)]. Calculations of the single-particle correlations using this method yield values of the velocity for correlation spreading that match well with exact methods in one dimension and experiments in one and two dimensions. However, the 2PISC method is less accurate for determining the amplitude of correlations, especially in the regime where interactions are not very strong. Viewing the calculation of the single-particle correlations as an image correction problem, we train a neural network (NN) to take input from the 2PISC approach to reproduce the output of exact diagonalization calculations. We show that the predictions of the NN improve on 2PISC results for parameters outside the training region. Our approach is not specific to the Bose-Hubbard model and may find application to the out-of-equilibrium dynamics of other quantum many-body systems.
Lieb-Robinson and related bounds set an upper limit on the rate of spreading of information in non-relativistic quantum systems. Experimentally, they have been observed in the spreading of correlations in the Bose-Hubbard model after a quantum quench. Using a recently developed two particle irreducible (2PI) strong coupling approach to out-of-equilibrium dynamics in the Bose-Hubbard model we calculate both the group and phase velocities for the spreading of single-particle correlations in one, two and three dimensions as a function of interaction strength. Our results are in quantitative agreement with measurements of the velocities for the spreading of single particle correlations in both the one and two dimensional Bose-Hubbard model realized with ultra-cold atoms. They also agree with the claim that the phase velocity rather than the group velocity was observed in recent experiments in two dimensions. We demonstrate that there can be large differences between the phase and group velocities for the spreading of correlations and also explore the variation of the anisotropy in the velocity at which correlations spread across the phase diagram of the Bose-Hubbard model. Our results establish the 2PI strong coupling approach as a powerful tool to study out-of-equilibrium dynamics in the Bose-Hubbard model in dimensions greater than one.
We study correlations in space and time after a quantum quench in the Bose Hubbard model. We obtain a two particle irreducible effective action within the Schwinger-Keldysh formalism that allows us to study dynamics in both the superfluid and Mott insulating phases. We derive equations of motion for both the superfluid order parameter and two-point correlation functions. We discuss the numerical solutions of these equations and calculate the quasimomentum distribution for quenches in the Mott phase. We demonstrate light-cone like spreading of correlations in the Mott phase in both one and two dimensions and calculate the propagation velocity in each dimension.
Lieb-Robinson and related bounds set an upper limit on the speed at which information propagates in nonrelativistic quantum systems. Experimentally, light-cone-like spreading has been observed for correlations in the Bose-Hubbard model (BHM) after a quantum quench. Using a two-particle irreducible (2PI) strong-coupling approach to out-of-equilibrium dynamics in the BHM we calculate both the group and phase velocities for the spreading of single-particle correlations in one, two, and three dimensions as a function of interaction strength. Our results are in quantitative agreement with measurements of the speed of spreading of single-particle correlations in both the one- and two-dimensional BHM realized with ultracold atoms. They also are consistent with the claim that the phase velocity rather than the group velocity was observed in recent experiments in two dimensions. We demonstrate that there can be large differences between the phase and group velocities for the spreading of correlations and explore how the anisotropy in the velocity varies across the phase diagram of the BHM. Our results establish the 2PI strong-coupling approach as a powerful tool to study out-of-equilibrium dynamics in the BHM in dimensions greater than one.