COMMUNICATIONS

ARTICLE

Received 11 Aug 2015 | Accepted 3 Dec 2015 | Published 27 Jan 2016 | Updated 26 Feb 2016^^BDOI:iai038/ncomms10386l OPEN

Polaronic metal state at the LaAlO3/SrTiO3 interface

C. Cancellieri1'2'*, A.S. Mishchenko3'*, U. Aschauer4'5, A. Filippetti6, C. Faber4, O.S. Barisic7, V.A. Rogalev1, T. Schmitt1, N. Nagaosa3,8 & V.N. Strocov1*

Interplay of spin, charge, orbital and lattice degrees of freedom in oxide heterostructures results in a plethora of fascinating properties, which can be exploited in new generations of electronic devices with enhanced functionalities. The paradigm example is the interface between the two band insulators LaAlO3 and SrTiO3 that hosts a two-dimensional electron system. Apart from the mobile charge carriers, this system exhibits a range of intriguing properties such as field effect, superconductivity and ferromagnetism, whose fundamental origins are still debated. Here we use soft-X-ray angle-resolved photoelectron spectroscopy to penetrate through the LaAlO3 overlayer and access charge carriers at the buried interface. The experimental spectral function directly identifies the interface charge carriers as large polarons, emerging from coupling of charge and lattice degrees of freedom, and involving two phonons of different energy and thermal activity. This phenomenon fundamentally limits the carrier mobility and explains its puzzling drop at high temperatures.

1 Swiss Light Source, Paul Scherrer Institute, Villigen CH-5232, Switzerland. 2EMPA, Swiss Federal Laboratories for Materials Science & Technology,

Ueberlandstrasse 129, Duebendorf 8600, Switzerland. 3 RIKEN Center for Emergent Matter Science, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan.

4 Materials Theory, ETH Zurich, Wolfgang-Pauli-Strasse 27, Zurich CH-8093, Switzerland. 5 Department of Chemistry and Biochemistry, University of Bern,

Freiestrasse 3 3012 Bern, Switzerland. 6CNR-IOM, Istituto Officina dei Materiali, Cittadella Universitaria, Cagliari, Monserrato 09042-I, Italy. 7 Institute of

Physics, Bijenicka 46, 10000 Zagreb, Croatia. 8 Department of Applied Physics, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8656, Japan.

* These authors contributed equally to this work. Correspondence and requests for materials should be addressed to V.N.S. (email: vladimir.strocov@psi.ch).

Coupling of the electron and lattice degrees of freedom in solids through electron-phonon interaction (EPI) is a key concept in electron transport and many other phenomena of condensed matter physics. An electron moving in the lattice can displace atoms from their equilibrium positions in response to the EPI. Such an electron (or hole) dragging behind a local lattice distortion—or phonon 'cloud'—forms a composite charge carrier known as polaron1'2. The increased effective mass m* of this quasi-particle (QP) fundamentally limits its mobility mp1/m* beyond incoherent scattering processes. The polarons are key players in many technological devices' a widespread example of which are high-electron mobility transistors (HEMTs) utilized in high-frequency devices such as mobile phones. In typical HEMTs, a donor layer of n-doped AlGaAs injects electrons into the channel layer of intrinsic GaAs where, escaping scattering on the dopant impurities, the electrons are limited in their mobility only by the polaronic coupling enhanced by spatial confinement in the GaAs quantum well (QW)3.

