Next Article in Journal
Co3O4 Nanopetals Grown on the Porous CuO Network for the Photocatalytic Degradation
Next Article in Special Issue
Recent Advances in Thermal Interface Materials for Thermal Management of High-Power Electronics
Previous Article in Journal
Micro-/Nanofibrillated Cellulose-Based Coating Formulations: A Solution for Improving Paper Printing Quality
Previous Article in Special Issue
Exploring the Impact of the Linker Length on Heat Transport in Metal–Organic Frameworks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

How Hydrodynamic Phonon Transport Determines the Convergence of Thermal Conductivity in Two-Dimensional Materials

Center for Phononics and Thermal Energy Science, China–EU Joint Lab for Nanophononics, MOE Key Laboratory of Advanced Micro-structured Materials, School of Physics Science and Engineering, Tongji University, Shanghai 200092, China
*
Author to whom correspondence should be addressed.
Nanomaterials 2022, 12(16), 2854; https://doi.org/10.3390/nano12162854
Submission received: 22 July 2022 / Revised: 13 August 2022 / Accepted: 15 August 2022 / Published: 18 August 2022
(This article belongs to the Special Issue Heat Transfer in Nanostructured Materials)

Abstract

:
The phonon Boltzmann transport equation combined with first-principles calculation has achieved great success in exploring the lattice thermal conductivity ( κ ) of various materials. However, the convergence of the predicted κ is a critical issue, leading to quite scattered results recorded in the literature, even for the same material. In this paper, we explore the origin for the convergence of thermal conductivity in two-dimensional (2D) materials. Two kinds of typical 2D materials, graphene and silicene, are studied, and the bulk silicon is also compared as a control system for a three-dimensional material. The effect of the cutoff radius ( r c ) in the third-order interatomic force constants on κ is studied for these three materials. It is found that that κ of these three materials exhibits diverse convergence behaviors with respect to r c , which coincides very well with the strength of hydrodynamic phonon transport. By further analyzing the phonon lifetime and scattering rates, we reveal that the dominance of the normal scattering process gives rise to the hydrodynamic phonon transport in both graphene and silicene, which results in long-range interaction and a large lifetime of low-frequency flexural acoustic phonons, while the same phenomenon is absent in bulk silicon. Our study highlights the importance of long-range interaction associated with hydrodynamic phonon transport in determining the thermal conductivity of 2D materials.

1. Introduction

As the size of novel electronic devices gradually shrinks to the nanoscale, thermal management and energy conversion at such length scale are facing huge challenges [1,2,3]. The diverse thermophysical properties make the two-dimensional (2D) materials suitable candidates for both heat dissipation [4,5,6,7,8] and thermoelectric applications [9,10,11]. Understanding the phonon thermal transport in 2D materials [12] is very important for designing new 2D materials and developing their applications. Important factors affecting the thermal transport properties of solids have been revealed, such as the phonon coherence [13,14,15,16,17], symmetry [18,19], four-phonon scattering [20,21,22,23], system size [24], and hydrodynamic phonon transport [25,26,27,28,29,30].
The phonon Boltzmann transport equation (BTE) is a powerful tool to predict the material’s lattice thermal conductivity ( κ ) and understand the underlying phonon transport mechanism [31,32,33,34,35,36,37,38]. Its prediction accuracy depends on a number of input parameters [39,40,41,42,43]. Among them, the cutoff radius ( r c ) for the nearest neighbors’ (NNs) interaction is critical in predicting the κ of 2D materials. For instance, Qin et al. [40] found that the predicted κ values for graphene by BTE are quite different with different r c , and a large r c is required to ensure the convergence of κ . Similarly, Sun et al. [44] found in four kinds of stable 2D honeycomb structures that the predicted value of κ varies greatly with different r c , which can differ by a factor of 2 if an insufficiently large r c is used, and a large r c up to the 14th NN should be used to reach convergence. In contrast, the predicted κ of three-dimensional (3D) materials is less sensitive to the choice of r c [45,46].
Possible origins for the convergence issue for κ of 2D materials have been identified. One is that the existence of long-range interatomic forces [40,47], such as resonant bonds, which were first found in rock-salt-like crystal structures, and the decoupling of the π -bond. Another is that the truncation errors in existing calculations can lead to the divergence of κ . Lindsay et al. [48] reported that a small supercell or a finite number of nearest neighbor atomic layers may break the principle of the translational invariance (TI) condition of the crystal. Similarly, Taheri et al. [35] found that the hybridization of flexural acoustic (ZA) modes and transverse/longitudinal acoustic (TA/LA) modes in 2D materials, as well as the limitation of the supercell size could break the rotational symmetry in second-order interatomic force constants (IFCs). However, the supercell used in previous works [40,44,49] is not sufficiently large, so that the r c for the NNs’ interaction exceeds half of the supercell length, which can greatly affect the prediction accuracy. Moreover, the influence of the NNs in the third-order IFCs, as well as the different sensitivity to the NNs for different dimensional systems are still unknow.
In this work, we investigate the relationship between κ and r c of the third-order IFCs in different materials, including graphene, silicene, and silicon. By calculating the value of κ with different r c , we find that the κ of graphene is the most difficult one to converge as r c increases, while silicon’s κ is not sensitive to the NNs and can converge within a small r c . Through the study of phonon–phonon scattering process, we discovered that the dominant normal scattering enables long-range interatomic force in hydrodynamic 2D materials, such as graphene. By further analysis of each phonon mode, we found that the underlying physical mechanism lies in the strong anharmonic ZA phonons, while TA and LA are insensitive to r c . Our work provides physical insights into the deeper understanding of heat transport and hydrodynamic phonon transport in different systems.

2. Computational Method

