# Spin Response of Electrons on the Surface of a Topological Insulator

###### Abstract

The surface of a topological insulator hosts a very special form of a quasi-two dimensional metallic system when it is embedded in a topologically trivial medium like the vacuum. The electronic properties of this unusual 2D metal are distinct in many aspects from both the conventional two-dimensional electron gas systems in quantum well heterostructures as well as those of a single layer graphene. In this paper, we study one of these distinct features i.e., the response of the electronic spins to an applied magnetic field perpendicular to the surface. We find an unusual behaviour of the spin magnetization and susceptibility as a function of both the magnetic field and the chemical potential for a generic topological surface. We propose that this behavior could be studied by the recently developed experimental technique called -NMR which is highly sensitive to the surface electron spins. We explain how this technique could be used to probe for spontaneous magnetic ordering caused by magnetic dopants or interactions discussed in the recent literature.

## I Introduction

Topologically protected electronic states residing on the surface of topological insulators (TI) hasan1 (); hasan2 (); shen1 (); hasan3 (); hasan4 (); fu-kane3D (); hasan_rev (); Zhanybek () form a unique 2D metal distinct from those so far realized in solid state systems. Similar to the low energy electronic states in a single layer grapheneNeto (), the robust 2D metal on an TI surface has conic branches touching (in the absence of the intrinsic gap) at high-symmetry points in the first Brillouin zone. This resemblance in the energy dispersion is responsible for some common properties between these metallic systems, such as the square-root dependence of Landau Levels (LLs) to the applied magnetic field. However, there are important differences between these systems as a result of the topological nature of the surface states in TIs.

The low energy Dirac-like surface bands of a topological insulator arise due to the mismatch of the bulk topological invariants on the two sides of the surface regardless of how constituent atoms have been arranged on the surface as long as the bulk remains in the topological phasefu-kane3D (); zhangprb (). In the case of a ‘strong’ topological insulator (STI) there is an odd number of Dirac points at any surface and, importantly, they are not spin degenerate. Furthermore, they exhibit a unique spin-momentum lockinghasan3 (); hasan4 () as a result of the strong spin-orbit interaction in the underlying STI which is essential for the formation of the topologically non-trivial insulator phase with time-reversal symmetry. These properties make electronic surface states in a STI very distinct from those of graphene with a pair of spin degenerate Dirac-like bands near two special points in the Brillouin zone which exist as a result of the unique arrangement of carbon atoms in the honeycomb lattice.

One of the interesting aspects in which these 2D metallic systems behave uniquely is their magnetic response.DiracMag2 () Here, motivated by the recent progres in the experimental methods and the importance of the spin susceptibility for our understanding of electronic systems, we study the spin response of electrons on the surface of a STI and show that it exhibits interesting features even in a simple non-interacting limit. We find that the characteristic spin-momentum locking of these electrons leads to a unique spin response that is distinct from spin-degenerate systems like graphene and 2DEG. Instead of the oscillatory behaviour of the susceptibility as a function of the chemical potential found in spin degenerate systems, in STIs our work predicts a plateau-type behaviour which arises from the strong correlation between spin and orbital degrees of freedom. In addition we find that the existence of a special LL with full spin polarization leads to a jump in the magnetization as the chemical potential crosses the energy of this LL. When an intrinsic magnetic ordering is present there is also a jump in the magnetization as a function of the applied magnetic field. Our results can therefore assist in detection of such intrinsic magnetizationMagDopped (); MagDopped2 () which is known to have profound consequences for the nature of the surface state; the magnetized surface state of a TI is predicted to become a quantum Hall liquid with half-integer quantized Hall conductivityzhangprb (); fu2 () and many unusualdhlee1 (); qi2 (); wormhole (); tse1 (); maciejko1 (); garate1 () and potentially usefulnagaosa1 (); garate2 () physical properties.