Angle-resolved photoelectron spectroscopy (ARPES) is a unique method to measure the single particle spectral function A(O'k) in crystalline solids resolved in electron energy o and momentum k. Containing all many-body (electron-electron, electron-phonon, and so on) interactions, A(o,k) reveals the formation of polarons by a characteristic peak-dip-hump (PDH) lineshape, where the sharp peak corresponds to a QP and the broad hump, extending to higher binding energies, corresponds to the cloud of entangled phonons with frequencies o0 (refs 1,2). However, the extreme surface sensitivity of conventional ARPES with photon energies hv below ~100eV sets the buried interfaces out of its reach. The crucial feature of our experiment is the use of soft-X-ray ARPES (SX-ARPES) operating in the hv range of hundreds eV (for a recent review see ref. 4). The longer photoelectron mean free path enables SX-ARPES to penetrate through the top layers and access A(o,k) at buried interfaces.

Complex oxide interfaces are presently at the forefront of fundamental research in view of their enhanced functionalities achieved by exploiting electron correlations5,6. The two-dimensional electron system (2DES) in LaAlO3/SrTiO3 (LAO/ STO)7 is confined within a narrow region of a few nanometres on the STO side8-10, where the mobile electrons populate the t2g-derived dxy-, dxz- and dyz-states of Ti ions acquiring reduced valence compared with bulk Ti4 +. Confinement in the interface QW further splits these states into a ladder of subbands8,11-14. This complex energy structure based on the correlated 3d orbitals, very different from conventional semiconductor heterostructures described as free particles embedded in the mean-field potential, is the source of a rich and non-trivial phenomenology. Here high 2DES mobility m2DES typical of uncorrelated ESs co-exists with superconductivity5,6, ferromagnetism15,16, large magneto-resistance17 and other phenomena typical of localized correlated electrons. Other intriguing puzzles in this intricate physics are why m2DES measured in transport falls short of estimates based on mean-field theories, and what causes the dramatic drop of m2DES with increase of temperature18.

Here, we directly access the nature of the LAO/STO interface carriers through their A(o,k) measured by SX-ARPES at ultrahigh-energy resolution. We discover that the LAO/STO interface forms a polaronic metal state involving at least two active phonons. While polaronic coupling to hard LO3 phonons fundamentally limits m2DES at low temperatures, coupling to soft TO1 phonons with increasing temperature provides the microscopic mechanism of m2DES drop observed in transport.

Results

Spectroscopic signatures of the polaronic metal state. Our

LAO/STO(001) samples with an LAO overlayer thickness of b 18 A corresponding to 5 unit cells (u.c.) were grown using Pulsed Layer Deposition (PLD), and subsequently annealed in oxygen atmosphere to minimize the concentration of oxygen vacancies (VOs) and the related extrinsic charge carriers7,19. SX-ARPES with its crucial advantage of enhanced probing depth is ideally suited to study this buried system where the 2DES only develops with the LAO layer thickness > 3 u.c. (ref. 20) The extremely small 2DES signal, however, has to be boosted using resonant photoemission21 locked to the interface Ti ions14,22,23. For details of the sample growth and SX-ARPES experiment, see Methods.

Our low-temperature experimental data set in Fig. 1 was measured at 12 K using s-polarized X-rays (the parallel p-polarization data are given in the Supplementary Note 1). The resonance map of (angle-integrated) photoemission intensity, Fig. 1a, was recorded under variation of hv across the Ti 2p absorption edge ~460eV. We identify there the 2DES signal at the Fermi level EF blowing up near the two Ti3 + L3- and L2-resonances and vanishing everywhere else (we notice that strong suppression of the in-gap states around — 1.2 eV, related to the VOs (refs 14,22,23), confirms the prevalence of the intrinsic interface charge carriers). Tuning hv onto the stronger L3-resonance produces the Fermi surface (FS) map in Fig. 1b where, by comparison with the superimposed theoretical FS contours, we recognize the manifold of merged circular dxy-derived FS sheets and the elliptical dyz-sheet extending in the kx-direction14,23.

The ARPES images measured along the TX (ky = 0) line of the square 2D Brillouin zone (BZ) at the L3- and L2-resonances are shown in Fig. 2a,b, respectively. With a high-energy resolution of 40 meV, these images resolve individual interface bands. The use of s-polarization selects the dxy- and dyz-derived states, which are antisymmetric relative to the TX line14 (although the selection rules are slightly relaxed by remnant structural distortions). By comparison with the overlaid E(k) dispersions calculated for LAO/STO interface with pseudo-self-interaction correction DFT (see Methods) one can recognize the lower dyz-band with its flat dispersion. Already at this point, we note signs of a QP interaction, which reduces the band dispersion compared with the overlaid DFT prediction as characterized by an effective mass ratio of m*/m0~2.5. The dxy-bands are not visible due to vanishing matrix elements, but the lowest dxy-band appears as two bright spots where it hybridizes with the dyz-band.

The most striking visual aspect of the experimental E(k) is, however, the two vertical waterfalls extending down from these

E-EF(eV) -6.0 -4.0 -2.0 0.0

kx (Â-1)

- 0.4-0.2 0.0 0.2 0.4

Int.max

468 , 465 ' 462 459 456

Int.min

Figure 1 | Experimental low-temperature (12 K) overview SX-ARPES data. (a) Resonance (angle-integrated) photoemission intensity map, identifying the 2DES signal at the L3 and L2 resonances of the interface Ti ions. (b) FS map at the L3-resonance, where the superimposed theoretical FS contours identify the dxy (pink) and dyz (green) sheets. The colour bar indicates the intensity (Int.) maximum and minimum.

high-intensity dxy spots. In Fig. 2c,d we show the energy distribution curves (EDCs)—that is, ARPES intensity as a function of binding energy for a given k—extracted from the images in Fig. 2a,b, respectively. These EDCs reveal, remarkably, a pronounced PDH structure of A(o,k), where the peak reflects the QP and the hump at ~ 118 meV below the peak its coupling to bosonic modes (such as magnon, plasmon, phonon) whose nature will be identified later on.

We note in passing that the EDC representing the whole dxy-band, Fig. 2c, exhibits a smaller but broader QP peak in comparison to its dyz-counterpart in Fig. 2d. Their nearly equal integral QP weight indicates that, non-trivially, the bosonic coupling is quite insensitive to different spatial distribution of the dxy- and dyz-states8. The larger broadening of dxy-EDC can reflect larger defect scattering, because the lowest dxy-state is localized closest to the interface where the concentration of defects generated by the non-equilibrium PLD growth is maximal, whereas the dyz-state extends deeper into the defect-free STO bulk. This observation is consistent with the recent Shubnikov-de

kx (Â-1)

kx (Â-1)

-0.4-0.2 0.0 0.2 0.4

-0.4 -0.2 0.0 0.2 0.4

Int.max

Int.min

-X r X^

0.3 Â-1

0 Â-1 -0.2 Â-1

-0.6 -0.4 -0.2 0.0 E-Ef (eV)

-0.6 -0.4 -0.2 0.0 E-Ef (eV)

Haas experiments15, which have found smaller dxy-mobility compared with dxz/yz-one.

We will now address the nature of the involved bosonic modes. Considering the extremely small ferromagnetic response of the LAO/STO interface15, the magnons can safely be ruled out. Plasmons can also be excluded since the energy of the hump does not depend on the interfacial carrier concentration ns varied via manipulation of VOs (see Supplementary Note 2) while the plasma frequency op is proportional to y/n~s. These bosonic modes can therefore be assigned to phonons coupling to electron excitations and forming polarons2 6. The hump apex, located at B118meV below the QP peak, identifies the main coupling phonon frequency oo0. The 2DES at the LAO/STO interface realizes therefore a polaronic metal state.

To identify the phonon modes forming the observed polaron, we used DFT to calculate the phonon dispersions (see Methods) for cubic bulk STO at different electron doping concentrations nv; the results are shown in Fig. 3 a. Although these initial calculations are for cubic STO, we will address shifts of the phonon frequencies at the cubic to tetragonal phase transition when discussing the temperature dependence. The response of the calculated dispersions quantifies the EPI strength as a function of phonon mode and q. With the actual electron density distribution8-10, our ns measured by the Luttinger area of the experimental FS roughly corresponds to nv = 0.12 electrons per u.c. Since the carriers reside on the STO side of the LAO/STO interface, we expect these results to be relevant for the interface as well. The polaron can be associated with the hard longitudinal optical phonon LO3, which is the only mode available in this high-energy range and, moreover, has the largest coupling constant l among all LO modes27'28 as observed by Raman29 and neutron spectroscopy30. At finite q vectors, the LO3 mode

q-vector

Figure 2 | Experimental low-temperature (12 K) high-resolution

SX-ARPES. (a,b) High-resolution ARPES images along the TX (ky = 0) line at the L3- and L2-edges, respectively, with the superimposed theoretical dxy (pink) and dyz (green) bands. The lower panels show the corresponding second derivative -d2//dE2>0 plots, which clearly show both the quasi-particle (QP) peak and the dispersive hump formed by the LO3 phonon. (c,d) A series of (normalized) EDCs extracted from a,b respectively, at the indicated kx-values through the occupied part of the BZ. The colour bar indicates the intensity maximum and minimum. The two curves at the bottom show EDCs integrated over the k-ranges indicated in c,d as well as the whole BZ in the kx-direction. The characteristic PDH spectral structure in c,d manifests a polaronic metal state formed by the hard LO3 phonon and renormalizing the dyz-band dispersion in a,b and the clear hump dispersion in b identifies a large polaron.

Figure 3 | Theoretical phonon modes in cubic phase doped STO.

(a) Phonon dispersion at various electron doping levels nv, with our case corresponding to nv~0.12. The arrows indicate the TO1 and TO2 modes shifting as a function of nv. The imaginary modes at the R- and M-points represent different octahedral rotation instabilities, whereas the one at the T-point in the undoped materials is the polar (quantum-paraelectric) instability. (b,c) Atomic displacements associated with the breathing LO3 mode at the R-point and the polar TO1 mode at the T-point, respectively.

represents a breathing distortion of the octahedral cage around a Ti site, Fig. 3b, which is typical of polaron formation driven by the (Holstein-type) short-range EPI, barely sensitive to electron screening in a metallic system. This character of EPI is confirmed by the increase of with increase of q away from the T-point. Recent optical studiesV on bulk STO31 have not only confirmed the involvement of the LO3 phonon in EPI but also found the corresponding effective mass renormalization m*/m0 B3.0 close to our value. This mode has also been observed by ARPES on (doped) bare STO(001)32-34 with the strength of the polaronic structure depending on ns. We note that, our theoretical LO3 energy of B100 meV perfectly matches that found experimentally for the bare STO bulk and surface, but it differs from our o0'~118meV measurement at the LAO/STO interface. To gain some insight on this discrepancy, we investigated the possible role of strain by performing additional calculations with the STO inplane lattice constants constrained to that of LAO. However, only a minor frequency shift to 102 meV was found. This implies that other interfacial factors such as the electric field, phonon coupling across the interface, propagation of the LAO distortion into STO35, structural changes due to electrostatic doping36 or weaker coupling to additional phonon modes could also play a role. We also note that the EPI at LAO/STO is enhanced by the tight 2D electron confinement in the interface QW2^, which was recently evidenced by huge oscillations of thermopower as a consequence of large phonon drag37.

Next, we perform a theoretical analysis of the PDH structure to estimate the strength of the EPI governing the polaron formation. With the total A(o,k) spectral weight normalized to unity, we define Z0 as the integral weight of the QP part. The observed A(o,k) of a polaron with momentum k can be expressed as a sum

of two terms1,2

A(o, k) = Z0d[o - E(k)} + AH (a, k)

The first term represents the sharp QP peak with the dispersion E(k) and effective mass m* resulting from the renormalization of the non-interacting single-particle band with dispersion e(k) and mass m0. The second term arises due to phonons coupling to the excited photohole. The EPI has a twofold effect on A(o,k): it reduces Z0 below unity, and builds up AH(o, k), with its dispersion following e(k), as a Frank-Condon series of phonon peak replicas at energy separations on = no0 from the QP peak, where o0 is the phonon frequency and n indexes the replicas.

For our analysis, we used the angle-integrated EDC at the L3-resonance (Fig. 2c), which is dominated by the dxy-intensity and thus insensitive to the dispersion effects in the dyz-band (certain admixture of the latter is not important because its dispersion range is much smaller than our experimental resolution of 40 meV and, moreover, the dxy- and dyz-bands have the same peak-to-hump relative weight). Gaussian fitting of its QP peak yields Z0~0.4. For the short-range EPI, forming the LO3 polaron, exact diagrammatic quantum Monte Carlo (QMC) calculations38 have shown that the relation m*/m0~1/Z0 holds with high accuracy. With the determined Z0, this yields m*/m0 = 2.5, which coincides with the value extracted above from the band dispersions, leaving not much room for electron correlations to contribute to the band renormalization. The latter is confirmed by theoretical analysis of tunnelling experiments39 where the local electron correlations were described with the Hubbard parameter U. According to (ref. 40), already at U = 4t (t is the hopping integral) the density of states significantly renormalizes and forms a pseudogap at EF. The absence of the latter in the tunnelling spectra39 indicates U<4t. Then the renormalization expe1cted from the correlation effects m*/m0 = [1 — (U/Uc)2] — (Uc = 8t is the critical Mott value) is < 1.3, which is much smaller than our experimental m*/m0.

Hence, the mass renormalization at the LAO/STO interface is dominated by EPI.

Furthermore, the coupling strength l for our short-range EPI can be estimated as l ~ 1 — Z0, which yields l e 0.6. This value is well-below the limit A~0.9 separating the weak- and strong-coupling regimes, which implies that we observe large polarons, where the lattice distortion extends over several unit cells1,2,41. This is perfectly consistent with the clear dispersion of the hump, which tracks the non-interacting e(k) dispersion of the dyz-band in Fig. 2d. Importantly, the small l also excludes that self-trapping of small polarons42,43 can be responsible for the discrepancy between the observed mobile charge and the 0.5 electrons per u.c. required for full compensation of the polar field in the LAO overlayer5,6,43. However, the EPI can assist charge trapping on shallow defects44 created by non-equilibrium PLD growth or VOs, in addition to the deep level trapping.

The polaronic reduction of m2DES fundamentally limits the application potential of the STO-based heterostructures. As illustrated in Supplementary Note 2, this limit can possibly be circumvented through manipulation of VOs. The VOs inject into the 2DES extrinsic charge carriers43, which increase the electron screening and thus reduce the EPI. This trend is consistent with the recent study at bare STO(100) surfaces33 and, furthermore, explains the recent results on the g-Al2O3/STO45, where the VOs have not only increased ns but also m2DES by almost 2 orders of magnitude. On the other hand, the VOs can assist the EPI due to charge trapping on shallow defects44 and also increase the defect scattering rate, both effects counteracting the above positive trend. Further experiments on oxygen-deficient LAO/STO will allow a better understanding of the role of VOs and ways to optimize m2DES.

Temperature dependence of polaronic effects. An intriguing peculiarity of the charge carriers at the LAO/STO interface is the drop of m2DES by more than 1 order of magnitude as the temperature increases above 200 K (ref. 18). To unveil the underlying microscopic mechanism, we measured the temperature dependence of the L3 angle-integrated EDC from Fig. 2c (the angle integration makes our analysis robust against extrinsic thermal scattering of high-energy photoelectrons, which averages the ARPES signal in k-space46). These results are shown in Fig. 4a in comparison with experimental temperature dependence of m2DES in Fig. 4b, characteristic of the LAO/STO samples18, which was derived from Hall effect measurements. The experimental A(o,k) in Fig. 4a shows that with increase of temperature, first, the polaronic hump broadens and propagates to higher binding energies. This trend is accompanied by increase of the QP peak width, Fig. 4d. These effects can be attributed to incoherent scattering on thermally populated phonon modes. A striking effect beyond this weight-conserving mechanism is, however, that the QP peak loses its integral weight Z0 (Fig. 4c) and completely dissolves in the phonon hump of A(o,k) towards 190 K (Fig. 4a). The decrease of Z0, following the trend of the temperature-dependent m2DES, unambiguously signals that the QP spectral weight leaks to soft phonon modes, different from the hard LO3, with concomitant further increase of m*. This phenomenon provides the microscopic mechanism of the puzzling m2DES drop with increasing temperature18,41. Interestingly, ARPES on the TiO2(001) surface shows similar temperature effects24 despite different crystallographic structure of TiO2. In line with our results, an optical study of bulk STO31 has also revealed reduction of the Drude weight with temperature. Surprisingly, the bare STO(100) surface at large nv (ref. 32) does not show any systematic

........ 1 1 1 T=

190 K 170 K 150 K

130 K 110 K

— 100 K

—'—~/ 70 K

30 K 12 K

1 . 1 . 1 . 1 . i i i

-0.6 -0 E-

4 -0.2 0.0 0.2 -EF(eV)

Figure 4 | Temperature dependence of polaronic effects. (a) Angle-integrated EDCs at the L3-resonance acquired for temperatures between 12 and 190 K. The QP peak dissolving into the hump towards ~190 K explains (b) the mobility drop with temperature measured in transport (the dashed line trends through the experimental points). (c) Temperature-dependent QP spectral weight Z0(T) fitted by the independent boson model of equation (2). The fit identifies a soft phonon mode (likely the TO1) with o0" = 18 meV (solid line) in the low-Tregion, which shifts to o0" = 14meV (dashed) in the high-T region through the tetragonal to cubic phase transition in STO. The fading colours represent the validity ranges of the fits. (d) QP peak width (including the instrumental resolution). The error bars in c,d correspond to the uncertainty in the Gaussian fit of the QP.

temperature effects in the QP weight; clearly, the LAO overlayer and 2DES significantly alter the EPI in the LAO/STO system.

Numerical analysis of the temperature-dependent ARPES spectral shape allows us to identify the frequency o000 of the soft phonon mode dominating the EPI. From the normalized EDCs, Fig. 4a, we evaluate the temperature dependence of the spectral weight Z0(T) by Gaussian fitting of the QP peak. The resulting Z0(T) in Fig. 4c can then be fitted with the analytic formula for the independent boson model1

Zo (T ) = e " 2g(2N +1) /0 [2g (2N + 1)],

where N = (effl°=T — 1) 1 is the Bose filling factor, which describes the population number as a function of the mode frequency o0 at T, and J0 is the modified Bessel function. This equation, neglecting momentum dependence and thus valid for our momentum-integrated EDCs, describes transfer of the spectral weight from the QP peak to the hump with temperature. The constant g-factor, measuring the EPI strength, was set to 0.95 to reproduce the experimental value Z0(12K)b0.4. Fitting Z0(T) in the temperature range below b120K yields o0" e 18 meV (solid line in Fig. 4c), while the high-T range yields o0" ~ 14 meV (dashed line in Fig. 4c). This crossover of o000 can be linked to the second-order tetragonal to cubic phase transition in STO at 105K (refs 32,47). One may therefore argue that the m2DES drop in LAO/STO has structural origin. Interestingly, the fastest change of the QP peak width in Fig. 4d falls into the same temperature region. All other possible temperature effects in the QP linewidth, including the Bloch-Gruneisen contribution due to acoustic phonons, are weaker for our system dominated by the optical phonon, and obscured by the phase transition and limited energy resolution. We note that, strictly speaking, our Z0(T) model (2) had to include EPI with

both low- and high-energy phonons, contributing to the total EPI. However, the result would not change significantly because the Bose filling factor N, determining the Z0(T) dependence, is most sensitive to small o0. In other words, the Z0(T) dependence is most sensitive to the low-energy sector of the phonon spectrum.

Next we return to our DFT calculations to identify the observed soft o000 phonon mode. The phonon dispersions of the cubic (high-T) phase, Fig. 3a, show a multitude of low-energy modes. The lowest energy LO mode (LO1), which was previously linked with kinks in ARPES dispersions28, has an energy of b 22 meV. This energy is considerably higher than our high-T fitted o000B14 meV and, moreover, our calculations do not show a significant frequency difference between the cubic and tetragonal phases of STO. These facts, together with the reported small coupling constant l (ref. 27), make the involvement of this LO1 mode unlikely. In Fig. 3 a we see, however, the TO modes strongly affected by the electron doping nv due to enhanced screening. The TO1 mode, sketched in Fig. 3c, is a polar mode whose instability in undoped STO leads to its quantum-paraelectric behaviour48,49. Its frequency rapidly increases and becomes real with increase of nv, stabilizing above our actual nv~0.12 electrons per u.c. at oTO1~ 13.7 meV at the T point. This oTO1 matches well with our o000B14 meV fitted in the high-T range. Turning to the low-T range, previous calculations on undoped STO have shown that the tetragonal phase transition increases oTO1 (ref. 50). Our computations with nv = 0.12 electrons per u.c. for the tetragonal phase indeed show oTO1 to shift to 15.3 meV for the doubly degenerate mode and 18.1 meV for the non-degenerate mode along the octahedral rotation axis, which is in good agreement with our low-T fitted value o000B18 meV. Octahedral rotations different from this bulk tetragonal distortion were observed at the LAO/STO interface due to electrostatic doping35,36, but since they are combinations of distortions similar to the one studied here, they are expected to also lead to an upwards shift of the TO1 frequency. Based on the t1u symmetry of the TO1 mode, its good agreement with the experimental o000 and its sensitivity to the cubic to tetragonal phase transition, we associate the experimental soft phonon with TO1. Moreover, significant polaronic coupling to this polar mode is consistent with gigantic dielectric constant e0 of STO caused by large polar ionic displacements under electric field in this material on the verge of a ferroelectric instability30. The TO1 mode is associated with long-range EPI, as evidenced by the increase of towards the T-point. The hard LO3 and soft TO1 modes inVvolved in the polaron formation are therefore associated with opposite types of EPI. The EPI strength in the latter case can be estimated from the Frohlich model which, although strictly valid only for the polar LO modes, represents the only viable model to find the long-range coupling constant a. Based on exact diagrammatic QMC calculations51 and Z0 at the top of our temperature range, we estimate aE2 for the TO1 mode. Similar to the above LO3 phonon, this value stays below the weak- to strong-coupling crossover at aE6, which implies that the TO1 phonon is also consistent with the large polaron scenario. We note that if LAO/STO superconductivity is driven by a phonon mechanism52, it can be related to the discovered polaronic activity. While the hard LO3 phonon energy much exceeds the energy scale of the superconducting transition at 0.3 K, involved in the electron pairing may be the soft TO1 phonon.

The 2DES at the LAO/STO interface realizes thus a polaronic metal state involving at least two phonons with different energies and thermal activity. The hard LO3 phonon at o00 ~ 118 meV, associated with short-range EPI, is directly resolved as the characteristic hump in the experimental A(o,k). It sets the fundamental limit of m2DES at low temperature, and exhaustively

accounts for the m* renormalization without notable effects of electron correlations. Another soft phonon, likely the TO1 one associated with long-range EPI, changes its frequency from O0/0~18 to 14meV across the phase transition in STO. This phonon causes a dramatic fading of the QP weight with temperature as the microscopic mechanism behind the m2DES drop observed in transport. The two phonons form a large polaron characterized by a lattice deformation extending over several unit cells. Physically, the involvement of both LO3 and TO1 phonons in the polaronic state means that the corresponding lattice distortion combines both breathing and polar atomic displacements. The polaronic activity is typical of oxide perovskites, reflecting their highly ionic character and easy structural transformations30,42. Our discovery may have implications for other related oxide systems, including LAO/STO interfaces with different crystallographic orientations. In a methodological perspective, we have demonstrated the power of the newly emerging experimental technique of ultrahighresolution SX-ARPES to retrieve information about polaronic effects at buried interfaces in the most direct way as embedded in one-electron A(o,k).

Methods

Sample preparation. Our LAO/STO samples were grown using PLD (for details of the growth procedure see refs 19,22) and subsequently annealed in oxygen atmosphere with a pressure of 200mBar at 550 0C for 1 h to minimize the concentration of VOs. These vacancies manifest themselves in ARPES spectra as characteristic dispersionless in-gap states at a binding energy of B1.2 eV (refs 14,22) which, similar to bare STO, grow with exposure to X-rays. Our resonant spectra in Fig. 1a show only traces of such spectral structures enhanced at the Ti3 + L3- and L2-resonances, indicating negligible concentration of Vqs. Measurements on samples with various ns varied through the VOs (Supplementary Note 2) demonstrated constant energy of the 118-meV spectral peak, excluding its plasmonic origin. The temperature-dependent mobility data were derived from standard Hall effect and conductivity measurements performed in the van der Paw geometry on square samples.

SX-ARPES experiment. SX-ARPES experiments were performed at the ADRESS beamline of the Swiss Light Source, Paul Scherrer Institute, Switzerland53. The experimental geometry allows symmetry analysis of the valence states using variable linear polarizations of incident X-rays. The experiments are normally performed at low sample temperatures B12 K to quench thermal scattering of high-energy photoelectron destructive for the coherent k-resolved spectral component46. The combined (beamline and analyser) energy resolution was set to 80 meV for measurements of the FS, and to 40 meV for high-resolution measurements of the band dispersions. Such resolution achieved for an interface buried behind a B18-A-thick overlayer presents the forefront of nowadays SX-ARPES instrumentation. The temperature dependence was measured with increasing temperature to avoid possible hysteresis effects54.

First-principles calculations. Band structure calculations for the LAO/STO interface were performed using the pseudo-self-interaction correction method55, with a plane wave basis set and ultrasoft pseudopotentials. This ab initio approach corrects the main deficiencies of basic density functional theory for a vast range of oxides8. The theoretical band structure and the FS were calculated for an interfacial carrier density of 0.115 electrons per u.c. This value is consistent with those determined using Hall effect measurements for the LAO/STO interface. The corresponding Fermi momentum kF of the dyz-states is in good agreement with the experimental value of B -0.29 A — 1 determined as the highest band dispersion point.

DFT-based phonon calculations under different electron doping were performed using the VASP code57. The electron count was adjusted while adding a compensating background charge. We used the PBEsol functional, which gives reliable lattice parameters and phonon frequencies56, and PAW potentials57 with Sr^s^p^s), Ti(3p,3d,4s) and O(2s,2p) valence shells. Phonons were then computed using the frozen phonon approach58. For further details, see ref. 50.

