tecktom.com

1つのモードの位相空間と、並べ方で分かれる3つの分布

2026年10月17日18分で読めますAI生成教科書『虎の巻』第5部 第6章 光の状態の地図:コヒーレント・スクイーズド・猫状態と P・Q 表示、位相空間の構成の限界
目次
  1. 複素の座標
  2. 特性関数の族
  3. 並べ方の定理と、Gauss 関数の畳み込み
  4. 木下の層

AI 生成・木下が未勉強

量子力学の様々な構成方法 第5部 第6章 §5.6.1 ★★★ 特性関数に es∣ξ∣2/2\mathrm{e}^{s\lvert \xi\rvert ^2/2} を掛けた族 FsF_s で PP(s=1s=1)・WW(s=0s=0)・QQ(s=−1s=-1)を1度に定義し、ss 順序の期待値を与えることと、ss を下げることが Gauss 関数での畳み込みであることを示す

第5部はここまで、1粒子の位置 xx と運動量 pp の相空間の上で量子力学を組み立ててきた。第5部第1章で Wigner 関数 WW を定義し、第5部第2章で密度行列と1対1に対応することを、第5部第3章で演算子の積が ⋆\star 積に写ることを示した。第5部第4章は WW の負の値の意味を調べ、第5部第5章は共通例題①を位相空間の道具だけで解いた。

この章は部の型(出発点の要請 → 基本の対象 → 基準との同値 → 共通例題 → 得意な問題 → 限界)の最後の2つ、「得意な問題」と「限界」を受け持つ。位相空間の構成が最も使われている分野は量子光学で、そこでは1つの光のモードの状態を、WW のほかに Glauber–Sudarshan の PP 表示と Husimi の QQ 関数という2つの分布で書く。3つの分布は、第二量子化(第6部)の生成消滅演算子 a^,a^†\hat a,\hat a^\dagger をどの順に並べた期待値を与えるかで分かれる。PP は a^†\hat a^\dagger を左に集めた正規順序、WW は並べ方を平均した対称順序、QQ は a^\hat a を左に集めた反正規順序の期待値を、古典の確率分布と同じ積分で与える。この対応が計画の地図の「第二量子化 ↔ 位相空間」の橋で、5.6.6節で両向きに示す。

出発点に置くのは前の章までの結果と、第2部の生成消滅演算子 a^\hat a、[a^, a^†]=1\left[\hat a,\,\hat a^\dagger\right]=1 である。5.6.1節で3つの分布を1つの族として定義し、互いに Gauss 関数の畳み込みでつながることを示す。5.6.2節で共通例題⑤(熱的な光)を3つの分布で解き、5.6.3節でコヒーレント状態とスクイーズド状態、5.6.4節で猫状態と損失による縞の消え方を扱う。5.6.5節は、実験で測られる分布がどれかと、位相空間の構成の限界をまとめる。1つのモードの x,px,p は粒子の位置と運動量ではなく、光の電場の2つの直交位相成分(余弦の振幅と正弦の振幅)である。式は第5部の1粒子の式と同じ形をしているので、前の章までの結果はそのまま使うことができる。

複素の座標 α\alpha

角周波数 ω\omega のモード(調和振動子と同じ H^=ℏω(a^†a^+12)\hat H=\hbar\omega(\hat a^\dagger\hat a+\frac12))の相空間に、複素の座標

α=12ℏ(mω x+ipmω)=X1+iX2,d2α≔dX1 dX2=dx dp2ℏ(5.6.1)\alpha=\frac{1}{\sqrt{2\hbar}}\Bigl(\sqrt{m\omega}\,x+\frac{\mathrm{i} p}{\sqrt{m\omega}}\Bigr)=X_1+\mathrm{i} X_2,\qquad \mathrm{d}^2\alpha\coloneqq\mathrm{d} X_1\,\mathrm{d} X_2=\frac{\mathrm{d} x\,\mathrm{d} p}{2\hbar}\tag{5.6.1}

