# Quantum superconductor-insulator transition: Implications of BKT-critical behavior

## Abstract

We explore the implications of Berezinskii-Kosterlitz-Thouless (BKT) critical behavior on the two dimensional (2D) quantum superconductor-insulator (QSI) transition driven by the tuning parameter . Concentrating on the sheet resistance BKT behavior implies: an explicit quantum scaling function for along the superconducting branch ending at the nonuniversal critical value ; a BKT-transition line where is the dynamic and the exponent of the zero temperature correlation length; independent estimates of , and from the dependence of the nonuniversal parameters entering the BKT expression for the sheet resistance. To illustrate the potential and the implications of this scenario we analyze data of Bollinger et al. taken on gate voltage tuned epitaxial films of LaSrCuO that are one unit cell thick. The resulting estimates and point to a 2D-QSI critical point where hyperscaling, the proportionality between and and the correspondence between quantum phase transitions in D and classical ones in (D+z) dimensions are violated and disorder is relevant.

###### pacs:

74.40.Kb, 74.62.-c, 74.78.-w^{1}

## I Introduction

A variety of different materials undergo a quantum superconductor- insulator (QSI) transition in the limit of two dimensions (2D) and zero temperature by variation of a tuning parameter including film thickness, disorder, applied magnetic field, and gate voltage.(1); (2); (3); (4) A widespread observable to study this behavior is the temperature dependence of the sheet resistance taken at various values of the tuning parameter . The curves of ) at different resemble the flow to a nearly temperature independent separatrix between superconducting and insulating phase with sheet resistance . This behavior implies a crossing point of the isotherms at different temperatures at which is a characteristic feature of a quantum phase transition and in the present case of a QSI transition. Traditionally the interpretation of experimental data taken close to the 2D-QSI transition is based on the quantum scaling relation with .(5); (6); (7); (8) is the dynamic, the critical exponent of the zero temperature correlation length and a nonuniversal coefficient of proportionality.(1); (2); (3) Given sheet resistance data fits to this scaling form yield estimates for the critical value of the tuning parameter and the exponent product , properties which are insufficient to distinguish different models, to fix the universality class to which the QSI transition belongs, or to clarify the relevance of disorder .

The nature of the 2D-QSI transition has been intensely debated.(1); (2); (3); (4) The scenarios can be grouped into two classes,
fermionic and bosonic. In the fermionic case the reduction of and
the magnitude of the order parameter is attributed to a combination of
reduced density of states, enhanced Coulomb interaction and depairing due to
an increase of the inelastic electron-electron scattering rate.(9); (10) The bosonic approach assumes that the fermionic degrees of
freedom can be integrated out, the mean square of the order parameter does
not vanish at , phase fluctuations dominate and the reduction of is attributable to quantum fluctuations and in disordered systems to
randomness in addition.(6); (7); (11) This scenario is closely
related to the suppression of ferroelectricity,(12) *e.g*. in
SrTiO.(13)

Here we adopt the bosonic scenario and concentrate on systems which undergo at finite temperature a normal state to superconductor transition with Berezinskii-Kosterlitz-Thouless (BKT) critical behavior,(14); (15) originally derived for the 2D xy-model with a two component order parameter and short range interactions. The occurrence of BKT criticality in 2D superconductors also implies that the mean square of the order parameter does not vanish at and with that there are, in analogy to He, condensed pairs (bosons) below and uncondensed ones above . In He and superconductors the order parameter is a complex scalar corresponding to the components in the xy-model. Supposing that in superconductors the interaction of Cooper pairs is short ranged and their effective charge is sufficiently small the critical properties at finite temperature are then those of the 3D-xy (bulk) and 2D-xy (thin films) models,(16); (17) reminiscent to the lamda transition in bulk He(18) and the BKT-transition in thin He films.(19); (20); (21); (22) In this context it is important to recognize that the existence of the BKT-transition (vortex-antivortex dissociation instability) in He films is intimately connected with the fact that the interaction energy between vortex pairs depends logarithmic on the separation between them. As shown by Pearl,(23) vortex pairs in thin superconducting films (charged superfluid) have a logarithmic interaction energy out to the characteristic length = , beyond which the interaction energy falls off as . Here is the magnetic penetration depth and the film thickness. As increases by approaching the diamagnetism of the superconductor becomes less important and the vortices in a clean and thin superconducting film become progressively like those in He films.(24); (25)

