Scaling analysis of negative differential thermal resistance

Transcription

1 Scaling analysis of negative differential thermal resistance Dahai He Department of Physics, Xiamen University, China Workshop on Nanoscale Heat Tehran April, 205

2 2

3 3

4 PhD Positions available! 4

5 Part I Outline. Introduction 2. Analysis 3. Examples 4. Conclusion 5

6 Negative differential electrical resistance in tunneling diode Leo Esaki (925- ) Nobel Prize in Physics(973) Illustrative figure for the work by L. Esaki, Phys. Rev. 09, 603 (958) 6

7 NDTR in thermal diode T L T R simulation analytical calculation j B. Li, L. Wang, G. Casati, PRL 93, 8430 (2004) B. Hu, DH, L.Yang, Y. Zhang, PRE 74, 0600 (2006) Negative differential thermal resistance (NDTR) 7

8 NDTR in thermal transistor B. Li, L. Wang, G. Casati, APL 88, 4350 (2006) 8

9 Mechanism of NDTR K int T L T R H L H R K 2 ( ) int 2 H = HL + x x0 + HR 2 pi 2 whre HLR, = + ( xi+ xi) + ULR, ( xi) 2 2 Difficulty of analysis : I. Nonequilibrium stationary state II. Nonlinearity III. Qusi-particle entity of phonons 9

10 Mechanism of NDTR K int T L T R H L K 2 ( ) int 2 H = HL + x x0 + HR H R 2 pi 2 whre HLR, = + ( xi+ xi) + ULR, ( xi) 2 2 j = σ T k B max where ( ) ωmin ω σ= αωdω 2π Transport process in the nonlinear response regime DH, S. Buyukdagli, B.Hu, Phys. Rev. B 80, (2009) 0

11 T + T NDTR region j (0-4 ) T Negetive differential thermal resistance Negetive in differential electrical the Frenkel- Kontorova (FK) model resistance in tunneling diode DH, B. Ai, H. Chan, and B. Hu, Phys. Rev. E 8, 043 (200) Illustrative figure for the work by L. Esaki, Phys. Rev. 09, 603 (958)

12 Size effect of NDTR N=32 N=64 N=28 N=52 0. J J = Nj T Shrinkage of the NDTR regime for increasing N: FK model DH, B. Ai, H. Chan, and B. Hu, Phys. Rev. E 8, 043 (200) 2

13 Size effect of NDTR N=32 N=64 N=28 N=256 N=52 N=024 N=2048 continuum J J = Nj 0 00 T Heat current approaches to saturation as N increases: ϕ 4 model DH, B. Ai, H. Chan, and B. Hu, Phys. Rev. E 8, 043 (200) 3

14 Q: what is the necessary conditions for the occurring of NDTR?. Spatially asymmetric structure? ( ) 2. Nonlinearity? ( ) (Φ 4 vs. FPU-β?) 3. Temperature? 4. System size? Q2: Is it possible to give a prediction of the occurring of NDTR? 4

15 A general theoretical analysis for NDTR jt ( +, T ) = jt (, T) 2 jt (, T) jt (, T) dj( T, T ) = dt + d T T T T T T = T + + T T + T T [ 0,2T ] For a particular constraint F( T, T ) = 0 Along this curve, negative differential thermal resistance (NDTR) corresponds to jt (, T) T F( T, T) = 0 < 0 H.-K. Chan, DH, B.Hu, Phys. Rev. E 89, (204) 5

16 A general theoretical analysis for NDTR NDTR: n ( T, T) n ( T, T) < [ + n ( T, T)] 2 3 where n ( T, T) n ( T, T) 2 dlnt d T = ln F( T, T) 0 ln κe( T, T) lnt a particular way of varying the temperature difference T Dependence of the effective thermal conductivity κ e on T ln κe( T, T) n3 ( T, T) ln T (, ) Nj κe T T T T Dependence of the effective thermal conductivity κ e on ΔT Example: T T dt = 0, T = + const, n ( T, T ) = [0,] 2 2T 6

17 NDTR: nt (, Tn ) ( T) < 2 3 n (T, T) NDTR 2 PDTR PDTR PDTR - n 2 (T) PDTR -2-3 NDTR Note: one can always ensure the inequality is satisfied by choosing a suitable value of n ( T, T) (, ) + H.-K. Chan, DH, B.Hu, Phys. Rev. E 89, (204) 7

18 κ e φ 4 model FPU-β model Fixing the value of T + or T - at T + t ( ) Nj κ ( ) e T CT+ t T γ 8

19 Example Φ 4 model 2 pi H = + ( x x ) + x i 2 4 i+ i i Constraint: fix T - N=32 N=64 N=28 N=256 N=52 N=024 N=2048 γ κeff ( T) = CT ( + t), where γ [0.95,.24], t 0 NDTR: nn 2< κ eff 0. γ >, T > T* _ T where: T* = 2( T + t) γ 9

