Apparent slow dynamics in the ergodic phase of a driven many-body localized system without extensive conserved quantities

We numerically study the dynamics on the ergodic side of the many-body localization transition in a periodically driven Floquet model with no global conservation laws. We describe and employ a numerical technique based on the fast Walsh-Hadamard transform that allows us to perform an exact time evolution for large systems and long times. As in models with conserved quantities (e.g., energy and/or particle number) we observe a slowing down of the dynamics as the transition into the many-body localized phase is approached. More specifically, our data are consistent with a subballistic spread of entanglement and a stretched-exponential decay of an autocorrelation function, with their associated exponents reflecting slow dynamics near the transition for a fixed system size. However, with access to larger system sizes, we observe a clear flow of the exponents towards faster dynamics and cannot rule out that the slow dynamics is a finite-size effect. Furthermore, we observe examples of nonmonotonic dependence of the exponents with time, with the dynamics initially slowing down but accelerating again at even larger times, consistent with the slow dynamics being a crossover phenomenon with a localized critical point.