tecktom.com

この章のまとめ

2026年10月17日12分で読めますAI生成教科書『虎の巻』第3部 第3章 自由粒子と調和振動子を Gauss 積分で解く

AI 生成・木下が未勉強

量子力学の様々な構成方法 第3部 第3章 §3.3.7 ★★ 2次の作用では「古典の作用+揺らぎの行列式」で経路積分が厳密に閉じ、共通例題①③が出そろう

作用が経路の2次式なら、経路を古典の経路と両端が 00 の揺らぎに分けたとき、作用は Scl+S2[y]S_{\mathrm{cl}}+S_2[y] と厳密に分かれ、伝播関数は K=F(t) eiScl/ℏK=F(t)\,\mathrm{e}^{\mathrm{i} S_{\mathrm{cl}}/\hbar} になる。端の位置への依存は古典の作用だけが担い、量子力学の計算は端によらない揺らぎの因子 F(t)F(t) に集まる。

自由粒子では揺らぎの積分が3重対角行列の行列式 NN になり、刻みによらず厳密に共通例題③の伝播関数 K0K_0 を得た。古典の作用の勾配が運動量であることから、時刻 tt の位置は初めの x^+p^t/m\hat x+\hat pt/m と同じ分布になり、波束の広がりは古典の経路の束の広がりと運動量の量子の広がりの和として読むことができる。遠方では位置の分布が運動量の分布の写しになる。

調和振動子では、古典の作用が境界項 m2[xx˙]0t\frac m2[x\dot x]_0^t で書け、揺らぎの因子は Gel'fand–Yaglom の初期値問題の解 ϕ(t)=sin⁡ωt/ω\phi(t)=\sin\omega t/\omega で決まり、共通例題①の Mehler 核を得た。橋の節では、ϕ(t)\phi(t) が Heisenberg 描像の交換子 (im/ℏ)[x^H(t), x^](\mathrm{i} m/\hbar)\left[\hat x_{\mathrm H}(t),\,\hat x\right] に等しいこと、伝播関数の跡からスペクトル ℏω(n+12)\hbar\omega(n+\frac12) と分配関数が出ることを示した。

自由粒子(共通例題③)調和振動子(共通例題①)
古典の作用 SclS_{\mathrm{cl}}m(xf−xi)2/(2t)m(x_f-x_i)^2/(2t)式 (3.3.21)
揺らぎの関数 ϕ(t)=m b(t)\phi(t)=m\,b(t)ttsin⁡ωt/ω\sin\omega t/\omega
揺らぎの因子 F(t)F(t)m/(2πiℏt)\sqrt{m/(2\pi\mathrm{i}\hbar t)}mω/(2πiℏsin⁡ωt)\sqrt{m\omega/(2\pi\mathrm{i}\hbar\sin\omega t)}
刻み NN のずれなし(厳密)O(1/N2)O(1/N^2)

2次でない作用では展開が2次で終わらず、この章の分け方は ħ が小さいときの近似(停留位相)になる。それが第3部第4章の主題である。

総合問題 一定の力を受ける調和振動子 ★★★

オリジナル質量 mm・角振動数 ω\omega の調和振動子に一定の力 ff がかかる系 L=m2x˙2−mω22x2+f xL=\frac m2\dot x^2-\frac{m\omega^2}2x^2+f\,x を考える。0<ωt<π0<\omega t<\pi とする。

  1. (1) 伝播関数が Kf=F(t) eiScl(f)/ℏK_f=F(t)\,\mathrm{e}^{\mathrm{i} S_{\mathrm{cl}}^{(f)}/\hbar} と書け、揺らぎの因子 F(t)F(t) が f=0f=0 の調和振動子と同じであることを示せ。(使う道具:3.3.1節・3.3.5節)
  2. (2) x0=f/(mω2)x_0=f/(m\omega^2) だけずらした座標 u=x−x0u=x-x_0 を使い、古典の作用が

