We investigate the phase diagram of two-dimensional (2D) Hubbard model by employing the optimization variational Monte Carlo method. The 2D Hubbard model is the most simple electronic model for cuprate high-temperature superconductors. The phase diagram consists of three regions; they are antiferromagnetic insulator (AFI) region, superconducting (SC) region and the coexistent region of superconductivity and antiferromgantism. The phase diagram obtained by numerical calculations well agrees with the experimental phase diagram for high-temperature cuprates. We mainly focused on the effect of
t′ on the antiferromagnetic (AF) correlation and the AFI region. The area of the AF phase increases when we include
t′ and thus the pure
d-wave SC phase decreases. The AFI phase near half filling decreases as
∣t′∣ increases.
It is well known that the metal-insulator transition in two dimensions for non-interacting fermions takes place at infinitesimal disorder. In contrast, the superconductor-to-insulator transition takes place at a finite critical disorder (on the order of Vc ~ 2t), where V is the typical width of the distribution of random site energies and t is the hopping scale. In this article we compare the localization/delocalization properties of one and two particles. Whereas the metal-insulator transition is a consequence of single-particle Anderson localization, the superconductor-insulator transition (SIT) is due to pair localization – or, alternatively, fluctuations of the phase conjugate to pair density. The central question we address is how superconductivity emerges from localized single-particle states. We address this question using inhomogeneous mean field theory and quantum Monte Carlo techniques and make several testable predictions for local spectroscopic probes across the SIT. We show that with increasing disorder, the system forms superconducting blobs on the scale of the coherence length embedded in an insulating matrix. In the superconducting state, the phases on the different blobs are coherent across the system whereas in the insulator long-range phase coherence is disrupted by quantum fluctuations. As a consequence of this emergent granularity, we show that the single-particle energy gap in the density of states survives across the transition, but coherence peaks exist only in the superconductor. A characteristic pseudogap persists above the critical disorder and critical temperature, in contrast to conventional theories. Surprisingly, the insulator has a two-particle gap scale that vanishes at the SIT despite a robust single-particle gap.
We numerically study the superconductor-insulator phase transition in a model disordered two-dimensional (2D) superconductor as a function of applied magnetic field. The calculation involves quantum Monte Carlo calculations of the
(2+1)D XY model in the presence of both disorder and magnetic field. The
XY coupling is assumed to have the form $\ensuremath{-}J\text{ }\text{cos}({\ensuremath{\theta}}_{i}\ensuremath{-}{\ensuremath{\theta}}_{j}\ensuremath{-}{A}_{ij})$, where
Aij has a mean of zero and a standard deviation $\ensuremath{\Delta}{A}_{ij}$. In a real system, such a model would be approximately realized by a 2D array of small Josephson-coupled grains with slight spatial disorder and a uniform applied magnetic field. The different values $\ensuremath{\Delta}{A}_{ij}$ then correspond to an applied field such that the average number of flux quanta per plaquette has various integer values
N: larger
N corresponds to larger $\ensuremath{\Delta}{A}_{ij}$. For any value of $\ensuremath{\Delta}{A}_{ij}$, there appears to be a critical coupling constant ${K}_{c}(\ensuremath{\Delta}{A}_{ij})=\sqrt{{[J/(2U)]}_{c}}$, where
U is the charging energy, below which the system is a Mott insulator; there is also a corresponding critical conductivity ${\ensuremath{\sigma}}^{\ensuremath{\ast}}(\ensuremath{\Delta}{A}_{ij})$ at the transition. For $\ensuremath{\Delta}{A}_{ij}=\ensuremath{\infty}$, the order parameter of the transition is a renormalized coupling constant
g. Using a numerical technique appropriate for disordered systems, we show that the transition at this value of $\ensuremath{\Delta}{A}_{ij}$ takes place from an insulating (I) phase to a Bose glass (BG) phase, and that the dynamical critical exponent characterizing this transition is $z\ensuremath{\sim}1.3$. By contrast,
z=1 for this model at $\ensuremath{\Delta}{A}_{ij}=0$. We suggest that the superconductor-to-insulator transition is actually of this I to BG class at all nonzero $\ensuremath{\Delta}{A}_{ij}$'s, and we support this interpretation by both numerical evidence and an analytical argument based on the Harris criterion [A. B. Harris, J. Phys. C 7, 1671 (1974)].
Kc is found to be a monotonically increasing function of $\ensuremath{\Delta}{A}_{ij}$. For certain values of
K, a disordered Josephson array may undergo a transition from an ordered, Bose glass phase to an insulator with increasing $\ensuremath{\Delta}{A}_{ij}$.