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

Universal semiclassical dynamics in disordered two-dimensional systems

Łukasz Iwanek, Marcin Mierzejewski, Anatoli Polkovnikov, Dries Sels, Adam S. Sajna

Published Mar 20, 2025• 0 citations
arXivHTMLPDF

Abstract

The dynamics of disordered two-dimensional systems is much less understood than the dynamics of disordered chains, mainly due to the lack of appropriate numerical methods. We demonstrate that a single-trajectory version of the fermionic truncated Wigner approximation (fTWA) gives unexpectedly accurate results for the dynamics of one-dimensional (1D) systems with moderate or strong disorder. Additionally, the computational complexity of calculations carried out within this approximation is small enough to enable studies of two-dimensional (2D) systems larger than standard fTWA. Using this method, we analyze the dynamics of interacting spinless fermions propagating on disordered 1D and 2D lattices. We find for both spatial dimensions that the imbalance exhibits a universal dependence on the rescaled time $t/{\ensuremath{\xi}}_{W}$, where in two dimensions the timescale ${\ensuremath{\xi}}_{W}$ follows a stretched-exponential dependence on disorder strength.

DOI
10.1103/physrevb.111.094207

Related papers

Semiclassical bounds on the dynamics of two-dimensional interacting disordered fermions

Łukasz Iwanek, Marcin Mierzejewski, Anatoli Polkovnikov +2 more

Using the truncated Wigner approximation (TWA) we study quench dynamics of two-dimensional lattice systems consisting of interacting spinless fermions with potential disorder. First, we demonstrate that the semiclassical dynamics generally relaxes faster than the full quantum dynamics. We obtain this result by comparing the semiclassical dynamics with exact diagonalization and Lanczos propagation of one-dimensional chains. Next, exploiting the TWA capabilities of simulating large lattices, we investigate how the relaxation rates depend on the dimensionality of the studied system. We show that strongly disordered one-dimensional and two-dimensional systems exhibit a transient, logarithmic-in-time relaxation, which was recently established for one-dimensional chains. Such relaxation corresponds to the infamous $1/f$ noise at strong disorder.

arXiv:2209.15062•Physical review. B./Physical review. B•Feb 8, 2023•6 citations

Semiclassical dynamics of a disordered two-dimensional Hubbard model with long-range interactions

Adam S. Sajna, Anatoli Polkovnikov

Quench dynamics in a two-dimensional system of interacting fermions is analyzed within the semiclassical truncated Wigner approximation (TWA). The models with short-range and long-range interactions are considered. We show that in the latter case, the TWA is very accurate, becoming asymptotically exact in the infinite-range limit, provided that the semiclassical Hamiltonian is correctly identified. Within the TWA, different dynamical timescales of charges and spins can be clearly distinguished. Interestingly, for a weak and moderate disorder strength, we observe subdiffusive behavior of charges, while spins exhibit diffusive dynamics. At strong disorder, the quantum Fisher information shows logarithmic growth in time with a slower increase for charges than for spins. It is shown that in contrast to the short-range model, strong inhomogeneities such as domain walls in the initial state can significantly slow down thermalization dynamics, especially at weak disorder. This behavior can put additional challenges in designing cold-atom experimental protocols aimed to analyze possible many-body localization in such systems. While within this approach we cannot make any definite statements about the existence of a many-body localized phase, we see a very fast crossover as a function of disorder strength from rapidly thermalizing to a slow glassylike regime both for the short-range and long-range models.

arXiv:2002.05549•Physical Review A•Sep 28, 2020•13 citations

Slow semiclassical dynamics of a two-dimensional Hubbard model in disorder-free potentials

Aleksander Kaczmarek, Adam S. Sajna

