Giant anisotropic magnetoresistance in a quantum anomalous Hall insulator.

When a three-dimensional ferromagnetic topological insulator thin film is magnetized out-of-plane, conduction ideally occurs through dissipationless, one-dimensional (1D) chiral states that are characterized by a quantized, zero-field Hall conductance. The recent realization of this phenomenon, the quantum anomalous Hall effect, provides a conceptually new platform for studies of 1D transport, distinct from the traditionally studied quantum Hall effects that arise from Landau level formation. An important question arises in this context: how do these 1D edge states evolve as the magnetization is changed from out-of-plane to in-plane? We examine this question by studying the field-tilt-driven crossover from predominantly edge-state transport to diffusive transport in Crx(Bi,Sb)(2-x)Te3 thin films. This crossover manifests itself in a giant, electrically tunable anisotropic magnetoresistance that we explain by employing a Landauer-Büttiker formalism. Our methodology provides a powerful means of quantifying dissipative effects in temperature and chemical potential regimes far from perfect quantization.

In the absence of magnetic doping, surface electrons in three-dimensional (3D) topological insulators such as (Bi,Sb)2Te3 are helical Dirac fermions protected from backscattering by time-reversal symmetry12. When such a 3D topological insulator is made ferromagnetic by doping with magnetic ions3456789 or by integration into magnetic heterostructures101112, time-reversal symmetry is broken. If the magnetization M is perpendicular to the surface, the surface state dispersion is gapped at the Dirac point1314, while M parallel to the surface is expected to keep the Dirac point intact with the origin shifting in momentum-space15. In a ferromagnetic topological insulator thin film with out-of-plane magnetic anisotropy16, the two surfaces are gapped, albeit with opposite mass. The mass domain wall created along the side walls of films then hosts a massless, one-dimensional (1D), dissipation-free edge state. If the chemical potential is tuned inside the magnetic gap of the surface states, transport is solely restricted to these 1D edge states, and characterized by a quantized Hall resistivity ρ=h/e2 (where h and e are the Planck constant and elementary charge, respectively) and (ideally) vanishing longitudinal resistivity ρ=0. This is known as the quantum anomalous Hall (QAH) state and was recently realized experimentally in thin films of Cr(Bi,Sb)2−Te3 (refs 6, 8, 9) at dilution fridge temperatures.

Anisotropic magnetoresistance (AMR)1718 is the response of the longitudinal resistance of a ferromagnet to the angle between the magnetization and direction of current flow. In this work, by tilting the magnetization of QAH insulator in-plane, we observe a novel, giant AMR effect that serves as a powerful quantitative probe of edge transport, even in temperature and chemical potential regimes away from perfect quantization (for instance, at high temperatures compared with the energy scale of the magnetic gap). Tilting the magnetization in-plane closes the magnetic gap and consequently destroys the edge states, leading to a crossover from dissipationless to diffusive transport. We demonstrate that this crossover yields a sharp change in the longitudinal resistivity (ρ) whose angular dependence is unlike classical AMR and show that modelling this dependence can be used to quantify edge-state contributions to transport, and the extent of its inter-mixing with dissipative channels.

Results

Quantum anomalous Hall effect

The magnetic topological insulator devices used in this study were fabricated from epitaxial films of Cr(Bi,Sb)2−Te3 grown by molecular beam epitaxy on 〈111〉 SrTiO3 (STO) substrates and measured in a vector magnet cryostat down to 280 mK (Fig. 1a, see Methods for details). The large dielectric constant of the STO substrate at cryogenic temperatures allows effective electrical back gating for tuning the position of the chemical potential19. Three thin-film devices—A, B and C, were studied with varying proximity to the QAH state (Table 1). The transport properties of these devices are detailed in Supplementary Figs 1–10 and Supplementary Notes 1–3. Unless otherwise mentioned, the data presented in the main text of this manuscript is from device A—a Hall bar of channel dimensions 1 mm × 0.5 mm fashioned from a 10 quintuple layer thin film of Cr0.1(Bi0.5,Sb0.5)1.9Te3. The ferromagnetic Curie temperature of this film is TC∼8 K, obtained from the temperature dependence of the transverse resistivity ρ. As shown in Fig. 1b,c, even at 280 mK, the AHE in device A takes a maximum value ρ∼0.95h/e2 (σ∼0.99h/e2) and the longitudinal resistivity its minimum value ρ∼0.25 h/e2, both at zero field, when the gate voltage is tuned to its value at charge neutrality VG0. This non-vanishing ρ, accompanied by ρ taking near-quantized values, is indicative of dissipative channels coexisting with the 1D chiral edge modes. ρ also displays a prominent hysteretic peak at magnetization reversal (Fig. 1b,c), and takes very large values approaching 260%. This large effect has been interpreted as enhanced backscattering between edge channels in a manner similar to the sharp ρ rise in the plateau transition of the quantum Hall effect. Below 3 K, at charge neutrality, ρ displays a metallic behaviour as ρ approaches quantization (Fig. 1d,e). This metallic behaviour is another signature of edge-state transport and ρ recovers its insulating behaviour when the chemical potential is tuned far out of the magnetic gap (Supplementary Fig. 2). Similarly, in the gate voltage dependence, as we tune the chemical potential towards the magnetic gap, ρ drops sharply, as ρ rises and finally flattens close to its maximum value.

