A strategy to create spin-split metallic bands on silicon using a dense alloy layer.

To exploit Rashba effect in a 2D electron gas on silicon surface for spin transport, it is necessary to have surface reconstruction with spin-split metallic surface-state bands. However, metals with strong spin-orbit coupling (e.g., Bi, Tl, Sb, Pt) induce reconstructions on silicon with almost exclusively spin-split insulating bands. We propose a strategy to create spin-split metallic bands using a dense 2D alloy layer containing a metal with strong spin-orbit coupling and another metal to modify the surface reconstruction. Here we report two examples, i.e., alloying reconstruction with Na and Tl/Si(111)1 × 1 reconstruction with Pb. The strategy provides a new paradigm for creating metallic surface state bands with various spin textures on silicon and therefore enhances the possibility to integrate fascinating and promising capabilities of spintronics with current semiconductor technology.

The Rashba spin splitting12 in the two-dimensional electron gas systems on semiconductors is considered to be the key concept for many promising spintronics applications. To combine Rashba-effect based spintronics with a silicon technology metallic surface states with a strong spin-orbit coupling on a silicon surface34 are demanded. Unfortunately, the most of the metal/silicon systems with spin-split surface state bands (e.g., 567, Tl/Si(111)1 × 1 8910, 11 and Pt/Si(110)“6” × 5 12) are semiconducting. The only known exceptions are the Tl/Si(111)1 × 1 modified by adsorption of additional Tl13 and modified by adsorption of In, Tl, Cs or Na4. It is worth noting that in both cases adsorbates do not alter the basic atomic arrangement of the pristine surface and their effect resides just in eliminating or generating surface defects (as for Au/Si(111) and Tl/Si(111), respectively) and doping electrons to the available surface-state bands. In the present study, we propose a novel more universal strategy for tailoring the spin-split metallic surface states. The main concept resides in formation of the 2D alloy layers containing a metal with a strong spin-orbit coupling and another suitable metal to obtain dense reconstructions with spin-split metallic bands. Validity of the approach is demonstrated by the formation of the Bi-Na and Tl-Pb 2D alloys on Si(111) surface possessing the required properties.

Results

Figure 1 (left panel) summarizes the main structural and electronic properties of the pristine surface. It contains 1.0 ML (monolayer, 7.8 × 1014 cm−2) of Bi arranged into the trimers centered at the T4 sites on the bulk-truncated Si(111) surface14. The main features of its electronic structure are spin-split non-metallic bands, S1 (S1′) and S2 (S2′)6. Room-temperature (RT) deposition of 0.33 ML Na onto this surface results in the formation of the ordered Bi-Na alloy layer whose structural and electronic properties are shown in Figure 1 (right panel). Comparison of its atomic structure with that of the original surface demonstrates that the Bi trimers are preserved but they increase in size by ~8% (Bi-Bi bond length changes from 3.13 to 3.34 Å) and twisted by ±6.8°. The resultant arrangement of Bi atoms at the surface can be thought then as a honeycomb network of the chained Bi trimers with Na atoms occupying the T4 sites in the center of each honeycomb unit.

Figure 1

Structural and electronic properties of the (left panel) and (right panel).

(a) and (e) 30 × 30 Å2 STM images. Upper and lower halves in (e) present experimental and simulated empty-state (+0.3 V) STM images, respectively. (b) and (f) Structural models where Bi atoms are shown by violet circles, Na atoms by green circles, and Si atoms by gray circles. Positions of the outlined unit cells in (e) and (f) correspond to each other. The inset in (b) shows a sketch of reciprocal space geometry with boundaries of the and 1 × 1 SBZs given by solid black and dashed blue lines, respectively. (c) and (g) Experimental ARPES spectra and (d) and (h) calculated electronic band structures. In the calculated band structures, the size of green circles corresponds to the strength of the surface character summed over all orbitals at a particular k|| value.

Changes in atomic structure of the surface cause modification of its electronic band structure as revealed by angle-resolved photoelectron spectroscopy (ARPES) measurements and density-functional theory (DFT) calculations (Figure 1g and h). We would like to remark that a perfect agreement between the calculated band structure and that determined in the experiment as well as coincidence of the experimental and simulated STM images can serve an indication of the proper (Bi, Na)/Si(111) structural model. The most essential feature in the Bi-Na alloy band structure is the appeared metallic spin-split surface-state band denoted by Σ1 (Σ1′).

Consider now the Pb-induced modification of the Tl/Si(111)1 × 1 surface. Figure 2 (left panel) presents the main structural and electronic properties of the pristine Tl/Si(111)1 × 1 reconstruction. It contains 1.0 ML of Tl atoms occupying every T4 site on the bulk-truncated Si(111) surface151617. The band structure of the Tl/Si(111)1 × 1 surface includes the insulating S1 surface-state band whose occurrence was proved before both experimentally815 and theoretically915. Spin-splitting of this band has been concluded from the spin-resolved ARPES results8 and from previous9 and current theoretical calculations (Figure 2d) and is clearly resolved in the spin-unpolarized ARPES spectrum in Figure 2c. Additional feature can be noted, namely, the shallow metallic band S2 around the point. It has recently been recognized that electron filling of this band is associated with extra Tl atoms on the Tl/Si(111)1 × 1 surface13 in the form of specific surface defects18.

