tecktom.com

トレースの公式と、どの関数が状態か

2026年10月17日13分で読めますAI生成教科書『虎の巻』第5部 第2章 Wigner 変換と Weyl 対応
目次
  1. トレースの公式
  2. どの関数が状態の Wigner 関数か
  3. 木下の層

AI 生成・木下が未勉強

量子力学の様々な構成方法 第5部 第2章 §5.2.3 ★★★ Tr⁡(A^B^)=∫AWBW/(2πℏ)\operatorname{Tr}(\hat A\hat B)=\int A_WB_W/(2\pi\hbar)。これで期待値が位相空間の積分になり、密度行列との同値が両向きに閉じる

トレースの公式

定理 5.2.5(トレースの公式 ★★★)

2つの演算子 A^\hat A、B^\hat B について(A^B^\hat A\hat B のトレースが有限なとき)

Tr⁡(A^B^)=∫dx dp2πℏ AW(x,p) BW(x,p).(5.2.10)\operatorname{Tr}(\hat A\hat B)=\int\frac{\mathrm{d} x\,\mathrm{d} p}{2\pi\hbar}\,A_W(x,p)\,B_W(x,p) .\tag{5.2.10}

とくに B^=ρ^\hat B=\hat\rho と置くと、(5.2.3) から

⟨A^⟩=Tr⁡(ρ^A^)=∫W(x,p) AW(x,p) dx dp.(5.2.11)\langle \hat A \rangle =\operatorname{Tr}(\hat\rho\hat A)=\int W(x,p)\,A_W(x,p)\,\mathrm{d} x\,\mathrm{d} p .\tag{5.2.11}

証明

位置の完全系を2回挿入して Tr⁡(A^B^)=∫dx′ dx′′⟨x′∣A^∣x′′⟩⟨x′′∣B^∣x′⟩\operatorname{Tr}(\hat A\hat B)=\int\mathrm{d} x'\,\mathrm{d} x''\langle x' | \hat A | x'' \rangle \langle x'' | \hat B | x' \rangle。両方の因子に (5.2.6) を入れる。B^\hat B の因子は x′x' と x′′x'' を入れ替えた形なので位相の符号が逆になり、

Tr⁡(A^B^)=∫dx′ dx′′∫dp dp′(2πℏ)2 AW(X,p) BW(X,p′) ei(p−p′)(x′−x′′)/ℏ,X=x′+x′′2.\operatorname{Tr}(\hat A\hat B)=\int\mathrm{d} x'\,\mathrm{d} x''\int\frac{\mathrm{d} p\,\mathrm{d} p'}{(2\pi\hbar)^2}\,A_W(X,p)\,B_W(X,p')\,\mathrm{e}^{\mathrm{i}(p-p')(x'-x'')/\hbar},\qquad X=\frac{x'+x''}{2}.

