Evidence for charge-vortex duality at the LaAlO3/SrTiO3 interface

Abstract

The concept of duality has proved extremely powerful in extending our understanding in many areas of physics. Charge-vortex duality has been proposed as a model to understand the superconductor to insulator transition in disordered thin films and Josephson junction arrays. In this model, on the superconducting side, one has delocalized Cooper pairs but localized vortices; while on the insulating side, one has localized Cooper pairs but mobile vortices. Here we show a new experimental manifestation of this duality in the electron gas that forms at the interface between LaAlO₃ and SrTiO₃ . The effect is due to the motion of vortices generated by the magnetization dynamics of the ferromagnet that also forms at the same interface, which results in an increase in resistance on the superconducting side of the transition, but an increase in conductance on the insulating side.

ARTICLE

Received 28 Mar 2012 | Accepted 18 Jun 2012 | Published 17 Jul 2012 DOI: 10.1038/ncomms1959

M.M. Mehta1, D.A. Dikin1, C.W. Bark2, S. Ryu2, C.M. Folkman2, C.B. Eom2 & V. Chandrasekhar1

The concept of duality has proved extremely powerful in extending our understanding in many areas of physics. Chargevortex duality has been proposed as a model to understand the superconductor to insulator transition in disordered thin lms and Josephson junction arrays. In this model, on the superconducting side, one has delocalized Cooper pairs but localized vortices; while on the insulating side, one has localized Cooper pairs but mobile vortices. Here we show a new experimental manifestation of this duality in the electron gas that forms at the interface between LaAlO3 and SrTiO3. The effect is due to the motion of vortices generated by the magnetization dynamics of the ferromagnet that also forms at the same interface, which results in an increase in resistance on the superconducting side of the transition, but an increase in conductance on the insulating side.

Evidence for chargevortex duality at the LaAlO3/SrTiO3 interface

1 Department of Physics and Astronomy, Northwestern University, Evanston, Illinois 60208, USA. 2 Department of Materials Science and Engineering, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA. Correspondence and requests for materials should be addressed to V.C. (email: [email protected]).

NATURE COMMUNICATIONS | 3:955 | DOI: 10.1038/ncomms1959 | www.nature.com/naturecommunications

2012 Macmillan Publishers Limited. All rights reserved.

ARTICLE

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms1959

The two-dimensional superconductor to insulator transition (SIT) has been studied extensively13, as experimentally it is one of the few easily accessible manifestations of a quan

tum phase transition. The two-dimensional electron gas formed at the LaAlO3/SrTiO3 (LAOSTO) interface46 has been shown to undergo a SIT as the density of the carriers is tuned79. In addition to superconductivity, our earlier work showed evidence of ferromagnetism, which manifests itself as a hysteresis in measurements of the magnetoresistance (MR), Hall eect and superconducting phase boundary8. Two other recent studies have conrmed the coexistence of superconductivity and magnetism in similar samples10,11. Theoretical calculations suggest that the magnetism arises from localized moments on Ti ions at the interface12. Ferromagnetic order is stablized via RKKY interactions mediated by the conductance electrons. The experimental evidence obtained so far8,10 indicates that the eect of the ferromagnetism on the electrical transport properties of the system can be accounted for by considering only the magnetic eld produced by the ferromagnet.

In this paper, we show that the coexistence of ferromagnetism and superconductivity in this system leads to a novel manifestation of chargevortex duality in the SIT. On the superconducting side of the transition, the dynamics of the ferromagnet during magnetization reversal leads to a sharp peak in resistance whose magnitude increases with the sweep rate of the magnetic eld. On the insulating side of the transition, the same magnetization dynamics lead to a dip in resistance whose magnitude again increases with the magnetic eld sweep rate. This behaviour can be understood by modelling the superconductor as a random network of superconducting islands that responds to the dynamic magnetic eld arising from the ferromagnet.

ResultsMR in the superconducting state. The samples in this work had 10 unit cells (uc) of LAO grown by pulsed laser deposition on TiO2-terminated (001) STO single-crystal substrates13,14. Details of sample fabrication were discussed in an earlier publication8

(see also Methods). The interaction between ferromagnetism and superconductivity can be seen in MR measurements in the superconducting state. Figure 1 shows the MR in a magnetic eld parallel and perpendicular to the plane of the sample, respectively, at T=50 mK and Vg = 80 V, a gate voltage at which the sample is maximally superconducting. Consider rst the case of a parallel magnetic eld, H|| (Fig. 1a). As has been seen earlier for this sample in perpendicular magnetic elds8, the parallel eld MR is hysteretic due to the presence of ferromagnetic order (inset of Fig. 1a). As H|| is

in the plane of the two-dimensional superconductor, the critical eld is large (see Methods), so that one expects the MR to be small in this eld range. In fact, the overall background MR is indeed quite small: most of the contribution to the MR is due to a set of twin peaks at 12 mT. For clarity, we show in Fig. 1a the MR for only one sweep direction. As H|| is increased from negative values, there is a small decrease in the resistance of the sample until H|| = 0. The peak in

