On Behind the Physics of the Thermoelectricity of Topological Insulators

Topological Insulators are the best thermoelectric materials involving a sophisticated physics beyond their solid state and electronic structure. We show that exists a topological contribution to the thermoelectric effect that arises between topological and thermal quantum field theories applied at very low energies. This formalism provides us with a quantized topological mass proportional to the temperature T leading, through an electric potential V, to a Seebeck coefficient where we identify an anomalous contribution that can be associated to the creation of real electron-hole Schwinger’s pairs close to the topological bands. Finally, we find a general expression for the dimensionless figure of merit of these topological materials, considering only the electronic contribution, getting a value of 2.73 that is applicable to the Bi2Te3, for which it was reported a value of 2.4 after reducing its phononic contribution, using only the most basic topological numbers (0 or 1).

Nowadays topological insulators (TI) are the best thermoelectrics (TE) at room temperature [1][2][3][4][5] , specially if they are combined with nanotechnological structures able to reduce the phononic thermal conductivity. A good example is 6 bismuth telluride, Bi 2 Te 3 , which has a small band gap giving a good number of carriers at room temperature (300 K) and reaching 2.4 for its dimensionless figure of merit ZT for p-type using alternating layers in a superlattice with Sb 2 Te 3 . This is the highest value of thermoelectricity observed 7 so far at room temperature. Despite the fact that the electronic structure of these materials was exhaustively studied 8 in relation with their thermoelectricity 9,10 , there is still lacking in the literature a physical model 11 able to explicitly take into account their common topological and physical features. Its importance might appear obvious with a counterexample, Pb 1−x Sn x Te has a good electronic structure to be a topological insulator with thermoelectricity 12 , but due to have an even number of band inversions, this prevents it to have time-reversal symmetry T and thus to be a TI, although it is a good thermoelectric at higher temperatures. Hence, the whole topology of a TI is not fully necessary for having good thermoelectricity as we are going to see. Understanding the physics behind these phenomena is not an easy task because it links different scientific branches which were developed independently: particle physics, statistical mechanics, condensed matter and algebraic topology. This is a characteristic of materials which exhibit linear dispersion laws instead of quadratic ones, allowing a quantum field interpretation where the spinors play a fundamental role substituting the usual non-relativistic wave function.
The paper is organized as follows. Firstly, we examine the basic concepts of topology and physics for topological insulators trying to show how they are related in the same structure. After that, we show that the Riemann-Hurwitz formula plays an important role which was not considered in the literature so far. This allows us to find five topological regions which are connected by four bands defined using the periodicity of the instanton solutions associated to the non-Abelian Berry fields introduced within the bulk of the TI. In the case of no time-reversal symmetry T , we have only three topological regions, connected by two bands. Finally we present a straightforward relationship between the temperature T and the topological index μ with the scalar electric potential V. This leads to a Seebeck coefficient for which we identify two terms, one related to the topological electron pump, and another associated to a change in the topological index that might be associated to the creation of real electron-hole pairs as will be analysed. We end with the calculation of a general expression for the dimensionless figure of merit ZT of the edge states in TI, taking into account the electronic contribution and neglecting the phononic part 7,13 .