All first-principles calculations were performed based on density functional theory (DFT), as implemented in the Vienna ab initio Simulation Package (VASP) [50,51,52] with the projected augmented wave method. The Perdew–Burke–Ernzerhof (PBE) of generalized gradient approximation (GGA) was used as the exchange-correlation functional in our study, since it has already been shown to have an accurate description of a smoother potential, correct behavior under uniform scaling, and a linear response of the uniform electron gas [53]. A vacuum space of 20 Å was imposed along the c-axis only for 2D materials to exclude the interaction between adjacent imaging layers. In all three structures, the unit cell was optimized with a cutoff energy of 600 eV until the energy and the Hellmann–Feynman force converged to 10−8 eV and 10−4 eV 1 , respectively. Besides, the k-mesh grid used to sample the Brillouin zone (BZ) was set as 21 × 21 × 1 for 2D graphene and silicene and 15 × 15 × 15 for 3D silicon, respectively. The second-order IFC and the third-order IFCs were recorded by the finite displacement method as implemented in the PHONOPY code [54] and thirdorder_vasp.py code [55], respectively. We built a 5 × 5 × 1 supercell for graphene and silicene and a 3 × 3 × 3 supercell for silicon to obtain the second-order IFCs. In order to ensure a large NN in graphene, when calculating the third-order IFCs, we used a 7 × 7 × 1 supercell so that a large r c up to the 15th NN can be correctly modeled, while for another two materials, the 5 × 5 × 1 supercell and 4 × 4 × 4 supercell were used for silicene and silicon, respectively. We used a finite displacement of 0.01 Å in the PHONOPY package, which is sufficient to ensure a convergent dynamical matrix.
Previous studies revealed the important role of the four-phonon scattering process in determining lattice thermal conductivity [20,21,56]. However, due to the complexity in modeling four-phonon scattering, only a short cutoff distance up to the second nearest neighbors can be considered in the calculation of fourth-order interatomic force constants [57,58,59,60]. Since our objective was to study the effect of a long-range cutoff distance in computing thermal conductivity, we limited our study to only the three-phonon scattering process.
By iteratively solving the BTE in the framework of three-phonon scattering, the lattice thermal conductivity tensor κ α β can be obtained as implemented in the ShengBTE package [55]:
κ α β = 1 k B T 2 Ω λ f 0 ( f 0 + 1 ) ( ω λ ) 2 υ λ α F λ β ,
where λ , , ω , k B , T, Ω , f 0 , and υ denote the phonon mode index, reduced Planck constant, phonon frequency, Boltzmann constant, temperature, system volume, equilibrium Bose–Einstein distribution, and phonon group velocity, respectively. Here, α and β are the Cartesian coordinate directions and the term F λ β denotes the linear coefficient of the nonequilibrium phonon distribution function given by
F λ β = τ λ 0 ( υ λ β + Δ λ β ) ,
where τ λ 0 is the phonon lifetime obtained under single-mode relaxation time approximation (RTA). The term Δ λ β is obtained by the fully iterative solution of the BTE, which is a correction term to reflect the phonon distribution deviation from the RTA scheme. In the calculations, one can start with the RTA solution and obtain F λ β by iteratively solving the BTE.
Based on Matthiessen’s rule, the phonon lifetime τ λ 0 can be expressed as
1 τ λ 0 = 1 τ λ 3 p h + 1 τ λ i s o ,
where 1 τ λ i s o is the isotope scattering rate and 1 τ λ 3 p h   is the three-phonon scattering rate. The intrinsic scattering rate 1 τ λ 3 p h   can be calculated by Fermi’s golden rule as
1 τ λ 3 p h = 1 N ( λ 1 λ 2 + Γ λ λ 1 λ 2 + + λ 1 λ 2 1 2 Γ λ λ 1 λ 2 ) ,
where N is the total number of q points and Γ with the superscripts ( ± ) represents the scattering rates for three-phonon processes given by:
Γ λ λ 1 λ 2 + = π 4 ( f 0 λ 1 f 0 λ 2 ) ω λ ω λ 1 ω λ 2 | V λ λ 1 λ 2 + | 2 δ ( ω λ + ω λ 1 ω λ 2 ) ,  
Γ λ λ 1 λ 2 = π 4 ( f 0 λ 1 + f 0 λ 2 + 1 ) ω λ ω λ 1 ω λ 2 | V λ λ 1 λ 2 | 2 δ ( ω λ ω λ 1 ω λ 2 ) .  
The scattering matrix V λ λ 1 λ 2 ± is given by
V λ λ 1 λ 2 ± = i , j , k α β γ Φ i j k α β γ e λ i α e ± λ 1 j β e λ 2 k γ m i m j m k ,
where Φ i j k α β γ is the third-order IFCs, i, j, and k are the atomic indices, m denotes the atomic mass, and e λ i α is the α th component of the phonon eigenvector on atom i. More details about the BTE calculations can be found in our previous works [19,38].

3. Results and Discussion