References

1. Mahan, G. D. Many Particle Physics (Plenum, 1981).

2. Holstein, T. Studies of polaron motion. Ann. Phys. (N.Y.) 8, 343-389 (1959).

3. Das Sarma, S. & Stopa, M. Phonon renormalization effects in quantum wells. Phys. Rev. B 36, 9595-9603 (1987).

4. Strocov, V. N. et al. Soft-X-ray ARPES at the Swiss Light Source: from 3D materials to buried interfaces and impurities. Synchr. Rad. News 27, 31-40 (2014).

5. Hwang, H. Y. et al. Emergent phenomena at oxide interfaces. Nat. Mater. 11, 103-113 (2012).

6. Mannhart, J. & Schlom, D. G. Oxide interfaces—an opportunity for electronics. Science 327, 1607-1611 (2010).

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

8. Delugas, P. et al. Spontaneous 2-dimensional carrier confinement at the n-Type SrTiO3/LaAlO3 Interface. Phys. Rev. Lett. 106, 166807 (2011).

9. Son, W.-J., Cho, E., Lee, B., Lee, J. & Han, S. Density and spatial distribution of charge carriers in the intrinsic n-type LaAlO3-SrTiO3 interface. Phys. Rev. B 79, 245411 (2009).

10. Dubroka, A. et al. Dynamical response and confinement of the electrons at the LaAlO3/SrTiO3 interface. Phys. Rev. Lett. 104, 156807 (2008).

