Open Access Article
Polina M.
Sheverdyaeva
*a,
Gustav
Bihlmayer
b,
Silvio
Modesti
cd,
Vitaliy
Feyer
e,
Matteo
Jugovac
ae,
Giovanni
Zamborlini
efg,
Christian
Tusche
eh,
Ying-Jiun
Chen
eh,
Xin Liang
Tan
e,
Kenta
Hagiwara
e,
Luca
Petaccia
i,
Sangeeta
Thakur
ij,
Asish K.
Kundu
akl,
Carlo
Carbone
a and
Paolo
Moras
a
aCNR-Istituto di Struttura della Materia (CNR-ISM), Strada Statale 14 km 163.5, 34149 Trieste, Italy. E-mail: polina.sheverdyaeva@ism.cnr.it
bPeter Grünberg Institut (PGI-1), Forschungszentrum Jülich and JARA, D-52425 Jülich, Germany
cDipartimento di Fisica, Università di Trieste, 34127 Trieste, Italy
dCNR-Istituto Officina dei Materiali (CNR-IOM), 34149 Trieste, Italy
ePeter Grünberg Institute (PGI-6), Forschungszentrum Jülich, 52428 Jülich, Germany
fTU Dortmund University, Otto-Hahn-Straße 4, 44227 Dortmund, Germany
gInstitute of Physics, University of Graz, Universitätsplatz 5, 8010 Graz, Austria
hFakultät für Physik, Universität Duisburg-Essen, 47057 Duisburg, Germany
iElettra Sincrotrone Trieste, Strada Statale 14 km 163.5, 34149 Trieste, Italy
jInstitut für Experimentalphysik, Freie Universität Berlin, Arnimallee 14, D-14195 Berlin, Germany
kNational Synchrotron Light Source II, Brookhaven National Laboratory, Upton, New York 11973, USA
lInternational Center for Theoretical Physics (ICTP), 34151 Trieste, Italy
First published on 2nd August 2024
Bismuth produces different types of ordered superstructures on the InAs(100) surface, depending on the growth procedure and coverage. The (2 × 1) phase forms at completion of one Bi monolayer and consists of a uniformly oriented array of parallel lines of Bi dimers. Scanning tunneling and core level spectroscopies demonstrate its metallic character, in contrast with the semiconducting properties expected on the basis of the electron counting principle. The weak electronic coupling among neighboring lines gives rise to quasi one-dimensional Bi-derived bands with open contours at the Fermi level. Spin- and angle-resolved photoelectron spectroscopy reveals a giant Rashba splitting of these bands, in good agreement with ab initio electronic structure calculations. The very high density of the dimer lines, the metallic and quasi one-dimensional band dispersion and the Rashba-like spin texture make the Bi/InAs(100)-(2 × 1) phase an intriguing system, where novel transport regimes can be studied.
The RB effect can significantly influence the electronic structure and spin texture of one-dimensional (1D) systems, such as quantum wires, and favor the emergence of a specific spintronic functionality.17 In a 1D system subject to the RB interaction the spin degeneracy of the spin-split bands at the time-reversal symmetry point can be removed by an external magnetic field opening a gap. If this gap opens at the Fermi level (EF), a pure and non-dissipative spin current can be established by applying a voltage. The exploitation of this mechanism is not trivial due to the scarcity of systems with genuine 1D RB-like spin texture. As an example, the 1D bands of Bi chains on Ag(110) display giant RB splitting, but the density of states of the system at EF is dominated by bulk Ag states.18 Semiconductor substrates, on the other hand, can support the formation of 1D structures with RB-split metallic states, as experimentally demonstrated for Au chains on Si(557),19 and theoretically predicted for Bi-adsorbed In atomic chains on Si(111).20
Bi is known to form different types of superstructures on III–V semiconducting surfaces.21–34 Some of these are based on Bi stripes displaying quasi 1D bands with giant RB splittings.28,30–32,34 Bi-terminated III–V semiconductors are expected to be non-metallic, according to the electron counting model, in order to decrease the surface energy.35 A notable exception to this model is the metallic behavior of Bi dimer lines grown on GaAs(100) (the so-called (2 × 1) phase), which emerges from scanning tunneling spectroscopy (STS) measurements and finds confirmation in density functional theory (DFT) calculations (without SOC).23,24 This full-monolayer phase made of symmetric dimers is stabilized by stress relief and pseudogap formation23 and was predicted to maintain long range compositional uniformity even at high temperature.36 An in-depth analysis aimed at establishing the dimensionality of the Bi-derived states and the magnitude of the RB effect in this and similar systems (Bi/GaAsxN1−x(100)-(2 × 1)22 and Bi/InAs(100)-(2 × 1)26) is still missing.
The present study reports on the electronic structure of the Bi/InAs(100)-(2 × 1) phase, which is examined by STS, photoelectron spectroscopy with spin analysis and DFT calculations. The (2 × 1) phase occurs at the completion of one Bi monolayer and consists of an array of uniformly oriented Bi dimer lines. The electronic coupling between neighboring lines turns out to be much weaker than along the dimer lines, thus giving rise to Bi-derived bands with highly anisotropic in-plane dispersion. The Fermi surface presents open quasi 1D contours with RB-like spin texture and giant splitting with αR values up to 4.6 eV Å, owing to the large SOC of Bi and the low structural symmetry. These findings suggest that the Bi/InAs(100)-(2 × 1) phase could support the generation of non-dissipative and spin-polarized currents and find application in spin-to-charge conversion processes.
DFT calculations were performed in the local density approximation43 using the full potential linearized augmented plane-wave method, as implemented in the Fleur code.44 The Bi/InAs(100)-(2 × 1) system was simulated with a slab of 21 atomic layers based on the symmetric dimer model proposed for the Bi/GaAs(100)-(2 × 1) system:24 one-atom-thick Bi dimer lines running along the [011] direction on the top surface; 19 atomic layers of InAs(100), with outermost In planes, as the substrate; hydrogen-terminated bottom surface. Fig. 1(d and e) display the top view and the three-dimensional representation of this structural model. The in-plane lattice constant (4.27 Å) and interlayer distance between In (or As) planes (3.02 Å) were set to those of bulk InAs. Structural relaxation on the Bi-terminated side of the slab led to 2.96 Å Bi–In interlayer distance, 3.07 Å Bi–Bi distance in the dimers and 0.5% expansion of the interlayer distance between the topmost In planes. Compared to the 1 × 1 Bi monolayer, about 0.31 eV/Bi are gained by Bi dimerization. Hydrogen atoms were placed on the bottom surface to saturate the dangling bonds. An attractive potential term of 2.67 eV was applied to the p-states of In and As to get a better description of the band gap and SOC was included in a self-consistent manner. The In terminations on both sides of the InAs substrate led to an artificial hole doping of the system and to a difference of 0.22 eV between the EF position in the calculations and experiments.
:
1 ratio between the long side (aligned with the [0
1] substrate axis) and the short side (aligned with the [011] substrate axis) is a signature of the dimerization of Bi atoms, which are not resolved individually in the present image,26 in analogy to similar systems.22,23 The dimers form an array of parallel lines running along the [011] axis. The continuity of the Bi dimer lines is occasionally interrupted by missing dimers (dark gray areas), as reported for the Bi/GaAs(100)-(2 × 1) phase.23 The structural features observed in the STM image are reproduced by the model of the Bi/InAs(100)-(2 × 1) phase. Bi atoms forming dimers present a reduced Bi–Bi distance along [0
1] with respect to the surface lattice constant of the substrate (Fig. 1(d)). These Bi dimers give rise to the parallel lines seen in Fig. 1(e).
Fig. 2 compares STS and photoelectron spectra of the clean InAs(100) surface and Bi/InAs(100)-(2 × 1) phase. The dI/dV signal of clean InAs(100) (top spectrum of Fig. 2(a)) is close to zero over an interval of 0.30 eV including EF that can be identified with the bulk bandgap of the substrate.45 The dI/dV signal of the Bi/InAs(100)-(2 × 1) phase (collected away from missing dimer regions, bottom spectrum of Fig. 2(a)) is more than a factor 10 higher than the noise level at EF, thus attesting the metallic character of the (2 × 1) phase. The peaks observed at −0.2, −0.08 and 0.08 eV (black arrows) will be discussed later in connection with the ARPES and DFT analyses.
Survey spectra of the photoelectron signal at hν = 70 eV for the two systems (Fig. 2(b)) are useful to compare the intensities of the core level lines. In both cases, the much higher signal of In with respect to As has combined structural and electronic origin, as the sputtering/annealing procedure results in an excess of In at the surface and In 4d states have 5.6 times larger photoionization cross section than As 3d states.46 The zoom of the photoelectron signal near EF for InAs(100) (top spectrum) and Bi/InAs(100)-(2 × 1) (bottom spectrum) in Fig. 2(c) highlights the semiconducting vs. metallic properties of the two systems.
Fig. 2(d–f) show the analysis of the As 3d, Bi 5d and In 4d levels of InAs(100) (top spectra) and Bi/InAs(100)-(2 × 1) (bottom spectra), respectively. The As 3d spectrum of clean InAs(100) can be fitted with two doublets related to sub-surface As atoms (red line, As 3d5/2 at 40.50 eV) and bulk-like As atoms (green line, As 3d5/2 at 40.69 eV). The sub-surface component is suppressed due to a restructuring of the (4 × 2)/c(8 × 2) termination, occurring in correspondence with the formation of the (2 × 1) phase. The Bi 5d spectrum can be fitted with one doublet (blue line, Bi 5d5/2 at 23.55 eV) very similar to that of metallic Bi,47 at variance with the two doublets used in the literature.26 The slight asymmetry of the peaks (tail on the high binding energy side) is interpreted as a signature of the metallic character of the Bi/InAs(100)-(2 × 1) phase,48 rather than ascribed to the presence of another phase.26 This interpretation is in agreement with the data of Fig. 2(a and c) and will be strengthened by the ARPES and DFT analyses reported in Fig. 3 and 4. Two doublets corresponding to surface (red line, In 4d5/2 at 17.48 eV) and bulk (green line, In 4d5/2 at 17.23 eV) components are sufficient to fit the In 4d spectrum of clean InAs(100). After the formation of the (2 × 1) phase, the surface doublet is quenched and another doublet appears on the low binding energy side of the bulk component. This new doublet has an asymmetric shape (blue line, In 4d5/2 at 16.93 eV) that is compatible with metallic In49 in contact with surface metallic Bi.
Fig. 3(a–c) show a constant energy cut of the photoelectron signal of the Bi/InAs(100)-(2 × 1) phase acquired at 0.15 eV below EF with the momentum microscope. This instrument allows to scan simultaneously an area of the (kx,ky) space including several surface Brillouin zones (SBZs) of the system. For clarity, the edges of the (2 × 1) SBZs (black dashed lines) are overlaid to the data. The central SBZ (rectangle with thick black edges,
and
) and its high symmetry points are reported in Fig. 3(b). The photoelectron signal is characterized by four bands S1–S4, which are elongated in the ky direction and cross the edges of neighboring SBZs. Fig. 3(c) displays the second derivative of the original data taken along the kx axis, to better visualize low intensity sections of S1–S4. All the features observed in Fig. 3(a–c) are absent in the clean substrate (Fig. S1 of the ESI†),50 reflect the periodicity of the (2 × 1) superstructure, and, therefore, can be identified as Bi-derived electronic states. The wavy S1/S2 and straight S3/S4 constant energy contours are schematically represented by red and blue continuous lines for positive kx values in Fig. 3(b), where the color is assigned on the basis of the spin analysis.
In the spin-resolved constant energy map of Fig. 3(d) the spin quantization axis (SQA), determined by the experimental geometry,40 is parallel to the ky axis. The intensity of red and blue colors for spin-up and spin-down states, respectively, is proportional to their spin projection along the SQA. The most evident feature of Fig. 3(d) is the full reversal of the spin texture with respect to the ky axis. Portions of the S1–S3 contours are highlighted by red/blue dashed lines for positive kx values to guide the eye. S1 and S2 display high and opposite spin polarizations, thus suggesting they form a RB pair. The spin analysis of S3 and S4 is hindered by their low intensity close to kx = 0 axis, while their high spin-polarization clearly emerges at larger |kx| values. The high spin-polarization of the S1–S4 states in Fig. 3(d) means that the spins of the Bi-derived bands are almost parallel or anti-parallel to the SQA.
The ARPES spectra of Fig. 3(e–i) show the energy–momentum dispersion of the S1–S4 bands along the segments marked by the green dashed lines in Fig. 3(b). Also in this case the data are presented in the second derivative form to enhance the sensitivity to weak features. The three ARPES maps taken along equivalent
directions (Fig. 3(e–g)) demonstrate the strong state-dependent modulation of the photoelectron signal in the (kx,ky) plane, due to matrix element effects. S1 and S2 clearly cross EF, thus making the Bi/InAs(100)-(2 × 1) phase metallic, in agreement with the analysis displayed in Fig. 2. This metallic character is a distinctive feature of the (2 × 1) phases that 1 ML Bi forms on GaAs(100)23,24 and GaAsxN1.x(100).22 All S1–S4 states can be observed simultaneously along the
–Ȳ–
direction of Fig. 3(h). The dispersion of S3 and S4 (Fig. 3(g and h) and Fig. S2 of the ESI†) closely reminds the RB-split states observed in the (2 × 1) phases of Bi on the (110) surfaces of III–V semiconductors.28,30–32 This suggests that also S3/S4 form a RB pair of bands, in analogy to S1/S2. Fig. 3(i) shows the ARPES spectra along the
direction. The flat feature represents the S3/S4 crossing point. The minimum at
(0.2 eV) and the flat dispersion about Ȳ (0.08–0.1 eV) can be directly correlated to the peaks observed at the same energies in the occupied part of the STS spectrum in Fig. 2(a). The STS peak at 0.08 eV above EF can be linked to the maxima of the S3 and S4 bands (see Fig. 4), which cross EF and, therefore, contribute to the metallicity of the system.
The data of Fig. 3 clearly demonstrate that all S1–S4 bands have an anisotropic in-plane dispersion, which can be associated to the structural properties of the (2 × 1) phase. Due to the formation of the Bi dimers, the separation between Bi atoms in neighboring dimer lines is larger than along the lines. Correspondingly, the electronic coupling of Bi-related states along the dimer lines is stronger that perpendicular to them, thus resulting in steeply dispersing and flat bands, respectively, and open elongated contours. Overall, the ARPES and spin polarization analyses allow to describe S1/S2 and S3/S4 as RB-split pairs of bands with quasi 1D and metallic character.
In order to interpret the experimental findings described above, the electronic structure of the Bi/InAs(100)-(2 × 1) system was computed by DFT using the model reported in Fig. 1(d and e), which is based on the symmetric dimer model proposed for the Bi/GaAs(100)-(2 × 1) system.24Fig. 4(a) shows the spin-resolved DFT band structure calculations along the high-symmetry directions of Bi/InAs(100)(2 × 1). The size of the symbols is proportional to the spin polarization of the states, with the SQA oriented along the ky axis. Red/blue colors correspond to up/down spin channels, in analogy to the experiment. The horizontal dashed line indicates the experimental position of EF, which lies 0.22 eV above EF in the calculations. This difference has been evaluated by aligning the experimental and calculated S3/S4 bands (Fig. 4(b)). It derives from the artificial hole doping due to the excess of In atoms in the model (10 In planes vs. 9 As planes). To ease the comparison between experiment and theory, the energy scale of the calculations will be referred to the experimental position of EF (right axis of Fig. 4(a)) from here onwards. The electronic structure in the proximity of EF is characterized by four spin-polarized and Bi-derived bands, which are labeled with S1–S4 in analogy to the experiment. The additional bands crossing EF along
and
and showing no spin-polarization originate from the H-terminated surface and are not relevant to the present discussion.
All S1–S4 bands display quasi 1D character, which emerges by comparing their steep dispersion along
and Ȳ–
(parallel to the Bi dimer lines) with their flat dispersion along
and
–
(perpendicular to the Bi dimer lines). The spin analysis reveals that S1/S2 and S3/S4 form two RB pairs. The strength of the RB effect can be evaluated in the proximity of
and Ȳ through αR = 2·ΔE/Δk, as shown in Fig. 4(b). The S1/S2 pair presents giant αR values of 4.6 eV Å along
–Ȳ–
and 3.5 eV Å along
. For the S3/S4 pair αR is 2.3 eV Å along
–Ȳ–
and 2.7 eV Å along
. In Fig. 4(b) the calculated bands are overlaid to the ARPES data of Fig. 3(h) to demonstrate the correspondence between theory and experiment. The slight offset of the experimental S1/S2 pair with respect to the calculations can indicate that the spin-splitting and, consequently, the αR value are larger than predicted. A similar effect is seen Fig. S2 of the ESI† for the S3/S4 pair. The overall good agreement allows to describe the Bi/InAs(100)-(2 × 1) phase as an array of Bi dimer lines with metallic character determined by two pairs of bands displaying quasi 1D dispersion and RB-like spin texture.
The RB states of the present system display a unique combination of metallic character, giant Rashba parameter and strong 1D anisotropy. The measured Rashba parameter of 4.6 eV A is larger than the giant values reported for the metallic states of BiTeI (3.85 eV A (ref. 51)), Bi/Ag(111) (3.05 eV A (ref. 6)), free-standing atomic Te chains (2.17 eV A (ref. 52)) and metallic nanowires on semiconducting substrates (0.8–3.5 eV A (ref. 20, 53 and 54)). It approaches the largest αR value so far observed in semiconducting Bi lines (5.5 eV A (ref. 28)). We emphasize that this giant value originates primarily from the large SOC of Bi and the 1D anisotropy of the dimer lines, and is rather insensitive to the choice of the potential and therefore of the substrate's gap (Fig. S3†). The unique properties of the Bi/InAs(100)-(2 × 1) phase can find application in spin-to-charge conversion via the inverse Edelstein effect.55 The Rashba-split surface states could carry a charge current in response to the nonzero spin density associated to spin injection from an adjacent ferromagnetic layer.12 Such an interface can be used to detect spin currents with an efficiency that is expressed by λ = αRτs/ħ, where τs is the spin-relaxation time.56 Assuming the average spin-relaxation time of 5 fs,28,56 one may expect an efficiency of about 3.5 nm, which is significantly larger than in 2D RB systems (0.3 nm (ref. 12)) and topological insulators (2.1 nm (ref. 57)). The decoupling of the S1/S2 surface bands from the bulk states of the semiconductor substrate near EF
15 and the reduced back-scattering associated to the quasi 1D dispersion58–60 could favor even larger τs values. From the fundamental point of view, the S1/S2 bands could be also used for studying exotic electronic phenomena, such as Majorana bound states61 and spin-dependent density waves.62 At variance with the S1/S2 bands, the maxima of the S3/S4 pair lie just above EF close to
. Most likely these maxima produce the peak located at 0.08 eV in the unoccupied part of the STS spectrum (Fig. 2(a)). In fact, S1 and S2 are expected to give rise to a featureless dI/dV signal, due to their steep dispersion and distance from
near EF. The vicinity of the S3/S4 crossing point to EF at
, which is a time-reversal symmetry point of the system, can realize the scenario depicted in ref. 17: an external magnetic field could open a gap between S3 and S4; if the gap includes EF, the application of a voltage would allow the flow of a non-dissipative, pure spin current through the Bi dimer lines. Notably, the exact location of the S3/S4 bands with respect to EF can be tuned by external doping or gating, thanks to the semiconducting nature of the substrate, to meet the requirements of ref. 17.
very close to EF, appears to be suitable to host pure spin-polarized and non-dissipative currents, upon the application of an external magnetic field. The experimental verification of these novel transport properties would open new perspectives for the exploitation of the RB effect in spintronic devices.
Footnote |
| † Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4nr01591f |
| This journal is © The Royal Society of Chemistry 2024 |