The occurrence of a 2D-QSI transition implies a line of BKT transition temperatures ending at the quantum critical point at where . It separates the superconducting from the insulating ground state. The BKT transition is rather special because the correlation length diverges above as and the sheet resistance tends to zero according to .(24); (25); (26) Approaching the 2D-QSI transition quantum phase fluctuations renormalize , , and . Indeed the BKT-transition line approaches the 2D-QSI transition as because the quantum scaling form exhibits at the universal value of the scaling argument a finite temperature singularity.(6); (8); (16) Noting that the amplitude of the BKT correlation length should match the divergence of the quantum counterpart, the exponents and should emerge from the amplitude in terms of .(26) Accordingly, the nonuniversal functions and entering the BKT expression for the sheet resistance exhibit close to the QSI transition quantum critical properties disclosing the quantum critical exponents , and . The exponents and are characteristic properties of the universality class to which the QSI transition belongs. In addition, given their values the relevance of disorder and the equivalence between quantum phase transitions in systems with spatial dimensions and the ones of classical phase transitions in dimensions can be checked.(5); (6); (16) The occurrence of a BKT transition line also implies: a nonuniversal critical sheet resistance because is nonuniversal; an explicit form of the superconductor branch of the quantum scaling function .

Even though BKT critical behavior is not affected by short-range correlated and uncorrelated disorder(27); (28) the observation of this behavior requires sufficiently homogeneous films and a tuning parameter which does not affect the disorder. Noting that sample inhomogeneity and vortex pinning are relevant in thickness and perpendicular magnetic field tuned transitions, electrostatic tuning of the 2D-QSI transition using the electric field effect appears to be more promising.(26); (29); (30); (31); (32) Indeed electrostatic tuning is not expected to alter physical or chemical disorder, but changes the mobile carrier density.

In Sec. II we sketch the theoretical background. To illustrate the potential and the implications of finite temperature BKT criticality on the 2D-QSI transition we analyze in Sec. II the sheet resistance data of Bollinger et al.(30) taken on gate voltage tuned epitaxial films of LaSrCuO that are one unit cell thick. Commenting on the difficulties in observing the BKT features in the magnetic penetration depth we consider the data of Bert et al.(32) taken on the superconducting LaAlO/SrTiO interface.

## Ii Theoretical background

Continuous quantum-phase transitions (QPT) are transitions at zero temperature in which the ground state of a system is changed by varying a parameter of the Hamiltonian.(5); (16); (17) The quantum superconductor- insulator transitions (QSI) in two-dimensional (2D) systems tuned by disorder, film thickness, magnetic field or with the electrostatic field effect are believed to be such transitions.(1); (2); (3); (26); (16); (30); (17) Traditionally the interpretation of experimental data taken close to the 2D-QSI transition is based on the quantum scaling relation,(5); (6); (7); (8)

(1) |

where is the resistance per square, the limiting ( and ) resistance. denotes the tuning parameter and a nonuniversal coefficient of proportionality. is a universal scaling function of its argument such that . In addition, is the dynamic and the critical exponent of the correlation length, supposed to diverge as

(2) |

The critical sheet resistance separating the superconducting and insulating ground states is determined from the isothermal sheet resistance at the crossing point in vs. tuning parameter at . The existence of such a crossing point, remaining temperature independent in the zero temperature limit, is the signature of a QPT. The data for plotted vs. should then collapse onto two branches joining at . The lower branch stems from the superconducting and the upper one from the insulating phase . This scaling form follows by noting that the divergence of the zero temperature correlation length, , is at finite temperature limited by the length .(5) Thus is a finite-size scaling function because . Supposing that there is a line of finite temperature phase transitions ending at the quantum critical point the quantum scaling form (1) exhibits at the universal value of the scaling argument a finite temperature singularity.(6); (8); (16) The phase transition line is then fixed by

(3) |

