tecktom.com

この章のまとめ

2026年10月17日10分で読めますAI生成教科書『虎の巻』第5部 第2章 Wigner 変換と Weyl 対応

AI 生成・木下が未勉強

量子力学の様々な構成方法 第5部 第2章 §5.2.7 ★★★ Wigner 変換は1対1で、期待値は ∫WAW\int WA_W。密度行列と位相空間の構成は同値である

演算子の位置表示を中点と差の座標に書き直し、差について Fourier 変換する Wigner 変換を、状態と物理量の両方に当てた。状態の変換が Wigner 関数 W=ρW/(2πℏ)W=\rho_W/(2\pi\hbar)、物理量の変換が Weyl 記号 AWA_W で、逆の写像(Weyl 量子化 QW\mathcal Q_{\mathrm W})は Fourier 逆変換として書くことができる。したがって WW から ρ^\hat\rho は元どおりに復元でき、エルミートは実数値に、トレースは ∫dx dp/(2πℏ)\int\mathrm{d} x\,\mathrm{d} p/(2\pi\hbar) に写る。トレースの公式 Tr⁡(A^B^)=∫AWBW/(2πℏ)\operatorname{Tr}(\hat A\hat B)=\int A_WB_W/(2\pi\hbar) から期待値は ∫WAW\int WA_W になり、状態を表す関数は「実数値・規格化・どの純粋状態の WφW_\varphi と重ねても非負」の3条件で特徴づけられる。これで密度行列と位相空間の構成の同値が両向きに示された。

Weyl 記号は古典の関数と一致するとは限らず、x^p^\hat x\hat p の記号は xp+iℏ/2xp+\mathrm{i}\hbar/2、12(x^2p^2+p^2x^2)\tfrac12(\hat x^2\hat p^2+\hat p^2\hat x^2) の記号は x2p2−ℏ2/2x^2p^2-\hbar^2/2 である。指数関数 ei(ax+bp)/ℏ\mathrm{e}^{\mathrm{i}(ax+bp)/\hbar} が ei(ax^+bp^)/ℏ\mathrm{e}^{\mathrm{i}(a\hat x+b\hat p)/\hbar} に写ることから、Weyl 対応は x^\hat x と p^\hat p の並べ方を平均する Weyl 順序であり、WW は量子の特性関数の Fourier 逆変換で、1次結合 ax^+bp^a\hat x+b\hat p の分布は WW の直線に沿った積分になる。ハミルトニアンが2次以下なら WW は古典の流れで運ばれ、共通例題③の波束の広がり σ(t)2=σ02+(ℏt/(2mσ0))2\sigma(t)^2=\sigma_0^2+(\hbar t/(2m\sigma_0))^2 は WW の剪断として出る。積が何に写るか(⋆\star 積)と一般のハミルトニアンでの時間発展は 第5部第3章で扱う。

総合問題 Schrödinger の猫型の重ね合わせを位相空間で読む ★★★

オリジナルφ(x)=(2πσ2)−1/4e−x2/(4σ2)\varphi(x)=(2\pi\sigma^2)^{-1/4}\mathrm{e}^{-x^2/(4\sigma^2)}(∣φ∣2\lvert \varphi\rvert ^2 の分散 σ2\sigma^2)とし、a>0a>0 として、2か所に置いた波束の重ね合わせ

ψ±(x)=N±[φ(x−a)±φ(x+a)],N±2=12(1±ϵ),ϵ=e−a2/(2σ2)\psi_\pm(x)=N_\pm\bigl[\varphi(x-a)\pm\varphi(x+a)\bigr],\qquad N_\pm^2=\frac{1}{2(1\pm\epsilon)},\quad \epsilon=\mathrm{e}^{-a^2/(2\sigma^2)}

