Abstract
Excess charge on polar surfaces of ionic compounds is commonly described by the two-dimensional electron gas (2DEG) model, a homogeneous distribution of charge, spatially-confined in a few atomic layers. Here, by combining scanning probe microscopy with density functional theory calculations, we show that excess charge on the polar TaO2 termination of KTaO3(001) forms more complex electronic states with different degrees of spatial and electronic localization: charge density waves (CDW) coexist with strongly-localized electron polarons and bipolarons. These surface electronic reconstructions, originating from the combined action of electron-lattice interaction and electronic correlation, are energetically more favorable than the 2DEG solution. They exhibit distinct spectroscopy signals and impact on the surface properties, as manifested by a local suppression of ferroelectric distortions.
Similar content being viewed by others
Introduction
Excess charge emerges spontaneously at the surfaces and interfaces of polar ionic compounds and dictates the physical and chemical material properties1. Usually, this intrinsic charge is described as a homogeneous two dimensional electron gas (2DEG), charge carriers spatially confined in the near-surface atomic layers2,3 that can trigger different phase transitions 4,5,6,7,8 and can be functionalized for a variety of applications9,10,11,12. On the other hand, signatures of ferromagnetic domains and superconductivity at the interfaces of these materials suggest the possibility for inhomogeneous distributions of the interfacial excess charge13,14,15,16,17,18. Moreover, non-intrinsic excess charge on neutral and polar surfaces, generated by chemical or lattice defects, is known to form localized polaron states19,20, excess holes or electrons coupled to local lattice distortions21. Localized polarons cause deep modification of surface properties, fundamentally different from those induced by a dispersed 2DEG, and play a key role in all applications involving charge transport and catalysis22,23,24,25. One question that has remained unanswered to date is whether intrinsic excess charge on polar surfaces can spontaneously form localized polaron states or charge inhomogeneities on the defect-free termination, thus providing an alternative scenario beyond the commonly accepted, delocalized 2DEG picture.
Here, by combining scanning probe microscopy experiments with density-functional theory (DFT) simulations, we examine the (001) surface of the cubic quantum-paraelectric perovskite KTaO326. This is a prototypical polar surface expected to host a homogeneous 2DEG on the ferroelectrically-polarized TaO2 termination27. Our data show that the intrinsic excess charge on KTaO3 is instead accommodated through the formation of charge-density waves (CDW, periodic charge modulations occurring in a variety of materials)28,29,30,31 and highly-localized small electron polarons and bipolarons (singly and doubly reduced Ta atoms, respectively). To weaken their mutual repulsive interaction, polarons and bipolarons form ordered patterns, and are found to coexist with the CDW. The coupling between charge localization and the crystal lattice can induce significant structural distortions that disrupt the surface ferroelectric polarization.
The formation of ordered patterns of (bi)polarons and a CDW on the ferroelectrically active KTaO3(001) is expected to impact the functionalities of the material as compared to the simple 2DEG picture. From a more fundamental point of view, our study demonstrates much richer physics of polar surfaces than considered so far, revealing charge modulation as an intrinsic property of the bare surface. In the following, we start by showing the experimental evidence for the coexistence of different types of electronic states on the surface. We then analyze these states on the basis of first-principles DFT results.
Results and discussion
Figure 1 shows the results of our surface-sensitive experiments on cleaved KTaO3(001) samples, obtained from slightly n-type-doped crystals with substitutions of Ca atoms on K sites to ensure bulk electrical conductivity for the measurements (an extended data set on variously doped crystals is shown in the Supplementary Figs. 1–7). Figure 1a shows a constant-height, non-contact atomic force microscopy (ncAFM) image of the bulk-terminated surface with atomically resolved KO terraces alternating with regions of TaO227. (The TaO2 regions lie lower by half a unit cell and are unresolved in the image.) Figure 1b–f show local conductance (LDOS) maps and scanning tunneling spectroscopy (STS) curves measured in the same region.
The occupied LDOS shows a complex structure confined on the TaO2 regions, while the vicinity of the steps (i.e., borders of the TaO2 terraces) as well as the KO terraces are free of in-gap states27. The electronic structure observed near the Fermi level (Fig. 1b) shows a wave pattern with a wavelength (peak-peak) of 1.55 ± 0.15 nm, typically rotated by ≈45 ± 10∘ with respect to the 〈100〉 directions (see Supplementary Figs. 1, 2). This has been previously interpreted as a simple 2DEG27,32,33, but detailed analysis of STS data brings a solid evidence for the presence of a charge density wave: There is a small band gap at the Fermi level (Fig. 1f), the wavelength does not disperse in energy (Supplementary Fig. 1), the empty states show a phase reversal with respect to the filled states (Supplementary Fig. 4)31, and the waves disappear above ≈80 K (Supplementary Fig. 6). At energies below ≈−0.1 V, where the LDOS is not dominated by the CDW, localized spots appear in the LDOS maps (Fig. 1c–e), indicating the presence of small polarons. These polaronic states exhibit a partially ordered pattern and are located in the minima of the occupied CDW states, as expected for repulsion of negative charges.
To gain insight into the experimental observations, we analyzed different electronic surface states (2DEG, CDW and (bi)polarons) individually, by performing DFT calculations modeling the TaO2 termination of KTaO3 with \(2\sqrt{2}\times 2\sqrt{2}\)-large slabs, and including an on-site interaction parameter of U = 4 eV to account for electronic correlation (see “Methods” section). Figure 2 highlights the essential structural and electronic features of the surface. By adopting a simple oxidation-state model, KTaO3 can be considered as built by negatively (KO)− and positively (TaO2)+ charged layers, stacking along [001] with large interlayer charge transfer 27. Broken bonds at the TaO2 surface boundary interrupt the charge transfer leading to an accumulation of 1/2 electron per 1 × 1 surface unit cell (as sketched in Fig. 2a)34. These uncompensated electrons constitute the intrinsic excess charge of the polar surface.
The conventional homogeneous 2DEG solution is shown in Fig. 2b, c. The 2DEG spreads over the surface and shallow sub-surface TaO2 layers, and it is associated with metallic states at the bottom of the conduction band (see the density of states (DOS) in Fig. 2b); conversely, deeper TaO2 layers are not affected by the metallic 2DEG and retrieve the bulk-like semiconducting DOS. The top view and the simulated scanning tunneling microscopy (STM) images in Fig. 2c further highlight the spatial homogeneity of the 2DEG on the Ta sites of the 1 × 1 uniform surface lattice. We note that the excess electrons here carry parallel spin moments, forming a spin-polarized 2DEG with the conduction band of the minority spin channel completely empty (a spin-degenerate 2DEG would be 50 meV per excess electron less stable, see Supplementary Fig. 8). Importantly, the surface undergoes a ferroelectriclike (FE) distortion: While atoms on bulk KO and TaO2 layers are perfectly aligned on their respective planes (as expected for a quantum paraelectric crystal as KTaO3)35, O atoms on the surface appear displaced outwards along [001] by δFE = 0.21 Å from the Ta plane (see also Table 1, and Supplementary Fig. 9 for further details on the surface structure). This is a typical feature of bulk-terminated polar surfaces, as ferroelectric displacements create local dipole moments opposed to the internal electric field, in order to stabilize the surface against the diverging internal electrostatic potential and the polar catastrophe27,36,37,38.
Figure 3 collects the alternative solutions (CDW and (bi)polarons), all significantly more stable than the homogeneous 2DEG (see Table 1). A spin-polarized CDW can form on the surface (ΔE = −50 meV per excess electron as compared to the homogeneous 2DEG), showing a weak modulation of the charge with \(\sqrt{2}\times \sqrt{2}\) periodicity, revealed by magnetic moments of 0.5 and 0.2 μB on Ta ions in a checkerboard-ordered pattern (Fig. 3a). We note that the formation of this CDW does not require any additional structural distortion with respect to the homogeneous 2DEG phase, preserving the FE distortions on the surface (see Table 1). The electronic fingerprint of the CDW phase consists of a shallow peak right below the Fermi level, and an almost vanishing LDOS signal at EF (see DOS in Fig. 3a), in agreement with the experimental observations (Fig. 1f).
Figure 3 b,c show the most representative localized (bi)polaron solutions. The most favorable phase is formed by an ordered arrangement of bipolarons and single-electron polarons (Fig. 3b, ΔE = −250 meV). Alternative (bi)polaronic solutions are shown in Fig. 3c (bipolarons on the surface aligned along the [110] direction, ΔE = −180 meV), Supplementary Fig. 10 (sub-surface single-electron polarons, ΔE = −230 meV) and Supplementary Fig. 11 ([100]-ordered bipolarons, ΔE = −100 meV). The bipolaronic Ta3+ surface sites posses a magnetic moment of 1.6 μB, in contrast to 0.8 μB in surface and subsurface single-electron polaron Ta4+ sites. The (bi)polaronic states are insulating and associated to in-gap peaks (see DOS in Fig. 3b, c) similar to the peaks at deep energies observed in the experiment (Fig. 1f). The two electrons forming the bipolaron in Fig. 3c give rise to degenerate in-gap peaks, which are split in the mixed phase in Fig. 3b due to the repulsive interaction with the surrounding single-electron polarons.
Due to the particularly strong charge confinement, bipolarons affect the lattice structure considerably, causing breathing-out displacements of oxygen atoms surrounding the Ta3+ site (Ta–O bond length increased by ΔbTa–O = +0.09 Å), and a local quenching of the ferroelectric distortions (δFE = 0.08 Å). Single-electron polarons instead induce smaller lattice distortions (ΔbTa–O = +0.03 Å) and preserve the surface ferroelectricity (δFE = 0.23 and 0.15 Å for subsurface and surface polarons, respectively).
In order to acquire a detailed understanding of the various states, we show in Fig. 4 the corresponding band structures. In the 2DEG phase the t2g bands (dxy, dxz, and dyz) cross the Fermi level in the majority spin channel, whereas the corresponding minority-spin states are unoccupied (Fig. 4a). The observed band splitting is mainly caused by the internal electric field of the polar surface32,39. In the CDW phase, the charge modulation is associated to low-dispersion dxz and dyz bands lying below the Fermi level, while the dxy bands are completely empty (Fig. 4b). The strong charge localization and the coupling with lattice distortions associated to the formation of polarons (Fig. 4c) and bipolarons (Fig. 4d) tends to form essentially dispersionless in-gap bands: deviations from completely flat bands reveal the repulsive interactions perturbing the (bi)polaronic states40.
Summing up at this point, we found that the uncompensated surface charge can be modeled with three different electronic phases: 2DEG, CDW, and (bi)polarons. The stability of these phases is controlled by the degree of charge ordering (facilitated by electron localization associated with sizable electron correlation) and local lattice distortions (Ta–O bond length and ferroelectriclike distortions), mainly originating from electron-phonon interactions and phonon instabilities. By removing the on-site U, charge localization is suppressed and the homogeneous 2DEG becomes the dominant phase (see Supplementary Fig. 12)32,33, in disagreement with our experimental observations. To decode the effect of structural distortions, we performed a series of calculations from the uniform 2DEG limit to the distorted mixed polaronic solution (i.e., the global energy minimum configuration shown in Fig. 3b) by interpolating the ion coordinates between these two end states (Fig. 5): In the FE-distorted surface (in the absence of polaronic-like oxygen breathing distortions), the CDW phase is more stable than the homogeneous 2DEG. By switching on polaron breathing displacements and modeling the evolution towards the fully bipolaron state, the system passes through an intermediate phase characterized by the coexistence of CDW and single-electron polarons (Fig. 5a). The transition from the bipolaron to the mixed polarons+bipolaron state is achieved by splitting one bipolaron into two single-electron polarons, as shown in the insets of Fig. 5 and corresponding energy profile. The relative stability of one phase over the other is controlled by the degree of lattice distortion that plays the role of the order-parameter in analogy to similar quantum-critical ferroelectric and multiferroics transitions26,41.
Our computational data show that polaronic states can coexist with the CDW, in agreement with the experiment. We note that a precise characterization of the relative stability of the electronic surface states is hampered in the DFT calculations by the dependence of the free energy on the computational setup (see Supplementary Fig. 12)42; moreover, the infinitely large TaO2 terminations in the calculations neglect the effects of the surface electrostatic field arising from the neighboring KO terraces, known to alter the surface electronic states27; also, the limited size of the DFT unit cell constrains the electronic phases, the CDW in particular. Nevertheless, the small differences in the formation energy of the different phases (−50 to −250 meV), including the CDW+polaron phase, support the possibility of a coexistence of these states as observed on the experimental samples. In fact, the DOS in Fig. 3 and band structures in Fig. 4 allow for interpretation of the experimental STS results from Fig. 1. The STS shoulder at −60 mV is at the bottom of the conduction band, separated from the empty states by a small gap induced by the CDW. The simulated CDW exhibits a wavelength of about 5.7 Å (see Fig. 3a), increased to 11.4 Å in the mixed CDW+polaron phase (see Supplementary Fig. 13), slightly shorter than the experimental one (15.5 ± 1.5 Å; a precise modeling of the experimentally observed CDW would require a prohibitively large supercell size). The peaks measured at −0.6 and ≈−1.3 V can be attributed to polaron or bipolaron states. While LDOS measurements do not allow us for unambiguous assignment of these deep states with specific polaronic types, support for predominance of Ta3+ bipolarons in KTaO3 over Ta4+ single-electron polarons can be found for instance in photoluminiscence experiments43,44. Further STS maxima at more negative bias depend on the tip; possibly these (as well as the negative tunneling conductance sometimes observed) are related to changes of the charge configuration by the electric field of the tip.
In summary, we have shown the emergence of coexisting charge density waves, polarons, and bipolarons originating from the intrinsic uncompensated charge of the polar TaO2-terminated KTaO3(001) surface. We have obtained a direct, real-space view of these states, in contrast to similar phenomena at interfaces of heterostructured materials17,45,46,47, where the charge distribution is not directly accessible. The competition between charge-ordering (CDW and (bi)polaronic states) and homogeneous distribution (2DEG) arises from the contrasting action of spatially extended 5d orbitals of Ta (favoring charge delocalization) and strong electronic correlation (favoring charge localization). The entanglement of the electronic degree of freedom with a phonon-active lattice enabled by electron-phonon coupling is crucial for the stabilization of charge trapping in the form of (bi)polarons, associated with breathing-out displacements of oxygen atoms surrounding Ta. The charge ordering originating from the large amount of uncompensated charge available on the polar surface is an interesting feature, which might explain the long lifetimes of polarized microdomains induced by photocarriers43,48. The possibility to control surface ferroelectric polarization by charge trapping through the formation of bipolarons represents a novel physical effect that could be important for applications such as piezocatalysis or pyrocatalysis49, where the surfaces turns polar once the bulk material switches between the ferroelectric and paraelectric states.
Methods
Computational methods
DFT calculations were performed by using the Vienna ab initio simulation package (VASP)50,51,52. We adopted the strongly constrained and appropriately normed meta-generalized gradient approximation (SCAN)53, with the inclusion of an on-site effective U54 of 4.0 eV on the d orbitals of Ta atoms (see Supplementary Fig. 12). The TaO2 termination of KTaO3(001) was modeled by using a slab with six TaO2 layers alternating with five KO layers, in a symmetric setup (mirror symmetry on the central (001) KO layer), including a vacuum region of more than 30 Å. To study the penetration of the 2DEG into deep layers, we used slabs with up to ten and nine TaO2 and KO layers, respectively. In order to enable spatial symmetry breaking we used \(2\sqrt{2}\times 2\sqrt{2}\) large slabs55,56. All the atomic coordinates (except those on the central KO layer) were relaxed using standard convergence criteria (residual forces smaller than 0.01 eV/Å), with a plane-wave energy cutoff of 500 eV, and a 3 × 3 × 1 grid for the integration in the reciprocal space, while electronic self consistence and densities of states were calculated by using a finer 7 × 7 × 1 grid. The various solutions for the intrinsic, uncompensated charge were obtained by adopting different initial conditions for the electronic density and wavefunction in unconstrained self-consistent calculations (see also Supplementary Information). We used VESTA57 to show the spatial extension of the excess charge and the Tersoff–Hamann approximation58 for the STM simulations, including all in-gap states and states at the bottom of the conduction band up to the Fermi energy level. In order to analyze the band structure of our large unit cells, we adopted the unfolding technique as implemented in ref. 40, by calculating the coherence (or Bloch character) of the supercell eigenstates unfolded in the Brillouin zone of a 1 × 1 surface slab.
The stability ΔE of charge density waves and (bi)polarons was calculated by comparing the total free energy E with the total free energy E2DEG obtained for the reference spin-polarized 2DEG: ΔE = (E − E2DEG)/8. The factor 1/8 scales the energy to the \(\sqrt{2}\times \sqrt{2}\) surface cell that contains one excess electron (note that our \(2\sqrt{2}\times 2\sqrt{2}\) slab contains 8 intrinsic excess electrons, i.e., 4 excess electrons per side of the symmetric slab); in the specific case of single-electron polaronic states, ΔE is equivalent to the polaron formation energy EPOL21.
Experimental methods
Combined STM/AFM measurements were performed at a temperature T of 4.8 K in a UHV chamber with a base pressure of 10−9 Pa, equipped with a commercial ScientaOmicron q-Plus LT head. Tuning-fork-based AFM sensors with a separate wire for the tunneling current (k = 1900 N/m, f0 = 30,500 Hz, Q ≃ 30,000)59, as well as a custom-design cryogenic preamplifier60, were used for the AFM measurements. Electrochemically etched W tips were glued to the tuning fork and cleaned in situ by field emission and self-sputtering and treated on a Cu(001) surface to ensure a metallic character of the tip. STS spectra were measured at open-loop conditions using a lock-in amplifier, with a modulation frequency of 123 Hz and an amplitude of 10 mV. The KTaO3 samples were prepared by solidification from a nonstoichiometric melt. Samples with 0.2% Ca doping were used, which ensures enough bulk electrical conductivity for performing the STS measurements. The cleaving was performed in situ by a tungsten carbide blade in a temperature range between 250 and 300 K, as described in ref. 27.
Data availability
Data generated or analyzed during this study are included in this published article (and its supplementary information files).
References
Santander-Syro, A. F. et al. Two-dimensional electron gas with universal subbands at the surface of SrTiO3. Nature 469, 189–193 (2011).
Wang, Z. et al. Anisotropic two-dimensional electron gas at SrTiO3(110). Proc. Natl Acad. Sci. USA 111, 3933–3937 (2014).
Wander, A. et al. Stability of polar oxide surfaces. Phys. Rev. Lett. 86, 3811 (2001).
Liu, C. et al. Two-dimensional superconductivity and anisotropic transport at KTaO3 (111) interfaces. Science 371, 716–721 (2021).
Chen, Z. et al. Two-dimensional superconductivity at the LaAlO3/KTaO3(110) heterointerface. Phys. Rev. Lett. 126, 026802 (2021).
Gattinoni, C. et al. Interface and surface stabilization of the polarization in ferroelectric thin films. Proc. Natl Acad. Sci. USA 117, 28589–28595 (2020).
Brinkman, A. et al. Magnetic effects at the interface between non-magnetic oxides. Nat. Mater. 6, 493–496 (2007).
Gariglio, S., Caviglia, A. D., Triscone, J.-M. & Gabay, M. A spin-orbit playground: Surfaces and interfaces of transition metal oxides. Rep. Prog. Phys. 82, 012501 (2019).
Sopiha, K. V., Malyi, O. I., Persson, C. & Wu, P. Suppression of surfaces states at cubic perovskite (001) surfaces by CO2 adsorption. Phys. Chem. Chem. Phys. 20, 18828 (2018).
Vaz, D. C. et al. Mapping spin-charge conversion to the band structure in a topological oxide two-dimensional electron gas. Nat. Mater. 18, 1187–1193 (2019).
Chen, Y. & Green, R. J. Progress and perspectives of atomically engineered perovskite oxide interfaces for electronics and electrocatalysts. Adv. Mater. Interfaces 6, 1900547 (2019).
Chen, J. & Park, N. G. Materials and methods for interface engineering toward stable and efficient perovskite solar cells. ACS Energy Lett. 5, 2742–2786 (2020).
Li, L., Richter, C., Mannhart, J. & Ashoori, R. C. Coexistence of magnetic order and two-dimensional superconductivity at LaAlO3/SrTiO3 interfaces. Nat. Phys. 7, 762–766 (2011).
Ariando, A. et al. Electronic phase separation at the LaAlO3/SrTiO3 interface. Nat. Commun. 2, 188 (2011).
Ristic, Z. et al. Nanoscale modulation of the density of states at the conducting interface between LaAlO3 and SrTiO3 band insulators. EPL (Europhys. Lett.) 93, 17004 (2011).
Bucheli, D., Caprara, S. & Grilli, M. Pseudo-gap as a signature of inhomogeneous superconductivity in oxide interfaces. Supercond. Sci. Technol. 28, 045004 (2015).
Bovenzi, N. et al. Density inhomogeneities and Rashba spin–orbit coupling interplay in oxide interfaces. J. Phys. Chem. Solids 128, 118–129 (2019).
Zhang, B. et al. Re-entrance to a ferromagnetic insulator with oxygen-vacancy ordering in the La0.7Sr0.3MnO3/SrTiO3 superlattice. J. Mater. Chem. A 9, 26717–26726 (2021).
Papageorgiou, A. C. et al. Electron traps and their effect on the surface chemistry of TiO2(110). Proc. Natl Acad. Sci. USA 107, 2391–2396 (2010).
Gerosa, M., Gygi, F., Govoni, M. & Galli, G. The role of defects and excess surface charges at finite temperature for optimizing oxide photoabsorbers. Nat. Mater. 17, 1122–1127 (2018).
Franchini, C., Reticcioli, M., Setvin, M. & Diebold, U. Polarons in materials. Nat. Rev. Mater. 6, 560–586 (2021).
Reticcioli, M. et al. Polaron-driven surface reconstructions. Phys. Rev. X 7, 031053 (2017).
Albanese, E., Di Valentin, C. & Pacchioni, G. H2O adsorption on WO3 and WO3−x (001) surfaces. ACS Appl. Mater. Interfaces 9, 23212–23221 (2017).
Rousseau, R., Glezakou, V. A. & Selloni, A. Theoretical insights into the surface physics and chemistry of redox-active oxides. Nat. Rev. Mater. 5, 460–475 (2020).
Zhang, D., Han, Z. K., Murgida, G. E., Ganduglia-Pirovano, M. V. & Gao, Y. Oxygen-vacancy dynamics and entanglement with polaron hopping at the reduced CeO2(111) surface. Phys. Rev. Lett. 122, 096101 (2019).
Rowley, S. E. et al. Ferroelectric quantum criticality. Nat. Phys. 10, 367–372 (2014).
Setvin, M. et al. Polarity compensation mechanisms on the perovskite surface KTaO3(001). Science 359, 572–575 (2018).
Chen, C.-W., Choe, J. & Morosan, E. Charge density waves in strongly correlated electron systems. Rep. Prog. Phys. 79, 084505 (2016).
Balandin, A. A., Zaitsev-Zotov, S. V. & Grüner, G. Charge-density-wave quantum materials and devices – New developments and future prospects. Appl. Phys. Lett. 119, 170401 (2021).
Zhu, X., Cao, Y., Zhang, J., Plummer, E. W. & Guo, J. Classification of charge density waves based on their nature. Proc. Natl Acad. Sci. USA 112, 2367–2371 (2015).
Spera, M. et al. Insight into the charge density wave gap from contrast inversion in topographic STM images. Phys. Rev. Lett. 125, 267603 (2020).
King, P. D. et al. Subband structure of a two-dimensional electron gas formed at the polar surface of the strong spin-orbit perovskite KTaO3. Phys. Rev. Lett. 108, 117602 (2012).
Santander-Syro, A. F. et al. Orbital symmetry reconstruction and strong mass renormalization in the two-dimensional electron gas at the surface of KTaO3. Phys. Rev. B 86, 121107 (2012).
Nakagawa, N., Hwang, H. Y. & Muller, D. A. Why some interfaces cannot be sharp. Nat. Mater. 5, 204��209 (2006).
Tyunina, M. et al. Evidence for strain-induced ferroelectric order in epitaxial thin-film KTaO3. Phys. Rev. Lett. 104, 227601 (2010).
Goniakowski, J., Finocchi, F. & Noguera, C. Polarity of oxide surfaces and nanostructures. Rep. Prog. Phys. 71, 016501 (2008).
Springborg, M., Zhou, M., Molayem, M. & Kirtman, B. Surfaces, shapes, and bulk properties of crystals. J. Phys. Chem. C 122, 11926–11932 (2018).
Hwang, J. et al. Tuning perovskite oxides by strain: Electronic structure, properties, and functions in (electro)catalysis and ferroelectricity. Mater. Today 31, 100–118 (2019).
Shanavas, K. V. & Satpathy, S. Electric field tuning of the Rashba effect in the polar perovskite structures. Phys. Rev. Lett. 112, 086802 (2014).
Dirnberger, D., Kresse, G., Franchini, C. & Reticcioli, M. Electronic state unfolding for plane waves: Energy bands, Fermi surfaces, and spectral functions. J. Phys. Chem. C 125, 12921–12928 (2021).
Narayan, A., Cano, A., Balatsky, A. V. & Spaldin, N. A. Multiferroic quantum criticality. Nat. Mater. 18, 223–228 (2019).
Reticcioli, M., Diebold, U. & Franchini, C. Modeling polarons in density functional theory: Lessons learned from TiO2. J. Phys.: Condensed Matter 34, 204006 (2022).
Grenier, P., Bernier, G., Jandl, S., Salce, B. & Boatner, L. A. Fluorescence and ferroelectric microregions in KTaO3. J. Phys.: Condensed Matter 1, 2515 (1989).
Vikhnin, V., Eden, S., Aulich, M. & Kapphan, S. The origin of the red luminescence in incipient ferroelectric KTaO3. Solid State Commun. 113, 455–460 (2000).
Cancellieri, C. et al. Polaronic metal state at the LaAlO3/SrTiO3 interface. Nat. Commun. 7, 10386 (2016).
Strocov, V. N., Cancellieri, C. & Mishchenko, A. S. Electrons and Polarons at Oxide Interfaces Explored by Soft-X-Ray ARPES. In: (eds Cancellieri, C. & Strocov, V.) Spectroscopy of Complex Oxide Interfaces. Springer Series in Materials Science, vol 266. (Springer, Cham, 2018).
Kong, W. et al. Formation of two-dimensional small polarons at the conducting LaAlO3/SrTiO3 interface. Phys. Rev. B 100, 085413 (2019).
Kapphan, S. E., Gubaev, A. I. & Vikhnin, V. S. UV light-induced absorption in pure KTaO3 at low temperature. Phys. Status Solidi C 2, 128–131 (2005).
Wang, M., Wang, B., Huang, F. & Lin, Z. Enabling PIEZOpotential in PIEZOelectric semiconductors for enhanced catalytic activities. Angew. Chem. 131, 7606–7616 (2019).
Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 54, 11169–11186 (1996).
Kresse, G. & Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. Mater. Sci. 6, 15–50 (1996).
Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B 59, 1758 (1999).
Sun, J. et al. Accurate first-principles structures and energies of diversely bonded systems from an efficient density functional. Nat. Chem. 8, 831–836 (2016).
Dudarev, S. & Botton, G. Electron-energy-loss spectra and the structural stability of nickel oxide: An LSDA+U study. Phys. Rev. B 57, 1505 (1998).
Varignon, J., Bibes, M. & Zunger, A. Origin of band gaps in 3d perovskite oxides. Nat. Commun. 10, 1658 (2019).
Malyi, O. I. & Zunger, A. False metals, real insulators, and degenerate gapped metals. Appl. Phys. Rev. 7, 041310 (2020).
Momma, K. & Izumi, F. VESTA 3 for three-dimensional visualization of crystal, volumetric, and morphology data. J. Appl. Crystallogr. 44, 1272–1276 (2011).
Tersoff, J. & Hamann, D. R. Theory and application for the scanning tunneling microscope. Phys. Rev. Lett. 50, 1998 (1983).
Giessibl, F. J. The qPlus sensor, a powerful core for the atomic force microscope. Rev. Sci. Instrum. 90, 011101 (2019).
Huber, F. & Giessibl, F. J. Low noise current preamplifier for qPlus sensor deflection signal detection in atomic force microscopy at room and low temperatures. Rev. Sci. Instrum. 88, 073702 (2017).
Acknowledgements
This work was supported by the Austrian Science Fund (FWF) projects P32148-N36 (SUPER, funding received by C.F. and M. Se) and SFB TACO (F81, by U.D. and C.F.), and by the joint Austrian (BMBWF CZ15/2021) Czech (MSMT 8J21AT004) project (by M.R. and M. Se). Martin Setvin acknowledges the support of the Czech Science Foundation (GACR 20-21727X and GAUK Primus/20/SCI/009). Calculations were performed at the Vienna Scientific Cluster.
Author information
Authors and Affiliations
Contributions
M. Se and C.F. designed the project. M.R. and Z.W. contributed equally to this work. M.R. performed the calculations, Z.W. conducted the experiments. All authors (M.R, Z.W., M. Sc, D.W., L.B., U.D., M. Se, and C.F.) contributed to a substantial analysis and discussion of the results, and to the manuscript writing.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Communications thanks Chi Ming Yim, and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Reticcioli, M., Wang, Z., Schmid, M. et al. Competing electronic states emerging on polar surfaces. Nat Commun 13, 4311 (2022). https://doi.org/10.1038/s41467-022-31953-6
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-022-31953-6
This article is cited by
-
Direct observation of strong surface reconstruction in partially reduced nickelate films
Nature Communications (2024)
-
Spin-orbital Jahn-Teller bipolarons
Nature Communications (2024)