マーナハンの状態方程式のソースを表示

{{混同|[[バーチ・マーナハンの状態方程式]]}}
'''マーナハンの状態方程式'''（マーナハンのじょうたいほうていしき、{{lang-en-short|Murnaghan equation of state}}）は、物体の体積と加えられた圧力との関係を表す[[状態方程式 (熱力学)|状態方程式]]である。この方程式は、[[地球科学]]や衝撃物理学において、物質が高圧条件下でどのように振る舞うかをモデル化するために使用される多くの状態方程式の1つである。1944年に[[ジョンズ・ホプキンス大学]]の{{仮リンク|フランシス・ドミニク・マーナハン|en|Francis Dominic Murnaghan (mathematician)}}によって提案され<ref name="Murnaghan1944">{{Citation | last = F.D. |first=Murnaghan | title = The Compressibility of Media under Extreme Pressures | journal = Proceedings of the National Academy of Sciences of the United States of America | volume = 30 |issue=9 | pages = 244–247 | year = 1944 | pmc=1078704 | doi=10.1073/pnas.30.9.244 | pmid=16588651|bibcode=1944PNAS...30..244M |doi-access=free }}</ref>、固体が圧縮されるほどさらに圧縮することが難しくなるという実験的事実を反映している。

この方程式は、特定の仮定のもとで[[連続体力学]]の方程式から導出される。調整可能な2つのパラメータ、すなわち体積弾性率<math>K_0</math>と圧力に対するその一次導関数<math>K'_0</math>（いずれも常圧条件下で測定）が含まれる。これらの係数は通常、体積<math>V</math>と圧力<math>P</math>の関係を示す実験的なデータに基づいて[[回帰分析]]により決定される。これらのデータは、X線回折や衝撃試験によって得られるほか、[[ab initio]]計算や[[分子動力学法]]で得られる、体積に対するエネルギーの値を用いて回帰分析が行われることもある。

マーナハンの状態方程式は、通常次のように表される。

<math display="block">
P(V) = \frac{K_0}{K_0'} \left[\left(\frac{V}{V_0}\right)^{-K_0'} - 1\right] \,
</math>

体積の圧縮率が低い場合、つまり<math>V/V_0</math>が約90%以上である場合、この方程式は実験データを満足のいく精度でモデル化できる。また、他の多くの状態方程式と異なり、圧力の関数として体積<math>V(P)</math>を明示的に表している。しかし、この方程式の有効範囲は限定的であり、物理的解釈が不十分な場合がある。それでも、この方程式は固体爆薬のモデルで広く使用されている。さらに詳しい状態方程式として、地球物理学では[[バーチ・マーナハンの状態方程式]]が、金属や合金の衝撃物理学では[[ミー・グリュナイゼンの状態方程式]]が広く使用されている。

== 背景 ==
地球内部の構造を、内層を構成する物質の機械的特性を通じて研究するには、極端な条件を扱う必要がある。圧力は数百[[パスカル (単位)|ギガパスカル]]、温度は数千度に達することがある。このような条件下での物質特性の研究は、[[ダイヤモンドアンビルセル]]のような装置を用いて静的圧力をかけたり、[[衝撃波]]を物質に与えるといった実験的手法にてできる。また、物質の状態を定義するさまざまなパラメータ（体積、密度、温度、圧力）の関係式を決定する理論的な研究も行われた。

状態方程式を定めるためには、次の2つのアプローチがある。

* [[原子間ポテンシャル]]、または[[ab initio]]計算に基づくアプローチ

* 状態方程式力学および熱力学の一般関係に基づくアプローチ（マーナハンの状態方程式はこのカテゴリである。）

これまでに、多くの研究者によって数十種類の状態方程式が提案されてきた<ref>{{Citation |first = P.T.|last= Wedepohl |title = Comparison of a simple two-parameter equation of state with the Murnaghan equation |journal = Solid State Communications |volume = 10 |issue= 10 |pages = 947–951 |year = 1972 |doi=10.1016/0038-1098(72)90228-1|bibcode= 1972SSCom..10..947W }}</ref>。これらの方程式は経験的な関係式であり、その質や適用の有用性は、関係式に含まれるパラメータの数、パラメータに割り当てられる物理的な意味、実験データの質 、理論的仮定の一貫性、および高圧縮状態での固体の挙動を外挿できる能力といった基準で判断できる<ref name=Stacey1981/>。

== 状態方程式の導出 ==
一定温度での体積弾性率は、以下のように定義される。

<math display="block"> K = -V \left(\frac{\partial P}{\partial V}\right)_T</math>