11. Salluzzo, M. et al. Orbital reconstruction and the two-dimensional electron gas at the LaAlO3/SrTiO3 interface. Phys. Rev. Lett. 102, 166804 (2014).

12. Santander-Syro, A. F. et al. Two-dimensional electron gas with universal subbands at the surface of SrTiO3. Nature 469, 189-193 (2011).

13. King, P. D. C. et al. Quasiparticle dynamics and spin-orbital texture of the SrTiO3 two-dimensional electron gas. Nat. Commun. 5, 3414 (2014).

14. Cancellieri, C. et al. Doping-dependent band structure of LaAlO3/SrTiO3 interfaces by soft x-ray polarization-controlled resonant angle-resolved photoemission. Phys. Rev. B 89, R121412 (2014).

15. Fete, A. et al. Large modulation of the Shubnikov-de Haas oscillations by the Rashba interaction at the LaAlO3/SrTiO3 interface. New J. Phys. 16, 112002 (2014).

16. Banerjee, S., Erten, O. & Randeria, M. Ferromagnetic exchange, spin-orbit coupling and spiral magnetism at the LaAlO3/SrTiO3 interface. Nat. Phys. 9, 626-630 (2013).

17. Caviglia, A. D. et al. Tunable Rashba spin-orbit interaction at oxide interfaces. Phys. Rev. Lett. 104, 126803 (2010).

