Anomalous Rashba Effect of Bi Thin Film Studied by Spin-Resolved ARPES Anomalous Rashba Effect of Bi Thin Film Studied by Spin-Resolved ARPES

The Rashba effect is a momentum‐dependent splitting of spin bands in two‐dimen‐ sional systems such as surface, interface and heterostructure. The effect is caused by broken space‐inversion symmetry and spin‐orbit coupling and allows to manipulate and generate the spin by the electric fields, that is, without the magnetic field. It means that the devices applied to the Rashba effect have many advantages. Bismuth is known as a promising candidate to investigate the surface Rashba effect, and the spin struc‐ ture of Bi surface has also been intensively discussed. However, it is unclear to what extent the so far believed simple vortical spin structure is adequate. To understand the surface properties of the Rashba system is particularly important when utilizing the Rashba effect to the spintronic devices, since it is desirable to control the spin polariza‐ tion when developing new types of devices. In this chapter, we report that the surface spin states of the Bi thin film exhibit unusual characteristics unlike the conventional Rashba splitting by using a spin‐ and angle‐resolved photoemission spectroscopy measurement.


Introduction
As we know and use, spintronic devices to use a spin-polarized electrons have actualized. The magnetic storage technology uses giant magneto-resistance [1]. A more advanced approach is to control spin-polarized electrons without the aid of a ferromagnetism nor to apply the magnetic field [2]. Spin-orbit coupling (SOC) makes it possible to generate and manipulate spin-polarized electrons only by the electric field, since the electric field acts on a moving charge carrier as an effective magnetic field. Thus, it is regarded as an essential ingredient for further development of next-generation spintronic devices such as the spin-field-effect transistor [3]. In nonmagnetic solids, the electronic states with opposite spin have the same energy (Kramers degeneracy) because of the time-reversal and the space-inversion symmetries (TRS and SIS). In the strong SOC environment with the broken space-inversion symmetry (typically at the surface or interface), the energy band splits in the momentum (k) space (Rashba effect [4]), leading to a spin helical structure of surface bands. This Rashba effect leads to the vortical spin structure of surface bands where the spin vector points parallel to the surface and perpendicular to the measured momentum. To be more specific, here we postulate the model of two-dimensional free electron gas at the surface. In the surface, as shown in Figure 1a, there is asymmetry of the potential in the direction perpendicular to the two-dimensional plane [∇V = (0, 0, E z )]. When an electron moves with momentum (p), the ∇V × p term acts as an effective magnetic field B eff , which is orthogonal to p (=ħk) and ∇V. As a result, the electron spin is quantized along the direction perpendicular to k in the surface plane. The energy of free electron gives the following, which is well known as the Rashba effect, where α R is a so-called Rashba parameter. Figure 1b shows the band dispersion of this model. As described above, two-dimensional (2D) system with strong SOC has provided a useful platform for realizing novel quantum phenomena applicable to advanced spintronic devices [2,5,6]. To discuss and understand Rashba effect, spin-and angle-resolved photoemission spectroscopy (spin-resolved ARPES) is a powerful experimental technique, which can simultaneously determine all key quantum parameters of electrons in solids, that is, momentum, energy and spin. Furthermore, by high-resolution measurements, we will be able to discuss the spin-related physical phenomena not only qualitatively but also quantitatively, which would certainly lead to the deeper understanding of the condensed matter physics. Elucidation of the electron spin is very important to understand physical property in solid state and its surface as well as the possible applications to spintronic devices.
The first observation of the surface Rashba effect by ARPES is the Au(1 1 1) surface [7]. After that, various materials such as group-V semimetals and their alloy surfaces [8][9][10][11][12][13][14][15][16], as well as heavy-atom adsorbed semiconductor surfaces [17][18][19][20][21][22] and so on, are studied. Among them, the group-V semimetal bismuth (Bi) is a prime candidate to investigate the surface Rashba effect and many experiments and theoretical calculations studied in order to clarify the fundamental properties of the Rashba effect [8,[10][11][12][13][14]. However, in previous researches, although the band structure and Fermi surfaces of Bi distinctly show strong anisotropy, the spin structure of Bi was argued by assuming an isotropic two-circular Fermi-surface model like Au(1 1 1) [7,[20][21][22]. The reason is because the energy and momentum resolutions of the previous spin-resolved ARPES machine are insufficient. To understand whole aspect of the Rashba effect, it is necessary to clarify that the spin structure of Bi consists with the conventional Rashba model or not. Moreover, in the film, electronic states originating from the bulk Bi are quantized, and they are connected to the surface bands continuously. So it would be necessary to take into account the relationship between bulk and surface states. On the other hand, intensive attempts have been made to extend the investigations on 2D Rashba systems to quasi one-dimensional (1D) system like artificially grown nanowires and quantum wires, because of the merits in downsizing of devices.1D Rashba effect in utilization of vicinal surfaces such as in Au chains on vicinal Si [23] and the vicinal Bi surface is also reported [24]. In such a case, breaking the TRS by applying magnetic field or adding magnetic impurities would create an energy gap at the Kramers point, and when the chemical potential is tuned to be located in the spin-orbit gap, the dissipation-less spin transport and the quantized conductance [25,26] may be realized (Figure 1c). However, the Rashba effect in edge state is not known because the signal from the edge is extremely faint. To understand the proposed novel properties of a true 1D system, we may be able to apply them to advanced spintronic devices.
In this chapter, we introduce the electronic structure of Bi thin film to elucidate the details of the Rashba effect by utilizing the high-resolution spin-resolved ARPES spectrometer equipped with a highly efficient mini-Mott detector. We show three novel Rashba effects of Bi thin film: (i) anisotropic Rashba effect from momentum-dependent measurement [27], (ii) the interface Rashba effect between metal-semiconductor from thickness-dependent [28] and (iii) 1D Rashba effect of edge state [29]. The present finding provides a useful platform to study the Rashba effect and at the same time opens a pathway to utilize the novel properties to advanced spintronic devices. Figure 2 shows a schematic diagram of the ultra-high-resolution spin-resolved ARPES spectrometer with a highly efficient mini-Mott detector [30]. This spectrometer consists of mainly four parts: (i) a photoemission measurement system including a hemispherical electron energy analyzer and an ultra-high-vacuum measurement chamber, (ii) a spin-detection system based on a mini-Mott detector, (iii) an intense xenon/helium plasma discharge lamp and (iv) a surface chamber to prepare the thin-film samples. We explain the detail of each part.

Spin-and angle-resolved photoemission spectromete
We have improved a MBS-A1 electron energy analyzer to achieve both spin-resolved and regular (non-spin-resolved) ARPES measurement. The spectrometer has two detectors: one is a multichannel plate for ARPES measurement, and the other is mini-Mott detector for spinresolved ARPES. To determine the three-dimensional spin polarization, an electron deflector has been placed between the analyzer and the Mott detector. The Mott detector observes the spin polarization of essentially two independent axes by using four channeltrons, enabling us to determine the in-plane and out-of-plane spin component. The scattering efficiency of the Mott detector is as high as 2.3 × 10 -2 . The optical system consists of helium (He) and xenon (Xe) plasma discharge lamps and a monochromator with the gratings, in which we can select photon energy if necessary. In this study, we used one of the Xe I lines (hν = 8.437 eV) to excite photoelectrons. In order to fabricate high-quality samples, a surface chamber has been constructed and is connected to the spin-resolved ARPRS spectrometer. The spectrometer achieves the energy resolutions of 0.9 and 8 meV for non-spin-resolved and spin-resolved modes, respectively. Figure 3 shows a schematic view of the surface chamber. A surface chamber contains heating systems for a semiconductor sample, a few kinds of dispensers, a quartz crystal microbalance and a low-energy electron diffraction (LEED) system for checking the quality of the sample surface. The vacuum of the surface chamber is basically kept 1 × 10 -10 Torr to prepare a high-quality sample surface. It allows the in situ preparation of the sample and its transfer to the spin-resolved spectrometer. It is particularly useful for elucidating the electronic states of the samples containing a clean well-ordered surfaces required for the accurate ARPES measurements. Thus, we can prepare the Bi thin-film samples by the evaporation of Bi on semiconducting substrate. We expect that the performance of the spectrometer is demonstrated by the observation of a clear Rashba splitting of the surface states in Bi. The energy and momentum resolutions for the regular (spin-integrated) ARPES were 5-20 meV and 0.3°, and for spin-resolved ARPES measurements were 40 meV and 3°, respectively. The Sherman function value was set at 0.07. The measurement temperature was 300 and 30 K.

Sample fabrication
To get a high-quality Bi(1 1 1) thin film, we prepare a clean surface of the Si(1 1 1) substrate. We use a commercially available Si wafer (n-type, As-doped: 0.001-0.005 Ω cm, Sb-doped: 0.01-0.02 Ω cm), and the surface of Si forms native SiO 2 in the air. So we must remove it by electrical heating in high vacuum. Figure 4 displays the design of a holder for the electrical heating system (Figure 4). The holder is made of molybdenum, and a main holder is insulated from a subholder by the alumina. We can control a current and temperature with 5 mA and 1°C, respectively. The sample size is 13 × 3 mm, and we have selected four different kinds of crystal orientations on the basis of orientation flat as shown in Figure 5. The sample geometry has been confirmed by the brightness symmetry of the LEED spots. The most stable structure of the Si(1 1 1) surface is the 7 × 7 reconstructed surface, as shown in Figure 6a [31,32]. Now, we describe the method to pre-  pare a Si(1 1 1)-7 × 7 reconstructed surface. First, the Si wafers were outgassed for more than 12 h below 750°C. After the outgassing enough, we have carried out flash annealing. Figure 6b shows a flash-annealing process to get a Si(1 11)-7 × 7 reconstructed surface; sample was (i) heated at 750°C to 1050-1200°C for a few seconds, (ii) keep temperature to maximum for 5 s, (iii) cooled down to 850°C for a few seconds and (iv) cooled to 750°C in 30 s [33]. We have repeated this cycle, and all of the above processes should be performed under the ultra-high vacuum of ~1 × 10 -10 Torr. As shown in Figure 7a, we have obtained the LEED pattern of the well-ordered 7 × 7 surface.
Next, Bi atoms are evaporated on Si substrate, which is called as a molecular beam epitaxy (MBE). Bi atoms are deposited at room temperature on Si(1 1 1)-7 × 7 reconstructed surface. Then, the Bi thin film was annealed at 150°C. The deposition rate is estimated by the quartz oscillator thickness monitor, and the film thickness was controlled by varying the deposition time with keeping the constant deposition rate. We can also estimate the film thickness from the energy position of the quantum well states (QWSs) in the ARPES spectra [28]. It is noted that 1 bilayer (BL) Bi is defined as 1.14 × 10 15 atoms/cm 2 , and the thickness is 0.39 nm [13]. In this study, we prepared several thickness sample . After the deposition of Bi atoms, the LEED pattern shows the 1 × 1 surface structure as shown in Figure 7b and c, and the intensity of the LEED spot has threefold symmetry.
When the sample has a multi-domain structure, the LEED pattern shows circular features surrounding the 1 × 1 spots (Figure 7d). We could repeatedly use one Si substrate by a flashing. Here we describe a structure of Bi thin film on Si substrate. Details of the thin-film growth process are shown in Figure 8a, which is reported in previous works such as STM and LEED [34]. Although the lattice constant of Bi (4.538 Å) is very different from Si (5.43 Å), it is possible to fabricate the Bi/Si thin film due to the existence of disordered layer called "wetting layer" between Bi and Si substrate. As shown in Figure 8b, the structural transition from {0 1 2} direction to (1 1 1) direction more than 8.4 ML suddenly takes place upon Bi deposition.

Anisotropic Rashba effect of Bi thin film
At first, we have performed normal (non-spin-resolved) ARPES measurement of the Bi thin film in order to check the sample quality and geometry, since even a subtle misalignment of the sample orientation would cause a significant error in determining the spin polarization. Figure 9a and b shows the band dispersion along the Γ ¯ ¯ M line and Fermi surface of Bi/Si (1 1 1). Areas with hatched lines in Figure 9a are bulk band projection. We distinctly see that several bands cross E F and three kinds Fermi surfaces exist; a hexagonal Fermi surface centered at the Γ ¯ point (S 1 ), surrounding elongated pockets (S 2 ) and the ellipsoidal pocket (S 3 ) near the ¯ M point. Judged from the band dispersion in Figure 9a, the S 1 and S 3 are attributed to the electron pockets, while the S 2 to a hole pocket, all of which arise from the spin-split surface states [8,[10][11][12][13][14].
This section focuses on the spin structure of the S 2 band. In a conventional 2D Rashba model, the in-plane spin has a vortical structure and isotropic magnitude as denoted by black arrows in Figure 10a. In Figure 10b and e, we display the near-E F spin-resolved energy distribution curves (EDCs) for the in-plane (y/x) and out-of-plane (z) spin components measured in various k regions, A-J as shown in Figure 10a. As displayed in Figure 10b and e, the spin-resolved  EDC mainly consists of two components: a slope-like feature, which rapidly increases its intensity at the binding energy (E B ) higher than 0.15 eV corresponding to the tail of the QWS, and a weaker broad feature at E F -0.1 eV assigned as the S 2 band.
First, we take a closer look at the spin polarization in regions A-D. In regions A and B, as seen in Figure 10b, the in-plane spin polarization of the S 2 band is dominated by the up spin. On the other hand, the down spin is barely superior in regions C and D. It means that the in-plane spin structure is qualitatively consistent with the Rashba picture [13][14]. However, it seems that the spin polarization is markedly suppressed in C and D. In fact, the magnitude of the spin polarization along the y direction |P y | is 0.5-0.7 in regions A and B, while it is 0.2-0.3 in regions C and D. This is unexpected since |P y | should keep the same value across the Γ ¯ point in the normal Rashba picture. Then, we focus on the out-of-plane (z) spin component. We immediately notice that there exists a sizable up-spin polarization, while theoretically predicted out-of-plane spin polarization for simple isotropic Rashba system is zero. More surprisingly, the magnitude of the z-axis spin polarization |P z | (0.4-0.7) is as large as that of |P y |.
We have also observed a finite in-plane and out-of-plane spin polarization in regions E-H. As shown in Figure 10c, the estimated maximum |P y | of the S 2 band in regions E and F The data in regions I and J are also the same trend as shown in Figure 10d. The difference in |P y | across the Γ ¯ point is also recognized by comparing the EDCs between regions I and J. In addition, the value of |P z |(~zero) is a good agreement in regions I and J where the up-and down-spin components almost overlap with each other because of the cancellation of two opposite spins. To see if a similar trend is observed along another high-symmetry line Γ ¯ ¯ K , we demonstrate in Figure 10e the spin-resolved EDCs measured in regions K-N. It is apparent that the sign of P z (also P x ) in region K (L) is the same as that in the region N (M), indicating that P z (P x ) does not switch the sign across the Γ ¯ ¯ K line. Figure 11 shows a schematical view of the spin polarization vectors of the S 2 band from Figure 10. We symmetrized the data by taking into account the threefold crystal symmetry due to the presence of a second bismuth layer. The in-plane spin component has a vortical structure, but the magnitude of the spin polarization is perpendicular to k, called here P θ . P z has a large component and switches the sign by every 60° step. These features would lead to the periodic oscillation of P θ and P z , unlike the general Rashba SOC where P θ = const and P z = 0. Now, we discuss the origin of anomalous Rashba effect of Bi thin film. In the conventional Rashba effect, the in-plane spin polarization of S 2 is symmetry with respect to the Γ ¯ point. In this study, the TRS is not broken by some magnetic impurities because the Bi and Si are nonmagnetic materials. Another possibility of causing the broken TRS is the local surface conditions or the final-state effect [35,36]. However, this possibility might be also unlikely by the reproducibility of data and comparing with previous studies. Here it is noted that the LEED pattern of Bi/Si(1 1 1) in Figure 7 shows the threefold symmetry due to the bilayer-Bi crystal structure, and we conjecture that the origin of the in-plane spin asymmetry would be related to the crystal symmetry at present. Next, we discuss the out-of-plane spin polarization. Unlike in-plane spin polarization, the threefold symmetry in the out-of-plane spin component is reported, and sign-switching behavior is in good agreement with the previous studies in Bi/Ag(1 1 1) and Bi 1-x Sb x [15,37]. The spin structure of Bi surface state is also similar to the spin-resolved ARPES experiments on the hexagonally warped Dirac-cone Fermi surface of the topological insulator as well as the prediction of the kÂ·p theory [38,39]. On the other hand, the absolute value of the out-of-plane spin polarization obtained by our experiment is 40-70%, while the theoretical value of previous studies is a few percent. In order to clarify the spin structure of Bi Rashba effect, we need to evolve theoretical and experimental studies.

Rashba effect at interface of a Bi Thin film on Si(1 1 1)
In this section, we focus on the spin structure around the ¯ M point of the Bi/Si(1 1 1) thin film. As shown in Figure 12a, the band dispersion of the QWSs dramatically alters around the ¯ M point as a function of the film thickness d, and we prepared five different thickness film for d = 8, 10, 15, 20 and 40 BL. All of the samples, the bottom energy of an electron-like dispersion is located at ~30 meV below E F . This band is well separated from the hole-like band at higher E B arising from the QWSs that originate from the confinement of the wave function in the direction perpendicular to the surface (z direction). There is no evidence for the Kramers degeneracy of the S 3 band or the QWSs at the ¯ M point (k y = 0). The absence of the Kramers degeneracy at the ¯ M point has also been observed in the previous research of the Bi, and it is explained by the hybridization of the surface states and the QWS [9,12]. As clearly visible, we identify the QWSs at E B = 0.1-0.25 eV for d = 40 BL. Upon decreasing d, the top of the highest QWS shifts downward, and the energy separation of each QWS becomes wider, enabling us to unambiguously identify the overall energy dispersion of the individual QWS in the band dispersion plots.
In conventional Rashba picture, the in-plane spin polarization of the S 3 band is dominated by the up spin independently of a thickness at the cut 1 shown in Figure 12b. Figure 12c shows the spin-resolved EDCs and spin polarization of P y for various d values of 8-40 BL. All of the spin-resolved EDCs consist of the tail of the QWSs at E B > 0.1 eV and the S 3 electron band at ~50 meV. As seen in Figure 12c, the in-plane spin-resolved EDCs of the S 3 band for 10-40 BL are dominated by the up spin, and this indicates that the in-plane spin direction is qualitatively consistent with the normal Rashba picture. However, our result clearly demonstrates that the observed P y value strongly depends on d. The obtained P y value at d = 40 BL has a maximum at ~0.7, then gradually reduces at d = 15 BL (P y ~0.3) and finally reaches ~0 at d = 8 BL. It is noted that the much smaller spin polarization in the above-E F region for d = 15 BL as compared to that for d = 40 BL might be due to much weaker peak weight in the original EDC and an intrinsic suppression of the spin polarization, which would lead to a relative enhancement of the background weight and less statistical reliability of the states above E F .
We discuss the physical mechanism behind the unusual thickness dependence of P. The coupling between the Bi film and the Si(1 1 1) substrate is fairly weak due to the presence of a disordered wetting layer at the interface, suggesting that the Bi film is nearly freestanding [34]. Namely, the Bi/Si interface is also broken SIS, and it can be thought of as another (bottom) surface as with vacuum-side (top) surface. These Rashba states of top and bottom surface should have opposite spin directions as illustrated in Figure 13. In case of bulk Bi, which is thick enough, the top and bottom Rashba states do not interfere with each other. On reducing thickness, the wave functions of the two Rashba states overlap and hybridize, which means that the up-and down-spin states merge and the P observed by spin-resolved ARPES would become small. This picture is a very simple model but explains in a nice way the thickness dependence of P. We also found that the decay length of the surface-state wave function is around at least 20 BL (~80 Å) since the reduction in experimental P is already started in 20 BL (Figure 12). A thin film of Bi 2 Se 3 topological insulator also exhibits a similar decay length, as inferred from the experimental fact that the hybridization gap starts to open at six quintuple layers (~60 Å) [40]. Our result also suggests that the Rashba effect of Bi/Si interface is useful for the spintronic devices because the interface is generally more stable than the surface.

1D edge state with Rashba effect in Bi thin film
As observed by the atomic force microscopy of our Bi thin film (Figure 14), triangular-shaped Bi BL islands with typically ~0.1 μm edge length are formed on the top surface of the Bi thin film as reported previously [34], and the edge of each island is perpendicular to the Γ ¯ ¯ M direction considering sample geometry; namely, the edge runs along the Γ ¯ ¯ K direction in the k space. According to previous experimental and theoretical studies, there are no bulk and surface states near the ¯ K point, and we might observe some electronic state if edge state exists along the Γ ¯ ¯ K line.
We demonstrate the band dispersion of Bi thin film along the Γ ¯ ¯ K line for d = 15 BL in Figure 15. With a careful look at the region between Γ ¯ ¯ K line (yellow rectangle in Figure 15a), one finds unexpected faint intensity displaying a finite energy dispersion, as better illustrated with enhanced color contrast in the inset. To establish the energy dispersion of this unexpected feature, we have measured the ARPES data along several cuts in the surface Brillouin zone (Figure 15b). A careful band searching with the enhanced intensity scale (Figure 15c) shows the band dispersion with a characteristic x-shape along cuts 1-3. The intersection of . This demonstrates that the x-shaped dispersion has a 1D character along the k y direction. Indeed, the observed band structure along cut 4 (perpendicular to cuts 1-3) shows no obvious dispersion, confirming the 1D nature. To experimentally clarify the spin-split nature of the edge band, we have performed a spin-resolved ARPES experiment. As shown in Figure 15d which plots the spin-resolved EDCs at a representative k y point (marked by a pink line in Figure 15c), we clearly find a difference between the up-and down-spin spectra in both the in-plane and out-of-plane components.
As for the origin of the unexpected 1D state, we have taken into account various possibilities such as the mixture of domains with different film thickness, surface reconstruction, slight isolation of the topmost Bi bilayer and surface stacking faults. A most natural and convincing explanation is that it originates from the edge states of Bi bilayer. To further strengthen our conclusion, we have carried out first-principles electronic band structure calculations for a specific crystal structure (Figure 16). Electronic band structure calculations were carried out by means of a first-principles density functional theory approach with the all-electron full-potential linearized augmented-plane-wave method in the scalar relativistic scheme. The spin-orbit coupling was included as the second variation in the self-consistent-field iterations. Thin-film systems were simulated by adopting periodic slab models with sufficiently thick vacuum layer. In this model, Bi atoms in 1D allay are removed from the topmost Bi 1BL (Figure 16a) so as to reproduce the infinitely long edge structure along the y direction. Assumption of such an idealized model crystal is turned out to be sufficient for reasonably simulating the edge band structure. As shown in Figure 16c, the Brillouin zone for this model crystal structure has a rectangular shape. The vertical length of Brillouin zone is the same as the ¯ M ¯ M interval since the size of unit cell along y-axis is the same for the edge structure and the Bi thin film. On the other hand, the horizontal Brillouin zone length has no important physical role because the unit-cell length for x-axis was simply chosen for the sake of calculations. Thus, the high-symmetry point ¯ Y now becomes the time-reversal-invariant momentum and is located exactly at the horizontal projection of the ¯ M point (k y = 0.69 Å -1 ), which actually coincides with the intersection of the x-shaped band. Figure 16c displays the calculated band dispersion along the Γ ¯ ¯ Y direction of Bi thin film for d = 10 BL. We identify two prominent dispersive bands, which cross E F and have the edge-state origin. Moreover, the edge bands show the Rashba spin splitting due to the strong spin-orbit coupling at the edge, as evidenced by the degeneracy at the ¯ Y point. We found that the overall spin-vector direction in the experiment, that is, sign of spin polarization for each component, is also consistent with the calculation, while the absolute magnitude of the spin polarization in the experiment (~0.2) is much smaller than that of the calculation (~0.26 to 0.64) in both the in-plane and out-of-plane components, likely due to a finite contribution from the angle-integrated-type background in the ARPES spectra and the spin-orbit entanglement effect [41]. Finally, we have estimated the Rashba parameter α R by numerical simulation for the 1D parabolic band with the SOC and obtained α R = 0.80 ± 0.05 eVÅ (Figure 17). This value is much larger than that for the 2D surface state (0.56 eVÅ) [11,15], and the difference could be explained in terms of the presence of an in-plane potential gradient [15,42] at the edge in addition to the out-of-plane component, which already exists in the 2D film, as supported by observation of the out-of-plane spin polarization as large as the in-plane counterpart. Recently, STM study on the edge states of Bi crystal and thin film [43,44]. The STM study reported that one of two types of edges, where Bi atoms are terminated at close to vacuum, which corresponds to our calculation, has a 1D character. The band dispersion of the edge state observed in ARPES shows a good agreement with the calculated bands, suggesting that both ARPES and STM observe the same edge state. It is true that 1D electronic state that observed our ARPES measurement is localized at edge.
The edge state in Bi would be also useful in realizing novel physical properties and new spintronic devices.

Conclusions
We have demonstrated anomalous Rashba effect of Bi thin film on Si(1 1 1) by using spinresolved ARPES. Major findings of the present work are the following three features: (i) the surface Rashba states of Bi exhibit the asymmetric in-plane spin polarization and the giant out-of-plane spin polarization, (ii) the spin polarization of the surface states is reduced on decreasing thickness and (iii) 1D band dispersion from the edge state of Bi islands on Si(1 1 1) exhibits large Rashba effect. These observed peculiar spin states are not explained in terms of Figure 17. Experimental band dispersion for cuts 1-3 extracted from the peak position of the EDCs in Figure 15c, compared with the numerical simulation for a simple 1D Rashba splitting (dashed black curve).
the conventional Rashba effect, and these results open a pathway for realizing exotic physical properties at the strong SOC systems.