Jocienne N Nelson1, Nathaniel J Schreiber2, Alexandru B Georgescu3,4, Berit H Goodge5, Brendan D Faeth1, Christopher T Parzyck1, Cyrus Zeledon2, Lena F Kourkoutis5,6, Andrew J Millis3,7, Antoine Georges3,8,9,10, Darrell G Schlom2,6,11, Kyle M Shen1,6. 1. Laboratory of Atomic and Solid State Physics, Department of Physics, Cornell University, Ithaca, NY 14853, USA. 2. Department of Materials Science and Engineering, Cornell University, Ithaca, NY 14853, USA. 3. Center for Computational Quantum Physics, Flatiron Institute, New York, NY 10010, USA. 4. Department of Materials Science and Engineering, Northwestern University, Evanston, IL 60208, USA. 5. School of Applied and Engineering Physics, Cornell University, Ithaca, NY 14853, USA. 6. Kavli Institute at Cornell for Nanoscale Science, Ithaca, NY 14853, USA. 7. Department of Physics, Columbia University, New York, NY 10027, USA. 8. Collège de France, 11 place Marcelin Berthelot, 75005 Paris, France. 9. CPHT, CNRS, Ecole Polytechnique, IP Paris, F-91128 Palaiseau, France. 10. DQMP, Universitè de Genéve, 24 quai Ernest Ansermet, CH-1211 Genéve, Suisse. 11. Leibniz-Institut für Kristallzüchtung, Max-Born-Str. 2, 12489 Berlin, Germany.
Abstract
Interface quantum materials have yielded a plethora of previously unknown phenomena, including unconventional superconductivity, topological phases, and possible Majorana fermions. Typically, such states are detected at the interface between two insulating constituents by electrical transport, but whether either material is conducting, transport techniques become insensitive to interfacial properties. To overcome these limitations, we use angle-resolved photoemission spectroscopy and molecular beam epitaxy to reveal the electronic structure, charge transfer, doping profile, and carrier effective masses in a layer-by-layer fashion for the interface between the Dirac nodal-line semimetal SrIrO3 and the correlated metallic Weyl ferromagnet SrRuO3. We find that electrons are transferred from the SrIrO3 to SrRuO3, with an estimated screening length of λ = 3.2 ± 0.1 Å. In addition, we find that metallicity is preserved even down to a single SrIrO3 layer, where the dimensionality-driven metal-insulator transition typically observed in SrIrO3 is avoided because of strong hybridization of the Ir and Ru t2g states.
Interface quantum materials have yielded a plethora of previously unknown phenomena, including unconventional superconductivity, topological phases, and possible Majorana fermions. Typically, such states are detected at the interface between two insulating constituents by electrical transport, but whether either material is conducting, transport techniques become insensitive to interfacial properties. To overcome these limitations, we use angle-resolved photoemission spectroscopy and molecular beam epitaxy to reveal the electronic structure, charge transfer, doping profile, and carrier effective masses in a layer-by-layer fashion for the interface between the Dirac nodal-line semimetal SrIrO3 and the correlated metallic Weyl ferromagnet SrRuO3. We find that electrons are transferred from the SrIrO3 to SrRuO3, with an estimated screening length of λ = 3.2 ± 0.1 Å. In addition, we find that metallicity is preserved even down to a single SrIrO3 layer, where the dimensionality-driven metal-insulator transition typically observed in SrIrO3 is avoided because of strong hybridization of the Ir and Ru t2g states.
Ultrathin quantum materials present a unique platform for the control of electronic, magnetic, and topological properties. A commonly observed phenomenon in many ultrathin complex oxides and topological semimetals is that a crossover from a metallic to insulating state occurs below a critical thickness (–). This presents a potential challenge for realizing ultrathin heterostructures of quantum materials where metallic properties are desirable. Here, we show that this thickness-driven metal-insulator transition can be avoided even in the single unit cell limit in certain ultrathin oxide heterostructures.We investigate the interface between SrIrO3 and SrRuO3, which has been of recent interest due to the observation of the topological Hall effect (), suggesting the presence of skyrmions induced by Dzyaloshinskii-Moriya interactions (), and may have potential for spintronic devices () due to the large spin Hall effect in SrIrO3 (). SrRuO3 is a correlated Weyl ferromagnet (), which gives rise to k-space magnetic monopoles in the Berry curvature () and for which a chiral anomaly has been recently reported (). SrIrO3 is a 5d Dirac nodal-line semimetal (, ), formed through the confluence of strong spin-orbit interactions, Coulomb interaction, and octahedral rotations. The topological nature of SrIrO3 and SrRuO3 suggest that their properties will be sensitive to the location of the Fermi energy (EF) with respect to the Dirac nodal line in SrIrO3 or the Weyl points in SrRuO3. Hole doping SrIrO3, in particular, is an attractive target, because the nodal points are only 50 meV below EF but has been challenging to achieve in iridates through chemical substitution (), thus interfacial charge transfer could be a promising alternative. Furthermore, the interplay of magnetism, spin-orbit coupling, electron correlations, and octahedral rotations make accurate theoretical predictions about the SrIrO3/SrRuO3 interface challenging.In this work, heterostructures of n layers of SrIrO3 are synthesized on bulk-like 20–unit cell–thick SrRuO3 films by molecular beam epitaxy (MBE), where n is varied from 1 to 8. Using angle-resolved photoemission spectroscopy (ARPES), we investigate the evolution of the electronic structure of the topmost SrIrO3 layer as it progresses away from the interface. Notably, a bulk-like metallic electronic structure is preserved down to the n = 1 limit, allowing us to isolate the impact of interfacial doping, where we find that approximately 0.5 electrons per Ir are transferred from the interfacial SrIrO3 layer to SrRuO3, shifting the EF of SrIrO3 downward by approximately 80 meV and moving the Dirac node from below to above EF. This charge transfer is confined within the first two unit cells of SrIrO3, with an estimated screening length of λ = 3.2 ± 0.1 Å. Last, the effective mass of the hole-like bands of SrIrO3 are reduced with increased hole doping, revealing that the electronic structure of SrIrO3 does not evolve in a perfectly rigid-band fashion.
RESULTS
Sample synthesis and characterization
SrIrO3/SrRuO3 heterostructures were synthesized by MBE. SrRuO3 films were deposited on LSAT substrates, followed by n layers of SrIrO3, described in Materials and Methods. Scanning transmission electron microscopy (STEM) images of (SrIrO3)/(SrRuO3)20 heterostructures (n = 1 and 2) are shown in Fig. 1 (D and E), demonstrating precise control over the atomic thickness. After growth, samples were transferred through an ultrahigh vacuum manifold (pressure better than 10−10 torr) into an analysis chamber where ARPES, x-ray photoemission spectroscopy (XPS), and low-energy electron diffraction (LEED) measurements were performed. Figure 1A shows a LEED image of an n = 1 sample, demonstrating that the SrIrO3 surface layer exhibits (0.5, 0.5), (±0.5, 0), and (0, ±0.5) diffraction spots (referenced in pseudocubic notation). This indicates that it has the same Pnma space group and octahedral rotation pattern (a−b+a−, Glazer notation) as that of bulk (fig. S4) (, ), although not necessarily with equal magnitude of rotations. In Fig. 1 (B and C), we show valence band and core level (XPS) measurements of the oxygen 2p and Ir 4f states, respectively. The difference in the relative position of the O 2p states has been used to estimate band alignment and charge transfer in oxide perovskite heterostructures (), but here, the positions of the O 2p states in SrIrO3 and SrRuO3 appear to be nearly identical. Furthermore, there is little observable shift in the Ir 4f core levels as a function of the SrIrO3 layer thickness. Therefore, inferring charge transfer between SrRuO3 and SrIrO3 appears to be beyond the capability of conventional core level and valence band spectroscopy used for estimating charge transfer in semiconductors (). This suggests that only direct measurements of the near-EF electronic structure by ARPES can detect charge transfer at the heterointerface.
Fig. 1.
Interface between SrIrO3 and SrRuO3.
(A) LEED measurement of (SrIrO3)1/(SrRuO3)20 (n = 1) with 200-eV electrons at 35 K, with the blue circle highlighting the main cubic peaks, purple circle indicating peaks consistent with R45° octahedral rotations, and yellow circle indicating peaks corresponding to c2 × 2 octahedral rotations. (B) Valence band spectra from bulk-like 24-nm-thick SrRuO3 (yellow) and 8-nm-thick SrIrO3 (blue) films. (C) XPS measurements of the Ir 4f core levels from n = 1 (black) and 3 (yellow), with dashed lines in (B) and (C) as guides to the eye. (D and E) High-angle annular dark-field scanning transmission electron microscopy (HAADF-STEM) images of n = 1 and 2 heterostructures measured along the [100] direction, where p indicates the pseudocubic basis. This shows precise heterostructures with smooth and abrupt interfaces and with the SrIrO3 layers terminated by the expected SrO layer. Samples were capped with either amorphous Al2O3 or crystalline SrTiO3 following ARPES measurements but before removal from UHV.
Interface between SrIrO3 and SrRuO3.
(A) LEED measurement of (SrIrO3)1/(SrRuO3)20 (n = 1) with 200-eV electrons at 35 K, with the blue circle highlighting the main cubic peaks, purple circle indicating peaks consistent with R45° octahedral rotations, and yellow circle indicating peaks corresponding to c2 × 2 octahedral rotations. (B) Valence band spectra from bulk-like 24-nm-thick SrRuO3 (yellow) and 8-nm-thick SrIrO3 (blue) films. (C) XPS measurements of the Ir 4f core levels from n = 1 (black) and 3 (yellow), with dashed lines in (B) and (C) as guides to the eye. (D and E) High-angle annular dark-field scanning transmission electron microscopy (HAADF-STEM) images of n = 1 and 2 heterostructures measured along the [100] direction, where p indicates the pseudocubic basis. This shows precise heterostructures with smooth and abrupt interfaces and with the SrIrO3 layers terminated by the expected SrO layer. Samples were capped with either amorphous Al2O3 or crystalline SrTiO3 following ARPES measurements but before removal from UHV.
ARPES measurements
In Fig. 2, we present the evolution of the electronic structure of (SrIrO3)/(SrRuO3)20 as a function of n, as measured by ARPES. Because of the short photoelectron mean free path (<1 nm) (), these measurements are predominantly sensitive to only the topmost SrIrO3 layer, as evidenced by the fact that no distinct SrRuO3 bands (fig. S5) () can be resolved in the near-EF (EB
<1 eV) region, even for n = 1. The Fermi surface map of the n = 8 sample shows small electron pockets composed of Dirac fermions at (π/2,π/2) and heavy hole pockets at (π,0) (again defined in pseudocubic notation), effectively identical to previous ARPES measurements of thick, bulk-like films of SrIrO3 () and in good agreement with density functional theory (DFT) calculations [k = 0.5π/c, consistent with ()]. We note that the features centered at (0,0) and (π,π) shows weak photoemission intensity likely due to matrix element effects but are present in the ARPES data. With increasing binding energy, the electron pockets shrink and disappear, while the hole pockets expand as determined from intermediate energies and the E versus k dispersion. At EB = 80 meV, because of the weak spectral intensity from the bands that are back-folded because of the bulk in-plane octahedral rotations, this approximates a circular contour. However, as highlighted in white (Fig. 2), the map actually consists purely of square hole pockets. This qualitative behavior is reproduced in the DFT calculations at higher binding energies when taking into account a mass renormalization (m*/mb) of 2, consistent with prior dynamical mean field theory calculations (). This correction is estimated on the basis of the velocity of the linear Dirac dispersion, which forms the electron pockets centered at (π/2,π/2), where vF,DFT ∼2vF. We compare DFT and ARPES using the empirical formula , which allows for reasonable agreement between DFT and ARPES on the size of the electron and hole pockets at each binding energy. Similar features are also observed for n = 3 (also n = 5; see fig. S6), suggesting that an undoped bulk-like electronic structure persists down to 3 SrIrO3 pseudocubic unit cells from the SrRuO3 interface. In contrast, ARPES spectra from n = 1 exhibits a notable departure from the bulk-like n ≥ 3 heterostructures, with a distinctive FS at EF that is nearly identical to the intensity maps of bulk-like n ≥ 3 at a binding energy of EB = 80 meV, indicating a downward shift of the chemical potential, ∆μ, as the SrRuO3 interface is approached with preservation of the bulk-like metallic SrIrO3 state. Notably, no evidence of any gap or suppression of the spectral weight is observed in the n = 1 sample (Fig. 4E), as might be expected for quasi–two-dimensional iridates, such as the Mott-insulating Sr2IrO4. The n = 2 heterostructure exhibits behavior intermediate between n = 1 and 3. This shifting ∆μ with decreasing n is also supported by the similarity of the intensity pattern at even higher binding energies (bottom row), where the intensity pattern at EB = 370 meV for the bulk-like n = 8 sample is closely reproduced at EB = 370 − ∆μ = 290 meV for n = 1. We also note that the extremely strong similarity between the n = 1 and shifted n = 8 maps indicates that the electronic states for the single SrIrO3 layer for n = 1 retains a character very similar to the bulk-like, three-dimensional electronic structure seen in the n = 8 sample. This may appear counterintuitive, because it might be expected that ultrathin SrIrO3 should exhibit a quasi–two-dimensional insulating electronic structure similar to Sr2IrO4. As we will describe later, such quantum confinement is avoided because of hybridization between the SrIrO3 and SrRuO3
t2 states. Such hybridization between the Ru and Ir states is indicated in our orbitally projected DFT calculations of a SrIrO3/SrRuO3 superlattice. For example, orbitally projected FS calculations (fig. S9) clearly show features with clearly mixed Ru and Ir character. In particular, the small “dot” of spectral weight at (0,0) in the n = 1 data appears to be of predominantly Ru character, based on the comparison to the calculations. In addition, some of the spectral weight between (0,0) and (π,π) in the n = 1 sample at EB = 290 meV is not predicted to occur in bulk SrIrO3 and thus should likely arise from hybridization with the SrRuO3. Together, these measurements strongly suggest that the effect of the SrRuO3 is to hole dope the iridate layer and that the interface consists of electronically distinct layers and is not an alloyed system, also suggested by the abrupt interfaces seen by STEM in Fig. 1 (D and E). Because the ARPES measurements are sensitive to only the topmost SrIrO3 layer of the heterostructure, this series of measurements can be thought of as effectively a layer-by-layer tomographic reconstruction of the electronic structure and charge transfer at the SrIrO3-SrRuO3 interface.
Fig. 2.
Constant energy intensity ARPES maps of (SrIrO3)/(SrRuO3)20 heterostructures.
(A) Density functional theory (DFT) calculations of constant energy intensity maps for bulk SrIrO3 at a variety of binding energies, k = 0.5 π/c. The binding energies are renormalized by a factor of 2 to account for many body effects described in the text. The conversion between experiment and DFT is EB = EDFT/2 − 0.035 eV. (B to E) ARPES intensity maps for the n = 8, 3, 2, and 1 heterostructures, with the top row showing the Fermi surface at EB = EF highlighted in red. The square hole pockets are highlighted in white as a guide to the eye in the 80 meV maps of the DFT, n = 8 and the FS of the n = 1, demonstrating that the map consists purely of hole pockets. All data are presented in the pseudocubic Brillouin zone (i.e., 1 Ir per unit cell) and symmetrized about k = 0; unsymmetrized data are shown in fig. S6.
Fig. 4.
Evolution of quasiparticle dispersion and effective mass with n.
(A and B) Dispersion along (0,0)-(π,π) for n = 1 and 8 showing the shifting of the Dirac cone in bulk to above EF for n = 1. (C and D) Dispersion along (0,0)-(0,π) for n = 1 and 8 showing the evolution of the quasiparticle dispersion of the hole-like bands. (E) Energy distribution curves at the Fermi crossings in (A) to (D) as indicated by the black and blue triangles. (F) Extracted dispersion along (0,0)-(0,π) for n = 1, 2, and 8 together with a parabolic fit used to extract m* from the curvature; dispersions are plotted relative to the band maximum. (G) Evolution of m* with n.
Constant energy intensity ARPES maps of (SrIrO3)/(SrRuO3)20 heterostructures.
(A) Density functional theory (DFT) calculations of constant energy intensity maps for bulk SrIrO3 at a variety of binding energies, k = 0.5 π/c. The binding energies are renormalized by a factor of 2 to account for many body effects described in the text. The conversion between experiment and DFT is EB = EDFT/2 − 0.035 eV. (B to E) ARPES intensity maps for the n = 8, 3, 2, and 1 heterostructures, with the top row showing the Fermi surface at EB = EF highlighted in red. The square hole pockets are highlighted in white as a guide to the eye in the 80 meV maps of the DFT, n = 8 and the FS of the n = 1, demonstrating that the map consists purely of hole pockets. All data are presented in the pseudocubic Brillouin zone (i.e., 1 Ir per unit cell) and symmetrized about k = 0; unsymmetrized data are shown in fig. S6.In Fig. 3 we shown the shift of ∆μ as a function of thickness extracted from the ARPES maps in Fig. 2. By fitting to an exponential, we find that the chemical potential shift due to charge transfer between SrIrO3 and SrRuO3 can be well described by the empirical function ∆μ ∝ e−, where d is the thickness of the SrIrO3 layers and the screening length is λ = 3.2 ± 0.1 Å. While this precise number may be affected by nonrigid band effects that we will detail later, this clearly demonstrates the very short length scale of charge transfer in this correlated oxide system, which is relevant for designing ultrathin electronics. The amount of hole doping, ∆h+(n), can also be estimated by shifting the chemical potential in the DFT calculations by the experimentally determined ∆μ (taking into account the mass renormalization factor vF,DFT ∼2vF), shown in Fig. 3A. The rapid decay of ∆μ (and ∆h+) indicates that, above three layers away from the interface, there is no appreciable amount of charge transfer, as shown in Fig. 3A. Here, charge transfer is not used in typical sense of transfer across a buried interface; instead, it describes the estimated doping of the topmost SrIrO3 probed by ARPES. However, it is a closely related quantity, and while the exact ∆μ(n) may change slightly in a buried interface between bulk-like SrIrO3 and SrRuO3, the magnitude and direction of charge transfer should hold. Even the angular distribution of ARPES intensity around a circular path through the FS of n = 1 matches quantitatively to that of n = 8 at EB = 80 meV. This is shown in Fig. 3B, where the distinctive eight-lobed intensity pattern is plotted as a function of polar angle around the FS.
Fig. 3.
Charge transfer and chemical potential shift at SrIrO3-SrRuO3 interface.
(A) Dependence of the shift of the chemical potential (red) and carrier concentration calculated by DFT (blue) on thickness. This illustrates that the n ≥ 3 samples are approximately bulk like. (B) Intensity integrated along a circular path around the FS for the n = 8 sample (red solid lines), which is shifted vertically for clarity; this contrasts with the integrated intensity from the n = 1 FS (black) and the n = 8, EB = 80 meV constant energy map (pink dashed lines), which display a similar angular dependence. The inset shows the Brillouin zone with the angle φ defined.
Charge transfer and chemical potential shift at SrIrO3-SrRuO3 interface.
(A) Dependence of the shift of the chemical potential (red) and carrier concentration calculated by DFT (blue) on thickness. This illustrates that the n ≥ 3 samples are approximately bulk like. (B) Intensity integrated along a circular path around the FS for the n = 8 sample (red solid lines), which is shifted vertically for clarity; this contrasts with the integrated intensity from the n = 1 FS (black) and the n = 8, EB = 80 meV constant energy map (pink dashed lines), which display a similar angular dependence. The inset shows the Brillouin zone with the angle φ defined.In Fig. 4 (A to D), we show ARPES spectra along two high-symmetry directions for the doped n = 1 and undoped n = 8 heterostructures. Energy distribution curves at the Fermi crossings shown in Fig. 4E show strong spectral weight at EF, further demonstrating the persistent metallic nature of the n = 1 heterostructure. At a qualitative level, the change from n = 1 to 8 can be well described by a shift of ∆μ = 80 meV, consistent with Fig. 3. This chemical potential shift pushes the Dirac node at (π/2,π/2) (located 50 meV below EF in bulk) slightly above EF for n = 1. The ability to shift the Dirac node closer to EF without substitutional disorder could be a promising approach to achieve very high mobilities in SrIrO3, as has been demonstrated in its sister compound, CaIrO3 (), where mobilities greater than 60,000 cm2 V−1 s−1 have been reported, ostensibly because of the proximity of the Dirac node to EF (). Nevertheless, a more detailed analysis of the hole pocket clearly shows that the electronic structure does not evolve purely in a simple, rigid-band fashion. In Fig. 4F, the dispersions around (0,π)-(0,0) are plotted as a function of n, along with parabolic fits that are vertically shifted to match the top of the parabolas. This clearly demonstrates that the extracted effective mass, m* (E = E0 − ħ2k2/2m*), decreases by nearly half when going from the bulk-like n = 8 heterostructure (m* = 1.6 ± 0.2me) to the single layer n = 1 heterostructure , summarized in Fig. 4G. We note that the effective mass of this band along this in-plane direction changes only by 6% as a function of k throughout the Brillouin zone.
Evolution of quasiparticle dispersion and effective mass with n.
(A and B) Dispersion along (0,0)-(π,π) for n = 1 and 8 showing the shifting of the Dirac cone in bulk to above EF for n = 1. (C and D) Dispersion along (0,0)-(0,π) for n = 1 and 8 showing the evolution of the quasiparticle dispersion of the hole-like bands. (E) Energy distribution curves at the Fermi crossings in (A) to (D) as indicated by the black and blue triangles. (F) Extracted dispersion along (0,0)-(0,π) for n = 1, 2, and 8 together with a parabolic fit used to extract m* from the curvature; dispersions are plotted relative to the band maximum. (G) Evolution of m* with n.
DISCUSSION
As has now been well established, the rich electronic and magnetic properties of the iridates arise from a complex interplay of spin-orbit coupling, electron correlations, octahedral rotations, and dimensionality. Insulating behavior is generally observed in other quasi–two-dimensional iridates such as Sr2IrO4 and Sr3Ir2O7 (, ) as well as in thin-film heterostructures such as ultrathin SrIrO3 on SrTiO3 (, ) or SrIrO3-SrTiO3 superlattices (, ), raising the obvious question as to why our ultrathin SrIrO3 heterostructure retain a highly metallic character, even for n = 1. The key distinction is that, in the prior examples, the Ir 5d states are confined to two dimensions, either by double SrO rock salt layers or by an insulating SrTiO3 substrate/blocking layers. In contrast, here, the SrRuO3 “substrate” is highly metallic with similar near-EF
t2 states that hybridize with the SrIrO3 electronic states (, ). This prevents quantum confinement, which requires the substrate to have a gap in the available density of states with which the SrIrO3 states could hybridize, either in an insulator with a full bandgap (e.g., SrTiO3), a material with a gap along the film-substrate direction [(e.g., Ag/Au(111) heterostructures ())] or from orbital symmetry considerations (). In this way, although for n = 1 the SrIrO3 layer is only a single unit cell thick from a structural and chemical perspective, its electronic wave function nevertheless retains a three-dimensional bulk-like character due to its extending into the SrRuO3 layers beneath.The hybridization of the Ir 5d states with the Ru t2 states and subsequent preserved three-dimensionality not only explains why a dimensionality-driven metal-insulator transition is avoided but also explains the strong similarity of the ARPES spectra from the n = 1 heterostructure to the thicker, bulk-like samples. For instance, the FS map of the n = 1 sample is nearly identical to the ARPES intensity map from the bulk-like n = 8 sample at EB = 80 meV. However, if the Ir 5d states were confined to two dimensions, then the n = 1 or 2 samples should have a radically different electronic structure than bulk SrIrO3, perhaps more closely resembling doped Sr2IrO4 or Sr3Ir2O7 (, –). Our experimental observations are therefore entirely consistent with a three-dimensional bulk-like electronic structure being preserved even in the n = 1 limit. Furthermore, because the terminating layer is SrIrO3 for all heterostructures, the work function and thus the inner potential and probed k value (k = 0.5 π/c) should be the same across all samples, consistent with our experimental observations.In iridates, it has been established that the octahedral rotations and tilts are also strongly related to the electronic structure through the large spin-orbit coupling. Because SrRuO3 exhibits a similar lattice constant as SrIrO3, the strain is minimized, and the octahedral rotations in SrIrO3 should be much less perturbed than on materials that impart substantial strain (e.g., SrTiO3). This is evidenced by the LEED patterns for all samples remaining unchanged from bulk all the way down to n = 1 (Fig. 1 and fig. S4), suggesting that both bulk SrIrO3 and (SrIrO3)1/(SrRuO3)20 have the same pattern of octahedral rotations as bulk. This is likewise consistent with our DFT calculations for a (SrRuO3)1/(SrIrO3)1 superlattice, which displayed very minor changes in bond angles (<1°) relative to bulk SrIrO3 (note S3). Last, because both SrIrO3 and SrRuO3 are nonpolar, the possibility of a polar reconstruction at the interface (e.g., LaAlO3/SrTiO3) can be eliminated.Therefore, having ruled out quantum confinement/dimensionality, structural modification, and polar discontinuity effects as possibilities, we can then conclude that the observed interfacial charge transfer should be driven by the inherent band alignment between SrIrO3 and SrRuO3. This is consistent with comparisons of the Schottky barrier heights of the SrRuO3/Nb:SrTiO3(001) () and SrIrO3/Nb:SrTiO3(001) (), which indicate a small (≈30 meV) difference in work functions that would imply a transfer of electrons from the Ir to Ru layers. This is also consistent with the aforementioned DFT calculations on a (SrRuO3)1/(SrIrO3)1 superlattice (fig. S8), which also indicate a transfer of electrons from the Ir to the Ru layers. While these differences in work function and charge transfer by DFT are qualitatively consistent with our own conclusions, they are also sufficiently small to underscore the need for high-resolution ARPES to provide a clear picture of the evolution of the electronic structure at the interface.In summary, we have demonstrated a powerful approach for elucidating the evolution of the low-energy interfacial electron structure in a layer-by-layer fashion in correlated metallic heterostructures using ARPES. Our work reveals that electrons are transferred from SrIrO3 to SrRuO3, resulting in hole doping of the closest two SrIrO3 layers to the SrIrO3-SrRuO3 interface, which may have implications to their topological and low-energy properties by pushing the Dirac node of SrIrO3 close to EF. Furthermore, we reveal that ultrathin layers of SrIrO3 on SrRuO3 unexpectedly do not exhibit a crossover to insulating behavior, as has been observed in layered iridates (Sr2IrO4 and Sr3Ir2O7) or ultrathin heterostructures (SrIrO3 on SrTiO3 and SrIrO3-SrTiO3 superlattices), due to the strong hybridization of the Ir and Ru t2 states. The enhanced metallicity (higher carrier density and lower effective mass) of a single unit cell of SrIrO3 on SrRuO3 is notable, given that it runs counter to the typical trend where metallicity is suppressed in ultrathin complex oxides and topological semimetals on insulating substrates [e.g., LaNiO3 (), SrRuO3 (), and Cd3As2 ()]. The enhanced screening for n = 1, potentially arising both from the increased carrier density in the Ir 5d orbitals and the highly metallic SrRuO3, could be responsible for the reduced m*, and possibly a decrease in the effective strength of electron correlations, possibly explaining why the interfacial SrIrO3 layers do not evolve in a perfectly rigid-band fashion. This opens a new avenue for engineering the properties of complex oxides, topological semimetals, and other quantum materials in the ultrathin limit through the judicious selection of the appropriate interface properties and density of states.
MATERIALS AND METHODS
Thin-film growth
Epitaxial SrRuO3 thin films 20 layers (approximately 8 nm) thick were grown using a Veeco GEN10 MBE system on single crystalline LSAT substrates at a substrate temperature of 650°C as measured by an optical pyrometer with a measurement wavelength of 1550 nm in a background partial pressure of 4 × 10−6 torr of distilled ozone (80% O3 + 20% O2). Strontium was evaporated using a low-temperature effusion cell at a flux of 1.5 × 1013 atom cm−2 s−1, while ruthenium was evaporated using an electron beam evaporator. Strontium and ruthenium were codeposited with a Ru/Sr atomic ratio of 1.7-2 starting with an initial nucleation layer of one SrO monolayer before codepositing to avoid formation of ruthenium metal at the interface. Films were then transported in air to a second Veeco GEN10 MBE chamber, where a pristine surface was recovered by annealing the films at 300°C in a background partial pressure of 1 × 10−6 torr of distilled ozone (80% O3 + 20% O2) for 30 min (). Subsequently, the substrate temperature was raised to 660°C (measured by an optical pyrometer with λ = 980 nm) for the growth of the SrIrO3 layer. Strontium was supplied by a low-temperature effusion cell using a flux of 3 × 1012 to 1.3 × 1013 atom cm−2 s−1, and iridium was supplied using an electron beam evaporator. Strontium and iridium were codeposited with an Ir/Sr ratio of 1.2-1.7. All films were capped following measurement before removal from ultrahigh vacuum (UHV) by either amorphous Al2O3 or crystalline SrTiO3. Additional details about the growth and sample characterization may be found in (, ) and note S1.
Photoemission measurements
Following SrIrO3 growth, the samples were immediately transferred through a UHV manifold (pressures better than 10−10 torr) to an analysis chamber, where high-resolution ARPES measurements were performed using a VG Scienta R4000 electron analyzer and a VUV5000 helium plasma discharge lamp and monochromator. He Iα (hν = 21.2 eV) photons with an energy resolution of ∆E = 11 meV at a temperature of 15 K were used unless otherwise noted. XPS measurements were performed immediately after ARPES measurements in the same vacuum chamber using a Scienta Omicron DAR 400 twin anode Mg/Al x-ray source using Mg Kα hν = 1253.6 eV photons with an energy resolution of ∆E = 0.625 eV. LEED measurements were performed using a four-grid LEED optics (SPECS ErLEED 150) in a reflection geometry.
STEM images
Cross-sectional STEM specimens were prepared using the standard focused ion beam (FIB) lift-out process on Thermo Scientific Helios G4 UX FIB or an FEI Strata 400 FIB equipped with an OmniProbe AutoProbe 200 nanomanipulator. High-angle annular dark-field STEM (HAADF-STEM) images were acquired on an aberration-corrected 300 keV FEI Titan Themis with a probe convergence semiangle of 21.4 mrad and inner and outer collection angles of 68 and 340 mrad, respectively.
DFT calculations
The calculations were performed using Quantum ESPRESSO. The structures were relaxed imposing in-plane 1.8% compressive lattice strain along the xy direction relative to the freely relaxed bulk SrIrO3. The final electronic calculations are performed using projector augmented wave pseudopotentials, using the Perdew-Burke-Ernzerhof functional and including spin-orbit coupling, a k-mesh of 10 × 10 × 8, and 816 eV plane wave cutoff. We did not include an additional Hubbard U, as both previous work () and our own has found that it is not necessary to reproduce the main features of the FS of bulk SrIrO3. To obtain the layer-projected FS, we performed a Wannierization of the Ru and Ir t2 states around the Fermi level using Wannier90. Our calculated FSs are then obtained by projecting only the Ir t2 atoms onto the FS. To best simulate the FS, we found a k = 0.5 π/c to be optimal, in agreement with previous work (), and use an energy renormalization of . We found that the (SrRuO3)1/(SrIrO3)1 superlattice adequately reproduces the main features of the FS. To simulate the (SrRuO3)1/(SrIrO3)1 superlattice and bulk SrIrO3, we performed 20 atom c(2 × 2) × 2 unit cell calculations. While we did perform slab calculations including one layer of SrIrO3 followed by five layers of SrRuO3 and vacuum, these were not as easy to interpret in a layer-by-layer basis because of the large number of Wannier functions and their relatively delocalized nature.
Authors: J Matsuno; K Ihara; S Yamamura; H Wadati; K Ishii; V V Shankar; Hae-Young Kee; H Takagi Journal: Phys Rev Lett Date: 2015-06-18 Impact factor: 9.161
Authors: Alexandru B Georgescu; Oleg E Peil; Ankit S Disa; Antoine Georges; Andrew J Millis Journal: Proc Natl Acad Sci U S A Date: 2019-07-02 Impact factor: 11.205
Authors: Lin Hao; D Meyers; Clayton Frederick; Gilberto Fabbris; Junyi Yang; Nathan Traynor; Lukas Horak; Dominik Kriegner; Yongseong Choi; Jong-Woo Kim; Daniel Haskel; Phil J Ryan; M P M Dean; Jian Liu Journal: Phys Rev Lett Date: 2017-07-14 Impact factor: 9.161
Authors: R Yamada; J Fujioka; M Kawamura; S Sakai; M Hirayama; R Arita; T Okawa; D Hashizume; M Hoshino; Y Tokura Journal: Phys Rev Lett Date: 2019-11-22 Impact factor: 9.161
Authors: B M Wojek; M H Berntsen; S Boseggia; A T Boothroyd; D Prabhakaran; D F McMorrow; H M Rønnow; J Chang; O Tjernberg Journal: J Phys Condens Matter Date: 2012-10-17 Impact factor: 2.333
Authors: Julian L Schmehr; Thomas R Mion; Zach Porter; Michael Aling; Huibo Cao; Mary H Upton; Zahirul Islam; Rui-Hua He; Rajdeep Sensarma; Nandini Trivedi; Stephen D Wilson Journal: Phys Rev Lett Date: 2019-04-19 Impact factor: 9.161
Authors: Yue Cao; Qiang Wang; Justin A Waugh; Theodore J Reber; Haoxiang Li; Xiaoqing Zhou; Stephen Parham; S-R Park; Nicholas C Plumb; Eli Rotenberg; Aaron Bostwick; Jonathan D Denlinger; Tongfei Qi; Michael A Hermele; Gang Cao; Daniel S Dessau Journal: Nat Commun Date: 2016-04-22 Impact factor: 14.919