tecktom.com

案内方程式:位置を状態に加える

2026年10月17日12分で読めますAI生成教科書『虎の巻』第7部 第4章 Bohm の軌道:案内方程式と、構成としての同値と解釈の境目
目次
  1. 4つの要請
  2. 等変性:分布は のまま運ばれる
  3. 2階の形:量子ポテンシャルの力
  4. 木下の層

AI 生成・木下が未勉強

量子力学の様々な構成方法 第7部 第4章 §7.4.1 ★★★ この構成の出発点の要請と、要請どうしを整合させる等変性の定理

第7部第1章では、波動関数を極形式 ψ=ρ eiS/ℏ\psi=\sqrt\rho\,\mathrm{e}^{\mathrm{i} S/\hbar} に書いて、Schrödinger 方程式を密度 ρ\rho と速度場 v=∇S/m\bm{v}=\nabla S/m の流体の方程式に書き直した(第7部第1章)。流体の方程式は場の方程式で、各点で流れの速さを与えるだけである。流れの速さが各点で決まっているなら、その流れに乗って運ばれる点の道筋(流線)も決まる。L. de Broglie(1927)と D. Bohm(1952)は、この道筋を「粒子が実際にたどる軌道」と読み、量子力学を「波動関数と、その波に案内される粒子の位置」の2つで組み立て直した。この章ではこの構成をBohm の力学(de Broglie–Bohm 理論、パイロット波理論とも呼ぶ)と呼ぶ。

部の型(出発点の要請 → 基本の対象 → 基準との同値 → 共通例題 → 得意な問題 → 限界)のうち、この章は第7部の中で「軌道」という対象についての一巡りを受け持つ。出発点に置くのは、波動関数 ψ\psi と粒子の配置 XX の組、ψ\psi の従う Schrödinger 方程式、XX の従う案内方程式 X˙=v(X,t)\dot X=\bm{v}(X,t)、そして初めの配置が ∣ψ∣2\lvert \psi\rvert ^2 に従って散らばっているという量子平衡の仮説の4つである。第1節で要請を置き、配置の分布が時間がたっても ∣ψ∣2\lvert \psi\rvert ^2 のままである(等変性)ことを証明する。第2節で共通例題①③を軌道で解き、第3節で二重スリットの軌道を数値で描く。第4節で測定を軌道の言葉で書き、Born の規則がどの仮定から出るかを確かめる。第5節は多粒子の速度場の非局所性と、相対論にしにくい理由を扱う。第6節の橋では、Bohm の力学が予言としては基準と同値であること、その同値が「位置の分布」までで止まり、軌道そのものは観測で確かめられない追加の構造であることを、式を並べて示す。

4つの要請

質量 m1,…,mNm_1,\dots,m_N の NN 粒子を考え、配置空間の点を x=(r1,…,rN)∈R3Nx=(\bm{r}_1,\dots,\bm{r}_N)\in\mathbb{R}^{3N} と書く。第7部第1章の流体は配置空間の上の場だったので(第7部第1章 の落とし穴の箱)、軌道も配置空間で考える。

Bohm の力学の要請 ★★★

  1. (B1) 系の状態は、波動関数 ψ(x,t)\psi(x,t) と、実際の配置 X(t)=(X1(t),…,XN(t))X(t)=(\bm{X}_1(t),\dots,\bm{X}_N(t)) の組である。
  2. (B2) 波動関数は Schrödinger 方程式 iℏ ∂tψ=H^ψ\mathrm{i}\hbar\,\partial_t\psi=\hat H\psi に従う(第1部第2章の Q5 と同じ)。
  3. (B3) 配置は案内方程式

dXkdt=vk(X,t),vk(x,t)=ℏmkIm⁡∇kψ(x,t)ψ(x,t)=∇kS(x,t)mk(7.4.1)\frac{\mathrm{d} \bm{X}_k}{\mathrm{d} t}=\bm{v}_k(X,t),\qquad \bm{v}_k(x,t)=\frac{\hbar}{m_k}\operatorname{Im}\frac{\nabla_k\psi(x,t)}{\psi(x,t)}=\frac{\nabla_kS(x,t)}{m_k}\tag{7.4.1}

に従う(∇k\nabla_k は rk\bm{r}_k についての勾配)。4. (B4) 量子平衡の仮説:同じ ψ\psi を用意した系の集団で、初めの配置 X(0)X(0) は確率密度 ∣ψ(x,0)∣2\lvert \psi(x,0)\rvert ^2 に従って散らばっている。

