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:2502.04927 · Chinese Physics B

Scaling corrections in driven critical dynamics: Application to the two-dimensional dimerized quantum Heisenberg model

Jing-Wen 静雯 Liu 刘, Shuai Yin, Yu-Rong Shu

Published Mar 28, 2025• 1 citation
arXivHTMLPDF

Abstract

Abstract Driven critical dynamics in quantum phase transitions holds significant theoretical importance, and also has practical applications in fast-developing quantum devices. While scaling corrections have been shown to play important roles in fully characterizing equilibrium quantum criticality, their impact on nonequilibrium critical dynamics has not been extensively explored. In this work, we investigate the driven critical dynamics in a two-dimensional quantum Heisenberg model. We find that in this model the scaling corrections arising from both finite system size and finite driving rate must be incorporated into the finite-time scaling form in order to properly describe the nonequilibrium scaling behaviors. In addition, improved scaling relations are obtained from the expansion of the full scaling form. We numerically verify these scaling forms and improved scaling relations for different starting states using the nonequilibrium quantum Monte Carlo algorithm.

DOI
10.1088/1674-1056/adc672

Related papers

Imaginary-time relaxation quantum critical dynamics in the two-dimensional dimerized Heisenberg model

Jiaqi Cai, Yu-Rong Shu, Xue-Qing Rao +1 more

We study the imaginary-time relaxation critical dynamics of the N\'eel paramagnetic quantum phase transition in the two-dimensional dimerized $S=\frac{1}{2}$ Heisenberg model. We focus on the scaling correction in the short-time region. A unified scaling form including both short-time and finite-size corrections is proposed. According to this full scaling form, improved short-imaginary-time scaling relations are obtained. We numerically verify the scaling form and the improved short-time scaling relations for different initial states using projector quantum Monte Carlo algorithm.

arXiv:2403.09084•Physical review. B./Physical review. B•May 6, 2024•3 citations

Universal Driven Critical Dynamics of Entanglement Entropy

Chang-Yu Shen, Shuai Yin, Zi-Xiang Li

The Kibble-Zurek mechanism (KZM) and finite-time scaling (FTS) provide a foundational framework for driven critical dynamics, yet their predictive power has been largely confined to local observables. Here, we establish a universal finite-time scaling theory for the nonequilibrium dynamics of quantum entanglement. Using unbiased quantum Monte Carlo simulations, we investigate the corner entanglement entropy of (2+1)-dimensional interacting Dirac fermions driven from ordered phases toward a quantum critical point. We find that the corner entanglement accurately obeys a universal driven scaling governed by the driving rate and system size, persisting whether the initial ordered state is fully gapped or hosts gapless Goldstone modes. Crucially, this dynamical entanglement exhibits a logarithmic dependence on the driving rate, from which the universal corner coefficient of the underlying conformal field theory can be robustly extracted far from equilibrium. These results generalize the KZM from local observables to the intrinsic nonlocal quantum information measures, offering a practical blueprint for characterizing quantum criticality and entanglement on programmable quantum simulators.

cond-mat.str-elCross-listedarXiv:2609.03854•Sep 4, 2026•0 citations

Nonequilibrium phase transition in a driven-dissipative quantum antiferromagnet

Mona H. Kalthoff, Dante Marvin Kennes, Andrew J. Millis +1 more

A deeper theoretical understanding of driven-dissipative interacting systems and their nonequilibrium phase transitions is essential both to advance our fundamental physics understanding and to harness technological opportunities arising from optically controlled quantum many-body states. This paper provides a numerical study of dynamical phases and the transitions between them in the nonequilibrium steady state of the prototypical two-dimensional Heisenberg antiferromagnet with drive and dissipation. We demonstrate a nonthermal transition that is characterized by a qualitative change in the magnon distribution from subthermal at low drive to a generalized Bose-Einstein form including a nonvanishing condensate fraction at high drive. A finite-size analysis reveals static and dynamical critical scaling at the transition, with a discontinuous slope of the magnon number versus driving field strength and critical slowing down at the transition point. Implications for experiments on quantum materials and polariton condensates are discussed.

arXiv:2107.03841•Physical Review Research•May 12, 2022•9 citations

Dynamical critical behavior of the three-dimensional Heisenberg model

Pawel Peczak, David Paul Landau

