We introduce a semi-classical limit for many-body localization in the absence of global symmetries. Microscopically, this limit is realized by disordered Floquet circuits composed of Clifford gates. In d = 1, the resulting dynamics are always many-body localized with a complete set of strictly local integrals of motion. In d ≥ 2, the system realizes both localized and delocalized phases separated by a continuous transition in which ergodic puddles percolate. We argue that the phases are stable to deformations away from the semi-classical limit and estimate the resulting phase boundary. The Clifford circuit model is a distinct tractable limit from that of free fermions and suggests bounds on the critical exponents for the generic transition.PACS numbers: 05.30. Rt,72.15.Rn The central assumption of statistical mechanics is that interactions between particles establish local equilibrium. It is now believed that this assumption of ergodicity breaks down in isolated quantum systems in the presence of sufficient quenched disorder [1][2][3][4][5][6][7][8][9][10][11][12][13][14][15][16][17][18]. The fundamental feature of such many-body localized (MBL) systems is the slowness of the spread of entanglement [19,20] due to the presence of local integrals of motion (IOM) [8,12,[21][22][23][24]. These local observables remember initial conditions forever and thereby block both equilibration and transport. At lower disorder strengths, the system delocalizes and ergodicity is restored. There is now much theoretical and numerical evidence for this picture in one dimension [4][5][6][7][8][9][10][11][12][13][14][15][16][19][20][21][22][23][24][25][26] and the experimental search in quantum optical and atomic systems is underway [25,[27][28][29].The nature of the two phases and the dynamical transition between them is poorly understood in dimensions greater than one. In this article, we present a semiclassical theory of the delocalization transition in the absence of global symmetries in arbitrary spatial dimension. In the semi-classical limit, the transition maps precisely onto classical site percolation. Each lattice site either blocks the flow of quantum information or permits it to spread. Clusters of unblocked sites form 'ergodic puddles', as such regions locally thermalize. Blockage sites, on the other hand, host local integrals of motion. Thus, at large density of blockage sites p b , the system is localized. As p b decreases, the ergodic puddles grow and eventually percolate (p c b ) leading to delocalization. The sharp distinction between blocked and unblocked sites is a consequence of the lack of quantum tunneling in our semi-classical limit. More precisely, the localization length ξ q of the integrals of motion on blocked sites is identically zero. The delocalization transition is instead driven by the divergence of the percolation length ξ e characterizing the size of ergodic puddles. Perturbing away from the semi-classical limit, both length scales be- * Electronic address: achandran@perimeterinstitute.ca come imp...