Eigenstates hybridize on all length scales at the many-body localization transition

An interacting quantum system can transition from an ergodic to a many-body localized (MBL) phase under the presence of sufficiently large disorder. Both phases are radically different in their dynamical properties, which are characterized by highly excited eigenstates of the Hamiltonian. Each eigenstate can be characterized by the set of quantum numbers over the set of (local, in the MBL phase) integrals of motion of the system. In this work we study the evolution of the eigenstates of the disordered Heisenberg model as the disorder strength, $W$, is varied adiabatically. We focus on the probability that two `colliding' eigenstates hybridize as a function of both the range $R$ at which they differ as well as the strength of their hybridization. We find, in the MBL phase, that the probability of a colliding eigenstate hybridizing strongly at range $R$ decays as $Pr(R)\propto \exp [-R/\eta]$, with a length scale $\eta(W) = 1 / (B \log(W / W_c) )$ which diverges at the critical disorder strength $W_c$. This leads to range-invariance at the transition, suggesting the formation of resonating cat states at all ranges. This range invariance does not survive to the ergodic phase, where hybridization is exponentially more likely at large range, a fact that can be understood with simple combinatorial arguments. In fact, compensating for these combinatorial effects allows us to define an additional correlation length $\xi$ in the MBL phase which is in excellent agreement with previous works and which takes the critical value $1 / \log(2)$ at the transition, found in previous works to destabilize the MBL phase. Finally, we show that deep in the MBL phase hybridization is dominated by two-level collisions of eigenstates close in energy.