We demonstrate that a nonzero concentration n v of static, randomly placed vacancies in graphene leads to a density w of zero-energy quasiparticle states at the band center ϵ ¼ 0 within a tight-binding description with nearest-neighbor hopping t on the honeycomb lattice. We show that w remains generically nonzero in the compensated case (exactly equal number of vacancies on the two sublattices) even in the presence of hopping disorder and depends sensitively on n v and correlations between vacancy positions. For low, but not-too-low, jϵj=t in this compensated case, we show that the density of states ρðϵÞ exhibits a strong divergence of the form ρ Dyson ðϵÞ ∼ jϵj −1 =½logðt=jϵjÞ ðyþ1Þ , which crosses over to the universal low-energy asymptotic form (modified Gade-Wegner scaling) expected on symmetry grounds ρ GW ðϵÞ ∼ jϵj −1 e −b½logðt=jϵjÞ 2=3 below a crossover scale ϵ c ≪ t. ϵ c is found to decrease rapidly with decreasing n v , while y decreases much more slowly. DOI: 10.1103/PhysRevLett.117.116806 Static impurities, which give rise to random timeindependent terms in the single-particle Hamiltonian for quasiparticle excitations of a condensed matter system, can lead to the phenomenon of Anderson localization, whereby quasiparticle wave functions lose their plane-wave character and become localized [1]. Such localization transitions and universal low-energy properties of the localized phase have been successfully described in many cases using effective field theories [2,3] whose form depends on symmetry properties of the quasiparticle Hamiltonian in the presence of impurities. In some cases [4,5], it has also been possible to refine these field-theoretical predictions using real-space strong-disorder renormalization group ideas [6].In this Letter, we study the effects of a nonzero concentration n v of static, randomly located vacancies in graphene. We use a tight-binding description for electronic states of graphene, with hopping amplitude t between nearest-neighbor sites on a honeycomb lattice, and model vacancies by the deletion of the corresponding site in this tight-binding model [7][8][9][10][11]. We focus on the compensated case, i.e., exactly equal numbers of vacancies on the two sublattices of the honeycomb lattice, and demonstrate that vacancies generically lead to a nonuniversal density w of zero-energy quasiparticle states at the band center ϵ ¼ 0 even in this compensated case, including in the presence of hopping disorder. For low, but not-too-low, jϵj=t in this compensated case, the density of states (DOS) ρðϵÞ exhibits a strong divergence of the formfamiliar in the context of various random-hopping problems in one dimension [12][13][14][15][16][17][18][19][20]. At still lower energies, below a crossover scale ϵ c that is several orders of magnitude smaller than t even for moderately small values of n v (0.05-0.1), we show that the DOS crosses over to the low-energy asymptotic behavior [4-6,21] of the chiral orthogonal universality class (to which our tight-binding model belongs on symmetry grou...