積分変数を (x′,x′′)(x',x'') から (X,u)(X,u)、u=x′−x′′u=x'-x'' に替える。この変換のヤコビアンは 11 である。uu の積分は ∫du ei(p−p′)u/ℏ=2πℏ δ(p−p′)\int\mathrm{d} u\,\mathrm{e}^{\mathrm{i}(p-p')u/\hbar}=2\pi\hbar\,\delta(p-p') を与え、

Tr⁡(A^B^)=∫dX∫dp dp′2πℏ AW(X,p)BW(X,p′) δ(p−p′)=∫dX dp2πℏ AW(X,p) BW(X,p).\operatorname{Tr}(\hat A\hat B)=\int\mathrm{d} X\int\frac{\mathrm{d} p\,\mathrm{d} p'}{2\pi\hbar}\,A_W(X,p)B_W(X,p')\,\delta(p-p')=\int\frac{\mathrm{d} X\,\mathrm{d} p}{2\pi\hbar}\,A_W(X,p)\,B_W(X,p) .

この定理が位相空間の構成の中心の式である。期待値は「状態の関数 WW と物理量の関数 AWA_W の積を相空間で積分する」という、古典の統計力学と同じ形になった。古典と違うのは2点で、WW が負になりうること(前の章)と、AWA_W が古典の関数と一致するとは限らないこと(例題 5.2.1)である。2つ目の点を、正の演算子の期待値が負になる例で確かめる。

例題 5.2.3 基底状態での 12(x^2p^2+p^2x^2)\tfrac12(\hat x^2\hat p^2+\hat p^2\hat x^2) ★★★

オリジナル調和振動子の基底状態 ∣0⟩| 0 \rangle(W0=1πℏe−2H/(ℏω)W_0=\frac{1}{\pi\hbar}\mathrm{e}^{-2H/(\hbar\omega)})で、演算子 C^=12(x^2p^2+p^2x^2)\hat C=\tfrac12(\hat x^2\hat p^2+\hat p^2\hat x^2) を考える。

  1. (1) C^\hat C の Weyl 記号が CW=x2p2−ℏ2/2C_W=x^2p^2-\hbar^2/2 であることを、例題 5.2.5 の結果 (x^2p^2)W=x2p2+2iℏxp−ℏ2/2(\hat x^2\hat p^2)_W=x^2p^2+2\mathrm{i}\hbar xp-\hbar^2/2 から導け。
  2. (2) (5.2.11) で ⟨C^⟩\langle \hat C \rangle を計算せよ。
  3. (3) 古典の関数 x2p2x^2p^2 をそのまま W0W_0 で平均したもの、および ⟨x^2⟩⟨p^2⟩\langle \hat x^2 \rangle \langle \hat p^2 \rangle と比べよ。
解答

(1) p^2x^2=(x^2p^2)†\hat p^2\hat x^2=(\hat x^2\hat p^2)^\dagger なので、命題 5.2.4 から (p^2x^2)W=x2p2−2iℏxp−ℏ2/2(\hat p^2\hat x^2)_W=x^2p^2-2\mathrm{i}\hbar xp-\hbar^2/2。平均すると虚部が消えて CW=x2p2−ℏ2/2C_W=x^2p^2-\hbar^2/2。

(2) W0W_0 は xx と pp について独立な Gauss 関数で、∫x2W0=ℏ/(2mω)\int x^2W_0=\hbar/(2m\omega)、∫p2W0=mℏω/2\int p^2W_0=m\hbar\omega/2。よって

⟨C^⟩=ℏ2mω⋅mℏω2−ℏ22=ℏ24−ℏ22=−ℏ24.\langle \hat C \rangle =\frac{\hbar}{2m\omega}\cdot\frac{m\hbar\omega}{2}-\frac{\hbar^2}{2}=\frac{\hbar^2}{4}-\frac{\hbar^2}{2}=-\frac{\hbar^2}{4}.

位置表示で C^ψ0\hat C\psi_0 を微分して直接計算しても同じ値になる(検算で確かめた)。

(3) 古典の関数の平均は ∫x2p2W0=ℏ2/4\int x^2p^2W_0=\hbar^2/4 で、⟨x^2⟩⟨p^2⟩\langle \hat x^2 \rangle \langle \hat p^2 \rangle と同じく正である。x^2\hat x^2 と p^2\hat p^2 はどちらも正の演算子だが、交換しないのでその対称化した積は正の演算子でなく、基底状態での期待値は負になる。その差 −ℏ2/2-\hbar^2/2 がちょうど記号の ℏ2\hbar^2 の項である。位相空間で期待値を計算するときは、物理量を古典の関数のまま使ってよいのではなく、その演算子の Weyl 記号を使わなければならない。

どの関数が状態の Wigner 関数か

トレースの公式を2つの状態に当てると

Tr⁡(ρ^1ρ^2)=2πℏ∫W1W2 dx dp,Tr⁡ρ^2=2πℏ∫W2 dx dp(5.2.12)\operatorname{Tr}(\hat\rho_1\hat\rho_2)=2\pi\hbar\int W_1W_2\,\mathrm{d} x\,\mathrm{d} p,\qquad \operatorname{Tr}\hat\rho^2=2\pi\hbar\int W^2\,\mathrm{d} x\,\mathrm{d} p\tag{5.2.12}

となる。純粋状態 ρ^2=∣φ⟩⟨φ∣\hat\rho_2=| \varphi \rangle \langle \varphi | なら左辺は ⟨φ∣ρ^1∣φ⟩≥0\langle \varphi | \hat\rho_1 | \varphi \rangle \ge0 で、前の章の「純粋状態で ∫W2=1/(2πℏ)\int W^2=1/(2\pi\hbar)」は2つ目の式で Tr⁡ρ^2=1\operatorname{Tr}\hat\rho^2=1 と置いたものである。これで、位相空間の関数のうちどれが量子状態を表すかを言い切ることができる。

定理 5.2.6(密度行列と Wigner 関数の1対1対応 ★★★)

Wigner 変換 ρ^↦W\hat\rho\mapsto W は、密度行列の全体から次の3条件を満たす関数 W(x,p)W(x,p) の全体への1対1の写像で、逆は (5.2.8) で与えられる。

  1. (i) WW は実数値。
  2. (ii) ∫W dx dp=1\int W\,\mathrm{d} x\,\mathrm{d} p=1。
  3. (iii) どの純粋状態 φ\varphi についても 2πℏ∫W Wφ dx dp≥02\pi\hbar\int W\,W_\varphi\,\mathrm{d} x\,\mathrm{d} p\ge0(WφW_\varphi は φ\varphi の Wigner 関数)。
証明

(密度行列 → 関数)ρ^\hat\rho はエルミートなので命題 5.2.4 から (i)、Tr⁡ρ^=1\operatorname{Tr}\hat\rho=1 と (5.2.9) から (ii)、ρ^≥0\hat\rho\ge0 と (5.2.12) から 2πℏ∫WWφ=⟨φ∣ρ^∣φ⟩≥02\pi\hbar\int WW_\varphi=\langle \varphi | \hat\rho | \varphi \rangle \ge0 で (iii) が成り立つ。異なる ρ^\hat\rho は異なる位置表示の行列要素を持ち、Fourier 変換は1対1なので、異なる WW に写る。

(関数 → 密度行列)(i)〜(iii) を満たす WW に (5.2.8) で演算子 ρ^=QW[2πℏW]\hat\rho=\mathcal Q_{\mathrm W}[2\pi\hbar W] を作る。(5.2.7) からその Wigner 関数は WW に戻る。(i) と命題 5.2.4 から ρ^\hat\rho はエルミート、(ii) と (5.2.9) から Tr⁡ρ^=1\operatorname{Tr}\hat\rho=1、(iii) と (5.2.12) からすべての φ\varphi で ⟨φ∣ρ^∣φ⟩≥0\langle \varphi | \hat\rho | \varphi \rangle \ge0、つまり ρ^≥0\hat\rho\ge0 である。よって ρ^\hat\rho は密度行列である。

条件 (iii) が、古典の「至る所で非負」の代わりに入る条件である。WW の値そのものは負でもよいが、どの純粋状態の WφW_\varphi と重ねて積分しても負にならない。古典の分布の非負性が各点ごとの条件だったのに対し、量子の条件は全体にまたがる条件になっている。この条件のため、勝手な正の Gauss 関数が状態になるとは限らない。

例題 5.2.4 細すぎる Gauss 関数は状態でない ★★★

オリジナルWg(x,p)=12πσxσpexp⁡(−x22σx2−p22σp2)W_g(x,p)=\frac{1}{2\pi\sigma_x\sigma_p}\exp\bigl(-\frac{x^2}{2\sigma_x^2}-\frac{p^2}{2\sigma_p^2}\bigr) を考える。至る所で正で、積分は 11 である。

  1. (1) 2πℏ∫Wg2 dx dp2\pi\hbar\int W_g^2\,\mathrm{d} x\,\mathrm{d} p を求め、σxσp<ℏ/2\sigma_x\sigma_p<\hbar/2 なら WgW_g はどの密度行列の Wigner 関数でもないことを示せ。
  2. (2) σxσp=ℏ/4\sigma_x\sigma_p=\hbar/4 のとき、条件 (iii) を破る純粋状態 φ\varphi を1つ挙げ、「⟨φ∣ρ^∣φ⟩\langle \varphi | \hat\rho | \varphi \rangle」に当たる量を計算せよ。
解答

(1) ∫Wg2=1(2πσxσp)2⋅πσx⋅πσp=14πσxσp\int W_g^2=\frac{1}{(2\pi\sigma_x\sigma_p)^2}\cdot\sqrt{\pi}\sigma_x\cdot\sqrt\pi\sigma_p=\frac{1}{4\pi\sigma_x\sigma_p} なので

2πℏ∫Wg2 dx dp=ℏ2σxσp.2\pi\hbar\int W_g^2\,\mathrm{d} x\,\mathrm{d} p=\frac{\hbar}{2\sigma_x\sigma_p}.

WgW_g が密度行列 ρ^\hat\rho の Wigner 関数なら、これは Tr⁡ρ^2\operatorname{Tr}\hat\rho^2 で、11 以下でなければならない(第4部第1章の純粋度)。σxσp<ℏ/2\sigma_x\sigma_p<\hbar/2 なら 11 を超えるので矛盾する。等号 σxσp=ℏ/2\sigma_x\sigma_p=\hbar/2 で純粋度がちょうど 11 になり、そのとき WgW_g は Gauss 波束(調和振動子の基底状態を伸縮したもの)の WW である。不確定性関係の下限が、位相空間では「Gauss 関数の幅の積の下限」として現れる。

(2) 同じ σx\sigma_x を持つ最小の Gauss 波束 φ\varphi、つまり Wφ=1πℏexp⁡(−x22σx2−2σx2p2ℏ2)W_\varphi=\frac{1}{\pi\hbar}\exp\bigl(-\frac{x^2}{2\sigma_x^2}-\frac{2\sigma_x^2p^2}{\hbar^2}\bigr) を取る。σp=ℏ/(4σx)\sigma_p=\hbar/(4\sigma_x) では WgW_g の指数の p2p^2 の係数は 1/(2σp2)=8σx2/ℏ21/(2\sigma_p^2)=8\sigma_x^2/\hbar^2 である。積 WgWφW_gW_\varphi の指数は、x2x^2 の係数が1/(2σx2)+1/(2σx2)=1/σx21/(2\sigma_x^2)+1/(2\sigma_x^2)=1/\sigma_x^2、p2p^2 の係数が (8+2)σx2/ℏ2=10σx2/ℏ2(8+2)\sigma_x^2/\hbar^2=10\sigma_x^2/\hbar^2 で、

∫e−x2/σx2dx=π σx,∫e−10σx2p2/ℏ2dp=π ℏ10 σx.\int\mathrm{e}^{-x^2/\sigma_x^2}\mathrm{d} x=\sqrt\pi\,\sigma_x,\qquad \int\mathrm{e}^{-10\sigma_x^2p^2/\hbar^2}\mathrm{d} p=\frac{\sqrt\pi\,\hbar}{\sqrt{10}\,\sigma_x}.

前の係数は 2πℏ⋅12πσxσp⋅1πℏ=1πσxσp=4πℏ2\pi\hbar\cdot\frac{1}{2\pi\sigma_x\sigma_p}\cdot\frac{1}{\pi\hbar}=\frac{1}{\pi\sigma_x\sigma_p}=\frac{4}{\pi\hbar} なので、全体は

2πℏ∫WgWφ dx dp=4πℏ⋅πℏ10=410=1.2652\pi\hbar\int W_gW_\varphi\,\mathrm{d} x\,\mathrm{d} p=\frac{4}{\pi\hbar}\cdot\frac{\pi\hbar}{\sqrt{10}}=\frac{4}{\sqrt{10}}=1.265

である。

これは「確率」⟨φ∣ρ^∣φ⟩\langle \varphi | \hat\rho | \varphi \rangle に当たる量なのに 11 を超えている。ρ^\hat\rho が密度行列なら ⟨φ∣ρ^∣φ⟩≤Tr⁡ρ^=1\langle \varphi | \hat\rho | \varphi \rangle \le\operatorname{Tr}\hat\rho=1 なので(ρ^≥0\hat\rho\ge0 のとき、1つの規格化ベクトルでの値は固有値の和を超えない)、WgW_g は状態でない。QW[2πℏWg]\mathcal Q_{\mathrm W}[2\pi\hbar W_g] は、固有値の和が 11 で最大の固有値が 11 を超える演算子、つまり負の固有値を持つ演算子である。

木下の勉強ステータス

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

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

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