tecktom.com

この章のまとめ

2026年10月17日10分で読めますAI生成教科書『虎の巻』第7部 第5章 確率過程量子化:Nelson の拡散から Schrödinger 方程式へ

AI 生成・木下が未勉強

量子力学の様々な構成方法 第7部 第5章 §7.5.7 ★★★ Nelson の構成の出発点と同値性、共通例題①③の答え、熱平衡の扱い、Parisi–Wu との違いを1つにまとめ、総合問題で結ぶ

Nelson の確率過程量子化は、粒子の位置を拡散係数 ν=ℏ/(2m)\nu=\hbar/(2m) の Markov の拡散過程に置き、前向きと後ろ向きの平均微分を対称に組んだ加速度に Newton の法則を課す構成である。前向きと後ろ向きのドリフトの半差の浸透速度 u=ν∇ln⁡ρu=\nu\nabla\ln\rho は、その非線形項と拡散の項がまとめて量子ポテンシャルの力になり、半和の流れの速度 vv が第7部第1章の量子 Euler 方程式に従う。こうして Madelung 方程式が出て、循環の量子化を足せば Schrödinger 方程式と同値になる(定理 7.5.3)。量子化の条件が要請から出ないことは Wallstrom の指摘のとおりである。

共通例題①の基底状態は Ornstein–Uhlenbeck 過程 dX=−ωX dt+dW\mathrm{d} X=-\omega X\,\mathrm{d} t+\mathrm{d} W で、定常分散が ℏ/(2mω)\hbar/(2m\omega)、相関は σg2e−ω∣t∣\sigma_{\mathrm g}^2\mathrm{e}^{-\omega\lvert t\rvert } だった。励起状態の過程は節を越えず、コヒーレント状態の過程は古典の振動に Ornstein–Uhlenbeck の揺れを重ねる。③の波束では、t∗=2mσ02/ℏt_*=2m\sigma_0^2/\hbar より前の内向きのドリフトが拡散を打ち消し、後の外向きのドリフトが粒子を運んで、正答の σ(t)2\sigma(t)^2 が出る。熱平衡は1つの定常な過程では書けず(命題 7.5.4)、分布を合わせた過程は柔らかい振動子の基底状態を表してしまう。混合状態は過程の混合として書くが、道筋はアンサンブルの選び方で変わる。

Parisi–Wu の確率過程量子化は、仮想時間の Langevin 方程式の平衡で虚時間の経路積分の重みを作る計算法で、Nelson の実時間の過程とは時間の意味も確率変数も違う。ただし Nelson の基底状態の過程は、基底状態の変換を通して虚時間の時間発展と同じ確率法則を持つ(命題 7.5.5)。この構成の限界のうち、多時刻の相関が基準と食い違うこと・多粒子の過程が配置空間に住むこと・節の特異性は、Feynman–Kac の公式と合わせて第7部第6章で扱う。

総合問題 熱平衡の調和振動子を、3つの確率過程で比べる ★★★

オリジナル1次元の調和振動子(質量 mm、角周波数 ω\omega、σg2=ℏ/(2mω)\sigma_{\mathrm g}^2=\hbar/(2m\omega))の温度 TT の熱平衡を考え、nˉ=1/(eβℏω−1)\bar n=1/(\mathrm{e}^{\beta\hbar\omega}-1) と書く。

  1. (1) Parisi–Wu の方法の平衡分布は虚時間の経路積分の重み e−SE/ℏ\mathrm{e}^{-S_{\mathrm E}/\hbar} である。例題 7.5.5 のモード分解を使い、平衡での等虚時間の分散⟨x(τ)2⟩\langle x(\tau)^2 \rangle を Matsubara 和として求め、σg2coth⁡(βℏω/2)\sigma_{\mathrm g}^2\coth(\beta\hbar\omega/2) になることを示せ。(使う道具:この章の第5節・第3部第6章 の虚時間の経路積分)
  2. (2) 同じ分散を定常分布に持つ Ornstein–Uhlenbeck 過程 dX=−γX dt+dW\mathrm{d} X=-\gamma X\,\mathrm{d} t+\mathrm{d} W(⟨dW2⟩=2ν dt\langle \mathrm{d} W^2 \rangle =2\nu\,\mathrm{d} t、ν=ℏ/(2m)\nu=\hbar/(2m))の γ\gamma を求め、これが T>0T>0 で Nelson の法則 ma=−mω2xma=-m\omega^2x を満たさないことを示せ。(使う道具:第2節の (7.5.11)・第4節)
  3. (3) (2) の過程の「エネルギー」E=∫ρ[m2(v2+u2)+V] dxE=\int\rho[\tfrac m2(v^2+u^2)+V]\,\mathrm{d} x を計算し、熱平衡のエネルギー ⟨H^⟩=ℏω(nˉ+12)\langle \hat H \rangle =\hbar\omega(\bar n+\tfrac12) と比べよ。(使う道具:例題 7.5.2・第4部第4章)
  4. (4) コヒーレント状態の過程の混合(第4節の (B))で、1時刻の分散と、成分ごとの EE を P(α)P(\alpha) で平均したものを求め、どちらも熱平衡と一致することを示せ。(使う道具:第3節の (7.5.20)・第7部第4章 の総合問題の P 表示の混合)
  5. (5) kBTk_{\mathrm B} T を固定して ℏ→0\hbar\to0 とすると、(B) の過程の道筋は古典の熱平衡の振動子の集団(Boltzmann 分布で振幅がばらつく決定論的な振動)に近づくことを示せ。(使う道具:第6節の ν→0\nu\to0・第7部第2章)
