Inhomogeneous pairing in highly disordered s-wave superconductors

We study a simple model of a two-dimensional s-wave superconductor in the presence of a random potential as a function of disorder strength. We first use the Bogoliubov–de Gennes (BdG) approach to show that, with increasing disorder, the pairing amplitude becomes spatially inhomogeneous, and the system cannot be described within conventional approaches for studying disordered superconductors that assume a uniform order parameter. In the high-disorder regime, we find that the system breaks up into superconducting islands, with large pairing amplitude, separated by an insulating sea. We show that this inhomogeneity has important implications for the physical properties of this system, such as superfluid density and the density of states. We find that a finite spectral gap persists in the density of states, even in the weak-coupling regime, for all values of disorder, and we provide a detailed understanding of this remarkable result. We next generalize Anderson’s idea of the pairing of exact eigenstates to include an inhomogeneous pairing amplitude, and show that it is able to qualitatively capture many of the nontrivial features of the full BdG analysis. Finally, we study the transition to a gapped insulating state driven by quantum phase fluctuations about the inhomogeneous superconducting state.