We have carried out extensive Monte Carlo simulations of the dynamic critical behavior of L\ifmmode\times\else\texttimes\fi{}L\ifmmode\times\else\texttimes\fi{}L simple-cubic classical Heisenberg models with periodic boundary conditions. The equilibrium relaxation time \ensuremath{\tau} is determined from the long-time exponential decay of the time-displaced correlation function of the magnetization. Using finite-size scaling and corrections to scaling we find that at the critical point ${\mathit{T}}_{\mathit{c}}$=(0.6929${)}^{\mathrm{\ensuremath{-}}1}$, \ensuremath{\tau} grows like \ensuremath{\tau}\ensuremath{\simeq}${\mathit{L}}^{\mathit{z}}$ with z=1.96(6). In the paramagnetic region we show that the relaxation time is represented by a universal scaling function \ensuremath{\tau}${\mathit{L}}^{\mathrm{\ensuremath{-}}\mathit{z}}$ vs (T/${\mathit{T}}_{\mathit{c}}$-1)${\mathit{L}}^{1/\ensuremath{\nu}}$ with z=1.96(6). This estimate is consistent with theoretical predictions and results of experimental studies concerning the dynamical critical behavior of Heisenberg magnets perturbed by order-parameter nonconserving dipolar interactions. We also describe how a simple over-relaxation algorithm of Creutz and of Brown and Woch can be incorporated into the Metropolis sampling method for the classical Heisenberg model, thereby significantly reducing the critical slowing down; the resultant correlation time varies as \ensuremath{\tau}\ensuremath{\sim}(T/${\mathit{T}}_{\mathit{c}}$-1${)}^{\mathrm{\ensuremath{-}}0.77}$, which implies z\ensuremath{\simeq}1.10.

Physical review. B, Condensed matter•Jun 1, 1993•45 citations

Evidence for an Unconventional Universality Class from a Two-Dimensional Dimerized Quantum Heisenberg Model

S. C. Wenzel, L. Bogacz, Wolfhard Janke E

The two-dimensional $J\mathrm{\text{\ensuremath{-}}}{J}^{\ensuremath{'}}$ dimerized quantum Heisenberg model is studied on the square lattice by means of (stochastic series expansion) quantum Monte Carlo simulations as a function of the coupling ratio $\ensuremath{\alpha}={J}^{\ensuremath{'}}/J$. The critical point of the order-disorder quantum phase transition in the $J\mathrm{\text{\ensuremath{-}}}{J}^{\ensuremath{'}}$ model is determined as ${\ensuremath{\alpha}}_{c}=2.5196(2)$ by finite-size scaling for up to approximately 10 000 quantum spins. By comparing six dimerized models we show, contrary to the current belief, that the critical exponents of the $J\mathrm{\text{\ensuremath{-}}}{J}^{\ensuremath{'}}$ model are not in agreement with the three-dimensional classical Heisenberg universality class. This lends support to the notion of nontrivial critical excitations at the quantum critical point.

arXiv:0805.2500•Physical Review Letters•Sep 19, 2008•88 citations

Quantum phase transitions of 2-d dimerized spin-1/2 Heisenberg models with spatial anisotropy

Ming-Tso Kao, D.-R. Tan, Fu-Jiun Jiang

Motivated by the unexpected Monte Carlo results as well as the theoretical proposal of a large correction to scaling for the critical theory of the 2-d staggered-dimer spin-1/2 Heisenberg model on the square lattice, we study the phase transitions induced by dimerization of several dimerized quantum Heisenberg models with spatial anisotropy using first principles Monte Carlo method. Remarkably, while our Monte Carlo data for all the models considered here, including the herringbone- and ladder-dimer models on the square lattice, are compatible with the recently proposed scenario of an enhanced correction to scaling, we find it is likely that the enhanced correction to scaling manifests itself as amplification of the nonuniversal prefactors appeared in the scaling forms. In other words, our data are in consistence with the established numerical values for the critical exponents, including the confluent exponent, in the O(3) universality class. Convincing numerical evidence is provided to support this proposed scenario.

