tecktom.com

出てくるもの (2):波動関数と古典の分布

2026年10月15日8分で読めますAI生成教科書『虎の巻』第4部 第4章 調和振動子
目次
  1. 基底状態は1階の方程式で決まる
  2. 励起状態
  3. 古典の分布との比較
  4. 木下の層

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第4部 第4章 §4.4.4 ★★ a^ψ0=0\hat a\psi_0=0 が1階の微分方程式になり Gauss 関数が出る。nn が大きいと ∣ψn∣2\lvert \psi_n\rvert ^2 を平均したものは古典の滞在時間の分布に近づく

基底状態は1階の方程式で決まる

位置表示(q^→q\hat q\to q、p^→−iℏ ∂q\hat p\to-\mathrm{i}\hbar\,\partial_q。第4部第2章)で、a^∣0⟩=0\hat a| 0 \rangle =0 は

mω2ℏ(q+ℏmωddq)ψ0(q)=0,dψ0dq=−qℓ2 ψ0\sqrt{\frac{m\omega}{2\hbar}}\Bigl(q+\frac{\hbar}{m\omega}\frac{\mathrm{d} }{\mathrm{d} q}\Bigr)\psi_0(q)=0,\qquad \frac{\mathrm{d} \psi_0}{\mathrm{d} q}=-\frac{q}{\ell^2}\,\psi_0

という1階の方程式になる。2階の Schrödinger 方程式を解く代わりに、1階の方程式を積分するだけでよく、規格化して

ψ0(q)=1π1/4ℓ e−q2/(2ℓ2)(4.4.21)\psi_0(q)=\frac{1}{\pi^{1/4}\sqrt\ell}\,\mathrm{e}^{-q^2/(2\ell^2)}\tag{4.4.21}

を得る。∣ψ0∣2\lvert \psi_0\rvert ^2 は分散 ℓ2/2\ell^2/2 の Gauss 分布で、4.4.3節で不確定性から求めた最小の幅と一致する。基底状態は位置と運動量の不確定性の積が下限 ℏ/2\hbar/2 に等しい状態である。

励起状態

ξ=q/ℓ\xi=q/\ell と書くと a^†=12(ξ−ddξ)\hat a^\dagger=\frac{1}{\sqrt2}\bigl(\xi-\frac{\mathrm{d} }{\mathrm{d} \xi}\bigr) で、式 (4.4.17) の (a^†)n(\hat a^\dagger)^n を ψ0\psi_0 に掛けると

ψn(q)=1π1/42nn! ℓ Hn(ξ) e−ξ2/2,ξ=qℓ(4.4.22)\psi_n(q)=\frac{1}{\pi^{1/4}\sqrt{2^nn!\,\ell}}\,H_n(\xi)\,\mathrm{e}^{-\xi^2/2},\qquad \xi=\frac q\ell\tag{4.4.22}

を得る。HnH_n は nn 次の Hermite 多項式で、H0=1H_0=1, H1=2ξH_1=2\xi, H2=4ξ2−2H_2=4\xi^2-2, H3=8ξ3−12ξH_3=8\xi^3-12\xi である。a^†\hat a^\dagger を1回掛けるごとに Hn+1(ξ)=2ξHn(ξ)−Hn′(ξ)H_{n+1}(\xi)=2\xi H_n(\xi)-H_n'(\xi) の漸化式が1回使われ、次数が1つ上がり、節(ψn=0\psi_n=0 の点)が1つ増える。Hermite 多項式は、授業では級数の打ち切りとして出てくる(4.4.6節)が、この章の道筋では生成演算子を繰り返し掛けた結果である。

古典の分布との比較

古典の振動子をエネルギー EE で振らせ、ランダムな時刻に位置を測ると、位置の分布は粒子がそこに滞在する時間に比例する。速さ ∣q˙∣=ωA2−q2\lvert \dot q\rvert =\omega\sqrt{A^2-q^2}(AA は振幅)に反比例するので

Pcl(q)=1πA2−q2(∣q∣<A),A=2Jmω(4.4.23)P_{\mathrm{cl}}(q)=\frac{1}{\pi\sqrt{A^2-q^2}}\quad(\lvert q\rvert <A),\qquad A=\sqrt{\frac{2J}{m\omega}}\tag{4.4.23}

で、転回点 q=±Aq=\pm A で発散する。角変数で言えば、θ\theta は一様に進むので θ\theta について一様な分布を q=Asin⁡θq=A\sin\theta で写したものである。

量子の ∣ψn∣2\lvert \psi_n\rvert ^2 と比べるには、作用変数をそろえて J=Jn=(n+12)ℏJ=J_n=(n+\frac12)\hbar、つまり A=ℓ2n+1A=\ell\sqrt{2n+1} にとる。図 4.4.2 に n=10n=10 の場合を描いた。∣ψ10∣2\lvert \psi_{10}\rvert ^2 は 1111 個の山を持って激しく振動するが、山の高さの包絡線は古典の分布をなぞり、転回点の近くで高く、中央で低い。古典の転回点の外(古典では入れない領域)にも確率がしみ出しており、その確率は n=0n=0 で erfc(1)=0.157\mathrm{erfc}(1)=0.157、n=10n=10 で 0.0600.060 である。

図

