tecktom.com

Moyal 括弧と Moyal 方程式

2026年10月17日10分で読めますAI生成教科書『虎の巻』第5部 第3章 ⋆ 積と Moyal 方程式
目次
  1. Moyal 括弧
  2. Moyal 方程式
  3. の場合
  4. 木下の層

AI 生成・木下が未勉強

量子力学の様々な構成方法 第5部 第3章 §5.3.3 ★★★ 交換子は Moyal 括弧に、Heisenberg 方程式と von Neumann 方程式は Moyal 方程式に写る。Weyl 対応が1対1なので同値は両向きに成り立つ

Moyal 括弧

定義 5.3.4(Moyal 括弧)

{ ⁣{f,g} ⁣}=f⋆g−g⋆fiℏ\{\!\{f,g\}\!\}=\dfrac{f\star g-g\star f}{\mathrm{i}\hbar} を Moyal 括弧と呼ぶ。定理 5.3.2 から

(1iℏ[A^, B^])W={ ⁣{AW,BW} ⁣}.(5.3.8)\Bigl(\frac{1}{\mathrm{i}\hbar}\left[\hat A,\,\hat B\right]\Bigr)_W=\{\!\{A_W,B_W\}\!\} .\tag{5.3.8}

⋆\star 積の級数 (5.3.4) で ff と gg を入れ替えると、Π↔\overleftrightarrow\Pi の符号が変わるので nn 次の項に (−1)n(-1)^n が掛かる。差を取ると偶数次の項が消え、奇数次の項が2倍になる。iℏ\mathrm{i}\hbar で割ると

{ ⁣{f,g} ⁣}=2ℏ fsin⁡[ℏ2Π↔]g={f,g}−ℏ224(fxxxgppp−3fxxpgxpp+3fxppgxxp−fpppgxxx)+O(ℏ4)(5.3.9)\{\!\{f,g\}\!\}=\frac{2}{\hbar}\,f\sin\Bigl[\frac{\hbar}{2}\overleftrightarrow\Pi\Bigr]g =\{f,g\}-\frac{\hbar^2}{24}\bigl(f_{xxx}g_{ppp}-3f_{xxp}g_{xpp}+3f_{xpp}g_{xxp}-f_{ppp}g_{xxx}\bigr)+\mathcal{O}(\hbar^4)\tag{5.3.9}

となる。ℏ\hbar の1次の項 iℏ2{f,g}\frac{\mathrm{i}\hbar}{2}\{f,g\} が2倍されて iℏ\mathrm{i}\hbar で割られ、Poisson 括弧がちょうど係数1で残る。次の補正は ℏ2\hbar^2 で、ff と gg を合わせて6回微分する。展開の形は 第10部第1章で Groenewold の矛盾を読むときに使った式と同じで、ここで導いた。Moyal 括弧は Poisson 括弧と同じく反対称で、Jacobi 恒等式も満たす(交換子の Jacobi 恒等式を写したもの)。

2次以下なら Moyal 括弧=Poisson 括弧

ff と gg のどちらかが xx、pp の2次以下の多項式なら、3階以上の微分が消えるので { ⁣{f,g} ⁣}={f,g}\{\!\{f,g\}\!\}=\{f,g\} が厳密に成り立つ。ハミルトニアンが2次以下(自由粒子・調和振動子・一様な力)なら、量子の時間発展を位相空間で書いた式は古典の式と1文字も違わない。

Moyal 方程式

定理 5.3.5(Heisenberg 方程式・von Neumann 方程式 ⇔ Moyal 方程式 ★★★)

ハミルトニアン H^\hat H の Weyl 記号を H(x,p)H(x,p) とする。

  1. (a) 密度行列が von Neumann 方程式 ∂tρ^=[H^, ρ^]/(iℏ)\partial_t\hat\rho=\left[\hat H,\,\hat\rho\right]/(\mathrm{i}\hbar) に従うことと、Wigner 関数が

∂W∂t={ ⁣{H,W} ⁣}(5.3.10)\frac{\partial W}{\partial t}=\{\!\{H,W\}\!\}\tag{5.3.10}

(Moyal 方程式)に従うことは同値である。2. (b) Heisenberg 描像の演算子が dA^Hdt=[A^H, H^]/(iℏ)\frac{\mathrm{d} \hat A_{\mathrm H}}{\mathrm{d} t}=\left[\hat A_{\mathrm H},\,\hat H\right]/(\mathrm{i}\hbar) に従うことと、その記号 A(x,p,t)=(A^H(t))WA(x,p,t)=(\hat A_{\mathrm H}(t))_W が∂tA={ ⁣{A,H} ⁣}\partial_tA=\{\!\{A,H\}\!\} に従うことは同値である。