The atomic structures of graphene and silicene are shown in Figure 1a,b, respectively. Although silicene has a similar 2D honeycomb structure as graphene, there exists a buckling height Δ = 0.45 Å in silicene. The calculated lattice constants for graphene and silicene are 2.47 Å and 3.86 Å, respectively, in good agreement with the literature results [41,43,61]. In comparison, Figure 1c shows the atomic structure of 3D bulk silicon, with a lattice constant of 5.47 Å, which also agrees well with the previous study [55]. The phonon dispersion (Figure S1 in the Supporting Information) reveals that the quadratic dispersion of the ZA mode is only observed in 2D graphene and silicene, while all acoustic branches are linear in 3D bulk silicon.
We first verified the calculated κ of graphene versus the Q-grid size ( N × N × 1 ). Since a 7 × 7 × 1 supercell was used for graphene in our calculations, the largest r c should be less than half of the length of the supercell due to the periodic boundary condition. This means the upper limit of r c is 8.645 Å, which corresponds to the 15th NN ( r c = 8.61 Å). The detailed relation between r c and the nth NN can be found in Figure S2 in the Supporting Information. Figure 2a shows the convergence for the calculated κ of graphene at 300 K. When the cutoff radius is small (NNs = 10th, 11th, and 12th), the value of κ diverges with the increasing Q-grid. A similar divergent κ was also reported by Qin and Hu [40] when the r c is small. However, we noticed our calculation results are quantitatively different from the results of Qin and Hu [40], in which a converged κ for NNs = 10th was reported by Qin and Hu’s work (empty hexagon in Figure 2a). A 5 × 5 × 1 supercell was used in their study, which sets the upper limit of r c around 6.17 Å. However, r c = 6.35 Å (NNs = 9th) and r c   = 6.80 Å (NNs = 10th) were used in their calculations. This discrepancy between our result and the previous study highlights the importance of using a sufficiently large supercell when increasing the cutoff radius.
When the r c further increases beyond the 12th NN, a convergent κ with respect to the Q-grid size is observed in graphene. This is because the long-range interaction can affect the lifetime of low-frequency phonons, which correspond to acoustic phonon modes near the Γ point [40]. However, the converged κ values depend notably on the specific choice of r c . For instance, the predicted κ at 300 K with the 13th, 14th, and 15th NN is 5780 Wm−1K−1, 7300 Wm−1K−1, and 4867 Wm−1K−1, respectively, exhibiting a non-monotonic dependence on r c . A further increase of r c is not possible unless an even larger supercell than 7 × 7 × 1 is used, which is beyond our computational capacity. As shown in Figure 2a, the predicted κ value for the largest r c considered in this study (NNs = 15th) agrees reasonably well with the previous theoretical and experimental results [4,62,63,64]. This complex convergence behavior in graphene is intriguing and also demonstrates that the cutoff radius is very important in determining the κ of graphene.
As another 2D material, the non-monotonic dependence of κ on r c is also observed in silicene, as shown in Figure 2b. However, when the r c is beyond the 6th NN, the predicted κ value at 300 K in silicene reaches convergence much faster compared to that in graphene. In silicene, the converged κ at the large Q-grid limit for the 6th, 7th, and 8th NN is 22.77 Wm−1K−1, 22.85 Wm−1K−1, and 21.42 Wm−1K−1, respectively. These values agree quite well with the results recorded in the literature [65,66,67,68]. In contrast, when switching to a 3D material, this significant convergence issue with respect to r c observed in the 2D material vanishes in the bulk silicon (Figure 2c). In the large Q-grid limit, the predicted κ value at 300 K in silicon for the 3rd and 4th NN is 150 Wm−1K−1 and 154 Wm−1K−1, respectively, which is also in accordance with the literature results [39,69,70,71,72,73]. By comparing these three materials, we found that the κ of the 3D material converges more easily than the 2D material, while the κ of graphene is more difficult to converge in the 2D material. This difference in the convergence of κ in the above-mentioned three materials aroused our interest in the intrinsic phonon transport form.
The different convergence behaviors of κ between the 2D and 3D materials remind us of the unique ZA branches in 2D materials. As a special feature of 2D materials, the flexural ZA phonons not only lead to unique phonon transport regimes, but also facilitate the momentum-conserving normal phonon–phonon scattering (N-scattering) [30]. When the N-scattering dominates over other phonon scatterings, such as Umklapp scattering (U-scattering), hydrodynamic phonon transport will take place [27,29].
In order to quantitatively compare the degree of hydrodynamic phonon transport among these three materials, we computed the ratio between two thermal conductivity predictions, i.e., the iterative solution ( κ I ) and relaxation time approximation ( κ RTA ). The relaxation time approximation (RTA) incorrectly treats the N-scattering as a resistive process, the same as U-scattering. As a result, RTA overestimates the phonon–phonon scattering rates and, thus, leads to an underestimation of κ compared to the iterative solution, especially when the N-scattering is dominant [26,28]. Therefore, the difference between κ I and κ RTA can serve as an indicator of hydrodynamic phonon transport.
Figure 3 shows the ratio κ I / κ RTA versus temperature for three materials. For both graphene and silicene, this ratio is larger than unity, indicating the emergence of hydrodynamic phonon transport in these 2D materials. With decreasing temperature, this ratio gradually increases in both graphene and silicene, due to the suppressed U-scattering at a low temperature. Moreover, this ratio in graphene is larger than that in silicene at a given temperature, and the decreasing trend of this ratio with decreasing temperature is much more obvious in graphene. These results suggest a stronger hydrodynamic behavior occurring in graphene. In contrast, this ratio in bulk silicon is close to unity and is insensitive to temperature, confirming that there is no hydrodynamic phonon transport in 3D silicon.
In order to understand the origin for the different degrees of hydrodynamic phonon transport, we further calculated the phonon lifetime from both iterative and RTA solutions for these three materials. In the 2D system, there exists notable differences in the phonon lifetime between the iterative and RTA solutions for both graphene (Figure 4a) and silicene (Figure 4b), especially for the low-frequency phonons near the Γ point. These low-frequency phonons actually correspond to the ZA branch. Moreover, graphene has a mirror reflection symmetry due to the atomically flat surface, while silicene has a buckled surface. The mirror reflection symmetry in graphene imposes a symmetry selection rule on ZA phonons, causing the three-phonon processes involving an odd number of flexural phonons to be forbidden [63,74,75]. This unique selection rule in graphene gives rise to a surprisingly large phonon lifetime for ZA phonons near the Γ point, which is larger than that in silicene by more than three-orders of magnitude. In contrast, there exists a negligible difference in the phonon lifetime between the iterative and RTA solutions for the bulk silicon (Figure 4c). These results reveal that the ZA phonon is responsible for the hydrodynamic phonon transport observed in 2D graphene and silicene.
We demonstrated the convergence issue of κ and hydrodynamic phonon transport in 2D materials, while the same phenomenon is absent in 3D materials. In order to understand the physical origin for this difference, we finally computed the frequency-resolved phonon scattering rate from Fermi’s golden rule for both N-scattering ( Γ N ( ω ) ) and U-scattering ( Γ U ( ω ) ) in these three materials at 300 K. To further compare the ensemble-averaged scattering strength, the averaged scattering rate is computed as [27,28,29]
Γ ¯ s = C ( ω ) Γ s ( ω ) C ( ω ) ,
where the subscripts denote the N- or U-scattering and C ( ω ) is the mode specific heat given by
C ( ω ) = f 0 ( ω ) [ f 0 ( ω ) + 1 ] ( ω ) 2 k B T 2
Figure 5a,b show the computed Γ N ( ω ) and Γ U ( ω ) with different NNs in graphene, respectively. For a given NN, the Γ N ( ω ) is much larger than Γ U ( ω ) by two-orders of magnitude in the low-frequency region, which makes the most contribution to the thermal conductivity. The ensemble-averaged scattering rate for the acoustic phonons is Γ ¯ N   = 2.843 ps 1 and Γ ¯ U   = 0.015 ps 1 . In other words, Γ ¯ N accounts for 99.5% of the total scatterings, demonstrating the dominance of N-scattering over U-scattering. This gives rise to the hydrodynamic phonon transport in graphene, as witnessed in Figure 4. In the hydrodynamic transport regime, phonons experience little scattering inside the material and, thus, can propagate collectively, which leads to very long phonon mean free path (MFP) and fast transport of thermal energy [27]. Such a long MFP and collective motion of phonons result in a long-range interaction between atoms in graphene. This effect can be revealed by the obvious dependence of Γ N ( ω ) and Γ U ( ω ) on the cutoff radius (NNs). Therefore, the hydrodynamic-transport-induced long-range interaction is responsible for the significant convergence issue of κ for graphene observed in Figure 1a.
Similar to graphene, Figure 5d,e reveal that the Γ N ( ω ) in silicene is also larger than Γ U ( ω ) in the low-frequency region, and ensemble-averaged scattering rate Γ ¯ N accounts for 76.9% of the total scatterings, which is larger than Γ ¯ U . However, the difference between Γ ¯ N and Γ ¯ U in silicene is much smaller than that in graphene. Furthermore, the scattering rate spectrum in silicene is less sensitive to the choice of NNs compared to that in graphene. These results suggest a weaker degree of hydrodynamic phonon transport in silicene compared to graphene. This in turn leads to a faster convergence of κ at a smaller NN (Figure 2b) and a smaller difference in κ between the two BTE solutions (Figure 4) in silicene, compared to the case of graphene. In contrast, when it shifts to a 3D material, Figure 5g,h show no difference in the phonon scattering rate spectrum between different NNs in the bulk silicon, and the averaged scattering rate in Figure 5i confirms that the U-scattering is the dominant scattering process. This results in the vanishing of hydrodynamic phonon transport (Figure 4) and, consequently, no convergence issue of κ for the bulk silicon, even at very small NNs (Figure 2c).

