The imbalanced Hubbard model features a transition between dynamic regimes depending on the mass ratio and coupling strength between two different particle species. A slowdown of the lighter particle transport can be attributed to an emergent effective disorder induced by the heavy particles for high mass ratio and strong coupling. This subdiffusive regime has been interpreted as a Griffiths phase, linking the effect to the coexistence of metallic and insulating regions. Here, we investigate the dynamics of the heavy particles, which also reveals subdiffusive behavior, yet cannot be explained within the Griffiths picture. We demonstrate that heavy particles predominantly form small clusters, mainly pairs, during the dynamical process, which reduces their propagation speed and transiently shifts the time-dependent diffusion constant into the subdiffusive regime at late times. The necessary attraction between particles driving this process can be understood within the Born-Oppenheimer approximation. We introduce a classical one-dimensional random walk model that can quantitatively reproduce the subdiffusion dynamics in the strong coupling regime.
The imbalanced Hubbard model features a transition between dynamic regimes depending on the mass ratio and coupling strength between two different particle species. A slowdown of the lighter particle transport can be attributed to an emergent effective disorder induced by the heavy particles for high mass ratio and strong coupling. This subdiffusive regime has been interpreted as a Griffiths phase, linking the effect to the coexistence of metallic and insulating regions. Here, we investigate the dynamics of the heavy particles, which also reveals subdiffusive behavior, yet cannot be explained within the Griffiths picture. We demonstrate that heavy particles predominantly form small clusters, mainly pairs, during the dynamical process, which reduces their propagation speed and transiently shifts the time-dependent diffusion constant into the subdiffusive regime at late times. The necessary attraction between particles driving this process can be understood within the Born-Oppenheimer approximation. We introduce a classical one-dimensional random walk model that can quantitatively reproduce the subdiffusion dynamics in the strong coupling regime.
We consider two mutually interacting fermionic particle species on a one-dimensional lattice and study how the mass ratio $\ensuremath{\eta}$ between the two species affects the (equilibration) dynamics of the particles. Focusing on the regime of strong interactions and high temperatures, two well-studied points of reference are given by (i) the case of equal masses $\ensuremath{\eta}=1$, i.e., the standard Fermi-Hubbard chain, where initial nonequilibrium density distributions are known to decay, and (ii) the case of one particle species being infinitely heavy, $\ensuremath{\eta}=0$, leading to a localization of the lighter particles in an effective disorder potential. Given these two opposing cases, the dynamics in the case of intermediate mass ratios $0<\ensuremath{\eta}<1$ is of particular interest. To this end, we study the real-time dynamics of pure states featuring a sharp initial nonequilibrium density profile. Relying on the concept of dynamical quantum typicality, the resulting nonequilibrium dynamics can be related to equilibrium correlation functions. Summarizing our main results, we observe that diffusive transport occurs for moderate values of the mass imbalance and manifests itself in a Gaussian spreading of real-space density profiles and an exponential decay of density modes in momentum space. For stronger imbalances, we provide evidence that transport becomes anomalous on intermediate timescales, and in particular, our results are consistent with the absence of strict localization in the long-time limit for any $\ensuremath{\eta}>0$. Based on our numerical analysis, we provide an estimate for the ``lifetime'' of the effective localization as a function of $\ensuremath{\eta}$.
We study transport in a one-dimensional boundary-driven Anderson insulator (the XX spin chain with onsite disorder) with randomly positioned onsite dephasing, observing a transition from diffusive to subdiffusive spin transport below a critical density of sites with dephasing. This model is intended to mimic the passage of an excitation through (many-body) insulating regions or ergodic bubbles, therefore providing a toy model for the diffusion-subdiffusion transition observed in the disordered Heisenberg model by \ifmmode \check{Z}\else \v{Z}\fi{}nidari\ifmmode \check{c}\else \v{c}\fi{} et al. [Phys. Rev. Lett. 117, 040601 (2016)]. We also present the exact solution of a semiclassical model of conductors and insulators introduced by Agarwal et al. [Phys. Rev. Lett. 114, 160401 (2015)], which exhibits both diffusive and subdiffusive phases, and qualitatively reproduces the results of the quantum system. The critical properties of both models, when passing from diffusion to subdiffusion, are interpreted in terms of ``Griffiths effects.'' We show that the finite-size scaling comes from the interplay of three characteristic lengths: one associated with disorder (the localization length), one with dephasing, and the third with the percolation problem defining large, rare, insulating regions. We conjecture that the latter, which grows logarithmically with system size, may potentially be responsible for the fact that heavy-tailed resistance distributions typical of Griffiths effects have not been observed in subdiffusive interacting systems.
We study transport in the boundary-driven XX spin chain with onsite disorder and randomly positioned onsite dephasing, observing a transition from diffusive to subdiffusive spin transport below a critical density of sites with dephasing. We then present an exactly solvable semiclassical model of conductors and insulators, which exhibits both diffusive and subdiffusive phases, and qualitatively reproduces the results of the quantum system. The subdiffusion in these models is a consequence of rare insulating regions, and they therefore exhibit, in a more controlled setup, the physics of which have been conjectured to cause the subdiffusive transport observed in interacting many-body localizable systems. For the quantum model we show that finite-size effects come from the interplay of three characteristic lengths: one associated with disorder (the localization length), one with dephasing, and the third with the percolation problem defining large, rare, insulating regions. The latter grows logarithmically with system size, and we conjecture that this may be the reason why the heavy-tailed distributions typical of Griffiths effects have not been observed in subdiffusive interacting systems.
We develop an analytical theory of anomalous transport in a noninteracting Anderson model with heavy-tailed hopping amplitudes. The broad distribution of hopping amplitudes gives rise to an extended intermediate-time regime with a subdiffusive effective exponent, despite the absence of interactions or genuine many-body effects. By solving the transport equations analytically, we derive the time dependence of the mean-square displacement and identify a continuous crossover from an intermediate anomalous regime to asymptotically diffusive transport. As the localization transition is approached, the spatial extent of the subdiffusive window diverges parametrically, while the crossover to conventional diffusion remains finite in units of Γ0−1. This produces an increasingly broad anomalous transport regime in space that can closely resemble Griffiths-type transport observed near the many-body localization transition. Our results demonstrate that rare hopping processes alone provide a microscopic single-particle mechanism for robust transport anomalies, establishing an analytical benchmark for distinguishing interaction-induced effects from disorder-driven dynamics.
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.
The effects of mass imbalance in a three-component Hubbard model are studied by the dynamical mean-field theory combined with exact diagonalization. The model describes a fermion-fermion mixture of two different particle species with a mass imbalance. One species is two-component fermion particles, and the other is single-component ones. The local interaction between particle species is considered isotropically. It is found that the mass imbalance can drive the mixture from insulator to metal. The insulator-metal transition is a species-selective-like transition of lighter mass particles and occurs only at commensurate particle densities and moderate local interactions. For weak and strong local interactions the mass imbalance does not change the ground state of the mixture.