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:2302.13347 · Physical review. B./Physical review. B

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

Published Apr 24, 2023• 15 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physrevb.107.144510

Related papers

Dynamics of the quantum phase transition in the one-dimensional Bose-Hubbard model: Excitations and correlations induced by a quench

Bartłomiej Gardas, Jacek Dziarmaga, Wojciech Hubert Zurek

The ground state of the one-dimensional Bose-Hubbard model at unit filling undergoes the Mott-superfluid quantum phase transition. It belongs to the Kosterlitz-Thouless universality class with an exponential divergence of the correlation length in place of the usual power law. We present numerical simulations of a linear quench both from the Mott insulator to superfluid and back. The results satisfy the scaling hypothesis that follows from the Kibble-Zurek mechanism (KZM). In the superfluid-to-Mott quenches there is no significant excitation in the superfluid phase despite its gaplessness. Since all critical superfluid ground states are qualitatively similar, the excitation begins to build up only after crossing the critical point when the ground state begins to change fundamentally. The last process falls into the KZM framework.

arXiv:1612.05084•Physical review. B./Physical review. B•Mar 24, 2017•45 citations

Quench in the 1D Bose-Hubbard mode

Bartłomiej Gardas, Jacek Dziarmaga, Wojciech Hubert Zurek

The ground state of the one-dimensional Bose-Hubbard model at unit filling undergoes the Mott-superfluid quantum phase transition. It belongs to the Kosterlitz-Thouless universality class with an exponential divergence of the correlation length in place of the usual power law. We present numerical simulations of a linear quench both from the Mott insulator to superfluid and back. The results satisfy the scaling hypothesis that follows from the Kibble-Zurek mechanism (KZM). In the superfluid-to-Mott quenches there is no significant excitation in the superfluid phase despite its gaplessness. Since all critical superfluid ground states are qualitatively similar, the excitation begins to build up only after crossing the critical point when the ground state begins to change fundamentally. The last process falls into the KZM framework.

•Dec 15, 2016•1 citation

Kibble-Zurek dynamical scaling hypothesis in the Google analog-digital quantum simulator of the X X model

Yintai Zhang, Francis A. Bayocboc, Jacek Dziarmaga

State-of-the-art tensor networks are employed to simulate the Hamiltonian ramp in the analog-digital quantum simulation of the quantum phase transition to the quasi-long-range ordered phase of the two-dimensional square-lattice $XX$ model [T. I. Andersen et al. Nature (London) 638, 79 (2025)]. We focus on the quantum Kibble-Zurek (KZ) mechanism near the quantum critical point. Using the infinite projected entangled pair state, we simulate an infinite lattice and demonstrate the KZ scaling hypothesis for the $XX$ correlations across a wide range of ramp times. We use the time-dependent variational principle algorithm to simulate a finite $8\ifmmode\times\else\texttimes\fi{}8$ lattice, similar to the one in the quantum simulation, and find that adiabatic finite-size effects dominate for longer ramp times, where the correlation length's growth with increasing ramp time saturates and the excitation energy's dependence on the ramp time crosses over to a power-law decay characteristic of adiabatic transitions. This finding contradicts the quantum simulation data where the correlation length seems to obey KZ-like power laws, although with modified exponents.

Physical review. B./Physical review. B•Sep 29, 2025•2 citations

Kibble-Zurek dynamical scaling hypothesis in the Google analog-digital quantum simulator of the $XX$ model

Yintai Zhang, Francis A. Bayocboc, Jacek Dziarmaga

State-of-the-art tensor networks are employed to simulate the Hamiltonian ramp in the analog-digital quantum simulation of the quantum phase transition to the quasi-long-range ordered phase of the two-dimensional square-lattice $XX$ model [T.I. Andersen \textit{et al.}, Nature (London) \textbf{638}, 79 (2025)]. We focus on the quantum Kibble-Zurek (KZ) mechanism near the quantum critical point. Using the infinite projected entangled pair state, we simulate an infinite lattice and demonstrate the KZ scaling hypothesis for the $XX$ correlations across a wide range of ramp times. We use the time-dependent variational principle algorithm to simulate a finite $8\times 8$ lattice, similar to the one in the quantum simulation, and find that adiabatic finite-size effects dominate for longer ramp times, where the correlation length's growth with increasing ramp time saturates and the excitation energy's dependence on the ramp time crosses over to a power-law decay characteristic of adiabatic transitions. This finding contradicts the quantum simulation data where the correlation length seems to obey KZ-like power laws, although with modified exponents.

