Voir la source de Équation d'état de Murnaghan

{{En-tête label|BA|année=2010}}
{{Voir homonymes|Murnaghan}}
{{Confusion|Équation d'état de Birch-Murnaghan}}

L''''équation d'état de Murnaghan''' est une relation qui lie le volume d'un corps et la pression à laquelle il est soumis.  C'est une des nombreuses [[Équation d'état|équations d'état]] qui ont été utilisées en [[sciences de la Terre]] pour [[Modèle mathématique|modéliser]] le comportement de la matière dans les conditions de hautes pressions qui règnent à l'intérieur du globe terrestre. Elle doit son nom à [[Francis Dominic Murnaghan]]<ref name="Murnaghan1944">{{Article | auteur = F.D. Murnaghan | titre = The Compressibility of Media under Extreme Pressures | journal = Proceedings of the National Academy of Sciences of the United States of America | vol = 30 | pages = 244-247 | date = 1944 | lang = en | texte = http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=1078704}}.</ref> qui l'a proposée en 1944 afin de rendre compte sur une gamme de pressions aussi large que possible d'un fait expérimentalement établi : plus on comprime un solide, plus il devient difficile de le comprimer.

L'équation de Murnaghan est déduite, moyennant certaines hypothèses, des équations de la [[mécanique des milieux continus]]. Elle fait intervenir deux paramètres ajustables qu'on identifie au [[Module d'élasticité isostatique|module d'incompressibilité]] <math>K_0</math> et à sa [[dérivée]] première par rapport à la pression, <math>K'_0</math>, tous deux pris à pression nulle. En général, on détermine ces deux coefficients par une [[Régression (statistiques)|régression]] sur les valeurs du volume <math>V</math> en fonction de la pression <math>P</math> obtenues expérimentalement, le plus souvent par [[Diffractométrie de rayons X|diffraction des rayons X]]. La régression peut également être effectuée sur les valeurs de l'énergie en fonction du volume obtenues par [[Méthode ab initio de chimie quantique|calcul ''ab initio'']].

Tant que la réduction de volume reste faible, c'est-à-dire pour <math>V/V_0</math> supérieur à 90 % environ, l'équation de Murnaghan permet de modéliser avec une précision satisfaisante les données expérimentales. De plus, à la différence de la plupart des nombreuses équations d'état proposées, elle donne une expression explicite du volume en fonction de la pression <math>V(P)</math>. Mais son domaine de validité reste trop limité et son interprétation physique insatisfaisante. On a plus souvent recours, pour l'analyse des données de compression, à d'autres équations d'état plus élaborées dont la plus utilisée est l'[[équation d'état de Birch-Murnaghan]].

== Contexte général ==

L'étude de la [[structure interne de la Terre]] passe par la connaissance des propriétés mécaniques des constituants des couches internes de la planète, dans lesquelles se rencontrent des conditions extrêmes : les pressions s'y comptent en centaines de [[Pascal (unité)|gigapascals]] et les températures en milliers de degrés. L'étude des propriétés de la matière dans ces conditions peut se faire de manière expérimentale grâce à des dispositifs comme la [[cellule à enclumes de diamant]] pour les pressions statiques, ou en soumettant la matière à des [[Onde de choc|ondes de choc]]. Elle a également donné lieu à des travaux théoriques visant à déterminer des [[équations d'état]], c'est-à-dire des relations liant les différents paramètres qui définissent dans ce cas l'état de la matière : le volume (ou la densité), la température et la pression.

On peut distinguer deux approches :
* les équations d'état dérivées de [[Potentiel interatomique|potentiels interatomiques]], ou éventuellement de [[Méthode ab initio de chimie quantique|calculs ''ab initio'']] ;
* les équations d'état dérivées des relations générales de la [[mécanique (science)|mécanique]] et de la [[thermodynamique]]. L'équation de Murnaghan appartient à cette seconde catégorie.

