arxiv: v2 [cond-mat.str-el] 4 Nov 2010

Transcription

1 The Kelvin Formula for Thermopower Michael R. Peterson Conense Matter Theory Center, Department of Physics, University of Marylan, College Park, MD 2742, USA B. Sriram Shastry Physics Department, University of California, Santa Cruz, CA 9564, USA (Date: August 29, 218) arxiv: v2 [con-mat.str-el] 4 Nov 21 Thermoelectrics are important in physics, engineering, an material science ue to their useful applications an inherent theoretical ifficulty. Recent experimental interest has shifte to strongly correlate materials, where the calculations become particularly ifficult. Here we reexamine the framework for calculating the thermopower, inspire by ieas of Lor Kelvin from We fin an approximate but concise expression, which we term as the Kelvin formula for the the Seebeck coefficient. Accoring to this formula, the Seebeck coefficient is given as the particle { number N erivative of the entropy S, at constant volume V an temperature T, S Kelvin = 1 S }. This q e N V,T formula is shown to be competitive compare to other approximations in various contexts incluing strongly correlate systems. We finally connect to a recent thermopower calculation for non-abelian fractional quantum Hall states, where we point out that the Kelvin formula is exact. PACS numbers: Jf, a, C I. INTRODUCTION A complete unerstaning of thermoelectric effects is important in the physical sciences where wie ranging applications utilize materials with large thermoelectric power S (Seebeck coefficient). Thermoelectrics of strongly correlate materials are of funamental interest since they present an important an challenging problem. Recent experiments have reveale that some materials, such as soium cobalt oxie Na x CoO 2 (NCO), possess unusually large thermopower 1, ue in part to strong electron interactions 2. Frustrate systems 3, like NCO, might prouce further surprises in enhance thermopower in some situations 2,4. In aition, emerging work 5 from the fractional quantum Hall effect (FQHE) is revitalizing thermopower as a tool to investigative the topological non-abelian quasiparticles 6 thought to exist at filling factor 5/2 7. Here we present the Kelvin formula for thermopower, S Kelvin. This is a formula inspire by Lor Kelvin s thermoynamic treatment of this variable in It is foun by reconsiering the sequence of taking the thermoynamic an uniform limits, an is a valuable approximation to the exact, but computationally intractable result, obtaine via Onsager an Kubo s treatments 9,1. For strongly correlate systems, such as the t-j moel, S Kelvin isfountopossessanaccuracybetweentherather coarse Mott-Heikes formulation, an a better argue high frequency limit formulation ue to Shastry 4 an stuie in Refs. 2 an 11. For intermeiate couplings, such as the Hubbar moel, we argue that S Kelvin provies one of the best available approximations, it is better than the high frequency limit. In certain issipation-less situations, suchas the FQHE,S Kelvin is exact, therebyproviing an elegant an simple erivation for the thermopower formula use in Ref. 5 (erive originally in Refs. 12.) S Kelvin isobtainebycompletingshastry sargument 11 for the absolute thermopower, i.e. S of an isolate system. Kelvinoriginallystuie 8 this objectusing the then available techniques, later he an others emphasize relative thermopower between two materials. Let us revert to the absolute thermopower as a starting point an imagine a long isolate cyliner of material of length L subject to a time epenent electric fiel Φ an temperature graient T. Φ couples to the ipole moment an T couples to the moment of the energy ensity (cf. Luttinger 13 ). These fiels iniviually generate a ipole moment linear in the fiels to lowest orer, an the conition for the cancellation of the two contributions, i.e., the zero ipole moment (or zero current) conition, leas to the thermopower S for a finite system size L at finite frequencies ω as S( L,ω) = Φ T { L,ω}. The thermoynamic limit, L, an the static limit, ω, must both be taken, as known from Onsager 9 an others 13,14. Kubo s exact formulas obtain in the fast or transport limit, where L before ω. Taking the static limit ω before L leas to the slow, where Kelvin s approximate formula arises an is expressible solely in terms of equilibrium thermoynamic variables. We transcribe this iscussion to a more convenient perioic system, by traing the length scale L for a wave vector q x = 2π/ L an the L limit by the uniform limit q x. The slow limit correspons to lim{q x, ω }, an the fast limit correspons to lim{ω, q x }. Thethermopowermeasurestheinuce thermoelectric voltage ue to a temperature graient an, as such, a useful an general formula for thermopower is given by the ratio between the thermoelec-