There are various ways of measuring the weak magnetization produced by electrons on a metallic surface. The SQUID scanning magnetometry is a highly sensitive probe which can detect tiny magnetization on the surfaces. However, it is a challenge to use this device to probe magnetization when there is a large applied magnetic field which interferes with the superconducting part of the SQUID and causes noise. A variant of this method using a superconducting pickup coil has been developed in order to study deHaas-van Alphen oscillations of the 2D metallic systems in a large perpendicular magnetic field.Squid1 (); Squid2 () High-sensitity micro-mechanical cantilever magnetometry MagExp1 () is another way of measuring electronic magnetization, however, this method measures the magnetization of the whole sample and it is difficult to isolate the contribution from the electrons on a single surface.

The nuclear magnetic resonance (NMR) is another powerfull experimental technique which can be used to study electronic spin magnetic response in a bulk metal.nmrbook () Through the so-called Knight shift in the nuclear resonance peak, it is possible to probe the spin part of the magnetic susceptibility of the electrons as they interact with the resonating nuclei in their proximity. Unfortunately, this method, in its conventional form, fails to be useful in very thin films and 2D metallic systems due to the limitation in the number of available nuclei and the resulting weakness of the signal. Thanks to the progresses made by experimentalists in controlling and implementing high energy beams of unstable ions, an experimental technique, the so-called -NMR, has been recently developed to overcome the above limitations. Briefly, in this exotic variety of NMR, unstable radioactive ions such as Li and Be are implanted in the sample surface. The nuclear spin precession signal is then detected through the products of the beta decay of the radioactive nucleus. Since the ion implantation depth can be controlled by tuning the beam energy it is possible to acquire information about the behaviour of the electronic spins in very thin metallic films.betaNMR1 (); betaNMR2 () Experiments are currently underway to study the surface magnetic response of STI crystals BiSe and BiTe using -NMR.farlane ()

The paper is organized as follows. First we introduce a simple model known to describe the electrons on the surface of a STI in the presence of a perpendicular magnetic field. We review the exact solutions of its eigenvalue problem which has been studied previously in various contexts.DiracMag1 (); DiracMag2 () In section III, we calculate the spin magnetization and susceptibility assuming that the chemical potential lies inside the gap between positive and negative energy eigenstates. We discuss the magnetic response of the surface as one tunes the magnetic field and the chemical potential for various values of the intrinsic gap. In closing, we explain how a -NMR experiment might be able to detect these effects through Knight shift measurements.

## Ii The Model

Many of the interesting features of the surface states in a topological insulator can be captured, at least qualitatively, using a simple non-interacting Dirac Hamiltonian

(1) |

where and is the fermionic creation operator of the spin up(down) states with wave vector . is the intrinsic gap in the surface spectrum which might be nonzero when the time-reversal symmetry is spontaneously broken due to magnetic ordering. The latter can arise due to the presence of magnetic dopants with spin in the proximity of the surface exchange-coupled to the electronic spinsMagDopped (); MagDopped2 (); Rosenberg () or as a result of electron-electron interactionstern1 ().

Turning on a perpendicular magnetic field adds two terms to the above Hamiltonian. One is the minimal coupling of the magnetic vector potential, (electron charge ). The other is the coupling of the spins to the magnetic field, the Zeeman effect, expressed as , where is the effective electron gyromagnetic constant. In the bulk BiSe crystal (Ref. ZhangFP, ) although much smaller values have been reported for electrons near the surface.ong1 () In our calculations below we use two representative values, and , which yield qualitatively similar results with some interesting differences.

In the continuum limit, the leading order Hamiltonian describing these surface states in the presence of the applied magnetic field, , in the Landau gauge , can be written in the following form

(2) |

where is the creation operator for the spinor mode extended along the direction and localized at . for is defined as

(3) |

where and

(4) |

The term in proportional to the magnetic field is the Zeeman contribution. is the one-dimensional harmonic oscillator bosonic operator defined as

(5) |