(B3) の速度場は、第7部第1章の速度場 v=∇S/m\bm{v}=\nabla S/m を配置空間に広げたものである。1粒子なら確率流密度 j=ρv\bm{j}=\rho\bm{v} を密度で割った v=j/ρ\bm{v}=\bm{j}/\rho で、流体の各点の流れの速さに等しい。つまり Bohm の軌道は、第7部第1章の確率の流体の流線そのものである。ここで新しく加わったのは方程式でなく、「流線のうち1本が実際に起きている」という主張と、(B4) の初めの分布だけである。

速度場は ψ\psi に定数を掛けても変わらない(∇ψ/ψ\nabla\psi/\psi の比で決まる)。全体の位相を変えても、規格化を変えても軌道は同じなので、第1部第2章の Q1(状態は射線)と矛盾しない。ψ\psi の節(ψ=0\psi=0 の点)では速度場が定まらないが、節は配置空間で余次元2(実部と虚部が両方 00)の集合なので、典型的な初めの位置から出た軌道は節に有限の時間でぶつからない(自己共役なハミルトニアンの大きな族でこれを示したのが Berndl, Dürr, Goldstein, Peruzzi, Zanghì(1995、Commun. Math. Phys. 173, 647、arXiv: quant-ph/9503013)である)。

等変性:分布は ∣ψ∣2\lvert \psi\rvert ^2 のまま運ばれる

(B3) は1階の常微分方程式なので、初めの配置を決めれば軌道はただ1本に決まる。(B4) で初めの配置を ∣ψ∣2\lvert \psi\rvert ^2 に従って散らすと、時刻 tt の配置もある分布 P(x,t)P(x,t) に従う。この PP が ∣ψ(x,t)∣2\lvert \psi(x,t)\rvert ^2 と一致するかどうかは自明でない。一致しなければ、ある時刻に Born の規則を満たしていた集団が、別の時刻には満たさなくなる。次の定理が、要請 (B2)(B3)(B4) が互いに整合していることを保証する。

定理 7.4.1(等変性 ★★★)

速度場 v\bm{v} が滑らかで、軌道が有限の時間で節や無限遠に達しないとする。初めの配置が P(x,0)=∣ψ(x,0)∣2P(x,0)=\lvert \psi(x,0)\rvert ^2 に従うなら、すべての時刻でP(x,t)=∣ψ(x,t)∣2P(x,t)=\lvert \psi(x,t)\rvert ^2 である。

証明

第1段:分布 PP が連続の式に従うことを示す。任意の滑らかで台が有界な関数 f(x)f(x) について、期待値の時間微分は軌道の式から

ddt∫f(x)P(x,t) dx=ddt⟨f(X(t))⟩=⟨∇f(X)⋅X˙⟩=∫∇f⋅v P dx=−∫f ∇⋅(Pv) dx\frac{\mathrm{d} }{\mathrm{d} t}\int f(x)P(x,t)\,\mathrm{d} x=\frac{\mathrm{d} }{\mathrm{d} t}\langle f(X(t)) \rangle =\langle \nabla f(X)\cdot\dot X \rangle =\int\nabla f\cdot\bm{v}\,P\,\mathrm{d} x=-\int f\,\nabla\cdot(P\bm{v})\,\mathrm{d} x

(最後は部分積分)。左辺は ∫f ∂tP dx\int f\,\partial_tP\,\mathrm{d} x なので、ff が任意であることから

∂P∂t+∇⋅(Pv)=0(7.4.2)\frac{\partial P}{\partial t}+\nabla\cdot(P\bm{v})=0\tag{7.4.2}

が成り立つ(∇⋅(Pv)=∑k∇k⋅(Pvk)\nabla\cdot(P\bm{v})=\sum_k\nabla_k\cdot(P\bm{v}_k))。

第2段:ρ=∣ψ∣2\rho=\lvert \psi\rvert ^2 も同じ式に従う。第7部第1章の Madelung 方程式の連続の式 ∂tρ+∇⋅(ρv)=0\partial_t\rho+\nabla\cdot(\rho\bm{v})=0 がそれで、そこで使った速度場と (B3) の速度場は同じものである。

第3段:同じ線形の式に同じ初期値なので両者は等しい。ρ>0\rho>0 の領域で P=hρP=h\rho と置くと

