tecktom.com

Schrödinger 方程式を実部と虚部に分ける

2026年10月17日7分で読めますAI生成教科書『虎の巻』第7部 第1章 波動関数を密度と位相に分ける:Madelung 変換と量子ポテンシャル
目次
  1. 代入して整理する
  2. 量子ポテンシャルの形
  3. 木下の層

AI 生成・木下が未勉強

量子力学の様々な構成方法 第7部 第1章 §7.1.2 ★★★ 同値性の橋の前半。連続の式と、量子ポテンシャルつきの Hamilton–Jacobi 方程式を導く

代入して整理する

Schrödinger 方程式

iℏ∂ψ∂t=−ℏ22m∇2ψ+V(r)ψ(7.1.5)\mathrm{i}\hbar\frac{\partial \psi}{\partial t}=-\frac{\hbar^2}{2m}\nabla^2\psi+V(\bm{r})\psi\tag{7.1.5}

に極形式 (7.1.1) を入れる。書きやすいように R=ρR=\sqrt\rho と置き、ψ=R eiS/ℏ\psi=R\,\mathrm{e}^{\mathrm{i} S/\hbar} とする。時間微分は

iℏ∂ψ∂t=(iℏ∂R∂t−R∂S∂t)eiS/ℏ.\mathrm{i}\hbar\frac{\partial \psi}{\partial t}=\Bigl(\mathrm{i}\hbar\frac{\partial R}{\partial t}-R\frac{\partial S}{\partial t}\Bigr)\mathrm{e}^{\mathrm{i} S/\hbar}.

空間の微分は、まず ∇ψ=(∇R+iℏR∇S)eiS/ℏ\nabla\psi=\bigl(\nabla R+\tfrac{\mathrm{i}}{\hbar}R\nabla S\bigr)\mathrm{e}^{\mathrm{i} S/\hbar}、もう1回発散を取って

∇2ψ=(∇2R+2iℏ∇R⋅∇S+iℏR∇2S−1ℏ2R∣∇S∣2)eiS/ℏ.\nabla^2\psi=\Bigl(\nabla^2R+\frac{2\mathrm{i}}{\hbar}\nabla R\cdot\nabla S+\frac{\mathrm{i}}{\hbar}R\nabla^2S-\frac{1}{\hbar^2}R\lvert \nabla S\rvert ^2\Bigr)\mathrm{e}^{\mathrm{i} S/\hbar}.

両辺から共通の eiS/ℏ\mathrm{e}^{\mathrm{i} S/\hbar} を除き、実部と虚部に分ける。RR と SS は実数なので、実部と虚部がそれぞれ等しくなければならない。実部からは

−R∂S∂t=−ℏ22m∇2R+12mR∣∇S∣2+VR,-R\frac{\partial S}{\partial t}=-\frac{\hbar^2}{2m}\nabla^2R+\frac{1}{2m}R\lvert \nabla S\rvert ^2+VR,

虚部(ℏ\hbar で割ったもの)からは

∂R∂t=−12m(2∇R⋅∇S+R∇2S)\frac{\partial R}{\partial t}=-\frac{1}{2m}\bigl(2\nabla R\cdot\nabla S+R\nabla^2S\bigr)

が出る。実部の式を −R-R で割り、虚部の式に 2R2R を掛けて ρ=R2\rho=R^2 に直すと、次の2本になる。

命題 7.1.2(Madelung 方程式 ★★★)

ρ>0\rho>0 の領域で、Schrödinger 方程式 (7.1.5) は次の2本の実数の式と同値である。

∂ρ∂t+∇⋅(ρ ∇Sm)=0,∂S∂t+∣∇S∣22m+V+Q=0,Q≔−ℏ22m∇2ρρ.\begin{align*}&\frac{\partial \rho}{\partial t}+\nabla\cdot\Bigl(\rho\,\frac{\nabla S}{m}\Bigr)=0,\tag{7.1.6}\\&\frac{\partial S}{\partial t}+\frac{\lvert \nabla S\rvert ^2}{2m}+V+Q=0,\qquad Q\coloneqq-\frac{\hbar^2}{2m}\frac{\nabla^2\sqrt\rho}{\sqrt\rho}.\tag{7.1.7}\end{align*}

1本目は確率の連続の式、2本目は量子 Hamilton–Jacobi 方程式で、QQ を量子ポテンシャルと呼ぶ。

証明

上の計算で、Schrödinger 方程式 × e−iS/ℏ\times\,\mathrm{e}^{-\mathrm{i} S/\hbar} は

−R(∂S∂t+∣∇S∣22m+V+Q)+i ℏ2R(∂ρ∂t+∇⋅ρ∇Sm)=0-R\Bigl(\frac{\partial S}{\partial t}+\frac{\lvert \nabla S\rvert ^2}{2m}+V+Q\Bigr)+\mathrm{i}\,\frac{\hbar}{2R}\Bigl(\frac{\partial \rho}{\partial t}+\nabla\cdot\frac{\rho\nabla S}{m}\Bigr)=0

