tecktom.com

橋:基準・位相空間・解析力学との対応

2026年10月17日8分で読めますAI生成教科書『虎の巻』第7部 第1章 波動関数を密度と位相に分ける:Madelung 変換と量子ポテンシャル
目次
  1. 地図の上の位置
  2. 基準の言葉で:密度と流れは演算子の期待値
  3. 第5部の言葉で:Wigner 関数の運動量のモーメント
  4. 解析力学と授業の形
  5. 木下の層

AI 生成・木下が未勉強

量子力学の様々な構成方法 第7部 第1章 §7.1.6 ★★★ この章の結果が同値性の地図の「波動関数 ↔ 流体」の橋に当たること、基準と第5部の言葉で同じ量がどう書けるかを式で並べる

地図の上の位置

第1部第1章の同値性の地図(第1部第1章)では、第7部は「波動関数 ──Madelung 変換──▶ 流体」の矢印から始まる。この章で示したのはその矢印の両向きで、基準から流体へは第2節と第3節、流体から基準へは第4節の定理 7.1.4 が受け持つ。図 7.1.3 に、この章の結果がほかの部とつながる所をまとめた。

図

図 7.1.3 この章の結果と同値性の地図。太枠がこの章の流体。上下の2本の矢印がこの章で証明した両向きの橋で、下の3本は第7部の後の章、右の1本は第5部へ続く。

基準の言葉で:密度と流れは演算子の期待値

第1部の基準では、ρ\rho と j\bm{j} は演算子の期待値として書ける。1次元で書くと、δ(x^−x)=∣x⟩⟨x∣\delta(\hat x-x)=| x \rangle \langle x | を使って

ρ(x)=⟨ψ∣δ(x^−x)∣ψ⟩,j(x)=⟨ψ∣12m{p^, δ(x^−x)}∣ψ⟩(7.1.19)\rho(x)=\langle \psi | \delta(\hat x-x) | \psi \rangle ,\qquad j(x)=\langle \psi | \frac{1}{2m}\bigl\{\hat p,\ \delta(\hat x-x)\bigr\} | \psi \rangle\tag{7.1.19}

である({⋅,⋅}\{\cdot,\cdot\} は反交換子)。2つ目は、⟨ψ∣p^∣x⟩=(⟨x∣p^∣ψ⟩)∗=(−iℏψ′(x))∗\langle \psi |\hat p| x \rangle =\bigl(\langle x |\hat p| \psi \rangle \bigr)^*=(-\mathrm{i}\hbar\psi'(x))^* から12m[ψ∗(−iℏψ′)+ψ(iℏψ′∗)]=ℏmIm⁡(ψ∗ψ′)\frac{1}{2m}\bigl[\psi^*(-\mathrm{i}\hbar\psi')+\psi(\mathrm{i}\hbar\psi'^*)\bigr]=\frac{\hbar}{m}\operatorname{Im}(\psi^*\psi') となり、(7.1.3) と一致する。連続の式は、演算子 δ(x^−x)\delta(\hat x-x) の Heisenberg の運動方程式(第2部第3章)の期待値である。運動量の保存則 (7.1.11) を全空間で積分したものは ddt⟨p^⟩=−⟨V′(x^)⟩\frac{\mathrm{d} }{\mathrm{d} t}\langle \hat p \rangle =-\langle V'(\hat x) \rangle で、Heisenberg 方程式 p^˙=−V′(x^)\dot{\hat p}=-V'(\hat x) の期待値と同じ式になる。流体の側では「量子応力は表面の項なので積分すると消える」、基準の側では「[p^,p^2]=0[\hat p,\hat p^2]=0 なので運動エネルギーの項は p^\hat p を変えない」という別の理由で同じ式になっている。

運動量の分散も2つに分かれる。(7.1.12) を1次元で書き直し、⟨p^⟩=∫ρ mv dx\langle \hat p \rangle =\int\rho\,mv\,\mathrm{d} x を引くと