Scl(f)=mω2sin⁡ωt[(uf2+ui2)cos⁡ωt−2ufui]+f2t2mω2,uf,i=xf,i−x0S_{\mathrm{cl}}^{(f)}=\frac{m\omega}{2\sin\omega t}\Bigl[(u_f^2+u_i^2)\cos\omega t-2u_fu_i\Bigr]+\frac{f^2t}{2m\omega^2},\qquad u_{f,i}=x_{f,i}-x_0

となることを示せ。(使う道具:3.3.4節)3. (3) 伝播関数の跡をとり、エネルギー準位を求めよ。基準の言葉で Hamiltonian の平方完成から出した値と比べよ。(使う道具:3.3.6節・第1部第2章の Q5)4. (4) ω→0\omega\to0 の極限で、一様な力の古典の作用 m(xf−xi)22t+ft(xf+xi)2−f2t324m\frac{m(x_f-x_i)^2}{2t}+\frac{ft(x_f+x_i)}2-\frac{f^2t^3}{24m} を得ることを示せ。(2) の第2項 f2t/(2mω2)f^2t/(2m\omega^2) は ω→0\omega\to0 で発散するが、どこで打ち消されるか。(使う道具:3.3.2節・3.3.4節・例題 3.3.4)5. (5) Heisenberg 描像の x^H(t)\hat x_{\mathrm H}(t) を求め、式 (3.3.33) の b(t)b(t) が ff によらないことから (1) を基準の側から確かめよ。(使う道具:3.3.6節・第2部第3章の Heisenberg 方程式)6. (6) Paul トラップの 40Ca+^{40}\mathrm{Ca}^+ イオン(m=6.636×10−26 kgm=6.636\times10^{-26}\,\mathrm{kg}、ω/2π=1 MHz\omega/2\pi=1\,\mathrm{MHz})に 10 mV m−110\,\mathrm{mV\,m^{-1}} の電場をかけた。x0x_0・準位の下がり f2/(2mω2)f^2/(2m\omega^2) を ℏω\hbar\omega 単位で・基底状態の広がり ℏ/(2mω)\sqrt{\hbar/(2m\omega)} を求めよ。(使う道具:(2)(3))

解答

(1) Lagrangian は xx と x˙\dot x の2次式なので、3.3.1節の議論がそのまま当てはまる。古典の経路 xclx_{\mathrm{cl}} からのずれ yy を入れると、1次の項は Euler–Lagrange 方程式 mx¨cl=−mω2xcl+fm\ddot x_{\mathrm{cl}}=-m\omega^2x_{\mathrm{cl}}+f で消え、2次の項は式 (3.3.4) の S2[y]S_2[y] で ff を含まない。ff は xx について1次なので、2次の揺らぎの作用に入りようがない。したがって F(t)F(t) は3.3.5節の mω/(2πiℏsin⁡ωt)\sqrt{m\omega/(2\pi\mathrm{i}\hbar\sin\omega t)} そのものである。

(2) Lagrangian を平方完成すると

−mω22x2+fx=−mω22(x−x0)2+f22mω2-\frac{m\omega^2}2x^2+fx=-\frac{m\omega^2}2(x-x_0)^2+\frac{f^2}{2m\omega^2}

である。u˙=x˙\dot u=\dot x なので L=m2u˙2−mω22u2+f22mω2L=\frac m2\dot u^2-\frac{m\omega^2}2u^2+\frac{f^2}{2m\omega^2} となり、最後の定数は作用に f2t/(2mω2)f^2t/(2m\omega^2) を足すだけである。uu についての作用は式 (3.3.21) の xx を uu に置き換えたものなので、示すべき式を得る。古典の経路は xcl=x0+uclx_{\mathrm{cl}}=x_0+u_{\mathrm{cl}} で、つり合いの位置 x0x_0 のまわりの振動である。

(3) 跡では xf=xi=xx_f=x_i=x とし xx で積分する。u=x−x0u=x-x_0 に変数を変えても積分範囲は変わらないので、f=0f=0 の跡 (3.3.34) に定数の位相が掛かるだけで

