tecktom.com

この章のまとめ

2026年10月17日8分で読めますAI生成教科書『虎の巻』第7部 第2章 ℏ→0 で Hamilton–Jacobi へ:流体から見た古典極限

AI 生成・木下が未勉強

量子力学の様々な構成方法 第7部 第2章 §7.2.6 ★★★ QQ を無視した流れが古典の HJ 方程式と連続の式で、その成り立つ範囲が焦点と定在波で区切られることを振り返る

Madelung の流体で ℏ\hbar を含むのは量子ポテンシャル QQ だけで、ℏ→0\hbar\to0 とは ρ\rho と SS を固定したまま de Broglie 波長を密度の変わる長さより短くする極限である。QQ を無視すると、位相 SS は古典の Hamilton–Jacobi 方程式に、密度は連続の式に従い、その流れは各点で運動量が ∇S\nabla S に決まった古典粒子の集団を表す。特性曲線は Hamilton 方程式の解で、密度は軌道の Jacobi 行列式で運ばれる。逆向きに、古典の解から作った ρ eiS/ℏ\sqrt\rho\,\mathrm{e}^{\mathrm{i} S/\hbar} は Schrödinger 方程式を残差 −Qψ-Q\psi で満たし、本物の解との差は焦点より前の時刻で O(ℏt)\mathcal{O}(\hbar t) に収まる。時間によらない流れでは ρ∝1/p\rho\propto1/p で、これは WKB 近似の主要項であり、無視した QQ と運動エネルギーの比が一様な力では (5/16π2)(dλ/dx)2(5/16\pi^2)(\mathrm{d}\lambda/\mathrm{d} x)^2 になる。

極限が一様でないのは2つの場合である。軌道が交わる焦点では古典の密度が発散し、量子ポテンシャルのエネルギーが古典のエネルギーと同じ大きさになる。束縛状態の定在波では2つの流れが重なって密度が波長の長さで振動するので、Q=En−VQ=E_n-V は小さくならず、局所的に平均した密度だけが古典の滞在時間の密度に近づく。位相空間の言葉では、単一運動量の分布の Liouville 方程式が流体の古典極限と同値で(定理 7.2.3)、経路積分の言葉では Van Vleck の前因子の2乗が古典の密度になる。3つの見方は同じ古典軌道の族を使うので、同じ焦点で崩れる。焦点を越えて流れを続ける処方は経路積分の停留位相にあり、流体の2本の式には無い。

総合問題 1点に集まる速度場を持たせた原子の波束 ★★★

オリジナル自由な原子(質量 mm、V=0V=0)の1次元の波束を、時刻 00 にψ0(x)=(2πσ02)−1/4exp⁡(−x2/(4σ02)−imx2/(2ℏT0))\psi_0(x)=(2\pi\sigma_0^2)^{-1/4}\exp\bigl(-x^2/(4\sigma_0^2)-\mathrm{i} mx^2/(2\hbar T_0)\bigr)(T0>0T_0>0)と用意する。位相は S0=−mx2/(2T0)S_0=-mx^2/(2T_0) で、速度場 v0=−x/T0v_0=-x/T_0 は原点へ向かう流れである。ε=ℏT0/(2mσ02)\varepsilon=\hbar T_0/(2m\sigma_0^2) と置く。

  1. (1) QQ を無視した流れについて、S(x,t)=−mx2/(2(T0−t))S(x,t)=-mx^2/(2(T_0-t)) が Hamilton–Jacobi 方程式と初期条件を満たすことを確かめ、軌道と古典の幅 σcl(t)\sigma_{\mathrm{cl}}(t) を求めよ。(使う道具:第2節)
  2. (2) 量子の幅は σ¨=ℏ2/(4m2σ3)\ddot\sigma=\hbar^2/(4m^2\sigma^3) に σ(0)=σ0\sigma(0)=\sigma_0、σ˙(0)=−σ0/T0\dot\sigma(0)=-\sigma_0/T_0 を与えた解である。σ(t)2=σ02(1−t/T0)2+(ℏt/(2mσ0))2\sigma(t)^2=\sigma_0^2(1-t/T_0)^2+(\hbar t/(2m\sigma_0))^2 がその解であることを確かめよ。(使う道具:第7部第1章の量子 Euler 方程式・第2節の例題)
  3. (3) 幅が最小になる時刻 t∗t_* と最小の幅 σmin⁡\sigma_{\min} を ε\varepsilon で表せ。87{}^{87}Rb 原子(m=1.443×10−25 kgm=1.443\times10^{-25}\,\mathrm{kg})で σ0=1 mm\sigma_0=1\,\mathrm{mm}、T0=0.1 sT_0=0.1\,\mathrm s のときの値を求めよ。(使う道具:第4節)
  4. (4) 運動エネルギーの期待値は流れのエネルギーと ∫ρQ=ℏ2/(8mσ2)\int\rho Q=\hbar^2/(8m\sigma^2) の和である。時刻 00 と t∗t_* で、全運動エネルギーのうち ∫ρQ\int\rho Q が占める割合を求めよ。(使う道具:第1節・第4節・第7部第1章の Fisher 情報)
  5. (5) 古典粒子の集団を、位置 x0x_0 を平均 00・分散 σ02\sigma_0^2 のガウス分布で、運動量を p0=−mx0/T0+ηp_0=-mx_0/T_0+\eta(η\eta は x0x_0 と独立で平均 00・分散 ℏ2/(4σ02)\hbar^2/(4\sigma_0^2))で用意し、自由に運動させる。位置の分散を求めて (2) と比べ、この初めの分布が ψ0\psi_0 の Wigner 関数であることを確かめよ。結果から、QQ を無視することを位相空間の言葉で言い直せ。(使う道具:第5節・第5部第3章)
