Interfacial thermal resistance

Last updated

Interfacial thermal resistance, also known as thermal boundary resistance, or Kapitza resistance, is a measure of resistance to thermal flow at the interface between two materials. While these terms may be used interchangeably, Kapitza resistance technically refers to an atomically perfect, flat interface whereas thermal boundary resistance is a more broad term. [1] This thermal resistance differs from contact resistance (not to be confused with electrical contact resistance) because it exists even at atomically perfect interfaces. Owing to differences in electronic and vibrational properties in different materials, when an energy carrier (phonon or electron, depending on the material) attempts to traverse the interface, it will scatter at the interface. The probability of transmission after scattering will depend on the available energy states on side 1 and side 2 of the interface.

Contents

Assuming a constant thermal flux is applied across an interface, this interfacial thermal resistance will lead to a finite temperature discontinuity at the interface. From an extension of Fourier's law, we can write

where is the applied flux, is the observed temperature drop, is the thermal boundary resistance, and is its inverse, or thermal boundary conductance.

Understanding the thermal resistance at the interface between two materials is of primary significance in the study of its thermal properties. Interfaces often contribute significantly to the observed properties of the materials. This is even more critical for nanoscale systems where interfaces could significantly affect the properties relative to bulk materials.

Low thermal resistance at interfaces is technologically important for applications where very high heat dissipation is necessary. This is of particular concern to the development of microelectronic semiconductor devices as defined by the International Technology Roadmap for Semiconductors in 2004 where an 8 nm feature size device is projected to generate up to 100000 W/cm2 and would need efficient heat dissipation of an anticipated die level heat flux of 1000 W/cm2 which is an order of magnitude higher than current devices. [2] On the other hand, applications requiring good thermal isolation such as jet engine turbines would benefit from interfaces with high thermal resistance. This would also require material interfaces which are stable at very high temperature. Examples are metal-ceramic composites which are currently used for these applications. High thermal resistance can also be achieved with multilayer systems.

As stated above, thermal boundary resistance is due to carrier scattering at an interface. The type of carrier scattered will depend on the materials governing the interfaces. For example, at a metal-metal interface, electron scattering effects will dominate thermal boundary resistance, as electrons are the primary thermal energy carriers in metals.

Two widely used predictive models are the acoustic mismatch model (AMM) and the diffuse mismatch model (DMM). The AMM assumes a geometrically perfect interface and phonon transport across it is entirely elastic, treating phonons as waves in a continuum. On the other hand, the DMM assumes scattering at the interface is diffusive, which is accurate for interfaces with characteristic roughness at elevated temperatures.

Molecular dynamics (MD) simulations are a powerful tool to investigate interfacial thermal resistance. Recent MD studies have demonstrated that the solid-liquid interfacial thermal resistance is reduced on nanostructured solid surfaces by enhancing the solid-liquid interaction energy per unit area, and reducing the difference in vibrational density of states between solid and liquid. [3]

Theoretical models

The primary model that has historically described Kapitza resistance is the phonon gas model. [4] [5] [6] Within this model are the acoustic mismatch and diffuse mismatch models (AMM and DMM respectively). For both models the interface is assumed to behave exactly as the bulk on either side of the interface (e.g. bulk phonon dispersions, velocities, etc.), with hybrid vibrational modes and the phonons that occupy them being completely neglected. In addition, the AMM and DMM models are based only on elastic phonon transport, usually ignoring electrical contributions, although it is possible to take electron contributions into account within the phonon gas model. [7] The AMM and DMM models should apply for interfaces where at least one of the materials is electrically insulating. The thermal resistance then results from the transfer of phonons across the interface. Energy is transferred when higher energy phonons which exist in higher density in the hotter material propagate to the cooler material, which in turn transmits lower energy phonons, creating a net energy flux. [8]