⟨p^2⟩−⟨p^⟩2=[∫ρ (mv)2dx−(∫ρ mv dx)2]+ℏ24J[ρ] ≥ ℏ24J[ρ] ≥ ℏ24 (⟨x^2⟩−⟨x^⟩2)(7.1.20)\langle \hat p^2 \rangle -\langle \hat p \rangle ^2=\Bigl[\int\rho\,(mv)^2\mathrm{d} x-\Bigl(\int\rho\,mv\,\mathrm{d} x\Bigr)^2\Bigr]+\frac{\hbar^2}{4}J[\rho] \ \ge\ \frac{\hbar^2}{4}J[\rho]\ \ge\ \frac{\hbar^2}{4\,(\langle \hat x^2 \rangle -\langle \hat x \rangle ^2)}\tag{7.1.20}

となる。最初の括弧は流れの速さのばらつき(分散なので 00 以上)、2つ目の不等号は位置の Fisher 情報についての Cramér–Rao の不等式 J[ρ] Var(x)≥1J[\rho]\,\mathrm{Var}(x)\ge1(相対エントロピーから情報と物理を組み立てる 第1部第7章)である。両端を結ぶと Heisenberg の不確定性関係 Δx Δp≥ℏ/2\Delta x\,\Delta p\ge\hbar/2 が出る。流体の言葉では、不確定性関係の下限は「流れが一様で、密度が Fisher 情報の下限を達成するガウス分布」のときに等号になる。

第5部の言葉で:Wigner 関数の運動量のモーメント

第5部の Wigner 関数 W(x,p)=1πℏ∫dy ψ(x+y)ψ∗(x−y)e−2ipy/ℏW(x,p)=\frac{1}{\pi\hbar}\int\mathrm{d} y\,\psi(x+y)\psi^*(x-y)\mathrm{e}^{-2\mathrm{i} py/\hbar} を運動量で積分すると、∫pkW dp\int p^kW\,\mathrm{d} p は ψ(x+y)ψ∗(x−y)\psi(x+y)\psi^*(x-y) を y=0y=0 で kk 回微分したものに (ℏ/(2i))k(\hbar/(2\mathrm{i}))^k を掛けた形になる。k=0,1,2k=0,1,2 を極形式で計算すると

∫W dp=ρ,∫p W dp=ρ mv,∫p2 W dp=ρ (mv)2+m P(7.1.21)\int W\,\mathrm{d} p=\rho,\qquad\int p\,W\,\mathrm{d} p=\rho\,mv,\qquad\int p^2\,W\,\mathrm{d} p=\rho\,(mv)^2+m\,P\tag{7.1.21}

となる(P=PxxP=P_{xx}。3本目は一般の R,SR,S について検算で確かめた)。つまり、流体の密度・運動量密度・運動量の流れは Wigner 関数の運動量についての0次・1次・2次のモーメントで、量子応力 PP は「その点での運動量の条件付き分散」に当たる。Moyal 方程式(第5部第3章)を pp で積分していくとモーメントの方程式が無限に続くが、純粋状態では高次のモーメントも ρ\rho と SS で決まるので、2次で閉じて Madelung の方程式になる。混合状態ではこの閉じ方が成り立たない。第5節の終わりで熱平衡を1つの流体で書けないと述べたのは、この事情による。

解析力学と授業の形

解析力学の巻の第1部第13章(解析力学から様々な理論分野を組み立てる 第1部第13章)は、同じ極形式から量子 HJ 方程式と連続の式を導き、ℏ→0\hbar\to0 の古典極限の条件を出した。その章は「古典力学の Hamilton–Jacobi 方程式が量子力学の位相の式の極限になる」向きに読み、この章は「量子力学を流体として書き直す構成の出発点」として読む。ℏ→0\hbar\to0 の極限そのものは次の章で扱う。