arXiv:2506.10771•Jun 12, 2025•0 citations

Quench in the 1D Bose-Hubbard model: Topological defects and excitations from the Kosterlitz-Thouless phase transition dynamics

Jacek Dziarmaga, Wojciech Hubert Zurek

Kibble-Zurek mechanism (KZM) uses critical scaling to predict density of topological defects and other excitations created in second order phase transitions. We point out that simply inserting asymptotic critical exponents deduced from the immediate vicinity of the critical point to obtain predictions can lead to results that are inconsistent with a more careful KZM analysis based on causality - on the comparison of the relaxation time of the order parameter with the "time distance" from the critical point. As a result, scaling of quench-generated excitations with quench rates can exhibit behavior that is locally (i.e., in the neighborhood of any given quench rate) well approximated by the power law, but with exponents that depend on that rate, and that are quite different from the naive prediction based on the critical exponents relevant for asymptotically long quench times. Kosterlitz-Thouless scaling (that governs e.g. Mott insulator to superfluid transition in the Bose-Hubbard model in one dimension) is investigated as an example of this phenomenon.

Scientific Reports•Aug 5, 2014•31 citations

Quench in 1D Bose-Hubbard Model: Topological Defects and Excitations from Kosterlitz-Thouless Phase Transition Dynamics

Jacek Dziarmaga, Wojciech Hubert Zurek

Kibble-Zurek mechanism (KZM) uses critical scaling to predict density of topological defects and other excitations created in second order phase transitions. We point out that simply inserting asymptotic critical exponents deduced from the immediate vicinity of the critical point to obtain KZM predictions can lead to results that are inconsistent with a more careful KZM-like analysis based on causality -- on the comparison of the relaxation time of the order parameter with the time "distance" from the critical point. As a result, scaling of quench-generated excitations with quench rates can exhibit behavior that is locally (i.e., in the neighborhood of any given quench rate) well approximated by the power law, but with exponents that depend on that rate, and that are quite different from the naive prediction based on the critical exponents relevant for asymptotically long quench times. Kosterlitz-Thouless scaling (that governs e.g. Mott insulator to superfluid transition in the Bose-Hubbard model in one dimension) is investigated as an example of this phenomenon.

