  Symmetries play a central role in single-particle localization. Recent
research focused on many-body localized (MBL) systems, characterized by new
kind of integrability, and by the area-law entanglement of eigenstates. We
investigate the effect of a non-Abelian $SU(2)$ symmetry on the dynamical
properties of a disordered Heisenberg chain. While $SU(2)$ symmetry is
inconsistent with the conventional MBL, a new non-ergodic regime is possible.
In this regime, the eigenstates exhibit faster than area-law, but still a
strongly sub-thermal scaling of entanglement entropy. Using exact
diagonalization, we establish that this non-ergodic regime is indeed realized
in the strongly disordered Heisenberg chains. We use real-space renormalization
group (RSRG) to construct approximate excited eigenstates, and show their
accuracy for systems of size up to $L=26$. As disorder strength is decreased, a
crossover to the thermalizing phase occurs. To establish the ultimate fate of
the non-ergodic regime in the thermodynamic limit, we develop a novel approach
for describing many-body processes that are usually neglected by RSRG,
accessing systems of size $L>2000$. We characterize the resonances that arise
due to such processes, finding that they involve an ever growing number of
spins as the system size is increased. The probability of finding resonances
grows with the system size. Even at strong disorder, we can identify a large
lengthscale beyond which resonances proliferate. Presumably, this eventually
would drive the system to a thermalizing phase. However, the extremely long
thermalization time scales indicate that a broad non-ergodic regime will be
observable experimentally. Our study demonstrates that symmetries control
dynamical properties of disordered, many-body systems. The approach introduced
here provides a versatile tool for describing a broad range of disordered
many-body systems.