According to the AMM and DMM models, a crucial factor in determining the thermal resistance at an interface is the overlap of phonon states. Specifically, the models completely disregard the effects of inelastic scattering and multiple phonon interactions. For example, the models only allow for a phonon occupying a particular mode frequency to interact with another phonon occupying a mode of exactly the same frequency. In reality, however, this is not the case and the interaction probability of two phonons can be calculated using perturbation theory (quantum mechanics). As an example within the AMM and DMM models, given two materials A and B, if material A has a low population (or no population) of phonons with certain k value, there will be very few phonons of that wavevector (or equivalently, frequency) to propagate from A to B. Likewise, due to the principle of detailed balance, AMM and DMM predict that very few phonons of that wavevector will propagate in the opposite direction, from B to A, even if material B has a large population of phonons with that wavevector. Thus as the overlap between phonon dispersions is small, there are fewer modes to allow for heat transfer in the material, giving at a high thermal interfacial resistance relative to materials with a high degree of overlap. [9] Neither model is very effective for predicting the thermal interface resistance (with the exception of very low temperature), but rather for most materials they act as upper and lower limits for real behavior.

The AMM and DMM differ in the conditions they require for propagation across the interface, because the models differ greatly in their treatment of scattering at the interface. In AMM the interface is assumed to be perfect, resulting in no scattering, thus phonons propagate elastically across the interface. The wavevectors that propagate across the interface are determined by conservation of momentum. In DMM, the opposite extreme is assumed, a perfectly scattering interface. In this case the wavevectors that propagate across the interface are random and independent of incident phonons on the interface. For both models the detailed balance must still be obeyed.

For both models the basic equations of the phonon gas model apply. [4] [5] [6] The flux of energy from one material to the other in one dimension is just:

where is the group velocity which is approximated to be the speed of sound in the material for the AMM and DMM models, is the number of phonons at a given wavevector, E is the energy, and α is the probability of transmission across the interface. The net flux is thus the difference of the energy fluxes:

Since both fluxes are dependent on T1 and T2, the relationship between the flux and the temperature difference can be used to determine the thermal interface resistance based on:

where A is the area of the interface. These basic equations form the basis for both models. n is determined based on the dispersion relation for the materials (for example, the Debye model) and Bose–Einstein statistics. Energy is given simply by the De Broglie Wavelength Equation:

where . The main difference between the two models is the transmission probability, whose determination is more complicated. In each case it is determined by the basic assumptions that form the respective models. The assumption of elastic scattering makes it more difficult for phonons to transmit across the interface, resulting in lower probabilities. As a result, the acoustic mismatch model typically represents an upper limit for thermal interface resistance, while the diffuse mismatch model represents the lower limit. [10]

Examples

Liquid helium interfaces

Typical Interfacial Resistance of Liquid Helium with metals. Resistance has been multiplied by T to remove the expected T dependence. Adapted from InterfacialHeliumResistance.png
Typical Interfacial Resistance of Liquid Helium with metals. Resistance has been multiplied by T to remove the expected T dependence. Adapted from

The presence of thermal interface resistance, corresponding to a discontinuous temperature across an interface was first proposed from studies of liquid helium in 1936. While this idea was first proposed in 1936, [9] it wasn't until 1941 when Pyotr Kapitsa (Peter Kapitza) carried out the first systematic study of thermal interface behavior in liquid helium. [11] The first major model for heat transfer at interfaces was the acoustic mismatch model which predicted a T−3 temperature dependence on the interfacial resistance, but this failed to properly model the thermal conductance of helium interfaces by as much as two orders of magnitude. Another surprising behavior of the thermal resistance was observed in the pressure dependence. Since the speed of sound is a strong function of temperature in liquid helium, the acoustic mismatch model predicts a strong pressure dependence of the interfacial resistance. Studies around 1960 surprisingly showed that the interfacial resistance was nearly independent of pressure, suggesting that other mechanisms were dominant.