4. Conclusions

In summary, we calculated the room temperature thermal conductivity κ of graphene, silicene, and bulk silicon within different cutoff radius r c . In contrast to the fast convergence of κ at very small r c in the 3D bulk silicon, the κ of the 2D materials exhibits a persistent dependence on r c in both graphene and silicene, which is more pronounced in graphene even up to the 15th nearest neighbor. By comparing the κ values predicted from the iterative solution and relaxation time approximation, the hydrodynamic phonon transport was confirmed in both graphene and silicene, while it is absent in the bulk silicon. More interestingly, the sensitivity of κ with r c in these three materials correlates very well with the strength of the hydrodynamic transport behavior. The calculations of phonon lifetime revealed that the presence of a unique ZA phonon in 2D materials is a key factor for the hydrodynamic phonon transport observed in graphene and silicene. By further analyzing the scattering rates, we demonstrated that the dominance of the N-scattering process causes the hydrodynamic transport, so that phonons can propagate collectively, which results in a long-range interaction in a 2D material such as graphene compared to the 3D material. Therefore, the long-range interaction induced by the strong hydrodynamic phonon transport is responsible for the significant convergence behavior of κ observed in graphene. This work uncovers the underlying physical mechanism for the notable convergence behavior of κ in 2D systems.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/nano12162854/s1, Figure S1: The phonon dispersion of graphene, silicene, and silicon. Figure S2: The corresponding cutoff radius ( r c ) to the nth nearest neighbors (NNs) for graphene, silicene, and silicon.

Author Contributions

J.C. conceived the idea of the study and supervised the project; J.J. performed the simulations and conducted the data analysis; S.L. and Y.O. contributed to the theoretical data analysis; all authors contributed to the writing and editing of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This project is supported in part by grants from the National Natural Science Foundation of China (Grant Nos. 11,890,703 and 12,075,168), the Science and Technology Commission of Shanghai Municipality (Grant No. 21JC1405600), and the Fundamental Research Funds for the Central Universities (Grant No. 22120220060).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

κ thermal conductivity
r c cutoff radius
NNnearest neighbor
reduced Planck constant
ω phonon frequency
k B Boltzmann constant
Ttemperature
Ω system volume
f 0 equilibrium Bose–Einstein distribution
υ phonon group velocity
τ λ 0 phonon lifetime
C ( ω ) mode specific heat
α ,   β Cartesian coordinate directions
matomic mass
i, j, katomic indices
κ I thermal conductivities from the iterative solution
κ RTA thermal conductivities from the relaxation time approximation
Γ N normal scattering rates
Γ U Umklapp scattering rates
Γ ¯ N averaged normal scattering rates
Γ ¯ U averaged Umklapp scattering rates

