© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

arXiv:2210.01082 · Physical review. B./Physical review. B

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

Aleksander Kaczmarek, Adam S. Sajna

Published Oct 9, 2023• 5 citations

Abstract

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.

DOI: 10.1103/physrevb.108.134304

Related papers

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

Universal semiclassical dynamics in disordered two-dimensional systems

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

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.

arXiv:2409.12956•Physical review. B./Physical review. B•Mar 20, 2025•0 citations

Systematic large flavor fTWA approach to interaction quenches in the Hubbard model

Alexander Osterkorn, Stefan K. Kehrein

Abstract We study the nonequilibrium dynamics after an interaction quench in the two-dimensional Hubbard model using the recently introduced fermionic truncated Wigner approximation (fTWA). To assess the range of validity of the method in a systematic way, we consider the SU(N) Hubbard model with the fermion degeneracy N as a natural semiclassical expansion parameter. Using both a numerical and a perturbative analytical approach we show that fTWA is exact at least up to and including the prethermalization dynamics. We discuss the limitations of the method beyond this regime.

arXiv:2007.05063•Journal of Physics A Mathematical and Theoretical•Aug 14, 2024•3 citations

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

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

Correlated quantum dynamics of graphene clusters

François Rousse, Olle Eriksson, Magnus Ögren

Phase-space representations are a family of methods for dynamics of both bosonic and fermionic systems, that work by mapping the system's density matrix to a quasiprobability density and the Liouville-von Neumann equation of the Hamiltonian to a corresponding density differential equation for the probability. We investigate here the accuracy and the computational efficiency of one approximate phase-space representation, called the fermionic truncated Wigner approximation (fTWA), applied to the Fermi-Hubbard model. On a many-body 2D system, with hopping strength and Coulomb $U$ tuned to represent the electronic structure of graphene, the method is found to be able to capture the time evolution of first-order (site occupation) and second-order (correlation functions) moments significantly better than the mean-field, Hartree-Fock method. The fTWA was also compared to results from the exact diagonalization method for smaller systems, and in general the agreement was found to be good. The fully parallel computational requirement of fTWA scales in the same order as the Hartree-Fock method, and the largest system considered here contained 198 lattice sites.

arXiv:2212.04352•Physical review. B./Physical review. B•Apr 20, 2023•2 citations

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.

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 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.

Abstract We study the nonequilibrium dynamics after an interaction quench in the two-dimensional Hubbard model using the recently introduced fermionic truncated Wigner approximation (fTWA). To assess the range of validity of the method in a systematic way, we consider the SU(N) Hubbard model with the fermion degeneracy N as a natural semiclassical expansion parameter. Using both a numerical and a perturbative analytical approach we show that fTWA is exact at least up to and including the prethermalization dynamics. We discuss the limitations of the method beyond this regime.

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.

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.

Phase-space representations are a family of methods for dynamics of both bosonic and fermionic systems, that work by mapping the system's density matrix to a quasiprobability density and the Liouville-von Neumann equation of the Hamiltonian to a corresponding density differential equation for the probability. We investigate here the accuracy and the computational efficiency of one approximate phase-space representation, called the fermionic truncated Wigner approximation (fTWA), applied to the Fermi-Hubbard model. On a many-body 2D system, with hopping strength and Coulomb

U

tuned to represent the electronic structure of graphene, the method is found to be able to capture the time evolution of first-order (site occupation) and second-order (correlation functions) moments significantly better than the mean-field, Hartree-Fock method. The fTWA was also compared to results from the exact diagonalization method for smaller systems, and in general the agreement was found to be good. The fully parallel computational requirement of fTWA scales in the same order as the Hartree-Fock method, and the largest system considered here contained 198 lattice sites.