The quench dynamics of the Hubbard model in tilted and harmonic potentials is discussed within the semiclassical picture. Applying the fermionic truncated Wigner approximation (fTWA), the dynamics of imbalances for charge and spin degrees of freedom is analyzed and its time evolution is compared with the exact simulations in a one-dimensional lattice. Quench from charge or spin density waves is considered. We show that introduction of harmonic and spin-dependent linear potentials sufficiently validates fTWA for longer times. Such an improvement of fTWA is also obtained for the higher-order correlations in terms of quantum Fisher information for charge and spin channels. This allows us to discuss the dynamics of larger system sizes and connect our discussion to the recently introduced Stark many-body localization. In particular, we focus on a finite two-dimensional system and show that at intermediate linear potential strength, the addition of a harmonic potential and spin dependence of the tilt results in subdiffusive dynamics, similar to that of disordered systems. Moreover, for specific values of harmonic potential, we observed phase separation of ergodic and nonergodic regions in real space. The latter fact is especially important for ultracold-atom experiments in which harmonic confinement can be easily imposed, causing a significant change in relaxation times for different lattice locations.

arXiv:2210.01082•Physical review. B./Physical review. B•Oct 9, 2023•5 citations

Dynamics of correlations in two-dimensional quantum spin models with long-range interactions: a phase-space Monte-Carlo study

Johannes Schachenmayer, Alexander Pikovski, A M Rey

Interacting quantum spin models are remarkably useful for describing different types of physical, chemical, and biological systems.Significant understanding of their equilibrium properties has been achieved to date, especially for the case of spin models with short-range couplings.However, progress toward the development of a comparable understanding in long-range interacting models, in particular out-of-equilibrium, remains limited.In a recent work, we proposed a semiclassical numerical method to study spin models, the discrete truncated Wigner approximation (DTWA), and demonstrated its capability to correctly capture the dynamics of one-and two-point correlations in one-dimensional (1D) systems.Here we go one step forward and use the DTWA method to study the dynamics of correlations in two-dimensional (2D) systems with many spins and different types of long-range couplings, in regimes where other numerical methods are generally unreliable.We compute spatial and time-dependent correlations for spin-couplings that decay with distance as a power-law and determine the velocity at which correlations propagate through the system.Sharp changes in the behavior of those velocities are found as a function of the power-law decay exponent.Our predictions are relevant for a broad range of systems including solid state materials, atom-photon systems and ultracold gases of polar molecules, trapped ions, Rydberg, and magnetic atoms.We validate the DTWA predictions for small 2D systems and 1D systems, but ultimately, in the spirt of quantum simulation, experiments will be needed to confirm our predictions for large 2D systems.

arXiv:1501.06593•New Journal of Physics•Jun 10, 2015•69 citations

Parafermionic Truncated Wigner Approximation

Javad Vahedi, Martin Gärttner

We introduce the parafermionic truncated Wigner approximation ($p$TWA), a semiclassical phase-space framework for simulating the nonequilibrium dynamics of lattice systems with fractional exchange statistics. The method extends truncated Wigner approaches developed for bosonic and fermionic systems to $\mathbb{Z}_n$ Fock parafermions by expressing the Hamiltonian in terms of local Hubbard operators that form a closed Lie algebra. This representation leads to a Lie--Poisson phase-space formulation in which quantum dynamics is approximated by stochastic sampling of initial conditions followed by deterministic semiclassical evolution. We benchmark the approach in several settings, including single-site clock dynamics, the fully connected $\mathbb{Z}_n$ clock model, long-range $\mathbb{Z}_3$ clock chains, and disordered $\mathbb{Z}_3$ Fock parafermion chains. The method reproduces key features of the exact dynamics, including excitation spreading, disorder-induced suppression of transport, and the emergence of long-time imbalance plateaus. Our results demonstrate that $p$TWA provides a practical tool for exploring the dynamics of parafermionic systems in regimes where exact numerical methods are limited by Hilbert-space growth.

arXiv:2603.29344•Mar 31, 2026•0 citations

Quench dynamics in higher-dimensional Holstein models: Insights from truncated Wigner approaches

Eva Paprotzki, Alexander Osterkorn, Vibhu Mishra +1 more

