Many-body localization (MBL) is a phase of matter that is characterized by the absence of thermalization. Dynamical generation of a large number of local quantum numbers has been identified as one key characteristic of this phase, quite possibly the microscopic mechanism of breakdown of thermalization and the phase transition itself. We formulate a robust algorithm, based on Wegner-Wilson flow (WWF) renormalization, for computing these conserved quantities and their interactions. We present evidence for the existence of distinct fixed point distributions of the latter: a Gaussian white-noise-like distribution in the ergodic phase, a 1=f law inside the MBL phase, and scale-free distributions in the transition regime. DOI: 10.1103/PhysRevLett.119.075701 Recent progress on the theory of many-body localization (MBL) demonstrates clearly that the conventional quantum statistical description of interacting many-body problems is incomplete. Concrete analytic [1], numerical [2][3][4][5], and mathematical [6,7] results establish the existence and robustness of many-body localized phases in sufficiently strongly disordered and/or low-dimensional interacting models at finite extensive entropy. While the understanding of the transition between thermal and MBL phases is only beginning to emerge [8][9][10][11][12], several distinct new directions of inquiry related to MBL and the fundamental issue of ergodicity in quantum many-body systems have taken shape. These include the interplay of MBL with spontaneous symmetry breaking and topological order [13][14][15][16], selflocalization (glassiness) in translationally invariant quantum systems [17][18][19][20], and MBL in driven systems [21][22][23]. MBL has also stimulated considerable progress in developing tools for describing excited eigenstates of many-body systems [12,[24][25][26][27][28]. MBL has been realized in recent experiments [29,30] and may also have important implications for quantum engineering problems, e.g., quantum computing [31][32][33][34][35].One natural route to the breakdown of thermalization is via proliferation of a large number of conserved quasilocal quantities. The extreme version of such a proposal has gained considerable traction as a model phenomenology [36] of the so-called fully MBL regime, where the entire many-body spectrum is localized. Consider a generic system, e.g., the n-site spin-1=2 random-field Heisenberg chain [see Eq. (7)], which is diagonalized by a (nonunique) unitary matrix U.