を置く。α\alpha は消滅演算子 a^\hat a を古典の数にしたもので(第6部第1章の古典の複素振幅と同じ)、第5部第4章がコヒーレント状態の中心として使ったものと同じ規格化である。実部と虚部 X1,X2X_1,X_2 は、直交位相の演算子 X^θ=(a^e−iθ+a^†eiθ)/2\hat X_\theta=(\hat a\mathrm{e}^{-\mathrm{i}\theta}+\hat a^\dagger\mathrm{e}^{\mathrm{i}\theta})/2 の θ=0,π/2\theta=0,\pi/2 の値に当たる(光学と量子エレクトロニクス 第31章と同じ 12\frac12 の規格化で、真空の分散は 14\frac14)。

この章の分布は、すべて d2α\mathrm{d}^2\alpha について規格化する。dx dp\mathrm{d} x\,\mathrm{d} p について規格化した前の章の W(x,p)W(x,p) とは

W(α)=2ℏ W(x,p)(5.6.2)W(\alpha)=2\hbar\,W(x,p)\tag{5.6.2}

の関係にある。W(α)W(\alpha) と W(x,p)W(x,p) を引数の字で書き分ける。例えば基底状態(真空)は W(x,p)=1πℏe−2H/(ℏω)W(x,p)=\frac{1}{\pi\hbar}\mathrm{e}^{-2H/(\hbar\omega)}、H/(ℏω)=∣α∣2H/(\hbar\omega)=\lvert \alpha\rvert ^2 からW(α)=2πe−2∣α∣2W(\alpha)=\frac2\pi\mathrm{e}^{-2\lvert \alpha\rvert ^2} である。

特性関数の族

変位演算子 D^(ξ)=exp⁡(ξa^†−ξ∗a^)\hat D(\xi)=\exp(\xi\hat a^\dagger-\xi^*\hat a)(光学と量子エレクトロニクス 第31章の D^\hat D)を使う。ξa^†\xi\hat a^\dagger と −ξ∗a^-\xi^*\hat a の交換子は[ξa^†, −ξ∗a^]=∣ξ∣2\left[\xi\hat a^\dagger,\,-\xi^*\hat a\right]=\lvert \xi\rvert ^2 で数なので、Baker–Campbell–Hausdorff の公式 eA^+B^=eA^eB^e−[A^, B^]/2\mathrm{e}^{\hat A+\hat B}=\mathrm{e}^{\hat A}\mathrm{e}^{\hat B}\mathrm{e}^{-\left[\hat A,\,\hat B\right]/2} が厳密に成り立ち、

D^(ξ)=e−∣ξ∣2/2 eξa^†e−ξ∗a^=e+∣ξ∣2/2 e−ξ∗a^eξa^†(5.6.3)\hat D(\xi)=\mathrm{e}^{-\lvert \xi\rvert ^2/2}\,\mathrm{e}^{\xi\hat a^\dagger}\mathrm{e}^{-\xi^*\hat a}=\mathrm{e}^{+\lvert \xi\rvert ^2/2}\,\mathrm{e}^{-\xi^*\hat a}\mathrm{e}^{\xi\hat a^\dagger}\tag{5.6.3}

となる。中央の形では a^†\hat a^\dagger が左(正規順序)、右の形では a^\hat a が左(反正規順序)に並んでいる。並べ方の違いは、Gauss 型の因子 e∓∣ξ∣2/2\mathrm{e}^{\mp\lvert \xi\rvert ^2/2} だけに集まる。

定義 5.6.1(ss 順序の特性関数と分布)

実数 ss について、状態 ρ^\hat\rho の ss 順序の特性関数と分布を

χs(ξ)=Tr⁡(ρ^D^(ξ))es∣ξ∣2/2,Fs(α)=1π2∫χs(ξ) eαξ∗−α∗ξ d2ξ(5.6.4)\chi_s(\xi)=\operatorname{Tr}\bigl(\hat\rho\hat D(\xi)\bigr)\mathrm{e}^{s\lvert \xi\rvert ^2/2},\qquad F_s(\alpha)=\frac{1}{\pi^2}\int\chi_s(\xi)\,\mathrm{e}^{\alpha\xi^*-\alpha^*\xi}\,\mathrm{d}^2\xi\tag{5.6.4}

