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:2603.25659

Diffusion in interacting two-dimensional systems under a uniform magnetic field

Łukasz Iwanek, Marcin Mierzejewski, Adam S. Sajna

Published Mar 26, 2026• 0 citations
arXivHTMLPDF

Abstract

The dynamics of interacting particles in orbital magnetic fields are notoriously difficult to study, as this physics is inherently connected to electronic correlations in two-dimensional systems, for which no straightforward theoretical methods are available. Here, we report on the diffusive relaxation dynamics of two-dimensional interacting fermionic systems under a uniform magnetic field in the infinite temperature regime. We first show that the fermionic truncated Wigner approximation captures the equilibration dynamics unexpectedly well for intermediate interaction strengths when going beyond one dimension. This high accuracy holds at least for relatively small ladder systems, which are accessible to the Lanczos method that we use to benchmark the reliability of the Wigner approximation. We find that strong interactions, which exceed the hopping energy, suppress magnetic-field effects on diffusive transport. However, when the interactions are comparable to the kinetic energy, the diffusion is significantly reduced by the magnetic flux. This is observed for sufficiently large systems (above approximately 400 lattice sites), where finite-size effects weakly affect particle transport. We suggest that our results should be directly accessible on current optical lattice platforms.

DOI
10.48550/arxiv.2603.25659

Related papers

Stark Many-Body Localization in Interacting Infinite Dimensional Systems

Hristiana Atanasova, André Erpenbeck, Emanuel Gull +2 more

We study bulk particle transport in a Fermi-Hubbard model on an infinite-dimensional Bethe lattice, driven by a constant electric field. Previous numerical studies showed that one dimensional analogs of this system exhibit a breakdown of diffusion due to Stark many-body localization at least up to time that scales exponentially with the system size. Here, we consider systems initially in a spin density wave state using a combination of numerically exact and approximate techniques. We show that for sufficiently weak electric fields, the wave's momentum component decays exponentially with time in a way consistent with normal diffusion. By studying different wavelengths, we extract the dynamical exponent and the generalized diffusion coefficient at each field strength. Interestingly, we find a nonmonotonic dependence of the dynamical exponent on the electric field. As the field increases toward a critical value proportional to the Hubbard interaction strength, transport slows down, becoming subdiffusive. At large interaction strengths, however, transport speeds up again with increasing field, exhibiting superdiffusive characteristics when the electric field is comparable to the interaction strength. Eventually, at the large field limit, localization occurs and the current through the system is suppressed.

arXiv:2311.08893•Physical Review Letters•Apr 19, 2024•8 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

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

Slow dynamics in a two-dimensional Anderson-Hubbard model

Yevgeny Bar Lev, David R. Reichman

We study the real-time dynamics of a two-dimensional Anderson-Hubbard model using nonequilibrium self-consistent perturbation theory within the second-Born approximation. When compared with exact diagonalization performed on small clusters, we demonstrate that for strong disorder this technique approaches the exact result on all available timescales, while for intermediate disorder, in the vicinity of the many-body localization transition, it produces quantitatively accurate results up to nontrivial times. Our method allows for the treatment of system sizes inaccessible by any numerically exact method and for the complete elimination of finite-size effects for the times considered. We show that for a sufficiently strong disorder the system becomes nonergodic, while for intermediate disorder strengths and for all accessible timescales transport in the system is strictly subdiffusive. We argue that these results are incompatible with a simple percolation picture, but are consistent with the heuristic random resistor network model where subdiffusion may be observed for long times until a crossover to diffusion occurs. The prediction of slow finite-time dynamics in a two-dimensional interacting and disordered system can be directly verified in future cold-atoms experiments.

arXiv:1508.05391•Europhysics Letters (EPL)•Feb 1, 2016•63 citations

Nonequilibrium Green Functions Approach to Strongly Correlated Fermions in Lattice Systems

Niclas Schlünzen, M. Bönitz

Abstract Quantum dynamics in strongly correlated systems are of high current interest in many fields including dense plasmas, nuclear matter and condensed matter and ultracold atoms. An important model case are fermions in lattice systems that is well suited to analyze, in detail, a variety of electronic and magnetic properties of strongly correlated solids. Such systems have recently been reproduced with fermionic atoms in optical lattices which allow for a very accurate experimental analysis of the dynamics and of transport processes such as diffusion. The theoretical analysis of such systems far from equilibrium is very challenging since quantum and spin effects as well as correlations have to be treated non‐perturbatively. The only accurate method that has been successful so far are density matrix renormalization group (DMRG) simulations. However, these simulations are presently limited to one‐dimensional (1D) systems and short times. Extension of quantum dynamics simulations to two and three dimensions is commonly viewed as one of the major challenges in this field. Recently we have reported a breakthrough in this area [N. Schlünzen et al., Phys. Rev. B (2016)] where we were able to simulate the expansion dynamics of strongly correlated fermions in a Hubbard lattice following a quench of the confinement potential in 1D, 2D and 3D. The results not only exhibited excellent agreement with the experimental data but, in addition, revealed new features of the short‐time dynamics where correlations and entanglement are being build up. The method used in this work are nonequilibrium Green functions (NEGF) which are found to be very powerful in the treatment of fermionic lattice systems filling the gap presently left open by DMRG in 2D and 3D. In this paper we present a detailed introduction in the NEGF approach and its application to inhomogeneous Hubbard clusters. In detail we discuss the proper strong coupling approximation which is given by T ‐matrix selfenergies that sum up two‐particle scattering processes to infinite order. The efficient numerical implemen‐tation of the method is discussed in detail as it has allowed us to achieve dramatic performance gains. This has been the basis for the treatment of more than 100 particles over large time intervals. The numerical results presented in this paper concentrate on the diffusion in 1D to 3D lattices. We find that the expansion dynamics consist of three different phases that are linked with the build‐up of correlations. In the long time limit, a universal scaling with the particle number is revealed. By extrapolating the expansion velocities to the macroscopic limit, the obtained results show excellent agreement with recent experiments on ultracold fermions in optical lattices. Moreover we present results for the site‐resolved behavior of correlations and entanglement that can be directly compared with experiments using the recently developed atomic microscope technique. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)