Results
Solid state physics allow us to tackle the problem of a crystal with translation symmetry, reducing the analysis of its different physical properties to the first Brillouin zone. This is permitted thanks to Bloch theorem, i.e., for a periodic potential V(x + a) = V(x), where a is the spatial period or the lattice constant, the wave function associated to the electrons have also a periodicity ψ k (x + a) = e ika ψ k (x). Being the eigenvalues ξ k = ξ k + T , ξ k = ξ −k also periodic in time. These are the necessary conditions for calculating the energy bands in a solid. It is easy to see that parity P , or space inversion is straightforwardly followed, while time-reversal symmetry T is not so obviously fulfilled. For example, the Schrödinger equation is not a solution due to the first order time derivative. This can be solved if the T operator has also associated a complex conjugation K operator. In fact, we must define σ = −T i K y for spinless states where U is a unitary operator. In the case of having half-integer spin particles, the unitary operator can be written in function of the σ y Pauli matrix as . This, given that σ = 1 y 2 , allow us to write the time-reversal operator as σ = −T i K y , which acting on the multiparticle state gives T 2 = 1 for an even number of fermions or T 2 = −1 when the number is odd. More generally written, T 2 has eigenvalue (−1) 2s for a particle of spin s, and if |n> is an energy eigenstate, then T|n> is also an eigenstate of the system, sharing the same energy and being orthogonal to each other. Thus, if there is an odd number of electrons there must be (at least) a twofold degeneracy; known as Kramers degeneracy. This is Kramers theorem, which when it is completed with the previous Bloch theorem for the bands, should provide us the first tools 14 to examine topologically the first Brillouin Zone (BZ) of a TI.
In Fig. 1 we represent schematically how the 2D BZ evolves under different translation operations transforming the square BZ into two equivalent cylinders S 1 . Consequently, it is also easy to see that combining the two lattice periodicities available in two dimensions, a 2D torus S 1 × S 1 can arise. That is, the existence of translation symmetry allows us to transform a topologically trivial square into a non-trivial torus with genus g ≠ 0. Given that all crystals are able to develop these fundamental properties, non-trivial topology is not enough when we study TI and the presence of singularities which are usually counted on the band structure are needed. Obviously, without extra-information, these singularities break the translation symmetry or the periodicity of the Bloch states 15,16 . So, let's analyze a little bit the possibility of these translations symmetry breaking. Physically the infinitesimal translations are generated by the linear momentum operator 1  follows the Heisenberg's principle of uncertainty. By extending this operation over the whole crystal one gets the finite trans- . Taking into account the periodicity of the crystal potential V(x) = V(x + a), which involves the periodicity of the states, we can write L a |ψ(x)> = |ψ(x−a)> = exp( ika)|ψ(x)> where k is the wavelength number associated to the momentum p = ℏk. On the contrary, the presence of singularities makes that |ψ(x − a)> = exp(ika)exp(iγ)|ψ(x)>, being γ the Berry phase which gives to electrons an anomalous phase factor during a complete cycle in the order parameter. Thus, in order to maintain the periodicity of the Bloch states, preserving the non-trivial topology that results for having a torus, singularities must take place on the edges of the crystal, where translation symmetry is no longer satisfied. Since singularities manifest in energy bands by means of Fermi points, bands are good hosts for feeding the non-trivial topology of a material whose electronic structure is appropriate and which its surface shows the breaking of translation symmetry. In this way, the Berry phase joints the non-trivial topology of the crystal employing its curvature and connection on the bands 17 .