図 4.4.2 n=10n=10 の確率密度 ∣ψ10∣2\lvert \psi_{10}\rvert ^2(点線)と、同じ作用変数 J=10.5 ℏJ=10.5\,\hbar の古典の振動子の位置の分布(実線)。破線は古典の転回点 ξ=±21\xi=\pm\sqrt{21}。量子の分布は細かく振動するが、振動を均すと古典の分布に沿い、転回点の外へわずかにしみ出す

nn を大きくすると、山の間隔は ℓ/n\ell/\sqrt n 程度に細かくなり、決まった幅で平均した分布は古典の分布に近づく。ξ\xi の幅 0.1 ξ00.1\,\xi_0(ξ0=2n+1\xi_0=\sqrt{2n+1} は転回点)で均した分布と古典の分布の相対差の最大値(転回点の近くを除いた ∣ξ∣<0.8 ξ0\lvert \xi\rvert <0.8\,\xi_0 の範囲)は、n=10n=10 で 0.710.71、n=40n=40 で 0.210.21 と小さくなる。平均をとらない ∣ψn∣2\lvert \psi_n\rvert ^2 そのものは、どれだけ nn を大きくしても古典の分布に近づかない(節の点では 00 のまま)。古典の分布に戻るのは、粗く見た量である。

近似の裏付け

型 B(極限) | 小さいもの ℏ/J→0\hbar/J\to0(1/n→01/n\to0。量子数が大きい) | 証明 要確認(決まった幅で平均した ∣ψn∣2\lvert \psi_n\rvert ^2 が古典の分布に近づくことは Hermite 関数の漸近形から出るとされるが、原論文で前提を確かめていない。平均をとらない密度は収束しない) | 行き先 近似の数理 第7章

例題 4.4.4 位置の2乗平均と4乗平均を古典と比べる ★★

オリジナル

  1. (1) 式 (4.4.12) を使って ⟨n∣q^2∣n⟩\langle n | \hat q^2 | n \rangle を求めよ。
  2. (2) 古典の振動子の q2q^2 の時間平均を作用変数 JJ で表し、J=JnJ=J_n とおいて (1) と比べよ。
  3. (3) 同じことを q4q^4 で行う。⟨n∣q^4∣n⟩=34ℓ4(2n2+2n+1)\langle n | \hat q^4 | n \rangle =\frac34\ell^4(2n^2+2n+1) を用いて、n=3n=3 のときの量子の値と、J=J3J=J_3 の古典の値を比べよ。
解答

(1) q^2=ℓ22(a^+a^†)2=ℓ22(a^2+a^†2+a^a^†+a^†a^)\hat q^2=\frac{\ell^2}{2}(\hat a+\hat a^\dagger)^2=\frac{\ell^2}{2}(\hat a^2+\hat a^{\dagger2}+\hat a\hat a^\dagger+\hat a^\dagger\hat a)。a^2\hat a^2 と a^†2\hat a^{\dagger2} は nn を2つずらすので対角要素を持たず、a^a^†+a^†a^=2n^+1\hat a\hat a^\dagger+\hat a^\dagger\hat a=2\hat n+1 なので⟨n∣q^2∣n⟩=ℓ22(2n+1)=ℏmω(n+12)\langle n | \hat q^2 | n \rangle =\frac{\ell^2}{2}(2n+1)=\frac{\hbar}{m\omega}\bigl(n+\frac12\bigr)。

(2) q=2J/(mω)sin⁡θq=\sqrt{2J/(m\omega)}\sin\theta で θ\theta は一様に進むので、sin⁡2θ‾=12\overline{\sin^2\theta}=\frac12 から q2‾=J/(mω)\overline{q^2}=J/(m\omega)。J=Jn=(n+12)ℏJ=J_n=(n+\frac12)\hbar を入れると ℏmω(n+12)\frac{\hbar}{m\omega}(n+\frac12) で、(1) と完全に一致する。2乗平均はエネルギーと同じく作用変数だけで決まり、作用変数を JnJ_n にとれば量子と古典の差は無い。

(3) 量子は 34ℓ4(2⋅9+6+1)=754ℓ4=18.75 ℓ4\frac34\ell^4(2\cdot9+6+1)=\frac{75}{4}\ell^4=18.75\,\ell^4。古典は sin⁡4θ‾=38\overline{\sin^4\theta}=\frac38 からq4‾=4J2m2ω2⋅38=32J2m2ω2\overline{q^4}=\frac{4J^2}{m^2\omega^2}\cdot\frac38=\frac32\frac{J^2}{m^2\omega^2}、J=72ℏJ=\frac72\hbar で 32⋅494ℓ4=1478ℓ4=18.375 ℓ4\frac32\cdot\frac{49}{4}\ell^4=\frac{147}{8}\ell^4=18.375\,\ell^4。差は 38ℓ4\frac38\ell^4 で、nn によらない定数である。古典の値 34ℓ4(2n2+2n+12)\frac34\ell^4(2n^2+2n+\frac12) と比べると、量子の値は 34ℓ4⋅12\frac34\ell^4\cdot\frac12 だけ大きい。nn が大きいと主要項 32ℓ4n2\frac32\ell^4n^2 に比べて差は相対的に 1/n21/n^2 で小さくなるが、2乗平均のような厳密な一致は4乗では起きない。

ほかの巻では

木下の勉強ステータス

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

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

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