The acoustic mismatch theory predicted a very high thermal resistance (low thermal conductance) at solid-helium interfaces. This is problematic for researchers working at ultra-cold temperatures because it greatly impedes cooling rates at low temperatures such as in dilution refrigerators. [12] Fortunately such a large thermal resistance was not observed due to many mechanisms which promoted phonon transport. In liquid helium, Van der Waals forces actually work to solidify the first few monolayers against a solid. This boundary layer functions much like an anti-reflection coating in optics, so that phonons which would typically be reflected from the interface actually would transmit across the interface. This also helps to understand the pressure independence of the thermal conductance. The final dominant mechanism to anomalously low thermal resistance of liquid helium interfaces is the effect of surface roughness, which is not accounted for in the acoustic mismatch model. For a more detailed theoretical model of this aspect see the paper by A. Khater and J. Szeftel. [13] Like electromagnetic waves which produce surface plasmons on rough surfaces, phonons can also induce surface waves. When these waves eventually scatter, they provide another mechanism for heat to transfer across the interface. Similarly, phonons are also capable of producing evanescent waves in a total internal reflection geometry. As a result, when these waves are scattered in the solid, additional heat is transferred from the helium beyond the prediction of the acoustic mismatch theory. For a more complete review on this topic see the review by Swartz. [14]

Notable room temperature thermal conductance

In general there are two types of heat carriers in materials: phonons and electrons. The free electron gas found in metals is a very good conductor of heat and dominates thermal conductivity. All materials though exhibit heat transfer by phonon transport so heat flows even in dielectric materials such as silica. Interfacial thermal conductance is a measure of how efficiently heat carriers flow from one material to another. The lowest room temperature thermal conductance measurement to date is the Bi/Hydrogen-terminated diamond with a thermal conductance of 8.5 MW m−2 K−1. As a metal, bismuth contains many electrons which serve as the primary heat carriers. Diamond on the other hand is a very good electrical insulator (although it has a very high thermal conductivity) and so electron transport between the materials is nil. Further, these materials have very different lattice parameters so phonons do not efficiently couple across the interface. Finally, the Debye temperature between the materials is significantly different. As a result, bismuth, which has a low Debye temperature, has many phonons at low frequencies. Diamond on the other hand has a very high Debye temperature and most of its heat-carrying phonons are at frequencies much higher than are present in bismuth. [15]

Thermal Conductance Data adapted from, InterfacialThermalConductance.png
Thermal Conductance Data adapted from,

Increasing in thermal conductance, most phonon mediated interfaces (dielectric-dielectric and metal-dielectric) have thermal conductances between 80 and 300 MW m−2 K−1. The largest phonon mediated thermal conductance measured to date is between TiN (Titanium Nitride) and MgO. These systems have very similar lattice structures and Debye temperatures. While there are no free electrons to enhance the thermal conductance of the interface, the similar physical properties of the two crystals facilitate a very efficient phonon transmission between the two materials. [11]

At the highest end of the spectrum, one of the highest thermal conductances measured is between aluminum and copper. At room temperature, the Al-Cu interface has a conductance of 4 GW m−2 K−1. The high thermal conductance of the interface should not be unexpected given the high electrical conductivity of both materials. [16]

Interfacial resistance in carbon nanotubes

The superior thermal conductivity of Carbon nanotubes makes it an excellent candidate for making composite materials. But interfacial resistance impacts the effective thermal conductivity. This area is not well studied and only a few studies have been done to understand the basic mechanism of this resistance. [17] [18]

Related Research Articles

The thermal conductivity of a material is a measure of its ability to conduct heat. It is commonly denoted by , , or and is measured in W·m−1·K−1.

Conduction is the process by which heat is transferred from the hotter end to the colder end of an object. The ability of the object to conduct heat is known as its thermal conductivity, and is denoted k.

In physics, a phonon is a collective excitation in a periodic, elastic arrangement of atoms or molecules in condensed matter, specifically in solids and some liquids. A type of quasiparticle, a phonon is an excited state in the quantum mechanical quantization of the modes of vibrations for elastic structures of interacting particles. Phonons can be thought of as quantized sound waves, similar to photons as quantized light waves. However, photons are fundamental particles that can be individually detected, whereas phonons, being quasiparticles, are an emergent phenomenon.

<span class="mw-page-title-main">Thermoelectric materials</span> Materials whose temperature variance leads to voltage change

Thermoelectric materials show the thermoelectric effect in a strong or convenient form.

<span class="mw-page-title-main">Drude model</span> Model of electrical conduction

