Particle pairing causes subdiffusion of heavy particles in the imbalanced Hubbard model

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.