18. Gariglio, S., Reyren, N., Caviglia, A. D. & Triscone, J.-M. Superconductivity at the LaAlO3/SrTiO3 interface. J. Phys. Condens. Matter 21, 164213 (2009).

19. Cancellieri, C. et al. Influence of the growth conditions on the LaAlO3/SrTiO3 interface electronic properties. Europhys. Lett. 91, 17004 (2010).

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

21. Molodtsov, S. L. et al. Angle-resolved resonant photoemission as a probe of spatial localization and character of electron states. Phys. Rev. Lett. 78, 142 (1997).

22. Cancellieri, C. et al. Interface Fermi states of LaAlO3/SrTiO3 and related heterostructures. Phys. Rev. Lett. 110, 137601 (2013).

23. Berner, G. et al. Direct k-space mapping of the electronic structure in an oxideoxide interface. Phys. Rev. Lett. 110, 247601 (2013).

24. Moser, S. et al. Tunable polaronic conduction in anatase TiO2. Phys. Rev. Lett. 110, 196403 (2013).

25. Mishchenko, A. S. & Nagaosa., N. Electron-phonon coupling and a polaron in the t — J model: from the weak to the strong coupling regime. Phys. Rev. Lett. 93, 036402 (2004).

26. Mishchenko, A. S. et al. Polaronic metal in lightly doped high-Tc cuprates. Europhys. Lett. 95, 57007 (2011).

