tecktom.com

量子の粒子は古典の輪:統計の対応

2026年10月17日6分で読めますAI生成教科書『虎の巻』第3部 第6章 虚時間と統計力学
目次
  1. 同じ式の読み替え
  2. 木下の層

AI 生成・木下が未勉強

量子力学の様々な構成方法 第3部 第6章 §3.6.5 ★★ 時間分割した ZPZ_P は、ばねでつないだ PP 個のビーズの輪の古典の配置積分と同じ式で、量子の統計力学を古典の統計力学の計算法で扱える

同じ式の読み替え

式 (3.6.5) の指数を β\beta でくくり直すと

ZP=(mP2πβℏ2)P/2∫∏k=0P−1dxk  exp⁡{−β∑k=0P−1[kP2(xk+1−xk)2+V(xk)P]},kP=mPβ2ℏ2(3.6.22)Z_P=\Bigl(\frac{mP}{2\pi\beta\hbar^2}\Bigr)^{P/2}\int\prod_{k=0}^{P-1}\mathrm{d} x_k\;\exp\Bigl\{-\beta\sum_{k=0}^{P-1}\Bigl[\frac{k_P}2(x_{k+1}-x_k)^2+\frac{V(x_k)}P\Bigr]\Bigr\},\qquad k_P=\frac{mP}{\beta^2\hbar^2}\tag{3.6.22}

となる(1ℏ⋅m2ϵ=mP2βℏ2=βkP2\frac1\hbar\cdot\frac{m}{2\epsilon}=\frac{mP}{2\beta\hbar^2}=\beta\frac{k_P}2、1ℏϵV=βVP\frac1\hbar\epsilon V=\frac{\beta V}P)。右辺は、温度 TT の古典の系で、PP 個の粒子(ビーズ)がばね定数 kPk_P のばねで輪につながり、各ビーズにポテンシャル V/PV/P がかかっているときの配置積分である(図 3.6.4)。量子の1粒子の熱平衡が、古典の輪のポリマーの熱平衡と同じ数になる。分子の系でこの読み替えを系統的に使ったのは Chandler–Wolynes (1981), J. Chem. Phys. 74, 4078 で、この輪の配置を Monte Carlo 法で標本化して量子の熱平衡を計算する方法(経路積分 Monte Carlo)は液体ヘリウムの計算で発展した(総説 Ceperley (1995), Rev. Mod. Phys. 67, 279)。

図

図 3.6.4 時間分割した虚時間の経路(P=12P=12)を古典の輪のポリマーとして描いたもの。ビーズ kk は虚時間 kβℏ/Pk\beta\hbar/P での位置 xkx_k を表し、隣の時刻どうしがばねでつながる。トレースの周期境界条件 xP=x0x_P=x_0 が輪を閉じる。輪の広がりは熱的波長の程度で、高温ほど縮んで1点(古典の粒子)に近づく。

この読み替えでは、量子の性質が輪の形に移る。高温(β\beta が小さい)ではばね kP∝1/β2k_P\propto1/\beta^2 が強く、輪は1点に縮んで古典の粒子になる。低温ではばねが弱くなり、輪は熱的波長の程度まで広がり、ビーズがポテンシャルの起伏の広い範囲にまたがる。前節の Feynman のぼかしは、この輪の広がりを Gauss 分布で置き換えたものである。空間 dd 次元の量子の粒子は、虚時間の方向が1つ増えた d+1d+1 次元の古典の物体(d=1d=1 なら輪=閉じた紐)に写る。

多粒子の Bose 粒子では、トレースを対称化した状態でとるので、終わりの位置が始めの位置の入れ替えであってもよい。輪の言葉では、ll 個の粒子の輪がつながって、周期 lβℏl\beta\hbar の長い1本の輪になる配置も数えることになる。入れ替えの輪の和から Bose 分布と Fermi 分布が出る計算は、解析力学の巻(解析力学から様々な理論分野を組み立てる 第5部第1章)と第6部第5章(第6部第5章)にある。

例題 3.6.5 自由な輪の大きさ ★★