Plusieurs dizaines d'équations ont ainsi été proposées par différents auteurs<ref>{{Article |auteur = P.T. Wedepohl |titre = Comparison of a simple two-parameter equation of state with the Murnaghan equation |revue = Solid State Communications |vol = 10 |pages = 947-951 |date = 1972 |lang = en |texte = https://dx.doi.org/10.1016/0038-1098(72)90228-1}}.</ref>. Ce sont des relations empiriques, dont la qualité et la pertinence dépendent de l'usage qui en est fait. On peut en juger selon différents critères : le nombre de paramètres indépendants qu'elles font intervenir, la signification physique qu'il est possible d'attribuer à ces paramètres, la qualité des affinements qu'elles permettent sur les données expérimentales, la cohérence des hypothèses théoriques qui les sous-tendent, leur capacité à extrapoler le comportement des solides aux fortes compressions<ref name="Stacey1981" />.

== Expressions de l'équation d'état ==

De manière générale, à température constante, le [[Module d'élasticité isostatique|module d'incompressibilité]] est défini par :
:<math> K = -V \left(\frac{\partial P}{\partial V}\right)_T</math>.
La manière la plus simple d'obtenir une équation d'état liant <math>P</math> et <math>V</math> consiste à supposer que <math>K</math> est constant, c'est-à-dire indépendant de la pression et de la déformation du solide ; on retrouve alors simplement la [[loi de Hooke]]. Dans cette hypothèse, le volume décroît exponentiellement avec la pression. Ce n'est pas un résultat satisfaisant car il est expérimentalement établi que plus on comprime un solide, plus il devient difficile de le comprimer. Pour aller plus loin, il faut donc prendre en compte les variations des propriétés élastiques du solide avec la compression.