The Drude model of electrical conduction was proposed in 1900 by Paul Drude to explain the transport properties of electrons in materials. Basically, Ohm's law was well established and stated that the current J and voltage V driving the current are related to the resistance R of the material. The inverse of the resistance is known as the conductance. When we consider a metal of unit length and unit cross sectional area, the conductance is known as the conductivity, which is the inverse of resistivity. The Drude model attempts to explain the resistivity of a conductor in terms of the scattering of electrons by the relatively immobile ions in the metal that act like obstructions to the flow of electrons.

<span class="mw-page-title-main">Electromigration</span> Movement of ions in an electrical field

Electromigration is the transport of material caused by the gradual movement of the ions in a conductor due to the momentum transfer between conducting electrons and diffusing metal atoms. The effect is important in applications where high direct current densities are used, such as in microelectronics and related structures. As the structure size in electronics such as integrated circuits (ICs) decreases, the practical significance of this effect increases.

In solid-state physics, the electron mobility characterises how quickly an electron can move through a metal or semiconductor when pulled by an electric field. There is an analogous quantity for holes, called hole mobility. The term carrier mobility refers in general to both electron and hole mobility.

<span class="mw-page-title-main">Seebeck coefficient</span> Measure of voltage induced by change of temperature

The Seebeck coefficient of a material is a measure of the magnitude of an induced thermoelectric voltage in response to a temperature difference across that material, as induced by the Seebeck effect. The SI unit of the Seebeck coefficient is volts per kelvin (V/K), although it is more often given in microvolts per kelvin (μV/K).

<span class="mw-page-title-main">Wiedemann–Franz law</span> Law of physics

In physics, the Wiedemann–Franz law states that the ratio of the electronic contribution of the thermal conductivity (κ) to the electrical conductivity (σ) of a metal is proportional to the temperature (T).

In condensed matter physics, a spin wave is a propagating disturbance in the ordering of a magnetic material. These low-lying collective excitations occur in magnetic lattices with continuous symmetry. From the equivalent quasiparticle point of view, spin waves are known as magnons, which are bosonic modes of the spin lattice that correspond roughly to the phonon excitations of the nuclear lattice. As temperature is increased, the thermal excitation of spin waves reduces a ferromagnet's spontaneous magnetization. The energies of spin waves are typically only μeV in keeping with typical Curie points at room temperature and below.

In mesoscopic physics, ballistic conduction is the unimpeded flow of charge carriers, or energy-carrying particles, over relatively long distances in a material. In general, the resistivity of a material exists because an electron, while moving inside a medium, is scattered by impurities, defects, thermal fluctuations of ions in a crystalline solid, or, generally, by any freely-moving atom/molecule composing a gas or liquid. Without scattering, electrons simply obey Newton's second law of motion at non-relativistic speeds.

In physics, thermal contact conductance is the study of heat conduction between solid or liquid bodies in thermal contact. The thermal contact conductance coefficient, , is a property indicating the thermal conductivity, or ability to conduct heat, between two bodies in contact. The inverse of this property is termed thermal contact resistance.

Helium atom scattering (HAS) is a surface analysis technique used in materials science. It provides information about the surface structure and lattice dynamics of a material by measuring the diffracted atoms from a monochromatic helium beam incident on the sample.

<span class="mw-page-title-main">Surface phonon</span>

In solid state physics, a surface phonon is the quantum of a lattice vibration mode associated with a solid surface. Similar to the ordinary lattice vibrations in a bulk solid, the nature of surface vibrations depends on details of periodicity and symmetry of a crystal structure. Surface vibrations are however distinct from the bulk vibrations, as they arise from the abrupt termination of a crystal structure at the surface of a solid. Knowledge of surface phonon dispersion gives important information related to the amount of surface relaxation, the existence and distance between an adsorbate and the surface, and information regarding presence, quantity, and type of defects existing on the surface.

There are a number of possible ways to measure thermal conductivity, each of them suitable for a limited range of materials, depending on the thermal properties and the medium temperature. Three classes of methods exist to measure the thermal conductivity of a sample: steady-state, time-domain, and frequency-domain methods.

As the devices continue to shrink further into the sub-100 nm range following the trend predicted by Moore’s law, the topic of thermal properties and transport in such nanoscale devices becomes increasingly important. Display of great potential by nanostructures for thermoelectric applications also motivates the studies of thermal transport in such devices. These fields, however, generate two contradictory demands: high thermal conductivity to deal with heating issues in sub-100 nm devices and low thermal conductivity for thermoelectric applications. These issues can be addressed with phonon engineering, once nanoscale thermal behaviors have been studied and understood.