解答

(1) ∂tS=−mx2/(2(T0−t)2)\partial_tS=-mx^2/(2(T_0-t)^2)、∂xS=−mx/(T0−t)\partial_xS=-mx/(T_0-t) で、∂tS+(∂xS)2/(2m)=0\partial_tS+(\partial_xS)^2/(2m)=0。t=0t=0 で S0S_0 に戻る。特性曲線は x˙=−x/(T0−t)\dot x=-x/(T_0-t) で、解は x(t)=x0(1−t/T0)x(t)=x_0(1-t/T_0) である(速度 −x0/T0-x_0/T_0 の等速運動)。Jacobi 行列式は J=1−t/T0J=1-t/T_0 なので、ガウス分布の幅はσcl(t)=σ0∣1−t/T0∣\sigma_{\mathrm{cl}}(t)=\sigma_0\lvert 1-t/T_0\rvert で、t=T0t=T_0 ですべての軌道が原点で交わり、密度が発散する。SS も t→T0t\to T_0 で発散し、その先で1価の解は続かない。

(2) σ2\sigma^2 を tt で微分すると 2σσ˙=−2σ02(1−t/T0)/T0+ℏ2t/(2m2σ02)2\sigma\dot\sigma=-2\sigma_0^2(1-t/T_0)/T_0+\hbar^2t/(2m^2\sigma_0^2) で、t=0t=0 で σ˙(0)=−σ0/T0\dot\sigma(0)=-\sigma_0/T_0 である。σ¨=ℏ2/(4m2σ3)\ddot\sigma=\hbar^2/(4m^2\sigma^3) を満たすことは、2回微分して代入し、検算で確かめた。第2節の例題の τ\tau を −T0-T_0 に置き換えた形で、σ2=σcl2+(ℏt/2mσ0)2\sigma^2=\sigma_{\mathrm{cl}}^2+(\hbar t/2m\sigma_0)^2 と古典の幅に量子圧力の項が足されている。

(3) dσ2/dt=0\mathrm{d}\sigma^2/\mathrm{d} t=0 から t∗(σ02/T02+ℏ2/(4m2σ02))=σ02/T0t_*\bigl(\sigma_0^2/T_0^2+\hbar^2/(4m^2\sigma_0^2)\bigr)=\sigma_0^2/T_0、つまり

t∗=T01+ε2,σmin⁡2=σ02ε21+ε2.t_*=\frac{T_0}{1+\varepsilon^2},\qquad \sigma_{\min}^2=\sigma_0^2\frac{\varepsilon^2}{1+\varepsilon^2}.

数値は ε=1.0546×10−34×0.1/(2×1.443×10−25×10−6)=3.65×10−5\varepsilon=1.0546\times10^{-34}\times0.1/(2\times1.443\times10^{-25}\times10^{-6})=3.65\times10^{-5} で、t∗t_* は T0T_0 と区別できず、σmin⁡=σ0ε/1+ε2=36.5 nm\sigma_{\min}=\sigma_0\varepsilon/\sqrt{1+\varepsilon^2}=36.5\,\mathrm{nm} である。幅 1 mm1\,\mathrm{mm} の原子の集団が、古典の流れなら1点に集まる時刻に、量子の幅は 36.5 nm36.5\,\mathrm{nm} で止まる。ε\varepsilon は端の運動量 mσ0/T0m\sigma_0/T_0 と長さ σ0\sigma_0 で作った ℏ/(2pL)\hbar/(2pL) で、(7.2.3) の小さい比に当たる。

