等エンタルピー定圧集団

等エンタルピー定圧集団（とうエンタルピーていあつしゅうだん、テンプレート:Lang-en-short）は統計力学的集団の1 つ。エネルギー、体積、テンプレート:仮リンクを熱力学的独立変数（多くの場合は定数）であるミクロカノニカル集団に対して、等エンタルピー定圧集団では、エンタルピー(テンプレート:Mvar)、圧力(テンプレート:Mvar)、粒子数(テンプレート:Mvar)を独立変数とし、そこからNPHアンサンブルとも呼ばれる。定圧系のモデル化としては最も基礎的なアンサンブルだが、専門文献ではほとんど言及されない。定圧下の熱力学系のシミュレーションに主に用いられる。

1980年に物理学者のH. C. Andersenによって発表された。

モデルの性質と特徴

等エンタルピー定圧系の図式

等エンタルピー定圧集団により説明される単純な物理系としては、熱を通さない断熱ピストンと一定圧力に維持された容積浴を結合させた系があげられる。

この系のエンタルピーの全微分は以下のように書ける。

$\mathrm{d}H=\mathrm{d}U+p\mathrm{d}V+V\mathrm{d}p=(\mathrm{d}q+\mathrm{d}w)+p\mathrm{d}V+V\mathrm{d}p=\mathrm{d}q+V\mathrm{d}p$

ピストンは断熱されており、また圧力は一定なので $\mathrm{d}H=\mathrm{d}q=\mathrm{d}p=0$ が成り立ち、したがって等エンタルピーが成り立つ。

位相空間上の体積

この系のハミルトニアンは以下のように与えられる。

$\mathscr{H}=\mathscr{H}(q_i,p_i,x_w,p_w,V,N)=H+PV$

ここで、テンプレート:Mvar, テンプレート:Mvarは各粒子の位置座標および運動量、テンプレート:Mvar, テンプレート:Mvarはピストンの位置座標および運動量、テンプレート:Mvarはシステムの体積、テンプレート:Mvarは粒子数、テンプレート:Mvarは系のエンタルピー、テンプレート:Mvarは浴の圧力である。系の位相空間全体を積分し、それを順列数でわることで位相空間上の体積が得られる。これによりエントロピーが導出できることが保証される。

$\mathrm{\Phi }(H,P,N)=\frac{1}{N!}\underset{\mathscr{H}\le H-PV}{\int }\mathrm{d}q\mathrm{d}p\mathrm{d}{x}_w\mathrm{d}{p}_w$

ここで定数倍は物理量に影響を与えないため省略した。巨視的な系ではテンプレート:Mathの影響は無視でき、テンプレート:Mathはテンプレート:Mathで置き換えることができる。

$\mathrm{\Phi }(H,P,N)=\frac{1}{N!}\underset{\mathscr{H}(q,p,V,N)\le H-PV}{\int }\mathrm{d}q\mathrm{d}p\mathrm{d}V$

簡単のために体積要素を$\mathrm{d}\tau =\mathrm{d}q\mathrm{d}p\mathrm{d}V$と定義し、ヘヴィサイドの階段関数を用いると上式は次のように書ける。

$\mathrm{\Phi }(H,P,N)=\frac{1}{N!}\int \mathrm{\Theta }[\mathscr{H}-(H-PV)]\mathrm{d}\tau$

状態密度と確率密度分布

状態密度は位相空間上の体積を微分することにより簡単に得ることができる。

$\omega (H,P,N)=(\frac{\partial \mathrm{\Phi }}{\partial H}{)}_{P,N}=\frac{1}{N!}\int \delta [\mathscr{H}-(H-PV)]\mathrm{d}\tau$

ここでテンプレート:Mathはディラックのデルタ関数である。

エネルギーが範囲$H-PV\le \mathscr{H}<H-PV+\delta H$にある確率は以下のように計算できる。

$W(q,p,V)=\frac{1}{\omega \delta H}=\frac{\delta [\mathscr{H}-(H-PV)]}{\omega }$

ここでテンプレート:Mvarはディラックのデルタ関数ではなくエンタルピーの微小変化を意味する。

確率密度分布が得られたので、これを用いて状態量を計算することができる。一般に関数テンプレート:Mvarの期待値は以下のように計算できる。

$\overline{f}=\int Wf\mathrm{d}\tau =\frac{1}{\omega }\int f\delta [\mathscr{H}-(H-PV)]\mathrm{d}\tau =\frac{1}{\omega }(\frac{\partial }{\partial H}\int f\mathrm{\Theta }[\mathscr{H}-(H-PV)]\mathrm{d}\tau {)}_{P,N}$

