量子力学の様々な構成方法 第7部 第5章 §7.5.4 ★★★ 第7部第1章から回ってきた共通例題①の熱平衡を扱い、Nelson の過程で混合状態を書く方法と、その方法がアンサンブルの選び方に依存することを示す
1つの定常な過程では書けない
共通例題①の熱平衡状態 ρ ^ = e − β H ^ / Z \hat\rho=\mathrm{e}^{-\beta\hat H}/Z ρ ^ = e − β H ^ / Z の位置の分布は、分散 σ g 2 coth ( β ℏ ω / 2 ) \sigma_{\mathrm g}^2\coth(\beta\hbar\omega/2) σ g 2 coth ( β ℏ ω /2 ) のガウス分布である(第4部第4章 ・第1部第1章の例題)。第7部第1章では、混合状態は1つの波動関数で書けないので、1つの流体 ( ρ , v ) (\rho,v) ( ρ , v ) でも書けないと述べた。Nelson の過程の言葉では、この主張は次の形で証明できる。
命題 7.5.4(1つの定常な Nelson の過程は固有状態しか表せない ★★★)
1次元のポテンシャル V V V の中で、密度 ρ > 0 \rho>0 ρ > 0 が時刻によらず、ρ v → 0 \rho v\to0 ρ v → 0 (∣ x ∣ → ∞ \lvert x\rvert \to\infty ∣ x ∣ → ∞ )を満たす Nelson の過程を考える。このとき v = 0 v=0 v = 0 で、ρ \sqrt\rho ρ はハミルトニアンの固有値 E E E の固有関数である。調和振動子では ρ = ∣ ψ 0 ∣ 2 \rho=\lvert \psi_0\rvert ^2 ρ = ∣ ψ 0 ∣ 2 (基底状態)に限られ、T > 0 T>0 T > 0 の熱平衡の分布を持つ Nelson の過程は存在しない。
証明 連続の式 (7.5.8) で ∂ t ρ = 0 \partial_t\rho=0 ∂ t ρ = 0 なので ρ v \rho v ρ v は x x x によらない定数で、無限遠で 0 0 0 だから v = 0 v=0 v = 0 。(N2) と補題 7.5.2 から 0 = − ∂ x ( V + Q ) / m 0=-\partial_x(V+Q)/m 0 = − ∂ x ( V + Q ) / m 、すなわち V + Q = E V+Q=E V + Q = E (定数)。Q Q Q の定義に戻すと− ℏ 2 2 m ( ρ ) ′ ′ + V ρ = E ρ -\frac{\hbar^2}{2m}(\sqrt\rho)''+V\sqrt\rho=E\sqrt\rho − 2 m ℏ 2 ( ρ ) ′′ + V ρ = E ρ で、ρ \sqrt\rho ρ は固有関数である。ρ > 0 \sqrt\rho>0 ρ > 0 は節を持たないので、1次元では基底状態に限られる(節の数が n n n の固有関数が第 n n n 励起状態)。調和振動子の基底状態の分散は σ g 2 \sigma_{\mathrm g}^2 σ g 2 で、T > 0 T>0 T > 0 の熱平衡の分散 σ g 2 coth ( β ℏ ω / 2 ) \sigma_{\mathrm g}^2\coth(\beta\hbar\omega/2) σ g 2 coth ( β ℏ ω /2 ) はこれより大きいので一致しない。
熱平衡のガウス分布を定常分布に持つ拡散過程そのものは作れる。分散 s 2 = σ g 2 coth ( β ℏ ω / 2 ) s^2=\sigma_{\mathrm g}^2\coth(\beta\hbar\omega/2) s 2 = σ g 2 coth ( β ℏ ω /2 ) を持つ Ornstein–Uhlenbeck 過程d X = − γ X d t + d W \mathrm{d} X=-\gamma X\,\mathrm{d} t+\mathrm{d} W d X = − γ X d t + d W は、ν / γ = s 2 \nu/\gamma=s^2 ν / γ = s 2 から
γ = ω tanh β ℏ ω 2 (7.5.23) \gamma=\omega\tanh\frac{\beta\hbar\omega}{2}\tag{7.5.23} γ = ω tanh 2 β ℏ ω ( 7.5.23 )
と取ればよい。しかしこの過程の平均の加速度は a = − γ 2 x a=-\gamma^2x a = − γ 2 x で、− ω 2 x -\omega^2x − ω 2 x と違う。Nelson の法則で読むと、この過程は角周波数 γ < ω \gamma<\omega γ < ω の柔らかい振動子の基底状態 を表していて、もとの振動子の熱平衡ではない。分布を合わせると力学が合わなくなる。
混合状態は過程の混合として書く
Nelson の構成で混合状態を書くには、ρ ^ \hat\rho ρ ^ を純粋状態の混合 ρ ^ = ∑ k p k ∣ ψ k ⟩ ⟨ ψ k ∣ \hat\rho=\sum_kp_k| \psi_k \rangle \langle \psi_k | ρ ^ = ∑ k p k ∣ ψ k ⟩ ⟨ ψ k ∣ に分け、各 ψ k \psi_k ψ k の過程を確率 p k p_k p k で選ぶ。各成分で等変性が成り立つので、混合した過程の各時刻の密度は ∑ k p k ∣ ψ k ∣ 2 = ⟨ x ∣ ρ ^ ∣ x ⟩ \sum_kp_k\lvert \psi_k\rvert ^2=\langle x | \hat\rho | x \rangle ∑ k p k ∣ ψ k ∣ 2 = ⟨ x ∣ ρ ^ ∣ x ⟩ で、基準の予言と一致する。熱平衡の調和振動子では、次の2つの分け方が自然である(第7部第4章の総合問題と同じ分け方)。
(A) 固有状態 ∣ n ⟩ | n \rangle ∣ n ⟩ を重み p n = ( 1 − e − β ℏ ω ) e − n β ℏ ω p_n=(1-\mathrm{e}^{-\beta\hbar\omega})\mathrm{e}^{-n\beta\hbar\omega} p n = ( 1 − e − β ℏ ω ) e − n β ℏ ω で混ぜる。各成分は定常な過程で、n ≥ 1 n\ge1 n ≥ 1 では節で半直線に分かれる。
(B) コヒーレント状態を重み P ( α ) = e − ∣ α ∣ 2 / n ˉ / ( π n ˉ ) P(\alpha)=\mathrm{e}^{-\lvert \alpha\rvert ^2/\bar n}/(\pi\bar n) P ( α ) = e − ∣ α ∣ 2 / n ˉ / ( π n ˉ ) で混ぜる(n ˉ = 1 / ( e β ℏ ω − 1 ) \bar n=1/(\mathrm{e}^{\beta\hbar\omega}-1) n ˉ = 1/ ( e β ℏ ω − 1 ) 。光の熱状態の P 表示と同じ形)。各成分は (7.5.20) の「古典の振動+Ornstein–Uhlenbeck の揺れ」で、古典の振幅が熱的にばらつく。
どちらでも1時刻の分布は熱平衡のガウス分布になる。(B) では X = X c + ξ X=X_{\mathrm c}+\xi X = X c + ξ で、X c X_{\mathrm c} X c の分散が 2 n ˉ σ g 2 2\bar n\sigma_{\mathrm g}^2 2 n ˉ σ g 2 、ξ \xi ξ の分散が σ g 2 \sigma_{\mathrm g}^2 σ g 2 、両者が独立なので和は σ g 2 ( 1 + 2 n ˉ ) = σ g 2 coth ( β ℏ ω / 2 ) \sigma_{\mathrm g}^2(1+2\bar n)=\sigma_{\mathrm g}^2\coth(\beta\hbar\omega/2) σ g 2 ( 1 + 2 n ˉ ) = σ g 2 coth ( β ℏ ω /2 ) である。2つの分け方は同じ ρ ^ \hat\rho ρ ^ を与え、第4部第3章の意味で区別できない(第4部第3章 )。ところが Nelson の道筋は、(A) と (B) で別物である。(A) では粒子は1つの固有状態の過程に閉じ込められ、(B) では古典の振動をする。1時刻の分布だけを見る限りこの違いは現れないが、多時刻の相関には現れる。この違いと、それが基準の予言とどう食い違うかは第7部第6章で扱う。
例題 7.5.4 熱平衡の分布に合わせた過程は、柔らかい振動子を表す ★★
オリジナル
β ℏ ω = 1 \beta\hbar\omega=1 β ℏ ω = 1 のとき、熱平衡の位置の分散を σ g 2 \sigma_{\mathrm g}^2 σ g 2 の倍数で求め、それを定常分布に持つ Ornstein–Uhlenbeck 過程の γ / ω \gamma/\omega γ / ω を求めよ。
例題 7.5.3 の 87 R b {}^{87}\mathrm{Rb} 87 Rb (ω = 2 π × 100 H z \omega=2\pi\times100\,\mathrm{Hz} ω = 2 π × 100 Hz )を T = 10 n K T=10\,\mathrm{nK} T = 10 nK の熱平衡に置いたときの同じ量を求めよ。
(1) の過程の平均の加速度が − ω 2 x -\omega^2x − ω 2 x からどれだけずれるかを、ある位置 x x x での比 a / ( − ω 2 x ) a/(-\omega^2x) a / ( − ω 2 x ) で答えよ。
解答 (1) coth ( 1 / 2 ) = 2.164 \coth(1/2)=2.164 coth ( 1/2 ) = 2.164 なので分散は 2.164 σ g 2 2.164\,\sigma_{\mathrm g}^2 2.164 σ g 2 。γ / ω = tanh ( 1 / 2 ) = 0.462 \gamma/\omega=\tanh(1/2)=0.462 γ / ω = tanh ( 1/2 ) = 0.462 。
(2) β ℏ ω = ℏ ω / ( k B T ) = 0.480 \beta\hbar\omega=\hbar\omega/( k_{\mathrm B} T)=0.480 β ℏ ω = ℏ ω / ( k B T ) = 0.480 、coth ( 0.240 ) = 4.25 \coth(0.240)=4.25 coth ( 0.240 ) = 4.25 で、標準偏差は 4.25 σ g = 1.57 μ m \sqrt{4.25}\,\sigma_{\mathrm g}=1.57\,\mu\mathrm m 4.25 σ g = 1.57 μ m (σ g = 0.763 μ m \sigma_{\mathrm g}=0.763\,\mu\mathrm m σ g = 0.763 μ m )。γ / ω = tanh ( 0.240 ) = 0.235 \gamma/\omega=\tanh(0.240)=0.235 γ / ω = tanh ( 0.240 ) = 0.235 。
(3) a = − γ 2 x a=-\gamma^2x a = − γ 2 x なので比は ( γ / ω ) 2 (\gamma/\omega)^2 ( γ / ω ) 2 で、(1) では 0.214 0.214 0.214 。熱平衡の分布を1つの過程で再現すると、Newton の法則の力が ω 2 \omega^2 ω 2 から約5分の1に弱まった振動子の基底状態になる。T → 0 T\to0 T → 0 で tanh → 1 \tanh\to1 tanh → 1 となり、比は 1 1 1 に戻る。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com