学部の量子力学の授業では、確率流密度 j=ℏ2mi(ψ∗∇ψ−ψ∇ψ∗)\bm{j}=\frac{\hbar}{2m\mathrm{i}}(\psi^*\nabla\psi-\psi\nabla\psi^*) と連続の式 ∂t∣ψ∣2+∇⋅j=0\partial_t\lvert \psi\rvert ^2+\nabla\cdot\bm{j}=0 を、Schrödinger 方程式に ψ∗\psi^* を掛けて複素共役と引き算する形で導く。この章の (7.1.6) はそれと同じ式で、授業の導き方は第2節の虚部だけを取り出したものに当たる。授業で通常省かれる実部の式が、この部の中心になる量子 HJ 方程式である。

例題 7.1.7 チャープしたガウス波束の Wigner 関数のモーメント ★★

オリジナルψ=(2πσ2)−1/4exp⁡(−x2/(4σ2)+icx2)\psi=(2\pi\sigma^2)^{-1/4}\exp\bigl(-x^2/(4\sigma^2)+\mathrm{i} cx^2\bigr)(cc は実数)の Wigner 関数を求め、(7.1.21) の3本を確かめよ。

解答

ψ(x+y)ψ∗(x−y)=(2πσ2)−1/2exp⁡(−(x+y)2+(x−y)24σ2+ic [(x+y)2−(x−y)2])=(2πσ2)−1/2e−x2/(2σ2) e−y2/(2σ2) e4icxy\psi(x+y)\psi^*(x-y)=(2\pi\sigma^2)^{-1/2}\exp\bigl(-\tfrac{(x+y)^2+(x-y)^2}{4\sigma^2}+\mathrm{i} c\,[(x+y)^2-(x-y)^2]\bigr) =(2\pi\sigma^2)^{-1/2}\mathrm{e}^{-x^2/(2\sigma^2)}\,\mathrm{e}^{-y^2/(2\sigma^2)}\,\mathrm{e}^{4\mathrm{i} cxy} である。yy についての Gauss 積分∫e−y2/(2σ2)e−2i(p−2ℏcx)y/ℏdy=2π σ e−2σ2(p−2ℏcx)2/ℏ2\int\mathrm{e}^{-y^2/(2\sigma^2)}\mathrm{e}^{-2\mathrm{i}(p-2\hbar cx)y/\hbar}\mathrm{d} y=\sqrt{2\pi}\,\sigma\,\mathrm{e}^{-2\sigma^2(p-2\hbar cx)^2/\hbar^2} を使うと

W(x,p)=1πℏ e−x2/(2σ2) e−2σ2(p−2ℏcx)2/ℏ2.W(x,p)=\frac{1}{\pi\hbar}\,\mathrm{e}^{-x^2/(2\sigma^2)}\,\mathrm{e}^{-2\sigma^2(p-2\hbar cx)^2/\hbar^2}.

各点 xx で、WW は pp について中心 2ℏcx2\hbar cx、分散 ℏ2/(4σ2)\hbar^2/(4\sigma^2) のガウス分布に ρ(x)\rho(x) を掛けたものになっている。したがって∫Wdp=ρ\int W\mathrm{d} p=\rho、∫pWdp=ρ⋅2ℏcx\int pW\mathrm{d} p=\rho\cdot2\hbar cx、∫p2Wdp=ρ [(2ℏcx)2+ℏ2/(4σ2)]\int p^2W\mathrm{d} p=\rho\,[(2\hbar cx)^2+\hbar^2/(4\sigma^2)] である。一方、流体の側では S=ℏcx2S=\hbar cx^2 から mv=2ℏcxmv=2\hbar cx、例題 7.1.3 の (1) から mP=ℏ2ρ/(4σ2)mP=\hbar^2\rho/(4\sigma^2) で、3本とも一致する。位相空間で見ると、チャープ cc は WW を p=2ℏcxp=2\hbar cx の直線に沿って傾けるだけで、その傾きが流体の速度場、直線の周りの pp の広がりが量子圧力である。

木下の勉強ステータス

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

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

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