解答

(1) 平衡で各モードは ⟨∣cn∣2⟩=ℏ/(m(ωn2+ω2))\langle \lvert c_n\rvert ^2 \rangle =\hbar/(m(\omega_n^2+\omega^2))。⟨x(τ)2⟩=(βℏ)−1∑n⟨∣cn∣2⟩\langle x(\tau)^2 \rangle =(\beta\hbar)^{-1}\sum_n\langle \lvert c_n\rvert ^2 \rangle(ei(ωn−ωn′)τ\mathrm{e}^{\mathrm{i}(\omega_n-\omega_{n'})\tau} の交差項は平均で 00)なので

⟨x(τ)2⟩=1mβ∑n=−∞∞1ωn2+ω2.\langle x(\tau)^2 \rangle =\frac{1}{m\beta}\sum_{n=-\infty}^{\infty}\frac{1}{\omega_n^2+\omega^2}.

ωn=2πn/(βℏ)\omega_n=2\pi n/(\beta\hbar)、a=βℏω/(2π)a=\beta\hbar\omega/(2\pi) と置くと和は (βℏ/2π)2∑n1/(n2+a2)(\beta\hbar/2\pi)^2\sum_n1/(n^2+a^2) で、公式 ∑n1/(n2+a2)=(π/a)coth⁡(πa)\sum_n1/(n^2+a^2)=(\pi/a)\coth(\pi a) を使うと

∑n1ωn2+ω2=(βℏ2π)22π2βℏωcoth⁡βℏω2=βℏ2ωcoth⁡βℏω2,⟨x2⟩=ℏ2mωcoth⁡βℏω2=σg2coth⁡βℏω2.\begin{aligned}\sum_n\frac{1}{\omega_n^2+\omega^2}&=\Bigl(\frac{\beta\hbar}{2\pi}\Bigr)^2\frac{2\pi^2}{\beta\hbar\omega}\coth\frac{\beta\hbar\omega}{2}=\frac{\beta\hbar}{2\omega}\coth\frac{\beta\hbar\omega}{2},\\ \langle x^2 \rangle &=\frac{\hbar}{2m\omega}\coth\frac{\beta\hbar\omega}{2}=\sigma_{\mathrm g}^2\coth\frac{\beta\hbar\omega}{2}.\end{aligned}

Parisi–Wu の方法は熱平衡の答えを正しく出す。これは虚時間の経路積分が混合状態 e−βH^\mathrm{e}^{-\beta\hat H} をそのまま扱えるからで、実時間の道筋は何も主張していない。

(2) 定常分散は ν/γ\nu/\gamma なので γ=ν/(σg2coth⁡)=ωtanh⁡(βℏω/2)\gamma=\nu/(\sigma_{\mathrm g}^2\coth)=\omega\tanh(\beta\hbar\omega/2)。v=0v=0、u=−γxu=-\gamma x で a=−uu′−νu′′=−γ2xa=-uu'-\nu u''=-\gamma^2x。T>0T>0 では tanh⁡<1\tanh<1 で γ<ω\gamma<\omega なので a≠−ω2xa\ne-\omega^2x。命題 7.5.4 の具体例である。

(3) ρ\rho は分散 s2=ν/γs^2=\nu/\gamma のガウス分布で、m2⟨u2⟩=m2γ2s2=mνγ2\tfrac m2\langle u^2 \rangle =\tfrac m2\gamma^2s^2=\tfrac{m\nu\gamma}{2}、⟨V⟩=mω2s22=mνω22γ\langle V \rangle =\tfrac{m\omega^2s^2}{2}=\tfrac{m\nu\omega^2}{2\gamma}。mν=ℏ/2m\nu=\hbar/2 なので

E=ℏ4(γ+ω2γ)=ℏω4(tanh⁡βℏω2+coth⁡βℏω2)=ℏω2coth⁡(βℏω)E=\frac{\hbar}{4}\Bigl(\gamma+\frac{\omega^2}{\gamma}\Bigr)=\frac{\hbar\omega}{4}\Bigl(\tanh\frac{\beta\hbar\omega}{2}+\coth\frac{\beta\hbar\omega}{2}\Bigr)=\frac{\hbar\omega}{2}\coth(\beta\hbar\omega)

(tanh⁡x+coth⁡x=2coth⁡2x\tanh x+\coth x=2\coth2x)。熱平衡は ℏω(nˉ+12)=ℏω2coth⁡(βℏω/2)\hbar\omega(\bar n+\tfrac12)=\tfrac{\hbar\omega}{2}\coth(\beta\hbar\omega/2) で、T>0T>0 では coth⁡(βℏω)<coth⁡(βℏω/2)\coth(\beta\hbar\omega)<\coth(\beta\hbar\omega/2) なので小さい。たとえば βℏω=1\beta\hbar\omega=1 で 0.657ℏω0.657\hbar\omega と 1.082ℏω1.082\hbar\omega。位置の分布を合わせても、運動量の広がり(浸透速度)が足りない。

(4) X=Xc+ξX=X_{\mathrm c}+\xi で、P(α)P(\alpha) のもとで Xc=2σgRe⁡αX_{\mathrm c}=2\sigma_{\mathrm g}\operatorname{Re}\alpha の分散は 2nˉσg22\bar n\sigma_{\mathrm g}^2、ξ\xi の分散は σg2\sigma_{\mathrm g}^2、独立なので和は σg2(1+2nˉ)=σg2coth⁡(βℏω/2)\sigma_{\mathrm g}^2(1+2\bar n)=\sigma_{\mathrm g}^2\coth(\beta\hbar\omega/2)。成分のエネルギーは v=Pc/mv=P_{\mathrm c}/m、u=−ωξu=-\omega\xi からEα=Pc22m+m2ω2σg2+mω22(Xc2+σg2)=Pc22m+mω2Xc22+ℏω2=ℏω(∣α∣2+12)E_\alpha=\frac{P_{\mathrm c}^2}{2m}+\frac m2\omega^2\sigma_{\mathrm g}^2+\frac{m\omega^2}{2}(X_{\mathrm c}^2+\sigma_{\mathrm g}^2)=\frac{P_{\mathrm c}^2}{2m}+\frac{m\omega^2X_{\mathrm c}^2}{2}+\frac{\hbar\omega}{2}=\hbar\omega\bigl(\lvert \alpha\rvert ^2+\tfrac12\bigr)(mω2σg2=ℏω/2m\omega^2\sigma_{\mathrm g}^2=\hbar\omega/2)。⟨∣α∣2⟩P=nˉ\langle \lvert \alpha\rvert ^2 \rangle _P=\bar n なので平均は ℏω(nˉ+12)\hbar\omega(\bar n+\tfrac12) で、熱平衡と一致する。1つの過程では合わなかった分布とエネルギーが、過程の混合では両方合う。

(5) ℏ→0\hbar\to0 で σg2=ℏ/(2mω)→0\sigma_{\mathrm g}^2=\hbar/(2m\omega)\to0 なので揺れ ξ\xi の幅は 00 になり、道筋は X=Xc(t)X=X_{\mathrm c}(t) の決定論的な古典の振動に近づく。振幅の分散は 2nˉσg2→(2kBT/ℏω)⋅ℏ/(2mω)=kBT/(mω2)2\bar n\sigma_{\mathrm g}^2\to(2 k_{\mathrm B} T/\hbar\omega)\cdot\hbar/(2m\omega)= k_{\mathrm B} T/(m\omega^2) で、XcX_{\mathrm c} と Pc/(mω)P_{\mathrm c}/(m\omega) は独立なガウス分布なので、(Xc,Pc)(X_{\mathrm c},P_{\mathrm c}) の分布は exp⁡[−(Pc2/2m+mω2Xc2/2)/kBT]\exp[-(P_{\mathrm c}^2/2m+m\omega^2X_{\mathrm c}^2/2)/ k_{\mathrm B} T] に比例する Boltzmann 分布になる。(A) の固有状態の混合では、どれほど ℏ\hbar が小さくても各成分は定常で節に区切られた過程のままで、古典の振動は現れない(第7部第4章の総合問題の Bohm の軌道と同じ結論)。

木下の勉強ステータス

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

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

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