2 2 trical an electrical conuctivities 11, where χâ, ˆB(ω) = i S(q x,ω) = χ ρ(q x), ˆK( q x) (ω) T χ ρ(qx),ρ( q x)(ω), (1) = n,m te (iω + )t [Â(t), ˆB()] (2) p n p m ε m ε n +ω n  m m ˆB n (3) is the susceptibility of any two operators  an ˆB, where ρ, ˆK = Ĥ µ ˆN, an Ĵx are the charge ensity, the (gran) Hamiltonian, an the charge current operator, respectively, at finite wave vectors; Ĥ, µ, an ˆN are the Hamiltonian, the chemical potential, an the total number operator, respectively. The susceptibility written in Eq. (3) is the Lehmann representation (where p n = exp( βε n )/Z is the probability of the quantum state n with energy ε n an Z is the partition function an β = 1/k B T with k B the Boltzmann constant) which we fin useful below. With Eq. (1), we can take ifferent limits an obtain various interesting formulas. II. THERMOPOWER FORMULAE A. The Kubo formula For narrow ban systems, such as NCO, high T C superconuctors, or heavy fermion systems, the socalle Mott-Heikes (MH) approximation introuce by Heikes (popularize by Mott 15 ) is written S MH = (µ µ(t))/q e T, where µ µ(t = ). S MH is obtaine by rather rastically replacing the first part of Eq. (4) by the zero temperature chemical potential µ to make the theory sensibly behave as T. From thermoynamics, we know that µ(t) T = ( S N ) E,V, an, hence, S MH relates thermopower to the partial erivative of entropy S with particle number N, at a fixe energy E an volume V. We see below that S Kelvin is similar, but with more natural hel variables, namely T an V. C. High frequency formula From Eq. (1), we can make a high frequency approximation, where ω ω c (ω c representing all finite characteristic energy scales), leaing to the object S. The formal expression an evaluation for S are iscusse elsewhere 11 an we only quote the results. We have argue that S is the best possible approximation to the exact Kubo formula for strongly correlate systems 2 such as the t-j moel, since the high frequency limit respects the single occupancy constraint an is closer to the DC limit than initially expecte. It is not specifically suite for Hubbar type moels, since the high frequency limit assumes ω U, an cannot capture the physics of correlations effectively 11. We will see that S Kelvin steps into this breach, an provies a very useful alternative for Hubbar type moels 16. Taking the fast limit an using the continuity equations to pass from ensities to current operators, Eq. (1) gives the exact Kubo result 1 S Kubo = 1 t β τ ĴE x (t iτ)ĵx() T t µ(t) β τ Ĵx(t iτ)ĵx() q e T, (4) where q e is the charge of the carriers an ĴE the energy current. B. The Mott-Heikes formula D. The Kelvin formula To obtain an approximate thermoynamical expression, we consier the slow limit of Eq. (1). S is among the few objects (along with Hall constant an Lorentz number) where this process gives finite an approximate results, unlike the electrical conuctivity where the slow limit gives meaningless results 11. This limit is ientifie with Kelvin since he essentially took the equilibrium limit of an interacting gas of particles. The slow limit (q x, ω ) is easiest to compute starting from Eq. (1) χ ρ(qx), S Kelvin = lim ˆK( q () x) q x T χ. (5) ρ(qx),ρ( q x)() To simplify we first consier the numerator of Eq. (5) which we rewrite by first using the Lehmann representation an then taking the q x limit. Note that ˆρ(q x ) tens to a conserve quantity q e N an cannot mix states of ifferent