Figure 1

Quantum anomalous Hall effect

(a) Cartoon of device and measurement scheme. Current flows in the x direction, and z is the direction perpendicular to plane. (b) Near-quantized anomalous Hall effect at charge neutrality (VG=VG0) and T=280 mK, with ρ taking a maximum value of 0.95 h/e2 at zero field. (b,c) Corresponding longitudinal magnetoresistance, displaying a minimum ρ of 0.25 h/e2, also at zero field. (d,e) Characteristic signatures of quantum anomalous Hall effect in the temperature and gate voltage dependence, respectively, of ρ (red) and ρ (blue), after magnetic training, at zero field and 200 mT, respectively. As edge-state transport dominates over diffusive transport, the decrease in ρ is complemented by a rise in ρ.

Table 1

Comparison of transport properties and AMR fitting parameters of three different Cr(Bi,Sb)2−Te3 thin-film devices with varying degrees of dissipation.

Device	ρxy (h/e2)	ρxx (0)(h/e2)	η		AMR (%)
A	0.95	0.25	0.66	0.62	140.9
B	0.72	0.55	0.16	0.54	75.6
C	0.48	2.79	0.013	0.025	100.2

AMR, Anisotropic magnetoresistance.

The proximity of devices A, B and C to the ideal QAH regime is seen by a comparison of their zero-field ρ and ρ, and is indicative of the strength of dissipation. The effect of dissipative channels is quantified by the edge-state transmission coefficient η, obtained by modelling the out-of-plane AMR of the different devices.

Metal–insulator transition

Having established the signatures of edge-state transport, we now examine their evolution as the magnetization is tilted in-plane by an external 1 T radial field (Fig. 2). The temperature dependence of ρ in Fig. 2a shows a metal–insulator transition as the magnetization is tilted from out-of-plane to in-plane, irrespective of the azimuthal angle (also see Supplementary Fig 5). Similarly, in the gate voltage dependence (Fig. 2b), as the chemical potential is tuned into the magnetic gap, the drop in ρ is replaced by a sharp rise when the magnetization is in-plane. This is indicative of a crossover from dominant 1D edge transport to diffusive transport, as the system transitions from a QAH regime to a gapless ferromagnetic topological insulator (when magnetized in-plane). Figure 2a,b already reveals a giant AMR effect when the chemical potential is in the magnetic gap, and this crossover is the origin of the novel angular dependence of the AMR that we now discuss.

Figure 2

Magnetization tuned edge-state transport.

An external radial field of 1 T is used for all the data in this figure. (a) Transition from metallic to insulating behaviour in the temperature dependence of ρ at VG=VG0, as the magnetization is tilted from out-of-plane to in-plane. The red, green and blue curves correspond to the x, y and z directions, respectively. (b) Gate voltage dependence of ρ at T=280 mK when the film is magnetized along the x (red),y (green) and z (blue) directions. (c) Out-of-plane anisotropic magnetoresistance of device A (light red, circles) at charge neutrality, as the field is rotated in the x−z plane. The angular dependence is unexplained by a conventional cos2θ fit (red, dotted line) and is understood by fits (red, bold line) that use an edge-dissipation mixing model. AMR data (circles) and corresponding fits (bold lines) from device B (light blue) and device C (light green), described in the text. (d) Cartoon of the device geometry highlighting field rotation from out-of-plane to in-plane, in the x−z plane. (e) In-plane anisotropic magnetoresistance of device A (light red, circles) at charge neutrality, as the field is rotated in the x−y plane. The angular dependence is explained by a typical cos2ϕ fit (red, bold line). (f) Cartoon of the device geometry highlighting field rotation in the x−y plane.

Giant anisotropic magnetoresistance

