Identification of helicity-dependent photocurrents from topological surface states in Bi2Se3 gated by ionic liquid.

Dirac-like surface states on surfaces of topological insulators have a chiral spin structure with spin locked to momentum, which is interesting in physics and may also have important applications in spintronics. In this work, by measuring the tunable helicity-dependent photocurrent (HDP), we present an identification of the HDP from the Dirac-like surface states at room temperature. It turns out that the total HDP has two components, one from the Dirac-like surface states, and the other from the surface accumulation layer. These two components have opposite directions. The clear gate tuning of the electron density as well as the HDP signal indicates that the surface band bending and resulted surface accumulation are successfully modulated by the applied ionic liquid gate, which provides a promising way to the study of the Dirac-like surface states and also potential applications in spintronic devices.

Three dimensional topological insulators are insulators with time-reversal symmetry and strong spin-orbit coupling (SOC)12. It is a new quantum phase of matter with insulating bulk states and Dirac-like helical surface states whose spin is perpendicularly locked to its momentum234. The helical surface states reside in the bulk insulating gap, and are protected from backscattering by time-reversal symmetry34. That makes the topological insulators potentially promising to serve as a platform for spintronics and quantum computing applications, based on the crucial spin coherence. Bi2Se3 is currently most applied one of these fascinating materials, with topologically non-trivial energy gap about 0.3 eV and only a single Dirac cone in the surface energy spectrum5. The static electronic properties of the surface states have been comprehensively imaged by angle-resolved photoemission spectroscopy (ARPES) and other surface sensitive techniques5678. Transport measurements at low temperature have also been carried out to study the new physics of these novel massless two-dimensional electron systems91011. However, the reported transport studies have been bothered by its inevitably metallic bulk1011. ARPES studies have also shown that the cleaved sample surfaces progressively become more electron doped by either gas contamination, or formation/migration of defects, e.g. vacancies1213. Although it has been reported that the chemical potential in Bi2Se3 could be tuned by surface deposits, bulk doping, and back gates, these methods are nevertheless cumbersome or restricted101114.

Spin-polarized helicity-dependent photocurrent (HDP) generated by circularly polarized light in the helical Dirac cone (HDC) is a useful technique to study the novel properties of topological surface states15. It had been studied theoretically as well as experimentally1617. Different from and advantaged than transport measurements, the contribution from the bulk states to HDP can be neglected, since the bulk states are spin degenerate517 However, clear Rashba spin-splitting of the conduction band in samples cleaved in air or exposure to water vapor has been reported by several groups through ARPES measurements1819. Furthermore, quantized sub-bands existed in the surface-electron-accumulation layer (SEAL) at the surface of the samples. Thus, one could not directly make the conclusion that the HDP is generated in HDC only. Therefore, it is very important to find a convenient and direct way to modulate the SEAL, and thereby identify the HDP truly stemmed from the HDC.

Here we report a study of the topological surface states in HDC by measuring the HDP generated by illumination of circularly polarized light. The SEAL in the samples was modulated by ionic liquid gate (ILG)202122. The clear gate-voltage dependence of the density of electrons and HDP signal shows that the band bending as well as the Rashba spin-orbit splitting of SEAL are successfully modulated by the ILG. These results demonstrate that the measured HDP has two components with different directions, one of which indeed stems from the HDC while the other from the SEAL. This work provides a promising way to the study of the Dirac-like surface states and also potential applications in spintronic devices.

Results

Photocurrents with different planes of incidence were measured, as shown in Figure 1b–1e. In the measurement, the total photocurrent j can be fitted by the formula: where j is the magnitude of the HDP, j the photocurrents induced by linear polarized light, and j0 the background photocurrent which is polarization-independent15. Only the HDP j is of concern in this study. The helicity dependence evidently indicates that j here is generated through a spin-related process, as in a spin-orbit coupling system the electrons with spin polarizations aligned or anti-aligned to the wave vector would be preferentially excited by left- or right-circularly polarized light, respectively. The results here show that j is extremely small when the light is obliquely incident in the y-z plane, in comparison with that in the x-z plane with the same incident angle. Although no clear signal can be observed when light is normally incident, j increases when the incident angle becomes larger in the x-z plane (Figure 1f).

Figure 1

(a) SEM image of the sample. (b) Sketch of HDP measurement at oblique incidence (45°) with y-z the plane of incidence. (c) Sketch of HDP measurement at oblique incidence (45°) with x-z the plane of incidence. (d) HDP measured at the setup of (b). (e) HDP measured at the setup of (c). In (d) and (e), solid line is a fit to equation (1), and dashed line indicates j. (f) HDP measured at different incident angles.

To clarify the origins of the measured HDP, we modulated the surface band-bending curvature by ILG. This is the key advance of this work along with the multiple advantages of ILG over the conventional solid gates202122. The ionic liquid together with the sample and the gate electrode forms an electric double-layer transistor geometry, where the interfaces could be regarded as nano-gap capacitors with huge capacitance, resulting in very high electron-density tunability23. Electric charges will accumulate at the interface when a gate voltage is applied. The electric charge accumulation will then induce an electric field at the interface to modulate the density of electrons and band-bending curvature in the SEAL. As sketched in Figure 2a, when a positive gate voltage is applied, electrons will accumulate at the interface, and the band-bending curvature will increase.

