tecktom.com

逆向き:位相空間の関数から演算子を作る

2026年10月17日8分で読めますAI生成教科書『虎の巻』第5部 第2章 Wigner 変換と Weyl 対応
目次
  1. 逆の写像
  2. 実数とエルミート、トレースと積分
  3. 木下の層

AI 生成・木下が未勉強

量子力学の様々な構成方法 第5部 第2章 §5.2.2 ★★★ Wigner 変換は Fourier 変換なので逆を書くことができる。WW から ρ^\hat\rho が元どおりに復元できる

逆の写像

Wigner 変換 (5.2.2) は、中点 xx を止めると差 yy についての Fourier 変換である。Fourier 変換は逆変換を持つので

⟨x−y2∣A^∣x+y2⟩=∫dp2πℏ AW(x,p) e−ipy/ℏ.\langle x-\tfrac y2 | \hat A | x+\tfrac y2 \rangle =\int\frac{\mathrm{d} p}{2\pi\hbar}\,A_W(x,p)\,\mathrm{e}^{-\mathrm{i} py/\hbar} .

左辺を元の座標 x′=x−y/2x'=x-y/2、x′′=x+y/2x''=x+y/2 に戻すと y=x′′−x′y=x''-x'、x=(x′+x′′)/2x=(x'+x'')/2 で、次の式を得る。

定義 5.2.3(Weyl 量子化)

位相空間の関数 f(x,p)f(x,p) に対し、演算子 QW[f]\mathcal Q_{\mathrm W}[f] を位置表示の行列要素

⟨x′∣QW[f]∣x′′⟩=∫dp2πℏ f(x′+x′′2,p) eip(x′−x′′)/ℏ(5.2.6)\langle x' | \mathcal Q_{\mathrm W}[f] | x'' \rangle =\int\frac{\mathrm{d} p}{2\pi\hbar}\,f\Bigl(\frac{x'+x''}{2},p\Bigr)\,\mathrm{e}^{\mathrm{i} p(x'-x'')/\hbar}\tag{5.2.6}

で定める。f↦QW[f]f\mapsto\mathcal Q_{\mathrm W}[f] を Weyl 量子化、この対応 A^↔AW\hat A\leftrightarrow A_W を Weyl 対応と呼ぶ。

作り方から、QW\mathcal Q_{\mathrm W} と Wigner 変換は互いに逆である。

QW[AW]=A^,(QW[f])W=f.(5.2.7)\mathcal Q_{\mathrm W}[A_W]=\hat A,\qquad \bigl(\mathcal Q_{\mathrm W}[f]\bigr)_W=f .\tag{5.2.7}

名前の由来は H. Weyl が 1927 年の論文(H. Weyl (1927)、Z. Phys. 46、1–46)で、古典の関数から演算子を作る規則としてこの対応を(次の 5.2.4 節の Fourier 展開の形で)導入したことにある。Wigner の 1932 年の関数がその逆向きの写像を状態に当てたものであることは、H. J. Groenewold の学位論文(H. J. Groenewold (1946)、Physica 12、405–460)で明確にされた。

密度行列に当てると、(5.2.3) から ρ^=QW[2πℏ W]\hat\rho=\mathcal Q_{\mathrm W}[2\pi\hbar\,W]、つまり

⟨x′∣ρ^∣x′′⟩=∫dp W(x′+x′′2,p) eip(x′−x′′)/ℏ(5.2.8)\langle x' | \hat\rho | x'' \rangle =\int\mathrm{d} p\,W\Bigl(\frac{x'+x''}{2},p\Bigr)\,\mathrm{e}^{\mathrm{i} p(x'-x'')/\hbar}\tag{5.2.8}

である。純粋状態ならこれが ψ(x′)ψ∗(x′′)\psi(x')\psi^*(x'') に等しく、x′′x'' を1点に固定すれば ψ(x′)\psi(x') が全体の位相を除いて読み取ることができる。WW を知っていれば ρ^\hat\rho のすべての行列要素が分かるので、WW は状態の情報を何も失っていない。

実数とエルミート、トレースと積分

Weyl 対応は、演算子の基本的な性質を関数の性質へ写す。

命題 5.2.4(エルミート共役と複素共役 ★★)

(A^†)W=(AW)∗(\hat A^\dagger)_W=(A_W)^*。とくに A^\hat A がエルミートであることと AWA_W が実数値であることは同値である。

証明

⟨x′∣A^†∣x′′⟩=⟨x′′∣A^∣x′⟩∗\langle x' | \hat A^\dagger | x'' \rangle =\langle x'' | \hat A | x' \rangle ^* なので

(A^†)W=∫⟨x+y2∣A^∣x−y2⟩∗ eipy/ℏ dy=(∫⟨x+y2∣A^∣x−y2⟩ e−ipy/ℏ dy)∗.(\hat A^\dagger)_W=\int\langle x+\tfrac y2 | \hat A | x-\tfrac y2 \rangle ^*\,\mathrm{e}^{\mathrm{i} py/\hbar}\,\mathrm{d} y =\Bigl(\int\langle x+\tfrac y2 | \hat A | x-\tfrac y2 \rangle \,\mathrm{e}^{-\mathrm{i} py/\hbar}\,\mathrm{d} y\Bigr)^* .

括弧の中で y→−yy\to-y と替えると (5.2.2) の AWA_W に戻る。後半は、Wigner 変換が1対1であることからA^†=A^  ⟺  (A^†)W=AW  ⟺  AW∗=AW\hat A^\dagger=\hat A\iff(\hat A^\dagger)_W=A_W\iff A_W^*=A_W。

トレースは対角成分の和なので、y=0y=0 の成分だけが寄与する。(5.2.6) で x′=x′′=xx'=x''=x と置いて xx で積分すると

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

を得る。演算子のトレースは、記号を相空間で積分して面積 2πℏ2\pi\hbar で割ったものである。ρ^\hat\rho に当てれば Tr⁡ρ^=∫W dx dp=1\operatorname{Tr}\hat\rho=\int W\,\mathrm{d} x\,\mathrm{d} p=1 で、前の章の規格化 (W1) が出る。2πℏ2\pi\hbar は前の章で「状態1つが占める面積」として現れた量で、学部の統計力学で古典の分配関数に入れる 1/h1/h と同じものである。

例題 5.2.2 熱平衡の WW から密度行列を復元する ★★★

オリジナル共通例題①の熱平衡の Wigner 関数は、s=βℏωs=\beta\hbar\omega として、中心が原点で相関の無い Gauss 関数

W(x,p)=12πσxσpexp⁡(−x22σx2−p22σp2),σx2=ℏ2mωcoth⁡s2,σp2=mℏω2coth⁡s2W(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),\qquad \sigma_x^2=\frac{\hbar}{2m\omega}\coth\frac s2,\quad \sigma_p^2=\frac{m\hbar\omega}{2}\coth\frac s2

である(第5部第1章の総合問題)。

  1. (1) (5.2.8) で位置表示の密度行列 ⟨x′∣ρ^∣x′′⟩\langle x' | \hat\rho | x'' \rangle を求めよ。
  2. (2) 結果が、⟨x′∣e−βH^∣x′′⟩=mω2πℏsinh⁡sexp⁡{−mω2ℏsinh⁡s[(x′2+x′′2)cosh⁡s−2x′x′′]}\langle x' | \mathrm{e}^{-\beta\hat H} | x'' \rangle =\sqrt{\frac{m\omega}{2\pi\hbar\sinh s}}\exp\Bigl\{-\frac{m\omega}{2\hbar\sinh s}\bigl[(x'^2+x''^2)\cosh s-2x'x''\bigr]\Bigr\}を Z=1/(2sinh⁡(s/2))Z=1/(2\sinh(s/2)) で割ったものに一致することを、s=0.7s=0.7、ℏ=m=ω=1\hbar=m=\omega=1、(x′,x′′)=(0.3,−0.5)(x',x'')=(0.3,-0.5) で数値で確かめよ。