27. Verbist, G., Peeters, F. M. & Devreese, J. T. Extended stability region for large bipolarons through interaction with multiple phonon branches. Ferroelectrics 130, 27-34 (1992).

28. Meevasana, W. et al. Strong energy-momentum dispersion of phonon-dressed carriers in the lightly doped band insulator SrTiO3. New J. Phys. 12, 023004 (2010).

29. Cardona, M. Optical properties and band structure of SrTiO3 and BaTiO3. Phys. Rev. A 140, 651 (1965).

30. Choudhury, N., Walter, E. J., Kolesnikov, A. I. & Loong, C.-K. Large phonon band gap in SrTiO3 and the vibrational signatures of ferroelectricity in ATiO3 perovskites: First-principles lattice dynamics and inelastic neutron scattering. Phys. Rev. B 77, 134111 (2008).

31. van Mechelen, J. L. M. et al. Electron-phonon interaction and charge carrier mass enhancement in SrTiO3. Phys. Rev. Lett. 100, 226403 (2010).

32. Chang, Y. J., Bostwick, A., Kim, Y. S., Horn, K. & Rotenberg, E. Structure and correlation effects in semiconducting SrTiO3. Phys. Rev. B 81, 235109 (2010).

33. Wang, Z. et al. Tailoring the nature and strength of electron-phonon interactions in the SrTiO3(001) two-dimensional electron liquid. Preprint at http://arxiv.org/abs/1506.01191 (2015).