∂(hρ)∂t+∇⋅(hρv)=h(∂ρ∂t+∇⋅(ρv))+ρ(∂h∂t+v⋅∇h)\frac{\partial (h\rho)}{\partial t}+\nabla\cdot(h\rho\bm{v})=h\Bigl(\frac{\partial \rho}{\partial t}+\nabla\cdot(\rho\bm{v})\Bigr)+\rho\Bigl(\frac{\partial h}{\partial t}+\bm{v}\cdot\nabla h\Bigr)

で、左辺と右辺の第1項はともに 00 だから ∂th+v⋅∇h=0\partial_th+\bm{v}\cdot\nabla h=0、つまり hh は軌道に沿って一定である。初めに h=1h=1 なので、軌道が届く所ではずっと h=1h=1、すなわち P=ρP=\rho である。

証明の第1段は「軌道の束が運ぶ分布は、その速度場の連続の式に従う」という一般的な事実で、量子力学をまったく使っていない。量子力学が入るのは第2段で、∣ψ∣2\lvert \psi\rvert ^2 が同じ速度場の連続の式に従うことである。案内方程式の速度場を ∇S/m\nabla S/m に選んだ理由はここにある。同じ連続の式 ∂tρ+∇⋅j=0\partial_t\rho+\nabla\cdot\bm{j}=0 を満たす流れは j\bm{j} に発散のない場を足しても作れるので、等変性だけからは速度場は1つに決まらない。1つに決める理由は第6節の橋で、弱測定で測れる流れという形で述べる。

軌道に沿った言い方にすると、初めの点 x0x_0 から出た軌道を X(t;x0)X(t;x_0) と書き、その Jacobi 行列式を J(x0,t)=det⁡(∂X/∂x0)J(x_0,t)=\det(\partial X/\partial x_0) とすれば

ρ(X(t;x0),t) J(x0,t)=ρ(x0,0)(7.4.3)\rho\bigl(X(t;x_0),t\bigr)\,J(x_0,t)=\rho(x_0,0)\tag{7.4.3}

である。流れに乗った小さな体積の中の確率が保たれる、という意味で、流体力学の Lagrange 描像での質量の保存と同じ式である(第7部第2章が ℏ→0\hbar\to0 の流れで使った式と同じ形。第7部第2章)。

2階の形:量子ポテンシャルの力

案内方程式は速度を与える1階の式だが、時間で微分すると Newton の運動方程式の形になる。1粒子で、軌道に沿った加速度は

X¨=ddtv(X(t),t)=(∂v∂t+(v⋅∇)v)∣r=X\ddot{\bm{X}}=\frac{\mathrm{d} }{\mathrm{d} t}\bm{v}(\bm{X}(t),t)=\Bigl(\frac{\partial \bm{v}}{\partial t}+(\bm{v}\cdot\nabla)\bm{v}\Bigr)\Big|_{\bm{r}=\bm{X}}

で、右辺の括弧は流体力学の物質微分である。第7部第1章の量子 Euler 方程式はこれが −∇(V+Q)/m-\nabla(V+Q)/m に等しいと言っているので、

mX¨=−∇(V+Q)∣r=X,Q=−ℏ22m∇2ρρ(7.4.4)m\ddot{\bm{X}}=-\nabla\bigl(V+Q\bigr)\big|_{\bm{r}=\bm{X}},\qquad Q=-\frac{\hbar^2}{2m}\frac{\nabla^2\sqrt\rho}{\sqrt\rho}\tag{7.4.4}

が成り立つ。Bohm(1952)は理論をこの形で書いた。ただし (7.4.4) は2階の式なので初めの速度を別に与える必要があり、それを mX˙(0)=∇S(X(0),0)m\dot{\bm{X}}(0)=\nabla S(\bm{X}(0),0) に限ったときに限って (7.4.1) と同じ軌道になる。初めの速度を勝手に選べる理論にすると等変性が成り立たず、Born の規則を再現しない。この意味で、1階の案内方程式のほうが要請として素直である。(7.4.4) の利点は、古典力学との違いが「量子ポテンシャル QQ の力が加わる」の1点に集まって見えることにある。

文献:L. de Broglie, J. Phys. Radium 8, 225 (1927);D. Bohm, Phys. Rev. 85, 166 および 180 (1952)。

例題 7.4.1 水素原子の電子は止まっているか ★★★