Charge-density wave phases in quantum materials stem from the complex interplay of electronic and lattice degrees of freedom. Nowadays, various time-resolved spectroscopy techniques allow to actively manipulate such phases and monitor their dynamics in real time. Modeling such nonequilibrium dynamics theoretically is a great challenge and exact methods can usually only treat a small number of atoms and finitely many phonons. We approach the melting of charge-density waves in a Holstein model after a sudden switch-on of the electronic hopping from two perspectives. We prove that in the noninteracting and in the strong-coupling limit, the CDW order parameter on high-dimensional hypercubic lattices obeys a factorization relation for long times, such that its dynamics can be reduced to the one-dimensional case. Secondly, we present numerical results from semiclassical techniques based on the truncated Wigner approximation for two spatial dimensions. A comparison with exact data obtained for a Holstein chain shows that a semiclassical treatment of both the electrons and phonons is required in order to correctly describe the phononic dynamics. This is confirmed, in addition, for a quench in the electron-phonon coupling strength.

arXiv:2312.12291•Physical review. B./Physical review. B•May 1, 2024•8 citations
The dynamics of disordered two-dimensional systems is much less understood than the dynamics of disordered chains, mainly due to the lack of appropriate numerical methods. We demonstrate that a single-trajectory version of the fermionic truncated Wigner approximation (fTWA) gives unexpectedly accurate results for the dynamics of one-dimensional (1D) systems with moderate or strong disorder. Additionally, the computational complexity of calculations carried out within this approximation is small enough to enable studies of two-dimensional (2D) systems larger than standard fTWA. Using this method, we analyze the dynamics of interacting spinless fermions propagating on disordered 1D and 2D lattices. We find for both spatial dimensions that the imbalance exhibits a universal dependence on the rescaled time t/\ensuremathξWt/{\ensuremath{\xi}}_{W}t/\ensuremathξW​, where in two dimensions the timescale \ensuremathξW{\ensuremath{\xi}}_{W}\ensuremathξW​ follows a stretched-exponential dependence on disorder strength.
Using the truncated Wigner approximation (TWA) we study quench dynamics of two-dimensional lattice systems consisting of interacting spinless fermions with potential disorder. First, we demonstrate that the semiclassical dynamics generally relaxes faster than the full quantum dynamics. We obtain this result by comparing the semiclassical dynamics with exact diagonalization and Lanczos propagation of one-dimensional chains. Next, exploiting the TWA capabilities of simulating large lattices, we investigate how the relaxation rates depend on the dimensionality of the studied system. We show that strongly disordered one-dimensional and two-dimensional systems exhibit a transient, logarithmic-in-time relaxation, which was recently established for one-dimensional chains. Such relaxation corresponds to the infamous 1/f1/f1/f noise at strong disorder.
Quench dynamics in a two-dimensional system of interacting fermions is analyzed within the semiclassical truncated Wigner approximation (TWA). The models with short-range and long-range interactions are considered. We show that in the latter case, the TWA is very accurate, becoming asymptotically exact in the infinite-range limit, provided that the semiclassical Hamiltonian is correctly identified. Within the TWA, different dynamical timescales of charges and spins can be clearly distinguished. Interestingly, for a weak and moderate disorder strength, we observe subdiffusive behavior of charges, while spins exhibit diffusive dynamics. At strong disorder, the quantum Fisher information shows logarithmic growth in time with a slower increase for charges than for spins. It is shown that in contrast to the short-range model, strong inhomogeneities such as domain walls in the initial state can significantly slow down thermalization dynamics, especially at weak disorder. This behavior can put additional challenges in designing cold-atom experimental protocols aimed to analyze possible many-body localization in such systems. While within this approach we cannot make any definite statements about the existence of a many-body localized phase, we see a very fast crossover as a function of disorder strength from rapidly thermalizing to a slow glassylike regime both for the short-range and long-range models.
The quench dynamics of the Hubbard model in tilted and harmonic potentials is discussed within the semiclassical picture. Applying the fermionic truncated Wigner approximation (fTWA), the dynamics of imbalances for charge and spin degrees of freedom is analyzed and its time evolution is compared with the exact simulations in a one-dimensional lattice. Quench from charge or spin density waves is considered. We show that introduction of harmonic and spin-dependent linear potentials sufficiently validates fTWA for longer times. Such an improvement of fTWA is also obtained for the higher-order correlations in terms of quantum Fisher information for charge and spin channels. This allows us to discuss the dynamics of larger system sizes and connect our discussion to the recently introduced Stark many-body localization. In particular, we focus on a finite two-dimensional system and show that at intermediate linear potential strength, the addition of a harmonic potential and spin dependence of the tilt results in subdiffusive dynamics, similar to that of disordered systems. Moreover, for specific values of harmonic potential, we observed phase separation of ergodic and nonergodic regions in real space. The latter fact is especially important for ultracold-atom experiments in which harmonic confinement can be easily imposed, causing a significant change in relaxation times for different lattice locations.
Interacting quantum spin models are remarkably useful for describing different types of physical, chemical, and biological systems.Significant understanding of their equilibrium properties has been achieved to date, especially for the case of spin models with short-range couplings.However, progress toward the development of a comparable understanding in long-range interacting models, in particular out-of-equilibrium, remains limited.In a recent work, we proposed a semiclassical numerical method to study spin models, the discrete truncated Wigner approximation (DTWA), and demonstrated its capability to correctly capture the dynamics of one-and two-point correlations in one-dimensional (1D) systems.Here we go one step forward and use the DTWA method to study the dynamics of correlations in two-dimensional (2D) systems with many spins and different types of long-range couplings, in regimes where other numerical methods are generally unreliable.We compute spatial and time-dependent correlations for spin-couplings that decay with distance as a power-law and determine the velocity at which correlations propagate through the system.Sharp changes in the behavior of those velocities are found as a function of the power-law decay exponent.Our predictions are relevant for a broad range of systems including solid state materials, atom-photon systems and ultracold gases of polar molecules, trapped ions, Rydberg, and magnetic atoms.We validate the DTWA predictions for small 2D systems and 1D systems, but ultimately, in the spirt of quantum simulation, experiments will be needed to confirm our predictions for large 2D systems.
We introduce the parafermionic truncated Wigner approximation (pppTWA), a semiclassical phase-space framework for simulating the nonequilibrium dynamics of lattice systems with fractional exchange statistics. The method extends truncated Wigner approaches developed for bosonic and fermionic systems to Zn\mathbb{Z}_nZn​ Fock parafermions by expressing the Hamiltonian in terms of local Hubbard operators that form a closed Lie algebra. This representation leads to a Lie--Poisson phase-space formulation in which quantum dynamics is approximated by stochastic sampling of initial conditions followed by deterministic semiclassical evolution. We benchmark the approach in several settings, including single-site clock dynamics, the fully connected Zn\mathbb{Z}_nZn​ clock model, long-range Z3\mathbb{Z}_3Z3​ clock chains, and disordered Z3\mathbb{Z}_3Z3​ Fock parafermion chains. The method reproduces key features of the exact dynamics, including excitation spreading, disorder-induced suppression of transport, and the emergence of long-time imbalance plateaus. Our results demonstrate that pppTWA provides a practical tool for exploring the dynamics of parafermionic systems in regimes where exact numerical methods are limited by Hilbert-space growth.
Charge-density wave phases in quantum materials stem from the complex interplay of electronic and lattice degrees of freedom. Nowadays, various time-resolved spectroscopy techniques allow to actively manipulate such phases and monitor their dynamics in real time. Modeling such nonequilibrium dynamics theoretically is a great challenge and exact methods can usually only treat a small number of atoms and finitely many phonons. We approach the melting of charge-density waves in a Holstein model after a sudden switch-on of the electronic hopping from two perspectives. We prove that in the noninteracting and in the strong-coupling limit, the CDW order parameter on high-dimensional hypercubic lattices obeys a factorization relation for long times, such that its dynamics can be reduced to the one-dimensional case. Secondly, we present numerical results from semiclassical techniques based on the truncated Wigner approximation for two spatial dimensions. A comparison with exact data obtained for a Holstein chain shows that a semiclassical treatment of both the electrons and phonons is required in order to correctly describe the phononic dynamics. This is confirmed, in addition, for a quench in the electron-phonon coupling strength.