tecktom.com

この章のまとめ

2026年10月17日9分で読めますAI生成教科書『虎の巻』第3部 第6章 虚時間と統計力学

AI 生成・木下が未勉強

量子力学の様々な構成方法 第3部 第6章 §3.6.7 ★★★ 虚時間の経路積分は実時間の経路積分の解析接続で、厳密な測度を持ち、量子の熱平衡を古典の統計力学の形に書く

この章は、実時間の伝播関数を複素時間の下半平面へ延ばすことから出発した。ハミルトニアンが下に有界なら e−iH^z/ℏ\mathrm{e}^{-\mathrm{i}\hat Hz/\hbar} は下半平面で解析的で、虚軸上の値が e−βH^\mathrm{e}^{-\beta\hat H} になる(命題 3.6.1)。第3部第3章の Mehler 核への t=−iβℏt=-\mathrm{i}\beta\hbar の代入はこの命題で正当化され、平方根の枝は t−i0t-\mathrm{i}0 の約束で決まる。Trotter の積公式を e−βH^\mathrm{e}^{-\beta\hat H} に当てると、分配関数は Euclid 作用の重み e−SE/ℏ\mathrm{e}^{-S_{\mathrm E}/\hbar} をつけた周期 βℏ\beta\hbar の閉じた経路の和になる。重みが正なので、この和は Wiener 測度の上の期待値(Feynman–Kac の公式)として厳密に定義でき、実時間の経路積分にない数学的な足場を持つ。

閉じた経路を Matsubara 成分に分けると、調和振動子の ZZ は古典の 1/(βℏω)1/(\beta\hbar\omega) と無限積の積として出て、量子の補正は全部その無限積に含まれる。同じ Gauss 積分から虚時間の相関関数 G(τ)G(\tau) が出て、その減衰率が励起のエネルギーを与え、τ↦βℏ−τ\tau\mapsto\beta\hbar-\tau の対称性が KMS 条件に対応する。経路を重心と揺らぎに分けると、揺らぎの分散は βℏ2/(12m)\beta\hbar^2/(12m) で、そのぼかしで作った古典の ZFZ_{\mathrm F} が厳密な ZZ の下限になり、ℏ2\hbar^2 の次数まで正しい(第5部第1章の VeffV_{\mathrm{eff}})。高温の極限では古典の ∫dx dp2πℏe−βH\int\frac{\mathrm{d} x\,\mathrm{d} p}{2\pi\hbar}\mathrm{e}^{-\beta H} に戻る。時間分割した ZPZ_P は古典の輪のポリマーの配置積分と同じ式で、ビーズ1個の輪が古典の粒子、ビーズを増やした輪が量子の粒子である。橋の節では、熱平衡の核の対称性・合成則・初期条件から下に有界なハミルトニアンが一意に取り戻せることを示し、逆向きの同値を完成させた。

総合問題 非調和振動子の熱平衡を4つの方法で求める ★★★