オリジナル水素原子の電子(質量 mm、Bohr 半径 a=4πε0ℏ2/(me2)a=4\pi\varepsilon_0\hbar^2/(me^2))について、Bohm の軌道を調べる。

  1. 基底状態 ψ1s∝e−r/a\psi_{1s}\propto\mathrm{e}^{-r/a} で、速度場と量子ポテンシャル QQ を求め、V+QV+Q が rr によらないことを示せ。電子に働く力はいくらか。
  2. 2p2p の mℓ=+1m_\ell=+1 の状態 ψ∝r e−r/2asin⁡θ eiφ\psi\propto r\,\mathrm{e}^{-r/2a}\sin\theta\,\mathrm{e}^{\mathrm{i}\varphi} で速度場を求めよ。r=4ar=4a、θ=π/2\theta=\pi/2 にいる電子の速さと1周の時間を数値で答えよ。
解答

(1) ψ1s\psi_{1s} は実数なので SS は定数で、速度場は v=∇S/m=0\bm{v}=\nabla S/m=0。電子はどの位置にいても止まっている。量子ポテンシャルは、球対称な関数のラプラシアン ∇2R=r−2(r2R′)′\nabla^2R=r^{-2}(r^2R')' を R=e−r/aR=\mathrm{e}^{-r/a} に使って

∇2RR=1a2−2ar,Q=−ℏ22ma2+ℏ2mar.\frac{\nabla^2R}{R}=\frac{1}{a^2}-\frac{2}{ar},\qquad Q=-\frac{\hbar^2}{2ma^2}+\frac{\hbar^2}{mar}.

Coulomb ポテンシャルは V=−e2/(4πε0r)V=-e^2/(4\pi\varepsilon_0r) で、Bohr 半径の定義から ℏ2/(ma)=e2/(4πε0)\hbar^2/(ma)=e^2/(4\pi\varepsilon_0) なので、QQ の第2項は VV をちょうど打ち消し、

V+Q=−ℏ22ma2=−13.6 eVV+Q=-\frac{\hbar^2}{2ma^2}=-13.6\,\mathrm{eV}

が rr によらない。(7.4.4) の力 −∇(V+Q)-\nabla(V+Q) は 00 で、止まった電子が止まったままでいることと整合する。古典の電子なら Coulomb 力で核に落ちるところを、量子ポテンシャルの外向きの力がちょうど支えている。基準の言葉では ⟨p^2⟩/(2m)=13.6 eV\langle \hat p^2 \rangle /(2m)=13.6\,\mathrm{eV} の運動エネルギーがあるが、Bohm の力学ではそれは運動の形でなく ∫ρQ\int\rho Q の形で蓄えられている(第7部第1章の「運動エネルギー=流れ+形のエネルギー」の分け方で、流れの分が 00)。

(2) 極座標の勾配は ∇=er∂r+eθr−1∂θ+eφ(rsin⁡θ)−1∂φ\nabla=\bm{e}_r\partial_r+\bm{e}_\theta r^{-1}\partial_\theta+\bm{e}_\varphi(r\sin\theta)^{-1}\partial_\varphi である。∂rψ/ψ\partial_r\psi/\psi と ∂θψ/ψ\partial_\theta\psi/\psi は実数なのでvr=vθ=0v_r=v_\theta=0、∂φψ/ψ=i\partial_\varphi\psi/\psi=\mathrm{i} なので

v=ℏmrsin⁡θ eφ.\bm{v}=\frac{\hbar}{mr\sin\theta}\,\bm{e}_\varphi .

電子は zz 軸の周りを、半径 rsin⁡θr\sin\theta の円の上で一定の速さで回る。角運動量の zz 成分は m(rsin⁡θ)vφ=ℏm(r\sin\theta)v_\varphi=\hbar で、どの軌道でも mℓℏm_\ell\hbar に等しい。r=4ar=4a、θ=π/2\theta=\pi/2 では v=ℏ/(4ma)=5.47×105 m s−1v=\hbar/(4ma)=5.47\times10^5\,\mathrm{m\,s^{-1}}(ℏ/(ma)\hbar/(ma) は微細構造定数と光速の積 αc=2.19×106 m s−1\alpha c=2.19\times10^6\,\mathrm{m\,s^{-1}})、1周の時間は 2π(4a)/v=32πma2/ℏ=2.43×10−15 s2\pi(4a)/v=32\pi ma^2/\hbar=2.43\times10^{-15}\,\mathrm{s} である。同じ 2p2p でも mℓ=0m_\ell=0 の状態 ∝r e−r/2acos⁡θ\propto r\,\mathrm{e}^{-r/2a}\cos\theta は実数なので電子は止まる。どちらの状態も同じエネルギーを持つので、軌道が動くかどうかはエネルギーでなく波動関数の位相で決まる。

木下の勉強ステータス

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

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

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