Otherwise one expects that sufficiently homogeneous 2D superconductors exhibit at the superconductor to normal state transition BKT critical behavior.(26); (14); (15) Note that there is the Harris criterion,(27) stating that short-range correlated and uncorrelated disorder is irrelevant at the unperturbed critical point, provided that , where is the dimensionality of the system and the critical exponent of the finite-temperature correlation length. With and , appropriate for the BKT transition, (14); (15) this disorder should be irrelevant. Given a BKT superconductor to normal state transition the sheet resistance scales for as(24); (25); (26)

(4) |

allowing to probe the characteristic BKT correlation length(14); (15); (25); (26)

(5) |

While , , and depend on the tuning parameter and are subject to quantum fluctuations the characteristic BKT form of the correlation length and with that of the sheet resistance applies for any . Through standard arguments (see, e.g., Ref. (33)) quantum mechanics does not modify universal finite temperature properties. is given by(15); (26)

(6) |

with . The parameter is expected to remain constant in the low regime.(26) It is related to the vortex core energy in terms of (34) and controls below the Nelson-Kosterlitz jump the temperature dependence of the magnetic penetration depth in terms of .(35) The amplitude of the BKT correlation length is proportional to the vortex core radius(19) known to increase with reduced .(20); (22) Indeed, the zero temperature correlation length diverges as (Eq. (2) and combined with (Eq. (3)) we obtain . Noting that approaches from above the scaling relation

(7) |

should apply,(26) making the determination of the exponents and possible. Consistency requires that the resulting agrees with the estimate obtained from the critical BKT line . Other implications concern the universality class of the 2D-QSI transition and the relevance of disorder. Given estimates for and the equivalence between quantum phase transitions in clean systems with spatial dimensions and the ones of classical phase transitions in dimensions can be checked. The fate of a clean critical point under the influence of disorder is controlled by the Harris criterion:(27); (28) If the zero temperature correlation length critical exponent fulfills the inequality the disorder does not affect the critical behavior. If the Harris criterion is violated, , the generic result is a new critical point with conventional power law scaling but new exponents which fulfill the Harris criterion. Another option is that the disorder destroys the QSI transition.

Finally, given a BKT transition line expression (4) the sheet resistance given by Eq. (5) transforms with Eq. (6) and the scaling variable to

(8) |

valid foe . Close to quantum criticality where reduces to

(9) | |||||

These relations are explicit forms of the quantum scaling function applicable to the superconductor branch. They reveal that the critical sheet resistance is the endpoint of a nonuniversal function and accordingly nonuniversal. Another implication concerns the universality class of the 2D-QSI transition. Supposing that the equivalence between quantum phase transitions in clean systems with spatial dimensions and the ones of classical phase transitions in dimensions applies, the 2D-QSI transition at the endpoint of a BKT line should belong to the finite temperature universality class. denotes an order parameter with two components, including the complex scalar, , of a superconductor.(5); (16)

The BKT - theory of thermally-excited vortex-antivortex pairs also predicts a super-to-normal state phase transition marked by the Nelson-Kosterlitz jump, a discontinuous drop in superfluid density from(35),

(10) |

to zero. The numerical relationship applies for in cm and in K. denotes the thickness of the 2D system, the in-plane magnetic penetration depth, and . In addition there is the prediction that , a measure of the phase stiffness, scales near the endpoint of the BKT transition line as(16); (17); (6); (8)

(11) |

provided that is below the upper critical dimension where hyperscaling holds. Since in the -xy-model the validity of Eq. (11) is in restricted to . Relation (11) is the quantum counterpart of the Nelson-Kosterlitz relation (8). is a dimensionless critical amplitude with the lower bound(16); (17)

(12) |

dictated by the characteristic temperature dependence of below the Nelson-Kosterlitz jump (Eq. (10)).(16); (17) Combining Eqs. (10) and (11) we obtain

(13) |

where is the superfluid density. The superfluid transition temperature as a function of the superfluid density has been measured in He films for transition temperatures ranging from to K by Crowell et al.(21) They studied He films adsorbed in two porous glasses, aerogel and Vycor, using high-precision torsional oscillator and dc calorimetry techniques. The investigation focused on the onset of superfluidity at low temperatures as the He coverage is increased. Their data yields with in moles/m and in K. Combined with the BKT transition line we obtain .

## Iii Comparison With Experiment

To illustrate the potential and the implications of the outlined BKT scenario we analyze next the data of Bollinger et al.(30) taken on epitaxial films of LaSrCuO that are one unit cell thick. Very large electric fields and the associated changes in surface carrier density enabled shifts in the midpoint transition temperature by up to K. Hundreds of resistance versus temperature and carrier density curves were recorded and shown to collapse onto a single function, as the quantum scaling form (Eq. (1)) for a 2D-QSI transition predicts. The observed critical resistance is close to the quantum resistance for pairs, k. Our starting point is the temperature dependence of the sheet resistance taken at various gate voltages where a BKT transition is expected to occur. As an example we depicted in Fig. 1a the data for V. The observation of BKT-behavior requires that the data extend considerably below the mean-field transition temperature . We estimated it with the aid of the Aslamosov-Larkin (AL) expression(36) for the conductivity, , with , where Gaussian fluctuations are taken into account and is the mean-field transition temperature. The resulting temperature dependence is included in Fig. 1a. It clearly reveals that the data extend considerably below K. To establish and characterize BKT behavior below we invoke Eq. (4) and

(14) |

As indicated in Fig. 1b this relation is used to fix and while is estimated by adjusting Eq. (4) with given and to the sheet resistance data. Comparing the data with the respective lines we observe that the BKT regime is attained and that hat the BKT is almost an order of magnitude lower than the mean-field counterpart. This uncovers a BKT transition from uncondensed to condensed Cooper pairs driven by strong phase fluctuations. On the other hand it is important to recognize that agreement with BKT-criticality is established in a temperature window only. Its upper bound reflects the crossover from BKT- to AL-behavior while the lower bound stems from the rounded BKT-transition. Precursors of this phenomenon are clearly visible in Fig. 1 around K. Here the correlation length is prevented to grow beyond a limiting length , i.e. the linear extent of the homogeneous domains. As a result, a finite size effect and with that a rounded transition occurs.(26) Because the BKT correlation length does not exhibit the usual power law divergence of the correlation length as is approached, it is particularly susceptible to such finite-size effect.

To explore the quantum critical behavior latent in , and we performed this analysis of the temperature dependence of the sheet resistance for additional gate voltages. To fix the critical gate voltage we used the empirical gate voltage dependence of the number of mobile holes per one formula unit of Bollinger et al.(30) yielding down to V with and V. The results , are shown in Fig. 2a and Fig. 2b while is depicted in Fig. 3. As can be seen in Fig. 2a the BKT transition line differs substantially from the so called superconducting dome behavior observed in bulk cuprate superconductors. It is approximately given by where is the maximum .(37) Close to the QSI transition where even bulk cuprate superconductors become essentially 2D(17) it reduces to with and suggests . In any case it differs substantially from the BKT line shown in Fig. 2a yielding the estimate

(15) |

in agreement with the value derived by Bollinger et al.(30) using the quantum scaling approach. In contrast to this, from the nonuniversal parameters and shown in Fig. 2b we derive with Eq. (7) in addition

(16) |

As these exponents satisfy the inequality we use the correct scaling argument, because where denotes the schemical potential.(6) The agreement between these values confirms the applicability of the scaling relation (7). Similarly, the nonuniversal parameter exhibits according to Fig. 3 the expected dependence (Eq. (6)). Noting that exceeds the upper critical dimension the critical exponent of the zero temperature correlation length should adopt its mean-field value . However the fate of this clean critical point under the influence of disorder is controlled by the Harris criterion.(27); (28) If the inequality is fulfilled, the disorder does not affect the critical behavior. If the Harris criterion is violated ( ), the generic result is a new critical point with conventional power law scaling but new exponents which fulfill . Since violates this inequality disorder is relevant and drives the system from the mean-field to an other critical point with different critical exponents as our estimate , fulfilling , uncovers. The resulting 2D-QSI transition with and violates then the equivalence between quantum phase transitions in systems with spatial dimensions and the ones of classical phase transitions in dimensions. In addition the proportionality between and (Eq. (11)) valid below the upper critical dimension does no longer hold because is above . In fact magnetic penetration depth measurements taken on underdoped high quality YBaCuO single crystals revealed .(38) Another option is that the disorder destroys the QSI transition. Given the evidence for a rounded BKT transition (see Fig. 1) and the missing data close to the QSI transition (see Fig. 4) further studies are required to elucidate this option.

To complete the analysis of the data of Bollinger et al.(30) we depicted in Fig. 4 the plot vs.corresponding to the quantum scaling function (Eq. (8)) in terms of . Apparently, the data does not fall completely on the BKT curve indicated by the dashed line. It corresponds to the superconductor branch of the quantum scaling function. Instead we observe a flow to and away from the universal characteristics. As decreases for fixed the crossover to AL-behavior sets in, while the rounding of the transition leads with increasing to a flow away from criticality. The important lesson then is that the quality of the data collapse on a single curve heavily depends on the temperature range of the data entering the plot. Another striking feature is the extended scaling regime. Within the BKT scenario it simply follows from the fact that the scaling form (4) applies along the BKT transition line irrespective of the distance from the QSI transition. In the quantum scaling approach this property remains hidden and merely suggests an extended quantum critical regime. The excellent quality of the piecewise data collapse also reveals that the observance of the substantial variation of (see Fig. 2b) is essential, while in the quantum scaling approach it is fixed by the critical sheet resistance.

To demonstrate these features even more compelling we depicted in Fig. 5a vs.. Apparently the data do not fall even piecewise on a single curve. After all this is not surprising because the quantum scaling form (1) holds close to the critical sheet resistance only and varies substantially in the gate voltage regime considered here (see Fig. 2b). As shown in Fig. 5b and Fig. 6 this behavior allows to estimate from vs. by rescaling in terms of with to achieve piecewise a collapse of the data. The resulting is shown in Fig. 6 and agrees well with V derived from the BKT behavior of the sheet resistance (see Fig. 2b).

A BKT line with a QSI transition at its endpoint was also explored rather detailed at the interface between the insulating oxides LaAlO and SrTiO exhibiting a superconducting 2D electron system that can be modulated by a gate voltage.(26); (39); (40); (32) BKT behavior and with that a 2D electron system was established as follows: The the current-voltage characteristics(39) revealed at the BKT transition temperature the characteristic BKT form with . (25) Consistency with the characteristic temperature dependence of the sheet resistance (Eq. (4)) was established.(26); (39); (40) It was also shown that the effective thickness of the superconducting 2D system can be extracted from the magnetic field dependence of the conductivity at .(41) The gate voltage tuned BKT phase transition line, , derived from the temperature dependence of the sheet resistance at various gate voltages uncovered with Eq. (4) consistency with K pointing to quantum critical behavior (Eq. (3)) with and the critical sheet resistance k.(26); (40) Furthermore the estimates and , confirming , have been derived from . (26) This suggests that the gate voltage tuned QSI transition of the 2D electron system at the LaAlO/ SrTiO interface belongs to the -xy universality class where hyperscaling and with that the proportionality between and (Eq. (11)) applies. On the other hand because disorder is according to Harris theorem relevant as well.(27); (28)

To comment the difficulties in observing the BKT features in the magnetic
penetration depth as well as the relation between and , we reproduced in Fig. 7 the data of Bert
et al.(32) for a gate voltage tuned superconducting LaAlO/SrTiO interface in terms of vs. . is defined here as the temperature at
which the diamagnetic screening drops below the noise level corresponding to
a detectable of cm. A glance at Fig. 7 reveals that this is just the regime where the universal quantum
behavior applies, indicated by the dashed and solid lines, corresponding to
the lower bound (12) and the behavior derived from the He data
of Crowell et al.,(21) respectively. Nevertheless the data
reveal the flow to the 2D-QSI transition which is attained at much lower ’s. Otherwise the data points resemble the outline of a fly’s wing,(17) remarkably similar to the vs. plots of the bulk superconductors YCu-123, Tl-1212,(42) and Tl-2201(43), covering nearly the
doping regime of the so called superconducting dome extending from the
underdoped to the overdoped limit . According to the generic plot *vs*. , where is the
anisotropy and denote the in-plane and c-axis correlation
length, these cuprates become nearly 2D in the underdoped limit.(17)
In any case, Fig. 7 shows that in this plot the universal QSI
behavior is attained at comparatively low values only. Accordingly,
in the QSI regime of interest the Nelson-Kosterlitz jump given by Eq. (10) becomes very small and appears to be beyond present experimental
resolution.(32) On the contrary in the temperature dependence of the
sheet resistance is the BKT critical regime accessible because ,(41) even though small compared to that in
the LaSrCuO films where .

## Iv Summary and discussion

In sum, we sketched and explored the implications of Berezinskii-Kosterlitz-Thouless (BKT) critical behavior on the quantum critical properties of a two dimensional (2D) quantum superconductor-insulator transition (QSI) driven by the tuning parameter . It was shown that the finite temperature BKT scenario, implies in terms of the characteristic temperature dependence of the BKT correlation length an explicit quantum scaling function for the sheet resistance along the superconducting branch ending at the nonuniversal critical value . This scaling form fixes the BKT-transition line and provides estimates for the quantum critical exponent product . In addition, independent estimates of , and follow from the dependence of the nonuniversal parameters and entering the characteristic BKT expression for the sheet resistance . This requires that the BKT critical regime where phase fluctuations dominate is attained and the finite temperature BKT relation for the sheet resistance applies for any . The last condition is satisfied because the BKT expression for the sheet resistance is simply related to the characteristic temperature dependence of the BKT correlation length. Quantum fluctuations enter via the nonuniversal parameters and disclosing close to the respective quantum critical behavior. Even though BKT critical behavior is not affected by short-range correlated and uncorrelated disorder(27); (28) the observation of this requires sufficiently homogeneous films and a tuning parameter which does not affect the disorder. Noting that in the magnetic field tuned case there is no BKT-line, thickness and gate voltage tuned 2D-QSI transitions appear to be promising candidates. In any case the scenario outlined here requires a line of BKT-transitions with a quantum critical endpoint and sheet resistance data which attain the BKT critical regime.

To illustrate the potential and the implications of this scenario we analyzed data of Bollinger et al.(30) taken on gate voltage tuned epitaxial films of LaSrCuO that are one unit cell thick. Evidence for dominant phase fluctuations and BKT-critical behavior was established in terms of the temperature dependence of the sheet resistance revealing a large critical regime extending substantially above the lowest attained temperature K. From the nonuniversal parameters and disclosing the respective quantum critical properties we derived for the quantum critical exponents the estimates: from , and from , yielding . Thus, in contrast to the standard quantum scaling approach, providing an estimate for only, the BKT scenario uncovered , and from the quantum critical behavior disclosed in and . Additional evidence for was established from the comparison of the scaled data with the explicit scaling BKT scaling form of the superconductor branch. We observed that the scaled data does not fall entirely on the BKT curve. Instead a flow to and away from the universal characteristics occurred. As decreases for fixed a crossover to AL-behavior sets in, while the rounding of the transition leads with increasing to a flow away from criticality. The important lesson then is that the quality of the data collapse on a single curve heavily depends on the temperature range of the data entering the plot. Another striking feature is the extended scaling regime. Within the BKT scenario it follows from the fact that the explicit scaling form of the sheet resistance applies along the entire BKT transition line irrespective of the distance from the QSI transition. In the quantum scaling approach this property remains hidden and merely suggests an extended quantum critical regime. The piecewise excellent quality of the data collapse also reveals that the provision of the substantial variation of is essential, while in the quantum scaling approach it is fixed by the critical sheet resistance. Supposing that the equivalence between quantum phase transitions with spatial dimensions and the ones of classical phase transitions in dimensions applies, the 2D-QSI transition at the endpoint of a BKT line should belong to the finite temperature universality class. Our estimate and with that critical behavior where . However the fate of this clean critical point under the influence of disorder is controlled by the Harris criterion.(27); (28) If the inequality is fulfilled, the disorder does not affect the critical behavior. If the Harris criterion is violated ( ), the generic result is a new critical point with conventional power law scaling but new exponents which fulfill . Since violates this inequality disorder is relevant and drives the system from the mean-field to an other critical point with different critical exponents as our estimate, , consistent with , uncovers. The resulting 2D-QSI transition with and violates then the equivalence between quantum phase transitions in systems with spatial dimensions and the ones of classical phase transitions in dimensions. In addition the proportionality between and (Eq. (11)), valid below the upper critical dimension , does no longer hold because is above . In fact magnetic penetration depth measurements taken on underdoped high quality YBaCuO single crystals revealed .(38) Another option is that the disorder destroys the QSI transition. Given the evidence for a rounded BKT transition (see Fig. 1) and noting that the analyzed data do not extend very close to the QSI transition (see Fig. 4) further studies are required to elucidate this option. In any case, our estimate , the evidence for BKT behavior at finite temperature and the Harris theorem imply the presence and relevance of disorder at the QSI transition and its irrelevance at finite temperature. Further evidence for the importance of disorder follows from the temperature dependence of the sheet resistance in the insulating phase. Considering of Bollinger et al.(30) we observe consistency with the Mott variable range hopping model in 2D. The conductivity exhibits the characteristic temperature dependence which applies to strongly disordered systems with localized states.(44)

A previous and analogous analysis of the sheet resistance data of the superconducting LaAlO/ SrTiO interface revealed the critical sheet resistance k and the exponents , and .(26) In this case the gate voltage tuned QSI transition of the 2D electron system at the LaAlO/ SrTiO interface has a finite temperature counterpart in dimensions namely the -xy model where hyperscaling and with that the proportionality between and (Eq. (11)) applies. On the other hand implies according to Harris theorem(27); (28) that disorder is relevant as well.

To comment on the BKT-features in the temperature dependence of the magnetic penetration depth we considered the data of Bert et al.(32) for a gate voltage tuned superconducting LaAlO/SrTiO interface in terms of vs. . The data reveals the flow to the universal relationship (11) but much lower must be attained to reach quantum critical regime. However in this low regime is the Nelson-Kosterlitz jump given by Eq. (10) very small and appears to be beyond present experimental resolution.(32) Otherwise the data points resemble the outline of a fly’s wing,(17) remarkably similar to the vs. plots of the bulk superconductors YCu-123, Tl-1212,(42) and Tl-2201(43), covering nearly the entire doping regime in the so called superconducting dome extending from the underdoped to the overdoped limit.

### Footnotes

- preprint: PREPRINT (March 6, 2018)

### References

- A. M. Goldman, Physica E 18, 1 (2003).
- N. Marković, C. Christiansen, A. Mack, and A. M. Goldman, Phys. Status Solidi B 218, 221 (2000).
- V. F. Gantmakher, V. T. Dolgopolov, Phys. Usp. 53, 1 (2010).
- A. M. Goldman, Int. J. Mod. Phys. B 24, 408 (2010).
- S. L. Sondhi, S. M. Girvin, J. P. Carini, and D. Shahar, Rev. Mod. Phys. 69, 315 (1997).
- M. P. A. Fisher, P. B. Weichman, G. Grinstein, and D. S. Fisher, Phys. Rev. 40, 546 (1989).
- M. P. A Fisher, G. Grinstein, and S. M. Girvin, Phys. Rev. Lett. 64, 587 (1990).
- Kihong Kim and Peter B. Weichman, Phys. Rev. B 43, 13583 (1991).
- S. Maekawa and H. Fukuyama, J. Phys. Soc. Jpn. 51, 1380 (1981).
- A. M. Finkel’shtein, Sov. Phys. JETP. 45, 46 (1987).
- P. B. Weichman and R. Mukhopadhyay, Phys. Rev. B 77, 214516 (2008).
- R. Morf, T. Schneider, and E. Stoll, Phys. Rev. B 16, 462 (1977).
- K. A. Müller and H. Burkhard, Phys. Rev. B 19, 3593 (1979).
- V. L. Berezinskii, Sov. Phys. JETP 32, 493 (1971).
- J. M. Kosterlitz and D. J. Thouless, J. Phys. C 6, 1181 (1973) .
- T. Schneider and J. M. Singer, Phase Transition Approach to High Temperature Superconductivity (Imperial College Press, London, 2000).
- T. Schneider, in The Physics of Superconductors, edited by K. Bennemann and J. B. Ketterson (Springer, Berlin, 2004) , p. 111.
- F. M. Gasparini, M. O. Kimball, K. P. Mooney, and M. Diaz-Avila, Rev. Mod. Phys. 80, 1009 (2008).
- D. Finotello and F. M. Gasparini, Phys. Rev. Lett. 55, 2156 (1985).
- H. Cho and G. Williams, Phys. Rev. Lett. 75, 1562 (1995).
- P. A. Crowell, F. W. Van Keuls, and J. D. Reppy, Phys. Rev. B 55, 12620 (1997).
- T. Oda, M. Hieda, R. Toda, T.Matsushita, and N. Wada, J. Low Temp. Phys. 158, 262 (2010).
- J. Pearl, Appl. Phys. Lett. 5, 65 (1964).
- P. Minnhagen, Rev. Mod. Phys. 59, 1001 (1987).
- B.I. Halperin, D.R. Nelson, J. Low Temp. Phys. 36, 599 (1979).
- T. Schneider, A. D. Caviglia, S. Gariglio, N. Reyren, and J.-M. Triscone, Phys. Rev. B 79, 184502, (2009).
- A. B. Harris, J. Phys. C 7, 1671 (1974).
- J. T. Chayes, L. Chayes, D. S. Fisher, and T. Spencer, Phys. Rev. Lett. 57, 2999 (1986) .
- Kevin A. Parendo, K. H. Sarwa B. Tan, and A. M. Goldman, Phys. Rev. B 73, 174527 (2006).
- A. T. Bollinger, G. Dubuis, J. Yoon, D. Pavuna, J. Misewich, and I. Božoviĉ, Nature, 472, 458 (2011).
- Xiang Leng, Javier Garcia-Barriocanal, Shameek Bose, Yeonbae Lee, and A. M. Goldman, Phys. Rev. Lett. 107, 027001 (2011).
- Julie A. Bert, Katja C. Nowack, Beena Kalisky, Hilary Noad, John R. Kirtley, Chris Bell, Hiroki K. Sato, Masayuki Hosoda, Yasayuki Hikita, Harold Y. Hwang, and Kathryn A. Moler, Phys. Rev B 86, 060503(R) (2012).
- J. A. Hertz, Phys. Rev. B 14 (1976).
- Lindsay M. Steele, Ch. J. Yeager, and D. Finotello, Phys. Rev. Lett. 71, 3673 (1993).
- D. R. Nelson and J. M. Kosterlitz, Phys. Rev. Lett. 39, 1201 (1977).
- L. G. Aslamosov and A. I. Larkin, Sov. Phys. Solid State 10,875 (1968).
- M. R. Presland, J. L. Tallon, R. G. Buckley, R. S. Liu, and N. E. Flower, Physica C 176, 95 (1991).
- R. Liang, D. A. Bonn, W. N. Hardy, and D. Bonn, Phys. Rev. Lett. 91, 117001 (2005).
- N. Reyren, S. Thiel, A. D. Caviglia, L. Fitting Kourkoutis, G. Hammerl, C. Richter, C. W. Schneider, T. Kopp, A.-S. Rüetschi, D. Jaccard, M. Gabay, D. A. Muller, J.-M. Triscone, and J. Mannhart, Science, 317, 1196 (2007).
- A. D. Caviglia, S. Gariglio, N. Reyren, D. Jaccard, T. Schneider, M. Gabay, S. Thiel, G. Hammerl, J. Mannhart, and J.-M. Triscone, Nature 456, 624 (2008).
- T. Schneider, Phys. Rev. B 80 , 214507 (2009).
- C. Bernhard, J. L. Tallon, Th. Blasius, A. Golnik, and Ch. Niedermayer, Phys. Rev. Lett. 86, 1614 (2001).
- Ch. Niedermayer, C. Bernhard, U. Binninger, H. Glückler, J. L. Tallon, E. J. Ansaldo, and J. I. Budnick, Phys. Rev. Lett. 71, 1764 (1993).
- N. F. Mott, Phil. Mag. 19, 835 (1969).