証明

(a) 前の章の命題から W=ρW/(2πℏ)W=\rho_W/(2\pi\hbar)。Weyl 記号を取る操作は線形で時間を含まないので、∂t\partial_t と交換する。von Neumann 方程式の両辺の記号を取り、(5.3.8) を使うと ∂tρW={ ⁣{H,ρW} ⁣}\partial_t\rho_W=\{\!\{H,\rho_W\}\!\}、2πℏ2\pi\hbar で割って (5.3.10) を得る。逆に (5.3.10) が成り立てば、両辺に 2πℏ QW2\pi\hbar\,\mathcal Q_{\mathrm W} を当てて (5.3.8) を逆に読むと von Neumann 方程式に戻る。Weyl 対応は1対1なので(前の章の定理)、どちらの向きでも情報は失われない。(b) も同じ議論で、交換子の順序が [A^H, H^]\left[\hat A_{\mathrm H},\,\hat H\right] であることに応じて括弧の順序が { ⁣{A,H} ⁣}\{\!\{A,H\}\!\} になる。

(a) と (b) は、期待値が描像によらないこと

ddt∫WAW dx dp=∫{ ⁣{H,W} ⁣}AW dx dp=∫W{ ⁣{AW,H} ⁣} dx dp\frac{\mathrm{d} }{\mathrm{d} t}\int W A_W\,\mathrm{d} x\,\mathrm{d} p=\int\{\!\{H,W\}\!\}A_W\,\mathrm{d} x\,\mathrm{d} p=\int W\{\!\{A_W,H\}\!\}\,\mathrm{d} x\,\mathrm{d} p

の2通りの書き方に当たる。2つ目の等号は命題 5.3.3 (iii) と結合律から ∫(H⋆W)⋆A=∫W⋆(A⋆H)\int(H\star W)\star A=\int W\star(A\star H) などとして、積分の中で ⋆\star 積を巡回させてよいことによる。この式で AWA_W を xjpkx^jp^k に取れば、位相空間のモーメントの時間変化が Moyal 括弧1つで書ける。

H=p2/(2m)+V(x)H=p^2/(2m)+V(x) の場合

ハミルトニアンが H^=p^2/(2m)+V(x^)\hat H=\hat p^2/(2m)+V(\hat x) なら、記号は H=p2/(2m)+V(x)H=p^2/(2m)+V(x) そのものである(前の章:xx だけ、pp だけの関数の記号は古典の関数と同じ)。(5.3.4) を運動エネルギーとポテンシャルに分けて当てる。V(x)V(x) は pp で微分できないので、V⋆WV\star W では k=0k=0 の項だけ、W⋆VW\star V では k=nk=n の項だけが残り

V⋆W−W⋆V=∑n1n!(iℏ2)n(1−(−1)n)V(n) ∂pnW.V\star W-W\star V=\sum_{n}\frac{1}{n!}\Bigl(\frac{\mathrm{i}\hbar}{2}\Bigr)^n\bigl(1-(-1)^n\bigr)V^{(n)}\,\partial_p^nW .

奇数の nn だけが残り、iℏ\mathrm{i}\hbar で割ると係数は 1n!(iℏ/2)n−1\frac{1}{n!}(\mathrm{i}\hbar/2)^{n-1} になる。運動エネルギーの側も同じ計算で、p2p^2 は3階以上の微分が消えるので n=1n=1 だけが残る。まとめて

∂W∂t=−pm∂W∂x+V′(x)∂W∂p−ℏ224V′′′(x)∂3W∂p3+ℏ41920V(5)(x)∂5W∂p5−⋯(5.3.11)\frac{\partial W}{\partial t}=-\frac pm\frac{\partial W}{\partial x}+V'(x)\frac{\partial W}{\partial p}-\frac{\hbar^2}{24}V'''(x)\frac{\partial^3W}{\partial p^3}+\frac{\hbar^4}{1920}V^{(5)}(x)\frac{\partial^5W}{\partial p^5}-\cdots\tag{5.3.11}