the MR develops just as H|| becomes positive for this eld sweep direction, but then rapidly dies out as H|| is increased further.

The most surprising fact about the MR of Fig. 1a is that the amplitude of the peak depends on the sweep rate of H||: the amplitude increases as the eld sweep rate is increased. For very slow sweep rates, the peak is almost absent, and the amplitude of the overall MR is very small, as expected from our earlier discussion.

The rate dependence of the MR can be understood if we consider the inuence of the domain wall of the ferromagnet on the superconductor, as shown schematically in Fig. 1b. In the simplest model, as H|| is swept from a large negative to a large positive value, the magnetization of the ferromagnetic layer reverses by forming a domain wall that propagates from one end of the sample to the other once H|| > 0. If the domain wall is a Bloch wall15, the

magnetization points out of the plane of the interface, and the conducting layer sees a magnetic eld perpendicular to the interface that is localized near the domain wall, and travels along with it. We emphasize that for the real sample, the magnetization dynamics are likely far more complicated, but the end result of a localized magnetic eld perpendicular to the plane of the interface that moves as H|| is swept is the same. At the eld sweep rates in the experiment, the motion of the domains is quasi-static. Consequently, the motion of the localized perpendicular magnetic eld is determined only by the sweep rate of the external magnetic eld. The moving perpendicular component of the eld in turn generates moving vortices in the superconductor. For our samples, which can be thought of as granular superconductors (see below), moving vortices will result in a change in resistance16 as the vortices cross the weak links between superconducting grains.

Similar behaviour is also observed for the MR in a perpendicular external magnetic eld H (Fig. 1c). Given the two-dimensional nature of the ferromagnet12, it is unlikely that the perpendicular magnetic elds applied can cant the moments appreciably out of plane. However, a Bloch domain wall may form at the interface, and it is domains within this Bloch wall that are reoriented in an external magnetic eld, as shown in Fig. 1d. Such behaviour has been seen before in thin magnetic lms17. (We stress again that this picture is highly simplied, and the real reversal process is likely much more complicated.) H itself leads to a large overall MR, with a minimum at zero eld. The eld from the ferromagnet gives rise to a sweep ratedependent peak at H~ 15 mT, and a second, less prominent peak at H~ 48 mT, which is likely due to more complex magnetization dynamics in perpendicular eld. Further evidence that the hysteresis and the dips are due to magnetization dynamics in the ferromagnet can be seen from the MR in crossed parallel and perpendicular elds. When one measures the background MR in a perpendicular eld in the presence of a large constant parallel eld that aligns the magnetization of the ferromagnet, no hysteresis or rate dependence is observed (see Methods).

MR in the insulating state. In the LAO/STO interface system, the superconducting transition can be tuned by Vg (refs 7 and 8). Early studies of LAO/STO interfaces identied this as a SIT7, which has been studied extensively in the past in thin superconducting lms as a function of the lm thickness or an applied magnetic eld1,2. However, in the LAO/STO system, it appears difficult to tune the system deep into the insulating regime with a gate voltage9. Nevertheless, measurements of the current-dependent dierential resistance (Methods) indicate that the transition is to a weakly insulating state where regions of superconductivity get increasingly isolated from one another as Vg is decreased. Thus, it makes sense to model our sample as a granular superconductor consisting of a network of superconducting islands whose phase is coupled through the Josephson eect, parametrized by the Josephson energy EJ and a Coulomb charging energy Ec, which is associated with the energy cost of adding a Cooper pair to an island3,18. (For a random network of islands, these quantities are averages over the network.) Changing Vg in our sample tunes the transition by modifying the ratio between EJ and Ec. The superconducting and insulating regimes are duals of each other, and this duality is predicted to manifest itself through an interchange between specic measurable quantities1821. For example, the current I and voltage V are interchanged across the SIT, so that the IV curve in the superconducting phase is similar in shape to the VI curve in the insulating state. Such dual IV characteristics have been observed in Josephson junction arrays and disordered thin lms3,22,23.

In our samples, as noted above, it appears difficult to tune the system deep into the insulating regime, so that observing this signature of chargevortex duality is difficult6,8. However, the interplay between ferromagnetism and superconductivity at the LAO/STO

NATURE COMMUNICATIONS | 3:955 | DOI: 10.1038/ncomms1959 | www.nature.com/naturecommunications

2012 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms1959

ARTICLE

a

25

c

25 20 15 10

5 0

600

20

15

10

R S ()

R S ()

60 30 0 30 60 30 T s1 60 T s1

120 T s1

240 T s1

300 T s1

180 T s1

60 T s1

120 T s1

360 T s1

480 T s1

600 T s1

240 T s1

5

0

80

60 40

20

0 20 40 60 80

500

400

300

200

100

0 150

100 50 0 50 100 150

H|| (mT)

H||

H||

H||

H (mT)

b

d

H

H

H