で定める。s=1s=1 を Glauber–Sudarshan の PP 表示、s=0s=0 を Wigner 関数 WW、s=−1s=-1 を Husimi の QQ 関数と呼ぶ。

(5.6.3) から χ1(ξ)=Tr⁡(ρ^ eξa^†e−ξ∗a^)\chi_1(\xi)=\operatorname{Tr}(\hat\rho\,\mathrm{e}^{\xi\hat a^\dagger}\mathrm{e}^{-\xi^*\hat a})、χ−1(ξ)=Tr⁡(ρ^ e−ξ∗a^eξa^†)\chi_{-1}(\xi)=\operatorname{Tr}(\hat\rho\,\mathrm{e}^{-\xi^*\hat a}\mathrm{e}^{\xi\hat a^\dagger}) で、χ0\chi_0 は並べ替えない Tr⁡(ρ^D^(ξ))\operatorname{Tr}(\hat\rho\hat D(\xi)) である。この族は Cahill と Glauber が導入した(K. E. Cahill and R. J. Glauber (1969), Phys. Rev. 177, 1857;同 1882)。PP 表示は R. J. Glauber (1963), Phys. Rev. 131, 2766 と E. C. G. Sudarshan (1963), Phys. Rev. Lett. 10, 277 が、QQ 関数は K. Husimi (1940), Proc. Phys.-Math. Soc. Japan 22, 264–314 が最初に書いた。

χ0\chi_0 が前の章の特性関数と同じものであることを確かめておく。ξ=ξ1+iξ2\xi=\xi_1+\mathrm{i}\xi_2 と書き、(5.6.1) で a^\hat a を x^,p^\hat x,\hat p に戻すとξa^†−ξ∗a^=i(ax^+bp^)/ℏ\xi\hat a^\dagger-\xi^*\hat a=\mathrm{i}(a\hat x+b\hat p)/\hbar、a=−ξ22ℏmωa=-\xi_2\sqrt{2\hbar m\omega}、b=ξ12ℏ/(mω)b=\xi_1\sqrt{2\hbar/(m\omega)} となる。したがって χ0(ξ)\chi_0(\xi) は 第5部第2章の量子の特性関数 χ(a,b)=Tr⁡(ρ^ ei(ax^+bp^)/ℏ)\chi(a,b)=\operatorname{Tr}(\hat\rho\,\mathrm{e}^{\mathrm{i}(a\hat x+b\hat p)/\hbar}) そのもので、同じ章で示した「WW は特性関数の Fourier 逆変換」から、(5.6.4) の F0F_0 は (5.6.2) の W(α)W(\alpha) に一致する。

並べ方の定理と、Gauss 関数の畳み込み

定理 5.6.2(ss 順序の期待値 ★★★)

FsF_s が関数として存在するとき、a^†\hat a^\dagger を mm 個、a^\hat a を nn 個含む ss 順序の積 {(a^†)ma^n}s\{(\hat a^\dagger)^m\hat a^n\}_s について

Tr⁡(ρ^{(a^†)ma^n}s)=∫Fs(α) (α∗)mαn d2α.(5.6.5)\operatorname{Tr}\bigl(\hat\rho\{(\hat a^\dagger)^m\hat a^n\}_s\bigr)=\int F_s(\alpha)\,(\alpha^*)^m\alpha^n\,\mathrm{d}^2\alpha .\tag{5.6.5}

s=1s=1 では {⋅}1\{\cdot\}_1 は正規順序 (a^†)ma^n(\hat a^\dagger)^m\hat a^n、s=−1s=-1 では反正規順序 a^n(a^†)m\hat a^n(\hat a^\dagger)^m、s=0s=0 では m+nm+n 個の演算子の並べ方すべての平均(対称順序)である。

証明

