Critical angle for interfacial phonon scattering:results from ab initio lattice dynamics calculations

Critical angle for interfacial phonon scattering:
results from ab initio lattice dynamics calculations

A. Alkurdi [    S. Pailhès Institut Lumière Matière, UMR5306 Université Claude Bernard Lyon 1-CNRS, Université de Lyon 69622 Villeurbanne Cedex, France    S. Merabia
July 15, 2019

Thermal boundary resistance is a critical quantity that controls heat transfer at the nanoscale, and which is primarily related to interfacial phonon scattering. Here, we combine lattice dynamics calculations and inputs from first principles ab initio simulations to predict phonon transmission at Si/Ge interface as a function of both phonon frequency and phonon wavevector. This technique allows us to determine the overall thermal transmission coefficient as a function of the phonon scattering direction and frequency. Our results show that, the thermal energy transmission is highly anisotropic while thermal energy reflection is almost isotropic. In addition, we found the existence of a global critical angle of transmission beyond which almost no thermal energy is transmitted. This critical angle around . is found to be almost independent on the interaction range between Si and Ge, the interfacial bonding strength, and the temperature above K. We interpret these results by carrying out a spectral and angular analysis of phonon transmission coefficient and differential thermal boundary conductance.

Thermal boundary resistance, interfacial phonon scattering, lattice dynamics
Thermal conduction in metals and alloys and semiconductors, 66.70.Df, -Lattice dynamics crystals (see 63), Phonons scattering by, 72.10.Di
preprint: AIP/123-QED

Also at ] Department of Physics, Al-Baath University, Homs, Syria

I Introduction

Thermal transport at the nanoscale is at the heart of modern microelectronics. With the continuous reduction of the size of the microprocessors, there is a strong need to understand heat flow across interfaces from a fundamental standpoint. Indeed, in the modern microelectronic era, transistors reach nanometer length scales, and their relevant dimensions become comparable with the Kapitza length where is the thermal conductivity and the thermal boundary conductance (TBC) characterizing the interface. Typical values for the Kapitza length are between and nm, showing the importance of thermal interface resistance for nanostructured materials Pop2010 ().

In conventional metal/dielectrics interfaces, heat transport is mediated by both phonon-phonon and electron-phonon scattering Lombard2015 (). For interfaces between semiconductors, only phonon-phonon scattering is relevant, but even this simpler situation is far from being understood. In particular, the thermal energy transmission coefficient has not been measured experimentally and thermal process at interfaces remain poorly understood. Notable exceptions may be found where phonon reflection and internal reflection have been evidenced experimentally Northtrop1984 (); Höss1990 (), but these manifestations concern sub-Thz frequencies, a situation where continnum acoustics and thus Snell-Descartes laws should be at work. In contrast to phonon scattering, reflection and transmission coefficients for photon interfacial scattering are well described by Fresnel coefficients Hua2015 (). Even the involved problem of the intensity scattered by a rough interface has been adressed theoretically by Maradudin and coworkers Simonsen2010 (). Such a milestone has not been reached for phonon scattering and for thermal transport at crystalline interfaces in general. The reasons are manifold: i) crystals are anisotropic; ii) phonon dispersion is not as simple as electromagnetic dispersion; iii) phonon transmission coefficient should depend on the interfacial bonding strength between the two solids.

In order to rationalize early experiments on thermal boundary resistance at low and moderate temperatures, phenomenological models have been proposed, including the popular acoustic mismatch model (AMM) Little1959 (); Chen1999 () and the diffuse mismatch model (DMM) Swartz1989 (). However, these two models generally fail in describing accurately interfacial thermal transport as they do not include any information relevant to the interface, and usually consider phonon dispersion in a simplified way. In the quest for more accurate models, several computational methods have been proposed during the last decade, including molecular dynamics Landry2009 (), Green’s function Zhang2007 () and lattice dynamics Young1989 (); Zhao2005 ().

Molecular dynamics is a flexible tool to describe bulk phonon properties and interfacial thermal transport, and its accuracy would be considerably improved by the use of interatomic potentials deriving from ab-initio calculations. Green’s function calculations may challenge these limitations, but they have not been devised to compute the frequency and wavevector dependence of the transmission coefficient, at least for an interface between two 3D materials Ong2015 ().

In this Letter, we combine lattice dynamics calculations with ab initio interatomic force constants in order to predict phonon transmission as a function of both phonon frequency and scattered angle, for a Si/Ge interface. Our results show that phonon transmission is highly anisotropic, while phonon reflection is almost isotropic. We demonstrate that almost all the thermal energy is transmitted at the interface for scattered angles smaller than a critical value which is found to be typically for the system considered here. Through the consideration of different numbers of interacting atoms at the interface, we systematically studied the effect of the interfacial interaction range and the interfacial bonding strength, and conclude that these factors mildly affect the angular distribution of phonon scattering. In contrast, temperature is found to have more influence, as a result of the phonon occupation statistics.