Riemann-Hurwitz formula and its application to topological insulators. The topological elements
of the TI have been found, but now it remains to see how they work together employing their associated invariants. The Riemann-Hurwitz formula 18 , which generalizes the Euler topological invariant, enable us to construct www.nature.com/scientificreports www.nature.com/scientificreports/ an equation relating the genus g and g' of two compact surfaces, i.e., whose boundary is zero. Actually this formula establishes the conditions for a map f:T → S being surjective and holomorphic, reducing the several topological invariants introducing the genus of T, the genus of S, N the degree of the map f and the number of ramifications e f (k). Riemann used the mentioned formula in the case that S were zero, i.e., spheres. Much later the general proofs were obtained by Zeuthen and Hurwitz. The formula is In our case we have a 2D torus T 2 with genus g T = 1 mapped in a S 3 sphere with genus g S = 0 and the degree of the map N = 2 due to the Kramers double degeneracy, see Fig. 2. Hence we obtain five ramifications branches, i.e., e f (k) = 5. This can be directly related to the second Chern number, playing a fundamental role in the transformation of heat in electricity. The number of ramification branches diminishes to three when there are not Kramer pairs, that is, N = 1 with e f (k) = 3. But it is very remarkable to observe that this formula doesn't depend on the dimensions of the involved sphere or torus, which justifies us to work with a T 2 torus instead of a T 3 or T 4 without being worry about new results.
In the presence of singularities in the band structure the map c of Fig. 2 can be interpreted as the d a map within the Hamiltonian H introduced 19 to study TI in (4 + 1)-D,  which is no more than the winding number resulting from the map of a four-dimensional torus T 4 to a sphere S 4 , having the mass m associated to the spinor's Γ's five critical values given by the condition allowing us to identify, in the same way as the Riemann-Hurwitz theorem, five different branches for the second Chern number where p is taken as a mass parameter which must be equal to the background kinetic energy of the particles k T v 2 / B F 2 for keeping its physical meaning. In this way, p is mainly associated with the quantization of the temperature 20 . It is easy to see that for non-T symmetric Hamiltonians there are only three ramifications available for the second Chern number. As we saw, this can be interpreted, in one straightforward form, into the Riemann-Hurwitz formalism as two different maps, where we change the degree of the map N from 2 to 1.

