Exploring the vortex phase diagram of Bogoliubov-de Gennes disordered superconductors

Abstract

We study the interplay of vortices and disorder in a two-dimensional disordered superconductor at zero temperature described by the Bogoliubov-de Gennes (BdG) self-consistent formalism for lattices of sizes up to 100×100 where the magnetic flux is introduced by the Peierls substitution. We model substantially larger lattice size than in previous approaches (≤36×36) which has allowed us to identify a rich phase diagram as a function of the magnetic flux and the disorder strength. For sufficiently weak disorder, and not too strong magnetic flux, we observe a slightly distorted Abrikosov triangular vortex lattice. An increase in the magnetic flux leads to an unexpected rectangular vortex lattice. A further increase in disorder, or flux, gradually destroys the lattice symmetry though strong vortex repulsion persists. An even stronger disorder leads to deformed single vortices with an inhomogeneous core. As the number of vortices increases, vortex overlap becomes more frequent. Finally, we show that global phase coherence is a feature of all these phases and that disorder enhances substantially the critical magnetic flux with respect to the clean limit with a maximum on the metallic side of the insulating transition.