tecktom.com

橋:同じ量をほかの構成で書く

2026年10月17日9分で読めますAI生成教科書『虎の巻』第5部 第1章 量子力学を「確率分布」で書けるか
目次
  1. 期待値の書き方を並べる
  2. 分配関数を4つの構成で並べる
  3. 学部の統計力学の書き方との対応
  4. 木下の層

AI 生成・木下が未勉強

量子力学の様々な構成方法 第5部 第1章 §5.1.6 ★★★ WW の1次のモーメントは確率流、ZZ は4つの構成で同じ ℏ2\hbar^2 補正を持つ。第5部の位置を地図の上で確かめる

この章の結果は、計画の同値性の地図で「密度行列 ↔ 位相空間」の橋の手前にある。橋そのもの(WW と ρ^\hat\rho が互いに一意に決まること、Tr⁡(ρ^A^)=∫WAW\operatorname{Tr}(\hat\rho\hat A)=\int W A_W)は次の章で両向きに証明する。この節では、この章で出した量がほかの構成でどう書けるかを式で並べる。

期待値の書き方を並べる

位置だけ・運動量だけの関数については、周辺分布の性質 (W3) からすでに次の等式が成り立っている。

Tr⁡(ρ^f(x^))=∫f(x) W dx dp,Tr⁡(ρ^g(p^))=∫g(p) W dx dp.(5.1.15)\operatorname{Tr}\bigl(\hat\rho f(\hat x)\bigr)=\int f(x)\,W\,\mathrm{d} x\,\mathrm{d} p,\qquad \operatorname{Tr}\bigl(\hat\rho g(\hat p)\bigr)=\int g(p)\,W\,\mathrm{d} x\,\mathrm{d} p .\tag{5.1.15}

xx と pp が混ざる量で一番簡単なのは pp の1次のモーメントである。∫p e−2ipy/ℏdp=(iℏ/2)∂y∫e−2ipy/ℏdp=(iℏ/2)πℏ δ′(y)\int p\,\mathrm{e}^{-2\mathrm{i} py/\hbar}\mathrm{d} p=(\mathrm{i}\hbar/2)\partial_y\int\mathrm{e}^{-2\mathrm{i} py/\hbar}\mathrm{d} p=(\mathrm{i}\hbar/2)\pi\hbar\,\delta'(y) を使って部分積分すると