L'hypothèse de Murnaghan est de supposer que le module d'incompressibilité est une fonction linéaire de la pression<ref name="Murnaghan1944"/> :
:<math> K = K_0 + P\ K_0'</math>
L'équation de Murnaghan est le résultat de l'intégration de cette [[équation différentielle]] : 
:<math>
P(V) = \frac{K_0}{K_0'} \left[\left(\frac{V}{V_0}\right)^{-K_0'} - 1\right]
</math>
On peut également exprimer le volume en fonction de la pression : 
:<math>
V(P) = V_0 \left[1+ P \left(\frac{K'_0}{K_0}\right)\right]^{-1/K'_0}
</math>

Cette présentation simplifiée est toutefois dénoncée par Poirier comme manquant de rigueur<ref>{{Harvsp|réf = Poirier2002 |texte = Poirier, 2002 |p = 65}}.</ref>. La même relation peut être démontrée de manière différente en partant du fait que le produit du module d'incompressibilité par le [[Dilatation thermique|coefficient de dilatation thermique]] ne dépend pas de la pression pour un matériau donné<ref name="Kumar1995">{{Article |auteur = M. Kumar |titre = High pressure equation of state for solids |journal = Physica B: Condensed Matter |vol = 212 |pages = 391-394 |date = 1995 |lang = en |texte = https://dx.doi.org/10.1016/0921-4526(95)00361-C}}.</ref>.

Dans certaines circonstances, en particulier en lien avec les calculs ''ab initio'', on préfèrera l'expression de l'énergie en fonction du volume<ref>{{Harvsp |réf = Silvi1997 |texte = Silvi, 1997 |p = 122}}.</ref>, qui peut s'obtenir en intégrant l'équation précédente selon la relation <math>P = -\mathrm dE/\mathrm dV</math>. Elle peut s'écrire, pour <math>K'_0</math> différent de 1,
:<math>
E(V) = E_0 + K_0\,V_0\left[\frac{1}{K_0'(K_0'-1)}\left(\frac{V}{V_0}\right)^{1-K_0'} + \frac{1}{K_0'}\frac{V}{V_0} - \frac{1}{K_0'-1}\right].
</math>

== Avantages et limites ==

En dépit de sa simplicité, l'équation de Murnaghan est en mesure de reproduire les données expérimentales dans une gamme de pression qui peut être assez large, de l'ordre de <math>K_0/2</math> selon Anderson<ref name="Anderson1995">{{Harvsp | réf = Anderson1995 |texte = Anderson, 1995 |p = 179}}.</ref>. On peut également retenir qu'elle reste satisfaisante tant que le rapport <math>V/V_0</math> reste supérieur à 90 % environ<ref>{{Ouvrage | langue=en | auteur1=R.J. Angel | titre=High-Pressure Crystallography | éditeur= | année= | passage=21-36 | isbn= | titre chapitre=Some pratical aspects of studying equations of state and structural phase transitions at high pressure}}.</ref>. Dans cette gamme, l'équation de Murnaghan présente un avantage par rapport aux autres équations d'état si l'on souhaite pouvoir exprimer le volume en fonction de la pression<ref name="Holzapfel1996">{{Article | auteur = W.B. Holzapfel | titre = Physics of solids under strong compression | revue = Reports on Progress in Physics | vol = 59 | pages = 29-90 | date = 1996 | lang = en | texte = http://www.iop.org/EJ/abstract/0034-4885/59/1/002 |id = Holzapfel1996}}.</ref>.

Néanmoins, plusieurs exemples ont montré que d'autres équations permettaient d'obtenir des affinements de meilleure qualité<ref name="Kumar1995" />. De plus, plusieurs arguments théoriques et expérimentaux montrent que l'équation de Murnaghan est insatisfaisante. Ainsi, dans la limite où le rapport <math>V/V_0</math> devient très faible, la théorie prédit que <math>K'</math> tend vers <math>5/3</math> ; c'est la [[limite de Thomas-Fermi]]<ref>La théorie de Thomas-Fermi considère schématiquement le solide très fortement comprimé comme un gaz d'électrons dégénéré ([[gaz de Fermi]]) auquel on ajoute un terme d'[[écrantage]] pour tenir compte de la présence des noyaux atomiques.</ref>{{,}}<ref name="Holzapfel1996" />. Or, dans l'équation de Murnaghan, <math>K'</math> est constant et fixé à sa valeur initiale. Sauf à poser systématiquement <math>K'_0=5/3</math>, on a l'assurance que l'équation de Murnaghan devient alors incompatible avec la théorie. De fait, lorsqu'il est extrapolé, le comportement prédit par l'équation de Murnaghan devient assez rapidement invraisemblable<ref name="Holzapfel1996" />.

Indépendamment de cet argument théorique, l'expérience montre clairement que <math>K'</math> décroît avec la pression, c'est-à-dire que la [[dérivée seconde]] du module d'incompressibilité <math>K''</math> est strictement négative. Un développement à l'ordre 2 basé sur le même principe (cf. [[#Extensions et généralisations|section suivante]]) permet de tenir compte de cette observation, mais cette solution reste encore insatisfaisante. En effet, elle conduit à un module d'incompressibilité qui devient négatif dans la limite où la pression tend vers l'infini. De fait, c'est une contradiction inévitable quel que soit le développement polynomial choisi, car on aura toujours un terme dominant, divergeant à l'infini<ref name="Stacey1981">{{Article |auteur = F.D. Stacey, B.J. Brennan et R.D. Irvine |titre = Finite strain theories and comparison with seismological data |revue = Surveys in Geophysics |vol = 4 |pages = 189-232 |date = 1981 |lang = en |texte = http://www.springerlink.com.chimie.gate.inist.fr/content/x0254w7q10k44w52/?p=8afe65621b9d4a79a3f0d9fd2466b263&pi=1}}.</ref>.

Ces limites importantes ont conduit à l'abandon de l'équation de Murnaghan, que W. Holzapfel qualifie d'{{citation|expression mathématique utile, mais dépourvue de justification physique}}<ref>{{Harvsp |texte = Holzapfel (2001) |réf = Holzapfel2001}}.</ref>. Dans la pratique, l'analyse des données de compression se fait à l'aide d'équations d'état plus élaborées. La plus utilisée au sein de la communauté des sciences de la terre est l'[[équation d'état de Birch-Murnaghan]], au deuxième ou troisième ordre selon la [[qualité des données]] collectées<ref>{{Ouvrage |langue=en |nom1=E. Boldyreva et P. Dera |directeur1=X |nom2=T. Boffa Ballaran |titre=High-Pressure Crystallography : From Fondamental Phenomena to Technological Applications |éditeur=Springer |année= |passage=135-145 |isbn= |titre chapitre=Equations of state and their applications in geosciences}}.</ref>.

Enfin, une limitation tout à fait générale de ce type d'équation d'état est leur incapacité à prendre en compte les [[Transition de phase|transitions de phases]] induites par la pression et la température : la fusion, mais aussi les multiples transitions solide-solide qui peuvent provoquer des changements abrupts de la densité et du module d'incompressibilité en fonction de la pression<ref name="Stacey1981" />.

== Exemples d'utilisation ==

Dans la pratique, l'équation de Murnaghan est utilisée pour effectuer une [[Régression (statistiques)|régression]] sur un [[jeu de données]], d'où on tire les valeurs des coefficients <math>K_0</math> et <math>K'_0</math>. Ces coefficients obtenus, et connaissant la valeur du volume à pression nulle, on est alors en principe capable de calculer le volume, la densité et le module d'incompressibilité pour n'importe quelle pression.

Le jeu de données est la plupart du temps une série de mesures du volume pour différentes valeurs de pression appliquée, obtenue la plupart du temps par [[Diffractométrie de rayons X|diffraction des rayons X]]. Il est également possible de travailler sur des données théoriques, en calculant l'énergie pour différentes valeurs du volume par des [[Méthode ab initio de chimie quantique|méthodes ''ab initio'']] puis en effectuant la régression sur ces résultats. On obtient alors une valeur théorique du module d'élasticité qui peut être comparée à un résultat expérimental.

Le tableau qui suit répertorie quelques résultats obtenus sur différents matériaux, avec pour seul objectif d'illustrer de quelques valeurs numériques les analyses qui ont pu être faites à l'aide de l'équation de Murnaghan, sans préjuger de la qualité des modèles obtenus. Étant donné les critiques qui ont été formulées dans la section précédente sur la signification physique de l'équation de Murnaghan, ces résultats doivent être considérés avec précaution.

{| class="wikitable centre"
! scope=col | Composé
! scope=col | <math>K_0</math> (GPa)
! scope=col | <math>K'_0</math>
|- 
|[[Fluorure de sodium|NaF]]<ref name="Kumar1995" />
|46,5
|5,28
|- 
|[[chlorure de sodium|NaCl]]<ref name="Kumar1995" />
|24,0
|5,39
|- 
|[[Bromure de sodium|NaBr]]<ref name="Kumar1995" />
|19,9
|5,46
|- 
|[[Iodure de sodium|NaI]]<ref name="Kumar1995" />
|15,1
|5,59
|- 
|[[Oxyde de magnésium|MgO]]<ref name="Anderson1995" />
|156
|4,7
|- 
|[[Calcite]] (CaCO<sub>3</sub>)<ref name="Silvi1997">{{Harvsp |id = Silvi1997 |texte = Silvi, 1997 |p = 123}}.</ref>
|75,27
|4,63
|- 
|[[Magnésite]] (MgCO<sub>3</sub>)<ref name="Silvi1997" />
|124,73
|3,08
|- 
|[[Carbure de silicium]] (3C-SiC)<ref>{{Article|auteur = K. Strössner, M. Cardona and W. J. Choyke |titre = High pressure X-ray investigations on 3C-SiC|journal = Solid State Communications |vol = 63 |pages = 113-114 |date = 1987 |doi = 10.1016/0038-1098(87)91176-8 |lang = en}}.</ref>
|248
|4,0
|}

== Extensions et généralisations ==

Afin d'améliorer les modèles ou échapper aux critiques signalées plus haut, plusieurs généralisations de l'équation de Murnaghan ont été proposées. Elles consistent habituellement en l'abandon d'une hypothèse simplificatrice et l'ajout d'un paramètre ajustable supplémentaire. Ceci peut améliorer les qualités des affinements mais complique les expressions. La question de la signification physique de ces paramètres supplémentaires reste également posée.

Une première sophistication consiste à ajouter un terme en <math>P^2</math> au développement précédent<ref>{{Article |auteur = J.R. MacDonald et D.R. Powell |titre= ''Discrimination Between Equations of State'' |revue = Journal of Research of the National Bureau of Standards, Section A |vol = 75 |pages = 441 |date = 1971 |lang = en}}.</ref>{{,}}<ref>{{Harvsp|réf = MacDonald1969 |texte = MacDonald, 1969 |p = 320}}.</ref> en posant <math> K = K_0 + PK_0' + P^2K_0''</math>. La résolution de cette équation différentielle donne l'équation de Murnaghan au second-ordre :
:<math>
P(V) = 2 \frac{K_0}{K_0'} \left[\frac{\Gamma}{K_0'}\,\frac{\left(\frac{V_0}{V}\right)^{\Gamma}+1}{\left(\frac{V_0}{V}\right)^{\Gamma}-1} - 1\right]^{-1}
</math>
pour <math>\Gamma^2 = {K_0'}^2 - 2 K_0 K_0'' > 0</math>. On retrouve naturellement l'équation au premier ordre en prenant <math>K_0''=0</math>. Des développements à un ordre supérieur à 2 sont possibles en principe<ref>{{Article | auteur = Kazuhiro Fuchizaki | titre = Murnaghan equation of state revisited | journal = Journal of the Physical Society of Japan | vol = 75 | date = 2006 | pages = 034601 | lang = en | texte = http://jpsj.ipap.jp/link?JPSJ/75/034601/}}.</ref> mais au prix de l'ajout d'un paramètre ajustable pour chaque terme.

D'autres généralisations peuvent également être citées :
* Kumari et Dass ont proposé une généralisation en abandonnant la condition <math>K''=0</math> mais en supposant le rapport <math>K''/K'</math> indépendant de la pression<ref>{{Article |auteur = M. Kumari et N. Dass |titre = An equation of state applied to sodium chloride and caesium chloride at high pressures and high temperatures |journal = Journal of Physics: Condensed Matter |vol = 2 |pages = 3219 |date = 1990 |lang = en |texte = http://iopscience.iop.org/0953-8984/2/14/006}}.</ref> ;
* Kumar a proposé une généralisation en prenant en compte une dépendance du paramètre d'Anderson en fonction du volume<ref name="Kumar1995" />. Il a été montré par la suite que cette équation généralisée n'était pas nouvelle mais se ramenait à l'[[équation de Tait]]<ref>{{Article |auteur = J. Shanker, B. Singh and S.S. Kushwah |titre = On the high-pressure equation of state for solids |journal = Physica B: Condensed Matter |vol = 229 |pages = 419-420 |date = 1997 |lang = en |texte = https://dx.doi.org/10.1016/S0921-4526(96)00528-5}}.</ref>.

== Notes et références ==
{{Références}}

== Annexes ==
=== Bibliographie ===
==== Ouvrages ====
* {{Ouvrage |langue=en |auteur1=J.P. Poirier |titre=Introduction to the physics of the Earth's interior |éditeur=[[Cambridge University Press]] |année=2002 |isbn= |lire en ligne=https://books.google.fr/books?id=_YOCl86XPHoC&pg=PP1 |id=Poirier2002}}.
* {{Ouvrage |langue=en |auteur1=B. Silvi et P. d'Arco |titre=Modelling of Minerals and Silicated Materials |éditeur=Kluwer Academic Publishers |année=1997 |isbn= |lire en ligne=https://books.google.fr/books?id=60wSAQAAIAAJ |id=Silvi1997}}.
* {{Ouvrage |langue=en |auteur1=O.L. Anderson |titre=Equations of state of solids for geophysics and ceramic science |éditeur=[[Oxford University Press]] |année=1995 |isbn= |lire en ligne=https://books.google.fr/books?id=qSqiz2NK7TIC&pg=PR4 |id=Anderson1995}}.

==== Articles de revue ====
* {{Article | auteur = W.B. Holzapfel | titre = Equations of state for solids under strong compression | revue = Zeitschrift für Kristallografie | vol = 216 | pages = 473-488 | date = 2001 | lang = en | texte = http://www.oldenbourg-link.com/doi/abs/10.1524/zkri.216.9.473.20346 |id = Holzapfel2001}}.
* {{Article |auteur = J.R. MacDonald |titre = Review of Some Experimental and Analytical Equations of State |revue = Review of Modern Physics |vol = 41 |pages = 316-349 |date = 1969 |lang = en |texte = http://rmp.aps.org/abstract/RMP/v41/i2/p316_1 |id = MacDonald1969}}.

=== Articles connexes ===
* [[Équation d'état]]
* [[Équation d'état de Birch-Murnaghan]]
* [[Équation d'état de Vinet]]

=== Liens externes ===
* [http://www.crystal.vt.edu/crystal/software.html EosFit], un programme permettant l'affinement de données expérimentales et le calcul de relations <math>P(V)</math> pour différentes équations d'état, dont l'équation de Murnaghan.

{{Portail|Physique}}
{{Bon article|vote=BA|oldid=59439976|date=25 novembre 2010}}

{{DEFAULTSORT:Équation d'etat de Murnaghan}}
[[Catégorie:Équation d'état|Murnagan]]