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}$.
The scaling behavior of the current-voltage characteristics of chiral and gauge-glass models of disordered superconductors, are studied numerically, in two dimensions. For both models, the linear resistance is nonzero at finite temperatures and the scaling analysis of the nonlinear resistivity is consistent with a phase transition at
T=0 temperature characterized by a diverging correlation length $\ensuremath{\xi}\ensuremath{\propto}{T}^{\ensuremath{-}{\ensuremath{\nu}}_{T}}$ and thermal critical exponent ${\ensuremath{\nu}}_{T}.$ The values of ${\ensuremath{\nu}}_{T},$ however, are found to be different for the chiral and gauge-glass models, suggesting different universality classes, in contrast to the result obtained recently in three dimensions.
The zero-temperature universal conductivity of two-dimensional (2D) films at the supeconductor-insulator transition is studied. The existence of a finite conductivity at T=0 and the universality class for this transition is discussed. Neglecting disorder as a first approximation, so the transition is from a commensurate Mott-Hubbard insulator to a superconductor, we calculate analytically the universal conductivity for the 2D pure boson Hubbard model up to the first order in a large-N expansion and numerically by Monte Carlo simulation of the (2+1)-D XY model. From the Monte Carlo results we find the universal conductivity to be ${\mathrm{\ensuremath{\sigma}}}^{\mathrm{*}}$=(0.285\ifmmode\pm\else\textpm\fi{}0.02)${\mathrm{\ensuremath{\sigma}}}_{\mathit{Q}}$, where ${\mathrm{\ensuremath{\sigma}}}_{\mathit{Q}}^{\mathrm{\ensuremath{-}}1}$==
RQ==h/(2e
)2\ensuremath{\approxeq}6.45 k\ensuremath{\Omega}. An analysis in 1D suggests that in the presence of disorder, the universal conductivity in films might be somewhat smaller than this value. The possible existence of universal dissipation in
4He films is also discussed briefly.
In contrast to the case of free fermions which have all states localized in two dimensions, dirty, repulsively-interacting bosons exhibit a zero-temperature localization phase transition between superconducting and insulating states as disorder is increased. The conductivity σ* = g (2 e ) 2 / h at the quantum critical point is a universal critical amplitude analogous to the universal jump in superfluid density at the classical Kosterlitz-Thouless transition. Using numerical Monte Carlo methods we compute the value of g for disorder-tuned transitions with and without long-range Coulomb interactions. We argue that this model is relevant to experiments on superconducting thin films and Josephson junction arrays.
We investigate the zero-temperature superfluid to insulator transitions in a diluted two-dimensional quantum rotor model with particle-hole symmetry. We map the Hamiltonian onto a classical
(2+1)-dimensional
XY model with columnar disorder which we analyze by means of large-scale Monte Carlo simulations. For dilutions below the lattice percolation threshold, the system undergoes a generic superfluid-Mott glass transition. In contrast to other quantum phase transitions in disordered systems, its critical behavior is of conventional power-law type with universal (dilution-independent) critical exponents $z=1.52(3), \ensuremath{\nu}=1.16(5), \ensuremath{\beta}/\ensuremath{\nu}=0.48(2), \ensuremath{\gamma}/\ensuremath{\nu}=2.52(4)$, and $\ensuremath{\eta}=\ensuremath{-}0.52(4)$. These values agree with and improve upon earlier Monte Carlo results [Phys. Rev. Lett. 92, 015703 (2004)] while (partially) excluding other findings in the literature. As a further test of universality, we also consider a soft-spin version of the classical Hamiltonian. In addition, we study the percolation quantum phase transition across the lattice percolation threshold; its critical behavior is governed by the lattice percolation exponents in agreement with recent theoretical predictions. We relate our results to a general classification of phase transitions in disordered systems, and we briefly discuss experiments.
We study the vortex glass transition in disordered high temperature superconductors using Monte Carlo simulations. We use a random pinning model with strong point-correlated quenched disorder, a net applied magnetic field, long-range vortex interactions, and periodic boundary conditions. From a finite size scaling study of the helicity modulus, the rms current, and the resistivity, we obtain critical exponents at the phase transition. The new exponents differ substantially from those of the gauge glass model, but are close to those of the pure three-dimensional XY model.