Ii Theoretical Model

The model we developed is inspired by the lattice dynamics model for Si/Ge interface, as described in Zhao2005 (). We have extended this latter model basically to employ interatomic force constants inputs from ab initio calculations. In this perspective, it was necessary to modify the equations of motion detailed in Zhao2005 () so as to account for the interaction with th neighbouring atoms, while the original work by Zhao and Freund considers an empirical Stillinger-Weber potential which is short range and simpler to tackle, but which lacks the accuracy to describe phonon bulk dispersion curves especially for Ge. The interaction in each medium is supposed to be determined by interatomic force constants derived from ab initio calculations Aouissi2006 () that have been calculated employing density functional pertubation theory, a pseudopotential approach and the local density approximation (LDA) for the exchange-correlation potential. As shown in the supplementary material, the phonon spectra are in good agreement with experimental data for Si and Ge. Note that the phonons dispersions curves obtained with the force constants described above doesn’t not show negative phonon branch.

Under the harmonic approximation, the equation of motion for an atom of mass displaced by can be read as :  where the sum is over all atoms including atom itself and is the interatomic force constant matrix limited here to harmonic terms. The equation of motion can be rewritten as :  where is the wavevector, denotes the branch mode and is the dynamical matrix here having dimensions . The interface is represented in Fig. 1, where we divide the system into three regions: the left lead where the incident phonons come from, the central region corresponding to the location of the interface and the right lead where the transmitted phonons come out, denoted by L, C and R, respectively. The interface is oriented to be perpendicular to the -axis corresponding to the [001] crystalline direction.

Figure 1: Schematic of the perfect interface between Si/Ge. Coming from the left lead, the incident phonon () strikes the interface in the central region and is decomposed into transmitted phonons in the right lead () and reflected phonons in the left lead ().

The dynamical matrix can thus be decomposed into three or more region-depending dynamical matrices as Zhao2005 (): where and denotes the direction normal to the interface. We assume elastic scattering at the interface, an incident phonon can only excite phonons having the same frequency and the components of the wave vector parallel to the interface are conserved Pettersson1990 () . These equations have complex displacement vectors as solution, of the th atom in the th primitive unit cell in the left lead which can be expressed as a superposition of the incident mode and the excited modes in the left lead :


and similarly, the displacement vector in the outcoming medium (right lead) is


where , are the amplitudes of the reflected and transmitted waves, and are the numbers of allowed reflected and transmitted modes, respectively (here for the diamond like systems considered). By solving this set of equations at the interface, we obtain the amplitude of each excited mode leading to the knowledge of the individual mode-dependent energy transmission coefficient which is defined as the fraction of phonon energy transmitted across the interface :


where and are the mass density in the incidence and transmission medium, respectively. are the group velocities of the transmitted phonons projected along the direction perpendicular to the interface. Similarly, the individual reflection coefficient is :


The mode-dependent thermal conductance has the expression :


where is the volume of the system and denotes the phonon occupation density and the sum is over the transmitted modes. In this work, the ab-initio interatomic force constants are taken from Ref. Aouissi2006 (). Unless specified, we consider that the interatomic force constants at the interface are equal to the arithmetic mean of Si and Ge force constants, which can be summarized as : with obvious notations. We note that our approach neglects anharmonic effects as a result of the elastic interfacial scattering assumption, but as noticed in Tian2012 (); Landry2009 () anharmonic effects become important for Si/Ge at temperatures higher than K.

To critically assess the method, we have compared the total thermal conductance calculated using the transmission coefficient and found a value (MW/m2/K) at high temperatures. This value is in good agreement with the conductance reported in a DFT-Green’s function calculation for the same interface Tian2012 (), thus validating our methodology. Note also, that the resulting value of the TBC is smaller by a factor than the value predicted by the empirical Stillinger-Weber potential Landry2009 (); Zhao2005 (), showing the necessity to employ ab-initio interatomic force constants.

Iii Results and Discussion

