Résistance de Kapitza
La résistance de Kapitza est la résistance à un flux de chaleur par une interface géométriquement parfaite entre deux matériaux. Elle est le résultat de la diffusion de phonons due au changement de structure atomique à cette interface. Elle a été mise en évidence par Piotr Kapitsa en 1941[1]. Si le phénomène a été découvert sur les interfaces hélium liquide-solide il est présent sur une interface solide-solide.
Résistance thermique d'interface
La résistance thermique d'interface mesure la discontinuité de température à une interface entre deux matériaux. Elle est définie par l'expression où est la densité de flux de chaleur. est la conductance d'interface.
Cette quantité est d'importance majeure en microélectronique où l'on atteint des densités de flux de l'ordre de 103 MW/m2 nécessitant des conductances d'interface supérieures à 102 MW/m2/K pour pouvoir assurer un refroidissement efficace[5].
Phénoménologie
Le phénomène peut être décrit par deux voies :
- La dynamique du réseau atomique qui consiste à simuler un réseau atomique simple (cubique par exemple) partagé en deux parties dont les atomes ont des propriétés différentes[6]. Il s'agit là de la généralisation de l'approche classique de description des phonons conduisant à un système linéaire qui, dans ce cas, est résolu numériquement. La représentativité du phénomène est jugée suffisante pour mettre en évidence les phénomènes malgré ses défauts, en particulier la juxtaposition des réseaux aux interfaces sans réarrangement local ;
- La dynamique moléculaire, beaucoup plus lourde à mettre en œuvre, permet de mieux représenter l'interface de deux matériaux, solides, liquides ou gaz, à l'échelle atomique et d'expliquer les mécanismes mis en jeu à cette échelle[7]. Dans le cas d'un liquide on y voit de notables modifications des propriétés vibrationnelles des atomes ou molécules lorsqu'on se rapproche des interfaces, typiquement sur une longueur de 4 à 5 fois la dimension de l'atome ou de la molécule[8].
Dans les deux approches il est possible de générer un phonon unique et de mesurer son comportement à l'interface[6]. La probabilité de passage d'un milieu à l'autre pour ce phonon dépend du nombre d'onde, de la polarisation et du sens de propagation, la différence entre les sens de propagation étant plus marquée pour les phonons acoustiques. Toutefois il convient de prendre en compte les populations de phonons dans chaque milieu, très différentes suivant la nature du cristal à une même fréquence, ainsi que l'incompatibilité des bandes de propagation. Ce dernier point est outrepassé par des conversions acoustique-optique de phonons.
Modélisation cinétique
Dans la description cinétique donnée par l'équation de Boltzmann-Peierls la densité de flux de chaleur entre les milieux 1 de température et 2 de température est donnée par :
où est le nombre de phonons de nombre d'onde k, leur énergie et c la vitesse de groupe assimilée à la vitesse du son. est la probabilité de passage de 1 vers 2.
s'exprime de façon identique et la densité de flux de chaleur résultante est , la conductance .
Il existe une relation entre les probabilités de passage liée à la microréversibilité. Un matériau à l'équilibre thermodynamique a une population de phonons donnée par la statistique de Bose-Einstein. Il émet une densité de flux égale à où est l'équivalent de la constante de Stefan-Boltzmann et la quantité traverse l'interface. La densité de flux totale est donc . La réversibilité impose pour d'où la relation entre coefficients :
- .
Tout le problème est contenu dans la connaissance de l'une des quantités , modélisée de différentes façons qui toutes ignorent les possibles transitions en fréquence[9].
Approximation acoustique
Dans cette approximation due à Isaak Khalatnikov (1952)[10] le phonon est traité comme une onde plane analogue au photon arrivant sur une interface et étant réfléchi ou transmis. Dans ce dernier cas on suppose vérifiées les lois de Snell-Descartes. Cela suppose de connaître la distribution angulaire des phonons ou de supposer l'isotropie.
Par souci de simplicité le coefficient de transmission est parfois calculé en utilisant une impédance acoustique où est la masse volumique du matériau :
Cette expression ne vérifie pas la relation de réversibilité.
Si de plus on utilise le modèle de Debye pour le calcul de , ce qui suppose une propagation isotrope des phonons dans un milieu lui-même isotrope, on obtient une dépendance en pour la conductance thermique thermique[9] :
où j est l'indice correspond à la polarisation.
Comparé à l'expérience, ce modèle surestime la résistance d'interface.
Approximation diffusive
Afin d'augmenter la probabilité de transmission l'idée introduite par E. T. Swartz (1987) est de remplacer la réflexion et la transmission à l'interface par un phénomène de diffusion élastique. Le phénomène de diffusion est mesurable sur une interface avec le vide pour les fréquences hautes (>100 GHz). Il est connu pour sa sensibilité à l'état de surface[11]. Il existe des modèles décrivant ce phénomène [12].
Pour une diffusion isotrope la probabilité de transmission d'un phonon est de 0,5. Ce chiffre est valable à l'interface entre deux matériaux identiques, alors qu'il est de 1 pour l'approximation acoustique.
La diffusion donne des résistances d'interface sous-évaluées.
Le choix du modèle
Un modèle plus général prenant en compte une diffusion inélastique anisotrope se heurte à difficulté de connaître l'interface à l'échelle atomique. De plus, si l'écart entre approximation acoustique et diffusion est important sur une interface solide-liquide, il est faible pour une interface solide-solide. D'une façon générale on se contente de répartir le phénomène : une part pour la diffusion, l'autre pour la réflexion/transmission.
Références
- (ru) Piotr Kapitsa, « The Study of Heat Transfer in Helium II », Journal of Experimental and Theoretical Physics, vol. 11,
- (en) Ho-Ki Lyeo et David G. Cahill, « Thermal conductance of interfaces between highly dissimilar materials », Physical Review B, vol. 73, no 14, , p. 144301
- (en) Ruxandra M. Costescu, Marcel A. Wall et David G. Cahill, « Thermal conductance of epitaxial interfaces », Physical Review B, vol. 67, no 5, , p. 054302
- (en) Bryan C. Gundrum, David G. Cahill et Robert S. Averback, « Thermal conductance of metal-metal interfaces », Physical Review B, vol. 72, no 24, , p. 245426
- (en) Ling Hu, Pawel Keblinski, Jian-Sheng Wang et Nachiket Raravikar, « Interfacial thermal conductance between silicon and a vertical carbon nanotube », Journal of Applied Physics, vol. 104, no 8, , p. 083503
- (en) Xiabo Li et Ronggui Yang, « Size-dependant phonon transmission across dissimilar material interfaces », Journal of Physics: Condensed Matter, vol. 24, , p. 155302
- (en) Zhi Liang et Ming Hu, « Tutorial: Determination of thermal boundary resistance by molecular dynamics simulations », Journal of Applied Physics, no 123, , p. 191101
- (en) Xiao Lu, Donatas Surblys, Yoshiaki Kawagoe, Abdul Rafeq Bin Salemen, Hiroki Matsubara, Gota Kikugawa et Taku Ohara, « A molecular dynamics study of thermal boundary resistance over solid interfaces with an extremely thin liquid film », International Journal of Heat and Mass Transfer, vol. 147, , p. 118949 (lire en ligne)
- (en) E. T. Swartz et R. O. Pohl, « Thermal boundary resistance », Reviews of Modern Physics, vol. 61, no 3,
- (en) I. M. Khalatnikov, An introduction to the theory of superfluidity, CRC Press, (lire en ligne)
- (en) W. Eisenmenger, Phonon Scattering at Surfaces and Interfaces, vol. 68, Springer Series in Solid-State Sciences book series, , 204-211 p.
- (en) M. Adamenko et I. N., Fuks, « Roughness and Thermal Resistance of the Boundary Between a Solid and Liquid Helium », Journal of Experimental and Theoretical Physics, vol. 32, , p. 1123-1129