Figure 2

Structural and electronic properties of the Tl/Si(111)1 × 1 (left panel) and (right panel).

(a) and (e) 30 × 30 Å2 empty-state (+0.4 V) STM images. Upper and lower halves in (e) present experimental and simulated empty-state (+0.4 V) STM images, respectively. (b) and (f) Structural models where Tl atoms are shown by violet circles, Pb atoms by green circles, and Si atoms by gray circles. Positions of the outlined unit cells in (e) and (f) correspond to each other. (c) and (g) Experimental ARPES spectra and (d) and (h) calculated electronic band structures.

Adding 0.33 ML of Pb to the Tl/Si(111)1 × 1 surface at RT produces a homogeneous Tl-Pb alloy having periodicity. Its atomic arrangement (Figure 2f) looks qualitatively akin that of the Bi-Na alloy (Figure 1f), namely the first adsorbate, Bi or Tl, form the chained trimers arranged into a honeycomb network, while atoms of the second adsorbate, Na or Pb, occupy the centers of the honeycomb units. The difference is that Pb atom residence and Tl trimers centers in the Tl-Pb alloy are in the T1 (on-top) sites while the positions of Na atoms and Bi trimers in the Bi-Na alloy are the T4 sites. Both Bi-Na and Tl-Pb alloys are confined essentially in single atomic layers with Na being by ~0.62 Å higher than Bi and Pb by 0.37 Å higher than Tl.

Substantial structural changes due to Tl/Pb alloying, including, in particular, the change in the surface periodicity, 1 × 1 to , leads to changing the band structure. The most essential newly developed features are the two spin-split metallic surface-state bands (denoted by Σ1 (Σ1′) and Σ2 (Σ2′)) which are clearly seen in Figure 2g and h.

The calculated Fermi surfaces of the alloys reproduce nicely the contours observed experimentally (Figure 3). The main features of the alloy spin texture predicted by DFT calculations are also illustrated in Figure 3. The planar spin components show in-plane helicity, suggestive of a Rashba scenario. For momentum vectors along the directions, the spin is fully aligned in-plane and perpendicular to the momentum vector. The out-of-plane spin component undulates between positive and negative values along the contours according to the C3 symmetry of the surfaces.

Figure 3

Calculated band structure of the (left panel) and (right panel) including spin-orbit coupling illustrating the (a, e) in-plane and (b, f) out-of-plane polarization components. (c, g) Experimental and (d, h) calculated Fermi contours. Inset in (g) shows the fine structure within the outlined area around the point with a greater contrast. Arrows adjacent to the calculated contours and their length denote the in-plane spin component. The out-of-plane spin component is indicated by the colour with red and blue corresponding to the upward and downward directions, respectively. White colour indicates fully in-plane spin alignment. The tic values in parentheses in (d) refer to the pocket contours around the .

In particular, the (Bi, Na)/Si(111) band structure contains spin-split Σ1 (Σ1′) band which shows up in the Fermi map (Figure 3c and d) as two hexagon-shaped contours with corners pointing in the directions (where the splitting of the band is maximal, momentum splitting Δk∥ = 0.044 Å−1 and energy splitting ΔEF = 210 meV). Besides, the Σ1 (Σ1′) band forms hole pockets having the shape of smoothed triangles around the points. The out-of-plane spin component for pocket contours is negligible (albeit finite and different for neighboring pockets), the in-plane component dominates and the counterclockwise spin helicity is the same for all six pocket contours.

As for the (Tl, Pb)/Si(111), it displays the two spin-split metallic bands, Σ1 (Σ1′) and Σ2 (Σ2′). In the Fermi map of the Σ1 (Σ1′), the outer contour has almost round shape, while the inner contour is a hexagon with corners pointing in the directions. The maximal splitting for the Σ1 (Σ1′), Δk∥ = 0.038 Å−1 and ΔEF = 250 meV, is along direction. The other spin-split metallic band Σ2 (Σ2′) shows up as hexagonal contours with corners along the direction, i.e. rotated by 30° compared to the inner contour of the Σ1 band. The maximal splitting for the Σ2 (Σ2′), Δk∥ = 0.050 Å−1 and ΔEF = 140 meV, is along direction.

Comparing the Fermi contours of the present alloys with those found for spin-split surface states of the reported phases with hexagonal symmetry, one can notice certain common features. For example, the contour with two hexagons having corners in the directions (as for the (Bi, Na)/Si(111) Σ1 (Σ1′) and (Tl, Pb)/Si(111) Σ2 (Σ2′) bands) was also observed at the surface3. The contour of (Tl, Pb)/Si(111) Σ1 (Σ1′) band with hexagon having corners in the directions inside a circle is akin that reported for 4 and 19.