arXiv:1202.1057•Bulletin of the American Physical Society•Feb 6, 2012•0 citations
Abstract Driven critical dynamics in quantum phase transitions holds significant theoretical importance, and also has practical applications in fast-developing quantum devices. While scaling corrections have been shown to play important roles in fully characterizing equilibrium quantum criticality, their impact on nonequilibrium critical dynamics has not been extensively explored. In this work, we investigate the driven critical dynamics in a two-dimensional quantum Heisenberg model. We find that in this model the scaling corrections arising from both finite system size and finite driving rate must be incorporated into the finite-time scaling form in order to properly describe the nonequilibrium scaling behaviors. In addition, improved scaling relations are obtained from the expansion of the full scaling form. We numerically verify these scaling forms and improved scaling relations for different starting states using the nonequilibrium quantum Monte Carlo algorithm.
We study the imaginary-time relaxation critical dynamics of the N\'eel paramagnetic quantum phase transition in the two-dimensional dimerized S=12S=\frac{1}{2}S=21​ Heisenberg model. We focus on the scaling correction in the short-time region. A unified scaling form including both short-time and finite-size corrections is proposed. According to this full scaling form, improved short-imaginary-time scaling relations are obtained. We numerically verify the scaling form and the improved short-time scaling relations for different initial states using projector quantum Monte Carlo algorithm.
The Kibble-Zurek mechanism (KZM) and finite-time scaling (FTS) provide a foundational framework for driven critical dynamics, yet their predictive power has been largely confined to local observables. Here, we establish a universal finite-time scaling theory for the nonequilibrium dynamics of quantum entanglement. Using unbiased quantum Monte Carlo simulations, we investigate the corner entanglement entropy of (2+1)-dimensional interacting Dirac fermions driven from ordered phases toward a quantum critical point. We find that the corner entanglement accurately obeys a universal driven scaling governed by the driving rate and system size, persisting whether the initial ordered state is fully gapped or hosts gapless Goldstone modes. Crucially, this dynamical entanglement exhibits a logarithmic dependence on the driving rate, from which the universal corner coefficient of the underlying conformal field theory can be robustly extracted far from equilibrium. These results generalize the KZM from local observables to the intrinsic nonlocal quantum information measures, offering a practical blueprint for characterizing quantum criticality and entanglement on programmable quantum simulators.
A deeper theoretical understanding of driven-dissipative interacting systems and their nonequilibrium phase transitions is essential both to advance our fundamental physics understanding and to harness technological opportunities arising from optically controlled quantum many-body states. This paper provides a numerical study of dynamical phases and the transitions between them in the nonequilibrium steady state of the prototypical two-dimensional Heisenberg antiferromagnet with drive and dissipation. We demonstrate a nonthermal transition that is characterized by a qualitative change in the magnon distribution from subthermal at low drive to a generalized Bose-Einstein form including a nonvanishing condensate fraction at high drive. A finite-size analysis reveals static and dynamical critical scaling at the transition, with a discontinuous slope of the magnon number versus driving field strength and critical slowing down at the transition point. Implications for experiments on quantum materials and polariton condensates are discussed.
We have carried out extensive Monte Carlo simulations of the dynamic critical behavior of L\ifmmode\times\else\texttimes\fi{}L\ifmmode\times\else\texttimes\fi{}L simple-cubic classical Heisenberg models with periodic boundary conditions. The equilibrium relaxation time \ensuremath{\tau} is determined from the long-time exponential decay of the time-displaced correlation function of the magnetization. Using finite-size scaling and corrections to scaling we find that at the critical point Tc{\mathit{T}}_{\mathit{c}}Tc​=(0.6929${)}^{\mathrm{\ensuremath{-}}1}$, \ensuremath{\tau} grows like \ensuremath{\tau}\ensuremath{\simeq}Lz{\mathit{L}}^{\mathit{z}}Lz with z=1.96(6). In the paramagnetic region we show that the relaxation time is represented by a universal scaling function \ensuremath{\tau}${\mathit{L}}^{\mathrm{\ensuremath{-}}\mathit{z}}$ vs (T/Tc{\mathit{T}}_{\mathit{c}}Tc​-1)${\mathit{L}}^{1/\ensuremath{\nu}}$ with z=1.96(6). This estimate is consistent with theoretical predictions and results of experimental studies concerning the dynamical critical behavior of Heisenberg magnets perturbed by order-parameter nonconserving dipolar interactions. We also describe how a simple over-relaxation algorithm of Creutz and of Brown and Woch can be incorporated into the Metropolis sampling method for the classical Heisenberg model, thereby significantly reducing the critical slowing down; the resultant correlation time varies as \ensuremath{\tau}\ensuremath{\sim}(T/Tc{\mathit{T}}_{\mathit{c}}Tc​-1${)}^{\mathrm{\ensuremath{-}}0.77}$, which implies z\ensuremath{\simeq}1.10.
The two-dimensional $J\mathrm{\text{\ensuremath{-}}}{J}^{\ensuremath{'}}$ dimerized quantum Heisenberg model is studied on the square lattice by means of (stochastic series expansion) quantum Monte Carlo simulations as a function of the coupling ratio $\ensuremath{\alpha}={J}^{\ensuremath{'}}/J$. The critical point of the order-disorder quantum phase transition in the $J\mathrm{\text{\ensuremath{-}}}{J}^{\ensuremath{'}}$ model is determined as ${\ensuremath{\alpha}}_{c}=2.5196(2)$ by finite-size scaling for up to approximately 10 000 quantum spins. By comparing six dimerized models we show, contrary to the current belief, that the critical exponents of the $J\mathrm{\text{\ensuremath{-}}}{J}^{\ensuremath{'}}$ model are not in agreement with the three-dimensional classical Heisenberg universality class. This lends support to the notion of nontrivial critical excitations at the quantum critical point.
Motivated by the unexpected Monte Carlo results as well as the theoretical proposal of a large correction to scaling for the critical theory of the 2-d staggered-dimer spin-1/2 Heisenberg model on the square lattice, we study the phase transitions induced by dimerization of several dimerized quantum Heisenberg models with spatial anisotropy using first principles Monte Carlo method. Remarkably, while our Monte Carlo data for all the models considered here, including the herringbone- and ladder-dimer models on the square lattice, are compatible with the recently proposed scenario of an enhanced correction to scaling, we find it is likely that the enhanced correction to scaling manifests itself as amplification of the nonuniversal prefactors appeared in the scaling forms. In other words, our data are in consistence with the established numerical values for the critical exponents, including the confluent exponent, in the O(3) universality class. Convincing numerical evidence is provided to support this proposed scenario.