粗視化

粗視化そしか: coarse graining)とは、ある変数空間で定義された連続的な物理量を、その変数を任意の単位スケールで離散化し、単位スケール内で物理量の平均を取ることで、その物理量そのものも離散化し情報量を減らす手法。単純に言うと、モザイクを掛ける事、解像度を下げる事がこれに相当する。解像度の高い空間から低い空間への射影操作であるとも言える。可逆な基礎方程式(ニュートン方程式など)から不可逆な現象(エントロピー増大則など)を説明するためのトリックの一つだと考えられている。

分布関数の粗視化

物理学ではよく粗視化された分布関数が用いられる。厳密な(粗視化されていない)分布関数 f {\displaystyle \left.f\right.} に対し、粗視化された分布関数 F {\displaystyle \left.F\right.} は次で定義される。

F = 1 A δ Γ δ Γ f d Γ {\displaystyle F={\frac {1}{A_{\delta \Gamma }}}\int _{\delta \Gamma }fd\Gamma }

ここで δ Γ {\displaystyle \left.\delta \Gamma \right.} は粗視化の領域であり、2次元ならば A δ Γ {\displaystyle \left.A_{\delta \Gamma }\right.} はその面積、 d Γ {\displaystyle \left.d\Gamma \right.} は面積素である。

例えば、統計力学エントロピーの計算をする際に位相空間を細胞に分割する操作は

F ( p , x ) = 1 Δ p Δ x p j p j + Δ p x i x i + Δ x f ( p ¯ , x ¯ ) d p ¯ d x ¯ {\displaystyle F(p,x)={\frac {1}{\Delta p\Delta x}}\int _{p_{j}}^{p_{j}+\Delta p}\int _{x_{i}}^{x_{i}+\Delta x}f({\bar {p}},{\bar {x}})d{\bar {p}}d{\bar {x}}}  但し  p j p < p j + 1 ,   x i x < x i + 1 {\displaystyle p_{j}\leq p<p_{j+1},\ x_{i}\leq x<x_{i+1}}

等の様に書けるだろう。この粗視化された分布関数に対してエントロピー S [ F ] {\displaystyle \left.S[F]\right.} が定義される。

S [ F ] = k B F ln F   d p d x {\displaystyle S[F]=-k_{B}\int F\ln F\ dpdx}

ここで k B {\displaystyle \left.k_{B}\right.} ボルツマン定数である。

粗視化とエントロピー増加則

厳密な分布関数 f {\displaystyle \left.f\right.} に対して同様な「エントロピー」を定義する。

S [ f ] = k B f ln f   d p d x {\displaystyle S[f]=-k_{B}\int f\ln f\ dpdx}

ここで k B {\displaystyle \left.k_{B}\right.} ボルツマン定数 ln {\displaystyle \left.\ln \right.} 自然対数。この時間変化は次の様になる。

d S [ f ] d t = k B ( 1 + ln f ) d f d t   d p d x {\displaystyle {\frac {dS[f]}{dt}}=-k_{B}\int \left(1+\ln f\right){\frac {df}{dt}}\ dpdx}

ところが、リウヴィル方程式より、

d f d t = f t + f p p ˙ + f x x ˙ = 0 {\displaystyle {\frac {df}{dt}}={\frac {\partial f}{\partial t}}+{\frac {\partial f}{\partial p}}{\dot {p}}+{\frac {\partial f}{\partial x}}{\dot {x}}=0}

であるから、

d S [ f ] d t = 0 {\displaystyle {\frac {dS[f]}{dt}}=0}

つまり、厳密な分布関数に対するエントロピーは時間依存性を持たない事がわかる。一般的によく言われる「エントロピーの増加法則」は、何らかの情報に対して粗視化された分布関数に対して初めて成り立つ。

簡単な例で考えてみよう。

「配置数」w 配置数W
1     2 1
2     2 2
3     2 2
4     2 2
5     2 2
6     2 1

4つの箱イロハニと2つのボール●がある。2つのボールは区別可能であり、1つの箱にボールは一つずつしか入れられない。そしてボールは4つの箱のいずれかに等確率で入っているとする(そのメカニズムは問わない)。2つのボールの違いには興味が無いので、その配置パターンは全部で6種類考えられる(上図)。

それぞれのパターンに対し、「配置数」wを計算する。

w = N ! i n i ! {\displaystyle w={\frac {N!}{\prod _{i}n_{i}!}}}

ここでNはボールの総数2であり、 n i {\displaystyle n_{i}} (i=イロハニ)はそれぞれの箱に入っているボールの数である。ところがこの場合「配置数」wは全て2となることがわかる。例えばこの系の初期状態がパターン1であったとした場合、それから時間発展して他のパターンへと移り変わったとしても「配置数」wは変化せず、それに対応した「エントロピー」sも初期状態から変化しない事がわかる。

s = k B ln w {\displaystyle \left.s=k_{B}\ln w\right.}

k B {\displaystyle \left.k_{B}\right.} ボルツマン定数 ln {\displaystyle \left.\ln \right.} 自然対数

ここで箱を2種類に分ける。イロを白グループ、ハニを桃グループとする。そして箱単位ではなくこれらのグループ単位で配置数を考え、それをWとする。

W = N ! c g n c g ! {\displaystyle W={\frac {N!}{\prod _{cg}n_{cg}!}}}

ここでNはボールの総数2であり、 n c g {\displaystyle n_{cg}} (cg=or)はそれぞれのグループに入っているボールの数である。この箱のグループ分け操作は、箱の違いに関する詳細な情報を破棄するという粗視化操作に相当する。するとこの場合、パターンごとに配置数Wに変化が現れる事がわかる(上図)。例えばこの系の初期状態がパターン1であったとした場合、それから時間発展して他のパターンへと移り変わったとすると、多くの場合(パターン6以外)において配置数Wは増加する(1→2)。それに対応したエントロピーSも初期状態から増加する。

S = k B ln W {\displaystyle \left.S=k_{B}\ln W\right.}

k B {\displaystyle \left.k_{B}\right.} ボルツマン定数 ln {\displaystyle \left.\ln \right.} 自然対数

このモデルでの「1つの箱にボールは一つずつしか入れられない。」という仮定は、例えば多粒子系を考えた時、複数の粒子の座標が厳密に重なり合う確率はほとんど無い、という事を表している。

参考文献

  • G.M.ザスラフスキー 「カオス-古典および量子力学系-」 現代工学社 ISBN 4-87472-140-0