Discussion

The shown examples of 2D alloying provide a hint for elucidating the pathways for converting the originally semiconducting surface into the metallic one. Note that both original surfaces, and Tl/Si(111)1 × 1, contain 1.0 ML of metal adsorbate. Taking into account that the surface lattice constant of Si(111) is 3.84 Å, such a density of adsorbate atoms is not sufficient to ensure the overlapping of their electron wave functions. This is clearly seen in the calculated electron density maps in Figure 4a and c: electron density is localized either around Bi trimers or single Tl atoms, respectively. Note that it is typical for metal-induced Si(111) reconstructions where the atomic layers with metallic properties develop usually at metal coverage beyond 1.0 ML (e.g., with 1.2 and 2.4 ML In20 or with 1.33 ML Pb21). Another possibility is realized when the metal layer of 1.0 ML is compressed into the interconnected network due to Si atoms incorporated within the same reconstructed layer (e.g., 2223 or 2425, each containing 1.0 ML metal and 1.0 ML Si). One can see that adding 1/3 ML Na to the produces a similar effect: Na adsorption causes originally isolated Bi trimers to rotate and grow in size until they form a chained-trimer structure. As a result, electron density associated with Bi atoms becomes arranged into the interconnected network (Figure 4b) responsible for appearance of the surface metallic properties. The (Tl, Pb)/Si(111) system looks seemingly similar, as originally isolated Tl atoms also form chained-trimer structure. However, the electron density map (Figure 4d) shows that actually the Tl-associated electron density is interconnected through that of the Pb atoms. Hence, appearance of the surface metallic properties is controlled here just by increasing metal-atom coverage at the surface.

Figure 4

Calculated electron density maps for (a) , (b) , (c) Tl/Si(111)1 × 1, and (d) surfaces. The multi-colour (blue-yellow) scale for electron density is used to highlight the interconnected electron density network at metallic alloy surfaces.

Thus, our results demonstrate that formation of 2D alloys containing a metal with a strong spin-orbit coupling allows to obtain spin-split metallic surface-state bands on silicon. However, finding an appropriate metal pair is actually a challenging task. For example, while Na is a suitable component for making 2D alloy with Bi, it appears that the other alkali metals, Li and Cs, are not alloyed with Bi and grow as 3D islands on surface. Nevertheless, the list of the prospective alloys displaying variety of spin-related properties is believed to be very vast.

Methods

Sample preparation

The STM and ARPES experiments were performed in an ultra-high-vacuum Omicron MULTIPROBE system with a base pressure better than ~2.0 × 10−10 Torr. Atomically-clean Si(111)7 × 7 surfaces were prepared in situ by flashing to 1280°C after the samples were first outgassed at 600°C for several hours. Pristine and Tl/Si(111)1 × 1 reconstructions were formed by depositing 1 ML of the corresponding species onto Si(111)7 × 7 surface held at ~450°C and ~300°C, respectively. Sodium was deposited from a well-outgassed commercial SAES chromate dispenser, lead from a heated Pb-stuffed Mo tube.

STM

STM images were acquired using Omicron variable-temperature STM-XA operating in a constant-current mode. Electrochemically-etched W tips and mechanically cut PtIr tips were used as STM probes after annealing in vacuum.

ARPES

ARPES measurements were conducted in the ultrahigh vacuum chamber Omicron MULTIPROBE using VG Scienta R3000 electron analyzer and high-flux He discharge lamp (hν = 21.2 eV) with toroidal-grating monochromator as a light source.

Ab initio random structure search

Atomic structure of the new alloyed and reconstructions were elucidated using ab initio random structure search technique26 after it had been successfully tested for the known and Tl/Si(111)1 × 1 structures.

First-principles calculations

Our calculations were based on DFT as implemented in the Vienna ab initio simulation package VASP2728, using a planewave basis and the projector-augmented wave approach29 for describing the electron-ion interaction. The generalized gradient approximation (GGA) of Perdew, Burke, and Ernzerhof (PBE)30 has been used for the exchange correlation (XC) potential. The Hamiltonian contains the scalar relativistic corrections, and the spin-orbit interaction (SOI) was taken into account by the second variation method as has been implemented in VASP by Kresse and Lebacq31. To simulate the reconstructions we use a slab consisting of 12 bilayers (BL). Hydrogen atoms were used to passivate the Si dangling bonds at the bottom of the slab. Both bulk Si lattice constant and the atomic positions within the three BLs of the slab were optimized including SOI self-consistently. The silicon atoms of deeper layers were kept fixed at the bulk crystalline positions.

Author Contributions

D.V.G., L.V.B. and A.V.M. carried out ARPES and STM under the support of A.A.Y. and A.Y.T.. J.P.C., C.M.W., S.V.E. and E.V.C. carried out the theoretical calculation. D.V.G., Y.L.W., A.V.Z. and A.A.S. analyzed the data and wrote the manuscript with input from J.P.C., M.Y.L., S.V.E. and E.V.C. and conceived and coordinated the project.