3 energy so ε m ε n. Thus, lim χ q ρ(q x x), ˆK( q x) () = lim q x n,m p n p m ε m ε n n ˆρ(q x ) m m ˆK( q x ) n p n p m = lim n ˆρ(q x ) m m ε m ε n ε n,m m ε ˆK( q x ) n n 1 e β(εn εm) = lim p n n ˆρ(q x ) m m ε n,m m ε ˆK( q x ) n n = q e βp n δ εn,ε m n ˆN m m ˆK n ε m ε n n,m = q e β[ ˆN ˆK ˆN ˆK ] [ = q e µ Ĥ µ ] µ ˆN. (6) The erivative with respect to µ in Eq. (6) is within the gran canonical ensemble an performe with a fixe V an T. The enominator of Eq. (5) is treate similarly yieling q 2 eβ[ ˆN 2 ˆN 2 ] = q 2 e ˆN /µ. Combining it with Eq. (6), yiels S Kelvin = 1 q e T µ Ĥ µ µ ˆN µ ˆN. (7) 3 To further simplify Eq. (7) we note a relation foun in textbooks on thermoynamics in the gran canonical ensemble: Ĥ E = Ω + TS + µn (Ω the gran potential), so that ( E µ ) T,V = µ( N µ ) T,V + T( S µ ) T,V an hence S Kelvin = 1 q e ( S µ ) T,V ( N µ ) T,V = 1 ( ) S q e N = 1 ( ) µ q e T T,V N,V (8) (9) where we use, to go from the secon equality to the last equality, a Maxwellrelationobtaine with F = ST pv + µn, an equating 2 F T N = 2 F N T. We refer to the last two equivalent equations (Eqs. (8) an (9)) as the Kelvin formula for the thermopower. This formula is unknown in the literature as far as we are aware. Note that S MH is similar to S Kelvin. The istinction is that in S Kelvin, the number erivative of the entropy is taken at constant T rather than at constant E. Thus, in the low T limit of a metal, where µ(t) T 2, they iffer in the linear T coefficient by a significant factor of two. We show below that for non-interacting electrons, scattere by impurities, S Kelvin is closer to the exact result than S MH. Further, we see that the approximation of exchanging the slow an fast limits has some justification in issipation-less systems, such as in the FQHE where S Kelvin is ientical to that foun by several workers (see below). III. APPLICATIONS OF THERMOPOWER FORMULAE A. Free electrons To gain insight into the strengths an weaknesses of the various thermopower formulations iscusse above we consier non-interacting egenerate electrons treate within the limit of elastic scattering at the Born level with an energy momentum epenent relaxation time τ(p,ω). This is a moestly issipative system, but at such a simple level that the Boltzmann-Bloch equation is an aequate escription. The solution for S is available in textbooks an a useful benchmark for various approximations. In the low temperature limit 11, to O(T 3 ), S Mott = T π2 kb 2 3q e µ ln[ ρ (µ) (vp) x 2 ] τ(p,µ) µ µ µ, (1) a formula often ascribe to Mott an ρ (µ) is the single particle ensity of states per unit volume per spin. In this non-interacting electron context, S Kelvin gives (to O(T 3 )), S Kelvin = T π2 k 2 B 3q e µ ln[ρ (µ)] µ µ, (11) which iffers from the exact answer (Eq. (1)) in the neglect of the relaxation time τ an particle velocity v x p in the logarithm. S MH, to the same orer, gives S MH = T π2 k 2 B 6q e µ ln[ρ (µ)] µ µ (12)

