tecktom.com

逆向きの問題と、短い時間の核

2026年10月17日12分で読めますAI生成教科書『虎の巻』第3部 第5章 経路積分から Schrödinger 方程式を導く
目次
  1. 手元にあるもの
  2. 核の幅は
  3. 平面波では展開がいらない
  4. 木下の層

AI 生成・木下が未勉強

量子力学の様々な構成方法 第3部 第5章 §3.5.1 ★★★ 経路積分から微分方程式を出すには、1区間の核が「小さい幅の Gauss 関数」であることだけを使えばよい

第3部第2章(第3部第2章)は、基準(第1部の状態ベクトルと演算子)から出発し、時間分割と完全系の挿入と Trotter の積公式によって、伝播関数 ⟨xf∣e−iH^t/ℏ∣xi⟩\langle x_f | \mathrm{e}^{-\mathrm{i}\hat Ht/\hbar} | x_i \rangle が時間分割した経路積分に等しいことを示した。同値性の地図の矢印「状態ベクトル → 経路積分」はそれで証明済みである。この章は第3部の型(出発点の要請 → 基本の対象 → 基準との同値 → 共通例題 → 得意な問題 → 限界)のうち「基準との同値」の後半、逆向きの矢印「経路積分 → 状態ベクトル」を受け持つ。

逆向きの議論では、出発点に置くものを入れ替える。手元にあるのは第3部第1章(第3部第1章)の Feynman の要請だけで、「各経路に振幅 eiS[x]/ℏ\mathrm{e}^{\mathrm{i} S[x]/\hbar} を割り当てて足すと伝播関数 KK になり、波動関数は ψ(x,t)=∫dx′ K(x,t;x′,t′)ψ(x′,t′)\psi(x,t)=\int\mathrm{d} x'\,K(x,t;x',t')\psi(x',t') で伝わり、確率密度は ∣ψ∣2\lvert \psi\rvert ^2 である」という内容である。ヒルベルト空間も演算子もハミルトニアンも、まだ持っていない。この章で示すのは、これだけの材料から次の4つが出てくることである。

  1. 波動関数は Schrödinger 方程式 iℏ∂tψ=[−ℏ22m∇2+V]ψ\mathrm{i}\hbar\partial_t\psi=\bigl[-\frac{\hbar^2}{2m}\nabla^2+V\bigr]\psi を満たす(3.5.2節)。ベクトルポテンシャルがあれば (p^−qA)2/(2m)(\hat{\bm{p}}-q\bm{A})^2/(2m) の形になる(3.5.3節)。
  2. 時間発展はユニタリで、時間について群をなす。Stone の定理によりその生成子が自己共役なハミルトニアン H^\hat H で、Chernoff の積公式により時間分割の極限は e−iH^t/ℏ\mathrm{e}^{-\mathrm{i}\hat Ht/\hbar} そのものである(3.5.4節)。
  3. 位置と運動量の交換関係 [x^, p^]=iℏ\left[\hat x,\,\hat p\right]=\mathrm{i}\hbar と Heisenberg の運動方程式が、経路積分の部分積分の恒等式として出てくる(3.5.5節)。
  4. 1と2と3をあわせると基準の公理 Q5(ユニタリな時間発展)と正準交換関係が復元され、第3部第2章の命題と組んで両向きの同値が完成する(3.5.6節)。

議論の中心は、1区間の振幅 KϵK_\epsilon が「幅 ℏϵ/m\sqrt{\hbar\epsilon/m} の振動する Gauss 関数」であるという1つの事実である。第3部第2章の例題で求めた変位のモーメント⟨η2⟩=iℏϵ/m\langle \eta^2 \rangle =\mathrm{i}\hbar\epsilon/m が、2階微分(運動エネルギー)・評価点の問題・交換関係の3つをまとめて決める。記号は第3部第1・2章に揃え、1次元の位置を xx、区間の幅を ϵ\epsilon、1区間の変位を η=xk+1−xk\eta=x_{k+1}-x_k と書く。ポテンシャル VV は区間の前の点 xkx_k で評価するのを既定とする。

手元にあるもの

Feynman の要請(第3部第1章)では、伝播関数は時間分割した積分の極限