(4) 運動エネルギーの期待値は自由粒子では保存し、時刻 00 で ⟨p^2⟩/2m=mσ02/(2T02)+ℏ2/(8mσ02)\langle \hat p^2 \rangle /2m=m\sigma_0^2/(2T_0^2)+\hbar^2/(8m\sigma_0^2) である(流れの分と ∫ρQ\int\rho Q の分)。時刻 00 の割合は

ℏ2/(8mσ02)mσ02/(2T02)+ℏ2/(8mσ02)=ε21+ε2≈1.3×10−9\frac{\hbar^2/(8m\sigma_0^2)}{m\sigma_0^2/(2T_0^2)+\hbar^2/(8m\sigma_0^2)}=\frac{\varepsilon^2}{1+\varepsilon^2}\approx1.3\times10^{-9}

である。t∗t_* では σ˙=0\dot\sigma=0 で流れが止まるので、全運動エネルギーが ∫ρQ\int\rho Q になり、割合は 11 である。実際、ε2=ℏ2T02/(4m2σ04)\varepsilon^2=\hbar^2T_0^2/(4m^2\sigma_0^4) を使うと

ℏ28mσmin⁡2=ℏ2(1+ε2)8mσ02ε2=mσ02(1+ε2)2T02\frac{\hbar^2}{8m\sigma_{\min}^2}=\frac{\hbar^2(1+\varepsilon^2)}{8m\sigma_0^2\varepsilon^2}=\frac{m\sigma_0^2(1+\varepsilon^2)}{2T_0^2}

で、全運動エネルギーに等しい。初めは 10−910^{-9} の割合で無視できた量子ポテンシャルが、焦点では運動エネルギーのすべてを占める。第4節の調和ポテンシャルの例と同じく、焦点の近くでは QQ を無視できない。

(5) x(t)=x0+p0t/m=x0(1−t/T0)+ηt/mx(t)=x_0+p_0t/m=x_0(1-t/T_0)+\eta t/m で、x0x_0 と η\eta が独立なので

⟨x(t)2⟩=σ02(1−tT0)2+ℏ24σ02t2m2\langle x(t)^2\rangle=\sigma_0^2\Bigl(1-\frac{t}{T_0}\Bigr)^2+\frac{\hbar^2}{4\sigma_0^2}\frac{t^2}{m^2}

となり、(2) の量子の σ(t)2\sigma(t)^2 と一致する。初めの分布は

f0(x,p)=e−x2/(2σ02)2π σ0⋅e−2σ02(p+mx/T0)2/ℏ22π ℏ/(2σ0)=1πℏ e−x2/(2σ02) e−2σ02(p+mx/T0)2/ℏ2f_0(x,p)=\frac{\mathrm{e}^{-x^2/(2\sigma_0^2)}}{\sqrt{2\pi}\,\sigma_0}\cdot\frac{\mathrm{e}^{-2\sigma_0^2(p+mx/T_0)^2/\hbar^2}}{\sqrt{2\pi}\,\hbar/(2\sigma_0)} =\frac{1}{\pi\hbar}\,\mathrm{e}^{-x^2/(2\sigma_0^2)}\,\mathrm{e}^{-2\sigma_0^2(p+mx/T_0)^2/\hbar^2}

で、前の章の例題のチャープしたガウス波束の Wigner 関数(c=−m/(2ℏT0)c=-m/(2\hbar T_0) として 2ℏc x=−mx/T02\hbar c\,x=-mx/T_0)と同じである。自由粒子では Moyal 方程式が Liouville 方程式と一致するので、この一致は偶然でなく、Wigner 関数を古典の流れで運んだ結果が厳密に量子の答えになっている。位相空間で見ると、ψ0\psi_0 は各点 xx で運動量 −mx/T0-mx/T_0 のまわりに幅 ℏ/(2σ0)\hbar/(2\sigma_0) の広がりを持つ。QQ を無視することはこの広がり η\eta を 00 にして、分布を単一運動量の ρ δ(p+mx/T0)\rho\,\delta(p+mx/T_0) に置き換えることで(定理 7.2.3)、そうすると全粒子が t=T0t=T_0 で1点に集まる。位相空間で見ると、量子圧力は各点での運動量の広がりに当たる。

木下の勉強ステータス

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

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

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