Note that varying shifts the position of the localized state produced by the application of on the vacuum state along the direction. In the case, can be obtained from the one given for in Eq. (3) by exchanging the off-diagonal elements and replacing by its time reversed counterpart .

The eigenstates for and are given by DiracMag1 (); DiracMag2 ()

(6) |

(7) |

while for we have

(8) |

In the above and are the one-dimensional harmonic oscillator eigenstates, i.e., . We remark that the eigenstate in Eq. (8) is very special since it is fully spin polarized. This will have important consequences for the magnetic response discussed below.

## Iii Spin Susceptibility and magnetization

The electronic magnetic moment due to spin is proportional to the spin operator

(11) |

Therefore, to calculate the spin part of the magnetic moment for the eigenstates given in the previous section we only need to find the expectation value of the spin operator. A straightforward evaluation using Eqs. (6-8) shows that all the electronic states within the same LL contribute equally to the magnetization. For each of them we have

(12) |

(13) |

(14) |

The total magnetization for each Landau level can be obtained by multiplying the above quantities by the Landau level degeneracy , representing the total number of states with characteristic length that the surface area can accommodate. Since the magnetization contribution computed above for each eigenstate is an explicit function of its energy, we can perform the following integral to find the total magnetization density due to the electronic spins

(15) |

where is the magnetization of the eigenstate with energy and is the Fermi-Dirac distribution function. The electronic density of states associated with the surface states is . For the Hamiltonian we used in the previous section it takes the following form

(16) |

where is the Landau level index beyond which the energy exceeds the cutoff energy . The cutoff can be chosen to be the energy where the surface band becomes degenerate with the bulk bands and here we assume meV. Using this density of states function to perform integration in Eq. (15) yields

(17) |

where .

For a constant magnetic field the magnetization of a single TI surface given by Eq. (17) shows an interesting behaviour as a function of chemical potential illustrated in Fig. 2. In order to avoid ambiguity associated with the high-energy cutoff we choose to display , i.e. spin magnetization relative to the neutrality point at . For the negative values of magnetization initially decreases reflecting the fact that the negative-energy surface states exhibit negative spin polarization, as can be seen from Eq. (13). The large jump in near results from electrons filling the fully spin-polarized Landau level. Further increase in results in increase of now reflecting the fact that the positive-energy surface states exhibit positive spin polarization.

The above behavior is unique to topological insulators for it results from the Landau level structure in a single Dirac point (or more generally an odd number thereof). In a TI the ‘other’ Dirac point is located on the opposite surface where the magnetic field points in the direction opposite relative to the surface normal, see Fig. 1. The contribution of this surface to would be the same. We emphasize that for a relatively thick TI slab -NMR will be sensitive to a single surface facing the beam and the behaviour predicted here is in principle observable, except that continuous tuning of the chemical potential in a crystal might be difficult to achieve.

A much more feasible experiment involves varying magnetic field while keeping constant. We now show that a unique signature of the fully polarized Landau level still exists in samples with intrinsic magnetic ordering, when the chemical potential resides in the gap (i.e. both the bulk and the surface are insulating in the absence of the field). In his situation we can set in our model. Now consider the effect of the applied magnetic field. Since the gap between the adjacent energy levels in which the Fermi energy is located is fairly large ( K for a 1T field), at sufficiently low temperatures we can replace the Fermi function in Eq. (17) by the step function and write

(18) |

We have assumed that the chemical potential remains pinned at zero energy and does not change as we tune the magnetic field. With these assumptions the magnetization has an interesting discontinuity at a finite magnetic field. The discontinuity in occurs since tuning the magnetic field forces the fully spin polarized Landau level energy to evolve according to Eq. (4) and cross the chemical potential. The critical magnetic field at which the jump happens is given by

(19) |

Assuming that the intrinsic gap remains independent of we plot the resulting magnetization in Fig. 4. This predicted behaviour could be employed to experimentally detect the intrinsic magnetization gap in the surface of a magnetically doped TI and measure the size of through Eq. (19). Although this signature only occurs when the chemical potential lies very close to the surface state Dirac point we note that magnetically doped samples satisfying this requirement have been grown and studied.MagDopped ()