K(xf,t;xi,0)=lim⁡N→∞∫dx1⋯∫dxN−1∏k=0N−1Kϵ(xk+1,xk),ϵ=tN,x0=xi, xN=xf(3.5.1)K(x_f,t;x_i,0)=\lim_{N\to\infty}\int\mathrm{d} x_1\cdots\int\mathrm{d} x_{N-1}\prod_{k=0}^{N-1}K_\epsilon(x_{k+1},x_k),\qquad \epsilon=\frac tN,\quad x_0=x_i,\ x_N=x_f\tag{3.5.1}

で、1区間の振幅は

Kϵ(x′,x)=m2πiℏϵ  exp⁡{iℏ[m(x′−x)22ϵ−ϵV(x)]}(3.5.2)K_\epsilon(x',x)=\sqrt{\frac{m}{2\pi\mathrm{i}\hbar\epsilon}}\;\exp\Bigl\{\frac{\mathrm{i}}{\hbar}\Bigl[\frac{m(x'-x)^2}{2\epsilon}-\epsilon V(x)\Bigr]\Bigr\}\tag{3.5.2}

である。指数の中身は、区間 [tk,tk+1][t_k,t_{k+1}] を等速で進む直線の経路の作用 ∫(m2x˙2−V)dt\int(\frac m2\dot x^2-V)\mathrm{d} t を、VV だけ前の点で評価したものである。前因子 m/(2πiℏϵ)\sqrt{m/(2\pi\mathrm{i}\hbar\epsilon)} は第3部第1章では「ϵ→0\epsilon\to0 で KϵK_\epsilon がデルタ関数に近づく」ように選んだ定数で、この章では3.5.2節でもう一度その条件から決め直す。

波動関数は、この伝播関数で過去の波動関数から作られる量として定義される。式 (3.5.1) の積は区間ごとに分かれているので、最後の1区間だけを切り離すと

ψ(x,t+ϵ)=∫dy  Kϵ(x,y) ψ(y,t)(3.5.3)\psi(x,t+\epsilon)=\int\mathrm{d} y\;K_\epsilon(x,y)\,\psi(y,t)\tag{3.5.3}

が成り立つ(極限 N→∞N\to\infty の前の段階で、tt と t+ϵt+\epsilon がどちらも区切りの時刻であるとする)。この式は経路積分の合成則(第3部第1章)の、片方の区間を1区間にした特別な場合である。逆向きの問題は、式 (3.5.3) を ϵ\epsilon について展開し、ψ(x,t+ϵ)−ψ(x,t)\psi(x,t+\epsilon)-\psi(x,t) の ϵ\epsilon の1次の係数を ψ(x,t)\psi(x,t) とその微分で書くことに帰着する。それが書ければ、ϵ→0\epsilon\to0 で ∂tψ\partial_t\psi を与える微分方程式が得られる。

核の幅は ϵ\sqrt\epsilon

y=x−ηy=x-\eta と変数を替えると、核は η\eta の関数として m/(2πiℏϵ) eimη2/(2ℏϵ)\sqrt{m/(2\pi\mathrm{i}\hbar\epsilon)}\,\mathrm{e}^{\mathrm{i} m\eta^2/(2\hbar\epsilon)} に e−iϵV/ℏ\mathrm{e}^{-\mathrm{i}\epsilon V/\hbar} を掛けたものになる。前者の位相 mη2/(2ℏϵ)m\eta^2/(2\hbar\epsilon) は ∣η∣≲ℏϵ/m\lvert \eta\rvert \lesssim\sqrt{\hbar\epsilon/m} ではほとんど変わらず、それより外では急に振動する(図 3.5.1)。ψ(y)\psi(y) がこの幅より長い尺度でしか変化しなければ、外側の振動する部分は ψ\psi を掛けて積分すると打ち消し合い、積分に残るのは幅 ℏϵ/m\sqrt{\hbar\epsilon/m} の内側の寄与だけになる。内側では ψ(x−η)\psi(x-\eta) を η\eta について Taylor 展開でき、各項は核のモーメント

⟨ηn⟩≔∫dη m2πiℏϵ eimη2/(2ℏϵ) ηn:⟨1⟩=1,⟨η⟩=⟨η3⟩=0,⟨η2⟩=iℏϵm,⟨η4⟩=3(iℏϵm)2(3.5.4)\langle \eta^n \rangle \coloneqq\int\mathrm{d}\eta\,\sqrt{\frac{m}{2\pi\mathrm{i}\hbar\epsilon}}\,\mathrm{e}^{\mathrm{i} m\eta^2/(2\hbar\epsilon)}\,\eta^n:\qquad \langle 1 \rangle =1,\quad \langle \eta \rangle =\langle \eta^3 \rangle =0,\quad \langle \eta^2 \rangle =\frac{\mathrm{i}\hbar\epsilon}{m},\quad \langle \eta^4 \rangle =3\Bigl(\frac{\mathrm{i}\hbar\epsilon}{m}\Bigr)^2\tag{3.5.4}