4 4 which is off by an important factor of two from S Kelvin (Eq. (11)). The formulations (Mott-Heikes an Kelvin) woul be ientical if µ T, which occurs if the system possesses a groun state egeneracy, an in the classical regime. The high frequency approximation gives a better result than all these an, again in the low temperature limit, to O(T 3 ), S = T π2 k 2 B 3q e µ ln[ ρ (µ) (v x p) 2 µ ] µ µ. (13) Other than the neglect of the energy erivative of τ, this is the same as the exact result. Hence, ranking the thermopower approximations for non-interacting electrons we have, from worst to best, S MH, S Kelvin, an S with the exact result being S Mott. B. Hubbar moel For intermeiate coupling moels, the relative rankings of the various approximations can be ifferent. In particular, S Kelvin can be superior to S, since the effect of correlations is ilute in the latter by making the assumption of ω U, whereas S Kelvin retains ω U. The sign of the true (i.e. transport) thermopower an the transport Hall constant are expecte to flip as we approach half filling in the Hubbar or t-j moels ue to the onset of correlations (carriers become holes measure from half filling rather than from a completely fille ban). In the case of the t-j moel, the high frequency Hall constant R H an S o isplay this behavior 11. However, for the Hubbar moel, R H an S o not isplay a sign change 16,17. S Kelvin on the other han, oes appear to showthe expectechangeinsign 16,18. Furtheriscussion concerning the relative merits of S Kelvin an S will be reporte later 16. C. NCO an the t-j moel To show the usefulness of S Kelvin, we apply it to NCO since (i) we have previously investigate 2 this system while benchmarking S, (ii) the system is intrinsically interesting 1, an (iii) we can compare ifferent thermopower formulations on equal footing. As iscusse 2, the action in NCO takes place primarily in the cobalt oxie planes where -shell spin-1/2 electrons live on a triangular lattice an these strongly interacting 2D electrons canbemoelewiththet-j moel. Hence, weexactlyiagonalizethe t-j moel on a L = 12 site two-imensional triangular lattice with perioic bounary conitions (cf. Fig.1e). Notethatweonlyshowresultsforthet-J moel with zero super-exchange interaction (J = ), as the results only weakly epen on J. To map the t-j moel to NCO we follow Refs. 2 an 4 an give results as a function of electron oping x = 1 n away from half filling (n is electron number ensity). Thermopower [µv/k] Thermopower [µv/k] (a) x =.67, t >, J = S * S Kelvin S MH (c) x =.75, t >, J = (b) x =.67, t <, J = () x =.75, t <, J = FIG. 1: Thermopower vs. T for the t-j moel (with J = ) corresponing to NCO in the Curie-Weiss metallic phase near x.7. (a) an (c) correspon to x =.67 an x =.75 for the NCO system (t > ) while in (b) an () the sign of the hopping has been switche to investigate the enhancement expecte for frustrate systems. The black, re (light gray), an blue (ark gray) lines are S, S Kelvin, an S MH. Finite size effects at low T are treate in the spirit escribe previously 2. At each x, for T below an appropriately chosen cutoff temperature T =.5 t, the thermopower is fit to S(T) at +bt 2 where a an b are obtaine from the compute S(T ) an S (T ) proviing a sensible extrapolation to low T an plotte as ashe lines. The inset figure (e) epicts the 12-site unit cell. S aequately escribes the physics of NCO for x >.5 an, in particular the so-calle Curie-Weiss metallic phase 2 near x.7. The subject of this work, however, is S Kelvin. We see in Figs. 1a an c an 2a, similar to S MH, S Kelvin oes a goo job capturing the physics with minimal computational effort. However, S Kelvin oes seem to overestimate the thermopower for intermeiate temperatures an high opings as compare to S MH. Near x.7, S Kelvin an S MH are similar but as x is ecrease the two formulae iverge an for low opings, S Kelvin better captures the physics as it is closer to the more accurate high frequency limit S. An interesting property of the triangular lattice unerlying the physics of NCO is its geometrical frustration 3, cf. inset Fig. 1e. It was preicte 2,4 that if the sign of the hopping amplitue were flippe to t < the thermopower woul be enhance at low to intermeiate T. We have consiere this situation in Figs. 1b an an 2b. Since the thermopower enhancement for t < compare to t > is largely a consequence of electronelectron interaction it is important to etermine whether this effect is capture by S Kelvin. We see this enhancement is capture to some extent by S Kelvin an S Kelvin is better than S MH in the large oping region where the enhancement is the greatest, but is missing some of the electron-electron physics at very low T that is capture by S (as is S MH ). (e) t?

