A Series Circuit of Thermal Rectifiers: An Effective Way to Enhance Rectification Ratio

Transcription

1 Thermal Rectification A Series Circuit of Thermal Rectifiers: An Effective Way to Enhance Rectification Ratio Shiqian Hu, Meng An, Nuo Yang,* and Baowen Li* Phonon, the quantized mode of vibration, is responsible for heat transport in dielectric and semiconducting crystals. The understanding of phonon properties provides a theoretical foundation for thermal devices, such as thermal rectifiers, [1 13] optomechanical crystals, [14,15] thermal cloaking, [16 20] thermoelectrics, [21 25] thermocrystals, [26 28] and phonon nanocapacitor, [29] etc., which could lead to a new technological revolution. Graphene, a 2D structure of a single atomic layer of carbon arranged in a honeycomb lattice, [30] exhibits many intriguing properties. In addition to its unique electronic and optical properties, graphene is a promising material for thermal management of future nanoelectronic devices due to its superior thermal conductivity at room temperature. [5,31 33] In addition, due to the outstanding phonon confinement [34,35] and unique phonon coherent interference effect, [26,36,37] the porous graphene has shown potential applications in many fields, such as gas separation and purification, [38] water desalination, [39] ion channel, [40] and DNA sequencing. [41] Similar to electronic counterpart, the thermal rectifier (TR) also plays a vital role in phononics/thermal circuits. [1] The thermal rectification means that the magnitude of heat flux changes when the temperature gradient is reversed in direction. Massive efforts have been devoted to investigate a high-performance thermal rectifier. Up to now, most theoretical studies and experimental tests are realized based on S. Hu Center for Phononics and Thermal Energy Science School of Physics Science and Engineering Tongji University Shanghai , P. R. China M. An, Prof. N. Yang State Key Laboratory of Coal Combustion Huazhong University of Science and Technology (HUST) Wuhan , P. R. China nuo@hust.edu.cn M. An, Prof. N. Yang Nano Interface Center for Energy (NICE) School of Energy and Power Engineering Huazhong University of Science and Technology (HUST) Wuhan , P. R. China Prof. B. Li Department of Mechanical Engineering University of Colorado Boulder, CO 80309, USA Baowen.Li@Colorado.Edu DOI: /smll nonlinear lattice models, [6] graded mass density, [2,7] asymmetric structure, [3 5,9,42] molecular junction, [10] or an interface between two different materials. [11 13,43] Very recently, based on phase change materials, [44,45] the thermal rectification can be enhanced significantly. However, available reports to enhance the thermal rectification ratio only focus on optimizing the performance of single thermal rectifier. In this letter, a series circuit of TRs is proposed to enhance the thermal rectification. Based on our previous work, the different temperature dependence of thermal conductivity in graphene and graphene phononic crystal (GPnC), [46] we construct an asymmetric graphene/gpnc structure (Figure 1a) and found its thermal rectification. Furthermore, inspired by the series effect in electronic circuits, the graphene/gpnc 1 / GPnC 2 structure (Figure 2a) is used to investigate the series effect of thermal rectification. The dependence of the thermal rectification on the geometric size of structure is also discussed in details. A schematic picture of the asymmetric graphene/gpnc structure is shown in Figure 1a, which is comprised of graphene and GPnC with the same length. The size of pores in GPnC is characterized by L 1. The neck width (L n ) is fixed as 0.71 nm in all structures studied here. In our simulations, the lattice constant and thickness of graphene are used as and nm, respectively. The fixed (periodic) boundary condition is adopted along longitudinal (transversal) direction. In nonequilibrium molecular dynamics simulations, the three layers atoms close to the two ends are coupled with the Langevin heat bath [47] at temperature T L = T 0 (1 + Δ) and T R = T 0 (1 Δ), respectively, where T 0 is the average temperature and Δ is the normalized temperature difference. Therefore, the atoms in the left end are at a higher (lower) temperature in the case of Δ > 0 (Δ < 0). The bonding interaction between carbon atoms is described by the energy potential of a Morse bond and a harmonic cosine angle, which include both two-body and threebody potential terms. [48] We use the velocity Verlet algorithm to integrate the discretized differential equations of motions, where the time step is set as 0.5 fs. The heat flux transferred across the system is defined as the energy transported along it in unit time. [3] The structure is relaxed long enough (3 ns) such that the system reaches the steady state. Then, both the temperature and heat flux are recorded, which are averaged over 4.5 ns. It is noted that the values of temperature difference (ΔT = T L T R ) are consistent when the heat flux (J+/J ) used for calculating the thermal rectification ratio is 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim (1 of 5)