で平均される(第3部第2章の例題で求めたもの。振動積分は時間に ϵ(1−iδ)\epsilon(1-\mathrm{i}\delta) の収束因子を入れて計算する)。η2j\eta^{2j} のモーメントは ϵj\epsilon^j の大きさなので、ϵ\epsilon の1次までを求めるには ψ\psi を η\eta の2次まで展開すれば足りる。これがこの章のすべての計算の骨組みである。

図

図 3.5.1 1区間の核の実部 cos⁡(u2/2)\cos(u^2/2)(細線)と、ゆっくり変わる波動関数(太線)。横軸 u=η/ℏϵ/mu=\eta/\sqrt{\hbar\epsilon/m} は変位 η\eta を核の幅で割ったもの。破線の内側(∣u∣<1\lvert u\rvert <1)では位相がほとんど変わらず、外側では振動が速くなって ψ\psi との積の積分が打ち消し合う。ϵ→0\epsilon\to0 で幅が縮むので、ψ\psi はいずれ内側で Taylor 展開できる。

平面波では展開がいらない

自由粒子(V=0V=0)なら、核は η\eta だけの関数なので式 (3.5.3) は畳み込みで、Fourier 変換すると掛け算になる。平面波 eikx\mathrm{e}^{\mathrm{i} kx} を入れると、Gauss 積分が厳密に実行できて

∫dη m2πiℏϵ eimη2/(2ℏϵ) eik(x−η)=e−iℏk2ϵ/(2m) eikx(3.5.5)\int\mathrm{d}\eta\,\sqrt{\frac{m}{2\pi\mathrm{i}\hbar\epsilon}}\,\mathrm{e}^{\mathrm{i} m\eta^2/(2\hbar\epsilon)}\,\mathrm{e}^{\mathrm{i} k(x-\eta)}=\mathrm{e}^{-\mathrm{i}\hbar k^2\epsilon/(2m)}\,\mathrm{e}^{\mathrm{i} kx}\tag{3.5.5}

となる。左辺の η\eta の積分は c/π∫dη e−cη2−ikη=e−k2/(4c)\sqrt{c/\pi}\int\mathrm{d}\eta\,\mathrm{e}^{-c\eta^2-\mathrm{i} k\eta}=\mathrm{e}^{-k^2/(4c)} に c=−im/(2ℏϵ)c=-\mathrm{i} m/(2\hbar\epsilon) を入れたもので、−k2/(4c)=−iℏk2ϵ/(2m)-k^2/(4c)=-\mathrm{i}\hbar k^2\epsilon/(2m) である。平面波は形を変えずに位相 e−iEkϵ/ℏ\mathrm{e}^{-\mathrm{i} E_k\epsilon/\hbar}(Ek=ℏ2k2/(2m)E_k=\hbar^2k^2/(2m))だけを受け取る。すなわち、経路積分の自由な1区間は Fourier 空間で e−iℏk2ϵ/(2m)\mathrm{e}^{-\mathrm{i}\hbar k^2\epsilon/(2m)} を掛ける操作で、これは基準の自由な時間発展 e−ip^2ϵ/(2mℏ)\mathrm{e}^{-\mathrm{i}\hat p^2\epsilon/(2m\hbar)} の運動量表示に他ならない。展開を使う次節の計算は、この厳密な式の ϵ\epsilon の1次をとることと同じである。

例題 3.5.1 平面波が受け取る位相と展開の条件 ★★

