Introduction

Oxide heterostructures have attracted much attention in recent years due to a plethora of novel or enhanced physical phenomena observed at interfaces1,2. For example, the two dimensional electron gas (2DEG) formed at the interface between two perovskite insulators, LaAlO3 (LAO) and SrTiO3 (STO), becomes superconducting3 at around 200 mK. This discovery triggered intense investigations of the superconductivity mechanism4,5,6,7,8 at interfaces. A number of intriguing properties emerge at the superconducting interface. Literature indicate on a pseudogap-like behavior9, on the electron pre-formed pairing10 and on the coexistence of superconductivity and ferromagnetism11,12,13. The phase diagram for an interface is similar to that of high temperature superconducting (HTS) cuprates7,9. The interface exhibits superconductivity for an extremely low carrier density14,15. Presented aspects may indicate on unconventional pairing mechanism. Nevertheless, the conventional electron-phonon coupling has been also considered5,16. Therefore, the microscopic mechanism responsible for interface superconductivity is still unclear and more research is of high interest.

Recently, it was reported that an energy gap persists above the superconducting dome defined as the variation of the critical temperature, Tc, vs. carrier doping. Furthermore, the energy gap increases monotonically as the gate voltage VG decreases. This correlation indicates on a pseudogap-like behavior of (001) LAO/STO interfaces9 and it is opposite to the VG dependence of the superfluid density14. For additional information about electron pairing in interfaces, it is useful to determine fundamental superconducting parameters such as the coherence length, ξ0 and the strength of pairing potential. These parameters can be extracted from upper critical field, Hc2; a higher Hc2 means a smaller coherence length and a stronger pairing potential17. Previous investigations of Hc2 have shown that superconductivity observed for (001) or (110) LAO/STO interfaces is of two dimensional nature3,18,19,20, but the VG dependence of Hc2 that is associated to pairing strength while being distinct from the superfluid density, has not been reported. Noteworthy is also that most studies are focused on (001) LAO/STO, while the (110) system is less explored although its study can bring an important contribution to understanding of superconductivity at interfaces20,21,22.

In this work we report the dependence of Hc2 on VG in a superconducting 2DEG (110) LAO/STO interface. We found that Hc2 decreases when VG increases. At the same time, the Tc(VG) curve shows a dome-like shape (with a maximum). Our results indicate that the Cooper pairing potential becomes stronger in the underdoped regime.

Results and Discussions

Critical temperature, critical current and normal state resistance as a function of gate voltage

In Fig. 1 are presented as examples, curves of (interface or sheet) resistance, Rs, versus temperature for three different gate voltages (−5 V, 0 V and 25 V) under zero magnetic field. Temperature for which Rs drops below the detection limit of the measuring equipment defines the critical temperature, (Fig. 1). The curve of as a function of VG (in zero magnetic field) has a dome-like shape with a maximum at 25 V (Fig. 1 inset). The dome-like shape of the (VG) curves was previously observed in LAO/STO and (001) LaTiO3/SrTiO3 (LTO/STO) interfaces7,20,23. A dome-like shape with a maximum at VG = 25 V is also obtained for the variation of the critical current Ic(VG). The critical current, Ic, was determined at 50 mK from I-V curves measurements (see Supplementary Information, Fig. S1). The (VG) and Ic(VG) similar dome-like dependencies indicate a strong influence of the carrier density on and Ic.

Figure 1
figure 1

The resistance Rs of the (110) LAO/STO interface as a function of temperature at three representative gate voltages and for B = 0 T.

The inset shows the gate voltage VG dependence of the critical temperature and of the normal state resistance Rn. The maximum of the (VG) curve is at 138 mK and at 25 V. Lines are guide to the eyes.