Phonons can scatter through several mechanisms as they travel through the material. These scattering mechanisms are: Umklapp phonon-phonon scattering, phonon-impurity scattering, phonon-electron scattering, and phonon-boundary scattering. Each scattering mechanism can be characterised by a relaxation rate 1/ which is the inverse of the corresponding relaxation time.

In heat transfer, thermal engineering, and thermodynamics, thermal conductance and thermal resistance are fundamental concepts that describe the ability of materials or systems to conduct heat and the opposition they offer to the heat current. The ability to manipulate these properties allows engineers to control temperature gradient, prevent thermal shock, and maximize the efficiency of thermal systems. Furthermore, these principles find applications in a multitude of fields, including materials science, mechanical engineering, electronics, and energy management. Knowledge of these principles is crucial in various scientific, engineering, and everyday applications, from designing efficient temperature control, thermal insulation, and thermal management in industrial processes to optimizing the performance of electronic devices.

For typical three-dimensional metals, the temperature-dependence of the electrical resistivity ρ(T) due to the scattering of electrons by acoustic phonons changes from a high-temperature regime in which ρ ∝ T to a low-temperature regime in which ρ ∝ T5 at a characteristic temperature known as the Debye temperature. For low density electron systems, however, the Fermi surface can be substantially smaller than the size of the Brillouin zone, and only a small fraction of acoustic phonons can scatter off electrons. This results in a new characteristic temperature known as the Bloch–Grüneisen temperature that is lower than the Debye temperature. The Bloch–Grüneisen temperature is defined as 2ħvskF/kB, where ħ is the Planck constant, vs is the velocity of sound, ħkF is the Fermi momentum, and kB is the Boltzmann constant.

Heat transfer physics describes the kinetics of energy storage, transport, and energy transformation by principal energy carriers: phonons, electrons, fluid particles, and photons. Heat is energy stored in temperature-dependent motion of particles including electrons, atomic nuclei, individual atoms, and molecules. Heat is transferred to and from matter by the principal energy carriers. The state of energy stored within matter, or transported by the carriers, is described by a combination of classical and quantum statistical mechanics. The energy is different made (converted) among various carriers. The heat transfer processes are governed by the rates at which various related physical phenomena occur, such as the rate of particle collisions in classical mechanics. These various states and kinetics determine the heat transfer, i.e., the net rate of energy storage or transport. Governing these process from the atomic level to macroscale are the laws of thermodynamics, including conservation of energy.

