Diamagnetic Response of Disordered Graphene to Nonuniform Magnetic Field

The diamagnetic susceptibility is calculated in disordered monolayer graphene in spatially varying magnetic field within a self-consistent Born approximation. The universality that the susceptibility is a function of q/2k(F) in ideal graphene is nearly valid in the weak-disorder case except that the susceptibility becomes nonzero even for slowly-varying field. In the case of strong disorder the universal behavior is completely destroyed. The mirror-image-like response to a magnet is still valid at the Dirac point as long as the disorder is sufficiently weak.