# Anomalous galvanomagnetism, cyclotron resonance and microwave spectroscopy of topological insulators

###### Abstract

The surface quantum Hall state, magneto-electric phenomena and their connection to axion electrodynamics have been studied intensively for topological insulators. One of the obstacles for observing such effects comes from nonzero conductivity of the bulk. To overcome this obstacle we propose to use an external magnetic field to suppress the conductivity of the bulk carriers. The magnetic field dependence of galvanomagnetic and electromagnetic responses of the whole system shows anomalies due to broken time-reversal symmetry of the surface quantum Hall state, which can be used for its detection. In particular, we find negative linear dc magnetoresistivity and a quadratic field dependence of the Hall angle, shifted rf cyclotron resonance, nonanalytic microwave transmission coefficient and saturation of the Faraday rotation angle with increasing magnetic field or wave frequency.

## I Introduction

Unlike ordinary band insulators, semiconductors or semimetals, a recently identified class of materials - topological insulators (TIs) Kane05 (); Bernevig06 (); Koenig07 (); Fu07 (); Moore07 (); Hasan10 () - exhibit unusual conducting states on sample boundaries. On the surface of a three-dimensional TI such a state is characterized by a nodal spectrum with a single Dirac cone (or, in general, with odd number of Dirac cones). If time-reversal symmetry (TRS) is broken, an energy gap is induced at the Dirac points, and the surface state exhibits the anomalous quantum Hall (QH) effect. Qi08 (); Essin09 (); Tse10a (); Tse10b () Although the generation of a sizable Dirac gap requires an effort, the surface quantum Hall state in TIs is of great interest because it gives rise to rich magneto-electric phenomena Qi08 (); Essin09 (); Tse10a (); Tse10b (); Maciejko10 (); Garate10 () specific to axion electrodynamics. Wilczek87 ()

There is however a serious obstacle for identifying the surface-related magneto-electric phenomena in three-dimensional TIs. It stems from dissipative bulk conductivity which generally cannot be ignored because of the complex band structure of three-dimensional TI where the Fermi level does not necessarily lie in the bulk band gap or crosses both the surface and bulk states. Chen09 (); Kulbachinskii99 (); Ayala10 (); Ren10 () For a TI film with thickness , bulk zero-field dc conductivity and surface QH conductivity , the contribution of the surface with respect to the bulk is characterized by parameter . Tse10b () It has been shown that the well-resolved surface magneto-electric effects, such as the Kerr or Faraday rotation, require sufficiently large values of . Qi08 (); Maciejko10 (); Tse10b () If, however, the bulk conductivity is much larger than the surface one, i.e. , is it still possible to resolve surface magneto-electric effects in TIs? In this paper we demonstrate such a possibility on several different examples of electrodynamic phenomena.

We show that the surface contribution to the electrodynamics of TIs becomes more pronounced when the bulk conductivity is suppressed by an external magnetic field and by finite frequency of an applied ac electromagnetic field. This can be seen from Boltzmann transport theory expressions for the longitudinal and transverse (Hall) bulk conductivities (see e.g. Refs. LAK56, ; Fawcett64, ):

(1) |

where and the cyclotron frequency are both assumed much smaller than the frequency associated with the surface Dirac gap :

(2) |

is the effective cyclotron mass and is the elastic scattering time. Clearly, with increasing and , the real parts of conductivities can be made comparable with , even though for the zero-field dc case . Under these conditions the TRS breaking on the TI surface leads to anomalous galvanomagnetic and electromagnetic responses of the whole system. In particular, we find (i) negative linear dc magnetoresistivity and Hall angle quadratic with , (ii) rf cyclotron resonance at shifted frequency

(3) |

(iii) nonanalytic -dependence of the microwave transmission coefficient and (iv) saturation of the Faraday rotation angle with increasing magnetic field or wave frequency. Below we explain in detail how these anomalies are related to the surface QH state and how they can be used for its experimental identification.

The paper is organized as follows. In Sec. II we formulate the main equations of electrodynamics of a TI film and discuss approximations used throughout the paper. Then we present the solutions of this electrodynamic problem for different physical situations: galvanomagnetic phenomena (Sec. III), cyclotron resonance (Sec. IV), and electromagnetic transmission and Faraday rotation effects (Sec. V). Finally, in Sec. VI we summarize.