2 communications Figure 1. a) Schematic picture of asymmetric graphene/gpnc structure. The length (width) of the GPnC is set as L (W). L1 (size of pores) and Ln (neck width) are used to characterize the GPnC structure. Ln is fixed as 0.71 nm. b) The temperature profiles correspond to the heat flux along the direction from GPnC to graphene and the inverse direction in the asymmetric graphene/gpnc structure. c) The thermal rectification ratio as a function of the average temperature T0 with Δ = 0.5. d) The dependence of thermal rectification ratio on ΔT. The corresponding parameters are L = 50 nm, W = 7 nm, L1 = 0.71 nm, T0 = 1500 K. recorded. All the results presented here are averaged from 12 independent simulations from different initial conditions. The error bar is the standard errors of the results of these simulations. To show the thermal rectification effect quantitatively, the thermal rectification ratio (Rect) is defined as Rect = r r+ temperature profiles in the direction from GPnC to graphene and the inverse direction in asymmetric graphene/gpnc structure are both plotted in Figure 1b. We have checked the (1) where r+ (r ) is the thermal resistance of the system when Δ < 0 (> 0). We define the generalized heat flux transmission function (τr) as[49] τ r± = J± k B T (2) where J+ (J-) is the heat current from the right to the left (the left to the right), corresponding to Δ < 0 (> 0), kb is the Boltzmann constant, and ΔT is the temperature difference. In the ballistic limit, the transmission function is simply equal to 3N (the number of propagating phonon modes). According the definition of thermal resistance (r = ΔT/J), we can easily obtain the Rect as Rect = r = τ r+ r+ τ r (3) First, the thermal rectification effect is studied in an asymmetric graphene/gpnc structure (Figure 1a). The (2 of 5) Figure 2. Schematic pictures of a) graphene/gpnc1/gpnc2 structure, b) asymmetric graphene/gpnc structure, and c) GPnC1/GPnC2 structure. L = Lgraphene + LGPnC, and PL = Lgraphene/LGPnC. For all the three structures, the widths (W) are set as 7 nm. The fixed (periodic) boundary condition is applied along the transverse (longitudinal) direction. The sizes of pores in the middle (right) region are characterized by L1 (L2) Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim

3 overall size effect of rectification ratio in the graphene/gpnc structure. Although there are obvious size dependences of thermal conductivity of graphene [35,50] and GPnC, [46] the results show that the rectification ratio is insensitive to the overall length of the structure (details shown in Figure S1 in the Supporting Information). The thermal rectification ratio is shown in Figure 1c as a function of average temperature (T 0 ). It is found that the heat flux from GPnC region to graphene region is higher than that from graphene region to GPnC region when the temperature difference is exchanged in direction and the thermal rectification ratio is insensitive to T 0. Therefore, we choose T 0 as 1500 K in the following studies. The dependence of thermal rectification ratio on the temperature difference (ΔT) at 1500 K is shown in Figure 1d. It is observed that the thermal rectification ratio increases when ΔT increases from 300 to 2700 K. In the following, we give a discussion on the mechanism of rectification in asymmetric graphene/gpnc structure, based on the difference of thermal resistance between GPnC and graphene. Compared with the graphene, the thermal resistance of GPnC has a less sensitive dependence on the temperature. [46] The difference between the total thermal resistances of the two directions originates mainly from the resistance of graphene region. When the heat transfers from GPnC to graphene, the graphene region has a lower temperature ( 800 K, Figure 1b), corresponding to a lower thermal resistance. For the inverse direction, from graphene to GPnC, the graphene region has a higher temperature ( 2000 K, Figure 1b), corresponding to a higher thermal resistance. The smaller thermal resistance from GPnC to graphene means the heat flux prefers from GPnC to graphene. Further, inspired by the series effect in electronic rectifiers, the basic element can be composed into complex structure to enhance the rectification ratio. We proposed a three section structure, graphene/gpnc 1 /GPnC 2 (shown in Figure 2a), where the size of pores in the GPnC 1 (GPnC 2 ) region is characterized by L 1 (L 2 ) and fixed as 0.71 nm (1.55 nm). The GPnC 1 and GPnC 2 have the same length. In the three section structure, the length of the graphene and GPnC are set as L graphene and L GPnC ( L GPnC 1 + L GPnC 2), respectively. The results of thermal rectification ratio of graphene/ GPnC 1 /GPnC 2 versus ΔT are shown in Figure 3a (triangles). Similar to asymmetric graphene/gpnc, the heat flux from GPnC 2 /GPnC 1 region to graphene region is also higher and the value of thermal rectification ratio increases with the ΔT increase from 300 to 2700 K. Meanwhile, the interesting finding is that the thermal rectification ratio of graphene/ GPnC 1 /GPnC 2 is enhanced significantly compared with that of asymmetric graphene/gpnc. For a comparison with graphene/gpnc 1 /GPnC 2, the structure of graphene/gpnc 2 / GPnC 1 is also studied. The thermal rectification ratio results versus ΔT are shown in Figure 3b (triangles). The surprising finding is that, instead of the enhancement, the value of R ect is reduced. To interpret this enhancement (reduction) of the thermal rectification ratio by the series effect in three section structures, the thermal rectification ratio of the graphene/gpnc 1, TR 1, (shown in Figure 2b) and the GPnC 1 /GPnC 2, TR 2, Figure 3. Rectification ratio versus ΔT in the a) graphene/gpnc 1 /GPnC 2 structure and b) graphene/gpnc 2 /GPnC 1 structure corresponding to upright triangles. Inverse triangles are series effect prediction through Equation (4). Squares (circles) are the rectification ratio in the (a) graphene/gpnc 1 (GPnC 1 /GPnC 2 ) and (b) graphene/gpnc 2 (GPnC 2 / GPnC 1 ) structure, in which the values of ΔT are kept the same as the corresponding part in graphene/gpnc 1 /GPnC 2 and graphene/ GPnC 2 /GPnC 1 structure. The parameters are L = 50 nm, W = 7 nm, L n = 0.71 nm, P L = 1, T 0 = 1500 K, L 1 = 0.71 nm, and L 2 = 1.55 nm. The error bar is standard errors of 12 MD simulations with different initial conditions. (shown in Figure 2c) are studied. It is found that in Figure 1c, the thermal rectification ratio is not sensitive to T 0. However, it depends on ΔT (shown in Figure 1d). To simplify the analysis, we assumed the ratio as a function of ΔT and independent on T 0. In order to show the series effect and comparison directly, the same values of ΔT are used in simulation for graphene/gpnc 1 or GPnC 1 /GPnC 2 as the corresponding part in graphene/gpnc 1 /GPnC 2. The thermal rectification ratio results of graphene/gpnc 1 and GPnC 1 /GPnC 2 structure are shown in Figure 3 (circles and squares). For the GPnC 1 /GPnC 2 structure, it is found that the heat flux prefers from GPnC 2 to GPnC 1. When three section structures are arranged as graphene/gpnc 1 / GPnC 2 (graphene/gpnc 2 /GPnC 1 ), these two TRs are linked on the same (inverse) thermal rectification direction, which is responsible for the thermal rectification enhancement (reduction). For the series of two thermal rectifiers, TR 1 and TR 2, we make a very rough but intuitive assumption that they are independent and not interact with each other. Combined with the Equation (3), we can easily arrive at the constitutive thermal rectification ratio of the series rectifiers τ τ se = r12 = r1, + r2, + Rect τ τ = R1 2 ectrect r (4) 12+ where ect 1 ( ect 2 ) r1, r2, R R is the rectification ratio of TR 1 (TR 2 ). There are two sets of results for three section structures. The first one from series effect prediction by Equation (4) (inverse triangles in Figure 3) based on MD simulation of graphene/gpnc 1 and GPnC 1 /GPnC 2. And the other one from MD simulation of three section structures (triangles in 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim (3 of 5)

4 communications Figure 3). The two sets of results agree well with each other as shown in Figure 3. It demonstrates that, through a subtle arrangement of GPnCs, the thermal rectification ratio can be enhanced with series effect strategy. Furthermore, in order to verify our results, we calculated another structure of graphene/gpnc 1 /GPnC 2 whose L 2 is set as 2.41 nm. The results show that the thermal rectification ratio of graphene/gpnc 1 / GPnC 2 is also enhanced by series effect. And the MD results and theoretical predictions match well (details are shown in Figure S2 in the Supporting Information). The series effect could be extended to a graded GPnC structure whose pore sizes are continuously tapered. The performance may be maximized in the graded structure if the series effect works. Although there are a few studies on the rectification effect on graded structure, [2,4,5] it is still not clear on its relationship with the series effect which deserves to explore further. Moreover, the series effect could work in other rectifier systems [10 13] and would improve the performance of rectifier. However, considering the relative high operating temperature difference in our system, a system with applicable temperature difference is needed to develop in future. Besides the series effect, the size effect should be another important factor of nanoscale devices. The ratio of L graphene /L GPnC is defined as P L. Then, the effect of length ratio P L in asymmetric graphene/gpnc on the rectification is investigated. When P L is changed from 1/3 to 3, it is found that the thermal rectification ratio increases with P L increase (shown in Figure 4a). The value of P L is restricted by the overall length (L) and the length of GPnC region (L GPnC ). The hardware limits the overall length of system simulated. It should be pointed out that the value of P L can be further increased by increasing L. As mentioned above, the difference of the total thermal resistances come mainly from the resistance of graphene region. For graphene, the thermal conductivity (κ G ) increases with its length (L), as κ G ~ logl, which has been shown in recent experiment study. [50] As the length ratio (P L ) increases, due to size effect, the graphene region will have a larger thermal conductivity and smaller resistance. As shown in Figure 4b, the heat flux increase as P L increases, which results in a larger difference of resistance and thermal rectification ratio. As both the series effect and size effect are all effective methods to enhance the thermal rectification ratio, we further calculate the thermal rectification of graphene/gpnc 1 / GPnC 2 structures, whose P L s are set as 3 (the optimized value in size effect). The results show that the combination of size effect and series effect have more significant effect on thermal rectification ratio than anyone of those alone (details are shown in Figure S3 in the Supporting Information). In summary, using nonequilibrium molecular dynamics simulations, we studied thermal rectification in both two section and three section asymmetric graphene/gpnc structures (Figure 2). The series effect is demonstrated by the consistence between the results of theoretical prediction and that of MD simulations (Figure 3). It is found that both the series effect and size effect are effective strategies to enhance the thermal rectification ratio. The series model advances previous thermal rectifier designs. More recently, the large scale synthesis [51] of porous graphene sheet with pore size ranging from atomic precision (<1 nm) to nanoscaled dimension ( nm) has been demonstrated experimentally, therefore, the high-performance thermal rectifier proposed in this paper could be realized experimentally in a foreseeable future. Supporting Information Supporting Information is available from the Wiley Online Library or from the author. Acknowledgements This project was supported in part by the grants from the National Natural Science Foundation of China: and The authors are grateful to Junichiro Shiomi for useful discussions. The authors thank the National Supercomputing Center in Tianjin (NSCC-TJ) for providing assistance in computations. Figure 4. a) Rectification ratio versus ΔT in the asymmetric graphene/ GPnC structures. The value of thermal rectification ratio increases with the increase of P L. b) Heat flux along the direction from GPnC to graphene and the inverse direction in asymmetric graphene/gpnc structure versus the length ratio (P L ). It shows that heat flux (thermal resistance) increases (decreases) with P L increase. The ΔT increases from 300 to 2700 K. P L is in the range from 1/3 to 3. The structure parameters are L = 50 nm, W = 7 nm, L n = 0.71 nm, L 1 = 0.71 nm, and T 0 = 1500 K. The error bar is standard errors of 12 MD simulations with different initial conditions. [1] N. Li, J. Ren, L. Wang, G. Zhang, P. Hanggi, B. Li, Rev. Mod. Phys. 2012, 84, [2] N. Yang, N. Li, L. Wang, B. Li, Phys. Rev. B 2007, 76, [3] N. Yang, G. Zhang, B. Li, Appl. Phys. Lett. 2008, 93, [4] N. Yang, G. Zhang, B. Li, Appl. Phys. Lett. 2009, 95, [5] J. Hu, X. Ruan, Y. P. Chen, Nano Lett. 2009, 9, [6] B. Li, L. Wang, G. Casati, Phys. Rev. Lett. 2004, 93, [7] C. W. Chang, D. Okawa, A. Majumdar, A. Zettl, Science 2006, 314, (4 of 5) Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim

5 [8] Y.-Y. Liu, W.-X. Zhou, L.-M. Tang, K.-Q. Chen, Appl. Phys. Lett. 2014, 105, [9] Y. Y. Liu, W. X. Zhou, K. Q. Chen, Sci. Rep. 2015, 5, [10] D. M. Leitner, J. Phys. Chem. B 2013, 117, [11] W. Kobayashi, Y. Teraoka, I. Terasaki, Appl. Phys. Lett. 2009, 95, [12] H. Acharya, N. J. Mozdzierz, P. Keblinski, S. Garde, Ind. Eng. Chem. Res. 2011, 51, [13] N. A. Roberts, D. G. Walker, J. Heat Transf. 2011, 133, [14] M. Eichenfield, J. Chan, R. M. Camacho, K. J. Vahala, O. Painter, Nature 2009, 462, 78. [15] A. H. Safavi-Naeini, J. T. Hill, S. Meenehan, J. Chan, S. Groblacher, O. Painter, Phys. Rev. Lett. 2014, 112, [16] T. Han, X. Bai, J. T. Thong, B. Li, C. W. Qiu, Adv. Mater. 2014, 26, [17] S. Narayana, S. Savo, Y. Sato, Appl. Phys. Lett. 2013, 102, [18] S. Guenneau, C. Amra, D. Veynante, Opt. Express 2012, 20, [19] T. Han, X. Bai, D. Gao, J. T. L. Thong, B. Li, C. W. Qiu, Phys. Rev. Lett. 2014, 112, [20] T. Han, T. Yuan, B. Li, C. W. Qiu, Sci. Rep. 2013, 3, [21] G. Chen, M. S. Dresselhaus, G. Dresselhaus, J. P. Fleurial, T. Caillat, Int. Mater. Rev. 2003, 48, 45. [22] M. S. Dresselhaus, G. Chen, M. Y. Tang, R. G. Yang, H. Lee, D. Z. Wang, Z. F. Ren, J. P. Fleurial, P. Gogna, Adv. Mater. 2007, 19, [23] C. J. Vineis, A. Shakouri, A. Majumdar, M. G. Kanatzidis, Adv. Mater. 2010, 22, [24] Z. Jin, Q. Liao, H. Fang, Z. Liu, W. Liu, Z. Ding, T. Luo, N. Yang, Sci. Rep. 2015, 5, [25] L. Yang, N. Yang, B. Li, Int. J. Heat Mass Transfer 2016, 99, 102. [26] M. Maldovan, Phys. Rev. Lett. 2013, 110, [27] M. Maldovan, Nature 2013, 503, 209. [28] J. K. Yu, S. Mitrovic, D. Tham, J. Varghese, J. R. Heath, Nat. Nanotechnol. 2010, 5, 718. [29] H. Han, B. Li, S. Volz, Y. A. Kosevich, Phys. Rev. Lett. 2015, 114, [30] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, M. I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, A. A. Firsov, Nature 2005, 438, 197. [31] D. L. Nika, S. Ghosh, E. P. Pokatilov, A. A. Balandin, Appl. Phys. Lett. 2009, 94, [32] A. A. Balandin, S. Ghosh, W. Bao, I. Calizo, D. Teweldebrhan, F. Miao, C. N. Lau, Nano Lett. 2008, 8, 902. [33] Z. Wang, R. Xie, C. T. Bui, D. Liu, X. Ni, B. Li, J. T. Thong, Nano Lett. 2011, 11, 113. [34] L. Yang, N. Yang, B. Li, Nano Lett. 2014, 14, [35] L. N. Yang, J. Chen, N. Yang, B. Li, Int. J. Heat Mass Transfer 2015, 91, 428. [36] S. Alaie, D. F. Goettler, M. Su, Z. C. Leseman, C. M. Reinke, I. El-Kady, Nat. Commun. 2015, 6, [37] N. Zen, T. A. Puurtinen, T. J. Isotalo, S. Chaudhuri, I. J. Maasilta, Nat. Commun. 2014, 5, [38] A. Du, Z. Zhu, S. C. Smith, J. Am. Chem. Soc. 2010, 132, [39] D. Cohen-Tanugi, J. C. Grossman, Nano Lett. 2012, 12, [40] K. Sint, B. Wang, P. Kral, J. Am. Chem. Soc. 2008, 130, [41] G. F. Schneider, S. W. Kowalczyk, V. E. Calado, G. Pandraud, H. W. Zandbergen, L. M. Vandersypen, C. Dekker, Nano Lett. 2010, 10, [42] G. Wu, B. Li, Phys. Rev. B 2007, 76, [43] M. Hu, J. V. Goicochea, B. Michel, D. Poulikakos, Appl. Phys. Lett. 2009, 95, [44] R. Chen, Y. Cui, H. Tian, R. Yao, Z. Liu, Y. Shu, C. Li, Y. Yang, T. Ren, G. Zhang, R. Zou, Sci. Rep. 2015, 5, [45] T. Zhang, T. F. Luo, Small 2015, 11, [46] S. Hu, M. An, N. Yang, B. Li, Nanotechnology 2016, 27, [47] N. Yang, G. Zhang, B. Li, Nano Lett. 2008, 8, 276. [48] R. E. Tuzun, D. W. Noid, B. G. Sumpter, R. C. Merkle, Nanotechnology 1996, 7, 241. [49] K. Sääskilahti, J. Oksanen, S. Volz, J. Tulkki, Phys. Rev. B 2015, 91, [50] X. Xu, L. F. Pereira, Y. Wang, J. Wu, K. Zhang, X. Zhao, S. Bae, C. Tinh Bui, R. Xie, J. T. Thong, B. H. Hong, K. P. Loh, D. Donadio, B. Li, B. Ozyilmaz, Nat. Commun. 2014, 5, [51] D. Zhou, Y. Cui, P. W. Xiao, M. Y. Jiang, B. H. Han, Nat. Commun. 2014, 5, Received: August 17, 2016 Revised: October 14, 2016 Published online: 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim (5 of 5)