オリジナルH^=p^22m+12mω2x^2+λx^4\hat H=\frac{\hat p^2}{2m}+\frac12m\omega^2\hat x^2+\lambda\hat x^4 の自由エネルギー FF を、λ\lambda の1次まで求める。λ=0\lambda=0 の熱平衡についての平均を ⟨⋅⟩0\langle \cdot \rangle _0、s=βℏωs=\beta\hbar\omega、nˉ=1/(es−1)\bar n=1/(\mathrm{e}^s-1) とする。1次の摂動論では F=F0+λ⟨x4⟩0+O(λ2)F=F_0+\lambda\langle x^4 \rangle _0+\mathcal{O}(\lambda^2) である(Z=Z0⟨e−λℏ∫x4dτ⟩0Z=Z_0\langle \mathrm{e}^{-\frac\lambda\hbar\int x^4\mathrm{d}\tau} \rangle _0 を1次まで展開し、τ\tau によらない平均を使う)。

  1. (1) Matsubara 成分を使って ⟨x2⟩0=G(0)\langle x^2 \rangle _0=G(0) を求め、G(0)=ℏ2mω(2nˉ+1)G(0)=\frac\hbar{2m\omega}(2\bar n+1) と書けることを示せ。(使う道具:3.6.3節、第4部第4章の Bose 分布)
  2. (2) 経路積分の重みが Gauss 型であることを使い、⟨x(τ)4⟩0=3G(0)2\langle x(\tau)^4 \rangle _0=3G(0)^2 を示せ。(使う道具:3.6.2・3.6.3節)
  3. (3) 演算子の側で ⟨n∣x^4∣n⟩=3(ℏ2mω)2(2n2+2n+1)\langle n | \hat x^4 | n \rangle =3\bigl(\frac\hbar{2m\omega}\bigr)^2(2n^2+2n+1) を示し、熱平衡の重み (1−e−s)e−sn(1-\mathrm{e}^{-s})\mathrm{e}^{-sn} で平均すると (2) と一致することを確かめよ。(使う道具:第2部第4章の a^,a^†\hat a,\hat a^\dagger、第4部第4章の幾何分布)
  4. (4) Feynman のぼかしで同じ量を近似する。x4x^4 のぼかし (x4)σ(x^4)_\sigma を求め、重心の分布の分散 c=1/(βmω2)c=1/(\beta m\omega^2) を使って λ\lambda の1次の補正が 3λ(c+σ2)23\lambda(c+\sigma^2)^2 になることを示せ。ℏ=m=ω=1\hbar=m=\omega=1、s=1s=1 で厳密な値と比べ、ずれの向きを説明せよ。(使う道具:3.6.4節、第5部第1章)
  5. (5) β→∞\beta\to\infty で FF の1次の補正が基底エネルギーのずれ ΔE0=3λ(ℏ/(2mω))2\Delta E_0=3\lambda(\hbar/(2m\omega))^2 になることを示せ。(4) の近似はこの極限で何を与えるか。(使う道具:3.6.6節、第3部第3章)
解答

(1) 式 (3.6.15) で τ=0\tau=0 とおくと G(0)=ℏ2mωcoth⁡s2G(0)=\frac\hbar{2m\omega}\coth\frac s2 である。coth⁡s2=es+1es−1=1+2es−1=2nˉ+1\coth\frac s2=\frac{\mathrm{e}^s+1}{\mathrm{e}^s-1}=1+\frac2{\mathrm{e}^s-1}=2\bar n+1 なので G(0)=ℏ2mω(2nˉ+1)G(0)=\frac\hbar{2m\omega}(2\bar n+1) となる。零点の揺らぎ ℏ2mω\frac\hbar{2m\omega} に、Bose 分布の占有数に比例する熱の揺らぎが足された形である。

(2) λ=0\lambda=0 の重み e−SE/ℏ\mathrm{e}^{-S_{\mathrm E}/\hbar} は経路について2次なので、各時刻の x(τ)x(\tau) は平均0の Gauss 変数(成分 xnx_n の1次結合)で、分散は G(0)G(0) である(τ\tau によらない)。平均0・分散 vv の Gauss 変数では ⟨x4⟩=3v2\langle x^4 \rangle =3v^2 なので ⟨x(τ)4⟩0=3G(0)2\langle x(\tau)^4 \rangle _0=3G(0)^2 を得る。したがって F=F0+3λG(0)2+O(λ2)F=F_0+3\lambda G(0)^2+\mathcal{O}(\lambda^2) である。

(3) x^=ℏ/(2mω)(a^+a^†)\hat x=\sqrt{\hbar/(2m\omega)}(\hat a+\hat a^\dagger) で、(a^+a^†)4(\hat a+\hat a^\dagger)^4 の対角要素は、生成と消滅を2回ずつ含む6通りの並びの和である。∣n⟩| n \rangle に当てると、a^a^a^†a^†→(n+1)(n+2)\hat a\hat a\hat a^\dagger\hat a^\dagger\to(n+1)(n+2)、a^†a^†a^a^→n(n−1)\hat a^\dagger\hat a^\dagger\hat a\hat a\to n(n-1)、残りの4つ(a^a^†a^a^†\hat a\hat a^\dagger\hat a\hat a^\dagger・a^a^†a^†a^\hat a\hat a^\dagger\hat a^\dagger\hat a・a^†a^a^a^†\hat a^\dagger\hat a\hat a\hat a^\dagger・a^†a^a^†a^\hat a^\dagger\hat a\hat a^\dagger\hat a)は (n+1)2(n+1)^2・n(n+1)n(n+1)・n(n+1)n(n+1)・n2n^2 を与え、合計は 6n2+6n+3=3(2n2+2n+1)6n^2+6n+3=3(2n^2+2n+1) である。熱平衡の重みで平均すると、幾何分布の ⟨n⟩=nˉ\langle n \rangle =\bar n、⟨n2⟩=2nˉ2+nˉ\langle n^2 \rangle =2\bar n^2+\bar n から ⟨2n2+2n+1⟩=4nˉ2+4nˉ+1=(2nˉ+1)2\langle 2n^2+2n+1 \rangle =4\bar n^2+4\bar n+1=(2\bar n+1)^2 となり、⟨x4⟩0=3(ℏ2mω)2(2nˉ+1)2=3G(0)2\langle x^4 \rangle _0=3\bigl(\frac\hbar{2m\omega}\bigr)^2(2\bar n+1)^2=3G(0)^2 で (2) と一致する。ℏ=m=ω=1\hbar=m=\omega=1、s=1s=1 では 3.5123.512 である。