∫p W(x,p) dp=−iℏ2∂∂y[ψ(x+y)ψ∗(x−y)]y=0=ℏIm⁡(ψ∗ψ′)=m j(x)(5.1.16)\int p\,W(x,p)\,\mathrm{d} p=-\frac{\mathrm{i}\hbar}{2}\frac{\partial }{\partial y}\bigl[\psi(x+y)\psi^*(x-y)\bigr]_{y=0} =\hbar\operatorname{Im}\bigl(\psi^*\psi'\bigr)=m\,j(x)\tag{5.1.16}

となる。j=(ℏ/m)Im⁡(ψ∗ψ′)j=(\hbar/m)\operatorname{Im}(\psi^*\psi') は基準の Schrödinger 方程式から出る確率流である。WW の各位置での運動量の平均は、確率流を密度で割った速度場 j/∣ψ∣2j/\lvert \psi\rvert ^2 に mm を掛けたものになっている。これは第7部の Madelung の流体(第7部第1章)の速度場 v=∂xS/mv=\partial_xS/m そのものである(次の例題)。(5.1.16) にさらに xx を掛けて積分すると、∫xp W dx dp=ℏ∫xIm⁡(ψ∗ψ′)dx=Re⁡⟨x^p^⟩=12⟨x^p^+p^x^⟩\int xp\,W\,\mathrm{d} x\,\mathrm{d} p=\hbar\int x\operatorname{Im}(\psi^*\psi')\mathrm{d} x=\operatorname{Re}\langle \hat x\hat p \rangle =\tfrac12\langle \hat x\hat p+\hat p\hat x \rangleで、例題 5.1.3 で見た対称化した積の期待値が一般に出る。並べると次のとおりである。

量基準(状態ベクトル・密度行列)位相空間(WW)
位置の分布⟨x∣ρ^∣x⟩\langle x | \hat\rho | x \rangle∫W dp\int W\,\mathrm{d} p
運動量の分布⟨p∣ρ^∣p⟩\langle p | \hat\rho | p \rangle∫W dx\int W\,\mathrm{d} x
確率流(ℏ/m)Im⁡(ψ∗ψ′)(\hbar/m)\operatorname{Im}(\psi^*\psi')(1/m)∫pW dp(1/m)\int pW\,\mathrm{d} p
対称化した積12⟨x^p^+p^x^⟩\tfrac12\langle \hat x\hat p+\hat p\hat x \rangle∫xpW dx dp\int xpW\,\mathrm{d} x\,\mathrm{d} p
原点での値⟨ψ∣Π^∣ψ⟩/(πℏ)\langle \psi | \hat\Pi | \psi \rangle /(\pi\hbar)W(0,0)W(0,0)

例題 5.1.6 極形式で書いた運動量の平均 ★★

オリジナルψ(x)=ρ(x) eiS(x)/ℏ\psi(x)=\sqrt{\rho(x)}\,\mathrm{e}^{\mathrm{i} S(x)/\hbar}(ρ≥0\rho\ge0、SS は実数値)と書く。

  1. (1) ∫pW dp=ρ ∂xS\int pW\,\mathrm{d} p=\rho\,\partial_xS を示せ。
  2. (2) 位置 xx での条件付きの運動量の平均 pˉ(x)=∫pW dp/∫W dp\bar p(x)=\int pW\,\mathrm{d} p\big/\int W\,\mathrm{d} p を求め、例題 5.1.2 の波束で確かめよ。
解答

(1) ψ′=(ρ′2ρ+iS′ℏρ)eiS/ℏ\psi'=\bigl(\tfrac{\rho'}{2\sqrt\rho}+\tfrac{\mathrm{i} S'}{\hbar}\sqrt\rho\bigr)\mathrm{e}^{\mathrm{i} S/\hbar} なのでψ∗ψ′=ρ′/2+iρS′/ℏ\psi^*\psi'=\rho'/2+\mathrm{i}\rho S'/\hbar。虚部は ρS′/ℏ\rho S'/\hbar で、(5.1.16) から ∫pW dp=ℏ⋅ρS′/ℏ=ρ ∂xS\int pW\,\mathrm{d} p=\hbar\cdot\rho S'/\hbar=\rho\,\partial_xS。

(2) ∫W dp=ρ\int W\,\mathrm{d} p=\rho なので pˉ(x)=∂xS\bar p(x)=\partial_xS。速度に直すと ∂xS/m\partial_xS/m で、これが Madelung の速度場である。例題 5.1.2 の波束では S=ℏκx2S=\hbar\kappa x^2 なので pˉ=2ℏκx\bar p=2\hbar\kappa x。例題 5.1.3 の WW は指数に(p−2ℏκx)2(p-2\hbar\kappa x)^2 を含み、各 xx で pp の分布の中心が 2ℏκx2\hbar\kappa x にあるので、一致する。

分配関数を4つの構成で並べる

5節の分配関数は、ほかの構成でも計算できる。並べると次のようになる。

Z=Tr⁡e−βH^密度行列(第4部第4章)=∫dx dp2πℏ e−βH w(x,p,β)位相空間(この章 5節)=∮Dx e−SE[x]/ℏ虚時間の経路積分(第3部第6章)→∫dx dp2πℏ e−βH古典(ℏ→0、解析力学 第3部第3章)\begin{align*}Z&=\operatorname{Tr}\mathrm{e}^{-\beta\hat H} &&\text{密度行列(第4部第4章)}\tag{5.1.17}\\&=\int\frac{\mathrm{d} x\,\mathrm{d} p}{2\pi\hbar}\,\mathrm{e}^{-\beta H}\,w(x,p,\beta) &&\text{位相空間(この章 5節)}\tag{5.1.18}\\&=\oint\mathcal Dx\,\mathrm{e}^{-S_{\mathrm E}[x]/\hbar} &&\text{虚時間の経路積分(第3部第6章)}\tag{5.1.19}\\&\to\int\frac{\mathrm{d} x\,\mathrm{d} p}{2\pi\hbar}\,\mathrm{e}^{-\beta H} &&\text{古典(}{\hbar\to0}\text{、解析力学 第3部第3章)}\tag{5.1.20}\end{align*}

経路積分 (5.1.19) は周期 βℏ\beta\hbar の閉じた経路の和である。経路を、その平均の位置 xˉ\bar x と、まわりの揺らぎに分ける。揺らぎは ℏβ/m\hbar\beta/m 程度の分散を持つ Gauss 型で、VV を揺らぎについて2次まで展開して揺らぎで平均すると、R. P. Feynman が導いた「ぼかしたポテンシャル」

Veff(xˉ)=V(xˉ)+ℏ2β24mV′′(xˉ)(5.1.21)V_{\mathrm{eff}}(\bar x)=V(\bar x)+\frac{\hbar^2\beta}{24m}V''(\bar x)\tag{5.1.21}

を古典の ZZ に入れたものが ℏ2\hbar^2 の次数まで正しい ZZ になる(揺らぎの平均 ⟨δx2⟩=ℏ2β/(12m)\langle \delta x^2 \rangle =\hbar^2\beta/(12m) から 12V′′⟨δx2⟩\tfrac12V''\langle \delta x^2 \rangle が出る。導出は第3部第6章で行う)。これを展開すると

Z≈Zcl(1−β⋅ℏ2β24m⟨V′′⟩cl)=Zcl(1−ℏ2β324m⟨V′2⟩cl)Z\approx Z_{\mathrm{cl}}\Bigl(1-\beta\cdot\frac{\hbar^2\beta}{24m}\langle V'' \rangle _{\mathrm{cl}}\Bigr) =Z_{\mathrm{cl}}\Bigl(1-\frac{\hbar^2\beta^3}{24m}\langle V'^2 \rangle _{\mathrm{cl}}\Bigr)

で、2つ目の等号は部分積分 (5.1.13)。定理 5.1.6 と同じ補正が出る。位相空間では補正が被積分関数の因子 ww に、経路積分では補正がポテンシャルのぼかしに現れ、積分した結果は一致する。図 5.1.3 はこの関係をまとめたものである。

図

図 5.1.3 分配関数の4つの書き方。上から下への矢印はどれも ℏ2\hbar^2 の補正 −ℏ2β3⟨V′2⟩cl/(24m)-\hbar^2\beta^3\langle V'^2 \rangle _{\mathrm{cl}}/(24m) を経て古典に至る。

学部の統計力学の書き方との対応

学部の統計力学の授業では、古典の分配関数を Zcl=1h∫e−βHdq dpZ_{\mathrm{cl}}=\frac1h\int\mathrm{e}^{-\beta H}\mathrm{d} q\,\mathrm{d} p と書き、1/h1/h は「量子力学との整合のために入れる因子」として説明されることが多い。この章の言葉では、1/h=1/(2πℏ)1/h=1/(2\pi\hbar) は(5.1.9) で平面波の完全性 ∫dp ∣p⟩⟨p∣=1\int\mathrm{d} p\,| p \rangle \langle p |=1 を位置の積分と組み合わせたときに出る因子で、4節の「1区画の面積 2πℏ2\pi\hbar」と同じものである。また、量子の調和振動子の比熱が低温で古典の値 kBk_{\mathrm B} から下がることは、学部では固有値の和 ∑ne−βEn\sum_n\mathrm{e}^{-\beta E_n} から計算する。位相空間の構成では同じことが、ww の補正が βℏω\beta\hbar\omega の大きい所で大きくなること(図 5.1.2)として見える。

木下の勉強ステータス

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

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

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