解答

(1) X=(x′+x′′)/2X=(x'+x'')/2、u=x′−x′′u=x'-x'' と置く。pp についての Gauss 積分 ∫dp e−p2/(2σp2)eipu/ℏ=2π σp e−σp2u2/(2ℏ2)\int\mathrm{d} p\,\mathrm{e}^{-p^2/(2\sigma_p^2)}\mathrm{e}^{\mathrm{i} pu/\hbar}=\sqrt{2\pi}\,\sigma_p\,\mathrm{e}^{-\sigma_p^2u^2/(2\hbar^2)} から

⟨x′∣ρ^∣x′′⟩=12π σxexp⁡(−X22σx2−σp2u22ℏ2).\langle x' | \hat\rho | x'' \rangle =\frac{1}{\sqrt{2\pi}\,\sigma_x}\exp\Bigl(-\frac{X^2}{2\sigma_x^2}-\frac{\sigma_p^2u^2}{2\hbar^2}\Bigr) .

密度行列は、中点 XX については幅 σx\sigma_x、差 uu については幅 ℏ/σp\hbar/\sigma_p の Gauss 関数の積になる。差の方向の幅は「干渉が残る距離」で、温度が上がって σp\sigma_p が大きくなるほど短くなる。

(2) ℏ=m=ω=1\hbar=m=\omega=1、s=0.7s=0.7 では coth⁡0.35=2.9729\coth0.35=2.9729 で、σx2=σp2=1.4864\sigma_x^2=\sigma_p^2=1.4864。X=−0.1X=-0.1、u=0.8u=0.8 を入れると、(1) の結果は

12π×1.4864exp⁡(−0.012×1.4864−1.4864×0.642)=0.3272×e−0.4790=0.2027.\frac{1}{\sqrt{2\pi\times1.4864}}\exp\Bigl(-\frac{0.01}{2\times1.4864}-\frac{1.4864\times0.64}{2}\Bigr)=0.3272\times\mathrm{e}^{-0.4790}=0.2027 .

右辺は sinh⁡0.7=0.7586\sinh0.7=0.7586、cosh⁡0.7=1.2552\cosh0.7=1.2552、Z=1/(2sinh⁡0.35)=1.3998Z=1/(2\sinh0.35)=1.3998 で、前の係数は 1/(2π×0.7586)/1.3998=0.3272\sqrt{1/(2\pi\times0.7586)}/1.3998=0.3272、指数は −12×0.7586[0.34×1.2552+0.30]=−0.4790-\frac{1}{2\times0.7586}\bigl[0.34\times1.2552+0.30\bigr]=-0.4790 となり、同じ 0.20270.2027 を与える。係数も指数も別々に一致しているのは偶然でなく、cosh⁡s−1=2sinh⁡2s2\cosh s-1=2\sinh^2\frac s2 などを使うと両辺ともmωπℏtanh⁡s2exp⁡[−mωℏtanh⁡s2 X2−mω4ℏcoth⁡s2 u2]\sqrt{\frac{m\omega}{\pi\hbar}\tanh\frac s2}\exp\bigl[-\tfrac{m\omega}{\hbar}\tanh\tfrac s2\,X^2-\tfrac{m\omega}{4\hbar}\coth\tfrac s2\,u^2\bigr] に整理することができる。

木下の勉強ステータス

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

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

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