オリジナルV=0V=0 として、式 (3.6.22) の輪を考える。重心を固定し、揺らぎだけを見る。

  1. (1) 隣のビーズの距離の分散が ⟨(xk+1−xk)2⟩=βℏ2mP(1−1P)\langle (x_{k+1}-x_k)^2 \rangle =\frac{\beta\hbar^2}{mP}\bigl(1-\frac1P\bigr) であることを示せ。
  2. (2) 重心からの距離の2乗平均(回転半径の2乗)が 1P∑k⟨(xk−xˉ)2⟩=βℏ212m(1−1P2)\frac1P\sum_k\langle (x_k-\bar x)^2 \rangle =\frac{\beta\hbar^2}{12m}\bigl(1-\frac1{P^2}\bigr) であることを示し、P→∞P\to\infty で式 (3.6.17) に戻ることを確かめよ。
  3. (3) ⁴He 原子の 2.17 K2.17\,\mathrm K と電子の 300 K300\,\mathrm K での βℏ2/(12m)\sqrt{\beta\hbar^2/(12m)} を求めよ。
解答

(1) 結合のベクトル Δk=xk+1−xk\Delta_k=x_{k+1}-x_k を独立変数にとると、重みは ∏ke−βkPΔk2/2\prod_k\mathrm{e}^{-\beta k_P\Delta_k^2/2} で、輪が閉じる条件 ∑kΔk=0\sum_k\Delta_k=0 がつく。条件なしなら各 Δk\Delta_k は分散 a=1/(βkP)=βℏ2/(mP)a=1/(\beta k_P)=\beta\hbar^2/(mP) の独立な Gauss 変数である。独立な Gauss 変数を和が0の超平面に制限すると、共分散は a(δjk−1/P)a(\delta_{jk}-1/P) になる(和の方向の成分だけを取り除く射影)。したがって ⟨Δk2⟩=a(1−1/P)=βℏ2mP(1−1P)\langle \Delta_k^2 \rangle =a(1-1/P)=\frac{\beta\hbar^2}{mP}(1-\frac1P) である。

(2) Fourier 成分で書くと、ばねの項は βkP∑k∣Δ∣2\beta k_P\sum_k\lvert \Delta\rvert ^2 を対角化して、波数 q=2πj/Pq=2\pi j/P(j=1,…,P−1j=1,\dots,P-1)の成分の分散が 1βkP⋅4sin⁡2(πj/P)\frac{1}{\beta k_P\cdot4\sin^2(\pi j/P)} になる。回転半径の2乗はこの分散の和を PP で割ったもので、1P∑j=1P−1βℏ2mP14sin⁡2(πj/P)\frac1P\sum_{j=1}^{P-1}\frac{\beta\hbar^2}{mP}\frac1{4\sin^2(\pi j/P)} である。和の公式 ∑j=1P−11/sin⁡2(πj/P)=(P2−1)/3\sum_{j=1}^{P-1}1/\sin^2(\pi j/P)=(P^2-1)/3 を使うと βℏ2mP2⋅P2−112=βℏ212m(1−1P2)\frac{\beta\hbar^2}{mP^2}\cdot\frac{P^2-1}{12}=\frac{\beta\hbar^2}{12m}(1-\frac1{P^2}) となる。P→∞P\to\infty で βℏ2/(12m)\beta\hbar^2/(12m)、すなわち連続な経路の揺らぎの分散 (3.6.17) に戻る。PP が有限の輪は、連続な経路の揺らぎのうち高い Matsubara 成分を落としたものである。

(3) ⁴He の 2.17 K2.17\,\mathrm K では ℏ2/(12mkBT)=0.068 nm\sqrt{\hbar^2/(12m k_{\mathrm B} T)}=0.068\,\mathrm{nm}(=λT/24π=\lambda_T/\sqrt{24\pi}、λT=0.592 nm\lambda_T=0.592\,\mathrm{nm})、電子の 300 K300\,\mathrm K では 0.496 nm0.496\,\mathrm{nm} である。⁴He の輪そのものは原子の間隔 0.36 nm0.36\,\mathrm{nm} より小さいが、隣の原子と入れ替わる配置の重みは e−πd2/λT2\mathrm{e}^{-\pi d^2/\lambda_T^2} 程度の因子で決まり、λT\lambda_T が間隔を上回るとこの入れ替えが無視できなくなる。

木下の勉強ステータス

木下はまだこの記事を勉強していません。

AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com

← AI生成教科書『虎の巻』の記事一覧へ