20 N=32 N=64 N=28 N=256 N=52 N=024 N=2048 continuum.25 NDTR: γ >, T > T* = 2( T + t) γ J.20.5 * NDTR: N < N γ T N* The vertical lines indicates the 0.90 theoretical 0 value of ΔT* N 20

21 Example 2 Φ 4 model Constraint: 2 pi H = + ( x x ) + x i 2 4 i+ i i T+ + T m0 T = m0 m = T + T Tmax = 2 2 /2 n m T γ =, n = T + t/ T 2 m NDTR: ( γ ) m > 0, T > T* = m0 + t ( γ ) m m + t/ m 0 > 2 γ + t/ m0 2

22 m 0.50 Presence of a NDTR regime 0.25 Absence of a NDTR regime m o Phase diagram of m and m o for the ϕ 4 model 22

23 0.08 m o = j T Values of m m0 + t The vertical lines indicate the theoretical values of T* = ( γ ) m 23

24 FPU-β model Example 3 2 pi H = + ( x x ) + ( x x ) i 2 4 i+ i i+ i Constraint: n T = m0 + m T m T γ =, n = T + t/ T 2 NDTR: ( γ ) m > 0, T > T* = m0 + t ( γ ) m m t, m + t/ m 0 0 > < γ 2 γ + t/ m0 24

25 (b) 00 m Absence of a NDTR regime Absence of a NDTR regime Possible presence of a NDTR regime t / γ m o Phase diagram of m and m o for the FPU-β model 25

26 m o = j T* T* T Values of m m Note that 0 T*/ Tmax = Tmax = γ /2 m NDTR regime is too narrow to be observed! 26

27 FPU-β model A very specific example Consider a change of temperatures of heat baths ( T, T ) = (5, 49) ( T, T ) = (4.005, ) + + T = 2 T = 2.0 Theory: n T dt = 44.56, nn 5.35 < T d T 2 Numerics: j = j =

28 Conclusion We develop a system-independent scaling analysis, and obtain the general condition for the occurrence of NDTR. Based on the condition, one can judge whether NDTR exist; If NDTR exist, ΔT* and N* can be predicted. The occurrence of NDTR can be manipulated for any nonlinear model by suitably choosing the way of varying T + and T -. 28

29 Thank You! 29

30 Part II Thermal expansion and its impacts on thermal transport in the FPU-α-β model X. Cao, DH, H. Zhao, and B. Hu, AIP Advances 5, (205) 30

31 Motivation I Recent controversy on the effect of asymmetric interaction potential on normal thermal conduction. (Hong Zhao s and Shunda Chen s talk) 3

32 Motivation II: application aspects With the rapid development of nanotechnique, thermal expansion plays an important role for thermal measurement, designing nanodevices with intriguing electronic, mechanical and thermal properties. Thermal Expansion Thermometry Appl. Phys. Lett. 66,

33 Motivation III: theoretical aspects Most of previous analytical studies used the perturbation approach, such as lattice-dynamics calculations, and nonequilibrium Green s function theory, which is incapable of dealing with strong anharmonicity for which some concerned intriguing properties occur. 33

34 Potential Profile of FPU-ab model 60 0 V(x) -60 α=0 α= α= 2 α= 3 α= 4 α= 5 α= 6 x

35 Quantify the asymmetry 20 0 (a) 0 - (b) V(x) -40 Σ Σ= S -S 2 /(S +S 2 ) x S x S x x c α 35

36 Temperature profile T.9 α= α= α= 2 α= 3 α= i 36

37 Thermal conductance 0.5 G γ ~ 0 γ < 0 γ > α The nonmonotonic behavior of G can be divided by three domains, corresponding to negative, positive and vanishing coefficient of thermal expansion γ, respectively. X. Cao, DH, H. Zhao, and B. Hu, AIP Advances 5, (205) 37

38 Self-consistent phonon theory (SCPT) Incorporating the nonlinearity into normal modes by renormalizing the harmonic frequency spectrum, which is realized by performing thermal average with respect to a trial Hamiltonian 2 eff pi f H = + ( ui+ ui) 2m 2 Where the effective harmonic potential coefficient f(t) can be obtained from the self-consistent equations: 2 V x x 2 ( ) V( x) = 0, = x f T. Dauxois, et al, Phys. Rev. E 47, 684(993) DH, S. Buyukdagli, and B. Hu, Phys. Rev. E 78, 0603 (2008) 38

39 Coefficient of thermal expansion 0 γ ~ 0 γ < 0 l γ > 0 0. α=- α=-3 α= T 39

40 Effect of nonlinearity on thermal expansion 6 5 MD SCPT 4 3 l 2 0 α X. Cao, DH, H. Zhao, and B. Hu, AIP Advances 5, (205) 40

41 Conclusion Three domains of thermal conductance with respect to α are identified, which is related to thermal expansion effect. Self-consistent phonon theory is developed to study the effect of thermal expansion, which agrees well with the numerical simulations. 4

42 Thank You! 42