ss 順序の積は母関数 eξa^†−ξ∗a^es∣ξ∣2/2=∑m,nξm(−ξ∗)nm! n!{(a^†)ma^n}s\mathrm{e}^{\xi\hat a^\dagger-\xi^*\hat a}\mathrm{e}^{s\lvert \xi\rvert ^2/2}=\sum_{m,n}\frac{\xi^m(-\xi^*)^n}{m!\,n!}\{(\hat a^\dagger)^m\hat a^n\}_s で定義する。s=1s=1 では (5.6.3) の中央の形から左辺は eξa^†e−ξ∗a^\mathrm{e}^{\xi\hat a^\dagger}\mathrm{e}^{-\xi^*\hat a} で、展開の係数は正規順序の積になる。s=−1s=-1 では右の形から反正規順序になる。s=0s=0 では (ξa^†−ξ∗a^)m+n(\xi\hat a^\dagger-\xi^*\hat a)^{m+n} を展開すると a^†\hat a^\dagger と a^\hat a のすべての並べ方が同じ重みで現れるので、対称順序になる。母関数の両辺のトレースを取ると、左辺は χs(ξ)\chi_s(\xi) である。一方 (5.6.4) の Fourier 変換を逆に解くと χs(ξ)=∫Fs(α)eξα∗−ξ∗αd2α\chi_s(\xi)=\int F_s(\alpha)\mathrm{e}^{\xi\alpha^*-\xi^*\alpha}\mathrm{d}^2\alpha で(∫e(α−β)ξ∗−(α−β)∗ξd2ξ=π2δ2(α−β)\int\mathrm{e}^{(\alpha-\beta)\xi^*-(\alpha-\beta)^*\xi}\mathrm{d}^2\xi=\pi^2\delta^2(\alpha-\beta) を使う)、eξα∗−ξ∗α=∑m,nξm(−ξ∗)nm! n!(α∗)mαn\mathrm{e}^{\xi\alpha^*-\xi^*\alpha}=\sum_{m,n}\frac{\xi^m(-\xi^*)^n}{m!\,n!}(\alpha^*)^m\alpha^n と展開してξm(−ξ∗)n\xi^m(-\xi^*)^n の係数を比べると (5.6.5) を得る。

定理の要点は、量子の期待値が古典の確率分布と同じ積分 ∫F f d2α\int F\,f\,\mathrm{d}^2\alpha の形で書けるが、その代わりに演算子の並べ方を1つ選ばなければならないことである。選んだ並べ方に応じて分布が変わる。異なる ss の分布の関係は、(5.6.4) から直ちに読み取ることができる。

命題 5.6.3(ss を下げると Gauss 関数でならす ★★★)

Gσ(α)=1πσe−∣α∣2/σG_\sigma(\alpha)=\frac{1}{\pi\sigma}\mathrm{e}^{-\lvert \alpha\rvert ^2/\sigma}(σ>0\sigma>0、∫Gσd2α=1\int G_\sigma\mathrm{d}^2\alpha=1、各直交位相の分散 σ/2\sigma/2)とする。s>s′s>s' のとき

Fs′=Fs∗G(s−s′)/2,とくにW=P∗G1/2,Q=W∗G1/2=P∗G1.(5.6.6)F_{s'}=F_s*G_{(s-s')/2},\qquad\text{とくに}\quad W=P*G_{1/2},\quad Q=W*G_{1/2}=P*G_1 .\tag{5.6.6}

証明

