ハイパーネット連鎖方程式

統計力学において、ハイパーネット連鎖方程式は、直接相関関数と全相関関数を関連付けるオルンシュタイン・ゼルニケ方程式を解くための閉包関係である。流体理論では、例えば動径分布関数の式を得るためによく用いられる。これは次のように与えられる。

lnyr12)lngr12)βur12)ρ[hr13)lngr13)βur13)]hr23)dr3{\displaystyle \ln y(r_{12})=\ln g(r_{12})+\beta u(r_{12})=\rho \int \left[h(r_{13})-\ln g(r_{13})-\beta u(r_{13})\right]h(r_{23})\,d\mathbf {r_{3}} ,\,}

ここで、 は分子の数密度、 は動径分布関数、はペア間の直接相互作用です。は 熱力学温度、 はボルツマン定数ですρNV{\displaystyle \rho ={\frac {N}{V}}}hr)gr)1{\displaystyle h(r)=g(r)-1}gr){\displaystyle g(r)}ur){\displaystyle u(r)}β1kBT{\displaystyle \beta ={\frac {1}{k_{\rm {B}}T}}}T{\displaystyle T}kB{\displaystyle k_{\rm {B}}}

導出

直接相関関数は、N  − 2個の他の粒子を含む系における2つの粒子間の直接的な相関を表します。これは次のように表すことができます

cr)gtotalr)gdrctr){\displaystyle c(r)=g_{\rm {total}}(r)-g_{\rm {indirect}}(r)\,}

ここで(平均力のポテンシャル)と は、ペア間の直接的な相互作用を含まない動径分布関数である。つまり と書く。したがって、次のように 近似できる。gtotalr)gr)経験[βr)]{\displaystyle g_{\rm {total}}(r)=g(r)=\exp[-\beta w(r)]}r){\displaystyle w(r)}gdrctr){\displaystyle g_{\rm {indirect}}(r)}ur){\displaystyle u(r)}gdrctr)経験{β[r)ur)]}{\displaystyle g_{\rm {indirect}}(r)=\exp\{-\beta [w(r)-u(r)]\}}cr){\displaystyle c(r)}

cr)βr)β[r)ur)]{\displaystyle c(r)=e^{-\beta w(r)}-e^{-\beta [w(r)-u(r)]}.\,}

上記の式の間接部分を展開し、関数を導入すると、次のように近似できます。 gr){\displaystyle g(r)}yr)βur)gr)gdrctr)){\displaystyle y(r)=e^{\beta u(r)}g(r)(=g_{\rm {indirect}}(r))}cr){\displaystyle c(r)}

cr)βr)1β[r)ur)]gr)1lnyr)fr)yr)[yr)1lnyr)]HNC){\displaystyle c(r)=e^{-\beta w(r)}-1+\beta [w(r)-u(r)]\,=g(r)-1-\ln y(r)\,=f(r)y(r)+[y(r)-1-\ln y(r)]\,\,({\text{HNC}}),}

と。 fr)βur)1{\displaystyle f(r)=e^{-\beta u(r)}-1}

この方程式はハイパーネット連鎖方程式の本質である。

hr)cr)gr)1cr)lnyr){\displaystyle h(r)-c(r)=g(r)-1-c(r)=\ln y(r).}

この結果をオルンシュタイン・ゼルニケ方程式に代入すると

hr12)cr12)ρcr13)hr23)dr3{\displaystyle h(r_{12})-c(r_{12})=\rho \int c(r_{13})h(r_{23})d\mathbf {r} _{3},}

ハイパーネット連鎖方程式が得られる:

lnyr12)lngr12)βur12)ρ[hr13)lngr13)βur13)]hr23)dr3{\displaystyle \ln y(r_{12})=\ln g(r_{12})+\beta u(r_{12})=\rho \int \left[h(r_{13})-\ln g(r_{13})-\beta u(r_{13})\right]h(r_{23})\,d\mathbf {r_{3}} .\,}

参照