HYPOTHESIS AND THEORY article

Front. Energy Res., 13 June 2018

Sec. Nanotechnology for Energy Applications

Volume 6 - 2018 | https://doi.org/10.3389/fenrg.2018.00048

A Modified Theoretical Model to Accurately Account for Interfacial Roughness in Predicting the Interfacial Thermal Conductance

  • YZ

    Yingying Zhang 1,2†

  • DM

    Dengke Ma 1,3†

  • YZ

    Yi Zang 1,3

  • XW

    Xiaojia Wang 2

  • NY

    Nuo Yang 1,3*

  • 1. State Key Laboratory of Coal Combustion, Huazhong University of Science and Technology, Wuhan, China

  • 2. Department of Mechanical Engineering, University of Minnesota, Minneapolis, MN, United States

  • 3. Nano Interface Center for Energy, School of Energy and Power Engineering, Huazhong University of Science and Technology, Wuhan, China

Abstract

The acoustic mismatch model and the diffuse mismatch model (DMM) have been widely used to predict the thermal interface conductance. However, the acoustic mismatch model and the DMM are based on the hypothesis of a perfectly smooth interface and a completely disordered interface respectively. Here, we present a new modified model, named as the mixed mismatch model, which considers the roughness/bonding at the interface. By taking partially specular and partially diffuse transmission into account, the mixed mismatch model (MMM) can predict the thermal interface conductance with arbitrary roughness. The proportions of specular and diffuse transmission are determined by the interface roughness which is described by the interfacial density of states. We show that the predicted results of the MMM match well with the values of molecular dynamics simulation and experimental data.

Thermal transport across interfaces is an important issue for microelectronics, photonics, and thermoelectric devices, and has been studied both experimentally and theoretically recently (Swartz and Pohl, ; Zhu et al., ; Cahill et al., ; Wu et al., ; Zheng et al., ). Generally, the interfacial thermal conductance (ITC) is used to evaluate the physical properties of thermal transport in devices and materials (Yang et al., ), such as composites (Nan et al., ; Deng et al., ) superlattices (Dresselhaus et al., ), thin-film multilayers (Cahill et al., ), nanoscale devices (Goodson and Ju, ; Han et al., ), and nanocrystalline materials (Soyez et al., ). Therefore, a deep understanding and an accurate prediction of ITC are crucial to improve the design and performance of various devices and materials.

To date, two theoretical models, the acoustic mismatch model (AMM) and the diffuse mismatch model (DMM), have been widely used in predicting the ITC (Swartz and Pohl, ). The AMM makes an essential simplified assumption that phonons incident at an interface undergo specular transmission and are governed by continuum mechanics (Swartz and Pohl, ). Contrarily, the DMM assumes that the interface is completely disordered and phonons lose their memory after reaching the interface (Little, ; Swartz and Pohl, ). Although these two models have advanced the understanding of thermal transport across interface, both the AMM and DMM purely consider the bulk properties, and ignore the interfacial roughness (Tian et al., ; Merabia and Termentzidis, ; Ramiere et al., ; Shao et al., ), interfacial modes(Chalopin and Volz, ; Gong et al., ; Yang et al., ; Gordiz and Henry, ; Lee and Luo, ), electron-phonon coupling (Huang and Kaviany, ; Hohensee et al., ), and inelastic scattering at the interface (Xiong et al., ), which leads to the inaccuracy of the two models (Swartz and Pohl, ; Wu and Luo, ).

In the past decades, some works have improved the accuracy of the AMM and DMM, and sophisticated modifications to the original models that account for the inelastic scattering have been proposed. For example, the maximum transmission model (Dames and Chen, ), the higher harmonic inelastic model (Hopkins, ), the joint frequency DMM (Hopkins and Norris, ), the scattering-mediated acoustic mismatch model (Prasher and Phelan, ), and the anharmonic inelastic model (Hopkins et al., ). In addition, several literature studies have included the electron-phonon couplings at metal/non-metal interfaces (Majumdar and Reddy, ; Wang et al., ; Wu and Luo, ).