34. Chen, C., Avila, J., Frantzeskakis, E., Levy, A. & Asensio, M. C. Observation of a two-dimensional liquid of Fröhlich polarons at the bare SrTiO3 surface. Nat. Commun 6, 8585 (2015).

35. Jia, C. L. et al. Oxygen octahedron reconstruction in the SrTiO3/LaAlO3 heterointerfaces investigated using aberration-corrected ultrahigh-resolution transmission electron microscopy. Phys. Rev. B 79, 081405(R) (2009).

36. Doennig, D. & Pentcheva, R. Control of orbital reconstruction in (LaAlO3)M/(SrTiO3)N(001) quantum wells by strain and confinement. Sci. Rep. 20, 7909 (2015).

37. Pallecchi, I. et al. Giant oscillating thermopower at oxide interfaces. Nat. Commun. 6, 6678 (2015).

38. Mishchenko, A. S., Nagaosa, N. & Prokofev, N. Diagrammatic Monte Carlo method for many-polaron problems. Phys. Rev. Lett. 113, 166402 (2014).

39. Breitschaft, M. et al. Two-dimensional electron liquid state at LaAlO3-SrTiO3 interfaces. Phys. Rev. B 81, 153414 (2010).

40. Feng, G. S. & White, S. R. Numerical study of the spectral weight function for the half-filled two-dimensional Hubbard model. Phys. Rev. B 46, 8691(R) (1992).