Contributions to Plasma Physics•Jan 1, 2016•55 citations

Transport and dynamics of low-dimensional quantum spin systems

Yousef Rahnavard

In this thesis, I use a quantum Monte Carlo (QMC) method based on the stochastic series expansion (SSE) to study transport, quantum phase diagrams, and dynamics of several one- and quasi-one-dimensional spin systems, including spin-1/2 and spin-1 chains, a frustrated four-leg spin tube, and N-leg spin-1/2 ladders. I establish a diffusive channel for the finite temperature spin transport of the spin-1/2 antiferromagnetic Heisenberg chain. Extensive finite size studies are detailed, allowing to extract the diffusion kernel in the thermodynamic limit. The diffusion is found to exhibit no anomalous exponent, moreover it is shown to survive for finite values of exchange anisotropy and magnetic fields, however with a strong suppression away from the isotropic point and for increasing magnetic fields. In the case of the spin-1 chain, I obtain the quantum phase diagram and spin dynamics in the combined presence of finite single-ion anisotropies and external magnetic fields. I also prove the existence of magnon modes and multi-particle continua and determine their evolutions as a function of anisotropies and magnetic fields. Moreover, I calculate the magnetic field dependence of the nuclear magnetic relaxation as a function of the single ion anisotropy. For the quasi-one-dimensional spin system of a frustrated four-leg tube, I investigate the regime of large rung coupling, where QMC is applicable. I evaluate the uniform spin susceptibility and the dynamic structure factor and show that in this regime the spin-tube exhibts a crossover from a gapped effective antiferromagnetic spin-2 chain behavior with massive magnons to a gapless Luttinger liquid phase with spinon excitations. Finally, I investigate the spin dynamics of N-leg spin-1/2 ladders with spacially isotropic interactions for N=1,..,5. I obtain the low-lying excitations of the spin ladders as a function of the number of legs and show that they exhibt an even-odd oscillation which is consistent with Haldane's conjecture and other known results for the thermodynamics of spin ladders.

•Jul 29, 2014•0 citations
The dynamics of interacting particles in orbital magnetic fields are notoriously difficult to study, as this physics is inherently connected to electronic correlations in two-dimensional systems, for which no straightforward theoretical methods are available. Here, we report on the diffusive relaxation dynamics of two-dimensional interacting fermionic systems under a uniform magnetic field in the infinite temperature regime. We first show that the fermionic truncated Wigner approximation captures the equilibration dynamics unexpectedly well for intermediate interaction strengths when going beyond one dimension. This high accuracy holds at least for relatively small ladder systems, which are accessible to the Lanczos method that we use to benchmark the reliability of the Wigner approximation. We find that strong interactions, which exceed the hopping energy, suppress magnetic-field effects on diffusive transport. However, when the interactions are comparable to the kinetic energy, the diffusion is significantly reduced by the magnetic flux. This is observed for sufficiently large systems (above approximately 400 lattice sites), where finite-size effects weakly affect particle transport. We suggest that our results should be directly accessible on current optical lattice platforms.
We study bulk particle transport in a Fermi-Hubbard model on an infinite-dimensional Bethe lattice, driven by a constant electric field. Previous numerical studies showed that one dimensional analogs of this system exhibit a breakdown of diffusion due to Stark many-body localization at least up to time that scales exponentially with the system size. Here, we consider systems initially in a spin density wave state using a combination of numerically exact and approximate techniques. We show that for sufficiently weak electric fields, the wave's momentum component decays exponentially with time in a way consistent with normal diffusion. By studying different wavelengths, we extract the dynamical exponent and the generalized diffusion coefficient at each field strength. Interestingly, we find a nonmonotonic dependence of the dynamical exponent on the electric field. As the field increases toward a critical value proportional to the Hubbard interaction strength, transport slows down, becoming subdiffusive. At large interaction strengths, however, transport speeds up again with increasing field, exhibiting superdiffusive characteristics when the electric field is comparable to the interaction strength. Eventually, at the large field limit, localization occurs and the current through the system is suppressed.
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.
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.
We study the real-time dynamics of a two-dimensional Anderson-Hubbard model using nonequilibrium self-consistent perturbation theory within the second-Born approximation. When compared with exact diagonalization performed on small clusters, we demonstrate that for strong disorder this technique approaches the exact result on all available timescales, while for intermediate disorder, in the vicinity of the many-body localization transition, it produces quantitatively accurate results up to nontrivial times. Our method allows for the treatment of system sizes inaccessible by any numerically exact method and for the complete elimination of finite-size effects for the times considered. We show that for a sufficiently strong disorder the system becomes nonergodic, while for intermediate disorder strengths and for all accessible timescales transport in the system is strictly subdiffusive. We argue that these results are incompatible with a simple percolation picture, but are consistent with the heuristic random resistor network model where subdiffusion may be observed for long times until a crossover to diffusion occurs. The prediction of slow finite-time dynamics in a two-dimensional interacting and disordered system can be directly verified in future cold-atoms experiments.
Abstract Quantum dynamics in strongly correlated systems are of high current interest in many fields including dense plasmas, nuclear matter and condensed matter and ultracold atoms. An important model case are fermions in lattice systems that is well suited to analyze, in detail, a variety of electronic and magnetic properties of strongly correlated solids. Such systems have recently been reproduced with fermionic atoms in optical lattices which allow for a very accurate experimental analysis of the dynamics and of transport processes such as diffusion. The theoretical analysis of such systems far from equilibrium is very challenging since quantum and spin effects as well as correlations have to be treated non‐perturbatively. The only accurate method that has been successful so far are density matrix renormalization group (DMRG) simulations. However, these simulations are presently limited to one‐dimensional (1D) systems and short times. Extension of quantum dynamics simulations to two and three dimensions is commonly viewed as one of the major challenges in this field. Recently we have reported a breakthrough in this area [N. Schlünzen et al., Phys. Rev. B (2016)] where we were able to simulate the expansion dynamics of strongly correlated fermions in a Hubbard lattice following a quench of the confinement potential in 1D, 2D and 3D. The results not only exhibited excellent agreement with the experimental data but, in addition, revealed new features of the short‐time dynamics where correlations and entanglement are being build up. The method used in this work are nonequilibrium Green functions (NEGF) which are found to be very powerful in the treatment of fermionic lattice systems filling the gap presently left open by DMRG in 2D and 3D. In this paper we present a detailed introduction in the NEGF approach and its application to inhomogeneous Hubbard clusters. In detail we discuss the proper strong coupling approximation which is given by T ‐matrix selfenergies that sum up two‐particle scattering processes to infinite order. The efficient numerical implemen‐tation of the method is discussed in detail as it has allowed us to achieve dramatic performance gains. This has been the basis for the treatment of more than 100 particles over large time intervals. The numerical results presented in this paper concentrate on the diffusion in 1D to 3D lattices. We find that the expansion dynamics consist of three different phases that are linked with the build‐up of correlations. In the long time limit, a universal scaling with the particle number is revealed. By extrapolating the expansion velocities to the macroscopic limit, the obtained results show excellent agreement with recent experiments on ultracold fermions in optical lattices. Moreover we present results for the site‐resolved behavior of correlations and entanglement that can be directly compared with experiments using the recently developed atomic microscope technique. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
In this thesis, I use a quantum Monte Carlo (QMC) method based on the stochastic series expansion (SSE) to study transport, quantum phase diagrams, and dynamics of several one- and quasi-one-dimensional spin systems, including spin-1/2 and spin-1 chains, a frustrated four-leg spin tube, and N-leg spin-1/2 ladders. I establish a diffusive channel for the finite temperature spin transport of the spin-1/2 antiferromagnetic Heisenberg chain. Extensive finite size studies are detailed, allowing to extract the diffusion kernel in the thermodynamic limit. The diffusion is found to exhibit no anomalous exponent, moreover it is shown to survive for finite values of exchange anisotropy and magnetic fields, however with a strong suppression away from the isotropic point and for increasing magnetic fields. In the case of the spin-1 chain, I obtain the quantum phase diagram and spin dynamics in the combined presence of finite single-ion anisotropies and external magnetic fields. I also prove the existence of magnon modes and multi-particle continua and determine their evolutions as a function of anisotropies and magnetic fields. Moreover, I calculate the magnetic field dependence of the nuclear magnetic relaxation as a function of the single ion anisotropy. For the quasi-one-dimensional spin system of a frustrated four-leg tube, I investigate the regime of large rung coupling, where QMC is applicable. I evaluate the uniform spin susceptibility and the dynamic structure factor and show that in this regime the spin-tube exhibts a crossover from a gapped effective antiferromagnetic spin-2 chain behavior with massive magnons to a gapless Luttinger liquid phase with spinon excitations. Finally, I investigate the spin dynamics of N-leg spin-1/2 ladders with spacially isotropic interactions for N=1,..,5. I obtain the low-lying excitations of the spin ladders as a function of the number of legs and show that they exhibt an even-odd oscillation which is consistent with Haldane's conjecture and other known results for the thermodynamics of spin ladders.