In 1960, Ziman analyzed the partially specular and partially diffuse reflection of phonons at surface (Ziman, ). Then Chen applied this idea to interface scattering model to predict the thermal conductivity of superlattice structures (Chen, ). The idea is worth to be applied in calculating the ITC. Besides, it is found that the phonon interface states which are localized close to the interfacial region and different from bulk states have a large effect on ITC (Shin et al., ; Liang et al., ; Murakami et al., ; Gordiz and Henry, ; Polanco et al., ). However, all of the existing theoretical models do not consider the interface states. Therefore, there is a great demand for a theoretical model that accounts for interface states in predicting the ITC.

In this letter, we developed a new model to predict the ITC, named as the mixed mismatch model (MMM). Based on partially specular and partially diffuse transmissions, by taking interface states into account, the MMM can predict the ITC with arbitrary roughness and overcomes the shortcomings of AMM and DMM. Here, we considered the interface states in obtaining the proportions of diffuse and specular transmission. To validate the accuracy of the MMM, we compared the ITC values of MMM to those by MD simulations and experimental measurements.

At the interface between material A and B, an incident phonon with frequency ω and mode j can either scatter back or transmit (Swartz and Pohl, ). According to the Landauer formalism (Chen, ), we can predict the ITC, G, as

where ω is the frequency, ωv is the cutoff frequency, D is the phonon density of states, n(ω,T) is the Bose-Einstein distribution function, v is the phonon group velocity for the phonon mode of interest at frequency ω, and α is the phonon transmission coefficient. The subscript “j” refers to the phonon polarization.

In the AMM and DMM, the transmission coefficient can be calculated as (Swartz and Pohl, ; Cahill et al., )

where ρ is the mass density.

According to previous studies (Soffer, ; Chen, ; Chen and Neagu, ; Phelan, ; Bellis et al., ; Zhao and Freund, ), for solid-solid interfaces, both diffuse and specular transmissions are taking place at the interface. Hence, an assumed specular parameter, p, can be used to represent the proportion of phonons specularly transmitted across the interface (Soffer, ; Chen, ). Hence, in the MMM, we define the transmission coefficient by linearly combining AMM and DMM terms as

According to Ziman (Soffer, ), the p is related to the root mean square roughness (RMSR), η, and to phonon wavelength, λ. Ziman postulated the relationship as

Generally, for simplification the roughness is considered as a constant for a given surface. While the roughness of the surface is measured by the height of the asperities, and is a rather rapidly changing function at the point where the height and wavelength are about equal (Ziman, ). It means that constant value of η is too restrictive to be applied to real surfaces. The whole interface should be divided into many different small sampling areas whose size is relative to a few wavelengths in dimensions (Ziman, ). Regarding to different wavelength of phonons, the size of the sampling area is different. According to statistics from different sampling area, the roughness is different and depends on phonon frequency.

As we know, the η depends on the interfacial structure (e.g., the arrangement of atoms at the interface). Recently, it is found that the interfacial structure is related to interfacial density of states (IDOS), which may differ from that of the bulk materials (Yang et al., ; Gordiz and Henry, ). That is, the η depends on IDOS. For the specular interface corresponding to the AMM case, the interface is treated as a perfect plane and the η is zero (Swartz and Pohl, ). To make the interface as a perfect plane, the atoms close to interface stay at the equilibrium situation as in the bulk (Dbulk,A and Dbulk,B). Therefore, the atoms close to interface vibrate as the bulk and the IDOS should be the same as the DOS of the bulk. Generally, interface is consisting of two materials which have different force constant, they interact with each other and will definitely distort the interfacial structure. In this case, the specular plane does not exist in interface (Chalopin and Volz, ; Gordiz and Henry, ). On the contrary, for the completely diffuse interface corresponding to the DMM case, all the phonons are diffusely scattered and the acoustic correlations are completely destroyed (Swartz and Pohl, ). And the structure is similar to amorphous structures where each atom has a distinct local environment (Kaviany, ) and the acoustic correlations are completely destroyed. Thus, we speculate that the IDOS of a completely diffuse interface should also be similar to the DOS of an amorphous structure (Damo,A and Damo,B). Since both the real interface structure and its IDOS (Dint,A and Dint,B) are between the two extremities, there is a mapping between RMSR and IDOS. We choose to use a simple mapping function of RMSR on the density of states, expressed as:

