Disordered graphene Josephson junctions

A tight-binding approach based on the Chebyshev-Bogoliubov-de Gennes method is used to describe disordered single-layer graphene Josephson junctions. Scattering by vacancies, ripples, or charged impurities is included. We compute the Josephson current and investigate the nature of multiple Andreev reflections, which induce bound states appearing as peaks in the density of states for energies below the superconducting gap. In the presence of single-atom vacancies, we observe a strong suppression of the supercurrent, which is a consequence of strong intervalley scattering. Although lattice deformations should not induce intervalley scattering, we find that the supercurrent is still suppressed, which is due to the presence of pseudomagnetic barriers. For charged impurities, we consider two cases depending on whether the average doping is zero, i.e., existence of electron-hole puddles, or finite. In both cases, short-range impurities strongly affect the supercurrent, similar to the vacancies scenario.