Potentiel interatomique

Le potentiel est constitué par la superposition de deux parties, l'une répulsive dominante à faible distance, l'autre possiblement attractive[2] - [1] - [3]. L'origine de ces forces est la suivante :

• partie répulsive : pour celle-ci les noyaux sont suffisamment proches pour que les nuages d'électrons se superposent. La théorie des perturbations est inapplicable et le calcul, de nature quantique[4], est celui d'une quasi-molécule, les forces impliquées étant les forces coulombiennes et l'interaction d'échange.
• partie attractive : la théorie des perturbations est applicable, déclinée sous diverses méthodes[2]. Dans la partie ${\displaystyle r<0.5\,nm}$ l'interaction d'échange reste importante. C'est dans cette région que se situe le minimum du potentiel. Dans la région externe l'opérateur d'interaction peut être développé en série multipolaire. On parle alors de forces de van der Waals. Le premier terme du développement correspond au dipôle mutuellement induit (forces de London). Il est en ${\displaystyle r^{-6}}$. On retrouve ce terme dans de nombreux potentiels : Lennard-Jones, Buckingham, etc. Les termes d'ordre supérieur correspondent aux interactions dipôle induit / quadrupôle induit sont en ${\displaystyle r^{-8}}$ et généralement négligés.

Il est possible de calculer ab initio ce genre de potentiel mais ceux utilisés sont généralement des approximations plus ou moins grossières de la réalité basées sur des considérations physiques mais utilisant des paramètres de réglage pour représenter au mieux une propriété donnée : équation d'état, viscosité, etc. On trouve ainsi des lois à un paramètre comme le potentiel sphères dures, à deux paramètres comme le potentiel de Lennard-Jones ou le potentiel de Sutherland, à trois paramètres comme le potentiel de Buckingham modifié ou le potentiel de Morse.

À mi-chemin entre le calcul exact et le simple modèle on trouve des potentiels tels que celui de Rydberg, Klein et Rees pour les molécules diatomiques, basé sur le recalage d'un calcul de type quantification semi-classique sur un spectre optique mesuré. Le résultat peut être recherché sous forme purement numérique ou par une série limitée de termes analytiques. Le potentiel de Morse est un cas particulier de ce développement.

La stabilité est liée à la profondeur du puits de potentiel ${\displaystyle D_{e}}$ qu'il faut vaincre pour séparer les atomes. Les gaz rares forment des molécules mais la faiblesse de la liaison induite (${\displaystyle D_{e}/k_{B}\simeq 7\,K}$ pour l'hélium[1]) fait que de telles molécules ne résistent pas à une collision, d'où une durée de vie très faible sauf aux basses températures.

Le fait que les paramètres des potentiels sont issus de l'expérience fait que l'on peut utiliser ce type d'approche pour des molécules : le potentiel est alors interprétable comme une certaine moyenne statistique des potentiels[1] pour toutes les configurations possibles.

Ce qui précède a naturellement des limites. Pour les molécules longues des potentiels plus complexes prenant en compte l'orientation relative des molécules ont été développés, très simples comme le potentiel de Taro Kihara basé sur l'interaction d'ellipsoïdes élastiques infiniment durs ou celui le potentiel de J. Corner où chaque molécule est définie par deux points[1].

Les deux molécules peuvent être polaires. Il existe alors une interaction due aux forces de Keesom (interaction dipôle / dipôle), une interaction mutuelle dipôle / dipôle induit (forces de Debye) et une interaction dipôle induit / dipôle induit (forces de London). L'ensemble est connu sous le nom de potentiel de van der Waals conduisant à des termes en ${\displaystyle r^{-6}}$. Les termes d'ordre supérieur sont généralement négligés.

D'une façon générale le lien entre potentiel et propriétés marche dans les deux sens, le potentiel étant souvent calculé pour caler un potentiel, les autres propriétés s'en déduisant.

• On peut effectuer ce lien pour les propriétés spectroscopiques des gaz, par exemple par la méthode de Rydberg-Klein-Rees.
• On peut relier le potentiel à l'équation d'état (liquide / gaz) dont un exemple est l'équation d'état de van der Waals.
• Les propriétés de transport gazeux (viscosités, coefficients de diffusion, conductivités) entrent dans ce même schéma en utilisant les intégrales de collision.

• Les propriétés thermodynamiques peuvent être calculées à partir d'un potentiel simple de type oscillateur harmonique par le modèle de Debye.
• Les propriétés mécaniques telles que coefficient de dilatation ou paramètre de Grüneisen sont liées à l'équation d'état, donc au potentiel et plus précisément à ses termes anharmoniques.
• A contrario la conductivité thermique qui est en principe calculable à partir des termes anharmoniques du potentiel est dans les faits liée à tous les défauts de constitution du cristal.