(4) ∫dy e−y2/(2σ2)(xˉ+y)4/2πσ2=xˉ4+6σ2xˉ2+3σ4\int\mathrm{d} y\,\mathrm{e}^{-y^2/(2\sigma^2)}(\bar x+y)^4/\sqrt{2\pi\sigma^2}=\bar x^4+6\sigma^2\bar x^2+3\sigma^4 である(奇数次は消え、⟨y2⟩=σ2\langle y^2 \rangle =\sigma^2、⟨y4⟩=3σ4\langle y^4 \rangle =3\sigma^4)。λ=0\lambda=0 の ZFZ_{\mathrm F} の重みは e−βmω2xˉ2/2\mathrm{e}^{-\beta m\omega^2\bar x^2/2} なので重心は分散 cc の Gauss 分布で、λ\lambda の1次の補正は λ(3c2+6σ2c+3σ4)=3λ(c+σ2)2\lambda(3c^2+6\sigma^2c+3\sigma^4)=3\lambda(c+\sigma^2)^2 である。ℏ=m=ω=1\hbar=m=\omega=1、s=1s=1 では c=1c=1、σ2=1/12\sigma^2=1/12 で 3(1+1/12)2=3.5213(1+1/12)^2=3.521、厳密な値は 3.5123.512、古典(σ=0\sigma=0)は 33 である。厳密な値は 3(c+f)23(c+f)^2 の形に書ける。ここで f=G(0)−c=1mβω2[s2coth⁡s2−1]f=G(0)-c=\frac1{m\beta\omega^2}\bigl[\frac s2\coth\frac s2-1\bigr] は調和ポテンシャルの中の本当の揺らぎの分散で、s=1s=1 では f=0.0820<σ2=0.0833f=0.0820<\sigma^2=0.0833 である。ぼかしは自由粒子の揺らぎを使うので、ポテンシャルが揺らぎを抑える効果を落とし、⟨x4⟩\langle x^4 \rangle をわずかに大きく見積もる。高温では f≈βℏ2/(12m)=σ2f\approx\beta\hbar^2/(12m)=\sigma^2 で両者が一致する。

(5) β→∞\beta\to\infty では nˉ→0\bar n\to0 で G(0)→ℏ/(2mω)G(0)\to\hbar/(2m\omega) なので、FF の1次の補正は 3λ(ℏ/(2mω))23\lambda(\hbar/(2m\omega))^2 に近づく。β→∞\beta\to\infty で F→E0F\to E_0(3.6.6節の例題)なので、これは基底エネルギーのずれである。(3) の n=0n=0 の値 λ⟨0∣x^4∣0⟩=3λ(ℏ/(2mω))2\lambda\langle 0 | \hat x^4 | 0 \rangle =3\lambda(\hbar/(2m\omega))^2 と一致し、1次の摂動論の基底エネルギーと同じものが経路積分から出た。(4) の近似は c→0c\to0 でも σ2=βℏ2/(12m)→∞\sigma^2=\beta\hbar^2/(12m)\to\infty なので、補正が 3λσ4→∞3\lambda\sigma^4\to\infty となって意味を失う。自由粒子を基準にしたぼかしは低温では使えず、Feynman–Kleinert の改良はこの欠点を直すためのものである。

木下の勉強ステータス

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

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

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