## Ii Formulation of the problem

TRS breaking on the surface of a TI can be achieved by coating it with thin layers of a ferromagnetic (FM) material, Qi08 (); He10 (); Zhu10 () magnetized perpendicularly to the TI film plane by an external dc magnetic field (see, Fig. 1a). The FM magnetization acts on the electron spin, generating an energy gap at the Dirac point, which can be described by a Hamiltonian , Tse10a () where are the spin Pauli matrices, and and are the velocity and momentum for top (t) and bottom (b) surfaces. It is assumed that the Fermi level lies within the gap such that both surfaces of the TI have vanishing dissipative longitudinal conductivities and nonzero quantized Hall conductivities: Streda82 () , with half-integer filling factors . Qi08 (); Tse10b () If the variation of is restricted by Eq. (2), the surface states remain on the QH plateaus, and do not depend on . Strong ()

We also assume that the surface states respond to a time-dependent electromagnetic (EM) field () adiabatically. This is justified for low frequencies which can, at the same time, be much smaller than the plasma frequency (see below). In particular, the surface states remain dissipationless, Qi08 (); Tse10a () i.e. the surface current density induced by the electric field can be written as

(4) |

where is the unit vector perpendicular to the film. The use of the delta functions in Eq. (4) is justified if the penetration lengths of the surface states into the bulk is much smaller than the film thickness . Also, in this case there is no magnetically induced energy gap in the interior of the film. Therefore, the bulk conductivity tensor contains both the dissipative (longitudinal) and Hall components which are both -dependent and given for a single (lightest) carrier group by the Boltzmann transport theory expressions (1). The resulting bulk current density is

(5) |

To find the EM field inside the TI, we use the thin-film approximation , which for the upper frequency limit meV implies m. In addition, the film thickness, , should be smaller than the skin penetration depth, :

(6) |

where is the plasma frequency of the bulk carriers with density , Fermi velocity and momentum . For cm and ms the plasma frequency is s, yielding the lower bound for the skin depth m [see Eq. (6)]. Note also that the plasma frequency s is of the order of the frequency s related to the surface gap meV. Therefore, in addition to requirement we have .

Under condition (6), the electric field inside the film can be approximated by the average value . The equation for is obtained by averaging the Maxwell equation over the film thickness ( is the dielectric constant):

(7) |

The external EM perturbation enters via magnetic fields at the outer top and bottom surfaces of the TI. We specify in each concrete situation considered below.

## Iii Galvanomagnetic phenomena

We begin by considering dc galvanomagnetic phenomena in the standard four-contact geometry (see, e.g. Ref. Kulbachinskii99, and Fig. 1a) in the presence of perpendicular magnetic field and electric current . The current induces the jump of the magnetic field across the film, , so that Eq. (7) determines the longitudinal and Hall electric fields in terms of given : and where the longitudinal and Hall resistivities are given by

(8) | |||

(9) |

For we recover the usual results of the magnetotransport theory: and . The details of the method used in our calculations are very well described in literature (see, e.g. Refs. LAK56, ; Fawcett64, ). The -field independent resistivity reflects strong cyclotron drift in the direction of the current (), induced by the crossed magnetic and electric Hall fields. Moreover, this conclusion remains valid in the nonlinear electrodynamics where the dependence of the bulk conductivity on the magnetic field of the current (or of an external EM wave) is taken into account. Makarov98 ()

However, according to Eq. (9) the nonzero surface conductivity generates an additional Hall field that affects the cyclotron drift in the direction of the current. For this reason (8) acquires the magnetic field dependence. Moreover, since on the QH plateau does not change with , resistivity exhibits linear behavior, whereas the Hall angle (9) has the anomalous quadratic term [Figs. 1b and c]. We note that the anomalous terms in and remain identifiable even if several bulk carrier groups are taken into account, because in that case the bulk resistivities are still regular analytic functions of which can be subtracted from the total and .

## Iv Cyclotron resonance

Let us consider now the bulk cyclotron resonance. It can be realized in a contactless setup where the sample is placed in the maximum of the electric field of an rf resonator normal mode, which generates an antisymmetric magnetic field across the sample, i.e. . At the rf frequencies the displacement current contribution is usually smaller than conductivity , so that induces an average ac current density in a contactless way. The relevant observable is the longitudinal ac resistivity which we find from Eq. (7) as

(10) |

For it has a resonance in both and dependencies, shown in Fig. 2. Note that the linear on-resonance relation between and is the hallmark of the cyclotron resonance. However, the resonant frequency (3) is shifted with respect to because of the additional drift in the direction of the current, induced by the surface contribution to the Hall electric field in Eq. (9). For bulk Drude conductivity the frequency shift depends on the bulk carrier density , Fermi velocity and momentum as

(11) |

For cm, nm, ms and the resonance frequency shift is s, i.e. well below both s and s.

## V Electromagnetic transmission and Faraday rotation

We now turn to microwave spectroscopy which also allows one to probe the surface states in TIs Ayala10 (). We consider an EM wave incident normally at the bottom surface of a TI and will analyze both the transmission coefficient and the Faraday rotation of the EM field plane in the transmitted wave [see, also Fig. 3a]. This situation involves a new conductivity scale, viz. the inverse impedance of the dielectric media surrounding the TI, , where are the dielectric constants of the top and bottom materials (shown in gray in Fig. 3a). The new conductivity scale is important because it is much larger than the surface QH conductivity : the product is a small parameter proportional to the fine structure constant Qi08 (); Tse10a (); Tse10b ():

(12) |

Therefore, the response of the surface state to an EM wave is generically rather weak. To proceed we note that the electric and magnetic fields on the outer surfaces of the TI film are

(13) | |||

(14) |

where refers to the transmitted wave on the top surface, whereas and label the incident and reflected waves on the bottom () surface. In the thin film approximation the electric field on each surface equals to the average field. Eliminating the reflected field through , we express the magnetic field difference as insert this into Eq. (7), and solve it for . Since the incident electric field can be regarded real, we present the solution for the real part of the transmitted wave:

(15) | |||

(16) |

where is the transmission coefficient, is the rotation angle of the EM field plane with respect to the incident wave (Faraday angle), and and are real functions given by

(17) |

(18) |

Similar to resistivity (8), the TRS breaking leads to a nonanalytic linear dependence of the transmission coefficient (16). To illustrate this we extract the large zero-field value from and plot in Figs. 3b and c the difference for (solid curves) and for (dashed curves). The magnetic field range in which can be tuned by varying parameters and . This should help in finding the optimal regime for observation of the predicted anomalous magnetic-field dependence of .

Figure 3d shows the low-frequency Faraday angle (16) in units of the fine structure constant as a function of the magnetic field for zero and finite bulk conductivity . For the Faraday angle contains only the surface contribution Qi08 (); Tse10a (); Tse10b (). For the bulk contribution makes the dependence nonmonotonic with the following asymptotics:

(19) |

The low-field limit () agrees with the result of Ref. Tse10b, which found smaller than the surface contribution for nonzero bulk conductivity , e.g. for . However, for strong fields we find saturation of precisely at the surface value because of the suppression of the bulk conductivity via classical cyclotron motion. This way of extracting the topological surface contribution may have an advantage over the previously proposed low-field detection scheme Qi08 () because the magnetization of the FMs in strong fields leads to a more robust surface Dirac gap . As seen from Fig. 3e, the frequency dependence of the Faraday angle also saturates at the surface value , which can be used for its detection as well.

## Vi Summary

In summary, we have investigated galvanomagnetic and electromagnetic properties of topological insulators in which time-reversal symmetry is broken due to the surface quantum Hall effect. Our model includes both the dissipationless quantum Hall conductivity on the surface and the classical magnetoconductivity in the bulk of the system. Although the zero-field dc bulk conductivity may significantly exceed the surface one, the surface contribution can still be detected through anomalous magnetic field dependencies of electrodynamic responses, revealing the underlying broken time-reversal symmetry. With appropriate modifications our findings can be extended to HgTe quantum wells which also support single-valley Dirac fermions Buettner11 (); GT11a (); GT11b (); GT11c (); GT10 () and show a pronounced Faraday effect. Shuvaev11 ()

###### Acknowledgements.

We thank A. H. MacDonald, S.-C. Zhang, L.W. Molenkamp, A. Pimenov, A. M. Shuvaev and G. V. Astakhov for helpful discussions. The work was supported by DFG grant HA5893/1-1.## References

- (1) C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 226801 (2005).
- (2) B. A. Bernevig and T. L. Hughes and S. C. Zhang, Science 314, 1757 (2006).
- (3) M. König, S. Wiedmann, C. Brüne, A. Roth, H. Buhmann, L. W. Molenkamp, X.-L. Qi, and S.-C. Zhang, Science 318, 766 (2007).
- (4) L. Fu, C. L. Kane and E. J. Mele, Phys. Rev. Lett. 98, 106803 (2007).
- (5) J. E. Moore and L. Balents, Phys. Rev. B 75, 121306(R) (2007).
- (6) M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 (2010); X.-L. Qi and S.-C. Zhang, e-print arXiv:1008.2026 (unpublished).
- (7) X.-L. Qi, T. L. Hughes, and S.-C. Zhang, Phys. Rev. B 78, 195424 (2008).
- (8) A. M. Essin, J. E. Moore, and D. Vanderbilt, Phys. Rev. Lett. 102, 146805 (2009).
- (9) W.-K. Tse and A. H. MacDonald, Phys. Rev. Lett. 105, 057401 (2010).
- (10) W.-K. Tse and A. H. MacDonald, Phys. Rev. B 82, 161104(R) (2010).
- (11) J. Maciejko, X.-L. Qi, H. D. Drew, and S.-C. Zhang, Phys. Rev. Lett. 105, 166803 (2010).
- (12) I. Garate and M. Franz, Phys. Rev. Lett. 104, 146802 (2010).
- (13) F. Wilczek, Phys. Rev. Lett. 58, 1799 (1987).
- (14) Y. L. Chen, 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).
- (15) V.A. Kulbachinskii, N. Miura, H. Nakagawa, H. Arimoto, T. Ikaida, P. Lostak, and C. Drasar, Phys. Rev. B 59, 15733 (1999).
- (16) O. E. Ayala-Valenzuela, J. G. Analytis, J.-H. Chu, M. M. Altarawneh, I. R. Fisher, and R. D. McDonald, e-print arXiv:1004.2311 (unpublished).
- (17) Z. Ren, A. A. Taskin, S. Sasaki, K. Segawa, and Y. Ando, Phys. Rev. B 82, 241306(R) (2010).
- (18) I. M. Lifshitz, M. Ya. Azbel, and M. I. Kaganov, Zh. Eksp. Theor. Fiz. 31, 63 (1956) [Sov. Phys. JETP 4, 41 (1957)].
- (19) E. Fawcett, Adv. Phys. 13, 139 (1964).
- (20) H.-T. He, G. Wang, T. Zhang, I.-K. Sou, G. K. L. Wong, J.-N. Wang, H.-Z. Lu, S.-Q. Shen, and F.-C. Zhang, Phys. Rev. Lett. 106, 166805 (2011).
- (21) J.-J. Zhu, D.-X. Yao, S.-C. Zhang, and K. Chang, Phys. Rev. Lett. 106, 097201 (2011).
- (22) P. Streda, J. Phys. C 15, L717 (1982).
- (23) In contrast, Ref. Tse10b, considers strong orbital and Zeeman magnetic field effects on the surface and no magnetic field influence in the bulk.
- (24) N. M. Makarov, G. B. Tkachev, and V. E. Vekslerchik, J. Phys.: Condens. Matter. 10, 1033 (1998).
- (25) B. Büttner, C. X. Liu, G. Tkachov, E. G. Novik, C. Brüne, H. Buhmann, E. M. Hankiewicz, P. Recher, B. Trauzettel, S. C. Zhang and L. W. Molenkamp, Nature Phys. 7, 418 (2011).
- (26) G. Tkachov, C. Thienel, V. Pinneker, B. Büttner, C. Brüne, H. Buhmann, L. W. Molenkamp, and E. M. Hankiewicz, Phys. Rev. Lett. 106, 076802 (2011).
- (27) G. Tkachov and E. M. Hankiewicz, Phys. Rev. B 84, 035444 (2011).
- (28) G. Tkachov and E. M. Hankiewicz, Phys. Rev. B 83, 155412 (2011).
- (29) G. Tkachov and E. M. Hankiewicz, Phys. Rev. Lett. 104, 166803 (2010).
- (30) A. M. Shuvaev, G. V. Astakhov, A. Pimenov, C. Brüne, H. Buhmann, and L. W. Molenkamp, Phys. Rev. Lett. 106, 107404 (2011).