From the measurements of Rs with magnetic field (up to 1.6 T) applied perpendicular to the surface of the interface at a fixed temperature and gate voltage one can determine the normal state resistance, Rn. An example of Rs(B) curves at 50 mK and for VG from −25 to +200 V is presented in Fig. 2a. Inset to Fig. 2a shows for selected (T, VG) values how Rn is determined. Namely, Rn is the resistance of the cross point between the fitting lines of the steepest part of the Rs(B) experimental curve and of the region where Rs almost saturates. For a fixed temperature (50 mK), one Rn(VG) curve is plotted in Fig. 1 inset. Enhancement of VG decreases Rn. Curve is non-linear and the decrease rate is smaller for higher VG.

Figure 2
figure 2

(a) Curves of Rs(B) for −5 ≤ VG ≤ 200 V at 50 mK. The inset in the right corner shows determination of Rn, and . The inset in the left corner shows (VG) and (VG) curves at 50 mK; (b) Curves of dRs/dB vs. B at 50 mK. Arrows indicate points of maximum of these curves and dashed line is guide for eyes. Magnetic field of a maximum point is defined as . The inset shows variation of as a function of VG.

It is remarkable that VG can tune superconducting and normal state characteristics such as and Rn, respectively. Results suggest that VG influences 2DEG superconductivity features and the carrier density. Our results are consistent with literature20.

Upper critical field as a function of gate voltage

For a fixed temperature and VG, the upper critical field Hc2 is determined as the field where resistance is 10%, 50% of Rn in the Rs(B) curves (Fig. 2a inset, right corner). This methodology was used to determine Hc2 for different superconducting systems24,25,26. We note that magnetic field B is applied in this work only perpendicular to the interface surface. Reported articles19,20 indicate that Hc2 shows large anisotropy when B is applied in-plane and out-of-plane. At 50 mK, the (VG) or (VG) (Fig. 2 inset, left corner) curves are unexpectedly without a dome-like shape that is specific for the (VG) and Ic(VG) (Figs 1 and S1) curves. Namely, for a lower VG, Hc2 monotonically increases (Fig. 2 inset, left corner). At the same time, a decrease of VG below VG, max = 25 V counter-intuitively produces the decrease of both and Ic (Figs 1 and S1).

The differential curves dRs/dB as a function of B, obtained from Rs(B) data (Fig. 2a) for different VG and at 50 mK, are presented in Fig. 2b. Each curve displays a peak (marked with an arrow in Fig. 2b). The peak shifts to higher B and its intensity increases when VG decreases. The magnetic field of the peak, denoted , as a function of VG (at 50 mK) is shown in Fig. 2b, inset. A dome-like shape is not obtained and, as expected, this curve has a similar behavior as (VG) or . (VG) curves. Furthermore, the values of and of are close to each other, but we shall keep both parameters in the discussion because of their different background: Hc250% is arbitrary taken, while Hc2peak defines the inflexion point of the Rs(B) cur. The inflexion point may have a physical meaning and, hence, Hc2peak can be more sensitive to external factors than .

Figure 3 shows curves of (T) for different VG from −5 V up to 200 V. For all gate voltages, increases with temperature decrease. Curves of . (T) shift into the region of higher values of -T for decreasing VG. This enhancement occurs even for negative VG, where a lower VG induces a lower . Moreover, the slope (d/dT|Tc) of the (T) curve approaching (Tc) = 0 systematically increases when VG decreases. The overall observed tendency does not change, if we use in the analysis or instead of (Fig. S2 a,b). According to Ginzburg-Landau (GL) theory, near Tc the upper critical field is a linear function of (Tc − T) and smoothly saturates when lowering the temperature. Our results suggest an anomalous unconventional behavior of the upper critical field vs. temperature and it is no longer appropriate to describe this dependence using the GL theory. The (T = 0 K) values obtained by linear extrapolation of e data at low temperature are increasing with decreasing VG. The determination of Hc2 at T = 0 K lacks precision. However, if we estimate the relative Hc2-increase (δHc2) at a finite temperature (e.g. 60 mK) between the curves for VG = −5 V and VG = 0 V, a criterion closer to Rs = 0 shows a lower value (i.e. δ = 1.3%, while δ = 5.5%). This finds its understanding in the following: The emergence of Rs = 0 is due to the global superconductivity. With increasing magnetic field, a finite resistance occurs, and, hence, global superconductivity disappears. However, local superconductivity still exists before the system recovers to normal state. Thus, it is reasonable that we observe a larger Hc2 for a more negative VG. This indicates that Cooper pair can persist up to a higher magnetic field than the upper critical field corresponding to global superconductivity.