(5.6.4) から χs′(ξ)=χs(ξ)e−(s−s′)∣ξ∣2/2\chi_{s'}(\xi)=\chi_s(\xi)\mathrm{e}^{-(s-s')\lvert \xi\rvert ^2/2}。GσG_\sigma の特性関数は ∫Gσ(α)eξα∗−ξ∗αd2α=e−σ∣ξ∣2\int G_\sigma(\alpha)\mathrm{e}^{\xi\alpha^*-\xi^*\alpha}\mathrm{d}^2\alpha=\mathrm{e}^{-\sigma\lvert \xi\rvert ^2}(ξ=k1+ik2\xi=k_1+\mathrm{i} k_2 と書くと指数は 2i(k2X1−k1X2)2\mathrm{i}(k_2X_1-k_1X_2) で、各成分の Gauss 積分)なので、e−(s−s′)∣ξ∣2/2\mathrm{e}^{-(s-s')\lvert \xi\rvert ^2/2} は G(s−s′)/2G_{(s-s')/2} の特性関数である。特性関数の積は分布の畳み込みに対応するので (5.6.6) を得る。Gσ∗Gσ′=Gσ+σ′G_{\sigma}*G_{\sigma'}=G_{\sigma+\sigma'} も同じ理由で成り立つ。

QQ と PP には、コヒーレント状態 ∣α⟩| \alpha \rangle を使った別の書き方がある。反正規順序の期待値 Tr⁡(ρ^ a^n(a^†)m)\operatorname{Tr}(\hat\rho\,\hat a^n(\hat a^\dagger)^m) の a^n\hat a^n と (a^†)m(\hat a^\dagger)^m のあいだにコヒーレント状態の完全性 ∫∣α⟩⟨α∣ d2α/π=1\int| \alpha \rangle \langle \alpha |\,\mathrm{d}^2\alpha/\pi=1 を挟み、a^∣α⟩=α∣α⟩\hat a| \alpha \rangle =\alpha| \alpha \rangle を使うと

Tr⁡(ρ^ a^n(a^†)m)=∫⟨α∣ρ^∣α⟩π αn(α∗)m d2α\operatorname{Tr}\bigl(\hat\rho\,\hat a^n(\hat a^\dagger)^m\bigr)=\int\frac{\langle \alpha | \hat\rho | \alpha \rangle }{\pi}\,\alpha^n(\alpha^*)^m\,\mathrm{d}^2\alpha

となる。これがすべての m,nm,n で (5.6.5) の s=−1s=-1 と一致し、特性関数が一致するので

Q(α)=⟨α∣ρ^∣α⟩π ≥0(5.6.7)Q(\alpha)=\frac{\langle \alpha | \hat\rho | \alpha \rangle }{\pi}\ \ge0\tag{5.6.7}

である。逆に、ρ^\hat\rho をコヒーレント状態の射影の重ね合わせ

ρ^=∫P(α)∣α⟩⟨α∣ d2α(5.6.8)\hat\rho=\int P(\alpha)| \alpha \rangle \langle \alpha |\,\mathrm{d}^2\alpha\tag{5.6.8}

に書けたとすると、正規順序の期待値は Tr⁡(ρ^(a^†)ma^n)=∫P(α)⟨α∣(a^†)ma^n∣α⟩d2α=∫P (α∗)mαnd2α\operatorname{Tr}(\hat\rho(\hat a^\dagger)^m\hat a^n)=\int P(\alpha)\langle \alpha | (\hat a^\dagger)^m\hat a^n | \alpha \rangle \mathrm{d}^2\alpha=\int P\,(\alpha^*)^m\alpha^n\mathrm{d}^2\alpha となり、s=1s=1 の (5.6.5) と一致する。(5.6.8) の重み PP が PP 表示である。PP が非負の関数なら、(5.6.8) は「コヒーレント状態(古典の場)を確率 PP で混ぜたもの」と読むことができる。

第5部第1章の要請との関係

第5部第1章は、位相空間の分布への要請 (W1)〜(W5) のうち (W2) 線形性・(W3) 周辺分布・(W5) 非負の3つが両立しないことを示し、WW は (W5) を捨てて残りを守る選び方だった。QQ は (5.6.7) から (W5) を守り、代わりに (W3) を捨てる。QQ の X1X_1 の周辺分布は真の X1X_1 の分布を分散 14\frac14 の Gauss 関数でならしたもので、量子の揺らぎが二重に数えられる。PP はどちらも守らず、正規順序の期待値を与えることだけを要請に置く。3つの分布は、同じ密度行列の情報を、別々の要請の組で書いたものである。

例題 5.6.1 1光子の状態 ∣1⟩| 1 \rangle の3つの分布 ★★

オリジナル∣1⟩| 1 \rangle の QQ を (5.6.7) から求めよ。次に P(α)=(1+∂α∂α∗)δ2(α)P(\alpha)=\bigl(1+\partial_\alpha\partial_{\alpha^*}\bigr)\delta^2(\alpha) と置くと、(5.6.6) の Q=P∗G1Q=P*G_1 が正しい QQ を与えることを確かめ、同じ PP から W=P∗G1/2W=P*G_{1/2} を求めて、第5部第5章の W1W_1 と照合せよ。ここで α\alpha と α∗\alpha^* は独立な変数として微分する。

解答

⟨α∣1⟩=e−∣α∣2/2α∗\langle \alpha | 1 \rangle =\mathrm{e}^{-\lvert \alpha\rvert ^2/2}\alpha^* なので Q=∣α∣2e−∣α∣2/πQ=\lvert \alpha\rvert ^2\mathrm{e}^{-\lvert \alpha\rvert ^2}/\pi。

δ2\delta^2 の微分との畳み込みは、微分を Gauss 関数の方へ移せばよく、P∗Gσ=(1+∂α∂α∗)GσP*G_\sigma=(1+\partial_\alpha\partial_{\alpha^*})G_\sigma である。∂α∗e−αα∗/σ=−ασe−αα∗/σ\partial_{\alpha^*}\mathrm{e}^{-\alpha\alpha^*/\sigma}=-\frac{\alpha}{\sigma}\mathrm{e}^{-\alpha\alpha^*/\sigma}、もう一度 ∂α\partial_\alpha を当てて∂α∂α∗e−∣α∣2/σ=(∣α∣2σ2−1σ)e−∣α∣2/σ\partial_\alpha\partial_{\alpha^*}\mathrm{e}^{-\lvert \alpha\rvert ^2/\sigma}=\bigl(\frac{\lvert \alpha\rvert ^2}{\sigma^2}-\frac1\sigma\bigr)\mathrm{e}^{-\lvert \alpha\rvert ^2/\sigma} だから

(1+∂α∂α∗)Gσ=1πσ(1−1σ+∣α∣2σ2)e−∣α∣2/σ.(1+\partial_\alpha\partial_{\alpha^*})G_\sigma=\frac{1}{\pi\sigma}\Bigl(1-\frac1\sigma+\frac{\lvert \alpha\rvert ^2}{\sigma^2}\Bigr)\mathrm{e}^{-\lvert \alpha\rvert ^2/\sigma} .

σ=1\sigma=1 では ∣α∣2e−∣α∣2/π\lvert \alpha\rvert ^2\mathrm{e}^{-\lvert \alpha\rvert ^2}/\pi で、上の QQ と一致する。σ=12\sigma=\frac12 では

W(α)=2π(4∣α∣2−1)e−2∣α∣2.W(\alpha)=\frac2\pi\bigl(4\lvert \alpha\rvert ^2-1\bigr)\mathrm{e}^{-2\lvert \alpha\rvert ^2} .

u=2H/(ℏω)=2∣α∣2u=2H/(\hbar\omega)=2\lvert \alpha\rvert ^2 と (5.6.2) で x,px,p の単位に戻すと W(x,p)=1πℏ(2u−1)e−u=−1πℏe−uL1(2u)W(x,p)=\frac{1}{\pi\hbar}(2u-1)\mathrm{e}^{-u}=-\frac{1}{\pi\hbar}\mathrm{e}^{-u}L_1(2u) で、第5部第5章の W1W_1 と同じである。3つを並べると、QQ は原点で 00 になる非負の関数、WW は ∣α∣2<14\lvert \alpha\rvert ^2<\frac14 で負、PP は δ\delta 関数の2階微分を含む超関数で、関数としては存在しない。ss を下げるたびに Gauss 関数でならされ、特異さが減っていく。

木下の勉強ステータス

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

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

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