Topological Seebeck coefficient.
Physically the non-Abelian Berry phase takes into account the allowed bulk degenerate states which are directly related with the change of temperature, as we shall see soon. First, we take the definition of the non-Abelian Berry connection and second Chern number where the indixes of the Levi-Civita tensor stand by i, j, m, n = 1, 2, 3, 4 and α refers to the occupied bands. This can be written f Tf i i being the T i 's the generators of the inner symmetry group, which in our case will be SU (2) The first term of the integral is a total derivative which can be made zero in a manifold without boundary, while the second, is an integral written in short as w(U), which provide us topological information of the manifold as a winding number. This integral is actually an integer number coming from the homotopy group π 3 (S 3 ) for the SU(2) group and we can take the Abelian Chern-Simons gauge symmetry associated to the background electrodynamics, i.e., the U(1) as a subgroup of transformations close to the surface. In this case, we can rewrite the Chern-Simons transformation for both fields. Hence the transformation of the Chern-Simons can be rewritten as 2 But as the path integral e i W  must be invariant, hence it means that μ is an integer. These are going to be the dynamic fields meanwhile the electromagnetic gauge potential A α will appear as a background gauge field.
Let us go now to the thermodynamic part of the topological thermoelectricity, applying the instanton solution associated with the second Chern number developed above. Within this mathematical contest, it is possible to make a direct relationship between thermodynamics and quantum formalism. Feynman path integral give us an expression for the amplitude of probability to evolve a particle from x i (0) to x f (t) a time later by , where S is the classical action, which could be obtained by On the other hand, in Statistical Mechanics the partition function is defined in quite a similar form by There is a way to change from one to the other using the Wick trick, i.e., doing the time a pure complex variable and transforming the Minkowskian space-time in a Euclidean one. Thus, we have where S E is the new Euclidean action. This fixes the concept of temperature relating it directly with the time by τ = ℏβ making where we have assumed a cyclic motion x(0) = x(β), i.e, the particle must come back where it started after Euclidean time τ = ℏβ 22,23 . This is exactly how Berry's phase works and justify partially its introduction as a gauge field.
On the other hand, Berry's phase is directly related with the electric polarization P, i.e., where μ is a winding number quantizing the temperature. This allows to find the electric potential directly related with the temperature by B which means that we have transformed the electric potential into another, plus a μ π β e 2 thermal term. This turns out to be a fundamental result: the thermal energy appears quantized by the winding number being added to one electric potential under a gauge transformation. Notice that a Chern-Simons term was necessary since it is not gauge-invariant and appears as surface actions. This allows to have, in non-Abelian Chern-Simons, the coupling constant g, directly associated to the temperature, related to the Chern-Simons action transforming S CS → S′ CS = S CS + 2πℏ 2 μ/g 2 where μ is the topological mass.
Getting back to the topological electric potential Eq. (7), we can easy calculate the Seebeck coefficient choosing V′ = 0. We can identify two contributions, the first, that comes purely from the temperature gradient in a topological branch, and the second, which takes into account the contribution due to a change in the topological index μ, i.e, a jump between ramifications. As we will see, the first term is equivalent to that we find when we compute the surface Seebeck coefficient while the second, doesn't appear since there is only one ramification available on the surface. Once we have seen that the instanton solutions allow us to relate the electric potential V with the temperature T and hence the electric field E with ∇T, we can give a microscopic picture of how the thermal energy is invested in electricity. The main idea is that there is a creation (Fig. 3) of electron-hole Schwinger's pairs [23][24][25][26] , provided that the electric field is big enough. Calculating the critical electric field E c , following Heisenberg-Euler 27 , we notice that the limit for obtaining a great number of real electron-hole pairs is www.nature.com/scientificreports www.nature.com/scientificreports/ where we have considered that the energy gap is 0.21 eV and the Fermi velocity ×  v m s 6 10 / . This provides a critical electric field of almost ten orders of magnitude lower than the critical electric field in QED and with one equivalent temperature of ζ = 1.74 × 10 −6 Knm −1 , being these values accessible in these topological materials at so small scales as at hundredths of volt at distances of angstroms. In this way, we can rewrite the second term of the Seebeck coefficient considering that the dependence of μ with the temperature, as we show in Eq. (4), is represented by a Heaviside function, leading to the following expression for the Seebeck.
being p the different values where μ changes, that is, 0, ±m/2, ±m/4 and where the local increase Fig. 4, represented by the second term, can be interpreted as the contribution originated in the creation of real electron-hole pairs. This Eq. (9) of the Seebeck coefficient has two terms, the first one quantize S in integers due to the Chern-Simons topological mass, whereas the second corresponds to the variation of this mass respect to the temperature T. The relevant point is that the topological bands allow finding a new term in the Seebeck coefficient which can increase it depending on the value of the temperature. But what is more important, these strings of singularities separate thermally some regions from other, resulting in a temperature gradient (∇ → T ) which produce an electric field → E 2 capable to create Schwinger's pairs whose electrons and holes have different velocities depending on the level where they are situated respect to the Fermi level. Thus, unlike what happens with metals, the effect of temperature on electrons and holes can be quite different without cancelling each other, leading to a Seebeck coefficient that would be also higher than in a semimetal. It is fundamental to observe that without the topological ramifications we would have a homogeneous crystal without thermoelectricity.
Surface transport and dimensionless figure of merit ZT. Finally, let's calculate the topological contribution to the surface figure of merit in TI, limiting our calculations to (2 + 1) dimensions, where the topology is determined through the first Chern number n and the winding number n. We start considering that, in a given parameter m is far away from the boundary between two topological branches, that is, for example, 2p. When m is close to these value the jump in the topological mass μ generates a potential difference that is able to break the vacuum and to generate a real electron-hole pair representing a new mechanism to transform thermal energy into electric in solid state physics. (B) Representation of electron and holes creation at the hot (red) side of the TI where m ≈ 2p. Electrons and holes take place on different branches, leading to an electric field → E 2 . Due to charge carriers are responsible to thermalize the material, this originates electron J e and holes J h currents where each type of carriers have different mobilities since they below to different branches. www.nature.com/scientificreports www.nature.com/scientificreports/ direction x or y, exists a temperature gradient in our TI in such a way that eigenstates evolve adiabatically. We can define then the Seebeck coefficient through the entropy as where ψ b and ψ t are the final an initial states, θ is the angle between the states and n is the winding number which can only take odd integers values, otherwise the topological Seebeck coefficient would give zero. In fact, we are representing the entropy of a Riemann surface. Notice that although the von Neumann entropy would give zero considering pure states, this is not true considering that these states are in entanglement due to their global topological properties. Actually there are many local Wannier bands, i.e., , providing separate electrons which move coherently producing what is called Thouless charge pump 28 (see Fig. 5) and relates Berry phase with the electric polarization P [29][30][31] (see Supplementary Information). In fact, due to the previous definitions, it is easy to see that = being θ the angle directly associated to the Berry phase on the closed curve C. This leads to an expression very similar to that obtained in (4 + 1)-D, first term of Eq. (9), where the product μw(U) has been transformed into nn, where π is the Berry phase of a non trivial material, n counts the number of times we complete a cycle in a system (number of singularities) and n is the first Chern number which takes into account the whole topology on the Brillouin zone for our TI.
In order to complete the calculation of the surface figure of merit, we proceed to the calculation of the electronic thermal conductivity and the electrical conductivity. For the electronic thermal conductivity, we consider the 2D density of states in the semimetal region as which allows us to obtain the electronic thermal conductivity κ e as Being f(ξ) de Fermi-Dirac distribution function, ζ(3) the Riemann zeta function of dimension three (Apery's constant) and where we have supposed a temperature dependent mean free path l. On the other hand, given that we have a ballistic regime for electronic transport, its conductivity σ e appears given by the simple expression σ = n e e h 2 . Therefore, despite using so different expressions than the ones of the metals, where a quadratic dispersion equation is employed instead of the linear one of the semimetal, we obtain a good Wiedemann-Franz law yielded by where the Lorenz number ζ = L nk e 3 (3)/ ( / ) B 2 is one constant, as it ought to be, but divided by a Chern number n which tell us that this expression is only valid within the context of the non-trivial topological materials that we have considered. Finally, we can calculate the figure of merit Z for these topological insulators where we are not considering the phononic part of the thermal conductivity 13,32 . In this way, the dimensionless figure of merit turns out to be a simple expression 2 3 2 This is the extra pure topological figure of merit for the edge states, which is zero in the case of trivial topological materials. Although these conditions are quite ideal and transport constraints can diminish its efficiency under real physical features of each material, this result opens a great hope because it tells us how to improve www.nature.com/scientificreports www.nature.com/scientificreports/ highly the thermoelectricity associated to the topological materials. In the case of Bi 2 Te 3 7 , for the quantum numbers equal to one we obtain a value close to the one of its present maximum, i.e. ZT = 2.737.

Discussion
In summary, we have shown the relationship between topological insulators, such as the family Bi 2 Te 3 , as well as other topological related materials without time inversion protection as the Pb 1−x Sn x Te, and their associated thermoelectricity. We have also seen that the second Chern number obtained for the non-Abelian SU(2) field leads to a thermal topological mass on a Chern-Simons action. This is equivalent to have a quantized temperature working in a kind of topological bands that we define as the ramification branches using the Riemann-Hurwitz formula on a Euclidean spacetime where instanton solutions substitute Bloch oscillating states. Physically what we have is a pumped charge between bands connected by the non-Abelian Berry phase within the insulator bulk at low temperature with an electromagnetic background field on the surface. Therefore, close to the surface we have only an Abelian U(1) Chern-Simons term providing us with one transformation between electric and thermal energy because we have only one kind of states. Moreover, we show that the Schwinger's electron-hole pairs, close to the topological bands, produce an increase of the Seebeck coefficient contributing to the transformation of thermal into electric energy which is one of the key points of the model that we present in this paper. Finally, we calculate a general expression to the dimensionless figure of merit in terms of the Chern number and winding number, getting a value that coincides quite well with the one experimentally measure for the Bi 2 Te 3 , doing zero its phononic thermal contribution. It is open for future a new class of topological materials using topological indices higher than one which can cross what is considered nowadays the efficient critical value of four for the ZT figure of merit changing the physical conditions suggesting in the presented model.