It is important to note that is in general not independent of since magnetic moments of dopants will tend to align with the applied field. Thus, like in a ferromagnet, there will be a hysteresis effect whose features will depend on the material detailsMagDopped (). One possible scenario is shown schematically in Fig. 5. In the case when depends on the field the equation (19) continues to hold but must now be viewed as an implicit equation for the critical field .

Another question that arises has to do with the origin of the electrons that fill the Landau level upon changing the field through . One may wonder where the extra electrons come from in a fully gapped isolated system. The answer lies in the side surfaces which under generic conditions remain gapless and act as a reservoir of electrons. The model we consider here does not capture these states but they reflect themselves in solutions of the Hamiltonian which are not normalizable in an infinite system. These are fully spin polarized with opposite energy and spin direction. In fact they are the particle-hole conjugates of the fully spin polarized LL given in Eq. (8). Taking into account these electronic states and the fact that it is more favourable for electrons with higher energy to be transferred to the negative energy states the counting problem can be resolved.

## Iv The Knight Shift

We now outline how a -NMR experiment can in principle be used to probe some of the physics discussed in the previous section. The valence and conduction electrons in a metal posses magnetic moments arising from both their orbital motion and their spin degrees of freedom. Nuclear magnetic resonance technique can be used as a probe of the spin part of the total magnetization in the presence of the magnetic field by measuring the relative shift in the nuclear resonance peak with respect to the same resonance peak in a reference insulating system. This effect, which is due to the interaction between electronic spins and those of the nuclei, is known in the literature as the Knight shift and has been extensively studied in both theory and experiment. nmrbook ()

The mobile electrons in a metal interact with the nuclei in their proximity and the Knight shift in the resonance peak of these nuclei can be described by a local Fermi contact interaction term given by where is the position of the th nucleus and is its magnetic moment. The magnitude of , the gyromagnetic ratio, depends on the nucleus quantum state. The nucleus total spin, , couples to the applied magnetic field and therefore the position of the peak depends on the magnitude of the total magnetic field experienced by the nucleus which has a contribution due to the interaction with electrons. It turns out that this shift is proportional to the spin susceptibility. The constant of proportionality, known as the hyperfine coupling, can be computed using first principle calculations for the implanted nuclei. On the other hand, if we assume that the presence of the nuclei does not significantly alter the electronic states, then it is possible to approximate the shift for them by taking the expectation value of the aforementioned interaction term using the unperturbed electronic states. This is the lowest order approximation in the perturbative treatment of the interaction term. For metallic systems with spin degenerate bands this shift is proportional to the spin susceptibility as it can been seen from a simple calculation considering the fact that the spatial and spin degrees of freedom are uncorrelatednmrbook ().

The spin-momentum locking on the surface of a TI along with the energy dependence of the penetration depth can in principle change the above simple physics. Since it is not possible anymore to separate spin and orbital degrees of freedom, one might question the validity of the linear relation between the Knight shift and the spin susceptibility. We devote the rest of this section to addressing this issue by considering a very simple model. We assume that the nuclei do not alter the electronic states around them. It is important to note that this assumption may break down for the implanted nuclei if they modify the electronic states around them significantly and computing the Knight shift would then require a first principle calculations.

The field experienced by the th nucleus due to the interaction with the proximate electrons is given by

(20) |

where is the expectation value of over electronic states at temperature . Therefore, the effective Hamiltonian for the ensemble of nuclei takes the form

(21) |

This way, th nucleus would have a resonance peak . The Knight shift is then defined by comparing the resonance frequency with the frequency in a similar material without these electronic states

(22) |

The shift in the resonance peak of the nuclei ensemble is the average of the knight shift from each individual nucleus and is given by

(23) |

where is the number of the implanted nuclei. Using the electronic eigenstates given in Eq. (6-8) we get the following expression for