Figure 3
figure 3

Temperature dependence of for different VG.

Lines are guide for the eyes.

Superconductor-insulator transition

The magnetic-field-induced superconductor-insulator transition(SIT)27 shows a characteristic fan-shaped pattern of Rs(B) isotherms crossing at one point. For example, in Fig. 4a the SIT cross point for VG = −5 V is at 810 mT. The Rs(T) curves extracted from Rs(B) show a plateau for 810 mT (Fig. 4b). The plateau separates two regimes. Therefore, the magnetic field drives a continuous quantum phase transition from a superconducting 2DEG to a weakly insulating state. Further finite-size scaling analysis shows that the data can be collapsed onto a bi-value curve (Fig. 4a, inset). It results that the crossing point is a quantum critical point (QCP), at which the phase transition occurs. Literature often show magnetic-field-induced SIT for different 2D superconductors27,28, including for interfaces such as LAO/STO and (001) LTO/STO29,30.

Figure 4
figure 4

(a) SIM pattern: curves of Rs(B) at different temperatures for V= −5 V. Note the crossing point at 810 mT that defines the quantum critical point (QCP) with (Bc = 810 mT, Rc = 1297.54 Ω). The inset shows the bi-value curve obtained by collapse of the date by finite-size scaling analysis. (b) Curves of Rs(T) at different magnetic fields B and for VG = −5 V. Dashed line shows a plateau corresponding to QCP. Continuous lines are guide for the eyes.

As already noted, Hc2(T) for a fixed VG depends on the criterion adopted for the Hc2 determination in the case of a broad transition26,31 and the Hc2(T = 0 K) cannot be determined in a reliable manner. This situation questions the intrinsic nature of the Hc2 – enhancement (for a lower VG) in the underdoped region. Due to this, here we use another method to directly determine the zero - temperature upper critical field, . This method is independent of the criterion applied to Rs(B) curves. At the QCP, taken as a transition point between superconducting and normal states of the 2DEG, the corresponding magnetic field is defined as . The curve as a function of VG is plotted in Fig. 5. Although it is considered that this method does not give accurate values of 30, this barely affects our main results, the dependence on VG of being similar to that of (T = 0 K) (Fig. 5). One reason for a non-precise determination of vs. VG is: when VG ≠ −5 V, e.g., a VG = 75 V is used for the construction of the magnetic-field-induced SIT pattern, the crossing point of the Rs(B) isotherms transforms into a field domain centered at 0.55 T and extending over ±0.05 T (Fig. S3a). In this case, the conventional power-law scaling behavior fails to describe the quantum criticality (Fig. S3b). Multiple quantum criticality was also found and reported for the (001) LTO/STO interface29. The phenomenon of multiple critical exponents suggests an unconventional critical behavior of SIT in the (110) LAO/STO interface and it will be discussed elsewhere.

Figure 5
figure 5

The VG dependence of , (T = 0 K) (green diamonds) and .

(violet stars). The inset shows the VG dependence of the coherence length where is the quantum flux. Continuous lines are guide for the eyes.

Phase diagram

The superconducting phase diagram of the (110) LAO/STO interface is presented in Fig. 5. As already addressed in the previous sections, the upper critical field Hc2 independently of the criterion for its determination monotonically increases when VG decreases, while displays a dome-like curve with a maximum. In the inset to Fig. 5 is shown the GL coherence length ξ0 determined from Hc2 as a function of VG. The most striking anomalous result is enhancement of Hc2 accompanied by the decrease of ξ0 for decreasing VG, i.e. for the carrier depletion reflected by the dome-like curve in the underdoped region. Results indicate that the Cooper pairing potential is stronger in the underdoped region.