∫dx Kf(x,t;x,0)=eif2t/(2mω2ℏ)∑n=0∞e−iω(n+12)t=∑n=0∞exp⁡[−itℏ(ℏω(n+12)−f22mω2)]\int\mathrm{d} x\,K_f(x,t;x,0)=\mathrm{e}^{\mathrm{i} f^2t/(2m\omega^2\hbar)}\sum_{n=0}^{\infty}\mathrm{e}^{-\mathrm{i}\omega(n+\frac12)t} =\sum_{n=0}^{\infty}\exp\Bigl[-\frac{\mathrm{i} t}\hbar\Bigl(\hbar\omega\bigl(n+\tfrac12\bigr)-\frac{f^2}{2m\omega^2}\Bigr)\Bigr]

となる。したがって En=ℏω(n+12)−f2/(2mω2)E_n=\hbar\omega(n+\frac12)-f^2/(2m\omega^2) である。基準の側では H^=p^22m+mω22x^2−fx^=p^22m+mω22(x^−x0)2−f22mω2\hat H=\frac{\hat p^2}{2m}+\frac{m\omega^2}2\hat x^2-f\hat x=\frac{\hat p^2}{2m}+\frac{m\omega^2}2(\hat x-x_0)^2-\frac{f^2}{2m\omega^2} で、x^−x0\hat x-x_0 と p^\hat p は x^\hat x と p^\hat p と同じ交換関係を満たす(並進のユニタリ変換で移り合う)ので、スペクトルは調和振動子のものを f2/(2mω2)f^2/(2m\omega^2) だけ下げたものになる。両者は一致する。

(4) (2) の式を ω\omega について展開する。第1項を uf,i=xf,i−x0u_{f,i}=x_{f,i}-x_0 のまま展開すると、x0=f/(mω2)x_0=f/(m\omega^2) が ω−2\omega^{-2} に比例するので、uf2+ui2u_f^2+u_i^2 と ufuiu_fu_i の項から x02x_0^2 に比例する ω−4\omega^{-4}・ω−2\omega^{-2} の項が出る。u=x−x0u=x-x_0 を代入して x0x_0 の次数で整理すると

mω2sin⁡ωt[(uf2+ui2)cos⁡ωt−2ufui]=mω2sin⁡ωt[(xf2+xi2)cos⁡ωt−2xfxi]+mωx0(xf+xi)(1−cos⁡ωt)sin⁡ωt−mωx02(1−cos⁡ωt)sin⁡ωt\begin{aligned} \frac{m\omega}{2\sin\omega t}\Bigl[(u_f^2+u_i^2)\cos\omega t-2u_fu_i\Bigr] &=\frac{m\omega}{2\sin\omega t}\Bigl[(x_f^2+x_i^2)\cos\omega t-2x_fx_i\Bigr]\\ &\quad+\frac{m\omega x_0(x_f+x_i)(1-\cos\omega t)}{\sin\omega t}-\frac{m\omega x_0^2(1-\cos\omega t)}{\sin\omega t} \end{aligned}

となる(x0x_0 の1次と2次の項を集めた)。(1−cos⁡ωt)/sin⁡ωt=tan⁡(ωt/2)=ωt2+(ωt)324+⋯(1-\cos\omega t)/\sin\omega t=\tan(\omega t/2)=\frac{\omega t}2+\frac{(\omega t)^3}{24}+\cdots を使うと、第2項は mωx0(xf+xi)ωt2+O(ω2)=ft(xf+xi)2+O(ω2)\frac{m\omega x_0(x_f+x_i)\omega t}{2}+O(\omega^2)=\frac{ft(x_f+x_i)}2+O(\omega^2)、第3項は

−mωx02(ωt2+(ωt)324)+⋯=−f2t2mω2−f2t324m+O(ω2)-m\omega x_0^2\Bigl(\frac{\omega t}2+\frac{(\omega t)^3}{24}\Bigr)+\cdots=-\frac{f^2t}{2m\omega^2}-\frac{f^2t^3}{24m}+O(\omega^2)

である。第3項の主要部 −f2t/(2mω2)-f^2t/(2m\omega^2) が (2) の定数 +f2t/(2mω2)+f^2t/(2m\omega^2) とちょうど打ち消し合い、第1項は例題 3.3.4(1) で m(xf−xi)2/(2t)m(x_f-x_i)^2/(2t) に近づく。合わせて