References

  1. Chen, J.; Xu, X.; Zhou, J.; Li, B. Interfacial thermal resistance: Past, present, and future. Rev. Mod. Phys. 2022, 94, 025002. [Google Scholar] [CrossRef]
  2. Qian, X.; Zhou, J.; Chen, G. Phonon-engineered extreme thermal conductivity materials. Nat. Mater. 2021, 20, 1188–1202. [Google Scholar] [CrossRef] [PubMed]
  3. Ouyang, Y.; Zhang, Z.; Li, D.; Chen, J.; Zhang, G. Emerging Theory, Materials, and Screening Methods: New Opportunities for Promoting Thermoelectric Performance. Ann. Phys. 2019, 531, 1800437. [Google Scholar] [CrossRef]
  4. Balandin, A.A.; Ghosh, S.; Bao, W.; Calizo, I.; Teweldebrhan, D.; Miao, F.; Lau, C.N. Superior thermal conductivity of single-layer graphene. Nano Lett. 2008, 8, 902–907. [Google Scholar] [CrossRef]
  5. Cai, Q.; Scullion, D.; Gan, W.; Falin, A.; Zhang, S.; Watanabe, K.; Taniguchi, T.; Chen, Y. High thermal conductivity of high-quality monolayer boron nitride and its thermal expansion. Sci. Adv. 2019, 5, eaav0129. [Google Scholar] [CrossRef]
  6. He, J.; Ouyang, Y.; Yu, C.; Jiang, P.; Ren, W.; Chen, J. Lattice thermal conductivity of β12 and x3 borophene. Chin. Phys. B 2020, 29, 126503. [Google Scholar] [CrossRef]
  7. Yu, J.; Zhou, J.; Wan, X.; Li, Q. High intrinsic lattice thermal conductivity in monolayer MoSi2N4. New J. Phys. 2021, 23, 033005. [Google Scholar] [CrossRef]
  8. Ren, W.; Ouyang, Y.; Jiang, P.; Yu, C.; He, J.; Chen, J. The Impact of Interlayer Rotation on Thermal Transport Across Graphene/Hexagonal Boron Nitride van der Waals Heterostructure. Nano Lett. 2021, 21, 2634–2641. [Google Scholar] [CrossRef]
  9. Su, L.; Wang, D.; Wang, S.; Qin, B.; Wang, Y.; Qin, Y.; Jin, Y.; Chang, C.; Zhao, L.-D. High thermoelectric performance realized through manipulating layered phonon-electron decoupling. Science 2022, 375, 1385–1389. [Google Scholar] [CrossRef]
  10. Roychowdhury, S.; Ghosh, T.; Arora, R.; Samanta, M.; Xie, L.; Singh, N.K.; Soni, A.; He, J.; Waghmare, U.V.; Biswas, K. Enhanced atomic ordering leads to high thermoelectric performance in AgSbTe2. Science 2021, 371, 722–727. [Google Scholar] [CrossRef]
  11. He, J.; Hu, Y.; Li, D.; Chen, J. Ultra-low lattice thermal conductivity and promising thermoelectric figure of merit in borophene via chlorination. Nano Res. 2022, 15, 3804–3811. [Google Scholar] [CrossRef]
  12. Xu, X.; Chen, J.; Li, B. Phonon thermal conduction in novel 2D materials. J. Phys. Condens. Matter. 2016, 28, 483001. [Google Scholar] [CrossRef] [PubMed]
  13. Zhang, Z.; Guo, Y.; Bescond, M.; Chen, J.; Nomura, M.; Volz, S. Coherent thermal transport in nano-phononic crystals: An overview. APL Mater. 2021, 9, 081102. [Google Scholar] [CrossRef]
  14. An, M.; Chen, D.; Ma, W.; Hu, S.; Zhang, X. Directly visualizing the crossover from incoherent to coherent phonons in two-dimensional periodic MoS2/MoSe2 arrayed heterostructure. Int. J. Heat Mass Transf. 2021, 178, 121630. [Google Scholar] [CrossRef]
  15. Zhang, Z.; Guo, Y.; Bescond, M.; Chen, J.; Nomura, M.; Volz, S. Generalized decay law for particlelike and wavelike thermal phonons. Phys. Rev. B 2021, 103, 184307. [Google Scholar] [CrossRef]
  16. Zhang, Z.; Guo, Y.; Bescond, M.; Chen, J.; Nomura, M.; Volz, S. Heat Conduction Theory Including Phonon Coherence. Phys. Rev. Lett. 2022, 128, 015901. [Google Scholar] [CrossRef]
  17. Zhang, Z.; Guo, Y.; Bescond, M.; Chen, J.; Nomura, M.; Volz, S. How coherence is governing diffuson heat transfer in amorphous solids. NPJ Comput. Mater. 2022, 8, 96. [Google Scholar] [CrossRef]
  18. Gupta, R.; Dongre, B.; Bera, C.; Carrete, J. The Effect of Janus Asymmetry on Thermal Transport in SnSSe. J. Phys. Chem. C 2020, 124, 17476–17484. [Google Scholar] [CrossRef]
  19. Lu, S.; Ren, W.; He, J.; Yu, C.; Jiang, P.; Chen, J. Enhancement of the lattice thermal conductivity of two-dimensional functionalized MXenes by inversion symmetry breaking. Phys. Rev. B 2022, 105, 165301. [Google Scholar] [CrossRef]
  20. Feng, T.; Lindsay, L.; Ruan, X. Four-phonon scattering significantly reduces intrinsic thermal conductivity of solids. Phys. Rev. B 2017, 96, 161201. [Google Scholar] [CrossRef]
  21. Feng, T.; Ruan, X. Four-phonon scattering reduces intrinsic thermal conductivity of graphene and the contributions from flexural phonons. Phys. Rev. B 2018, 97, 045202. [Google Scholar] [CrossRef]
  22. Yu, C.; Hu, Y.; He, J.; Lu, S.; Li, D.; Chen, J. Strong four-phonon scattering in monolayer and hydrogenated bilayer BAs with horizontal mirror symmetry. Appl. Phys. Lett. 2022, 120, 132201. [Google Scholar] [CrossRef]
  23. Ouyang, Y.; Yu, C.; He, J.; Jiang, P.; Ren, W.; Chen, J. Accurate description of high-order phonon anharmonicity and lattice thermal conductivity from molecular dynamics simulations with machine learning potential. Phys. Rev. B 2022, 105, 115202. [Google Scholar] [CrossRef]
  24. Zhang, Z.; Ouyang, Y.; Cheng, Y.; Chen, J.; Li, N.; Zhang, G. Size-dependent phononic thermal transport in low-dimensional nanomaterials. Phys. Rep. 2020, 860, 1–26. [Google Scholar] [CrossRef]
  25. Guo, Y.; Wang, M. Phonon hydrodynamics and its applications in nanoscale heat transport. Phys. Rep. 2015, 595, 1–44. [Google Scholar] [CrossRef]
  26. Cepellotti, A.; Fugallo, G.; Paulatto, L.; Lazzeri, M.; Mauri, F.; Marzari, N. Phonon hydrodynamics in two-dimensional materials. Nat. Commun. 2015, 6, 6400. [Google Scholar] [CrossRef]
  27. Lee, S.; Broido, D.; Esfarjani, K.; Chen, G. Hydrodynamic phonon transport in suspended graphene. Nat. Commun. 2015, 6, 6290. [Google Scholar] [CrossRef]
  28. Zhang, Z.; Ouyang, Y.; Guo, Y.; Nakayama, T.; Nomura, M.; Volz, S.; Chen, J. Hydrodynamic phonon transport in bulk crystalline polymers. Phys. Rev. B 2020, 102, 195302. [Google Scholar] [CrossRef]
  29. Ding, Z.; Zhou, J.; Song, B.; Chiloyan, V.; Li, M.; Liu, T.H.; Chen, G. Phonon Hydrodynamic Heat Conduction and Knudsen Minimum in Graphite. Nano Lett. 2018, 18, 638–649. [Google Scholar] [CrossRef]
  30. Yu, C.; Ouyang, Y.; Chen, J. A perspective on the hydrodynamic phonon transport in two-dimensional materials. J. Appl. Phys. 2021, 130, 010902. [Google Scholar] [CrossRef]
  31. Lindsay, L.; Li, W.; Carrete, J.; Mingo, N.; Broido, D.A.; Reinecke, T.L. Phonon thermal transport in strained and unstrained graphene from first principles. Phys. Rev. B 2014, 89, 155426. [Google Scholar] [CrossRef]
  32. Jain, A.; McGaughey, A.J. Strongly anisotropic in-plane thermal transport in single-layer black phosphorene. Sci. Rep. 2015, 5, 8501. [Google Scholar] [CrossRef] [PubMed]
  33. Zhang, Z.; Chen, J.; Li, B. Negative Gaussian curvature induces significant suppression of thermal conduction in carbon crystals. Nanoscale 2017, 9, 14208–14214. [Google Scholar] [CrossRef]
  34. Lindsay, L.; Hua, C.; Ruan, X.L.; Lee, S. Survey of ab initio phonon thermal transport. Mater. Today Phys. 2018, 7, 106–120. [Google Scholar] [CrossRef]
  35. Taheri, A.; Pisana, S.; Singh, C.V. Importance of quadratic dispersion in acoustic flexural phonons for thermal transport of two-dimensional materials. Phys. Rev. B 2021, 103, 235426. [Google Scholar] [CrossRef]
  36. Liu, P.-F.; Bo, T.; Xu, J.; Yin, W.; Zhang, J.; Wang, F.; Eriksson, O.; Wang, B.-T. First-principles calculations of the ultralow thermal conductivity in two-dimensional group-IV selenides. Phys. Rev. B 2018, 98, 235426. [Google Scholar] [CrossRef]
  37. Taheri, A.; Da Silva, C.; Amon, C.H. Phonon thermal transport in β−NX(X = P,As,Sb) monolayers: A first-principles study of the interplay between harmonic and anharmonic phonon properties. Phys. Rev. B 2019, 99, 235425. [Google Scholar] [CrossRef]
  38. Ren, W.; Zhang, Z.; Chen, C.; Ouyang, Y.; Li, N.; Chen, J. Phononic Thermal Transport in Yttrium Hydrides Allotropes. Front. Mater. 2020, 7, 569090. [Google Scholar] [CrossRef]
  39. Jain, A.; McGaughey, A.J.H. Effect of exchange–correlation on first-principles-driven lattice thermal conductivity predictions of crystalline silicon. Comput. Mater. Sci. 2015, 110, 115–120. [Google Scholar] [CrossRef]
  40. Qin, G.; Hu, M. Accelerating evaluation of converged lattice thermal conductivity. NPJ Comput. Mater. 2018, 4, 3. [Google Scholar] [CrossRef]
  41. Xie, H.; Gu, X.; Bao, H. Effect of the accuracy of interatomic force constants on the prediction of lattice thermal conductivity. Comput. Mater. Sci. 2017, 138, 368–376. [Google Scholar] [CrossRef]
  42. Taheri, A.; Da Silva, C.; Amon, C.H. First-principles phonon thermal transport in graphene: Effects of exchange-correlation and type of pseudopotential. J. Appl. Phys. 2018, 123, 215105. [Google Scholar] [CrossRef]
  43. Qin, G.; Qin, Z.; Wang, H.; Hu, M. On the diversity in the thermal transport properties of graphene: A first-principles-benchmark study testing different exchange-correlation functionals. Comput. Mater. Sci. 2018, 151, 153–159. [Google Scholar] [CrossRef]
  44. Sun, Z.; Yuan, K.; Chang, Z.; Bi, S.; Zhang, X.; Tang, D. Ultra-low thermal conductivity and high thermoelectric performance of two-dimensional triphosphides (InP3, GaP3, SbP3 and SnP3): A comprehensive first-principles study. Nanoscale 2020, 12, 3330–3342. [Google Scholar] [CrossRef] [PubMed]
  45. Lindroth, D.O.; Erhart, P. Thermal transport in van der Waals solids from first-principles calculations. Phys. Rev. B 2016, 94, 115205. [Google Scholar] [CrossRef]
  46. Kundu, A.; Yang, X.; Ma, J.; Feng, T.; Carrete, J.; Ruan, X.; Madsen, G.K.H.; Li, W. Ultrahigh Thermal Conductivity of theta-Phase Tantalum Nitride. Phys. Rev. Lett. 2021, 126, 115901. [Google Scholar] [CrossRef]
  47. Lee, S.; Esfarjani, K.; Luo, T.; Zhou, J.; Tian, Z.; Chen, G. Resonant bonding leads to low lattice thermal conductivity. Nat. Commun. 2014, 5, 4525. [Google Scholar] [CrossRef]
  48. Lindsay, L.; Broido, D.A.; Reinecke, T.L. Ab initiothermal transport in compound semiconductors. Phys. Rev. B 2013, 87, 165201. [Google Scholar] [CrossRef]
  49. Wang, N.; Shen, C.; Sun, Z.; Xiao, H.; Zhang, H.; Yin, Z.; Qiao, L. High-Temperature Thermoelectric Monolayer Bi2TeSe2 with High Power Factor and Ultralow Thermal Conductivity. ACS Appl. Energ. Mater. 2022, 5, 2564–2572. [Google Scholar] [CrossRef]
  50. Kresse, G.; Furthmuller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. Mater. Sci. 1996, 6, 15–50. [Google Scholar] [CrossRef]
  51. Kresse, G.; Furthmuller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 1996, 54, 11169–11186. [Google Scholar] [CrossRef] [PubMed]
  52. Kresse, G.; Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B 1999, 59, 1758–1775. [Google Scholar] [CrossRef]
  53. Perdew, J.; Burke, K.; Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 1996, 78, 3865–3868. [Google Scholar] [CrossRef]
  54. Togo, A.; Tanaka, I. First principles phonon calculations in materials science. Scr. Mater. 2015, 108, 1–5. [Google Scholar] [CrossRef]
  55. Li, W.; Carrete, J.; Katcho, N.A.; Mingo, N. ShengBTE: A solver of the Boltzmann transport equation for phonons. Comput. Phys. Commun. 2014, 185, 1747–1758. [Google Scholar] [CrossRef]
  56. Feng, T.; Ruan, X. Quantum mechanical prediction of four-phonon scattering rates and reduced thermal conductivity of solids. Phys. Rev. B 2016, 93, 045202. [Google Scholar] [CrossRef]
  57. Tong, Z.; Yang, X.; Feng, T.; Bao, H.; Ruan, X. First-principles predictions of temperature-dependent infrared dielectric function of polar materials by including four-phonon scattering and phonon frequency shift. Phys. Rev. B 2020, 101, 125416. [Google Scholar] [CrossRef]
  58. Feng, T.; O’Hara, A.; Pantelides, S.T. Quantum prediction of ultra-low thermal conductivity in lithium intercalation materials. Nano Energy 2020, 75, 104916. [Google Scholar] [CrossRef]
  59. Yang, X.; Feng, T.; Li, J.; Ruan, X. Stronger role of four-phonon scattering than three-phonon scattering in thermal conductivity of III-V semiconductors at room temperature. Phys. Rev. B 2019, 100, 245203. [Google Scholar] [CrossRef]
  60. Gu, X.; Fan, Z.; Bao, H.; Zhao, C.Y. Revisiting phonon-phonon scattering in single-layer graphene. Phys. Rev. B 2019, 100, 064306. [Google Scholar] [CrossRef]
  61. Xie, H.; Ouyang, T.; Germaneau, É.; Qin, G.; Hu, M.; Bao, H. Large tunability of lattice thermal conductivity of monolayer silicene via mechanical strain. Phys. Rev. B 2016, 93, 075404. [Google Scholar] [CrossRef]
  62. Fugallo, G.; Cepellotti, A.; Paulatto, L.; Lazzeri, M.; Marzari, N.; Mauri, F. Thermal conductivity of graphene and graphite: Collective excitations and mean free paths. Nano Lett. 2014, 14, 6109–6114. [Google Scholar] [CrossRef] [PubMed]
  63. Lindsay, L.; Broido, D.A.; Mingo, N. Flexural phonons and thermal transport in graphene. Phys. Rev. B 2010, 82, 115427. [Google Scholar] [CrossRef]
  64. Ghosh, S.; Calizo, I.; Teweldebrhan, D.; Pokatilov, E.P.; Nika, D.L.; Balandin, A.A.; Bao, W.; Miao, F.; Lau, C.N. Extremely high thermal conductivity of graphene: Prospects for thermal management applications in nanoelectronic circuits. Appl. Phys. Lett. 2008, 92, 151911. [Google Scholar] [CrossRef]
  65. Gu, X.; Yang, R. First-principles prediction of phononic thermal conductivity of silicene: A comparison with graphene. J. Appl. Phys. 2015, 117, 025102. [Google Scholar] [CrossRef]
  66. Qin, G.; Qin, Z.; Yue, S.Y.; Yan, Q.B.; Hu, M. External electric field driving the ultra-low thermal conductivity of silicene. Nanoscale 2017, 9, 7227–7234. [Google Scholar] [CrossRef]
  67. Han, Y.; Dong, J.; Qin, G.; Hu, M. Phonon transport in the ground state of two-dimensional silicon and germanium. RSC Adv. 2016, 6, 69956–69965. [Google Scholar] [CrossRef]
  68. Liu, Z.; Wu, X.; Luo, T. The impact of hydrogenation on the thermal transport of silicene. 2D Mater. 2017, 4, 025002. [Google Scholar] [CrossRef]
  69. Parrish, K.D.; Jain, A.; Larkin, J.M.; Saidi, W.A.; McGaughey, A.J.H. Origins of thermal conductivity changes in strained crystals. Phys. Rev. B 2014, 90, 235201. [Google Scholar] [CrossRef]
  70. Broido, D.A.; Malorny, M.; Birner, G.; Mingo, N.; Stewart, D.A. Intrinsic lattice thermal conductivity of semiconductors from first principles. Appl. Phys. Lett. 2007, 91, 231922. [Google Scholar] [CrossRef]
  71. Fulkerson, W.; Moore, J.P.; Williams, R.K.; Graves, R.S.; McElroy, D.L. Thermal Conductivity, Electrical Resistivity, and Seebeck Coefficient of Silicon from 100 to 1300 °K. Phys. Rev. 1968, 167, 765–782. [Google Scholar] [CrossRef]
  72. Glassbrenner, C.J.; Slack, G.A. Thermal Conductivity of Silicon and Germanium from 3 °K to the Melting Point. Phys. Rev. 1964, 134, A1058–A1069. [Google Scholar] [CrossRef]
  73. Shanks, H.R.; Maycock, P.D.; Sidles, P.H.; Danielson, G.C. Thermal Conductivity of Silicon from 300 to 1400 °K. Phys. Rev. 1963, 130, 1743–1748. [Google Scholar] [CrossRef]
  74. Nika, D.L.; Pokatilov, E.P.; Askerov, A.S.; Balandin, A.A. Phonon thermal conduction in graphene: Role of Umklapp and edge roughness scattering. Phys. Rev. B 2009, 79, 155413. [Google Scholar] [CrossRef]
  75. Peng, B.; Zhang, D.; Zhang, H.; Shao, H.; Ni, G.; Zhu, Y.; Zhu, H. The conflicting role of buckled structure in phonon transport of 2D group-IV and group-V materials. Nanoscale 2017, 9, 7397–7407. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The atomic structure of (a) graphene, (b) silicene, and (c) bulk silicon. Graphene has an atomically flat surface, while the surface of silicene has a finite buckling height Δ .