Figure 1 | MR in the superconducting regime. (a) Parallel eld MR as a function of different rates in the superconducting regime at Vg = 80 V. Data for only one eld sweep direction are shown for clarity. Arrow indicates the direction of eld sweep. The inset shows the MR for forward and backwardeld sweep directions at a sweep rate of 300 T s 1. (b) Schematic of the system at different parallel eld values. The top layer is the ferromagnet, the superconductor is shown through a gradient as extending some distance into the STO. In this simplied picture, magnetization reversal occurs by means of domain wall propagation in the ferromagnet. The perpendicular component of the eld due to the domain wall induces vortices in the superconductor. (c) Perpendicular eld MR for elds swept from negative to positive values for different eld sweep rates for Vg = 80 V. (d) Schematic of the magnetization state of the system at different perpendicular elds. Owing to the shape anisotropy of the system, the majority of the moments lie in plane, but the external magnetic eld orients the direction of the perpendicular component of the magnetization of the domain wall. It should be emphasized that the magnetization conguration in the real system during reversal is denitely far more complicated, but would still give rise to a perpendicular component of the magnetic eld. All data were taken at T = 50 mK.

interface results in a unique manifestation of chargevortex duality associated with moving ferromagnetic domain walls that has not been observed in other systems.

Figure 2a,b demonstrates this new behaviour. Figure 2a shows the equivalent of the data of Fig. 1a, but at a gate voltage of Vg = 100 V, placing the sample on the insulating side of the transition. At the same magnetic elds at which a sweep ratedependent peak was seen in Fig. 1a, a sweep ratedependent dip is observed. As with

the peak, the magnitude of the dip increases with increasing sweep rate. Similar behaviour is observed in perpendicular magnetic eld (Fig. 2b). Note that in perpendicular eld, one even observes two sets of dips (at 15 and 48 mT), mirroring the behaviour seen in Fig. 1c. In the chargevortex duality model, current and voltage are switched in going across the SIT, so that conductance and resistance are also switched. Thus, where one observes peaks in resistance on the superconducting side, one should observe peaks in conductance

NATURE COMMUNICATIONS | 3:955 | DOI: 10.1038/ncomms1959 | www.nature.com/naturecommunications

2012 Macmillan Publishers Limited. All rights reserved.

ARTICLE

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms1959

a

a

b

2,540

2,500

2,500

2,530

2,400

2,400

R S ( ) R S ( )

2,300

2,520

2,300

2,200

R S ()

R S ( ) R S ( )

100 V 80 V 60 V 40 V 30 V 20 V 0 V20 V 40 V 60 V 80 V

2,200

2,510

2,000

0 150

30 T s1 60 T s1 120 T s1 180 T s1 240 T s1 300 T s1

2,500

1,500

80

60

40

20

1,500

2,490

750

375

1,000

250

2,480

125

500

2,470

0

80 60 40 20 0 20 40 60 80

40

60

80

100

50 0 50 100

150

H|| (mT)

H (mT)

H|| (mT)

b

R S ()

2,500

Figure 3 | SIT in the magnetoresistance. (a) Parallel eld MR for the fastest sweep rate (300 T s 1) as a function of gate voltage, tuningthe system from the superconducting to the insulating regime. The top and bottom panels are for insulating and superconducting regimes, respectively. Chargevortex duality manifests itself as the conversion of the peak in the superconducting regime to a dip in the insulating regime. The peak and dip occur at the external eld value of H||~ 12 mT. (There is an axis break on the y axis in the bottom panel.) The maximum change in the resistance occurs for Vg = 0 V, where the superconductivity isvery weak. (b) Similar behaviour in perpendicular eld. Note that the fastest sweep rate in this case is 600 T s 1. Additional structure inthe perpendicular eld at H~ 48 mT is due to the more complex magnetization dynamics in the perpendicular eld. T = 50 mK.

2,540

2,520

60 T s1 120 T s1 240 T s1 360 T s1 480 T s1 600 T s1

2,480

2,460

2,440

150 100 50 0 50 100 150

H (mT)

Figure 2 | MR in the insulating regime. (a) Parallel eld MR forVg = 100 V. Data for only one eld sweep direction are shown for clarity.

Arrow indicates the direction of eld sweep. (b) Perpendicular eld MR for Vg = 100 V. T = 50 mK.

(or dips in resistance) on the insulating side. This is exactly what is observed.

Gate voltage dependence of the peaks/dips. As Vg is changed from + 80 to 100 V, tuning the system through the SIT, the peaks in the MR change to dips, as shown in Fig. 3a,b. A more striking graphical demonstration of the SIT can be seen if we plot the rate dependence of the resistance at the peaks or dips at dierent gate voltages normalized to their values at the highest sweep rate (Fig. 4a,b). The transition from peaks to dips in Figs 3 and 4 occurs at exactly the same Vg as the transition from the superconducting to the insulating state in the R versus T dependence8. Thus the same phenomenon, that of magnetization dynamics in the ferromagnet, which causes an increase in resistance on the superconducting side, causes an increase in conductance on the insulating side. This interchange of resistance and conductance between the two dierent ground states of the system demostrates the duality of the two states. We note, however, that the resistance per square, RS, at which the transition occurs is approximately 2.1 k, less than the universal resistance value of h/4e2~6.45 k expected for the SIT20.