と書ける(検算で一般の R,S,VR,S,V について確かめた)。R>0R>0 なので、左辺が 00 であることと2つの括弧が両方 00 であることは同じである。

(7.1.6) の中の ρ∇S/m\rho\nabla S/m は (7.1.3) の j\bm{j} なので、連続の式は ∂tρ+∇⋅j=0\partial_t\rho+\nabla\cdot\bm{j}=0 と書ける。これは確率の保存を表し、全空間で積分すれば ∫ρ d3r\int\rho\,\mathrm{d}^3r が時間によらないことが出る。量子ポテンシャルを含むのは2本目だけである。

量子ポテンシャルの形

(7.1.7) の QQ を除くと、古典力学の Hamilton–Jacobi 方程式 ∂tS+∣∇S∣2/(2m)+V=0\partial_tS+\lvert \nabla S\rvert ^2/(2m)+V=0 そのものになる(解析力学から様々な理論分野を組み立てる 第1部第13章。そちらでは同じ極形式を ψ=AeiS/ℏ\psi=A\mathrm{e}^{\mathrm{i} S/\hbar} と書いて同じ2本を出している)。QQ は ℏ2\hbar^2 に比例し、密度 ρ\rho の形だけで決まり、位相 SS にも外のポテンシャル VV にもよらない。ρ\rho を直接使って書くと

Q=−ℏ24m[∇2ρρ−∣∇ρ∣22ρ2](7.1.8)Q=-\frac{\hbar^2}{4m}\Bigl[\frac{\nabla^2\rho}{\rho}-\frac{\lvert \nabla\rho\rvert ^2}{2\rho^2}\Bigr]\tag{7.1.8}

である(∇2ρ=∇⋅(∇ρ/(2ρ))\nabla^2\sqrt\rho=\nabla\cdot(\nabla\rho/(2\sqrt\rho)) を展開すれば出る)。この形から、QQ について次の3つが読み取れる。

1つ目に、ρ\rho を定数倍しても QQ は変わらない。量子ポテンシャルは密度の大きさでなく、密度の曲がり方だけを見ている。古典の流体の圧力が密度の大きさに比例するのと違う点である。2つ目に、密度が遠くまで薄く広がっていても、形が急に変わる所では QQ が大きくなる。この性質は第7部第4章の Bohm の軌道で非局所性の話につながる。3つ目に、ρ→0\rho\to0 になる節の近くでは ρ\sqrt\rho で割るので QQ が発散しうる。第5節で調和振動子の第1励起状態を例にこれを見る。

例題 7.1.2 ガウス分布の量子ポテンシャル ★★★

オリジナル1次元で、密度が平均 00・標準偏差 σ\sigma のガウス分布 ρ(x)=(2πσ2)−1/2e−x2/(2σ2)\rho(x)=(2\pi\sigma^2)^{-1/2}\mathrm{e}^{-x^2/(2\sigma^2)} のとき、量子ポテンシャル Q(x)Q(x) を求めよ。QQ が正になる範囲と、QQ が粒子に及ぼす力 −∂xQ-\partial_xQ の向きを答えよ。

解答

ρ∝e−x2/(4σ2)\sqrt\rho\propto\mathrm{e}^{-x^2/(4\sigma^2)} である。微分すると(ρ)′=−x2σ2ρ(\sqrt\rho)'=-\dfrac{x}{2\sigma^2}\sqrt\rho、(ρ)′′=(x24σ4−12σ2)ρ(\sqrt\rho)''=\Bigl(\dfrac{x^2}{4\sigma^4}-\dfrac{1}{2\sigma^2}\Bigr)\sqrt\rho なので、

Q(x)=−ℏ22m(x24σ4−12σ2)=ℏ24mσ2−ℏ2x28mσ4.Q(x)=-\frac{\hbar^2}{2m}\Bigl(\frac{x^2}{4\sigma^4}-\frac{1}{2\sigma^2}\Bigr)=\frac{\hbar^2}{4m\sigma^2}-\frac{\hbar^2x^2}{8m\sigma^4}.

QQ は上に凸の放物線で、∣x∣<2 σ\lvert x\rvert <\sqrt2\,\sigma で正、その外で負になる。力は −∂xQ=ℏ2x/(4mσ4)-\partial_xQ=\hbar^2x/(4m\sigma^4) で、x>0x>0 では正の向き、x<0x<0 では負の向きを向く。つまり量子ポテンシャルの力は分布の中心から外へ向かい、密度を広げる向きに働く。外のポテンシャルが無いとき(第5節の自由粒子)、波束が広がる理由はこの力で説明できる。規格化の定数 (2πσ2)−1/2(2\pi\sigma^2)^{-1/2} は QQ に現れず、QQ が密度の大きさによらないことの例にもなっている。

木下の勉強ステータス

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

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

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