Figure 2

(a) Sketch of electric double-layer transistor geometry formed by using ILG. (b) Oblique drawing diagram of the setup for gate modulation of the HDP with ILG. Gate voltage is applied between the flake and a gate pad nearby. (c) The dependence of the HDP and surface electric field shift on the gate voltage applied.

Hall measurements with different gate voltage applied on ILG were carried out to have a quantitative evaluation of the electron accumulation near the surface. From the measurement, the density of electrons in the sample is tuned from 4.94 × 1013 cm−2 to 7.48 × 1013 cm−2 in the gate-voltage range varying from −0.9 V to 1.2 V. Further increase of the magnitude of the gate voltage will cause saturation20. Surface electric field induced by the electron accumulation could be obtained according to the Poisson equation E = qN/ε0ε, where E is the surface electric field, q the elementary charge, N the surface electron density, ε0 the vacuum dielectric permittivity, and ε ~ 113 the dielectric permittivity of Bi2Se32024. From the change of the density of electrons given above, the surface electric field is tuned as much as 0.43 MV/cm.

We then measured the HDP with different gate voltages based on the setup shown in Figure 2b. The gate voltage is applied between the Au gate pad and Bi2Se3 flake. Figure 2c shows the dependence of the HDP and surface electric field shift on the applied gate voltage measured in one of the samples. The HDP decreases from 154 pA to 5.4 pA with increasing gate voltage from −0.4 V to 1 V. The saturation of the HDP with further increasing of the gate voltage may come from the effective screening of the electric field or saturation of the ILG, as in the case of electron density mentioned earlier202122. Note that an interesting happening is the decreasing of the HDP with increasing gate voltage.

Discussion

First of all, we should check whether or not the band bending near the surface could be effectively tuned by the ILG. From the Hall measurements, the bulk density of electron in the sample at zero gate voltage is 4.8 × 1018 cm−3, taking the thickness of the sample into account. It is an upper limit since we assume that there is no contribution from the surface during the calculation. The Fermi level of the bulk is about 71 meV above the bottom of the conduction band (CB), where the effective mass of Bi2Se3 is 0.13 ~ 0.15m. Then, electrons in the HDC and the SEAL are considered. From the formula in Ref. 24, the depth of the SEAL is estimated to be 10 nm. The total surface density with electrons in HDC and SEAL all included could be expressed as , where x is the distance between the Fermi level and the bottom of CB at the surface in the unit of eV, assuming the SEAL forms a triangle quantum well. If we reasonably suppose that the difference of the bottom of CB between the surface and bulk is 0.1 eV, the range of x modulated in our measurements is from 0.053 eV to 0.253 eV1819. The tuning range is then 0.2 eV, which is in agreement with the tuned surface electric field considering the depth of the SEAL and screen effect. Thus, the surface band bending and the Rashba spin-splitting in the surface could indeed be effectively modulated by the ILG as shown in the references1819.

We now consider the possible source to the HDP j, i.e. the spin of electrons in the HDC, which is locked perpendicular to their linear momentum, and lies in the plane with slight warping51425. Generally, the HDP can be quantitatively described by , where j is the HDP, χ the second-rank pseudotensor, E0 the complex amplitude of the electric field of the incident light, and P the degree of circular polarization of the incident light15. In the presence of a [111] surface, the symmetry of the crystal is reduced from D. From the view of group theory, the nonzero items of the second-rank pseudotensor χ in HDC are χ and χ, with χ = −χ, which indicates that the HDP is always perpendicular to the plane of incidence. When x-z plane is taken as the plane of incidence, the HDP could be expressed as , where γ is the SOC coefficient, θ the incident angle. It is also obvious that j should become larger when incident angle of the light increases, due to the enhanced in-plane spin polarization component. As such, the results in Figure 1 (d) to (f) reveal that the HDP j indeed arises from the asymmetric optical excitation of the HDC.

To understand the interesting decrease shown in Figure 2c, contributions from other possible origins should be checked carefully. HDP has been observed on many systems with SOC2627. As a system with strong SOC, there could be additional sources of HDP in the Bi2Se3 flakes. Firstly, the photon energy of 1.17 eV for the used laser light is greater than the band gap of the Bi2Se3 flakes (~0.3 eV). Inter-band absorption is thus possible, which means that the bulk may have a prominent contribution to the measured photocurrent. However, the contribution from the bulk should be helicity-independent, since the symmetry of the bulk Bi2Se3 is D indicating that the items in the second-rank pseudotensor are all zero. Some recent papers revealed the existence of new unoccupied topological states located in bulk band gap above the Fermi level2829. This new topological states locates about 1.5 eV above the HDC with quite small band gap and small momentum. Although the new topological states are similar to the HDC, they actually contribute quite little to the HDP comparing with the HDC. Secondly, due to the high thermoelectric power of Bi2Se3, laser generated heat gradients in the sample might induce a bulk thermoelectric current in additional to the photocurrent17. We argue however that this bulk thermoelectric current should be small and have nothing to do with the circular polarization of the incident light. The light spot used in the experiments is much larger (~mm2) than the sample itself, which means that the light intensity can be considered to be even in the whole area of the sample. Thus no significant heat gradient could be generated. On the other hand, as discussed above, the absorption of the light mainly takes place in the bulk and is therefore helicity-independent also, so that the induced thermoelectric current should contribute to the background current only, but not to the HDP.