を得る。初めの2項は古典の Liouville 方程式の右辺 {H,W}=HxWp−HpWx\{H,W\}=H_xW_p-H_pW_x と同じで、3項目以降が量子の補正である。補正の項はポテンシャルの3階以上の微分と運動量の3階以上の微分の組で、ポテンシャルが4次の多項式なら ℏ2\hbar^2 の項で厳密に止まる。

自由粒子では V=0V=0 で、(5.3.11) は ∂tW=−(p/m)∂xW\partial_tW=-(p/m)\partial_xW だけになる。Wt(x,p)=W0(x−pt/m,p)W_t(x,p)=W_0(x-pt/m,p) を入れると∂tWt=−(p/m)(∂xW0)(x−pt/m,p)=−(p/m)∂xWt\partial_tW_t=-(p/m)(\partial_xW_0)(x-pt/m,p)=-(p/m)\partial_xW_t で、確かに解である。共通例題③の剪断の解(第5部第2章の5節で特性関数から導いた Wt=W0(x−pt/m,p)W_t=W_0(x-pt/m,p) と、そこから出る σ(t)2=σ02+(ℏt/(2mσ0))2\sigma(t)^2=\sigma_0^2+(\hbar t/(2m\sigma_0))^2)が、Moyal 方程式の立場からも出てくる。前の章は Heisenberg 方程式が線形に閉じることを使ったが、ここでは Moyal 括弧が2次の HH で Poisson 括弧に一致することを使っている。

例題 5.3.3 3次のポテンシャルでのモーメントの時間変化 ★★

オリジナルV(x)=κx3/3V(x)=\kappa x^3/3 とする。

  1. (1) { ⁣{p,H} ⁣}\{\!\{p,H\}\!\}、{ ⁣{p2,H} ⁣}\{\!\{p^2,H\}\!\}、{ ⁣{p3,H} ⁣}\{\!\{p^3,H\}\!\} を求めよ。
  2. (2) ddt⟨p^3⟩\frac{\mathrm{d} }{\mathrm{d} t}\langle \hat p^3 \rangle に、古典の統計力学には無い項が、状態によらない定数として現れることを示せ。
解答

(1) pkp^k は xx で微分できないので、(5.3.9) の ℏ2\hbar^2 の項で残るのは −fpppgxxx-f_{ppp}g_{xxx} の項だけで、pkp^k と VV の Moyal 括弧は{pk,V}+ℏ224(pk)pppV′′′\{p^k,V\}+\frac{\hbar^2}{24}(p^k)_{ppp}V''' である(ℏ4\hbar^4 の項は VV の5階微分を含み、3次の VV では消える)。運動エネルギーとは可換なので

{ ⁣{p,H} ⁣}=−κx2,{ ⁣{p2,H} ⁣}=−2κx2p,{ ⁣{p3,H} ⁣}=−3κx2p2+ℏ224⋅6⋅2κ=−3κx2p2+κℏ22.\begin{gathered} \{\!\{p,H\}\!\}=-\kappa x^2,\qquad \{\!\{p^2,H\}\!\}=-2\kappa x^2p,\\ \{\!\{p^3,H\}\!\}=-3\kappa x^2p^2+\frac{\hbar^2}{24}\cdot6\cdot2\kappa=-3\kappa x^2p^2+\frac{\kappa\hbar^2}{2} . \end{gathered}

(2) 定理 5.3.5 (b) の後の式から ddt⟨p^3⟩=∫W{ ⁣{p3,H} ⁣}=−3κ∫Wx2p2+κℏ2/2\frac{\mathrm{d} }{\mathrm{d} t}\langle \hat p^3 \rangle =\int W\{\!\{p^3,H\}\!\}=-3\kappa\int Wx^2p^2+\kappa\hbar^2/2。古典の分布 ρ\rho が Liouville 方程式に従うなら ddt∫ρ p3=−3κ∫ρ x2p2\frac{\mathrm{d} }{\mathrm{d} t}\int\rho\,p^3=-3\kappa\int\rho\,x^2p^2 で、量子では κℏ2/2\kappa\hbar^2/2 が足される。この項は ∫W=1\int W=1 だけを使っていて、状態によらない。pp の1次と2次のモーメントの式は古典と同じ形なので、量子と古典の違いが最初に見えるのは3次のモーメントからである。ただし ∫Wx2p2\int Wx^2p^2 は Weyl 順序の演算子 QW[x2p2]\mathcal Q_{\mathrm W}[x^2p^2] の期待値で、古典の ∫ρ x2p2\int\rho\,x^2p^2 と同じ記号でも中身は違う。

木下の勉強ステータス

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

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

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