where C is a constant, which does not depend on the frequency. The Dbulk, Damo are the DOS in the bulk and in amorphous structure, respectively. It is worth noting that here we obtain the value of C by MD simulations, but not limited to MD simulations, the value of C can be obtained by experiment or any other simulation methods.

According to Equations (5) and (6), the specular parameter, p, is obtained

This is the key equation of this letter. By substituting Equation (7) into Equation (4), we can get the transmission coefficient α, then the ITC can be figured out by Equation (1). It is shown that the interfacial structure or IDOS will greatly affect the phonon transport. As shown in Figure 3, for the perfectly smooth interface, the specular reflection does not change the energy and the momentum components along the direction of the temperature gradient (Chen, ), thus the IDOS should be the same as the DOS of the bulk. Contrarily, for the extremely rough interface, the IDOS is similar to the DOS of amorphous structures. For a partially specular and partially diffuse interface, the IDOS is in between the two extreme situations.

The MD simulation and green function are being widely used to predict the ITC (Landry and McGaughey, ; Ma et al., ). While there is no assumption on the phonon scattering in MD (Ma et al., ). To demonstrate the accuracy of the MMM, we calculated the ITC of interfaces with different roughness, and then compared the results of MMM with MD values. In our MD simulations (The MD simulation details are shown in Supporting information (Figure 1d and Figure SI), we create two Al/Si interface structures with different roughness, structure (a) and (b) (shown in Figure SI). The structures after relaxation are shown in Figures 1a,b. The interface roughness of structure (b) is larger than that of structure (a). And the IDOS of structure (a) and (b) is calculated as the average phonon modes of atoms at the interface. And the interface thickness is set as 0.815 nm. The phonon group velocity of each mode is calculated by sorting the phonon dispersion relationship. (as shown in Figure SII) We also simulate an amorphous (completely disordered) structure (shown in Figure 1c) to obtain the Damo,Al and Damo,Si in Equation (8).

Figure 1

The main results are shown in Figure 2A. The ITC of AMM and DMM correspond to the ITC of the perfectly smooth interface and completely disordered interface, respectively. The interfaces of structure (a) and (b) are partially diffuse and partially specular. We calculated the ITC of the two interfacial structures [structure (a) with smaller roughness and structure (b) with larger roughness] by MD. We perform several simulations with different initial states for one structure, and the value shown here is the average value, where the error bar is the standard deviation. The results of MD show that both of them fall in AMM and DMM predictions, and the ITC of structure (b) is larger than that of structure (a). In the calculation of MMM, the value of C is unknown, but can be obtained by fitting to the MD data of structure (a), as 5.75 × 10−11. Then we applied Equation (7) to predict ITC of structure (b) by MMM which matches well with the data of MD (Figure 2A). It declares the validity of MMM. It is worth noting that here we obtain the value of C by MD, but not limited to MD, the value of C can be get by experiment or any other simulation method.

Figure 2

The specular parameter is the most important parameter in the MMM which presents the proportion of phonons specular transmitted across the interface. As shown in Figure 2B, the dependence of specular parameters on frequency for difference interface structures was calculated and compared. (The specular parameters of other structures are shown in Figure SV) There is an obvious decline in the p-curves of both structures (a) and (b). For lower frequency phonons in each p-curve, with longer wavelength, it shows that the values of p for both (a) and (b) are close to one. That is, the phonons with longer wavelength are not sensitive to the nanoscale interface roughness and transit across them similar to a specular interface. For higher frequency phonons, the shorter wavelength phonons transmit across the interface with more diffusive scatterings. Besides, the decline frequency of (a), ~6.5 THz, is higher than that of (b), ~3.5 THz, because the roughness of structure (a) is smaller than that of structure (b). As shown in Equation (5) the IDOS dominates the value of p, which makes the MMM consider the interfacial structure. To show the relationship between the IDOS and the interfacial structure, we calculated the IDOS of different interfacial structures. Figures 3A,B show the IDOS of a perfectly smooth interface in AMM and a completely disordered interface in DMM, corresponding to the DOS of bulk and amorphous structure [Figure S1 (f)], respectively. While Figures 3C,D show the IDOS of structures (a) and (b), respectively. (The IDOS of other structures are shown in Figures SIII, SIV) By comparing the difference of IDOS between Al and Si, we find that the IDOS of Al and Si become similar with each other as the interface roughness increases because the disorder structure at interface diminishes the bulk properties.

Figure 3

To further validate the MMM, we compare the transmission coefficient calculated by the MMM with the experimental results (Hua et al., ). Minnich et al. obtained the results by interpreting data from TDTR experiment with ab initio phonon transport model. In Figure 4, it shows that the transmission coefficients of the MMM match well with the experimental values. Besides, we can see that the DMM underestimates the transmission coefficient of lower-frequency phonons because the phonons with longer wavelengths are less likely to be scattered. Contrarily, the AMM overestimates the transmission coefficient of high frequency phonons due to the phonons with shorter wavelengths are more likely to be scattered.

Figure 4

).

