めもめも

このブログに記載の内容は個人の見解であり、必ずしも所属組織の立場、戦略、意見を代表するものではありません。

多体系の密度行列の非対角成分が0になる理由を系に対するノイズ（ランダムな摂動）の観点で説明する

何の話かというと

Statistical Physics: Volume 5

Statistical Physics: Volume 5

Landau の Statistical Physics の §6 では、密度行列の非対角成分が平衡状態では無視できる理由として、「時間依存性がないことが平衡状態の定義であり、したがって、時間依存性を持つ非対角成分が0の場合を考える」となかなかトートロジカルな説明がなされています。一応、脚注において、「系に含まれる粒子数が増加するなどして、相互作用の重要性が低下することにより非対角成分が0になる」と説明されていますが、この部分がイマイチ直感的に理解できません。というわけで、この点について、もう少し定量的な説明を試みたのが以下の内容です。

まず、他の粒子からの相互作用を1粒子のハミルトニアンに対する摂動パラメーター q と考えて、1粒子のエネルギー固有値が q に依存した関数になると考えます。

　\hat H=\hat H_0 + \hat V(q)\ \Rightarrow\ \omega_i = \omega_i(q)（\omega_i=\frac{E_i}{\hbar}）

したがって、1粒子のエネルギー固有状態は、時間依存を含めて、

　\mid \phi_i\rangle e^{-i\omega_i(q)t}

となります。（厳密には、\mid \phi_i\rangle も q に依存していますが、この点については後述します。）

この時、一般の重ね合わせ状態

　\mid\psi\rangle=\sum_i c_i{\mid\phi_i\rangle} e^{-i\omega_i(q)t}

について、Pure state としての密度行列を計算すると次が得られます。

　\rho = {\mid\psi\rangle}{\langle\psi\mid} = \sum_i |c_i|^2{\mid\phi_i\rangle}{\langle\phi_i\mid}+\sum_{i\ne j}c_ic_j^*{\mid\phi_i\rangle}{\langle\phi_j\mid} e^{-i\left\{\omega_i(q)-\omega_j(q)\right\}t} ---- (1)

一般に、\omega が大きい古典極限では非対角成分は時間変動が大きく、マクロな時間レベルで平均して0になると考えられます。ただし、これは、統計力学とは無関係の古典極限の話です。ここでは、非古典的な量子状態に対する統計的性質を考える必要があります。

今、（粒子が多数存在するなどの理由で）粒子間の相互作用が乱雑で、摂動パラメーター自体が時間に対して激しく変動すると仮定すると、\omega_i(q(t)) という時間依存により、\omega がそれほど大きくない領域においても、位相 -i\omega_i(q(t)) は激しく変動するため、古典極限と同様の近似として、非対角成分を0とみなすことができます。一般に、マクロな時間内に \Delta \omega(q) \sim 2\pi、すなわち、\Delta E(q) \sim h 程度の変動があれば十分です。

この話を一般化すると、系に対するランダムな摂動（ノイズ）がある場合は、非古典領域においても密度行列の非対角成分を0とする、（平衡系の）統計力学的な取り扱いが正当化されることになります。

なお、(1) でマクロな時間平均を取ると、上述の理由で非対角成分は0になりますが、同時に、\mid \phi_i\rangle の q 依存性についても平均が行われます。これは、他の粒子からの影響を平均場で近似して、平均場のもとに得られたエネルギー固有状態 \mid\phi_i\rangle に置き換える操作と理解できます。平均場が0になる場合（つまり、摂動項の平均が0になる場合）は、近似として、もとのハミルトニアンの固有状態をそのまま利用することができます。