簡単な例としては、エネルギー(ハミルトニアン)の期待値を計算することができ、$\overline{E}=H-P\overline{V}$のように定義式と一致する。ただし、エネルギーおよび体積はアンサンブル平均であり、系の独立変数ではない。

等配分の法則

ミクロカノニカルアンサンブルと同様に以下が成り立つ。

$⟨x_j\frac{\partial \mathscr{H}}{\partial x_k}⟩=\frac{\mathrm{\Phi }}{\omega }{\delta }_{jk}$

そして理想気体のハミルトニアンを代入すると、NPHアンサンンブルにおける温度が以下のとおり得られる。

$k_{B}T=\frac{\mathrm{\Phi }}{\omega }$

ここでテンプレート:Mvarは温度、テンプレート:Mathはボルツマン定数である。

断熱定理

ハミルトニアンが追加のパラメータテンプレート:Mvarに依存すると仮定する。このパラメータに共役な一般化力は$\frac{\partial \mathscr{H}}{\partial a}$であたえられる。一般化力のアンサンブル平均を計算すると、以下を得る。

$⟨\frac{\partial \mathscr{H}}{\partial a}⟩=-\frac{1}{\omega }(\frac{\partial \mathrm{\Phi }}{\partial a}{)}_{H,P,N}=(\frac{\partial H}{\partial a}{)}_{\mathrm{\Phi },P,N}$

統計力学と熱力学の関係

位相空間上の体積(実際には分配関数)を用いて、系のエントロピーを得ることができる。

$S(H,P,N)=k_B\ln \mathrm{\Phi }(H,P,N)$

ここから熱力学恒等式をもちいれば種々の状態量を導出できる。

熱力学量と統計的ゆらぎ

独立変数(粒子数、圧力、エンタルピー)以外の熱力学的状態量は統計的に定義され、その値は一定ではない。そのゆらぎの度合いは分散により定義される。

$(\mathrm{\Delta }f{)}^2=⟨(f-\overline{f}{)}^2⟩=⟨f^2⟩-⟨f{⟩}^2$

種々の平均値は前述のとおり計算できる。

重要な熱力学量は以下のように算出できる。

• 体積

$⟨V⟩=\frac{1}{\omega }(\frac{\partial \chi }{\partial H}{)}_P;⟨V^2⟩=-\frac{1}{\omega }(\frac{\partial \chi }{\partial P}{)}_H;(\mathrm{\Delta }V{)}^2=⟨V^2⟩-⟨V{⟩}^2$

ここで $\chi (H,P)=\frac{1}{N!}\int V\cdot \mathrm{\Theta }[\mathscr{H}-(H-PV)]\mathrm{d}\tau$である。マクスウェルの関係式と偏微分の性質を利用すると、体積の分散の別の表式が得られる。

$(\mathrm{\Delta }V{)}^2=-\frac{1}{\omega }(\frac{\partial \chi }{\partial P}{)}_S$

• エネルギー

$(\mathrm{\Delta }E{)}^2=P^2\cdot (\mathrm{\Delta }V{)}^2$

• 運動エネルギー（ただし、位置エネルギーが粒子の質量に依存しない場合にのみ正しい）

$(\mathrm{\Delta }K{)}^2=\frac{3}{2}N(k_{B}T{)}^2[1-\frac{3N{k}_B}{2C_P}]$

理想気体

自由古典的理想気体

ハミルトニアン

$\mathscr{H}=\sum _i\frac{p_i^2}{2m}$

に対して位相空間体積は次のように得られる。

$\mathrm{\Phi }(N,P,H)=\frac{(2\pi m{)}^{\frac{3N}{2}}}{[N!(\frac{3N}{2})!]}\underset{0}{\overset{\frac{H}{P}}{\int }}{V}^N(H-PV{)}^{\frac{3N}{2}}dV=\frac{(2\pi m{)}^{\frac{3N}{2}}}{(\frac{5N}{2}+1)!}\frac{H^{\frac{5N}{2}+1}}{P^{N+1}}$

したがって等エンタルピー・等圧集団における理想気体の状態方程式は以下のように得られる。

$PV=(N+1)k_{B}T$

熱力学的極限ではテンプレート:Mathであるから、この式は本質的には古典理想気体の標準的な状態方程式と同等である。