It is worth noting that in experiments, the oxidation of interface is a very common issue (Hopkins et al., , ; Lee et al., ), which will affect phonon transport. Generally, the native oxides will increase the interfacial roughness and induce more scattering, thus decrease the transmission coefficient (Hopkins et al., ; Mu et al., ). While if the oxidation strengthens the chemical bond of the interface, this effect will be much stronger than the induced roughness scattering, the transmission coefficient will increase, and the ITC will increase (Hopkins et al., ). In MMM, IDOS is used to describe the interfacial information, which includes the increasing of roughness and strengthening of bond caused by oxides. Therefore, the MMM is also valid to deal with interfaces which have oxidation.

In summary, we have developed a new theoretical model, named as the mixed mismatch model, which predicts the thermal interface conductance by considering the interfacial structure. The MMM takes into account partially specular and partially diffuse transmissions using the IDOS to determine the proportion of specular parameter. The value of specular parameter has an obvious dependence on the phonon wavelength and interfacial roughness. The MMM is validated by comparing its prediction with both MD results and measurements. Both ITC and the transmission coefficients calculated by the MMM agree well with those obtained from the direct MD simulation and from experiment.

Statements

Author contributions

YyZ and DM contributed equally to proposing the new model in this manuscript, and YZ helped validate the data. XW contributed to revising the manuscript and replying to reviewers' comments. NY provided essential guidance to this work.

Acknowledgments

NY was sponsored by National Natural Science Foundation of China (No. 51576076 and No. 51711540031), Hubei Provincial Natural Science Foundation of China (2017CFA046) and Fundamental Research Funds for the Central Universities (2016YXZD006). XW was sponsored by the National Science Foundation through the University of Minnesota MRSEC (DMR-1420013) and the Legislative-Citizen Commission on Minnesota Resources. The authors are grateful to Hua Bao, and Xin Qian for useful discussions. Besides, the authors thank the National Supercomputing Center in Tianjin (NSCC-TJ) and China Scientific Computing Grid (ScGrid) for providing assistance in computations.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fenrg.2018.00048/full#supplementary-material

References

Summary

Keywords

acoustic mismatch model, diffuse mismatch model, thermal interface conductance, mixed mismatch model, interfacial density of states

Citation

Zhang Y, Ma D, Zang Y, Wang X and Yang N (2018) A Modified Theoretical Model to Accurately Account for Interfacial Roughness in Predicting the Interfacial Thermal Conductance. Front. Energy Res. 6:48. doi: 10.3389/fenrg.2018.00048

Received

12 December 2017

Accepted

22 May 2018

Published

13 June 2018

Volume

6 - 2018

Edited by

Wey Yang Teoh, City University of Hong Kong, Hong Kong

Reviewed by

Tandeep Singh Chadha, Applied Particle Technology LLC, United States; Sreeram Vaddiraju, Texas A&M University, United States

Updates

Copyright

*Correspondence: Nuo Yang

This article was submitted to Nanoenergy Technologies and Materials, a section of the journal Frontiers in Energy Research

†These authors have contributed equally to this work.

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics