tecktom.com

逆向き:流体の方程式から波動関数を作り直す

2026年10月17日9分で読めますAI生成教科書『虎の巻』第7部 第1章 波動関数を密度と位相に分ける:Madelung 変換と量子ポテンシャル
目次
  1. 同値性の定理
  2. 条件 (iv) は流体の言葉の外にある
  3. 木下の層

AI 生成・木下が未勉強

量子力学の様々な構成方法 第7部 第1章 §7.1.4 ★★★ 同値性の橋の後半。流体の方程式だけでは足りず、循環の量子化条件を要請として加えたときに限って Schrödinger 方程式に戻ることを証明する

第2節で、Schrödinger 方程式から Madelung の流体方程式が出ることを示した。逆に、流体の方程式を満たす (ρ,v)(\rho,\bm{v}) を与えたとき、∣ψ∣2=ρ\lvert \psi\rvert ^2=\rho、j=ρv\bm{j}=\rho\bm{v} となる Schrödinger 方程式の解 ψ\psi が必ず作れるかを考える。答えは「そのままでは作れない」で、足りない条件を明示したのが次の定理である。

同値性の定理

定理 7.1.4(Madelung 方程式と Schrödinger 方程式の同値性 ★★★)

D⊂R3D\subset\mathbb{R}^3 を開集合とし、DD の上の滑らかな ρ>0\rho>0 と v\bm{v} が各時刻で次の4つを満たすとする。

  1. (i) 連続の式 ∂tρ+∇⋅(ρv)=0\partial_t\rho+\nabla\cdot(\rho\bm{v})=0。
  2. (ii) 量子 Euler 方程式 (7.1.9)(QQ は (7.1.7) の形)。
  3. (iii) 渦なし ∇×v=0\nabla\times\bm{v}=0。
  4. (iv) 循環の量子化条件:DD の中のどの閉曲線 CC についても

∮Cmv⋅dl∈2πℏ Z.(7.1.14)\oint_C m\bm{v}\cdot\mathrm{d}\bm{l}\in2\pi\hbar\,\mathbb{Z} .\tag{7.1.14}

このとき、DD の上の滑らかな関数 ψ\psi で、∣ψ∣2=ρ\lvert \psi\rvert ^2=\rho、(ℏ/m)Im⁡(∇ψ/ψ)=v(\hbar/m)\operatorname{Im}(\nabla\psi/\psi)=\bm{v} を満たし、Schrödinger 方程式 (7.1.5) を解くものが、全体の位相の定数倍を除いてただ1つある。逆に、Schrödinger 方程式の解から Madelung 変換で作った (ρ,v)(\rho,\bm{v}) は、ρ>0\rho>0 の領域で (i)〜(iv) を満たす。

証明

Schrödinger から流体へ。 (i) と (ii) は第2節と第3節で導いた。(iii) は v=∇S/m\bm{v}=\nabla S/m が局所的に勾配であることから出る。(iv) は、ψ\psi が DD の上で1価で 00 にならないので、閉曲線 CC に沿って1周したとき偏角 S/ℏS/\hbar の増え方は 2π2\pi の整数倍(ψ\psi の値が 00 の周りを回る回数)になり、∮C∇S⋅dl=2πℏn\oint_C\nabla S\cdot\mathrm{d}\bm{l}=2\pi\hbar n となることから出る。

流体から Schrödinger へ。 時刻 tt を止め、DD の連結成分の1つに基準点 r0\bm{r}_0 を選んで

S(r,t)=∫r0rmv(r′,t)⋅dr′+s(t)S(\bm{r},t)=\int_{\bm{r}_0}^{\bm{r}}m\bm{v}(\bm{r}',t)\cdot\mathrm{d}\bm{r}'+s(t)

と置く(s(t)s(t) は後で決める実数の関数。DD が複数の連結成分を持てば成分ごとに同じことをする)。(iii) より、積分の値は端点を固定したまま経路を連続的に変形しても変わらない。変形で移り合えない2本の経路では、値の差は閉曲線の上の積分なので、(iv) より 2πℏ2\pi\hbar の整数倍だけ違う。したがって SS そのものは経路によって2πℏ2\pi\hbar の整数倍だけ不定だが、eiS/ℏ\mathrm{e}^{\mathrm{i} S/\hbar} は経路によらず定まる。この SS は ∇S=mv\nabla S=m\bm{v} を満たす。

次に ψ=ρ eiS/ℏ\psi=\sqrt\rho\,\mathrm{e}^{\mathrm{i} S/\hbar} と置いて、量子 HJ 方程式が成り立つように s(t)s(t) を選ぶ。F(r,t)=∂tS+∣∇S∣2/(2m)+V+QF(\bm{r},t)=\partial_tS+\lvert \nabla S\rvert ^2/(2m)+V+Q と置くと、∂t∇S=m ∂tv\partial_t\nabla S=m\,\partial_t\bm{v} と (iii) から (v⋅∇)v=∇(∣v∣2/2)(\bm{v}\cdot\nabla)\bm{v}=\nabla(\lvert \bm{v}\rvert ^2/2) が使えるので、

∇F=m(∂v∂t+(v⋅∇)v)+∇(V+Q)=0\nabla F=m\Bigl(\frac{\partial \bm{v}}{\partial t}+(\bm{v}\cdot\nabla)\bm{v}\Bigr)+\nabla(V+Q)=0

が (ii) から出る。よって FF は r\bm{r} によらず F=f(t)F=f(t) である。s(t)s(t) を s′(t)s'(t) だけずらすと FF も同じだけずれるので、s(t)=−∫tf(t′) dt′s(t)=-\int^tf(t')\,\mathrm{d} t' と選び直せば F=0F=0、つまり量子 HJ 方程式 (7.1.7) が成り立つ。連続の式は (i) そのものである。命題 7.1.2 の証明の恒等式から、ψ\psi は Schrödinger 方程式を満たす。

一意性は、別の解 ψ′\psi' も同じ ρ\rho と v\bm{v} を与えるなら、ψ′/ψ\psi'/\psi は絶対値 11 で偏角の勾配が 00 なので各連結成分で eiα(t)\mathrm{e}^{\mathrm{i}\alpha(t)} の形になり、両方が Schrödinger 方程式を満たすことから α\alpha は時間にもよらない定数になることで示せる。

条件 (iv) は流体の言葉の外にある

(i)〜(iii) は古典の流体の方程式にもある種類の条件で、流体の言葉で自然に書ける。ところが (iv) は、循環を 2πℏ/m2\pi\hbar/m の整数倍に限るという条件で、流体の方程式から導くことはできない。DD が単連結(どの閉曲線も1点に縮められる)なら (iii) から循環は全部 00 になり (iv) は自動的に満たされるが、波動関数に節の線があると DD はその線を除いた領域になって単連結でなくなり、(iv) が独立な条件として必要になる。

T. C. Wallstrom(1994)は、この点を指摘して「Madelung の方程式は、量子化条件を手で加えない限り Schrödinger 方程式と同値でない」と論じた(T. C. Wallstrom, Phys. Rev. A 49, 1613 (1994))。Madelung の原論文は E. Madelung, Z. Phys. 40, 322 (1927) である。流体を出発点に置く構成では、(iv) を要請の1つとして最初から置くことになる。この章の要請を並べ直すと、第7部の出発点は次のようになる。

第7部の出発点の要請

状態は、配位空間の上の密度 ρ≥0\rho\ge0(∫ρ=1\int\rho=1)と速度場 v\bm{v} の組で表す。ρ>0\rho>0 の所で v\bm{v} は渦なしで、どの閉曲線の上の循環も 2πℏ/m2\pi\hbar/m の整数倍である。時間発展は連続の式 (7.1.6) と量子 Euler 方程式 (7.1.9) で決まる。定理 7.1.4 により、この要請は ρ>0\rho>0 の領域で第1部の基準(位置表示の Schrödinger 方程式)と同値である。

節そのものの上での扱いには、もう1つ注意が要る。定理は ρ>0\rho>0 の領域 DD の上だけで述べてあり、ψ=ρ eiS/ℏ\psi=\sqrt\rho\,\mathrm{e}^{\mathrm{i} S/\hbar} を節の上まで滑らかにつなげられるかは別に確かめなければならない。節の近くで ρ\sqrt\rho が距離に比例して 00 になり、位相の巻き数が循環の整数に一致していればψ\psi は節を越えて滑らかにつながる。第5節の調和振動子の第1励起状態は、1次元で節が1点だけの例である。

電磁場があるとき

ベクトルポテンシャル A\bm{A} の中の電荷 qq の粒子では、運動量の演算子が p^−qA\hat{\bm{p}}-q\bm{A} に置き換わるので、速度場は v=(∇S−qA)/m\bm{v}=(\nabla S-q\bm{A})/m になり、渦なしの条件は ∇×(mv)=−qB\nabla\times(m\bm{v})=-q\bm{B} に変わる。量子化条件は ∮C(mv+qA)⋅dl∈2πℏZ\oint_C(m\bm{v}+q\bm{A})\cdot\mathrm{d}\bm{l}\in2\pi\hbar\mathbb{Z} で、超伝導の磁束の量子化と同じ形になる。この形は第7部第3章(第7部第3章)と第10部第3章の前量子化の整数条件(第10部第3章)でもう一度現れる。

例題 7.1.4 輪の上を回る流れと、量子化条件の無い流体 ★★★

オリジナル半径 R0R_0 の細い輪の上に閉じ込めた質量 mm の粒子を考える。輪の上の位置を角度 φ\varphi で表し、外のポテンシャルは無いとする。

  1. 一様な密度 ρ=1/(2πR0)\rho=1/(2\pi R_0) と、輪に沿った一様な速さ v=νℏ/(mR0)v=\nu\hbar/(mR_0)(ν\nu は実数)の組が、(i)〜(iii) を満たすことを示し、そのエネルギーを求めよ。
  2. 量子化条件 (iv) を課すと ν\nu に何が起きるか。ν=1/2\nu=1/2 の流体の解を Madelung 変換で波動関数に戻そうとすると、何が起きるか。
解答

(1) ρ\rho も vv も φ\varphi と tt によらないので、連続の式の各項は 00 である。1次元の輪の上なので渦の条件 (iii) は意味を持たず自動的に満たされる。ρ\rho が一様なので Q=0Q=0、vv が一様で時間にもよらないので Euler 方程式の両辺は 00 である。エネルギーは流れの運動エネルギーだけでE=mv2/2=ℏ2ν2/(2mR02)E=mv^2/2=\hbar^2\nu^2/(2mR_0^2) となる。ν\nu は任意の実数でよいので、流体の方程式 (i)〜(iii) だけではエネルギーは連続的な値をとる。

(2) 輪を1周する循環は ∮mv dl=2πR0⋅νℏ/R0=2πℏν\oint mv\,\mathrm{d} l=2\pi R_0\cdot\nu\hbar/R_0=2\pi\hbar\nu なので、(iv) は ν∈Z\nu\in\mathbb{Z} を要求する。これで E=ℏ2n2/(2mR02)E=\hbar^2n^2/(2mR_0^2)(n∈Zn\in\mathbb{Z})という、輪の上の Schrödinger 方程式の固有値と同じ離散的な値になる。ν=1/2\nu=1/2 の解を戻すと、S=ℏφ/2−EtS=\hbar\varphi/2-Et より ψ∝eiφ/2\psi\propto\mathrm{e}^{\mathrm{i}\varphi/2} だが、φ=2π\varphi=2\pi で ψ=−ψ(φ=0)\psi=-\psi(\varphi=0) となって1周で値が戻らない。1価の波動関数にはならず、エネルギー ℏ2/(8mR02)\hbar^2/(8mR_0^2) は量子力学の 00 と ℏ2/(2mR02)\hbar^2/(2mR_0^2) の間の、許されない値である。エネルギーの量子化は、流体の方程式からでなく条件 (iv) から出てくる。

木下の勉強ステータス

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

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

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