We present now the results obtained by performing the LD calculations on Si/Ge interface described in the previous section. The interatomic interaction was limited to the th unit cell (corresponding to th nearest neighbours for diamond-like crystals). Fig. 2 shows the transmission and reflection of all phonon modes in polar coordinates, and for different interaction ranges that goes from the first unit cell (L1) up to the th unit cell (L4). It is important to mention that we have performed an intensive statistical treatment, and considered millions excited incident wavevectors. The scattering angle is defined using group velocities . Overall, the transmission coefficient displays a maximum value close to in the normal transmission direction corresponding to and it decreases when we go to the grazing direction, . From this analysis, we can see that the phonon transmission is anisotropic with relative low values lying in the interval of transmission angle between and . This transmission angle beyond which the transmission is small can be seen as a global critical angle for phonon scattering. Thermal phonon reflection is to be almost isotropic with maximum value about in the direction of of transmission angle. The transmission coefficient is found to be almost insensitive to the interaction range, but the reflection coefficient displays little variations when the interaction goes up to the th unit cell.

Figure 2: (Color online) (a) Phonon energy transmission at Si/Ge in polar coordinates for different interaction ranges. L1 refers to the interaction limited to only the nearest unit cell while L4 refers to the interaction range up to the 4th unit cell.

We have performed a spectral and angular analysis to quantify the contribution of each phonon mode in a given scattering direction, as depicted in Fig. 3. Again, this figure presents the result of million excited modes. We can split the phonon frequency range, which extends up to the Ge cut-off frequency, in three domains. For low frequency phonon modes (rad/ps) the transmission probability is small and we can clearly see the existence of a critical angle around for zero frequency modes. Intermediate frequency modes, (rad/ps), corresponding to longitudinal acoustic modes of Germanium have a large probability to be transmitted for small scattered angles only, . The corresponding contribution of heat flux is high, as can be inferred from Fig. 3 b). Lastly, Germanium optical modes having frequencies (rad/ps) are characterized by a large transmission coefficients for angles , but the amount of heat they carry is small. Hence, these modes only contribute to the heat capacity without transmitting significant heat flux.

Figure 3: (Color online) (a) : Colormap of the phonon energy transmission at Si/Ge interface as a function of both the angular frequency and the transmission angle. (b) : Same for the thermal boundary conductance (in units of MW ps/m2K rad). Note the log scale used in the color map. The interaction range was taken up to the nd unit cell here.

In Fig. 4, we plotted the cumulative TBC as a function of transmission angle for different interatomic interaction ranges (a) and for different interfacial interaction strengths (b). In both figures, we see again that the significant contribution to the energy transmission across the interface comes from the angular directions corresponding to transmission angles smaller than . The TBC exhibits a significant drop in magnitude when the interfacial interaction range goes beyond the first unit cell (L1). This behaviour may be surprising and contrary to the intuition, nevertheless, it demonstrates that increasing the range of the interaction does not necessarily imply enhanced phonon transmission. The conductance converges when the interfacial interaction range spans more than three unit cells. If we now examine the role of the interaction strength, taken as a percentage of the mean value of the bonding in the bulk of both media (Si and Ge), Fig. 4(b) shows that the TBC increases with the interfacial interaction strength up to a value where we have a maximum and then the TBC starts to decrease for larger interfacial bonding strengths. We are interested to know up to which direction 98% of thermal energy is transmitted to the outcoming medium (Ge) and how the corresponding critical angle changes with respect to the interfacial bonding strengths and the interatomic interaction range as well. To this end, we plotted the limiting angle in the insets of Fig. 4 (a) and (b). We found that the variations of this limiting angle with respect to the interatomic interaction range are not that large. Moreover, the limiting angle is found to be not so sensitive to the interfacial bonding strength, but when the relative bonding strength is larger than 100 % , the corresponding limiting angle slightly levels off.

Figure 4: (Color online) The cumulative TBC at room temperature for Si/Ge interface as a function of transmission angle for different interaction ranges (a) and different interfacial bonding strengths (b) taken as a percentage of the mean value of the bonding in the bulk of both media (Si and Ge). The insets show the critical angle corresponding to 98% of the total transmitted energy.

Finally, we illustrate in Fig. 5 the temperature dependence of the cumulative thermal boundary conductance as a function of the transmission angle. It is obvious that the TBC increases with respect to temperature and saturates at high temperatures. The inset shows the variations of the limiting angle with respect to the temperature. The limiting angle clearly increases with temperature and saturates for a temperature higher than K. The increase is related to the phonon statistics and the progressive population of high frequency phonon modes which may transfer energy at large scattering angles , as can be inferred from fig. 3.

Figure 5: (Color online) The cumulative TBC of Si/Ge interface as a function of transmission angle for different temperatures. The inset shows the critical angle corresponding to 98 % of the total transmitted energy as a function of temperature.

Iv Conclusion