(24) |

Here appears as a factor in the realistic 3D electronic wave functions of the surface electrons reflecting the fact that the electrons have an energy dependent penetration depth into the bulk. The nuclei implanted in the system have a spatial probablity distribution , which depends on the energy and diameter of the beam of the ions used in the -NMR experiment. If the distribution function is known, we can replace the above summation with a 3D integral over the crystal volume

(25) |

Assuming that the distribution is uniform in the plane of the surface, i.e., , we get

(26) |

where we have performed the integration over the in-plane degrees of freedom and replaced the summation over with the LL degeneracy. We have also defined the th LL weight, , as

(27) |

Now if we assume that different LL have the same penetration depth we have for all and we obtain

(28) |

We thus recover the linear proportionality of the Knight shift to the surface electronic spin susceptibility under reasonable assumptions. Note that if we relax the assumption that different LLs can now have different penetration depths, then the Knight shift would no longer be linearly proportional
to the total spin susceptibility. Instead, it would be a weighted superposition of contributions from each individual LL to the spin susceptibility. Nevertheless, the Knight shift will still display the interesting behavior discussed in this study as long as is a reasonably slowly varying function of .
We should emphasize once again that although above considerations
elaborate on the differences caused by the spin-momentum locking
and the energy dependent penetration depth, they do not take into account the fact that the electronic wave-functions could be altered by the
presence of the implanted nuclei and the hybridization with the adjacent nuclei.

## V Conclusions

The magnetic response of spins of the Dirac-like electrons on the surface of a topological insulator shows interesting features both in the absence and the presence of an intrinsic gap . When the surface states are gapped owing to the time-reversal breaking perturbation (i.e. due to magnetic doping) the Landau level, which is fully spin polarized, can have positive or negative energy depending on the sign of the intrinsic gap relative to that of the applied magnetic field. It will therefore be completely filled or empty when the chemical potential is tuned to zero energy. Our study shows that this structure results in an observable jump in the spin susceptibility, mesurable e.g. through the -NMR Knight shift, as one tunes the applied magnetic field through the critical value given in Eq. (19). The effect may be used as a means to measure the magnitude of the intrinsic gap on the surface of a magnetically doped topological insulator if -NMR or another surface-sensitive technique could capture the magnetic response of the electronic spins on the surface. This behaviour is a unique feature of the topological insulator exotic surface states closely related to the special form of the spin-momentum entanglement.

## Vi Acknowledgment

The authors are indebted to W.A. MacFarlane, R.F. Kiefl and K.A. Moler for illuminating discussions. This work was supported in part by NSERC and CIfAR.

## References

- (1) M.Z. Hasan, C.L. Kane, Rev. Mod. Phys. 82 3045 (2010)
- (2) D. Hsieh, D. Qian, L. Wray, Y. Xia, Y. S. Hor, R. J. Cava, and M. Z. Hasan, Nature 452, 970 (2008).
- (3) Y. Xia, L. Wray, D. Qian, D. Hsieh, A. Pal, H. Lin, A. Bansil, D. Grauer, Y.S. Hor, R.J. Cava, M.Z. Hasan, Nature Phys. 5, 398 (2009).
- (4) Chen, Y.L. J. G. Analytis, J.-H. Chu, Z. K. Liu, S.-K. Mo, X. L. Qi, H. J. Zhang, D. H. Lu, X. Dai, Z. Fang, S. C. Zhang, I. R. Fisher, Z. Hussain and Z.-X. Shen, Science 325, 178 (2009).
- (5) D. Hsieh, D. Qian, L. Wray, Y. Xia, Y. S. Hor, R. J. Cava, and M. Z. Hasan, Nature 460, 1101 (2009).
- (6) D. Hsieh, Y. Xia, L. Wray, D. Qian, A. Pal, J. H. Dil, J. Osterwalder, F. Meier, G. Bihlmayer, C. L. Kane, Y. S. Hor, R. J. Cava and M. Z. Hasan, Science 323, 919 (2009).
- (7) L. Fu, C. L. Kane, and E. J. Mele, Phys. Rev. Lett. 98 106803 (2007).
- (8) Z. Alpichshev, J. G. Analytis, J.-H. Chu, I. R. Fisher, Y. L. Chen, Z. X. Shen, A. Fang, and A. Kapitulnik, Phys. Rev. Lett. 104, 016401 (2010)
- (9) A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov and A. K. Geim, Rev. Mod. Phys. 81, 109 (2009)
- (10) X. L. Qi, T. L. Hughes, and S. C. Zhang, Phys. Rev. B78 195424 (2008)
- (11) Z. Wang, Z-G. Fu, S-Xi. Wang, and P. Zhang, Phys. Rev. B82 085429 (2010)
- (12) Y. L. Chen et al. Science 329, 659 (2010).
- (13) R. Yu et al. Science 329, 61 (2010).
- (14) L. Fu and C. L. Kane, Phys. Rev. B76, 045302 (2007).
- (15) D.-H. Lee, Phys. Rev. Lett. 103, 196804 (2009).
- (16) X.-L. Qi, R. Li, J. Zang, S.-C. Zhang, Science 323, 1184 (2009).
- (17) G. Rosenberg, H.-M. Guo, M. Franz, Phys. Rev. B82, 041104(R) (2010).
- (18) W.-K. Tse and A. H. MacDonald, Phys. Rev. Lett. 105, 057401 (2010).
- (19) J. Maciejko, X.-L. Qi, H.D. Drew, and S.-C. Zhang, Phys. Rev. Lett. 105, 166803 (2010).
- (20) I. Garate, M. Franz, Phys. Rev. B84, 045403 (2011).
- (21) T. Yokoyama, Y. Tanaka, and N. Nagaosa, Phys. Rev. B81, 121401(R) (2010).
- (22) I. Garate, M. Franz, Phys. Rev. Lett. 104, 146802 (2010).
- (23) I. Meinel, T. Hengstmann, D. Grundler, and D. Heitmann, Phys. Rev. Lett. 82 819 (1999)
- (24) I. Meinel, D. Grundler, S. B-Franke, C. Heyn, D. Heitmann, and B. David, Appl. Phys. Lett. 70 3305 (1997)
- (25) J. G. E. Harris, R. Knobel, K. D. Maranowski, A. C. Gossard, N. Samarth, and D. D. Awschalom, Phys. Rev. Lett. 86 4644 (2001)
- (26) Principles of Magnetic Resonance, C.P. Slichter (Springer, Berlin, 1990).
- (27) G. D. Morris, W. A. MacFarlane, K. H. Chow, Z. Salman, D. J. Arseneau, S. Daviel, A. Hatakeyama, S. R. Kreitzman, C. D. P. Levy, R. Poutissou, R. H. Heffner, J. E. Elenewski, L. H. Greene, and R. F. Kiefl Phys. Rev. Lett. 93 157601 (2004)
- (28) T.J. Parolin, J. Shi, Z. Salman, K. H. Chow, P. Dosanjh, H. Saadaoui, Q. Song, M. D. Hossain, R. F. Kiefl, C. D. P. Levy, M. R. Pearson, and W. A. MacFarlane, Phys. Rev. B80 174109 (2009)
- (29) W.A. MacFarlane and R.F. Kiefl (private communication).
- (30) J. W. McClure, Phys. Rev. 104, 666 (1956)
- (31) G. Rosenberg, M. Franz, unpublished.
- (32) Y. Baum, A. Stern, arXiv:1110.4037
- (33) C.-X. Liu, X.-L. Qi, H. Zhang, Xi Dai, Z. Fang, and S.-C. Zhang Phys. Rev. B82 045122 (2010)
- (34) J. Xiong, Y. Luo, Y.-H. Khoo, S. Jia, R.J. Cava, N.P. Ong, arXiv:1111.6031