We numerically study the interplay between superconductivity and disorder on the graphene honeycomb lattice with on-site Hubbard attractive interactions U using a spatially inhomogeneous self-consistent Bogoliubov-de Gennes (BdG) approach. In the abs
ence of disorder there are two phases at charge neutrality. Below a critical value Uc for attractive interactions there is a Dirac semimetal phase and above it there is a superconducting phase. We add scalar potential disorder to the system, while remaining at charge neutrality on average. Numerical solution of the BdG equations suggests that while in the strong attraction regime (U > Uc) disorder has the usual effect of suppressing superconductivity, in the weak attraction regime (U < Uc) weak disorder enhances superconductivity. In the weak attraction regime, disorder that is too strong eventually suppresses superconductivity, i.e., there is an optimal disorder strength that maximizes the critical temperature Tc. Our numerical results also suggest that in the weakly disordered regime, mesoscopic inhomogeneities enhance superconductivity significantly more than what is predicted by a spatially uniform mean-field theory a` la Abrikosov-Gorkov. In this regime, superconductivity consists of rare phase-coherent superconducting islands. We also study the enhancement of the superconducting proximity effect by disorder and mesoscopic inhomogeneities, and obtain typical spatial plots of the tunneling density of states and the superfluid susceptibility that can be directly compared to scanning tunneling miscroscopy (STM) experiments on proximity-induced superconductivity in graphene.