Scl(f)→m(xf−xi)22t+ft(xf+xi)2−f2t324mS_{\mathrm{cl}}^{(f)}\to\frac{m(x_f-x_i)^2}{2t}+\frac{ft(x_f+x_i)}2-\frac{f^2t^3}{24m}

を得る。発散はつり合いの位置 x0x_0 が無限遠へ逃げることによる見かけのもので、x02x_0^2 に比例する項と定数項を足し合わせると消える。揺らぎの因子は ω→0\omega\to0 で m/(2πiℏt)\sqrt{m/(2\pi\mathrm{i}\hbar t)} になるので、一様な力の伝播関数は m/(2πiℏt) eiScl(f)/ℏ\sqrt{m/(2\pi\mathrm{i}\hbar t)}\,\mathrm{e}^{\mathrm{i} S_{\mathrm{cl}}^{(f)}/\hbar} である(検算でこれが Schrödinger 方程式を満たすことを確かめた)。

(5) Heisenberg 方程式 x^˙H=p^H/m\dot{\hat x}_{\mathrm H}=\hat p_{\mathrm H}/m、p^˙H=−mω2x^H+f\dot{\hat p}_{\mathrm H}=-m\omega^2\hat x_{\mathrm H}+f は古典の運動方程式と同じ形で、初期条件 x^H(0)=x^\hat x_{\mathrm H}(0)=\hat x、p^H(0)=p^\hat p_{\mathrm H}(0)=\hat p の解は

x^H(t)=x^cos⁡ωt+p^mωsin⁡ωt+x0(1−cos⁡ωt)\hat x_{\mathrm H}(t)=\hat x\cos\omega t+\frac{\hat p}{m\omega}\sin\omega t+x_0(1-\cos\omega t)

である。ff は定数項 x0(1−cos⁡ωt)x_0(1-\cos\omega t) にしか入らず、交換子 [x^H(t), x^]=−iℏsin⁡ωt/(mω)\left[\hat x_{\mathrm H}(t),\,\hat x\right]=-\mathrm{i}\hbar\sin\omega t/(m\omega) は f=0f=0 と同じである。式 (3.3.33) から C(t)=F(t)C(t)=F(t) も ff によらない。定数項は式 (3.3.30) の左辺の xf−axix_f-a x_i から x0(1−cos⁡ωt)x_0(1-\cos\omega t) を引くだけで、指数の1次の項(つまり古典の作用)だけを変える。経路積分の側で「ff は1次の項にしか入らない」と言ったことが、行列力学の側では「ff は定数項にしか入らない」ことに対応する。

(6) 力は f=eE=1.602×10−19×10−2 N=1.602×10−21 Nf=eE=1.602\times10^{-19}\times10^{-2}\,\mathrm N=1.602\times10^{-21}\,\mathrm N、ω=2π×106 s−1\omega=2\pi\times10^6\,\mathrm s^{-1} である。

x0=fmω2=1.602×10−216.636×10−26×(6.283×106)2 m=6.12×10−10 mx_0=\frac{f}{m\omega^2}=\frac{1.602\times10^{-21}}{6.636\times10^{-26}\times(6.283\times10^6)^2}\,\mathrm m=6.12\times10^{-10}\,\mathrm m

で、準位の下がりは f2/(2mω2)=fx0/2=4.90×10−31 Jf^2/(2m\omega^2)=fx_0/2=4.90\times10^{-31}\,\mathrm J、ℏω=6.63×10−28 J\hbar\omega=6.63\times10^{-28}\,\mathrm J で割って 7.39×10−4 ℏω7.39\times10^{-4}\,\hbar\omega である。基底状態の広がりは ℏ/(2mω)=1.12×10−8 m\sqrt{\hbar/(2m\omega)}=1.12\times10^{-8}\,\mathrm m で、ずれ x0=0.61 nmx_0=0.61\,\mathrm{nm} はその約 1/181/18 にとどまる。(1) のとおり、この電場で変わるのは古典の作用(中心の位置と全体の位相)だけで、波束の形を決める揺らぎの因子は変わらない。

木下の勉強ステータス

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

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

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