Inhomogeneous
s-wave superconductivity is studied in the two-dimensional, square lattice attractive Hubbard Hamiltonian using the Bogoliubov-de Gennes (BdG) mean field approximation. We find that at weak coupling, and for densities mainly below half-filling, an inhomogeneous interaction in which the on-site interaction
Ui takes on two values,
Ui=0,2U results in a larger zero temperature pairing amplitude, and that the superconducting
Tc can also be significantly increased, relative to a uniform system with
Ui=U on all sites. These effects are observed for stripe, checkerboard, and even random patterns of the attractive centers, suggesting that the pattern of inhomogeneity is unimportant. Monte Carlo calculations which reintroduce some of the fluctuations neglected within the BdG approach see the same effect, both for the attractive Hubbard model and a Hamiltonian with
d-wave pairing symmetry.
We propose a highly efficient numerical method to describe inhomogeneous superconductivity by using the kernel polynomial method in order to calculate the Green's functions of a superconductor. Broken translational invariance of any type (impurities, surfaces, or magnetic fields) can be easily incorporated. We show that limitations due to system size can be easily circumvented and therefore this method opens the way for the study of scenarios and/or geometries that were unaccessible before. The proposed method is highly efficient and amenable to large scale parallel computation. Although we only use it in the context of superconductivity, it is applicable to other inhomogeneous mean-field theories.