The peaks in the MR on the superconducting side are due to the magnetic eld associated with a moving perpendicular component of the magnetization in the ferromagnet as the external magnetic eld is swept. Given the fact that the position of the peaks and dips

coincide for both H and H||, it is clear that the dips also have the same origin. How would this eld give rise to an increase in conductance on the insulating side? If the non-superconducting regime was simply a weakly localized metal, it would be hard to explain the dips in the MR that we observe (for a discussion, see Supplementary methods, which also contains alternate explanations of the MR dips, and why they do not account for our observations). These dips can be explained if one models the insulating side of the transition as one where the Cooper pairs are localized on isolated superconducting islands, and a Cooper pair requires an energy Ec to transfer between two neighbouring islands, even though Ec might be very small. Classically, the time-dependent magnetization of the ferro-magnet gives rise to a time-dependent magnetic eld, which in turn gives rise to an electric eld H = A/t through Faradays law; here

A is the vector potential associated with the eld. The voltage dierence V = Hd generated between two neighbouring superconducting islands separated by a distance d can exceed the charging energy

Ec, V > Ec/2e, so that the conductance of the system is momentarily increased. This explains the dips in MR in the insulating regime

(Fig. 3a,b) when the domain wall passes through the ferromagnet.

Discussion

As we saw above, modelling our system as a granular superconductor explains the MR data qualitatively in both the superconducting and the insulating regimes. A more quantitative understanding of the peaks and dips can be obtained if one considers a thermal activation model for vortex motion in the superconducting regime and for charge transport in the insulating regime. As before, we begin with a discussion of the MR in the superconducting state.

On the superconducting side of the transition, there is strong Josephson coupling between the islands (EJ Ec), and in the absence of motion of magnetic vortices, the sample is in the zero resistance state. Any magnetic vortices present exist in the interstitial regions between the islands. Motion of the magnetic vortices,

NATURE COMMUNICATIONS | 3:955 | DOI: 10.1038/ncomms1959 | www.nature.com/naturecommunications

2012 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms1959

ARTICLE

a

b

R d/R d at 300 T s 1

R p/R p at 300 T s 1

R d/R d at 600 T s 1

R p/R p at 600 T s 1

where the attempt frequency depends on microscopic parameters as well as external parameters of the circuit. Here F is the free energy barrier for phase jumps in the two directions26

F E I

e

2 2 .

J

p

1.02

1.02

0

1.015

1.015

100 V 80 V 60 V 40 V 30 V 20 V 0 V20 V 40 V 60 V 80 V

1.01

1.01

1.005

1.005

(2)(2)

(3)(3)

(4)(4)

(5)(5)

1

1

0.8

0.8

0.6

0.6

0.4

0.4

0.2

0.2

0

0 2 4 6 8 10 0 2 4 6 8 10

In the LA model, the application of a voltage V drives the generation of phase slips, and in steady state, the average voltage due to the generation of phase slips is equal to the applied voltage, giving the relation

V e e I ek T

E k T

= (2 / ) ( /2 )

2 /

p p

V J B B

sinh

|| Sweep rate ( 30 T s1) Sweep rate ( 60 T s1)

Figure 4 | Sweep rate dependence. (a) The evolution of the normalized resistance of the peaks and dips at H||~12 mT, with external eld sweep rate as a function of the gate voltage for parallel eld. The top and bottom panels are for insulating and superconducting regimes, respectively. (b) Similar plot for perpendicular eld. The resistance of the peaks quickly saturates in the superconducting regime due to the additional pair breaking effect of the external eld, unlike in the parallel eld case. T = 50 mK.

In the low current limit (I ek T

B /) the corresponding

resistance is

R V I e k T e E k T

= / = ( / ) .

2 2 2 /

p B J B

V

which in our case is induced by the magnetization dynamics in the ferromagnet, will give rise to a nite resistance as the vortices cross the junctions between the superconducting islands. To understand the mechanism of the generation of this resistance, we use the phase slip model due to Langer and Ambegaokar (LA)24. There are corrections to this model, for example by McCumber and Halperin25, but the basic concepts are all that are essential for our analysis.

For simplicity, we consider two adjacent superconducting islands that are Josephson coupled to each other at two points, so that the interstitial region between them can enclose a magnetic vortex: the sample consists of a network of such interconnected islands. (More generally, one probably has multiple interconnected islands enclosing vortices, but this does not change the physical picture.) These two islands can be thought of as a dc SQUID, whose energy prole is given by U() = EJcos, where is the phase dierence across the SQUID and EJ is the Josephson energy. EJ is proportional to the critical current Ic of the SQUID. This denes a periodic potential where the minima of the potential correspond to integral values of the phase = 2n. At nite temperature, the system can be thermally activated over the energy barrier represented by EJ from one potential minimum to an adjacent minimum which diers by a value = 2. Each such phase slip event will give rise to a voltage pulse according to the Josephson relation V e d dt