Figure 1. The atomic structure of (a) graphene, (b) silicene, and (c) bulk silicon. Graphene has an atomically flat surface, while the surface of silicene has a finite buckling height Δ .
Nanomaterials 12 02854 g001
Figure 2. The room temperature thermal conductivity κ with respect to the Q-grid for different nearest neighbors (NNs). (a) Graphene. The literature results were adapted from Qin et al. [40], Fugallo et al. [62], Lindsay et al. [63], Balandin et al. [4], and Ghosh et al. [64]. (b) Silicene. The literature results were adapted from Gu et al. [65], Qin et al. [66], Han et al. [67], and Liu et al. [68]. (c) Bulk silicon. The literature results were adapted from Parrish et al. [69], Broido et al. [70], Jain et al. [39], Fulkerson et al. [71], Glassbrenner et al. [72], and Shanks et al. [73].
Figure 2. The room temperature thermal conductivity κ with respect to the Q-grid for different nearest neighbors (NNs). (a) Graphene. The literature results were adapted from Qin et al. [40], Fugallo et al. [62], Lindsay et al. [63], Balandin et al. [4], and Ghosh et al. [64]. (b) Silicene. The literature results were adapted from Gu et al. [65], Qin et al. [66], Han et al. [67], and Liu et al. [68]. (c) Bulk silicon. The literature results were adapted from Parrish et al. [69], Broido et al. [70], Jain et al. [39], Fulkerson et al. [71], Glassbrenner et al. [72], and Shanks et al. [73].
Nanomaterials 12 02854 g002
Figure 3. Thermal conductivities from the iterative solution ( κ I ) and relaxation time approximation ( κ RTA ) of the BTE for graphene, silicene, and silicon.
Figure 3. Thermal conductivities from the iterative solution ( κ I ) and relaxation time approximation ( κ RTA ) of the BTE for graphene, silicene, and silicon.
Nanomaterials 12 02854 g003
Figure 4. The phonon lifetime calculated from the iterative solution and relaxation time approximation at 300 K in (a) graphene, (b) silicene, and (c) bulk silicon. Huge differences in the phonon lifetime for the low-frequency ZA mode are observed in graphene and silicene.
Figure 4. The phonon lifetime calculated from the iterative solution and relaxation time approximation at 300 K in (a) graphene, (b) silicene, and (c) bulk silicon. Huge differences in the phonon lifetime for the low-frequency ZA mode are observed in graphene and silicene.
Nanomaterials 12 02854 g004
Figure 5. The frequency-resolved N-scattering rate ( Γ N ( ω ) ) and U-scattering rate ( Γ U ( ω ) ) at 300 K with different NNs in (a,b) graphene, (d,e) silicene, and (g,h) silicon. The ensemble-averaged scattering rates ( Γ ¯ N and Γ ¯ U ) for three acoustic phonons are shown in (c) graphene, (f) silicene, and (i) silicon, respectively.
Figure 5. The frequency-resolved N-scattering rate ( Γ N ( ω ) ) and U-scattering rate ( Γ U ( ω ) ) at 300 K with different NNs in (a,b) graphene, (d,e) silicene, and (g,h) silicon. The ensemble-averaged scattering rates ( Γ ¯ N and Γ ¯ U ) for three acoustic phonons are shown in (c) graphene, (f) silicene, and (i) silicon, respectively.
Nanomaterials 12 02854 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jiang, J.; Lu, S.; Ouyang, Y.; Chen, J. How Hydrodynamic Phonon Transport Determines the Convergence of Thermal Conductivity in Two-Dimensional Materials. Nanomaterials 2022, 12, 2854. https://doi.org/10.3390/nano12162854

AMA Style

Jiang J, Lu S, Ouyang Y, Chen J. How Hydrodynamic Phonon Transport Determines the Convergence of Thermal Conductivity in Two-Dimensional Materials. Nanomaterials. 2022; 12(16):2854. https://doi.org/10.3390/nano12162854

Chicago/Turabian Style

Jiang, Jianhui, Shuang Lu, Yulou Ouyang, and Jie Chen. 2022. "How Hydrodynamic Phonon Transport Determines the Convergence of Thermal Conductivity in Two-Dimensional Materials" Nanomaterials 12, no. 16: 2854. https://doi.org/10.3390/nano12162854

APA Style

Jiang, J., Lu, S., Ouyang, Y., & Chen, J. (2022). How Hydrodynamic Phonon Transport Determines the Convergence of Thermal Conductivity in Two-Dimensional Materials. Nanomaterials, 12(16), 2854. https://doi.org/10.3390/nano12162854

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop