Topological delocalization of two-dimensional massless dirac fermions

The beta function of a two-dimensional massless Dirac Hamiltonian subject to a random scalar potential, which, e.g., underlies theoretical descriptions of graphene, is computed numerically. Although it belongs to, from a symmetry standpoint, the two-dimensional symplectic class, the beta function monotonically increases with decreasing conductance. We also provide an argument based on the spectral flows under twisting boundary conditions, which shows that none of the states of the massless Dirac Hamiltonian can be localized.