ここで、圧力<math>P</math>と体積<math>V</math>を結び付ける最も簡単な方法は、体積弾性率<math>K</math>が一定であり、圧力や固体の変形に依存しないと仮定し、フックの法則を発見する。この場合、体積は圧力に対して指数関数的に減少する。しかし、実際には、固体が圧縮されるほど、さらに圧縮することが難しくなることが知られており、この結果は実験的に満足のいくものではない。このため、圧縮に伴う固体の弾性特性の変化を考慮する必要がある。

マーナハンは、体積弾性率<math>K</math>が圧力<math>P</math>の一次関数であると仮定した<ref name="Murnaghan1944"/>。

<math display="block"> K = K_0 + P\ K_0'</math>

マーナハンの状態方程式は、以下の微分方程式を積分することで得られる。

<math display="block">
P(V) = \frac{K_0}{K_0'} \left[\left(\frac{V}{V_0}\right)^{-K_0'} - 1\right]
</math>

また、次のような体積<math>V</math>と圧力<math>P</math>の関係式が得られる。

<math display="block">
V(P) = V_0 \left[1+ P \left(\frac{K'_0}{K_0}\right)\right]^{-1/K'_0}
</math>

しかし、この説明は、ポワリエによって厳密性を欠くとして批判されている<ref>Poirier (2002), p. 65.</ref>。同じ関係式は、弾性率と熱膨張係数の積が圧力に依存しないという事実から、異なる方法で導き出すことができる<ref name=Kumar1995>{{Citation |first = M. |last=Kumar |title = High pressure equation of state for solids |journal = Physica B: Condensed Matter |volume = 212 |issue=4 |pages = 391–394 |year = 1995 |doi=10.1016/0921-4526(95)00361-C|bibcode=1995PhyB..212..391K }}</ref>。この状態方程式は、古い[[ポリトロープ]]関係の一般的な場合であり<ref name="mnras">Weppner, S. P., McKelvey, J. P., Thielen, K. D. and Zielinski, A. K., "A variable polytrope index applied to planet and material models", "Monthly Notices of the Royal Astronomical Society", Vol. 452, No. 2 (Sept. 2015), pages 1375–1393, Oxford University Press also found at [https://arxiv.org/abs/1409.5525 the arXiv]</ref>、一定の指数関係を持つ。

特定の状況下、特に''[[ab initio]]''計算に関連して、エネルギーを体積の関数として表すことが好まれる場合がある<ref>Silvi (1997), p. 122.</ref>。これは、以下の関係式を基に積分することで得られる。

{{math|1=''P'' = −d''E''/d''V''}}

体積弾性率<math>K</math>の一次微分<math>K'_0</math>が1でない場合、この関係式が用いられる。

<math display="block">
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>

{| class="toccolours collapsible collapsed" width="60%" style="text-align:left"
! マーナハンの状態方程式の導出
|-
|固体には特定の平衡体積<math>V_0</math>が存在し、その体積から膨張した場合も収縮した場合もエネルギーは増加する。このような条件をみたすエネルギーと体積の最も単純な依存関係は二次関数であり、以下のようになる。

<math display="block">
E = E_0 + \frac{1}{2} K_0 \frac{(V-V_0)^2}{V_0}
</math>

さらに、次に単純なモデルでは、体積弾性率が温度によらず一定とするモデルである。

<math display="block">K_0 = - V \left( \frac{\partial P}{\partial V} \right)_T</math>

これを積分すると以下の関係式群を得る。

<math display="block">\begin{align}
P &= K_0 \ln(V_0/V) \\
V &= V_0 \exp(-P/K_0) \\
E &= E_0 + K_0 \left( V_0 - V + V \ln(V/V_0) \right)
\end{align}</math>

上式よりもパラメータの多い状態方程式の導出のため、以下の式で定義される圧力および体積膨張率との関係式を考える。

{{NumBlk||<math display="block">P = - \left( \frac{\partial E}{\partial V} \right)_S</math>|{{EquationRef|1}}}}

{{NumBlk||<math display="block">K = -V\left( \frac{\partial P}{\partial V} \right)_T</math>|{{EquationRef|2}}}}

実験的には、体積弾性率の圧力微分

{{NumBlk||<math display="block">K' =  \left( \frac{\partial K}{\partial P} \right)_T</math>|{{EquationRef|3}}}}

は圧力によってあまり変化しないことが知られている。そこで、<math>K'=K'_0</math>のような定数関数を仮定し積分すると次を得る。

{{NumBlk||<math display="block">K = K_0 + K'_0 P \qquad</math>|{{EquationRef|4}}}}

ここで<math>K_0</math>は、<math>P = 0</math>のときの<math>K</math>の値である。これを({{EquationNote|2}})に代入して[[変数分離]]すると次を得る。 

{{NumBlk||<math display="block">\frac{\mathrm dV}{V} = -\frac{\mathrm d P}{K_0 + K'_0 P}</math>|{{EquationRef|5}}}}

この両辺を積分し整理すると、次のようになる。

{{NumBlk||<math display="block">P(V) = \frac{K_0}{K'_0} \left(\left(\frac{V_0}{V}\right)^{K'_0} - 1\right)</math>|{{EquationRef|6}}}}

また、次のようにも書ける。

{{NumBlk||<math display="block">V(P) = V_0 \left(1+K'_0 \frac{P}{K_0}\right)^{-1/K'_0}</math>|{{EquationRef|7}}}}

ここで、({{EquationNote|6}}) を <math>T = 0</math>([[熱力学第三法則|したがって]]<math>S = 0</math>)のとき、<math display="inline">E = E_0 - \int_{V_0}^V P \, \mathrm dV</math>に代入すると、エネルギーの状態方程式が得られる。

{{NumBlk||<math display="block">
E(V) = E_0
 + \frac{ K_0 V }{ K_0' } \left( \frac{ (V_0/V)^{K_0'} }{ K_0' - 1 } + 1 \right)
 - \frac{ K_0 V_0 }{ K_0' - 1 }
</math>|{{EquationRef|8}}}}

多くの物質では、<math>K'_0</math>はおよそ3.5でほとんど一定している。
|}

== 利点と限界 ==
マーナハンの状態方程式は、その単純さにもかかわらず、体積弾性率<math>K_0/2</math>のおおよそ半分までの圧力範囲で実験データを再現することが可能である<ref name=Anderson1995>{{Citation |first = O.L. |last=Anderson |title = Equations of state of solids for geophysics and ceramic science, p. 179 |year = 1995 |publisher = Oxford University Press |url = https://books.google.com/books?id=qSqiz2NK7TIC&pg=PR4|isbn=9780195345278 }}</ref>。また、体積比<math>V/V_0</math>が約90%以上であれば、比較的良好な結果を提供する<ref>{{Citation | title = High-Pressure Crystallography | first = R.J.|last= Angel | chapter = Some practical aspects of studying equations of state and structural phase transitions at high pressure | pages = 21–36}}</ref>。この範囲内では、圧力の関数として体積を表す必要がある場合、マーナハンの状態方程式は他の状態方程式と比べて有利である<ref name=Hol1996>{{Citation|first = W.B.|last=Holzapfel | title = Physics of solids under strong compression | journal = Reports on Progress in Physics | volume= 59 | pages = 29–90 |year = 1996 |issue=1 | doi=10.1088/0034-4885/59/1/002|bibcode=1996RPPh...59...29H |s2cid=250909120 }}</ref>。

しかし、多くの問題において、マーナハンの状態方程式は理論的・実験的な観点から満足のいく結果を提供しないことが示されている。特に、体積比<math>V/V_0</math>が非常に小さくなる状況、すなわち[[トーマス＝フェルミ模型|トーマス・フェルミ極限]]において、体積弾性率の一次微分<math>K'_0</math>は<math>5/3</math>に収束すると予測される<ref name=Hol1996/>{{Efn|トーマス・フェルミ理論は極度に圧縮された固体を[[縮退]]した[[電子ガス]]([[フェルミ気体]])と考え、これに原子核の存在による追加の遮蔽項をとりこんだモデルである。}}。しかし、マーナハンの状態方程式では、<math>K</math><math>'</math>は一定であり、初期値に設定される。 特に、<math>K'_0 = 5/3</math>が初期値に固定しているため、この予測に一致しない。マーナハンの状態方程式を外挿して使用すると、その予測は急速に現実とかけ離れた挙動を示す<ref name=Hol1996/>。

経験的には体積弾性率<math>K</math>の一次微分<math>K</math><math>'</math>は圧力が増加すると減少する。また、弾性率の二次微分<math>K</math><math>'</math><math>'</math>は常に負の値を持つ。この観測を説明するため、同じ原則に基づく二次理論が提案されている。しかし、このアプローチにも問題がある。まず、圧力が無限に増加する極限では、この理論に基づくと体積弾性率が負になる矛盾が生じる。次に、任意の多項式展開を採用した場合でも、必ず支配的な項が無限大に発散するため、同様の矛盾が生じる<ref name=Stacey1981>{{Citation |first1 = F.D.|last1= Stacey|first2= B.J. |last2=Brennan |first3=R.D. |last3=Irvine |title = Finite strain theories and comparison with seismological data |journal = Surveys in Geophysics |volume = 4 |issue= 3|pages = 189–232 |year = 1981 |url = http://www.springerlink.com.chimie.gate.inist.fr/content/x0254w7q10k44w52/?p=8afe65621b9d4a79a3f0d9fd2466b263&pi=1 |doi=10.1007/bf01449185|bibcode= 1981GeoSu...4..189S|s2cid= 129899060}}{{dead link|date=February 2020|bot=medic}}{{cbignore|bot=medic}}</ref>。

これらの重要な限界から、マーナハンの状態方程式は事実上廃止された。W・ホルツアプフェルはこの方程式を次のように評している<ref>{{Citation | first = W.B. |last =Holzapfel |title = Equations of state for solids under strong compression | journal = Zeitschrift für Kristallographie | volume = 216 |issue =9 | pages = 473–488 |year = 2001 | doi=10.1524/zkri.216.9.473.20346|bibcode =2001ZK....216..473H |s2cid =94908666 }}</ref>。実際の圧縮データ解析には、より洗練された状態方程式が使用されている。科学界で最も広く用いられているのは、[[バーチ・マーナハンの状態方程式]]（2次または3次）であり、特に高品質なデータの解析で一般的である<ref>{{Citation |first1 = E. |last1=Boldyreva|first2= P.|last2= Dera |first3= T. Boffa|last3= Ballaran |editor = Springer |title = High-Pressure Crystallography: From Fundamental Phenomena to Technological Applications |chapter = Equations of state and their applications in geosciences |pages = 135–145}}</ref>。

この種の状態方程式は、圧力や温度による融解や固体間の相転移による密度や体積弾性率の急激な変化に対応できない、圧力に応じた複数の固体間転移を正確に反映することが困難である、といった限界がある<ref name=Stacey1981/>。

== 例 ==
マーナハンの状態方程式は、データセットに対して回帰分析を行い、弾性率<math>K_0</math>とその圧力に対する1次導関数<math>K'_0</math>の値を得ることができる。これらの係数を取得し、常圧での体積の値を知ることで、体積、密度、体積弾性率を計算することが可能である。

データセットは主に、さまざまな圧力における体積の測定値（多くの場合、X線回折によるもの）で構成されている。また、理論的なデータも利用可能であり、ab initio計算によって異なる体積値におけるエネルギーを算出し、それを回帰分析に使用することができる。こうして得られる理論的な弾性率の値は、実験的な結果と比較される。

以下の表は、マーナハンの状態方程式により得られたさまざまな材料の数値例を示している。ただし、これは単なる数値解析の例として示されたものであり、モデルの質についての評価を含むものではない。これらの結果は、マーナハンの状態方程式を使用して得られた数値解析の例として示されており、前節で述べた物理的な限界点を踏まえ、慎重に扱う必要がある。

{| class="wikitable centre"
! scope=col | 材質
! scope=col | <math>K_0</math> (GPa)
! scope=col | <math>K'_0</math>
|- 
|[[フッ化ナトリウム]]<ref name="Kumar1995" />
|46.5
|5.28
|- 
|[[塩化ナトリウム]]<ref name="Kumar1995" />
|24.0
|5.39
|- 
|[[臭化ナトリウム]]<ref name="Kumar1995" />
|19.9
|5.46
|- 
|[[ヨウ化ナトリウム]]<ref name="Kumar1995" />
|15.1
|5.59
|- 
|[[酸化マグネシウム]]<ref name="Anderson1995" />
|156
|4.7
|- 
|[[方解石]] (CaCO<sub>3</sub>)<ref>Silvi,1997. p. 123.</ref>
|75.27
|4.63
|- 
|[[菱苦土鉱]] (MgCO<sub>3</sub>)<ref>Silvi, 1997.</ref>
|124.73
|3.08
|- 
|[[炭化ケイ素]] (3C-SiC)<ref>{{Citation|first1 = K. |last1=Strössner|first2= M. |last2=Cardona |first3= W. J. |last3=Choyke |title = High pressure X-ray investigations on 3C-SiC|journal = Solid State Communications |volume = 63 |issue=2|pages = 113–114 |year = 1987 |doi = 10.1016/0038-1098(87)91176-8|bibcode=1987SSCom..63..113S}}</ref>
|248
|4.0
|}

== 拡張と一般化 ==
前節のようなモデルの改良や、批判を回避するために、マーナハンの状態方程式の一般化が提案されている。これらの一般化は通常、単純化のための仮定を取り除き、新たな調整可能なパラメータを追加する形で行われる。このアプローチにより、モデルの精度を向上させることができるが、一方で式が複雑化する可能性もある。また、追加されたパラメータの物理的な意味についての問題も提起されている。

一つの可能な戦略として、先の展開に追加項<math>P^2</math>を含めることが挙げられる<ref>{{Citation |first1 = J.R.|last1= MacDonald|first2= D.R.|last2= Powell |title= Discrimination Between Equations of State |journal = Journal of Research of the National Bureau of Standards Section A |volume = 75 |issue= 5|pages = 441 |year = 1971 |doi= 10.6028/jres.075A.035|doi-access= free }}</ref><ref>MacDonald, 1969, p. 320</ref>。この場合、次の条件を課すことになる。

<math> K = K_0 + PK_0' + P^2K_0''</math>

この微分方程式を解くことで、2次のマーナハンの状態方程式が得られる。

<math display="block">
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>

ここで、<math>\Gamma^2 = K_0'^2 - 2 K_0 K_0'' > 0</math>であり、これは<math>K_0''=0</math>を取ることで1次方程式に自然に現れる形となる。また、2次より高次の展開も原理的には可能だが<ref>{{Citation | first = Kazuhiro|last= Fuchizaki | title = Murnaghan equation of state revisited | journal = Journal of the Physical Society of Japan | volume = 75 |issue= 3 | year = 2006 | pages = 034601 |url = http://jpsj.ipap.jp/link?JPSJ/75/034601/ | doi=10.1143/jpsj.75.034601|bibcode= 2006JPSJ...75c4601F }}</ref>、その場合は各項に対して調整可能なパラメータを追加する必要がある。

他の一般化としては、以下が挙げられる。

* クマリおよびダスは、<math>K=0</math>という条件を放棄し、<math>K/K'</math>が圧力に依存しないと仮定する一般化を提案した<ref>{{Citation|first1 = M. |last1= Kumari |first2= N.|last2= Dass |title = An equation of state applied to sodium chloride and caesium chloride at high pressures and high temperatures |journal = Journal of Physics: Condensed Matter |volume = 2 |issue= 14 |pages = 3219–3229 |year = 1990 |doi=10.1088/0953-8984/2/14/006|bibcode= 1990JPCM....2.3219K |s2cid= 250827859 }}</ref>。

* クマーは、アンダーソンパラメータが体積の関数として依存することを考慮に入れた一般化を提案した。その後、この一般化された方程式は新しいものではなく、実際には[[テイト方程式]]に還元可能であることが示された<ref name="Kumar1995"/><ref>{{Citation|first1= J. |last1=Shanker|first2= B. |last2=Singh|first3=S.S. |last3=Kushwah |title = On the high-pressure equation of state for solids |journal = Physica B: Condensed Matter |volume = 229 |issue=3–4|pages = 419–420 |year = 1997 |doi=10.1016/S0921-4526(96)00528-5|bibcode=1997PhyB..229..419S}}</ref>。

== 脚注 ==
{{脚注ヘルプ}}

=== 出典 ===
{{Reflist}}注釈
{{Notelist}}

== 参考文献 ==
* {{Citation |first = J.P.|last= Poirier |title = Introduction to the physics of the Earth's interior |year = 2002 |publisher = Cambridge University Press |url = https://books.google.com/books?id=_YOCl86XPHoC&pg=PP1|isbn= 9780521663922 }}
* {{Citation |first1 = B. |last1=Silvi |first2= P. |last2=d'Arco|title = Modelling of Minerals and Silicated Materials |publisher = Kluwer Academic Publishers |year = 1997 |url = https://books.google.com/books?id=60wSAQAAIAAJ|isbn=9780792343332 }}
* {{Citation |first = J.R.|last=MacDonald |title = Review of Some Experimental and Analytical Equations of State |journal = Reviews of Modern Physics |volume = 41 |issue=2 |pages = 316–349 |year = 1969 |doi=10.1103/revmodphys.41.316|bibcode=1969RvMP...41..316M }}

== 関連項目 ==
* [[状態方程式 (熱力学)]]
* [[バーチ・マーナハンの状態方程式]]
* [[ローズ・ヴィネットの状態方程式]]
* [[ポリトロープ]]

== 外部リンク ==
* [https://web.archive.org/web/20130606143226/http://www.crystal.vt.edu/crystal/software.html EosFit]、実験データの精密化と、マーナハンの状態方程式を含むさまざまな状態方程式に対するP (V)の計算関係のためのプログラム。

{{DEFAULTSORT:まあなはんのしようたいほうていしき}}
[[Category:状態方程式 (熱力学)]]
[[Category:数学に関する記事]]
[[Category:人名を冠した数式]]
[[Category:連続体力学]]