Recent reports demonstrated that the magnitude and sign of the surface dichroism of the insulating topological insulator may be oscillatory with the energy of the incident photon3031. However, the energy of the incident photons is the same in all of our measurements. According to the reports, surface dichroism depends significantly on the photon energy. What is more, the period of the oscillation indicated by this report is about 10 eV, which is almost a hundred times larger than the modulation range in our study. We thus could take the magnitude and the sign of the surface dichroism as the same in our analysis. The change of HDP from the HDC should mainly rely on the density of electron.

Another possible source of HDP is the Rashba spin-splitting states in the SEAL12131432. The Hamiltonian of the Rashba SOC can be expressed as , where α is the parameter of the SOC, the Pauli spin metrics, the electric field, and the electron momentum1520. It is noticed that the only nonzero component of the electric field in the SEAL is along z direction, which means that the spin of the electrons over there must lie in x-y plane. The nonzero items of the second-rank pseudotensor in the SEAL are also χ and χ. HDP stemmed from the SEAL would then have the same relationship with the plane of incidence and incident angle as that for the HDC. We thus cannot rule it out from our previous experimental results. The magnitude of HDP from SEAL depends on the strength of the SOC and electron density, hence it will increase (decrease) when a positive (negative) gate voltage is applied. The total HDP therefore has two components at oblique incidence and can be expressed as j = j + j, where j and j denote the HDP stemmed from the HDC, and SEAL, respectively.

As stated previously, the magnitude of the HDP from the SEAL (j) has a positive correlation with the applied gate voltage. Obviously, the magnitude of HDP from HDC should also become larger with increasing gate voltage, since the number of electrons in HDC increases. As such, the decreasing of HDP indicates that the two components must have different directions. Actually, by considering the band structure of Bi2Se3, the HDP from HDC is determined to be in opposite direction with the one from the SEAL. Thus, the measurements here make the identification that the total HDP comes mainly from the HDC in our measurements when no gate voltage is applied. The present technique thus provides an effective approach to determine which component accounts for the major part of the HDP by the direction of the HDP in a single measurement.

To summarize, we measured the HDP stemmed from the topological HDC, using generation by circularly polarized light in Bi2Se3. The observed HDP has two components. Besides the one from the HDC, the other one from the Rashba spin-splitting states in the SEAL is also significant. Fortunately, these two components in HDP have opposite directions, and can be distinguished from each other by the ILG modulation. HDP from the HDC dominate the total HDP in our measurements if no gate voltage is applied. In our study, the surface band-bending was successfully modulated by ILG, and the topological surface states as one of the sources was unambiguously identified from the total HDP. Both of the above results might be important for further study and possible applications in optical-spintronics33.

Methods

Materials

In our study, Bi2Se3 single crystals were grown by chemical vapor deposition (CVD)3435. The samples were fabricated by transferring Bi2Se3 flakes from the substrate onto a SiO2 surface. Suitable ones with smooth surface were picked up by using optical microscope and scanning electron microscope (SEM) (Figure 1a). The Bi2Se3 flakes were typically 5 μm wide and 10 μm long, and thicker than 150 nm. Electrodes were patterned by standard electron beam lithography procedures. Perfect Ohmic contacts with 50 nm Cr/100 nm Au were made by electron beam evaporation.

Measurements

In the measurements, a 1064 nm laser light was illuminated on the samples through a chopper. The induced photocurrent was pre-amplified by a voltage pre-amplifier before being measured by a lock-in amplifier. The circular polarization degree p of the incident light was shifted by rotating a quarter-wavelength plate, varying with the phase difference in terms of the rotation angle φ with a period of π. The polarization modulation is in a series from linearly to left-circular, and linearly again, to right circular, and then back to linearly polarized. The polarization-dependent photocurrents (Figure 1d and 1e) were measured at different polarization degree P ∝ sin2φ, where φ is between the polarization direction of the incident light and optical axis of the quarter-wave plate, with a step of π/18.

Author Contributions

J.X.D., N.T. and B.S. conceived and designed the experiments. J.X.D. and X.H. performed sample fabrication. Y.Y. grew the CVD Bi2Se3. J.X.D. conducted the measurements. S.Z. and X.D.Q. assisted with measurements and data analysis. J.X.D., N.T. and W.K.G. wrote the paper. X.Q.W., X.L.Y., F.J.X., Y.H.C. and B.S. gave scientific advice. All the authors contributed through scientific discussion.