Bi monocrystal formation on InAs(111)A and B substrates L. Nicola¨ı,1 J.-M. Mariot,2, 3 U. Djukic,4 W. Wang,5 O. Heckmann,4, 6 M. C. Richter,4, 6 J. Kanski,7 M. Leandersson,8 J. Sadowski,9, 10 T. Balasubramanian,8 I. Vobornik,11 J. Fujii,12 J. Braun,13 H. Ebert,13 J. Min´ar,1 and K. Hricovini4, 6
arXiv:1807.00306v1 [cond-mat.mtrl-sci] 1 Jul 2018
1
New Technologies-Research Center, University of West Bohemia, Univerzitni 8, 306 14 Pilsen, Czech Republic 2 Sorbonne Universit´e, CNRS (UMR 7614), Laboratoire de Chimie Physique–Mati`ere et Rayonnement, 4 place Jussieu, 75252 Paris Cedex 05, France 3 Synchrotron SOLEIL, L’Orme des Merisiers, Saint-Aubin, BP 48, 91192 Gif-sur-Yvette, France 4 Laboratoire de Physique des Mat´eriaux et des Surfaces, Universit´e de Cergy-Pontoise, 5 mail Gay-Lussac, 95031 Cergy-Pontoise, France 5 Link¨ oping University, Department of Physics, Chemistry and Biology / Surface and Semiconductor Physics, 581 83 Link¨ oping, Sweden 6 DRF, IRAMIS, SPEC – CNRS/UMR 3680, Bˆ at. 772, L’Orme des Merisiers, CEA Saclay, 91191 Gif-sur-Yvette Cedex, France 7 Chalmers University of Technology, Department of Physics, 412 96 Gothenburg, Sweden 8 MAX IV Laboratory, Lund University, P.O. Box 118, 221 00 Lund, Sweden 9 Lund University, MAX-lab, P.O. Box 118, 221 00 Lund, Sweden 10 Institute of Physics, Polish Academy of Sciences, al. Lotnik´ ow 32/46, 02-668 Warszawa, Poland 11 Istituto Officina dei Materiali, TASC Laboratory, CNR, Area Science Park - Basovizza, Strada Statale 14, km 163.5, 34149 Trieste, Italy 12 Istituto Officina dei Materiali, TASC Laboratory, CNR, 34014 Trieste, Italy 13 Ludwig-Maximilians-Universit¨ at M¨ unchen, Department Chemie, Butenandtstraße 11, 81377 M¨ unchen, Germany (Dated: July 3, 2018) The growth of Bi films deposited on both A and B faces of InAs(1 1 1) has been investigated by low-energy electron diffraction, scanning tunneling microscopy, and photoelectron spectroscopy using synchrotron radiation. The changes upon Bi deposition of the In 4d and Bi 5d5/2 photoelectron signals allow to get a comprehensive picture of the Bi/InAs(1 1 1) interface. From the initial stages the Bi growth on the A face (In-terminated InAs) is epitaxial, contrary to that on the B face (Asterminated InAs) that proceeds via the formation of islands. Angle-resolved photoelectron spectra show that the electronic structure of a ≈ 10 BL deposit on the A face is identical to that of bulk Bi, while more than ≈ 30 BL are needed for the B face. Both bulk and surface states are well accounted for by fully relativistic ab initio spin-resolved photoemission calculations. PACS numbers: 79.60.-i, 71.18.+y, 71.20.-b, 73.20.At
I.
INTRODUCTION
. A number of recently discovered materials, such as topological insulators (TIs) or Weyl semimentals, share a fundamental similarity. Rather than fermions obeying the usual Schr¨ odinger Hamiltonian, their low-energy fermionic excitations behave as massless Dirac particles. This emergent property of fermions in condensed matter systems defines the unifying framework for a class of materials we call “Dirac material”. The emergence of Dirac excitations is controlled by symmetries of the material, such as time-reversal, gauge, and spin-orbit coupling. The family of TIs, discovered a decade ago [1–3], has led to flourishing of novel and fascinating physics. The simplicity of their surface states, the robustness of both topological properties and surface metallicity under external perturbations, and the prediction of novel quantized states arising from the peculiar coupling between magnetic and electric fields make TIs a perfect venue for the next generation electronic and spintronic devices. Focusing on the spintronics aspect, one of the possibilities to obtain spin-polarized electronic currents is to use
the spin splitting in momentum space of surface states that results from the breaking of inversion symmetry at a surface (Rashba effect) [4]. This has recently stimulated studies of the electronic structure of Bi and Bi-based compounds, because of the large spin-orbit interaction in this heavy group-V semimetal (see, e.g., Ref. [5]). Along this line many investigations of the changes in the electronic structure of Bi thin films deposited on semiconductors of with respect to bulk Bi have been performed, e.g., Si(1 1 1) [6 and 7], SiC(0 0 0 1) [8], and highly-oriented pyrolytic graphite [9]. Such studies on Bi thin films are also of importance to the field of topological insulators [1, 2, 10–12]. As far as thin layers are concerned, it has been theoretically predicted that the Bi(1 1 1) surface could behave as a topological insulator [13–15]. Key factors leading to this property are tensile strength due to the mismatch between the lattice of the Bi deposit and that of the substrate, as well as the formation of chemical bonds at the Bi/substrate interface. Thus the investigation of such interfaces is a great relevance to explore the existence of new physical properties and, if any, to open ways to engineer the relevant interfaces. On the experimental side,
2 III.
RESULTS AND DISCUSSION
A.
Pristine InAs(1 1 1) surfaces
Although a (2 × 2) LEED pattern is observed for InAs(1 1 1)A [Fig. 1 (a)], the InAs(1 1 1)B surface is unreconstructed (not shown); this is due to the presence of In vacancies on the In-terminated A surface (see Refs. [24– 26], and the references cited therein).
(a)
<
In
4.28 Å
In
(e)
(d)
In
Bi >
In
Bi Bi
II.
In
EXPERIMENTAL AND COMPUTATIONAL METHODS
The experiments were performed at the APE beamline [17] of the Elettra synchrotron radiation facility (Trieste, Italy) and at the I3 beamline [18] of the MAX III synchrotron (MAX-lab, Lund, Sweden). The substrates were cut from 0.5 mm thick n-doped (3 × 1018 cm−3 ) InAs(1 1 1) wafers (Wafer Technology Ltd., UK) polished on both sides. After degassing, both InAs(1 1 1)A and B surfaces were prepared (simultaneously) by repeated cycles of ion bombardment (Ar+ , 600 eV) and annealing at 400◦ C until a sharp (2 × 2) or (1 × 1) LEED pattern for side A or side B, respectively, was obtained (see, e.g., Ref. [19] and the references cited therein). Bismuth was deposited from a Knudsen cell at a rate of about 0.5 BL min−1 under a pressure lower than 4 × 10−10 Torr. On the theoretical side, band structure calculations were performed using the SPR-KKR package [20] that is based on the Korringa–Kohn–Rostoker method which uses the Green’s function formalism within the multiple scattering theory. The package is based on the Dirac equation, thus fundamentally containing all relativistic effects, such as spin-orbit coupling. The calculations were performed within the local density and atomic sphere approximations for the bulk potential, that was used as an input to simulate the ARPES spectra within the one-step model of photoemission [21]. The model of Rundgren– Malmstr¨ om [22] was used to connect the bulk potential with the surface. The relaxation of the surface was also taken into account using the structural data given in Ref. [23].
>
In
Bi
4.55 Å
(c)
(b)
<
III-V substrates have not been explored up to now. We report here on a study of the electronic structure of Bi films deposited on both A (In-terminated) and B (Asterminated) faces of InAs(1 1 1). Although some results on the growth of Bi on InAs(1 1 1)B have been reported [16], to our knowledge that of Bi on InAs(1 1 1)A has not yet been investigated. We have followed the structural and electronic changes that occur during the deposition of Bi on both faces of InAs(1 1 1) using low-energy electron diffraction (LEED), scanning tunneling microscopy (STM), and core-level as well as valence-band photoelectron spectroscopy (PES) using synchrotron radiation. The growth mode is found to be dependent on the substrate side, an epitaxial growth being evidenced on the A side. As a consequence, ultra-thin Bi films deposited on InAs(1 1 1)A have an the electronic structure showing all the characteristics of bulk Bi. This does not occur when deposition is done on the B face because then the Bi growth proceeds via the formation of islands. The angle-resolved PES (ARPES) data on Bi/InAs(1 1 1)A are well explained by ab initio fully relativistic multiple scattering theory in the framework of the density functional theory.
Bi Bi
In
(f) (g)
(f)
50 nm
4.1 Å 4.1 Å
FIG. 1. a) and b) LEED patterns of pristine InAs(1 1 1)A showing the (2 × 2) reconstruction and of a 10 BL Bi deposit InAs(1 1 1)A (taken at 34 and 53 eV, respectively). Structures of (In-terminated) InAs(1 1 1)A c), Bi(1 1 1) d), and a possible arrangement of the Bi deposit on InAs e). (f and g) STM images of the InAs(1 1 1)A substrate covered by 10 BL Bi.
The In 4d PES spectra for both the A side and the B side of InAs(1 1 1) are in agreement with those reported previously [16, 27, and 28]. They are given in Fig. 2 to help comparison in the following discussion with those of the Bi-covered substrates. The origin of a ≈ 250 meV shift towards the high binding energy (BE) side of the main peak of the In 4d PES spectrum for the A side with respect to the B side is due to an effective charge transfer perpendicular to the surface resulting in a polarization of the As-In bi-layers of InAs(111)B (1 × 1) [29]. It was recognize early [27] that appart from the bulk component (denoted as B) In 4d core levels contain as well a surface contribution (denoted as S). Our experi-
3
Normalized amplitude -20 20
(b)
InAs(111) B In 4d
(a)
InAs(111) A In 4d
B
S B
S
-18 18
Binding energy (eV)
-16 16
18 -18
16 -16
Binding energy (eV)
FIG. 2. In 4d PES spectra of clean InAs(1 1 1): (a) A side; (b) B side. Raw data are shown as dots; the lines are the fit to the data and the components involved in the fit.
mental data are well reproduced fitting these two components by a Voigt profile with Gaussian and Lorentzian full width at half maximum of 0.18 eV and 0.35 eV, respectively. One can also note that the intensity of the S component is larger for the In-terminated A side than for the B side. The B–S energy shift is found to be 0.3 eV for the A side, in good agreement with [27].
B.
Bi/InAs(1 1 1): core levels
At the early stage of the deposition, the LEED pattern of the InAs(1 1 1)A (2 × 2) substrate [Fig. 1(a)] starts to blur and relatively broad spots showing to the formation of Bi(1 1 1) layers appear. A well-defined sharp spots of Bi(1 1 1) (1 × 1) being observed for a Bi coverage greater than approximately 10 BL, see Fig. 1(b). This is consistent with the fact that the In-terminated InAs(1 1 1)A surface exhibits an hexagonal pattern of In atoms separated by ≈ 4.28 ˚ A lattice [Fig. 1(c)] and the Bi(1 1 1) an honeycomb pattern in which the distance between a Bi atom and its second in-plane neighbor being ≈ 4.55 ˚ A [Fig. 1(d)]. This ≈ 6% mismatch is small enough to allow an epitaxial growth of Bi(1 1 1) on the InAs(1 1 1) substrate, as shown for instance in the manner depicted in Fig. 1(e) which is compatible with the LEED patterns in Figs. 1(a) and (b). The evolution of the intensity of the Bi 5d and In 4d PES signals as a function of the Bi evaporation duration on InAs(1 1 1)A (not shown) is made of a series of linear segments with break points, separated by a common evaporation duration, where a change in slope occurs (see top panel of Fig. 3). This is the well-known signature of a layer-by-layer (Frank–van der Merwe) growth mode. This is consistent with both the LEED patterns and the scanning tunneling microscopy (STM) images [Fig. 1(d)] of these surfaces which terraces with the six-fold symmetry [Fig. 1(c)] characteristic of the hexagonal arrangements of the Bi atoms on (1 1 1) terraces [Fig. 1(d)] [19]. On the InAs(1 1 1)B surface the plot of the intensity of the Bi 5d and In 4d PES signals as a function of
the Bi deposition duration shows a break in linear behavior only for the first 1 BL Bi. This suggests that the Bi growth on the B side of InAs(1 1 1) proceeds via a Stranski–Krastanow mode, contrary to what happens on InAs(1 1 1)A. This is confirmed by the fact that an amount of ∼ 50 BL Bi is necessary to clear the substrate In 4d PES signal, i.e., to ensure the coalescence of the Bi three-dimensional islands that build on the surface after one Bi BL has been deposited on the B side of InAs. Changes in the electronic structure of the Bi/InAs(1 1 1)A interface can be detected at the early stage of the Bi deposition in the In 4d PES spectrum (recorded at hν = 45 eV) [Fig. 3(a)]. For 0.2 BL coverage the whole spectrum shifts relative to that of the pristine substrate by ≈ 200 meV towards the high BE side; this shift reaches ≈ 300 meV for a coverage of ≈ 1 BL and stabilizes to this value for higher coverages. This shift is accompanied by a change in the overall shape of the In 4d spin-orbit doublet. Both effects are due to the progressive disappearance of the surface S component present in the spectrum of pristine InAs(1 1 1)A upon Bi deposition, see Fig. 2. The substrate In 4d PES signal disappears totally for a Bi coverage of ≈ 10 BL. The bottom spectrum in Fig. 3(a) presents the subtraction between 2 BL and 0.2 BL Bi coverages. This difference spectrum is very similar to that corresponding to the S component in the In 4d PES spectrum of pristine InAs(1 1 1)A. This is thus a confirmation of the disappearance of the S component at the early stage of Bi deposition on InAs(1 1 1)A. Note that we did not subtract the In 4d PES for 0.2 BL Bi coverage because of the band bending effect induced by the deposition of the first Bi atoms. Finally it is interesting to note that the subtraction reveals a faint structure at ≈ 17.4 eV BE, which corresponds to the BE value expected when Bi is chemically bound to In. Such a behavior is also observed for the Bi/GaAs(1 1 1)A and B interfaces [30 and 31] as well as for the Bi/InAs(1 0 0) interface [32 and 33]. We now turn to the Bi 5d5/2 PES signal [Fig. 3(b)]. When the amount of Bi deposited is increased its shape converges to that of bulk Bi at ≈ 10 BL of Bi. To shed light on the origin of the changes occurring at low Bi coverage we show at the bottom of Fig. 3(b) the difference spectrum obtained by subtracting the Bi 5d5/2 PES spectrum for 15 BL from that for 0.2 BL Bi coverage, removing in such a way the bulk component. This difference spectrum is made of two components: the one at low BE, of much higher intensity, can be attributed to Bi bound to In and that at high BE to Bi bound to As (see Ref. [16]). The latter can be present despite the InAs(1 1 1)A substrate is In-terminated because some defects due to surface preparation (ion bombardment and annealing) might allow As atoms to be present on the surface. The subtraction clearly indicates that the apparent change in the position of the Bi 5d5/2 PES upon Bi deposition, namely for 0.2 BL Bi coverage, arises from the intensity variation/decrease of the PES component
4 (b)
Bi 5d5/2
Normalized amplitude
Bi BL 10
Bi BL
10
2
2
1
1
0.2
0.2
0
(a)
In 4d
300 meV
Normalized amplitude
(a)
In 4d
S
(b)
Bi 5d5/2
In-Bi In-In
Bi BL 13 27
(27)-(0.2)
(13)-(0.2)
(0)-(0.2)
(0.2)-(15)
(2)-(0.2) -19 19
-18 18
-17 17
-16 -25 25 16
Binding energy (eV) (c) In 4d
-24 24
-23 23
Binding energy (eV) (d) Bi 5d5/2
Normalized amplitude
13
13
17 -17
16 25 -16-25
-24 24
23 -23
Binding energy (eV)
2.5
2.5
0.2
0.2
0 (0)-(0.2) -19 19
18 -18
Binding energy (eV)
FIG. 4. (Color online) PES difference spectra between Bicovered (in BL) and pristine InAs(1 1 1)B: (a) In 4d; (b) Bi 5d5/2 .
Bi BL 27
Bi BL
19 -19
(0.2)-(27)
-18 18
-17 17
Binding energy (eV)
-16 25 16 -25
-24 24
-23 23
Binding energy (eV)
FIG. 3. Changes in core level PES spectra upon Bi deposition (in BL) on InAs(1 1 1)A and B (top and bottom panels, respectively): (a) In 4d spectra. The difference spectrum between In 4d PES signals for 2.0 and 0.2 BL Bi coverage is shown (in green) at the bottom of the figure. (b) Bi 5d5/2 spectra. The difference spectrum between Bi 5d5/2 PES signals for 0.2 and 15 BL Bi coverage is shown (in green) at the bottom of the figure. (c) In 4d; the difference spectrum between In 4d PES signals for 0 and 0.2 BL Bi deposition is shown (in green) at the bottom of the figure; (d) Bi 5d5/2 spectra, the difference spectrum between Bi 5d5/2 PES signals for 0.2 and 27 BL Bi deposition is shown (in green) at the bottom of the figure.
corresponding to the Bi–In bond, similarly to what is observed for the In 4d PES signal upon Bi deposition. The evolution of In 4d and Bi 5d5/2 PES spectra on InAs(1 1 1)B is quite different from InAs(1 1 1)A [Fig. 3(c) and (d)]. The presence of Bi atoms on the B surface does not induce a shift as for InAs(1 1 1)A. This allows a subtraction of the In 4d spectrum of the bare substrate from that corresponding to 0.2 BL Bi deposition [”(0)-(0.2)”, at the bottom of Fig. 3(c)]. Similarly to InAs(1 1 1)A, the S component disappears at very early stages of Bi deposition. For thicker Bi layers a new low BE feature appears in the spectra. A better insight of its evolution is shown in Fig. 3(a) where we plot only difference spectra. The surface component, ”(0)-(0.2)” from Fig. 3(c), is sh for own for a reference. Subtracting In 4d spectrum corresponding to 0.2 BL Bi and containing only the bulk component, from that of 13 BL and 27 BL Bi depositions, middle and up-
per spectrum, respectively (only In 4d5/2 component is shown), reveals two new structures that can be attributed to Bi-In (higher BE) and In-In (lower BE) components. Such a feature has been previously suggested [16] but not resolved. Here, the In-In component clearly appears for higher Bi coverages. We believe that its formation is due to the fact that the Bi-As bond is stronger as compared to the Bi-In one, as deduced from BEs shifts [16]. At higher coverages Bi-As bonds being formed, releasing free In atoms on the surface. Similar chemical behavior has been observed when annealing Bi layer on top of InAs(111) crystal [19]. Evolution of Bi 5d5/2 PES spectra [Fig. 3(d)] are in agreement with this chemical interpretation. As in the case of InAs(111)A side we show the difference spectrum [bottom in Fig. 3(d)] obtained by subtracting the Bi 5d5/2 PES spectrum for 0.2 BL Bi coverage from that of 13 BL coverage. In the difference spectrum we find again two components, the one at high BE has now a larger intensity which is expected for the As-terminated B side. The presence of the low BE component, attributed to Bi-In, can be justified because the annealing during the sample preparation eliminates preferentially As atoms leaving In-terminated patches on the surface. This image is consistent with further evolution of Bi 5d5/2 spectra. For 2.5 BL Bi coverage, as all Bi-As bonds are saturated, the Bi-In is a dominant component. Upon increasing deposition the FHWM of Bi 5d5/2 decreases and reaches that of bulk Bi well above 27 BL [Fig. 3(b)], when Bi islands coalesce in larger terraces, in clear contrast with InAs(111)A surface. We suggest that the presence of strong Bi-As bonds prevents epitaxial growth at very initial stages of Bi deposition. Reminding that for such high depositions the In 4d PES signal is faint but still present, this is in line with the formation of 3D Bi islands during the growth. This situation is similar to that encountered on GaAs where Bi bands have been reported for a 2 BL deposit [30]. Indeed the STM images show the presence of such islands on
5 0 0
1
1
1
-2 2
-2Kinetic Energy [eV] 2
-2Kinetic Energy [eV] 2
3
3
3
0 0.0
0.2 0.2
k (Å-1)
-0.4 -0.4k^-1 0.4 0.4
-0.2 -0.2
0 0.0
0.2 0.2
(b)
0.0 0
0.2
0.2 0.2
-0.2
0.4 0.4
-0.4
0.4
0.4 0.4
k (Å-1)
FIG. 5. (Color online) ARPES spectra along the Γ–K line of InAs(1 1 1)A covered by (b) 2 BL Bi and (c) 12 BL Bi. (a) Normal emission EDC corresponding to (b) (red) and (c) (blue). Spectra were recorded at hν = 20 eV
InAs(1 1 1)B [26]. We deposited Bi films as well on MBE prepared surfaces. Interestingly, in spite of their higher quality, as compared to the ion-bombardment/annealing preparation, Bi growth on the InAs(111)A face is worse. In MBE systems, it is impossible to suppress instantaneously the As background pressure after the InAs substrate growth. So, even the perfect In-terminated A-face is covered by a fraction of monolayer of As preventing the good epitaxy. Paradoxically, Bi monocrystals prepared on bombarded/annealed A-face are of much higher quality.
C.
0.0 0
Binding energy (eV)
Binding energy (eV)
-0.4 -0.2 Intensity (a.u.) -0.4k^-1 -0.2
(a)
Binding energy (eV)
0 0 (c) (b)
(b)
0 0 (a)
Bi/InAs(1 1 1) valence band PES
The ARPES spectra recorded at hν = 20 eV along the Γ–K line for Bi coverages of 2 BL and 12 BL are shown in Fig. 5. A Bi coverage of 2 BL is sufficient for the ARPES spectrum to show main bands characteristic of bulk Bi [Fig. 5(b)], in particular Bi states near the Fermi level (EF ) start to appear in the gap of InAs. We observed that the bands sharpen progressively upon increasing Bi coverage, the whole ARPES spectrum becoming identical to that of bulk Bi [5] for 12 BL coverage [Fig. 5(c)]. Normal emission spectra in Fig. 5(a), extracted from Fig. 5(b) and (c), show namely strong intensity increase of the surface state resonance at binding energy of 2 eV that characterize high quality Bi crystal [5]. This is in complete consistency with the Bi epitaxial growth mode on the A side, as deduced from the analysis of the core level PES spectra (see Sec. III B). The measured spectrum was as well compared to calculations [see Fig. 6]. Bi band dispersion calculations were carried out with the SPR-KKR package. The following lattice parameters were used: a = 4,5332 ˚ A and c = 11,7967 ˚ A. The measured spectrum [Fig. 6(a)] in the vicinity of the Γ point along the Γ–K direction are well reproduced by the corresponding theoretical calculation [Fig. 6(b)]. The theoretical simulation was performed for a Bi semi-infinite crystal with the [111] direction. The following ARPES calculations are based on the one-step model of photo-emission applied to a bulk potential; this
-0.2
0.2k^-1
0.1 -0.1
0
0.0
-0.1 0.1
k (Å-1)
0.6Kinetic Energy [eV] 0.6 -0.2 0.2 0.1 0.2 -0.2 -0.1
0
0.0
0.1
-0.1
-0.6 0.6 0.2
-0.2
k (Å-1)
FIG. 6. (Color online) ARPES spectrum along the Γ–K direction of a Bi(111) monocrystal grown on an InAs(111)-A substrate (a) measured and (b) computed with the SPR-KKR package. The photon energies used were hν = 20 eV for the measurement and hν = 22 eV for the calculation with, for both cases, a circularly polarized light.
approximation appears to be sufficient to reproduce the measured spectra. As seen from Fig. 6 the calculated bands are in good agreement with the calculation. The electronic bands show a high sensitivity to the helicity of the light used to probe them as shown by the PES dichroic spectrum recorded at hν = 20 eV in the Γ–M direction [Fig. 7(a)]. The theoretical prediction at hν = 22 eV [Fig. 7(b)] reproduces nicely the experimental result. These results can be compared to the spin polarization to display the correlation between the observed circular dichroism and the fundamental spin polarization. In addition the theoretical analysis of the contribution to the dichroic signal of initial state having selected angular momentum is shown in Figs. 7(d–f). When only p [Fig. 7(d)] or d [Fig. 7(e)] initial states are considered not all the bands experimentally observed are recovered. The remaining states, circled in blue, are regained when both p and d symmetries are considered as seen in Fig. 7(d)] where only the s orbitals were deactivated. This suggests a partial hybridization of the p and d bands. At this stage it deserves to be reminded that circular dichroism in ARPES has been recently proposed as an alternative to spin-resolved ARPES because it is a much less experimentally experimental technique (see, e.g., Refs. [34–36]). However, extensive combined experimental and theoretical studies of ARPES using circularly polarized light have demonstrated that contrasted experimental findings between spin-resolved and dichroic ARPES arise from the spin dependence of the relativistic dipole matrix elements that also depend strongly on the final states reached in the PES, i.e., on the photon energy inducing the photoelectron. This underlines the fact that the final state involved the photoemission process play a role that must be more influential on the result than the spin texture of the initial state. Actually, sign inversions are observed rather regularly when increasing the photon energy by several eV and coincide only accidentally with the spin polarization [see Fig. 7(c)].
6 Measured Dichroism
Calculated Dichroism max.
a)
min.
0.6 0.6 -4
-2
0
2
θ (deg)
4
0.3
min.
0.6 -4
6
Binding energy (eV)
Binding energy (eV)
0.4 0.4
0 θ (deg)
0.3
0.6
4
-4
0 θ (deg)
4
0 θ (deg)
4
Different orbital l influences
e)
only p initial states
0.0
0.3
0.6 -4
0 θ (deg)
4
f)
only d initial states
0.0
Binding energy (eV)
Binding energy (eV)
0.0
Binding energy (eV)
d)
CONCLUSIONS
0.0
Intensity (arbitrary units)
Intensity (arbitrary units)
0.2 0.2
IV.
c)
0.0
0.0 0.0
Kinetic energy (eV)
Spin Polarisation max.
b)
0.3
no s initial states
0.3
We show that InAs(111) substrates prepared by two to three cycles of successive ion sputtering and annealing present distinctive results for the Bi growth depending on the atomic surface termination. Core-level analysis confirm that Bi growth is epitaxial on the InAs(111)-A surface even for very thin films. Deposition of 10 bilayers of Bi on InAs(111)-A resulted in a Bi monocrystal of very high quality as attested by reproduction of the band dispersion by theoretical calculations. Moreover, a close comparison of circular dichroism spectra between theory and experiment highly suggests a pd-orbitals hybridization.
0.6
0.6 -4
0 θ (deg)
4
-4
Figure 4.9: Experimentally measured (a) and theoretically calculated (b) circular dichroism spectra along the K direction with ~⌫ = 20 eV and ~⌫ = 22 eV , respectively. States that required an enhancement of the bulk effects are indicated by arrows in FIG. 7. (Color ARPES of (c) Bi(1 1 1):the relative (b). Further theoretical studies showonline) the corresponding computeddichroic perpendicularspectra spin polarisation and depict roles of each the orbitals l (d,e,f), without by the Fermi bands in blue are reproduced (a)of Experiment (hν correction = 20 eV) for function. 20 BLTheof Bicircled deposited on when both p and d orbitals are activated.
ACKNOWLEDGMENTS
InAs(1 1 1)A. (b) Theory (hν = 22 eV). (c) Calculated spin polarization. (d, e, and f) Theoretical analysis of the conThe research leading to these results has received fundof Rundgren-Malmström [153]angular is used, which connects the bulk muffin-tin potential the surface tribution of the momentum in the initial statewith (reing from the European Community’s Seventh Frameasymptotic one by a third order polynomial. The surface barrier value is here selected as 0.24, which spectively, ` = 1, 2, and in absence of ` = 0) to the dichroic work Programme (FP7/2007-2013) under grant agreedefines the bulk/surface interface.6 Finally, a relaxation of the surface was taken into account by using the response. ment n◦ 312284. parameters given by the work of Mönig et al. [79].
The measurements presented below for the dichroism study were performed at the APE beamline of Elettra. The measured spectra in the vicinity of the
point along the K direction are well reproduced by
the corresponding theoretical calculation (see respectively Fig. 4.8 (b,c)). The theoretical simulation was performed for a Bi semi-infinite crystal with the [111] direction. However, in order to obtain bands arising 1
M.bulkZ.effects, Hasan and C. L. Kane, “Colloquium : calculated Topological from primary the number of layers over which the photocurrent is N_LAYER had insulators,” Mod. Phys. 82,the3045 (2010). to be increased from 15 to 25.Rev. This consequently confirms good quality of the monocrystal of Bi(111). 2
11
Y. Ando, insulator J. light Phys. These electronic bands“Topological revealed a high sensitivity to thematerials,” polarisation of the used Soc. to probe them, Jpn.circular 82, 102001 (2013). giving a peculiar dichroism, as shown by the spectrum of Fig. 4.9 (a). On the theoretical side, 3
Xiao-Liang Qi spin andpolarisation Shou-Cheng “Topological insuall bands are reproduced. The in the Zhang, M direction (perpendicular to the analysed lators and Rev. circular Mod. dichroism Phys. 83, 1057 K direction), displays the superconductors,” correlation between the observed and the fundamental
12
(2011). spin-polarisation. Within the one-step model, it is possible to activate/deactivate the different orbitals: s, p 4
E. only I. Rashba, “Properties of initial semiconductors with an exand d. Keeping p (Fig. 4.9 (d)) or d (Fig. 4.9 (e)) states does not allow to recover all bands.7 tremum loop.in1.blue, Cyclotron combinational The remaining states, circled are regainedand if p and d orbitals are bothresonance activated: thisin is done in 6 7
a magnetic field perpendicular to the plane of the loop,”
an appropriate surface barrier is needed for the surface states not to be shifted in binding energy Sov. State 2,is 1109 Tela the influence fromPhys. the orbitalSolid s for circular dichroism negligible(1960), here (smaller [Fiz. by severalTverd. orders of magnitude)
5
6
7
8
9
10
(Leningrad) 2, 1224 (1960)]. Ph. Hofmann, “The surfaces of bismuth: Structural and electronic properties,” Prog. Surf. Sci. 81, 191 (2006). A. Takayama, T. Sato, S. Souma, and T. Takahashi, “Giant out-of-plane spin component and the asymmetry of spin polarization in surface Rashba states of bismuth thin film,” Phys. Rev. Lett. 106, 166401 (2011). A. Takayama, T. Sato, S. Souma, T. Oguchi, and T. Takahashi, “Tunable spin polarization in bismuth ultrathin film on Si(1 1 1),” Nano Lett. 12, 1776 (2012). H. Huang, S. L. Wong, Y. Wang, J.-T. Sun, X. Gao, and A. T. S. Wee, “Scanning tunneling microscope and photoemission spectroscopy investigations of bismuth on epitaxial graphene on SiC(0 0 0 1),” J. Phys. Chem. C 118, 24995 (2014). P. J. Kowalczyk, O. Mahapatra, D. N. McCarthy, W. Kozlowski, Z. Klusek, and S. A. Brown, “STM and XPS investigations of bismuth islands on HOPG,” Surf. Sci. 605, 659 (2011). C.-L. Gao, D. Qian, C.-H. Liu, Jia J.-F., and F. Liu, “Topological edge states and electronic structures of a 2D
13
14 71
15
16
17
18
topological insulator: Single-bilayer Bi(1 1 1),” Chin. Phys. B 22, 67304 (2013). I. K. Drozdov, A. Alexandradinata, S. Jeon, S. Nadj-Perge, H. Ji, R. J. Cava, A. B. Bernevig, and A. Yazdani, “Onedimensional topological edge states of bismuth bilayers,” Nat. Phys. 10, 664 (2014). A. Takayama, T. Sato, S. Souma, T. Oguchi, and T. Takahashi, “One-dimensional edge states with giant spin splitting in a bismuth thin film,” Phys. Rev. Lett. 114, 066402 (2015). L. Fu and C. L. Kane, “Topological insulators with inversion symmetry,” Phys. Rev. B 76, 045302 (2007). Z. Liu, C.-X. Liu, Y.-S. Wu, W.-H. Duan, F. Liu, and J. Wu, “Stable nontrivial Z2 topology in ultrathin Bi(1 1 1) films: A first-principles study,” Phys. Rev. Lett. 107, 136805 (2011). Z.-Q. Huang, F.-C. Chuang, C.-H. Hsu, Y.-T. Liu, H.-R. Chang, H. Lin, and A. Bansil, “Nontrivial topological electronic structures in a single Bi(1 1 1) bilayer on different substrates: A first-principles study,” Phys. Rev. B 88, 165301 (2013). K. Szamota-Leandersson, M. Leandersson, M. G¨ othelid, and U. O. Karlsson, “Correlated development of a (2 × 2) reconstruction and a charge accumulation layer on the InAs(1 1 1)-Bi surface,” Surf. Sci. 605, 12 (2011). G. Panaccione et al., “Advanced photoelectric effect experiment beamline at Elettra: A surface science laboratory coupled with synchrotron radiation,” Rev. Sci. Instrum. 80, 043105 (2009). T. Balasubramanian et al., in SRI 2009, 10th International Conference on Synchrotron Radiation Instrumentation, AIP Conf. Proc. Series, Vol. 1234, edited by R. Gar-
7
19
20
21
22
23
24
25
26
27
28
rett, I. Gentle, K. Nugent, and S. Wilkins (2010) p. 661. M. C. Richter et al., “Bi atoms mobility-driven circular domains at the Bi/InAs(1 1 1) interface,” Surf. Sci. 651, 147 (2016). H. Ebert, D K¨ odderitzsch, and J Min´ ar, “Calculating condensed matter properties using the KKR-Green’s function method–recent developments and applications,” Rep. Prog. Phys. 74, 096501 (2011). J. Braun, “The theory of angle-resolved ultraviolet photoemission and its applications to ordered materials,” Rep. Prog. Phys. 59, 1267 (1996). J. Rundgren and G. Malmstr¨ om, “Transmission and reflection of low-energy electrons at the surface barrier of a metal,” J. Phys. C: Solid State Phys. 10, 4671 (1977). H. M¨ onig, J. Sun, Yu. M. Koroteev, G. Bihlmayer, J. Wells, E. V. Chulkov, K. Pohl, and Ph. Hofmann, “Structure of the (1 1 1) surface of bismuth: LEED analysis and firstprinciples calculations,” Phys. Rev. B 72, 085410 (2005). A. Taguchi and K. Kanisawa, “Stable reconstruction and adsorbates of InAs(1 1 1)A surface,” Appl. Surf. Sci. 252, 5263 (2006). A. Taguchi, “First-principles investigations of surface reconstructions of an InAs(1 1 1)B surface,” J. Cryst. Growth 278, 468 (2005). E. Hilner, E. Lundgren, and A. Mikkelsen, “Surface structure and morphology of InAs(1 1 1)B with/without gold nanoparticles annealed under arsenic or atomic hydrogen flux,” Surf. Sci. 604, 354 (2010). ¨ Olsson, M. C. H˚ C. B. M. Andersson, L. O. akansson, L. Ilver, U. O. Karlsson, and J. Kanski, “Core level photoemission from (1 1 1)-type InAs surfaces,” J. Phys. IV France 4-C9, 209 (1994). K. Szamota-Leandersson, M. G¨ othelid, O. Tjernberg, and U. O. Karlsson, “Interaction between oxygen and
29
30
31
32
33
34
35
36
InAs(1 1 1) surfaces, influence of the electron accumulation layer,” Appl. Surf. Sci. 212-213, 589 (2003). S. Mankefors, P. O. Nilsson, and J. Kanski, “Semiconductor polar surfaces: mechanisms of the stability of nonreconstructed III-V (1 1 1) surfaces,” Surf. Sci. 443, L1049 (1999). C. McGinley, A. A. Cafolla, B. Murphy, D. Teehan, and P. Moriarty, “The interaction of bismuth with the GaAs(1 1 1)/B surface,” Appl. Surf. Sci. 152, 169 (1999). C. McGinley, A. A. Cafolla, E. McLoughlin, B. Murphy, D. Teehan, P. Moriarty, and D. A. Woolf, “Core-level photoemission study of the Bi-GaAs(1 1 1)/A interface,” Appl. Surf. Sci. 158, 292 (2000). K. Szamota-Leandersson, M. Leandersson, P. Palmgren, M. G¨ othelid, and U. O. Karlsson, “Electronic structure of bismuth terminated InAs(1 0 0),” Surf. Sci. 603, 190 (2009). M. Ahola-Tuomi, M. P. J. Punkkinen, P. Laukkanen, M. Kuzmin, J. L˚ ang, K. Schulte, A. Pietzsch, R. E. Per¨ al¨ a, N. R¨ as¨ anen, and I. J. V¨ ayrynen, “Properties of selfassembled Bi nanolines on InAs(1 0 0) studied by core-level and valence-band photoemission, and first-principles calculations,” Phys. Rev. B 83, 245401 (2011). Y. H. Wang, D. Hsieh, D. Pilon, L. Fu, D. R. Gardner, Y. S. Lee, and N. Gedik, “Observation of a warped helical spin texture in Bi2 Se3 from circular dichroism angle-resolved photoemission spectroscopy,” Phys. Rev. Lett. 107, 207602 (2011). C.-H. Park and S. G. Louie, “Spin polarization of photoelectrons from topological insulators,” Phys. Rev. Lett. 109, 097601 (2012). Y. Wang and N. Gedik, “Circular dichroism in angleresolved photoemission spectroscopy of topological insulators,” Phys. Status Solidi RRL 7, 64 (2013).