References

  1. ↑ Giri, A.; Hopkins, P.E. (2020). "A Review of Experimental and Computational Advances in Thermal Boundary Conductance and Nanoscale Thermal Transport across Solid Interfaces". Advanced Functional Materials. 30 (8): 1903857. doi:10.1002/adfm.201903857. S2CID   202037103.
  2. ↑ Hu, Ming; Keblinski, Pawel; Wang, Jian-Sheng; Raravikar, Nachiket (2008). "Interfacial thermal conductance between silicon and a vertical carbon nanotube". Journal of Applied Physics. 104 (8): 083503–083503–4. Bibcode:2008JAP...104h3503H. doi:10.1063/1.3000441.
  3. ↑ Hu, Han; Sun, Ying (2012). "Effect of nanopatterns on Kapitza resistance at a water-gold interface during boiling: A molecular dynamics study". Journal of Applied Physics. AIP Publishing. 112 (5): 053508–053508–6. Bibcode:2012JAP...112e3508H. doi:10.1063/1.4749393.
  4. 1 2 Stanley, Christopher M. (May 2018). Heat flow and Kapitza resistance across a Si|SiO2 interface: A first principles study (Thesis). hdl: 2346/73843 . OCLC   1045421642.
  5. 1 2 Gordiz, Kiarash; Henry, Asegun (16 March 2016). "Phonon Transport at Crystalline Si/Ge Interfaces: The Role of Interfacial Modes of Vibration". Scientific Reports. 6 (1): 23139. Bibcode:2016NatSR...623139G. doi: 10.1038/srep23139 . PMC   4793224 . PMID   26979787. S2CID   8893499.
  6. 1 2 Stanley, Christopher M.; Rader, Benjamin K.; Laster, Braxton H. D.; Servati, Mahsa; Estreicher, Stefan K. (December 2021). "The Role of Interface Vibrational Modes in Thermal Boundary Resistance". Physica Status Solidi A. 218 (23): 2100111. Bibcode:2021PSSAR.21800111S. doi:10.1002/pssa.202100111. S2CID   236408006.
  7. ↑ Monachon, Christian; Weber, Ludger; Dames, Chris (1 July 2016). "Thermal Boundary Conductance: A Materials Science Perspective". Annual Review of Materials Research. 46 (1): 433–463. Bibcode:2016AnRMS..46..433M. doi: 10.1146/annurev-matsci-070115-031719 .
  8. ↑ Swartz, Eric Thomas (1987). Solid-Solid Thermal Boundary Resistance (Thesis). Bibcode:1987PhDT.......170S. ProQuest   303581612.
  9. 1 2 3 Swartz, E. T.; Pohl, R. O. (1 July 1989). "Thermal boundary resistance". Reviews of Modern Physics. 61 (3): 605–668. Bibcode:1989RvMP...61..605S. doi:10.1103/revmodphys.61.605.
  10. ↑ Zeng, Taofang; Chen, Gang (1 April 2001). "Phonon Heat Conduction in Thin Films: Impacts of Thermal Boundary Resistance and Internal Heat Generation". Journal of Heat Transfer. 123 (2): 340–347. doi:10.1115/1.1351169.
  11. 1 2 Kapitza, P. L. (15 August 1941). "Heat Transfer and Superfluidity of Helium II". Physical Review. 60 (4): 354–355. Bibcode:1941PhRv...60..354K. doi:10.1103/PhysRev.60.354.
  12. ↑ Zu, H.; Dai, W.; de Waele, A.T.A.M. (2022). "Development of dilution refrigerators – A review". Cryogenics. 121. doi: 10.1016/j.cryogenics.2021.103390 . S2CID   244005391.
  13. ↑ Khater, A.; Szeftel, J. (1987-05-01). "Theory of the Kapitza resistance". Physical Review B. American Physical Society (APS). 35 (13): 6749–6755. Bibcode:1987PhRvB..35.6749K. doi:10.1103/physrevb.35.6749. PMID   9940924.
  14. 1 2 Lyeo, Ho-Ki; Cahill, David G. (2006-04-03). "Thermal conductance of interfaces between highly dissimilar materials". Physical Review B. American Physical Society (APS). 73 (14): 144301. Bibcode:2006PhRvB..73n4301L. doi:10.1103/physrevb.73.144301.
  15. 1 2 Costescu, Ruxandra M.; Wall, Marcel A.; Cahill, David G. (2003-02-27). "Thermal conductance of epitaxial interfaces". Physical Review B. American Physical Society (APS). 67 (5): 054302. Bibcode:2003PhRvB..67e4302C. doi:10.1103/physrevb.67.054302.
  16. 1 2 Gundrum, Bryan C.; Cahill, David G.; Averback, Robert S. (2005-12-30). "Thermal conductance of metal-metal interfaces". Physical Review B. American Physical Society (APS). 72 (24): 245426. Bibcode:2005PhRvB..72x5426G. doi:10.1103/physrevb.72.245426.
  17. ↑ Zhong, Hongliang; Lukes, Jennifer R. (2006-09-01). "Interfacial thermal resistance between carbon nanotubes: Molecular dynamics simulations and analytical thermal modeling". Physical Review B. American Physical Society (APS). 74 (12): 125403. Bibcode:2006PhRvB..74l5403Z. doi:10.1103/physrevb.74.125403.
  18. ↑ Estrada, David; Pop, Eric (2011-02-14). "Imaging dissipation and hot spots in carbon nanotube network transistors". Applied Physics Letters. AIP Publishing. 98 (7): 073102. arXiv: 1011.4551 . Bibcode:2011ApPhL..98g3102E. doi:10.1063/1.3549297. S2CID   40975236.
This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.