= ( /2 )( / )

f , where the phase change is 2. In the absence of an external current, phase jumps in either direction are equally likely, hence the average voltage measured across the junction is zero, although such phase slip events may give rise to voltage noise.

In the presence of an external measuring current I, the SQUIDs energy prole is modied to the tilted washboard potential U E e I

( ) = ( /2 )

In our case, a phase slip is generated whenever a magnetic vortex line crosses a weak link. The rate at which these phase slips are generated is then proportional to the time dependence of the eld generated by the magnetization of the ferromagnet. At the eld sweep rates in our experiment, the magnetization change is essentially quasi-static, hence this rate is proportional to the external magnetic eld sweep rate

B, resulting in an additional factor in equation (4) proportional to

B. The change in resistance due to the motion of vortices for the entire sample is an average of terms such as equation (4) with a distribution of EJs, but will be proportional to the magnetic eld sweep rate

B. On the other hand, the measured value of the critical current Ic of the sample will be determined by a parallel combination of random paths through the weak links connecting the superconducting islands. Nevertheless, it is reasonable to assume that this measured critical current Ic is proportional to the mean Josephson energy EJ for the superconducting network.

Thus, we can write the change in resistance due to the motion of magnetic vortices in the form

R ABe I k T

= /

a c B

where Ic is the measured critical current, and A and are numerical constants at a xed temperature. Thus, if one plots ln(R) as a function of Ic at a xed temperature, one should obtain a straight line. Ic in this system can be varied by changing the gate voltage

Vg. Figure 5a shows this plot. As can be seen, the expected exponential dependence is indeed observed for a number of dierent magnetic eld sweep rates. The slopes (equal to the exponential factor /kBT) of the linear ts in Fig. 5a are 2.7107, 2.22107 and 1.97107 A 1, for the three sweep rates, 180, 240 and 300 T s 1, respectively. The slopes (which should be equal for the three sweep rates) match to within about 25%.

Equation (5) also predicts that the resistance change should be proportional to the magnetic eld sweep rate

B. If one looks at the parallel eld data for Vg = 0 and 20 V shown in Fig. 4a, one can see that the dependence on sweep rate is indeed approximately linear.

However, for Vg = 40, 60 and 80 V, which are deeper in the superconducting regime, the dependence is clearly not linear. In fact, deep in the superconducting regime, one might expect very little change in resistance, as the exponential factor involving Ic in equation (5)

would suppress any resistance change. Thus, the data at the three slowest sweep rates are not surprising. However, at higher sweep rates, the peak amplitude for these gate voltage values does increase.

f f f

J cos . In this case, phase slip events corresponding to the system travelling down the washboard tilt are slightly more likely than those in the other direction, leading to a nite average voltage, and hence a nite resistance. However, the probability of such events is still quite small at low temperatures if I Ic.

LAs original paper focused on phase slips in a single weak link between two superconductors. They derived the rate at which phase slips occur in each direction

h

= /

Ve F k T

B (1)(1)

NATURE COMMUNICATIONS | 3:955 | DOI: 10.1038/ncomms1959 | www.nature.com/naturecommunications

2012 Macmillan Publishers Limited. All rights reserved.

ARTICLE

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms1959

a

9

b

0

0.16

8

7

100

0.15

Log (R)

6

R ()

200

5

180 T s1

240 T s1

300 T s1

0.14

4

300

0.13

3 1e07 2e07 3e07

0 2 4 6 8 10

|| Sweep rate ( 30 T s1)

Figure 5 | Thermal activation model for the peaks/dips. (a) Log(R) versus Ic in the superconducting regime for the three fastest sweep rates in parallel external eld. R is dened as the difference in resistance between the forward and reverse magnetic eld sweep traces at H||~12 mT. The ve gate voltages corresponding to the ve Ic are Vg = 0, 20, 40, 60 and 80 V from the smallest to the largest Ic. The three lines are the ts to the dataat the three different sweep rates. (b) R as a function of the parallel eld sweep rate in the insulating regime. The data for each gate voltage,

Vg = 30, 40, 60, 80 and 100 V, are tted to an exponentially decaying function. Fitting parameters are given in the Discussion section. T = 50 mK.

T c (K)

Ic (A)

0.12

0.11

0.1

600 400 200 0 200 400 600

H|| (mT)

Figure 6 | Parallel critical eld: Tc versus H|| for Vg = 80 V. A parabolic t (black curve) at high elds gives Hc||~1.2 T. The thickness of the superconducting layer estimated from this value is d = 13.6 nm.

We do not know the reason for this, although one possibility is that at higher sweep rates, the eective Josephson coupling between islands is modied. We note that the data for the perpendicular eld shown in Fig. 4b do not show a linear behaviour, but in this case, the situation is complicated by the fact that there is also a orbital contribution to the peak height due to the external perpendicular magnetic eld.

On the insulating side of the transition, a thermal activation model for charge transport can be developed in analogy with the vortex hopping model on the superconducting side. In this model, the rate at which charge can be transferred between islands is proportional to e E k T

c B

/ (ref. 27), and as with the superconducting regime above, charge transfers in both directions are equally likely, so there is no net average current. In the presence of a nite voltage dierence V between the two islands, the rate is modied to e E e V k T

( 2 )/

c B

,

and consequently favours charge transfer from the island at higher potential to the island at lower potential, resulting in a net average current between the islands. In our model, the potential V results from the electric eld that is generated by the moving magnetic eld generated by the magnetization dynamics in the ferromagnet through Faradays law. The time dependence of the changing magnetization, as we have pointed out above, is proportional to the sweep rate of the external magnetic eld

B, and hence V =

B,

where is a constant. Thus, the resistance of the dip should have a dependence on the magnetic sweep rate of the form

R A Be B k T

= .

/

As Vg is changed to bias the system closer to the transition from the insulating side, it is reasonable to expect that the eective spacing between islands decreases, resulting in a decrease in , as is observed in Fig. 5b.

In summary, the LAO/STO interface, with its unusual combination of superconductivity and ferromagnetism, provides a unique signature of chargevortex duality in the SIT. It would be interesting to see if similar signatures can be observed in articially fabricated hybrid ferromagnetsuperconductor structures.

Methods

Experimental techniques. The samples in this work had 10 unit cells (uc) of LaAlO3 grown by pulsed laser deposition on TiO2-terminated (001) SrTiO3 single-crystal substrates13,14. The electrical measurements were performed on a Hall bar dened by photolithography and etched using argon ion milling. Details of the lm preparation, characterization and sample fabrication have been discussed in detail in previous publications8,13,14. The samples were measured in an Oxford dilution refrigerator with a base temperature of 15 mK. This refrigerator was equipped with a two-axis magnet so that a magnetic eld could be applied both perpendicular and parallel to the LAO/STO interface. A gate voltage Vg was applied to the back of the 500-m-thick STO substrate.

Transport measurements were made by a standard ac lock-in detection method. The samples were current biased with an excitation current of Iac = 10 nA at a frequency of 11.3 Hz. For Vg = 20 V and higher, this excitation current is

much less than the critical current, Ic (Iac < 0.1 Ic). The voltage signals were rst amplied using a low-noise homemade instrumentation amplier and then

measured using a PAR 124 analogue lock-in amplier or an EG&G 7260 digital lock-in amplier.

Determination of lm thickness and parallel critical eld. From the continuous mapping of Tc versus H one can estimate the superconducting coherence length, , as discussed in ref. 8. We obtained a value for the coherence length of ~70 nm at Vg = 80 V. To determine the thickness of the superconducting layer, we performed a similar measurement of the dependence of Tc on H||. The system was biased at the mid-point of the resistive transition (RS = 288 at Vg = 80 V), under a proportional-integral-dierential (PID) feedback circuit. The PID output was used to drive the mixing chamber heater while the eld was ramped. This ensured that we always stayed at Tc, dened as the mid-point of the resistive transition.

A plot of Tc versus H|| measured in this way is shown in Fig. 6. As expected for a two-dimensional superconductor, the dependence of Tc on a parallel applied eldis quadratic, except at low elds where hysteresis due to the ferromagnet suppresses

Tc. From a t to this measurement at high elds, we extracted the parallel critical eld, Hc||~1.2 T. For a thin lm superconductor, Hc|| and the lm thickness are related by the formula28

H d

c o

+

B

(6)(6)

We note that this is a characteristically dierent dependence than that which is predicted and found in the superconducting regime.

Figure 5b shows a plot of the resistance change at the dip as a function of the sweep rate for dierent gate voltages on the insulating side: the lines are ts to the exponential form, equation (6). As can be seen, the ts are quite good. The t parameters for Vg = 30, 40, 60, 80 and 100 V are, respectively: A = 568.59, 332.77, 371.9, 381.05 and 373.59 ; B = 578.56, 360.572, 379.10, 424.60 and 391.51 ; and kBT/ = 25.11, 5.22, 5.97, 5.79 and 5.89 T s 1. For Vg = 100, 80, 60 and 40 V, the value of is almost the same. For Vg = 30 V, it is quite dierent. The time-dependent magnetic eld generates an electric eld; the resulting voltage V is proportional to the the electric eld, but also to the spacing a between the islands. Hence should also depend on a.

= 3F

p x

where o = h/2e is the superconducting ux quantum.

From a knowledge of and Hc||, we determined the thickness of the superconducting lm to be d = 13.6 nm, which is less than .

(7)(7)

NATURE COMMUNICATIONS | 3:955 | DOI: 10.1038/ncomms1959 | www.nature.com/naturecommunications

2012 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms1959

ARTICLE

a

2,400

400

2,200

300

R S ()

R S ()

R S ()

2,000

500 mT

250 mT

100 mT

20 mT

9 mT

0 mT

200

100 V 60 V 45 V 30 V

1,800

100

1 0.5 0 0.5 1

Idc (A)

0 150 100 50 0 50 100 150

Figure 8 | Insulating state IV characteristics; RS versus Idc for different gate voltages at T = 15 mK.

H (mT)

b

2,370

2,360

2,350

2,340

2,330

500 mT 0 mT

2,320

2,310

2,300

100 50 0 50 100

H (mT)

Figure 7 | Crossed parallel and perpendicular elds. (a) RS versus H for different values of parallel elds at Vg = 100 V. For clarity, eld sweeps in only one direction are shown, with the arrow indicating the directionof eld sweep. (b) RS versus H for H|| = 0 and 500 mT at Vg = 40 V.

For both curves, the eld is swept at the same rate of 120 T s 1. ForH|| = 500 mT, characteristic features due to the ferromagnet are absent.

T = 50 mK for both the curves.

1,600

MR in crossed parallel and perpendicular elds. As our system shows hysteresis due to the magnetization dynamics of the ferromagnet, it becomes inherently more complex to study the eect of an externally applied eld to the system. To eliminate the hysteresis we applied a persistent parallel eld, H||, which saturated the magnetization of the ferromagnet in one direction, and then performed perpendicular eld MR measurements. In Fig. 7a, we show data in which perpendicular MR is measured for various values of H|| for Vg = 100 V at T = 50 mK. It can be seen that for H|| = 20 mT and higher, the sharp resistance peak in the MR is absent and the magnetization of the ferromagnet is completely saturated. For H|| = 0 mT, the sharp peak is recovered. Note that the minimum of resistance and the eld at which this minimum occurs both increase as the parallel eld is increased. The additional pair breaking caused due to a constant parallel eld gives rise to the small increase in resistance. The shiing of the zero of the magnetic eld is due to a small misalignment of plane of the sample with respect to the axis of our split coil magnet. From the shi we calculate the misalignment to be 1.4.

In the insulating phase (for Vg = 40 V, say) the generation of ux due to domain wall motion is responsible for a rise in conductance. As discussed above, a changing vector potential,

A t

/ , generates an electric eld, H, which causes a potential, V, to develop between two adjacent islands. When this potential energy 2 eV exceeds the charging energy of the islands, Ec, one gets conduction between the islands. To see that the dips in resistance that we see on the insulating side are indeed due to such a mechanism, we show R versus H traces for Vg = 40 V, for

H|| = 0 and 500 mT in Fig. 7b. Both the traces in Fig. 7b are at the same sweep rate (120 T s 1). By applying a parallel eld of 500 mT, we have eliminated the eect of the magnetization dynamics in the ferromagnet. The sharp dips in resistance vanish completely, and only a small positive MR due to weak localization remains.

The shi in zero of the H|| = 500 mT curve is due to the small misalignment of the sample as discussed above.

Insulating state IV characteristics. The IV characteristics of the system were discussed in ref. 8. The IV curves in the insulating state were almost linear with no clear evidence of a Coulomb gap, indicating a weak insulating state. However, if dV/dI (RS) is plotted as a function of an applied bias current Idc, a peak in dV/dI is observed at low bias. This is shown in Fig. 8 for Vg = 30, 45, 60 and 100 V at T = 15 mK. For Vg = 30 and 45 V, a dip in dV/dI at zero bias signies the presence of some superconductivity; however, no evidence for that is seen in the temperature dependence of the resistance, as dR/dT < 0 at all temperatures in the temperature range of the current experiments. The dip in dV/dI probably arises due to localized superconducting islands weakly phase coupled by the Josephson eect. For Vg = 60 and 100 V, the system is in the insulating regime with a

peak in dV/dI at zero bias.

References

1. Haviland, D. B., Liu, Y. & Goldman, A. M. Onset of superconductivity in the two-dimensional limit. Phys.Rev.Lett. 62, 21802183 (1989).

2. Hebard, A. F. & Paalanen, M. A. Magnetic-eld-tuned superconductor-insulator transition in two-dimensional lms. Phys.Rev.Lett. 65, 927930 (1990).

3. Geerligs, L. J., Peters, M., de Groot, L. E. M., Verbruggen, A. & Mooij, J. E. Charging eects and quantum coherence in regular Josephson junction arrays. Phys.Rev.Lett. 63, 326329 (1989).

4. Ohtomo, A. & Hwang, H. Y. A high-mobility electron gas at the LaAlO3/SrTiO3 heterointerface. Nature 427, 423426 (2004).

5. Thiel, S., Hammerl, G., Schmehl, A., Schneider, C. W. & Mannhart, J. Tunable quasi-two-dimensional electron gases in oxide heterostructures. Science 313, 19421945 (2006).

6. Reyren, N. etal. Superconducting interfaces between insulating oxides. Science 317, 11961199 (2007).

7. Caviglia, A. D. etal. Electric eld control of the LaAlO3/SrTiO3 interface ground state. Nature 456, 624627 (2008).

8. Dikin, D. A. etal. Coexistence of superconductivity and ferromagnetism in two dimensions. Phys.Rev.Lett. 107, 056802 (2011).

