Anomalous transport of small polarons arises from transient lattice relaxation or immovable boundaries

Abstract: Elucidating transport mechanisms and predicting transport coefficients is crucial for advancing material innovation, design, and application. Yet, state-of-the-art calculations are restricted to exact simulations of small lattices with severe finite-size effects or approximate simulations that assume the nature of transport. We leverage recent algorithmic advances to perform exact simulations of the celebrated Holstein model that systematically quantify and eliminate finite-size effects to gain insights into small polaron formation and the nature and timescales of its transport. We perform the first systematic comparison of the performance of two distinct approaches to predict charge carrier mobility: equilibrium-based Green-Kubo relations and nonequilibrium relaxation methods. Our investigation uncovers when and why disparities arise between these ubiquitously used techniques, revealing that the equilibrium-based method is highly sensitive to system topology whereas the nonequilibrium approach requires bigger system sizes to reveal its diffusive region. Contrary to assumptions made in standard perturbative calculations, our results demonstrate that small polarons exhibit anomalous transport and that it manifests transiently, due to nonequilibrium lattice relaxation, or permanently, as a signature of immovable boundaries. These findings have consequences for applications including the utilization of organic polymers in organic electronics and transition metal oxides in photocatalysis.