PHYSICAL REVIEW RESEARCH 3, 013058 (2021)
Universal quantization of the magnetic susceptibility jump at a topological phase transition
Soshun Ozaki * and Masao Ogata Department of Physics, University of Tokyo, Bunkyo, Tokyo 113-0033, Japan
(Received 21 April 2020; accepted 16 November 2020; published 19 January 2021)
We examine the magnetic susceptibility of topological insulators microscopically and find that the orbital- Zeeman (OZ) cross term, the cross term between the orbital effect and the spin Zeeman effect, is directly related to the Berry curvature when the z component of spin is conserved. In particular, the OZ cross term reflects the spin Chern number, which results in the quantization of the magnetic susceptibility jump at the topological
phase transition. The magnitude of the jump is in units of the universal value 4|e|μB/h. We also apply the obtained formula to an explicit model and demonstrate the quantization. For this model, the physical origin of this quantization is clarified.
DOI: 10.1103/PhysRevResearch.3.013058
I. INTRODUCTION g0 = 2. (Here we neglect the relativistic correction of g0.) Thus, their conclusion means that the jump in χ is not Topological insulators (TIs) [1–13] show anomalous phe- OZ quantized in a universal value. nomena such as electric conduction on sample surfaces. In the present paper, we study χ microscopically based Experimentally, the search for candidate materials for TIs OZ on the Green’s function formalism and show that, in contrast is one of the most important problems. In particular, two- to the results of Nakai and Nomura, the jump in χ is exactly dimensional (2D) TIs are predicted to show unique phenom- OZ quantized in units of the universal value 4|e|μ /h [see Eq. (9) ena, such as the spin Hall effect and robust edge states against B below] even in the presence of SOI as long as the z component nonmagnetic impurities, only a few of which have been found of the spin is conserved. Here μ =|e|h¯/2m is the Bohr mag- [7–11]. So far, the confirmation of topological materials has B neton. When we study a model microscopically [as in Eq. (1)], been achieved by finding the edge state by angle-resolved the modification of the g factor does not occur explicitly, photoemission spectroscopy or from the transport coefficients. and instead the effect of SOI appears in the deformation of Since both methods detect anomalous electronic states at the the Bloch wave functions and the energy dispersion, which edge, it is desirable to develop some bulk-sensitive methods eventually leads to the orbital-dependent g factors. We show that enable us to confirm the topological nature of a material. below that the effect of SOI is exactly canceled out in χ , In this paper, we propose that the quantization of the bulk OZ which leads to the quantization of jump with a universal value. magnetic susceptibility jump can be used as strong evidence In the latter part of this paper, we apply the obtained formula for the topological phase transition in 2D TIs. to a low-energy effective model for a topological insulator Usually, the magnetic susceptibility is discussed in terms to show the validity of the present proposal. For this model, of the orbital effect of the magnetic field [14–30] and spin we clarify the physical origin of the quantization: It turns out Zeeman effect independently. In general, however, there can that the quantization is associated with the chiral edge current, be a cross term between the orbital and Zeeman effects which is characteristic of the topologically nontrivial state. [5,6,31–35], which we call the orbital-Zeeman (OZ) cross term χOZ in the following. Recently, Nakai and Nomura [34] II. GENERAL FORMALISM discussed the jump in χOZ at the topological phase transition using the formula of the orbital magnetization [36–42] and First, we develop microscopically a general formula for the Stredaˇ formula [43]. They calculated the OZ cross term magnetic susceptibility including orbital magnetism, Pauli in the Bernevig-Hughes-Zhang model [4] and concluded that paramagnetism, and the OZ cross term in terms of thermal the width of the jump depends on the g factors of the involved Green’s functions in the presence of SOI. Let H be the general orbitals introduced phenomenologically. In general, spin-orbit Hamiltonian derived from the Dirac equation in the presence interaction (SOI) modifies the g factor from its bare value of a general 3D periodic potential V (r) and a magnetic field, which is given by 1 eh¯ H = [p − eA(r)]2 − σ · B(r) + V (r) *[email protected] 2m 2m h¯2 h¯ Published by the American Physical Society under the terms of the + ∇2V + σ ·∇V × [p − eA(r)], (1) 2 2 2 2 Creative Commons Attribution 4.0 International license. Further 8m c 4m c distribution of this work must maintain attribution to the author(s) where A(r) is a vector potential, e < 0 for electrons, σ = and the published article’s title, journal citation, and DOI. (σx,σy,σz )are2× 2 Pauli matrices, and B(r) =∇×A(r)
2643-1564/2021/3(1)/013058(9) 013058-1 Published by the American Physical Society SOSHUN OZAKI AND MASAO OGATA PHYSICAL REVIEW RESEARCH 3, 013058 (2021) represents a magnetic field. The last term represents the SOI. perpendicular to the 2D x-y plane. Eventually, we obtain the It is to be noted that the second term representing the Zee- formula for 2D systems by substituting the 3D k summation man interaction has the bare value g0 = 2. As performed for a 2D k summation and by including trace over the subband by Fukuyama [20], we implement a perturbative calculation index for the z direction wave function. Hereafter, we use bold of the free energy in terms of the vector potential A(r)via for two-component vectors, e.g., k = (kx, ky ). the Luttinger-Kohn representation [44]. As a result, we ob- Then, to discuss the quantization, we rewrite Eq. (2c) tain the expression for each contribution as follows (the details in terms of the Bloch wave functions in a similar way to of the derivation are shown in the Appendix): Ref. [22]. The periodic part of the Bloch wave function t , = ↑ , , ↓ , 2 uˆlk(r z) (ulk (r z) ulk (r z)) satisfies χ = e kBT γ Gγ Gγ Gγ G, orbit Tr x y x y (2a) , = ε , , = −ikr ikr, 2¯h2 L3 Hkuˆlk(r z) l (k)ˆulk(r z) Hk e He (3) nk where l represents both the band index in the 2D Brillouin kBT χ =− Tr MsGMsG, (2b) zone and the subband index in the z direction. In the fol- Pauli L3 z z nk lowing, we abbreviate ulkσ (r, z)asulkσ . When the SOI is / 2 2 σ ∂ − − ∂ − i|e| k T given by (¯h 4m c ) z[( xV )(py eAy ) ( yV )(px eAx )], χ =− B sGγ Gγ G − sGγ Gγ G , OZ Tr M x y M y x the z component of spin is conserved. In this case, up- and h¯ L3 z z nk down-spin electrons are independent, and ulk↑ and ulk↓ are the (2c) eigenstates of σz with eigenvalues εl↑(k) and εl↓(k) (denoted ε ↑ ε ↓ G as l and l in the following). In this case, the matrices and where G is the thermal Green’s function G(k, iεn ), whose (ll ) s −1 M are diagonal and given by [Gσ ]ll = δll (iεn − εlσ + μ) component is the matrix element between the lth and lth z and bands. Each band index includes the pseudospin degrees of s ∗ freedom in the case with SOI. εn is the Matsubara frequency, = −μ σ σ Mz lσ,lσ ulkσ ( B z )ul k drdz γμ represents the current operator in the μ direction divided by / s −μ σ e h¯, and Mz is the matrix for the operator B z. The effect =−μBσδll δσσ , (4) of SOI is included in G and γμ. Tr is the trace over the band indices and the spin degrees of freedom. In Eq. (2), χ and respectively. Note that the normalization is fixed by the equa- orbit tion χPauli represent the orbital and Pauli magnetic susceptibility, respectively. These are the same expressions as were obtained ∗ u u σ drdz = δ δσσ , (5) before [20,45] even in the presence of SOI. On the other hand, lkσ l k ll χ OZ is the OZ cross term, which we focus on in this paper. where the spatial integral is over the entire crystal. On the other hand, the matrix γμ has off-diagonal matrix elements III. UNIVERSAL QUANTIZATION OF χ between the different bands and it becomes [20,46], OZ ∗ ∂Hk In the following, we consider the case where the z com- γ = [ μσ ]ll ulkσ ul kσ drdz ponent of spin is conserved even in the presence of SOI. ∂kμ This situation will be realized more easily in 2D systems ∂εlσ ∗ ∂ulkσ = δ + ε − ε . than in 3D systems. Therefore, we focus on 2D systems in ll ( l σ lσ ) ulkσ drdz (6) ∂kμ ∂kμ the following. As shown in Appendix, it is straightforward to extend Eq. (2) in 2D systems where V (r) gives a con- Substituting these quantities into Eq. (2c) and carrying out the finement potential in the z direction and the magnetic field Matsubara summation, we obtain