と、混合状態 ρ^mix=12(∣φ+⟩⟨φ+∣+∣φ−⟩⟨φ−∣)\hat\rho_{\mathrm{mix}}=\tfrac12\bigl(| \varphi_+ \rangle \langle \varphi_+ |+| \varphi_- \rangle \langle \varphi_- |\bigr)(φ±(x)=φ(x∓a)\varphi_\pm(x)=\varphi(x\mp a))を比べる。WG(x,p)=1πℏexp⁡(−x22σ2−2σ2p2ℏ2)W_G(x,p)=\frac{1}{\pi\hbar}\exp\bigl(-\frac{x^2}{2\sigma^2}-\frac{2\sigma^2p^2}{\hbar^2}\bigr) を原点の波束 φ\varphi の Wigner 関数とする。

  1. (1) ψ+\psi_+ の Wigner 関数が W+=N+2[WG(x−a,p)+WG(x+a,p)+2WG(x,p)cos⁡(2pa/ℏ)]W_+=N_+^2\bigl[W_G(x-a,p)+W_G(x+a,p)+2W_G(x,p)\cos(2pa/\hbar)\bigr] であることを示せ。(使う道具:この章 5.2.1 節)
  2. (2) W+(0,0)W_+(0,0) と、ψ−\psi_- の Wigner 関数の原点の値を求め、パリティの期待値と照合せよ。(使う道具:第5部第1章の W(0,0)=⟨ψ∣Π^∣ψ⟩/(πℏ)W(0,0)=\langle \psi | \hat\Pi | \psi \rangle /(\pi\hbar))
  3. (3) ψ+\psi_+ と ρ^mix\hat\rho_{\mathrm{mix}} の純粋度を 2πℏ∫W22\pi\hbar\int W^2 で求めよ。a=2σa=2\sigma での値を答えよ。(使う道具:この章 5.2.3 節・第4部第1章の純粋度)
  4. (4) ψ+\psi_+ と ρ^mix\hat\rho_{\mathrm{mix}} の ⟨p^2⟩\langle \hat p^2 \rangle を ∫p2W\int p^2W で求め、a=2σa=2\sigma で比べよ。位置の分布しか測らない実験で両者を区別できるか答えよ。(使う道具:この章 5.2.3 節・式 (5.2.5))
  5. (5) 自由粒子として時間発展させたとき、ψ+\psi_+ の ⟨x^2⟩t\langle \hat x^2 \rangle _t を求めよ。ℏ=m=σ=1\hbar=m=\sigma=1、a=2a=2、t=3t=3 での値を答えよ。(使う道具:この章 5.2.5 節・第2部第3章の Heisenberg 方程式)
解答

(1) Wigner 変換は ψ\psi と ψ∗\psi^* について双線形なので、N+2N_+^2 を除くと4つの項φμ(x+y)φν∗(x−y)\varphi_\mu(x+y)\varphi_\nu^*(x-y)(μ,ν=±\mu,\nu=\pm)の変換の和になる。対角の2項は平行移動の性質(前の章の (W4))から WG(x∓a,p)W_G(x\mp a,p)。交差の項 φ+(x+y)φ−(x−y)=φ(x+y−a)φ(x−y+a)\varphi_+(x+y)\varphi_-(x-y)=\varphi(x+y-a)\varphi(x-y+a) の指数は −[(x+y−a)2+(x−y+a)2]/(4σ2)=−[x2+(y−a)2]/(2σ2)-\bigl[(x+y-a)^2+(x-y+a)^2\bigr]/(4\sigma^2)=-\bigl[x^2+(y-a)^2\bigr]/(2\sigma^2) で、yy を y+ay+a にずらすと原点の波束の積分に e−2ipa/ℏ\mathrm{e}^{-2\mathrm{i} pa/\hbar} が掛かった形になり、WG(x,p)e−2ipa/ℏW_G(x,p)\mathrm{e}^{-2\mathrm{i} pa/\hbar} を与える。もう一方の交差の項はその複素共役で、足して 2WG(x,p)cos⁡(2pa/ℏ)2W_G(x,p)\cos(2pa/\hbar)。2つの波束の真ん中 x=0x=0 に、pp 方向の周期 πℏ/a\pi\hbar/a の縞が現れる。これが重ね合わせの干渉の位相空間での姿である。

(2) W+(0,0)=N+2[2ϵ+2]/(πℏ)=1/(πℏ)W_+(0,0)=N_+^2\bigl[2\epsilon+2\bigr]/(\pi\hbar)=1/(\pi\hbar)(WG(±a,0)=ϵ/(πℏ)W_G(\pm a,0)=\epsilon/(\pi\hbar))。ψ+\psi_+ は偶関数でパリティの期待値は 11 なので一致する。ψ−\psi_- では交差の項の符号が逆になり W−(0,0)=N−2[2ϵ−2]/(πℏ)=−1/(πℏ)W_-(0,0)=N_-^2[2\epsilon-2]/(\pi\hbar)=-1/(\pi\hbar)。ψ−\psi_- は奇関数でパリティ −1-1 と一致する。ψ−\psi_- は至る所で正の波束2つの重ね合わせなのに、その WW は原点で下限 −1/(πℏ)-1/(\pi\hbar) に達する。

(3) ψ+\psi_+ は純粋状態なので 2πℏ∫W+2=12\pi\hbar\int W_+^2=1 である(検算では数値積分で確かめた)。混合状態は Wmix=12[WG(x−a,p)+WG(x+a,p)]W_{\mathrm{mix}}=\tfrac12[W_G(x-a,p)+W_G(x+a,p)] で、(5.2.12) から 2πℏ∫WG(x−a,p)WG(x+a,p)=∣⟨φ+∣φ−⟩∣2=ϵ22\pi\hbar\int W_G(x-a,p)W_G(x+a,p)=\lvert \langle \varphi_+ | \varphi_- \rangle \rvert ^2=\epsilon^2。したがって

2πℏ∫Wmix2=14(1+1+2ϵ2)=1+e−a2/σ22.2\pi\hbar\int W_{\mathrm{mix}}^2=\frac14\bigl(1+1+2\epsilon^2\bigr)=\frac{1+\mathrm{e}^{-a^2/\sigma^2}}{2}.

