Thermoelectricity at nanoscale: theoretical models Dima Shepelyansky (CNRS, Toulouse) www.quantware.ups-tlse.fr/dima Joint work with Oleg Zhirov (BINP) Europhys. Lett. 03, 68008 (203) 2 3 3 T/Kc 4 2 0 0 2 3 4 K/K c Time evolution of figure of merit ZT (center) (from A.Majumdar Science 303, 777 (2004)) (Quantware group, CNRS, Toulouse) Inst Semicond Phys RAS Dec 0, 204 / 22
Early works from Y.Imry (Weizmann Inst) talk at Inst. H. Poincaré, Paris (202) (Quantware group, CNRS, Toulouse) Inst Semicond Phys RAS Dec 0, 204 2 / 22
Early works from Y.Imry (Weizmann Inst) talk at Inst. H. Poincaré, Paris (202) (Quantware group, CNRS, Toulouse) Inst Semicond Phys RAS Dec 0, 204 3 / 22
Main characteristics Seebeck coefficient: S = V/ T = π 2 k B 2 T[(d lnσ/de)] EF /e (Mott relation (958)) For 2DEG with Wiedemann-Franz law: S = 2πk B 2 Tm/(3eh 2 n e ); typical value S 0µV/K at T = 0.3K, n e = 4 0 0 cm 2 Thermoelectric figure of merit ZT = σs 2 T/κ, thermal conductivity κ = κ el +κ phonon (heat flux Q = κ T) (Quantware group, CNRS, Toulouse) Inst Semicond Phys RAS Dec 0, 204 4 / 22
Experiments on Seebeck coefficient for 2DEG R.Fletcher, V.M.Pudalov et al. Semicond. Sci. Technol. 5, 386 (200) (Quantware group, CNRS, Toulouse) Inst Semicond Phys RAS Dec 0, 204 5 / 22
Experiments on Seebeck coefficient for 2DEG V.Narayan, S.Goswami, M.Pepper et al. PRB 85, 25406 (202) In dimensionless units S = 0mV/K 00 (Quantware group, CNRS, Toulouse) Inst Semicond Phys RAS Dec 0, 204 6 / 22
ZT in various materials from http://www.thermoelectrics.caltech.edu/thermoelectrics/ (Quantware group, CNRS, Toulouse) Inst Semicond Phys RAS Dec 0, 204 7 / 22
ZT in p-type Bi 2 Te 3 /Sb 2 Te 3 superlattices R.Venkatasubramanian et al. (N.Carolina) Nature 43, 597 (200) (Quantware group, CNRS, Toulouse) Inst Semicond Phys RAS Dec 0, 204 8 / 22
ZT in SnSe Li-Dong Zhao et al. (Illinois) Nature 508, 373 (204) (Quantware group, CNRS, Toulouse) Inst Semicond Phys RAS Dec 0, 204 9 / 22
Research at CNRS Grenoble (O.Bourgeois Refs.3,4) (Quantware group, CNRS, Toulouse) Inst Semicond Phys RAS Dec 0, 204 0 / 22
Applications (Quantware group, CNRS, Toulouse) Inst Semicond Phys RAS Dec 0, 204 / 22
Thermoelectricity of Wigner crystal in a periodic potential Hamiltonian H = i ( p 2 i 2 + K cos x i + 2 j i ) x i x j Dynamic equations ṗ i = H/ x i + E dc ηp i + gξ i (t), ẋ i = p i Here the Langevin force is given by g = 2ηT, ξ i (t)ξ j (t ) = δ ij δ(t t ); n e = ν/2π, ν = ν g =.68... Fibonacci rational approximates. Aubry transition at K = K c = 0.0462. I.Garcia-Mata, O.Zhirov, DLS EPJD 4, 325 (2007) (Quantware group, CNRS, Toulouse) Inst Semicond Phys RAS Dec 0, 204 2 / 22
Time evolution of Wigner crystal 2 0 5 t 0 5 0 8 9 0 2 3 x/2π Electron density variation in space and time from one Langevin trajectory at K/K c = 2.6, T/K c = 0., η = 0.02, N = 34, M = L/2π = 2; density changes from zero (dark blue) to maximal density (dark red); only a fragment of x space is shown. (Quantware group, CNRS, Toulouse) Inst Semicond Phys RAS Dec 0, 204 3 / 22
Numerical fits. T e (x)/ T.0 0.9.64.62.6 ν(x) 0 7 4 2 x/2π 2.0.8.6.4.2 ν(x) 0 7 4 2 x/2π Left panels: dependence of electron temperature Te(x) (top, blue points) and rescaled density ν(x) (bottom, black points) on distance x along the chain placed on the Langevin substrate with a constant temperature gradient (it is shown by the blue line) at average temperature T = 0.0 and temperature difference T = 0.2 T; black line shows the fit of density variation in the bulk part of the sample. Right panel: density variation produced by a static electric field E dc = 4 0 4 at a constant substrate temperature T = 0.0; black line shows the fit of gradient in the bulk part of the sample. Here N = 34, M = 2, K =.52Kc, η = 0.02, averaging is done over time interval t = 0 7 ; S = 3.3 at T = 0.0 0.22Kc. (Quantware group, CNRS, Toulouse) Inst Semicond Phys RAS Dec 0, 204 4 / 22
Seebeck coefficient S 0 2 0 0 2 4 6 K/K c S 0 2 0 K = 0 0 T/K c Left panel: Dependence of the Seebeck coefficient S on rescaled potential amplitude K/Kc at temperatures T/Kc = 0.065, 0., 0.22 and 0.65 shown by black, blue, green and red colors, respectively from top to bottom. The full and open symbols correspond respectively to chains with N = 34, M = 2 and N = 55, M = 34. Right panel: Dependence of S on T/Kc at different K/Kc = 0, 0.75,.5, 2.2, 3 shown respectively by black, violet, blue, green and red points; N = 34, M = 2; the dashed gray line shows the case K = 0 for noninteracting particles. The stars show corresponding results from left plane at same N, M. Dotted curves are drown to adapt an eye. Here and in other Figs. the statistical error bars are shown when they are larger than the symbol size. Here η = 0.02. (Quantware group, CNRS, Toulouse) Inst Semicond Phys RAS Dec 0, 204 5 / 22
Conductivity and thermal conductivity σ σ 0 0 2 0 κ κ0 0 0 4 0 2 4 6 K/K c 0 2 0 2 4 6 K/K c Left panel: Rescaled electron conductivity σ/σ 0 as a function of K/K c shown at rescaled temperatures T/K c = 0.065, 0.22, 0.65 by black, green and red points respectively. Right panel: Rescaled thermal conductivity κ/κ 0 as a function of K/K c shown at same temperatures and colors as in left panel. Here we have N = 34, M = 2, η = 0.02, σ 0 = ν g /(2πη), κ 0 = σ 0 K c. (Quantware group, CNRS, Toulouse) Inst Semicond Phys RAS Dec 0, 204 6 / 22
ZT dependence on parameters.5 ZT.0 0.5 0.0.5.0 0.5 0.0.......... ZT.......................... 0 2 4 6 K/K c 2 ZT 0 3 2 0... ZT....... 0 0 T/K c Left panels: Dependence of ZT on K/K c at temperatures T/K c = 0. (top panel) and T/K c = 0.65 (bottom panel); the black points and open triangles correspond respectively to η = 0.02 and η = 0.05 at N = 34, M = 2. Right panels: Dependence of ZT on T/K c for K/K c = 0.75 at η = 0.02, N = 34, M = 2. Bottom right panel: Same as in top right panel at K/K c = 2.6 and N = 34, M = 2 (black points); N = 89, M = 55 (green circles); N = 44, M = 89 (red stars). (Quantware group, CNRS, Toulouse) Inst Semicond Phys RAS Dec 0, 204 7 / 22
ZT dependence on parameters 2 3 3 T/Kc 4 2 0 0 2 3 4 K/K c Dependence of ZT on K/K c and T/K c shown by color changing from ZT = 0 (black) to maximal ZT = 4.5 (light rose); contour curves show values ZT =, 2, 3, 4. Here η = 0.02, N = 34, M = 2. (Quantware group, CNRS, Toulouse) Inst Semicond Phys RAS Dec 0, 204 8 / 22
Other quantities dependences 0 0 κ/κ 0... 0 0 T/K c 2 0 S 2 σ/σ 0.......... 0 0 T/K c Left panel: Rescaled thermal conductivity κ/κ 0 as a function of rescaler temperature T/K c, to adapt an eye the straight dashed line shows the dependence κ/κ 0 = 0.6T/K c ; right panel: same as in left panel for S 2 σ/σ 0. Data are obtained at K/K c = 2.6, η = 0.02, N = 34, M = 2, σ 0 = ν g /(2πη), κ 0 = σ 0 K c. (Quantware group, CNRS, Toulouse) Inst Semicond Phys RAS Dec 0, 204 9 / 22
Physical parameters In physical units we can estimate the critical potential amplitude as U c = K c e 2 /(ǫd), where ǫ is a dielectric constant, x is a lattice period and d = ν x/2π is a rescaled lattice constant Ref.5. For values typical for a charge density wave regime Ref.6 we have ǫ 0, ν, x nm and U c 40mV 500K so that the Aubry pinned phase should be visible at room temperature. The obtained U c value is rather high that justifies the fact that we investigated thermoelectricity in the frame of classical mechanics of interacting electrons. (Quantware group, CNRS, Toulouse) Inst Semicond Phys RAS Dec 0, 204 20 / 22
Message from Novosibirsk GES (Quantware group, CNRS, Toulouse) Inst Semicond Phys RAS Dec 0, 204 2 / 22
References: R. A.F.Ioffe and L.S.Syil bans, Re. Prog. Phys. 22, 67 (969) R2. H.J. Goldsmid, Introduction to thermoelectricity, Springer, Berlin (2009). R3. Nanoelectronics : Concepts, Theory and Modeling. Network meeting and workshop on thermoelectric transport 2-27 October 202, Cargese, Corsica http://iramis.cea.fr/meetings/nanoctm/program.php R4. Slides at R3 by O.Bourgeois, L.W.Molenkamp, S.Voltz R5. I.Garcia-Mata, O.V.Zhirov, D.L.Shepelyansky, EPJD 4, 325 (2007) R6. S.Brazovskii et al., Phys. Rev. Lett. 08, 09680 (202) (Quantware group, CNRS, Toulouse) Inst Semicond Phys RAS Dec 0, 204 22 / 22