In summary, we have used LD calculations combined with ab initio interatomic force constants to investigate thermal transport at silicon/germanium interface. These calculations have allowed us to quantify the transmission and reflection coefficients as a function of scattered angle, enlighting the important contribution of phonons transmitted along the normal direction. A spectral analysis has been carried out to elaborate a complete view on the contribution of each phonon mode in a given scattering direction. We have also demonstrated the existence of a limiting scattering angle beyond which thermal transmission is negligible. This limiting angle has been shown to depend mildly on the interatomic interaction range and the interfacial bonding strength, but increases with temperature. Further improvement of the methodology, applied here to the case of the Si/Ge interface, will need to address more precisely the force constants at the interface by means of ab initio calculation in large supercells, especially for interfaces between materials that do not have the same bulk crystalline structure. This work may be considered as a first step towards the development of multi-scale theoretical analysis of phonon transport in real devices with nanostructured interfaces Merabia2014 ().

V Supplementary material

See supplementary material for the phonon spectra of Silicon and Germanium calculated with the ab initio force constants.

Vi Acknowledgements

The authors acknowledge fruitful discussions with P.-O. Chapuis, C. Adessi, R. Viennois and K. Termentzidis. Partial financial support from the ANR Mascoth is also acknowledged.



  • (1) E. Pop, Energy dissipation and transport in nanoscale devices, Nano Research 3, 147-169 (2010)
  • (2) J. Lombard, F. Detcheverry and S. Merabia, Influence of the electron-phonon interfacial conductance on the thermal transport at metal/dielectric interfaces, J. Phys.: Condens. Matter 27 (2015) 015007
  • (3) I. Simonsen, A. A. Maradudin, and T. A. Leskova, “Scattering of electromagnetic waves from two-dimensional randomly rough perfectly conducting surfaces: The full angular intensity distribution,” Phys. Rev. A 91, 013806(1-13) (2010)
  • (4) G. A. Northtrop and J. P. Wolfe, Phonon reflection imaging: a determination of specular versus diffuse boundary scattering Phys. Rev. Lett. 52, 2156 (1984)
  • (5) C. Hoss, J. P. Wolfe and H. Kinder, Total internal reflection of high frequency phonons: a test of specular refraction Phys. Rev. Lett. 64, 1134 (1990)
  • (6) Hua C., Chen X., Ravichandran N. K. and Minnich, A. J., Fresnel transmission coefficients for thermal phonons at solid interfaces, arXiv:1509.07806
  • (7) P.L. Kapitza, J. Phys. (USSR), 4 181 (1941)
  • (8) W.A. Little, The Transport of Heat Between Dissimilar Solids at Low Temperatures Can. J. Phys. 37 334 (1959)
  • (9) G. Chen, Phonon Wave Heat Conduction in Thin Films and Superlattices J. Heat Transf. 121 945 (1999)
  • (10) E.T. Swartz and R.O. Pohl, Thermal boundary resistance Rev. Mod. Phys. 61 605 (1989)
  • (11) D.G. Cahill, W.K. Ford, K.E. Goodson, G.D. Mahan, A. Majumdar, H.J. Maris, R. Merlin and S.R. Philpot, Nanoscale thermal transport J. App. Phys. 93 793 (2003)
  • (12) E.S. Landry and A.J.H. McGaughey, Phys. Rev. B 80 165304 (2009)
  • (13) D.A. Young, and H.J. Maris, Phys. Rev. B 40 3685 (1989)
  • (14) W. Zhang, T.S. Fisher and N. Mingo, AMSE J. Heat Transfer 129 (2007) 483
  • (15) L. Zhang, P. Keblinski, J.S. Wang, and B. Li, Phys. Rev. B 83 064303 (2011)
  • (16) R.J. Stoner and H.J. Maris, Phys. Rev. B 48 16373 (1993)
  • (17) M. Aouissi, I. Hamdi, N. Meskini and A. Qteish, Phonon spectra of diamond, Si, Ge, and : Calculations with real-space interatomic force constants Phys. Rev. B 74 054302 (2006)
  • (18) H. Zhao and J.B. Freund, J. Appl. Phys. 97 024903 (2005)
  • (19) S. Pettersson and G.D. Mahan, Phys. Rev. B 42 7386-7390 (1990)
  • (20) Z. Tian, K. Esfarjani, and G. Chen, Phys. Rev. B 86 235304 (2012)
  • (21) Z.-Y. Ong and G. Zhang, Efficient approach for modeling phonon transmission probability in nanoscale interfacial thermal transport Phys. Rev. B 91 174302 (2015)
  • (22) S. Merabia and K. Termentzidis, Thermal boundary conductance across rough interfaces probed by molecular dynamics Phys. Rev. B 89 054309 (2014)
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description