a=2σa=2\sigma では (1+e−4)/2=0.5092(1+\mathrm{e}^{-4})/2=0.5092。縞の項を消すと純粋度がほぼ半分に下がる。縞が「純粋な重ね合わせ」であるための情報を持っている。

(4) 運動量の分布は WW を xx で積分したもので、WGW_G の pp の周辺分布は分散 v=ℏ2/(4σ2)v=\hbar^2/(4\sigma^2) の正規分布である。∫p2cos⁡(kp) N(0,v)(p) dp=(v−k2v2)e−k2v/2\int p^2\cos(kp)\,\mathcal N(0,v)(p)\,\mathrm{d} p=(v-k^2v^2)\mathrm{e}^{-k^2v/2} に k=2a/ℏk=2a/\hbar(k2v=a2/σ2k^2v=a^2/\sigma^2)を入れると、交差の項からの寄与は 2v(1−a2/σ2)ϵ2v(1-a^2/\sigma^2)\epsilon で、

⟨p^2⟩+=v 1+(1−a2/σ2)ϵ1+ϵ,⟨p^2⟩mix=v.\langle \hat p^2 \rangle _+=v\,\frac{1+(1-a^2/\sigma^2)\epsilon}{1+\epsilon},\qquad \langle \hat p^2 \rangle _{\mathrm{mix}}=v .

a=2σa=2\sigma では ϵ=e−2=0.1353\epsilon=\mathrm{e}^{-2}=0.1353 で ⟨p^2⟩+=v(1−3×0.1353)/1.1353=0.5232 v\langle \hat p^2 \rangle _+=v(1-3\times0.1353)/1.1353=0.5232\,v。位置の分布は WW を pp で積分したもので、縞の項は ∫cos⁡(2pa/ℏ)WG dp=ϵ ∣φ(x)∣2\int\cos(2pa/\hbar)W_G\,\mathrm{d} p=\epsilon\,\lvert \varphi(x)\rvert ^2 を与える。a=2σa=2\sigma ではこれは小さいが 00 でなく、原点付近の位置の分布にも差が出る。a≫σa\gg\sigma では ϵ→0\epsilon\to0 で位置の分布の差は消え、2つを位置だけで区別することはできない。運動量の分布(縞の向きの射影)には、2つの波束の運動量表示を足して2乗した干渉の縞2N+2∣φ~(p)∣2(1+cos⁡(2pa/ℏ))2N_+^2\lvert \tilde\varphi(p)\rvert ^2\bigl(1+\cos(2pa/\hbar)\bigr) として違いが残る(φ~\tilde\varphi は φ\varphi の運動量表示)。これは前の章の定理(非負で双線形な分布は無い)の証明に使った干渉項である。

(5) 自由粒子では (5.2.20) から xt=x0+p0t/mx_t=x_0+p_0t/m の分布で、⟨x^2⟩t=⟨x2⟩0+2(t/m)⟨xp⟩0+(t/m)2⟨p2⟩0\langle \hat x^2 \rangle _t=\langle x^2 \rangle _0+2(t/m)\langle xp \rangle _0+(t/m)^2\langle p^2 \rangle _0。ψ+\psi_+ は実数値なので ∫pW dp=ℏIm⁡(ψ∗ψ′)=0\int pW\,\mathrm{d} p=\hbar\operatorname{Im}(\psi^*\psi')=0(前の章の確率流の式)で、⟨xp⟩0=0\langle xp \rangle _0=0。位置の2次モーメントは⟨x2⟩0=2N+2[(σ2+a2)+σ2ϵ]=σ2+a21+ϵ\langle x^2 \rangle _0=2N_+^2\bigl[(\sigma^2+a^2)+\sigma^2\epsilon\bigr]=\sigma^2+\dfrac{a^2}{1+\epsilon} なので

⟨x^2⟩t=σ2+a21+ϵ+t2m2 ℏ24σ2 1+(1−a2/σ2)ϵ1+ϵ.\langle \hat x^2 \rangle _t=\sigma^2+\frac{a^2}{1+\epsilon}+\frac{t^2}{m^2}\,\frac{\hbar^2}{4\sigma^2}\,\frac{1+(1-a^2/\sigma^2)\epsilon}{1+\epsilon}.

ℏ=m=σ=1\hbar=m=\sigma=1、a=2a=2、t=3t=3 では ⟨x2⟩0=1+4/1.1353=4.523\langle x^2 \rangle _0=1+4/1.1353=4.523、⟨p2⟩0=0.25×0.5232=0.1308\langle p^2 \rangle _0=0.25\times0.5232=0.1308 で、⟨x^2⟩t=4.523+9×0.1308=5.700\langle \hat x^2 \rangle _t=4.523+9\times0.1308=5.700。検算では、波動関数を Fourier 変換で時間発展させて直接計算した値と一致することを確かめた。

木下の勉強ステータス

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

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

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