arXiv:1312.5139•Dec 18, 2013•1 citation
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.
The ground state of the one-dimensional Bose-Hubbard model at unit filling undergoes the Mott-superfluid quantum phase transition. It belongs to the Kosterlitz-Thouless universality class with an exponential divergence of the correlation length in place of the usual power law. We present numerical simulations of a linear quench both from the Mott insulator to superfluid and back. The results satisfy the scaling hypothesis that follows from the Kibble-Zurek mechanism (KZM). In the superfluid-to-Mott quenches there is no significant excitation in the superfluid phase despite its gaplessness. Since all critical superfluid ground states are qualitatively similar, the excitation begins to build up only after crossing the critical point when the ground state begins to change fundamentally. The last process falls into the KZM framework.
The ground state of the one-dimensional Bose-Hubbard model at unit filling undergoes the Mott-superfluid quantum phase transition. It belongs to the Kosterlitz-Thouless universality class with an exponential divergence of the correlation length in place of the usual power law. We present numerical simulations of a linear quench both from the Mott insulator to superfluid and back. The results satisfy the scaling hypothesis that follows from the Kibble-Zurek mechanism (KZM). In the superfluid-to-Mott quenches there is no significant excitation in the superfluid phase despite its gaplessness. Since all critical superfluid ground states are qualitatively similar, the excitation begins to build up only after crossing the critical point when the ground state begins to change fundamentally. The last process falls into the KZM framework.
State-of-the-art tensor networks are employed to simulate the Hamiltonian ramp in the analog-digital quantum simulation of the quantum phase transition to the quasi-long-range ordered phase of the two-dimensional square-lattice XXXXXX model [T. I. Andersen et al. Nature (London) 638, 79 (2025)]. We focus on the quantum Kibble-Zurek (KZ) mechanism near the quantum critical point. Using the infinite projected entangled pair state, we simulate an infinite lattice and demonstrate the KZ scaling hypothesis for the XXXXXX correlations across a wide range of ramp times. We use the time-dependent variational principle algorithm to simulate a finite $8\ifmmode\times\else\texttimes\fi{}8$ lattice, similar to the one in the quantum simulation, and find that adiabatic finite-size effects dominate for longer ramp times, where the correlation length's growth with increasing ramp time saturates and the excitation energy's dependence on the ramp time crosses over to a power-law decay characteristic of adiabatic transitions. This finding contradicts the quantum simulation data where the correlation length seems to obey KZ-like power laws, although with modified exponents.
Kibble-Zurek dynamical scaling hypothesis in the Google analog-digital quantum simulator of the XXXXXX model
State-of-the-art tensor networks are employed to simulate the Hamiltonian ramp in the analog-digital quantum simulation of the quantum phase transition to the quasi-long-range ordered phase of the two-dimensional square-lattice XXXXXX model [T.I. Andersen \textit{et al.}, Nature (London) \textbf{638}, 79 (2025)]. We focus on the quantum Kibble-Zurek (KZ) mechanism near the quantum critical point. Using the infinite projected entangled pair state, we simulate an infinite lattice and demonstrate the KZ scaling hypothesis for the XXXXXX correlations across a wide range of ramp times. We use the time-dependent variational principle algorithm to simulate a finite 8×88\times 88×8 lattice, similar to the one in the quantum simulation, and find that adiabatic finite-size effects dominate for longer ramp times, where the correlation length's growth with increasing ramp time saturates and the excitation energy's dependence on the ramp time crosses over to a power-law decay characteristic of adiabatic transitions. This finding contradicts the quantum simulation data where the correlation length seems to obey KZ-like power laws, although with modified exponents.
Kibble-Zurek mechanism (KZM) uses critical scaling to predict density of topological defects and other excitations created in second order phase transitions. We point out that simply inserting asymptotic critical exponents deduced from the immediate vicinity of the critical point to obtain predictions can lead to results that are inconsistent with a more careful KZM analysis based on causality - on the comparison of the relaxation time of the order parameter with the "time distance" from the critical point. As a result, scaling of quench-generated excitations with quench rates can exhibit behavior that is locally (i.e., in the neighborhood of any given quench rate) well approximated by the power law, but with exponents that depend on that rate, and that are quite different from the naive prediction based on the critical exponents relevant for asymptotically long quench times. Kosterlitz-Thouless scaling (that governs e.g. Mott insulator to superfluid transition in the Bose-Hubbard model in one dimension) is investigated as an example of this phenomenon.
Kibble-Zurek mechanism (KZM) uses critical scaling to predict density of topological defects and other excitations created in second order phase transitions. We point out that simply inserting asymptotic critical exponents deduced from the immediate vicinity of the critical point to obtain KZM predictions can lead to results that are inconsistent with a more careful KZM-like analysis based on causality -- on the comparison of the relaxation time of the order parameter with the time "distance" from the critical point. As a result, scaling of quench-generated excitations with quench rates can exhibit behavior that is locally (i.e., in the neighborhood of any given quench rate) well approximated by the power law, but with exponents that depend on that rate, and that are quite different from the naive prediction based on the critical exponents relevant for asymptotically long quench times. Kosterlitz-Thouless scaling (that governs e.g. Mott insulator to superfluid transition in the Bose-Hubbard model in one dimension) is investigated as an example of this phenomenon.