9. Schneider, T., Caviglia, A. D., Gariglio, S., Reyren, N. & Triscone, J.- M. Electrostatically-tuned superconductor-metal-insulator quantum transition at the LaAlO3/SrTiO3 interface. Phys.Rev.B 79, 184502 (2009).

10. Bert, J. A. etal. Direct imaging of the coexistence of ferromagnetism and superconductivity at the LaAlO3/SrTiO3 interface. Nat.Phys. 7, 767771 (2011).

11. Li, L., Richter, C., Mannhart, J. & Ashoori, R. C. Coexistence of magnetic order and two-dimensional superconductivity at LaAlO3/SrTiO3 interfaces.

Nat.Phys. 7, 762766 (2011).12. Michaeli, K., Potter, A. C. & Lee, P. A. Superconducting and ferromagnetic phases in SrTiO3/LaAlO3 oxide interface structures: possibility of nite momentum pairing. Phys.Rev.Lett 108, 117003 (2012).

13. Park, J. W. etal. Creation of a two-dimensional electron gas at an oxide interface on silicon. Nat.Comm. 1, 94 (2010).

14. Bark, C. W. etal. Tailoring a two-dimensional electron gas at the LaAlO3/ SrTiO3 (001) interface by epitaxial strain. Proc.NatlAcad.Sci. 108, 47204724 (2011).

15. Kittel, C. IntroductiontoSolidStatePhysics, 7th edn (Wiley, 1995).

NATURE COMMUNICATIONS | 3:955 | DOI: 10.1038/ncomms1959 | www.nature.com/naturecommunications

2012 Macmillan Publishers Limited. All rights reserved.

ARTICLE

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms1959

16. Jaeger, H. M., Haviland, D. B., Orr, B. G. & Goldman, A. M. Onset of super-conductivity in ultrathin granular metal lms. Phys.Rev.B 40, 182196 (1989).

17. Kim, J. B., Akinaga, H. & Kim, J. Direct observation of the spin congurations of vertical Bloch line. Appl.Phys.Lett. 98, 052510 (2011).

18. Fazio, R. & Schon, G. Charge and vortex dynamics in arrays of tunnel junctions. Phys.Rev.B 43, 53075320 (1991).

19. Girvin, S. M. Duality in Perspective. Science 274, 524525 (1996).20. Fisher, M. P. A., Grinstein, G. & Girvin, S. M. Presence of quantum diusion in two dimensions: universal resistance at the superconductor-insulator transition. Phys.Rev.Lett. 64, 587590 (1990).

21. Fisher, M. P. A. Quantum phase transitions in disordered two-dimensional superconductors. Phys.Rev.Lett. 65, 923926 (1990).

22. van der Zant, H. S. J., Fritschy, F. C., Elion, W. J., Geerligs, L. J. & Mooij, J. E. Field-induced superconductor-to-insulator transitions in Josephson-junction arrays. Phys.Rev.Lett. 69, 29712974 (1992).

23. Rimberg, A. J. etal. Dissipation-driven superconductor-insulator transition in a two-dimensional Josephson-junction array. Phys.Rev.Lett. 78, 26322635 (1997).

24. Langer, J. S. & Ambegaokar, V. Intrinsic resistive transition in narrow superconducting channels. Phys.Rev. 164, 498510 (1967).

25. McCumber, D. E. & Halperin, B. I. Time scale of intrinsic resistive uctuations in thin superconducting wires. Phys.Rev.B 1, 10541070 (1970).

26. Halperin, B. I., Refael, G. & Demler, E. Resistance in superconductors. Preprint at arXiv:1005.3347v1 (2010).

27. Shimshoni, E., Auerbach, A. & Kapitulnik, A. Transport through quantum melts. Phys.Rev.Lett. 80, 33523355 (1998).

28. Tinkham, M. IntroductiontoSuperconductivity, 2nd edn (McGraw-Hill, 1996).

Acknowledgements

We thank P.A. Lee, A.M. Goldman and J.B. Ketterson for providing useful comments. Work at Northwestern was supported by a grant from the DOE Office of Basic Energy Sciences under grant no. DE-FG02-06ER46346. Work at the University of Wisconsin was supported by funding from the DOE Office of Basic Energy Sciences under award number DE-FG02-06ER46327 and the National Science Foundation under grant no. DMR-0906443.

Author contributions

M.M.M., D.A.D. and V.C. performed the experiments and analysed the data. C.W.B., S.R., C.M.F. and C.B.E. grew the samples. M.M.M. and V.C. wrote the paper. All authors contributed to the concept of the paper.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/ naturecommunications

Competing nancial interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/ reprintsandpermissions/

How to cite this article: Mehta, M.M. etal. Evidence for chargevortex duality at the LaAlO3/SrTiO3 interface. Nat.Commun. 3:955 doi: 10.1038/ncomms1959

(2012).

NATURE COMMUNICATIONS | 3:955 | DOI: 10.1038/ncomms1959 | www.nature.com/naturecommunications

2012 Macmillan Publishers Limited. All rights reserved.