量子力学の様々な構成方法 第7部 第6章 §7.6.2 ★★★ 消える Brown 運動を、正の固有関数で重みを付け直すと消えない拡散過程になり、それが Nelson の固有状態の過程である。経路積分と Nelson の過程の橋の本体
変換の定義
Feynman–Kac の公式の右辺の期待値は、粒子が消えていくので確率測度の期待値ではない(f = 1 f=1 f = 1 でも 1 1 1 にならない)。正の基底状態 ψ 0 \psi_0 ψ 0 を使うと、これを消えない確率過程に直せる。J. L. Doob が条件付きの Brown 運動の研究で導入した変換で、h h h 変換とも呼ぶ(J. L. Doob (1957), Bulletin de la Société Mathématique de France 85 , 431–458)。
命題 7.6.2(Doob の変換 ★★★)
H ^ ψ 0 = E 0 ψ 0 \hat H\psi_0=E_0\psi_0 H ^ ψ 0 = E 0 ψ 0 、ψ 0 > 0 \psi_0>0 ψ 0 > 0 とし、x x x を出た Brown 運動の経路の汎関数
M τ = ψ 0 ( W τ ) ψ 0 ( x ) exp ( E 0 τ ℏ − 1 ℏ ∫ 0 τ V ( W s ) d s ) (7.6.5) M_\tau=\frac{\psi_0(W_\tau)}{\psi_0(x)}\exp\Bigl(\frac{E_0\tau}{\hbar}-\frac1\hbar\int_0^\tau V(W_s)\,\mathrm{d} s\Bigr)\tag{7.6.5} M τ = ψ 0 ( x ) ψ 0 ( W τ ) exp ( ℏ E 0 τ − ℏ 1 ∫ 0 τ V ( W s ) d s ) ( 7.6.5 )
を考える。このとき E x [ M τ ] = 1 \mathrm E_x[M_\tau]=1 E x [ M τ ] = 1 で、M τ M_\tau M τ を重みとして付け直した経路の測度のもとで W W W は Markov 過程になり、その推移密度は
p τ ( x → y ) = ψ 0 ( y ) ψ 0 ( x ) e E 0 τ / ℏ ⟨ y ∣ e − H ^ τ / ℏ ∣ x ⟩ , (7.6.6) p_\tau(x\to y)=\frac{\psi_0(y)}{\psi_0(x)}\,\mathrm{e}^{E_0\tau/\hbar}\langle y | \mathrm{e}^{-\hat H\tau/\hbar} | x \rangle ,\tag{7.6.6} p τ ( x → y ) = ψ 0 ( x ) ψ 0 ( y ) e E 0 τ /ℏ ⟨ y ∣ e − H ^ τ /ℏ ∣ x ⟩ , ( 7.6.6 )
生成作用素は第7部第5章の L f = ν f ′ ′ + 2 ν ( ψ 0 ′ / ψ 0 ) f ′ \mathcal Lf=\nu f''+2\nu(\psi_0'/\psi_0)f' L f = ν f ′′ + 2 ν ( ψ 0 ′ / ψ 0 ) f ′ である。すなわち付け直した過程は、Nelson の基底状態の過程そのものである。
証明 (7.6.3) に f = ψ 0 f=\psi_0 f = ψ 0 を入れると、左辺は e − H ^ τ / ℏ ψ 0 = e − E 0 τ / ℏ ψ 0 \mathrm{e}^{-\hat H\tau/\hbar}\psi_0=\mathrm{e}^{-E_0\tau/\hbar}\psi_0 e − H ^ τ /ℏ ψ 0 = e − E 0 τ /ℏ ψ 0 なので、E x [ ψ 0 ( W τ ) e − ∫ V d s / ℏ ] = e − E 0 τ / ℏ ψ 0 ( x ) \mathrm E_x[\psi_0(W_\tau)\mathrm{e}^{-\int V\mathrm{d} s/\hbar}]=\mathrm{e}^{-E_0\tau/\hbar}\psi_0(x) E x [ ψ 0 ( W τ ) e − ∫ V d s /ℏ ] = e − E 0 τ /ℏ ψ 0 ( x ) 。両辺を e − E 0 τ / ℏ ψ 0 ( x ) \mathrm{e}^{-E_0\tau/\hbar}\psi_0(x) e − E 0 τ /ℏ ψ 0 ( x ) で割れば E x [ M τ ] = 1 \mathrm E_x[M_\tau]=1 E x [ M τ ] = 1 である。推移密度は、終点を y y y に固定した重み付きの期待値で、定理 7.6.1 の核の形を使えば (7.6.6) になる。Chapman–Kolmogorov の式 ∫ p τ 1 ( x → z ) p τ 2 ( z → y ) d z = p τ 1 + τ 2 ( x → y ) \int p_{\tau_1}(x\to z)p_{\tau_2}(z\to y)\mathrm{d} z=p_{\tau_1+\tau_2}(x\to y) ∫ p τ 1 ( x → z ) p τ 2 ( z → y ) d z = p τ 1 + τ 2 ( x → y ) は、真ん中の ψ 0 ( z ) \psi_0(z) ψ 0 ( z ) が約分され、e − H ^ τ / ℏ \mathrm{e}^{-\hat H\tau/\hbar} e − H ^ τ /ℏ の合成則に帰着するので成り立つ。これで Markov 過程の推移の族になる。生成作用素は τ → 0 \tau\to0 τ → 0 での ∂ τ \partial_\tau ∂ τ で、∂ τ p τ = ψ 0 ( y ) ψ 0 ( x ) − 1 e E 0 τ / ℏ ( E 0 − H ^ ) / ℏ \partial_\tau p_\tau=\psi_0(y)\psi_0(x)^{-1}\mathrm{e}^{E_0\tau/\hbar}(E_0-\hat H)/\hbar ∂ τ p τ = ψ 0 ( y ) ψ 0 ( x ) − 1 e E 0 τ /ℏ ( E 0 − H ^ ) /ℏ を f f f に当てるとL f = − ℏ − 1 ψ 0 − 1 ( H ^ − E 0 ) ( ψ 0 f ) \mathcal Lf=-\hbar^{-1}\psi_0^{-1}(\hat H-E_0)(\psi_0f) L f = − ℏ − 1 ψ 0 − 1 ( H ^ − E 0 ) ( ψ 0 f ) である。第7部第5章の基底状態の変換の命題(第7部第5章 )により、これは ν f ′ ′ + 2 ν ( ψ 0 ′ / ψ 0 ) f ′ \nu f''+2\nu(\psi_0'/\psi_0)f' ν f ′′ + 2 ν ( ψ 0 ′ / ψ 0 ) f ′ に等しい。
ドリフトの項 2 ν ψ 0 ′ / ψ 0 = ν ∂ x ln ψ 0 2 2\nu\psi_0'/\psi_0=\nu\partial_x\ln\psi_0^2 2 ν ψ 0 ′ / ψ 0 = ν ∂ x ln ψ 0 2 は第7部第5章の浸透速度 u u u で、基底状態では流れの速度が v = 0 v=0 v = 0 なので前向きのドリフト b = v + u b=v+u b = v + u に一致する。Feynman–Kac の公式では「ポテンシャルの高い所で消える」という形で入っていたポテンシャルが、Doob の変換のあとでは「ψ 0 \psi_0 ψ 0 の大きい方へ引くドリフト」という形に移る。消える割合 ( V − E 0 ) / ℏ (V-E_0)/\hbar ( V − E 0 ) /ℏ とドリフトの間の関係は、ψ 0 \psi_0 ψ 0 の満たす Schrödinger 方程式 ν ψ 0 ′ ′ = ( V − E 0 ) ψ 0 / ℏ \nu\psi_0''=(V-E_0)\psi_0/\hbar ν ψ 0 ′′ = ( V − E 0 ) ψ 0 /ℏ を通して保たれている。
例題 7.6.2 Mehler 核から Ornstein–Uhlenbeck 過程の推移密度を作る ★★
オリジナル 共通例題①の基底状態 ψ 0 ∝ e − x 2 / ( 4 σ g 2 ) \psi_0\propto\mathrm{e}^{-x^2/(4\sigma_{\mathrm g}^2)} ψ 0 ∝ e − x 2 / ( 4 σ g 2 ) 、E 0 = ℏ ω / 2 E_0=\hbar\omega/2 E 0 = ℏ ω /2 について、(7.6.6) に例題 7.6.1 の Euclid の Mehler 核を入れ、p τ ( x → y ) p_\tau(x\to y) p τ ( x → y ) が平均 x e − ω τ x\mathrm{e}^{-\omega\tau} x e − ω τ ・分散 σ g 2 ( 1 − e − 2 ω τ ) \sigma_{\mathrm g}^2(1-\mathrm{e}^{-2\omega\tau}) σ g 2 ( 1 − e − 2 ω τ ) のガウス分布の密度であることを示せ。この過程のドリフトを読み取れ。
解答 c = cosh ω τ c=\cosh\omega\tau c = cosh ω τ 、s = sinh ω τ s=\sinh\omega\tau s = sinh ω τ 、κ = m ω / ℏ = 1 / ( 2 σ g 2 ) \kappa=m\omega/\hbar=1/(2\sigma_{\mathrm g}^2) κ = mω /ℏ = 1/ ( 2 σ g 2 ) と書く。(7.6.6) の指数は
− κ 2 ( y 2 − x 2 ) + ω τ 2 − κ 2 s [ ( x 2 + y 2 ) c − 2 x y ] -\frac\kappa2\bigl(y^2-x^2\bigr)+\frac{\omega\tau}{2}-\frac{\kappa}{2s}\bigl[(x^2+y^2)c-2xy\bigr] − 2 κ ( y 2 − x 2 ) + 2 ω τ − 2 s κ [ ( x 2 + y 2 ) c − 2 x y ]
で、第1項が ψ 0 ( y ) / ψ 0 ( x ) \psi_0(y)/\psi_0(x) ψ 0 ( y ) / ψ 0 ( x ) 、第2項が e E 0 τ / ℏ \mathrm{e}^{E_0\tau/\hbar} e E 0 τ /ℏ 、第3項が Mehler 核の指数である。y y y の2次の係数は − κ 2 ( 1 + c / s ) = − κ 2 s + c s = − κ 2 e ω τ s -\frac\kappa2(1+c/s)=-\frac{\kappa}{2}\frac{s+c}{s}=-\frac{\kappa}{2}\frac{\mathrm{e}^{\omega\tau}}{s} − 2 κ ( 1 + c / s ) = − 2 κ s s + c = − 2 κ s e ω τ 、x y xy x y の係数は κ / s \kappa/s κ / s 、x 2 x^2 x 2 の係数は κ 2 ( 1 − c / s ) = − κ 2 e − ω τ s \frac\kappa2(1-c/s)=-\frac{\kappa}{2}\frac{\mathrm{e}^{-\omega\tau}}{s} 2 κ ( 1 − c / s ) = − 2 κ s e − ω τ である。2 s = e ω τ ( 1 − e − 2 ω τ ) 2s=\mathrm{e}^{\omega\tau}(1-\mathrm{e}^{-2\omega\tau}) 2 s = e ω τ ( 1 − e − 2 ω τ ) を使うと3つの項は
− κ e ω τ 2 s ( y − x e − ω τ ) 2 = − ( y − x e − ω τ ) 2 2 σ g 2 ( 1 − e − 2 ω τ ) -\frac{\kappa\,\mathrm{e}^{\omega\tau}}{2s}\bigl(y-x\mathrm{e}^{-\omega\tau}\bigr)^2=-\frac{(y-x\mathrm{e}^{-\omega\tau})^2}{2\sigma_{\mathrm g}^2(1-\mathrm{e}^{-2\omega\tau})} − 2 s κ e ω τ ( y − x e − ω τ ) 2 = − 2 σ g 2 ( 1 − e − 2 ω τ ) ( y − x e − ω τ ) 2
にまとまる(κ e ω τ / s = 2 κ / ( 1 − e − 2 ω τ ) \kappa\mathrm{e}^{\omega\tau}/s=2\kappa/(1-\mathrm{e}^{-2\omega\tau}) κ e ω τ / s = 2 κ / ( 1 − e − 2 ω τ ) )。前因子は κ / ( π s ) e ω τ / 2 = 2 κ e ω τ / ( 2 π s ) = [ 2 π σ g 2 ( 1 − e − 2 ω τ ) ] − 1 / 2 \sqrt{\kappa/(\pi s)}\,\mathrm{e}^{\omega\tau/2}=\sqrt{2\kappa\mathrm{e}^{\omega\tau}/(2\pi s)}=[2\pi\sigma_{\mathrm g}^2(1-\mathrm{e}^{-2\omega\tau})]^{-1/2} κ / ( π s ) e ω τ /2 = 2 κ e ω τ / ( 2 π s ) = [ 2 π σ g 2 ( 1 − e − 2 ω τ ) ] − 1/2 で、ちょうどガウス分布の規格化になる。したがって p τ p_\tau p τ は平均 x e − ω τ x\mathrm{e}^{-\omega\tau} x e − ω τ ・分散 σ g 2 ( 1 − e − 2 ω τ ) \sigma_{\mathrm g}^2(1-\mathrm{e}^{-2\omega\tau}) σ g 2 ( 1 − e − 2 ω τ ) のガウス分布である。
平均の短時間の変化は x ( e − ω τ − 1 ) ≈ − ω x τ x(\mathrm{e}^{-\omega\tau}-1)\approx-\omega x\tau x ( e − ω τ − 1 ) ≈ − ω xτ なのでドリフトは − ω x -\omega x − ω x 、分散の増え方は σ g 2 ⋅ 2 ω τ = ( ℏ / m ) τ = 2 ν τ \sigma_{\mathrm g}^2\cdot2\omega\tau=(\hbar/m)\tau=2\nu\tau σ g 2 ⋅ 2 ω τ = ( ℏ/ m ) τ = 2 ν τ で、第7部第5章の Ornstein–Uhlenbeck 過程 d X = − ω X d t + d W \mathrm{d} X=-\omega X\mathrm{d} t+\mathrm{d} W d X = − ω X d t + d W に一致する。ψ 0 \psi_0 ψ 0 による付け直しで、Mehler 核の中の量子力学の情報がドリフト − ω x -\omega x − ω x に移っている。
逆向き:過程からハミルトニアンを取り戻す
(7.6.6) は逆にも読める。基底状態の Nelson の過程の推移密度 p τ p_\tau p τ と定常分布 ρ 0 = ψ 0 2 \rho_0=\psi_0^2 ρ 0 = ψ 0 2 が分かれば、
⟨ y ∣ e − H ^ τ / ℏ ∣ x ⟩ = e − E 0 τ / ℏ ρ 0 ( x ) ρ 0 ( y ) p τ ( x → y ) (7.6.7) \langle y | \mathrm{e}^{-\hat H\tau/\hbar} | x \rangle =\mathrm{e}^{-E_0\tau/\hbar}\sqrt{\frac{\rho_0(x)}{\rho_0(y)}}\;p_\tau(x\to y)\tag{7.6.7} ⟨ y ∣ e − H ^ τ /ℏ ∣ x ⟩ = e − E 0 τ /ℏ ρ 0 ( y ) ρ 0 ( x ) p τ ( x → y ) ( 7.6.7 )
で虚時間の核が決まり、τ \tau τ で微分して H ^ \hat H H ^ が(定数 E 0 E_0 E 0 のずれを除いて)決まる。E 0 E_0 E 0 は ψ 0 \psi_0 ψ 0 の満たす方程式から決まるので、ずれも残らない。したがって、正の基底状態を持つ系では、基底状態の Nelson の過程と、虚時間の経路積分と、基準のハミルトニアンは、互いに作り直せる。第6節の橋はこの事実にもとづく。
励起状態:節で区切られた領域ごとの Doob の変換
励起状態 ψ n \psi_n ψ n は符号を変えるので、そのまま (7.6.5) に使えない。ψ n \psi_n ψ n が 0 0 0 にならない連結な領域(節の領域)D D D ごとに考える。D D D の上で ψ n \psi_n ψ n は符号を変えず、D D D の境界で 0 0 0 になる。D D D の境界で 0 0 0 になる条件(Dirichlet 条件)を課したハミルトニアン H ^ D \hat H_D H ^ D の固有関数のうち、符号を変えないものは最も低い固有値のものだけなので、ψ n ∣ D \psi_n|_D ψ n ∣ D は H ^ D \hat H_D H ^ D の基底状態で、固有値はもとの E n E_n E n である。H ^ D \hat H_D H ^ D の Feynman–Kac の公式は、(7.6.3) の Brown 運動を「D D D の境界に触れた瞬間に消す」ことで得られる。それを ψ n ∣ D \psi_n|_D ψ n ∣ D で Doob 変換すると、D D D の中で消えない過程になり、ドリフトは 2 ν ∇ ψ n / ψ n 2\nu\nabla\psi_n/\psi_n 2 ν ∇ ψ n / ψ n で、第7部第5章の固有状態の Nelson の過程に一致する。境界の近くで ψ n \psi_n ψ n は境界までの距離に比例して 0 0 0 になるので、ドリフトは境界から押し返す 2 ν / ( 距離 ) 2\nu/(\text{距離}) 2 ν / ( 距離 ) の形になり、過程は D D D の外へ出られない。第7部第5章では節を越えないことを Bessel 過程との比較で示したが、Doob の変換の言葉では「Dirichlet 条件で消える Brown 運動を、消えなかった経路だけに条件付けたもの」と読める。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com