オリジナル自由粒子の1区間の核を平面波 eikx\mathrm{e}^{\mathrm{i} kx} に当てる。

  1. (1) 式 (3.5.5) の右辺を ϵ\epsilon の1次まで展開し、iℏ (ψ(t+ϵ)−ψ(t))/ϵ→−ℏ22m∂x2ψ\mathrm{i}\hbar\,(\psi(t+\epsilon)-\psi(t))/\epsilon\to-\frac{\hbar^2}{2m}\partial_x^2\psi となることを確かめよ。
  2. (2) 一般の波束 ψ(x)=∫dk2πψ~(k)eikx\psi(x)=\int\frac{\mathrm{d} k}{2\pi}\tilde\psi(k)\mathrm{e}^{\mathrm{i} kx} について、1次で切った誤差の L2L^2 ノルムが ϵ22(ℏ2m)2∥∂x4ψ∥\frac{\epsilon^2}{2}\bigl(\frac{\hbar}{2m}\bigr)^2\lVert \partial_x^4\psi\rVert 以下であることを示せ。
  3. (3) 電子の運動エネルギーが 1 eV1\,\mathrm{eV} のとき、展開が使える時間の刻みの目安 ℏ/Ek\hbar/E_k を求めよ。
解答

(1) e−iℏk2ϵ/(2m)=1−iℏk2ϵ/(2m)+O(ϵ2)\mathrm{e}^{-\mathrm{i}\hbar k^2\epsilon/(2m)}=1-\mathrm{i}\hbar k^2\epsilon/(2m)+\mathcal{O}(\epsilon^2) なので ψ(x,t+ϵ)−ψ(x,t)=−iℏk2ϵ2meikx+O(ϵ2)\psi(x,t+\epsilon)-\psi(x,t)=-\mathrm{i}\frac{\hbar k^2\epsilon}{2m}\mathrm{e}^{\mathrm{i} kx}+\mathcal{O}(\epsilon^2) である。∂x2eikx=−k2eikx\partial_x^2\mathrm{e}^{\mathrm{i} kx}=-k^2\mathrm{e}^{\mathrm{i} kx} を使うと iℏ (ψ(t+ϵ)−ψ(t))/ϵ=ℏ2k22meikx+O(ϵ)=−ℏ22m∂x2eikx+O(ϵ)\mathrm{i}\hbar\,(\psi(t+\epsilon)-\psi(t))/\epsilon=\frac{\hbar^2k^2}{2m}\mathrm{e}^{\mathrm{i} kx}+\mathcal{O}(\epsilon)=-\frac{\hbar^2}{2m}\partial_x^2\mathrm{e}^{\mathrm{i} kx}+\mathcal{O}(\epsilon) となる。

(2) 実数 θ\theta について ∣e−iθ−1+iθ∣=∣∫0θ(e−is−1)(−i) ds∣≤∫0∣θ∣s ds=θ2/2\lvert \mathrm{e}^{-\mathrm{i}\theta}-1+\mathrm{i}\theta\rvert =\bigl\lvert\int_0^\theta(\mathrm{e}^{-\mathrm{i} s}-1)(-\mathrm{i})\,\mathrm{d} s\bigr\rvert\le\int_0^{\lvert \theta\rvert }s\,\mathrm{d} s=\theta^2/2 が成り立つ(∣e−is−1∣≤∣s∣\lvert \mathrm{e}^{-\mathrm{i} s}-1\rvert \le\lvert s\rvert による)。θ=ℏk2ϵ/(2m)\theta=\hbar k^2\epsilon/(2m) として各 kk の成分に当て、Plancherel の等式で xx の側に戻すと、誤差のノルムは 12(ℏϵ2m)2∥k4ψ~∥=ϵ22(ℏ2m)2∥∂x4ψ∥\frac12\bigl(\frac{\hbar\epsilon}{2m}\bigr)^2\lVert k^4\tilde\psi\rVert =\frac{\epsilon^2}{2}\bigl(\frac{\hbar}{2m}\bigr)^2\lVert \partial_x^4\psi\rVert 以下である。展開は ψ\psi が4回微分できて ∂x4ψ\partial_x^4\psi が2乗可積分であれば、各点でなく L2L^2 の意味で正当化される。

(3) ℏ/Ek=1.055×10−34/(1.602×10−19) s=6.58×10−16 s\hbar/E_k=1.055\times10^{-34}/(1.602\times10^{-19})\,\mathrm s=6.58\times10^{-16}\,\mathrm s である。ϵ\epsilon がこれより十分短ければ位相 Ekϵ/ℏE_k\epsilon/\hbar は小さく、1次の展開が使える。運動エネルギーの高い成分ほど小さい ϵ\epsilon を要するが、ϵ→0\epsilon\to0 の極限ではどの成分についても条件が満たされる。

木下の勉強ステータス

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

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

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