A systematic increase in Hc2 with decreasing doping has been reported both in high-Tc cuprates and iron-based superconductors31,32. This is usually considered as an evidence for the existence of a so-called ‘pseudogap’ state, in which the bosonic pairs form above Tc but cannot condensate into superconducting state due to dilution of pairs33. Recently, the planar tunneling spectroscopy study in 2DEG at a (001) LAO/STO interface has shown that the energy gap ∆ increases with charge carrier depletion in both underdoped and overdoped regions9. And, the coherence-peak-broadening parameter Γ derived from the Dynes fit, that is related to the strength of the superconducting pairing interaction, increases steeply with decreasing VG. GuangLei Cheng et al.10 recently reported that the Cooper pairs form at temperatures well above the superconducting transition temperature of the (001)LAO/STO superconducting system. In the experiments they used a superconducting single-electron transistor and they observed that pairs condensate at low magnetic fields and temperatures. The physical understanding of the processes at (110) LAO/STO interface is analogous to that of hole-doped cuprates31, namely, the pairing potential is stronger and the Ginzburg-Landau coherence length ξ0decreases in the underdoped region as VG decreases. The trend differs from that of superfluid density and superconducting transition temperature7,14. Namely, the superfluid density decreases with decreasing VG providing that the phase fluctuation is important14,34,35. The observation of ∆/Γ scaled with Tc in ref. 9 implies that the limited quasiparticle lifetime controls Tc effectively9 and the reduction in Tc versus ∆ was attributed to a competing order parameter or to a weak phase coherence. In addition, in LAO/STO system the spin-orbit coupling is non negligible and strongly depends on VG36,37. Both conventional and unconventional pairing mechanisms have been considered to describe superconductivity in interfaces5,6,16,38. For example, the spin-orbit coupling36,37 and the coexistence of superconductivity and ferromagnetism11,12,13 may indicate formation of possible exotic superconducting states such as finite momentum Cooper pairing39,40. One has also to consider a different orbital reconstruction between (001) and (110) systems20,41. Superconducting properties and Rashba spin-orbit coupling can be largely tuned by controlling selective orbital occupancy in different crystal orientations20. For the underdoped region, it has been reported that a Lishiftz transition is observed at (001) interface42,43. The (110) system is expected to be characterized only by a 3dxz/dyz filled electronic state; thus the superconducting properties of the (110) interface system can be substantially different from that of the (001) system.

Conclusions

We systematically investigated the upper critical field as a function of gate voltage by ultralow temperature magnetoresistance measurements in superconducting 2DEG of a (110) LAO/STO interface. We found that upper critical field increases as the gate voltage decreases. Two independent methods to determine the upper critical field give a similar trend. This implies that the pairing potential is stronger in the underdoped region. This observation is similar to recent reports that consider a pseudogap-like behavior at the (001) LAO/STO interface. Our results for an interface with a different orientation contribute to understanding of the pairing mechanism of superconductivity at LAO/STO interface.

Method

A five-unit-cell LaAlO3 thin film was grown on the (110) SrTiO3 substrate (500 μm thickness) by pulsed laser deposition. Details were described in ref. 19. A metallic back gate was evaporated and attached to the rear of the substrate. Leakage current was low (below the maximum value of 5 nA at VG = 200 V). Standard four-terminal resistance measurements were made using wedge-bonding contacts. The sample was cooled in a dilution refrigerator with a base temperature of 10 mK. The measurement current is sufficiently low (~50 nA) to avoid sample heating at ultralow temperatures. To ensure the reversible behavior of the superconductivity, the gate voltage was ramped up to 200 V after cooling down. Perpendicular magnetic field B was applied to the sample (interface) surface and the field direction is the same for all measurements.

Additional Information

How to cite this article: Shen, S. C. et al. Gate dependence of upper critical field in superconducting (110) LaAlO3/SrTiO3 interface. Sci. Rep. 6, 28379; doi: 10.1038/srep28379 (2016).