5 5 (a) t >, J= (b) t <, J=.8.6 x.6 x.4.4 S * S Kelvin S MH.2 S * S Kelvin S MH FIG. 2: The thermopower vs. T an x for the t-j moel (with J = ) for t > (a) an t < (b). Note that x.7 correspons to the Curie-Weiss metallic phase of NCO, cf. Fig. 1. The line-type an color-coing is the same as in Fig. 1. Note that for opings below.5 it is not clear whether the t-j moel aequately escribes the physics of NCO. D. FQHE at ν = 5/2 We now iscuss how S Kelvin is applie to issipationless systems such as the FQHE where thermopower can beuseasapossiblenon-abelianquasiparticleetector 5. For a weakly isorere electron system (from Eq. (1) an (11)) S Kelvin essentially gives the issipation-less thermopower where particle velocities are further approximate. If the system is issipation-less an the particle velocities are also energy inepenent, such as the FQHE, then we expect S Kelvin is exact. An expression for the thermopower in a 2D electron system in the presence of a perpenicular magnetic fiel (the FQHE system) has been erive 5,12, assuming zero impurities, as S q en (Eq. (6) in Ref. 5). Yang an Halperin show5 that S k B N log() where > 1 is the quantum imension of the quasiparticles for the FQHE at ν = 5/2 (provie they are non-abelian). Thus, a non-zero entropy linear in N is obtaine. From Eq. (8), we see that the thermopower is the erivative of the entropy with respect to the number of particles at constant T an V. When entropy is linear in particle number, as in non-abelian FQHE states, S/ N S/N an the formulas are ientical. Our erivation provies a simple an straightforwar insight into the formula given previously 5. E. High Temperature Superconuctors Before concluing, we point out an intriguing application of S Kelvin for high T c superconuctors. For ifferent families of high T c compouns, a universal curve of the thermopower, at T = 29K, as a function ofhole concentrationp 1 nhasbeen observe 19. The thermopower, in all families, vanishes near optimal oping (p.16) starting out positive at small p. Phillips et al. 2 appeal to the atomic limit of S MH as an explanation. Viewing this ata 19 more generally, through the prism of S Kelvin (Eq. (8)) we conclue that the optimal filling, i.e., maximum T c, aitionally correspons to a local maximum of the electronic entropy as a function of filling. This conclusion is powerful, since we avoie the ifficult issue of calculating either thermopower or entropy, merely using the link between them provie by S Kelvin. IV. CONCLUSION It is clear that S MH, which has serve as a virtual workhorse for years, has a new competitor in S Kelvin. This simple mine approximation can be written in close form an in many ifficult regimes, where the exact Kubo-Onsagerexpressions are not useful, an S Kelvin provies an excellent guie to the physics of the system. We acknowlege support from DOE-BES Grant No. DE-FG2-6ER46319 at UCSC, an MRP acknowleges support from LPS-CMTC at University of Marylan an Microsoft Project Q. 1 I. Terasaki, Y. Sasago, an K. Uchinokura, Phys. Rev. B 56, R12685 (1997); Y. Wang et al, Nature 423, 425 (23); M. L. Foo et al, Phys. Rev. Lett. 92, 2471 (24). 2 J. O. Haerter, M. R. Peterson, an B. S. Shastry, Phys. Rev. Lett. 97, (26); M. R. Peterson, B. S. Shastry, an J. O Haerter, Phys. Rev. B. 76, (27); J. O. Haerter, M. R. Peterson, an B. S. Shastry, ibi. 74, (26). 3 A. P. Ramirez, in More is Different, eite by N. P. Ong an R. N. Bhatt (Princeton University Press, New Jersey, 21). 4 B. S. Shastry, Phys. Rev. B. 73, (26). 5 K. Yang an B. I Halperin, Phys. Rev. B. 79, (29). 6 G. Moore an N. Rea, Nucl. Phys. B 36, 362 (1991). 7 R. Willett et al., Phys. Rev. Lett. 59, 1776 (1987). 8 W. L. K Thomson, Proc. Roy. Soc. Einburgh 123 Collecte Papers I, pp , (1854). 9 L. Onsager, Phys. Rev. 37, 45 (1931); ibi. 38, 2265 (1931). 1 R. Kubo, J. Phys. Soc. Jpn. 12, 57 (1957). 11 B. S. Shastry, Rep. Prog. Phys. 72, 1651 (29). 12 Y. N. Obraztsov, Fiz. Tver. Tela (Leningra) 7, 573 (1965) [Sov. Phys.Soli State 7, 455(1965)]; N. R.Cooper,

6 6 B. I. Halperin, an I. M.Ruzin, Phys. Rev. B 55, 2344 (1997). 13 J. M. Luttinger, Phys. Rev. 135, A155 (1964); ibi. 136, A1481 (1964). 14 S. Ewars, Phil. Mag. 3, 33, 12 (1958). 15 P. M. Chaikin, an G. Beni, Phys. Rev. B 13, 647 (1976); R. R. Heikes, Thermoelectricity (Wily-Interscience, New York, 1961). 16 L. Arsenault, B. S. Shastry, an A. M. Tremblay, In preparation.. 17 F. Assa, an M. Imaa, Phys. Rev. Lett. 74, 3868 (1995). 18 T. Pruschke, D. L. Cox, an M. Jarrell, Phys. Rev. B 47, 3553 (1993). 19 S. Obertelli,, J. R. Cooper,, an J. L Tallon. Phys. Rev. B 46, (1992); T. Honma an P. H. Hor, ibi. 77, (28). 2 P. Phillips, T.-P. Choy, an R. G. Leigh, Rep. Prog. Phys. 72, 3651 (29).