41. Mishchenko, A. S., Nagaosa, N., De Filippis, G., de Candia, A. & Cataudella, V. Mobility of Holstein polaron at finite temperature: an unbiased approach. Phys. Rev. Lett. 114, 146401 (2015).

42. Mannella, N. et al. Polaron coherence condensation as the mechanism for colossal magnetoresistance in layered manganites. Phys. Rev. B 76, 233102 (2007).

43. Hao, X., Wang, Z., Schmid, M., Diebold, U. & Franchini, C. Coexistence of trapped and free excess electrons in SrTiO3. Phys. Rev. B 91, 085204 (2015).

44. Mishchenko, A. S. et al. Localization-delocalization transition of a polaron near an impurity. Phys. Rev. B 79, 180301(R) (2009).

45. Chen, Y. Z. et al. A high-mobility two-dimensional electron gas at

the spinel/perovskite interface of g-Al2O3/SrTiO3. Nat. Commun. 4, 1371 (2013).

46. Braun, J. et al. Exploring the XPS limit in soft and hard x-ray angle-resolved photoemission using a temperature-dependent one-step theory. Phys. Rev. B 88, 205409 (2013).

47. Wang, F. F. Y. & Gupta, K. P. Phase transformation in the oxides. Metall. Trans. 4, 2767 (1973).

48. Muller, K. A. & Burkard., H. SrTiO3: an intrinsic quantum paraelectric below 4K. Phys. Rev. B 19, 3593 (1979).

49. Vanderbilt, D. & Zhong, W. First-principles theory of structural phase transitions for perovskites: competing instabilities. FerroeZectrics 206, 181-204 (1998).

50. Aschauer, U. & Spaldin, N. A. Competition and cooperation between antiferrodistortive and ferroelectric instabilities in the model perovskite SrTiO3. J. Phys. Condens. Matter 26, 122203 (2014).

51. Mishchenko, A. S., Prokofev, N. V., Sakamoto, A. & Svistunov, B. V. Diagrammatic quantum Monte Carlo study of the Fröhlich polaron. Phys. Rev. B 62, 6317 (2000).

52. Boschker, H., Richter, C., Fillis-Tsirakis, E., Schneider, C. W. & Mannhart., J. Electron-phonon coupling and the superconducting phase diagram of the LaAlO3-SrTiO3 interface. Sci. Rep. 5, 12309 (2015).

53. Strocov, V. N. et aZ. Soft-X-ray ARPES facility at the ADRESS beamline of the SLS: Concepts, technical realisation and scientific applications. J. Synchrotron Rad. 21, 32-44 (2014).

54. Schoofs, F., Egilmez, M., Fix, T., MacManus-Driscoll, J. L. & Blamire, M. G. Impact of structural transitions on electron transport at LaAlO3/SrTiO3 heterointerfaces. Appl. Phys. Lett. 100, 081601 (2012).

55. Filippetti, A. et aZ. Introducing a variational pseudo-self-interaction correction approach for solids and molecules. Phys. Rev. B 84, 195127 (2011).

56. Wahl, R., Vogtenhuber, D. & Kresse, G. SrTiO3 and BaTiO3 revisited using the projector augmented wave method: Performance of hybrid and semilocal functionals. Phys. Rev. B 78, 104116 (2008).

57. Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B 59, 1758-1775 (1999).

58. Togo, A., Oba, F. & Tanaka, I. First-principles calculations of the ferroelastic transition between rutile-type and CaCl2-type SiO2 at high pressures. Phys. Rev. B 78, 134106 (2008).

Acknowledgements

We thank P.R. Willmott, J.H. Dil, R. Claessen, J.-M. Triscone, F. Baumberger, F. Bisti, Z. Wang, M. Radovic, J. Mesot and N.A. Spaldin for fruitful discussions, and X. Wang for software support of the experiment. Parts of this research were supported by the ImPACT Program of the Council for Science, Technology and Innovation (Cabinet Office, Government of Japan), the Croatian QuantiXLie Center of Excellence, and the ETH Zurich and ERC Advanced Grant program, No. 291151.

Author contributions

C.C. and V.N.S. have performed the experiment supported by V.A.R. and T.S. and processed the data. C.C. has grown the samples. V.N.S. and A.S.M. have developed the scientific concept. A.S.M. has performed the polaronic analysis supported by O.S.B. and N.N. A.F. and U.A. have calculated the electron and phonon band structures, respectively. All authors have discussed the results and interpretations, as well as the manuscript written by V.N.S. and A.S.M. with contributions of C.F., U.A. and C.C.

Additional information

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

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

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

How to cite this article: Cancellieri, C. et al. Polaronic metal state at the LaAlO3/SrTiO3 interface. Nat. Commun. 7:10386 doi: 10.1038/ncomms10386 (2016).

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

COMMUNICATIONS

DOI: 10.1038/ncomms10920

Corrigendum: Polaronic metal state at the LaAlO3/SrTiO3 interface

C. Cancellieri, A.S. Mishchenko, U. Aschauer, A. Filippetti, C. Faber, O.S. Barisic, V.A. Rogalev, T. Schmitt, N. Nagaosa & V.N. Strocov

Nature Communications 7:10386 doi: 10.1038/ncomms10386 (2016); Published 27 Jan 2016; Updated 26 Feb 2016

In the original PDF version of this Article, which was published on 27 January 2016, the publication date was incorrectly given as 27 January 2015. This has now been corrected in the PDF; the HTML version of the paper was correct from the time of publication.

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/