Thermal transport in phononic crystals: The role of zone folding effect

Recent experiments [Yu et al., Nature Nanotech 5, 718 (2010); Tang et al., Nano Lett. 10, 4279 (2010); Hopkins et al., Nano Lett. 11, 107(2011)] on silicon based nanoscale phononic crystals demonstrated substantially reduced thermal conductivity compared to bulk Si, which cannot be explained by incoherent phonon boundary scattering within the Boltzmann Transport Equation (BTE). In this paper, partial coherent treatment of phonons, where phonons are regarded as either wave or particles depending on their frequencies, was considered. Phonons with mean free path smaller than the characteristic size of phononic crystals are treated as particles and the transport in this regime is modeled by BTE with phonon boundary scattering taken into account. On the other hand, phonons with mean free path longer than the characteristic size are treated as waves. In this regime, phonon dispersion relations are computed using the Finite Difference Time Domain (FDTD) method and are found to be modified due to the zone folding effect. The new phonon spectra are then used to compute phonon group velocity and density of states for thermal conductivity modeling. Our partial coherent model agrees well with the recent experimental results on in-plane thermal conductivity of phononic crystals. Our study highlights the importance of zone folding effect on thermal transport in phononic crystals. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.3699056]