We systematically probe the crossover from edge dominated to diffusive transport as we rotate a 1 T external magnetic field from perpendicular-to-plane to in-plane. For a 1 T field, the AMR, defined as (ρ−ρ)/ρ takes a maximum value ∼140%, and its angular dependence is completely inconsistent with the well-known cos2θ angular dependence seen in conventional ferromagnets. To explain these observations, we model the angular dependence using a simplified four-terminal Landauer–Büttiker formalism I=ΣG(V−V), where I and V are the current through and the voltage at the ith contact, and G is the conductance of the channel between the ith and jth contacts. Our model accounts for the coexistence of both dissipationless edge channels and dissipative bulk and surface channels. The dissipative channels likely arise from the gapless, two-dimensional surface Dirac states of the side walls when the chemical potential is in the magnetic gap20. We note that these surface states are not gapped because the external magnetic field and the magnetization are always parallel to the plane of the side walls. We cannot exclude other possible origins of dissipation, such as the bulk bands, mid-gap impurity bands, and even the gapped surface states. These details are however not important for the model. The conductance G is given by G=T(e2/h) where the transmission coefficient T takes the form T=ηδ+t. Here, η is the contribution from edge modes that only transmit between leads j and j+1, and t corresponds to dissipative contributions. From calculations explained in detail in Supplementary Note 4, we derive the following expression for the magnetoresistance:

where Rd is a measure of the dissipative resistance (and is the dissipative transmission coefficient). Physically, the first term in the denominator is the conductance of an inter-mixed channel, that constitues the edge states scattering into dissipative pathways. The second term though corresponds solely to the dissipative channels. The transmission coefficient for the edge modes η depends on the gap Δ induced by the magnetization. We expect η to be suppressed from its maximum value of 1 as the magnetic gap is reduced Δ: this arises from enhanced inter-mixing with dissipative channels that accompanies an increase in penetration length λ of the edge-state wave functions, where λ=ℏvF/Δ where vF is the Fermi velocity. We thus assume a phenomenological expression for the edge-state transmission as:

where l0 is a length scale that quantifies the strength of the inter-mixing, and Δ0 is the corresponding energy scale. This expression satisfies the physical requirement at the two extrema: Δ→0η→0 and Δ→∞η→1. The angle dependence of the magnetoresistance then arises via the dependence of the magnetic gap on the perpendicular component of field, that is, Δ=Δmcosθ. Δm is the size of the magnetic gap, when the film is magnetized out-of-plane. The physics behind the model is depicted in Fig. 3. We now use equation (1) to fit our angular data with two free parameters: Rd and Δ′ (defined as Δ′=Δ0/Δm). The model clearly fits the data very well, up to tilt angles ∼70°, as seen in Fig. 2b, confirming that the novel, giant AMR effect has its origins in the edge states. The model deviates at large tilt angles when the field approaches the in-plane orientation. In this regime, the dominant resistance is associated with the dissipative channels, whose angular magnetoresistance is not accounted for in our model (Rd is treated as a constant).

Figure 3

Edge-dissipation mixing model.

A pictorial representation of the physical origins of the giant AMR effect. The penetration length of edge states (depicted in dark blue) is dependent on the size of the magnetic gap Δ, which in turn depends on the perpendicular-to-plane component of magnetization. (a) As the magnetization (red arrows) is tilted in-plane from the intial perpendicular-to-plane configuration, the size of the magnetic gap decreases (illustrated by the energy-momentum dispersions, in red). Consequently, the penetration length of the edge states into the bulk (depicted in light blue) increases (b) eventually leading to complete mixing as the magnetization is tilted perfectly in-plane (c).

Beyond explaining the angular dependence of this novel AMR, our model also presents a powerful approach to quantify edge-state contributions to transport in an edge-dissipation mixed QAH insulator. The fitting parameter Δ′ can be used to extract the edge transmission coefficient η=1−e−1/Δ′. This is used to quantify and compare edge transport contributions in two other samples with varying degrees of dissipation, whose AMR is also shown in Fig. 2c. The strength of the dissipation is inferred by comparing the zero-field ρ and the magnitude of the AHE in these films. The angular dependence and fits to data from devices B and C, also serves as a test of reproducibility for the giant AMR. Table 1 summarizes the different device characteristics and the extracted fitting parameters, all at charge neutrality and at the same temperature. The extracted values of η consistently characterize the degree to which the samples are close to the ideal QAH regime. Also, it is important to point out that the magnitude of the AMR effect consistently depends on the relative values of η and (equation (1)).

In-plane anisotropic magnetoresistance

As we mentioned earlier, conventional ferromagnets show an AMR described by a cos2ϕ angular dependence. This behaviour is indeed observed for in-plane measurements (Fig. 2e). Compared the giant AMR, the in-plane AMR only has magnitude of ∼10.7%, at charge neutrality (Fig. 2e), and also shows a far weaker temperature dependence (difference between upper two curves in Fig. 2a). With the field and magnetization in-plane, the sample is expected to be in a diffusive regime, since tuning the surface states gapless destroys the edge states. The functional dependence of the in-plane AMR is unchanged with gate voltage (Supplementary Fig. 7).

Gate voltage and temperature dependence of giant AMR

Finally, we present a more detailed picture of dependence of the giant AMR on the gate voltage and temperature. This is shown for device A in Fig. 4. The effect is electrically tunable over an order of magnitude, and is largest in magnitude with the chemical potential in the magnetic gap, when edge-state contributions to transport are maximum. As the chemical potential is tuned out of the magnetic gap, surface state and bulk dissipative channels are populated, reducing edge-state contributions and consequently the magnitude of the AMR. This is captured in the gate dependence of the edge-state transmission coefficient η (Fig. 4b), extracted from our fits to the data of Fig. 4a. The qualitative behaviour of is also reasonable, since the dissipative transmission coefficient increases as more electrons are accumulated. Similarly, the temperature dependence of the AMR (Fig. 4c) also captures the competition between the edge states and dissipative channels since our model quantifies edge transport in temperature regimes away from perfect quantization (Fig. 4d). While perfect Hall quantization has required dilution fridge temperatures, our results here demonstrate that the temperature scale for observing this novel AMR, with its origins in the edge states can be significantly higher. Furthermore, both the gate and temperature dependence of the fitting parameters reveal the interplay between the edge states and dissipation channels: an increase in dissipative transmission is accompanied by a reduction in edge transmission and vice versa.

Figure 4

Giant anisotropic magnetoresistance.

(a) Out-of-plane AMR of device A at T=280 mK for a range of gate voltages. The AMR takes maximum values when the film is tuned to charge neutrality. (b) Gate voltage dependence of the edge-state transmission coefficient η (red) and the dissipative transmission coefficient (blue) at 280 mK, extracted from fits to the data of a. (c) Out-of-plane AMR at VG0 for a range of temperatures. (d) Temperature dependence of η (red) and (blue) at VG=VG0 extracted from fits to the data of c.

Conclusions

In summary, we have probed the evolution of edge states in a QAH insulator, as its magnetization is tilted in-plane, both by angular magnetoresistance measurements and by developing a model based on the Landauer–Büttiker formalism that accounts for coexisting dissipative channels. We demonstrate a crossover from dominant 1D edge transport to diffusive transport driven by the closing of the surface state magnetic gap. This crossover results in a novel, gate-tunable, giant AMR effect that serves as a quantitative probe of the interplay between the chiral edge modes and dissipative channels, even in non-ideal QAH regimes. This methodology is relevant to ongoing efforts to raise the temperature scale of the quantum anomalous Hall effect21. Furthermore, as with recent advances that exploit the interplay between spin polarized states in topological insulators and ferromagnetism222324, the observation of a giant AMR effect may provide a useful basis for proof-of-concept devices for ‘topological spintronics'.

Methods

Material synthesis and characterization

The Cr(Bi,Sb)2−Te3 thin films studied by transport were grown under ultra-high vacuum (low 10−10 torr) by molecular beam epitaxy on 〈111〉 STO. SQUID magnetometry was however performed on Cr doped thin films grown on <111> InP substrates. (See Supplementary Fig. 12, Supplementary Note 5). The source material was high purity Cr, Bi, Sb and Te thermally evaporated from conventional Knudsen cells. For device C however, e-beam evaporation of Cr was employed instead. Before growth, the STO substrates were subjected to an ex situ high temperature (925–950 °C) anneal in an oxygen atmosphere to obtain a smoother substrate morphology. The topological insulator film of device A was capped in situ with a thin layer of room-temperature-deposited Al, which naturally oxidized upon exposure to atmosphere. Devices B and C were uncapped. The thickness of the topological insulator films was estimated by comparison with the growth rate of calibration samples, whose thickness was measured using high resolution transmission electron microscopy and also confirmed by atomic force microscopy of the Hall bar step edge. The relative compositions of Cr, Bi, Sb and Te in the different films were obtained by high resolution energy dispersive X-ray spectroscopy and secondary ion mass spectroscopy.

Transport measurements

The thin-film devices were mechanically scratched into Hall bars using a probe station. Since the Hall bars were mechanically defined, our raw data were symmetrized to account for geometric offsets. The angular measurements were performed in an Oxford Triton He3 vector magnet system with a radial field limit of 1 T and a base temperature of 280 mK. All transport measurements typically employed standard low frequency (19 Hz) lock-in techniques with source currents in the range 10–20 nA, and a Keithley 6430 source meter for back gating.

Additional information

How to cite this article: Kandala, A. et al. Giant anisotropic magnetoresistance in a quantum anomalous Hall insulator. Nat. Commun. 6:7434 doi: 10.1038/ncomms8434 (2015).