tecktom.com

この章のまとめ

2026年10月17日9分で読めますAI生成教科書『虎の巻』第5部 第3章 ⋆ 積と Moyal 方程式

AI 生成・木下が未勉強

量子力学の様々な構成方法 第5部 第3章 §5.3.7 ★★★ 積は ⋆\star 積、交換子は Moyal 括弧、時間発展は Moyal 方程式に写り、ℏ→0\hbar\to0 で古典の Liouville 方程式に戻る

指数関数の演算子 ei(ax^+bp^)/ℏ\mathrm{e}^{\mathrm{i}(a\hat x+b\hat p)/\hbar} の積は、肩の交換子が数なので BCH の位相1つで閉じる。その位相を微分で書き、Fourier の重ね合わせで一般の演算子に広げると、積の Weyl 記号は ⋆\star 積 f⋆g=fexp⁡[iℏ2(∂x←∂p→−∂p←∂x→)]gf\star g=f\exp[\frac{\mathrm{i}\hbar}{2}(\overleftarrow{\partial_x}\overrightarrow{\partial_p}-\overleftarrow{\partial_p}\overrightarrow{\partial_x})]g で与えられる。⋆\star 積は結合的で可換でなく、積分の中では普通の積に戻せ、x⋆x\star と p⋆p\star は「半分は掛け算、半分は微分」の微分演算子になる。交換子の記号は Moyal 括弧 { ⁣{f,g} ⁣}={f,g}−ℏ224(⋯ )+O(ℏ4)\{\!\{f,g\}\!\}=\{f,g\}-\frac{\hbar^2}{24}(\cdots)+\mathcal{O}(\hbar^4) で、von Neumann 方程式と Heisenberg 方程式はMoyal 方程式 ∂tW={ ⁣{H,W} ⁣}\partial_tW=\{\!\{H,W\}\!\} と ∂tA={ ⁣{A,H} ⁣}\partial_tA=\{\!\{A,H\}\!\} に両向きに同値に写る。

ℏ→0\hbar\to0 で Moyal 括弧は Poisson 括弧に、Moyal 方程式は Liouville 方程式に戻り、量子力学の2つの描像が古典力学の2つの書き方(Hamilton・Liouville)に並ぶ。2次以下のハミルトニアンでは時間発展は古典と厳密に一致し、共通例題③の波束の剪断もこの形で出る。それ以外では ℏ2V′′′∂p3W\hbar^2V'''\partial_p^3W の補正が残り、V=x4/4V=x^4/4 では古典の集団が位相のばらけで減衰したまま戻らないのに対し、量子では ⟨x⟩\langle x \rangle が復活する。固有状態は H⋆W=EWH\star W=EW、純粋状態は W⋆W=W/(2πℏ)W\star W=W/(2\pi\hbar)、熱平衡は Bloch 方程式 ∂βΩ=−H⋆Ω\partial_\beta\Omega=-H\star\Omega で書け、調和振動子の熱平衡はこの方程式から共通例題①の ZZ と WW の形で厳密に出る。一般の VV では熱平衡の記号の ℏ2\hbar^2 補正(Wigner–Kirkwood)が実数で ℏ\hbar の偶数次だけになる。負の値を取る WW が古典の確率と矛盾しない理由は 第5部第4章、共通例題①②の Wigner 関数は 第5部第5章で扱う。

総合問題 4次の非調和振動子を位相空間で読む ★★★

オリジナルH=p2/(2m)+λx4H=p^2/(2m)+\lambda x^4(λ>0\lambda>0)の1粒子を、位相空間の言葉で調べる。

  1. (1) (使う道具:この章の5.3.3節)Moyal 方程式を書き下し、右辺が ℏ2\hbar^2 の項で厳密に止まることを示せ。
  2. (2) (使う道具:5.3.3節・第2部第3章の Heisenberg 方程式)ddt⟨x^⟩\frac{\mathrm{d} }{\mathrm{d} t}\langle \hat x \rangle、ddt⟨p^⟩\frac{\mathrm{d} }{\mathrm{d} t}\langle \hat p \rangle、ddt⟨p^3⟩\frac{\mathrm{d} }{\mathrm{d} t}\langle \hat p^3 \rangle を Moyal 括弧で求め、古典の Liouville 方程式から出る式と比べよ。
  3. (3) (使う道具:5.3.4節・第5部第2章の Gauss 関数の WW)ℏ=m=1\hbar=m=1 とし、W=1πe−(x−x0)2−p2W=\frac1\pi\mathrm{e}^{-(x-x_0)^2-p^2} のときの ddt⟨p^3⟩\frac{\mathrm{d} }{\mathrm{d} t}\langle \hat p^3 \rangle の量子の補正と古典の項の比を x0x_0 で表し、x0=1x_0=1 と x0=3x_0=3 で数値を出せ。
  4. (4) (使う道具:5.3.5節・第5部第1章の Wigner–Kirkwood)ℏ=m=λ=1\hbar=m=\lambda=1、β=0.5\beta=0.5 で、熱平衡の記号 Ω=e−βH(1+ℏ2χ)\Omega=\mathrm{e}^{-\beta H}(1+\hbar^2\chi) から ZZ を数値で求め、古典の ZclZ_{\mathrm{cl}} と厳密な値 1.146771.14677 と比べよ。
  5. (5) (使う道具:5.3.6節の例題 5.3.6・第5部第1章の基底エネルギー)固有状態で ⟨T⟩=2⟨V⟩\langle T \rangle =2\langle V \rangle を示し、ℏ=m=λ=1\hbar=m=\lambda=1 の基底エネルギー E0=0.66799E_0=0.66799 から ⟨V⟩\langle V \rangle と ⟨T⟩\langle T \rangle を求めよ。
解答

(1) (5.3.11) で V=λx4V=\lambda x^4 とすると V′=4λx3V'=4\lambda x^3、V′′′=24λxV'''=24\lambda x、V(5)=0V^{(5)}=0 なので

∂W∂t=−pm∂W∂x+4λx3∂W∂p−λℏ2x∂3W∂p3\frac{\partial W}{\partial t}=-\frac pm\frac{\partial W}{\partial x}+4\lambda x^3\frac{\partial W}{\partial p}-\lambda\hbar^2x\frac{\partial^3W}{\partial p^3}

で厳密に止まる。ℏ4\hbar^4 以降の項はすべて VV の5階以上の微分を含み、4次の多項式では 0 になる。

(2) 定理 5.3.5 (b) の後の式 ddt⟨A^⟩=∫W{ ⁣{AW,H} ⁣}\frac{\mathrm{d} }{\mathrm{d} t}\langle \hat A \rangle =\int W\{\!\{A_W,H\}\!\} を使う。{ ⁣{x,H} ⁣}=p/m\{\!\{x,H\}\!\}=p/m、{ ⁣{p,H} ⁣}=−4λx3\{\!\{p,H\}\!\}=-4\lambda x^3 は (5.3.9) の補正の項が消えるので

ddt⟨x^⟩=⟨p^⟩m,ddt⟨p^⟩=−4λ∫Wx3=−4λ⟨x^3⟩.\frac{\mathrm{d} }{\mathrm{d} t}\langle \hat x \rangle =\frac{\langle \hat p \rangle }{m},\qquad \frac{\mathrm{d} }{\mathrm{d} t}\langle \hat p \rangle =-4\lambda\int Wx^3=-4\lambda\langle \hat x^3 \rangle .

古典の集団でも同じ形 ddt⟨x⟩=⟨p⟩/m\frac{\mathrm{d} }{\mathrm{d} t}\langle x \rangle =\langle p \rangle /m、ddt⟨p⟩=−4λ⟨x3⟩\frac{\mathrm{d} }{\mathrm{d} t}\langle p \rangle =-4\lambda\langle x^3 \rangle になる。p3p^3 では ℏ2\hbar^2 の項のうち −fpppgxxx-f_{ppp}g_{xxx} が残り、

{ ⁣{p3,H} ⁣}=−12λx3p2+ℏ224⋅6⋅24λx=−12λx3p2+6λℏ2x,ddt⟨p^3⟩=−12λ∫Wx3p2+6λℏ2⟨x^⟩.\begin{gathered} \{\!\{p^3,H\}\!\}=-12\lambda x^3p^2+\frac{\hbar^2}{24}\cdot6\cdot24\lambda x=-12\lambda x^3p^2+6\lambda\hbar^2x ,\\ \frac{\mathrm{d} }{\mathrm{d} t}\langle \hat p^3 \rangle =-12\lambda\int Wx^3p^2+6\lambda\hbar^2\langle \hat x \rangle . \end{gathered}

古典では −12λ⟨x3p2⟩-12\lambda\langle x^3p^2 \rangle だけで、量子では 6λℏ2⟨x^⟩6\lambda\hbar^2\langle \hat x \rangle が足される。例題 5.3.3 と同じく、違いは3次のモーメントから現れる。

(3) この WW では xx と pp が独立で、xx は平均 x0x_0・分散 1/21/2、pp は平均 0・分散 1/21/2。⟨x3⟩=x03+3x0⋅12\langle x^3 \rangle =x_0^3+3x_0\cdot\frac12 なので古典の項は −12λ⟨p2⟩⟨x3⟩=−6λ(x03+32x0)-12\lambda\langle p^2 \rangle \langle x^3 \rangle =-6\lambda\bigl(x_0^3+\tfrac32x_0\bigr)、量子の補正は 6λx06\lambda x_0。比の大きさは

6λx06λ(x03+32x0)=1x02+32\frac{6\lambda x_0}{6\lambda(x_0^3+\tfrac32x_0)}=\frac{1}{x_0^2+\tfrac32}

で、x0=1x_0=1 で 0.400.40、x0=3x_0=3 で 0.0950.095 である。x02x_0^2 は振動の作用(ℏ\hbar を単位とした大きさ)の目安なので、作用が大きいほど補正は小さく、5.3.4節の見積もり (ℏ/(LΔp))2(\hbar/(L\Delta p))^2 と同じ向きになる。

(4) (5.3.16) で V′′=12x2V''=12x^2、V′2=16x6V'^2=16x^6。pp で積分すると χ→−β2V′′/12+β3V′2/24=−β2x2+23β3x6\chi\to-\beta^2V''/12+\beta^3V'^2/24=-\beta^2x^2+\tfrac23\beta^3x^6 なので

Z=12π2πβ∫e−βx4(1−β2x2+23β3x6)dx.Z=\frac{1}{2\pi}\sqrt{\frac{2\pi}{\beta}}\int\mathrm{e}^{-\beta x^4}\Bigl(1-\beta^2x^2+\frac23\beta^3x^6\Bigr)\mathrm{d} x .

β=0.5\beta=0.5 で数値積分すると Zcl=1.21628Z_{\mathrm{cl}}=1.21628(括弧の 1 だけ)、補正つきで Z=1.14361Z=1.14361。厳密な値 1.146771.14677 に対し、古典は +6.1%+6.1\%、補正つきは −0.28%-0.28\% のずれで、第5部第1章の表と同じ値になる。Kirkwood の位置表示から出した式と、Bloch 方程式の Weyl 記号から出した式が同じ ZZ を与えることの数値での確認である。

(5) 例題 5.3.6 から定常状態で 2⟨T⟩=⟨xV′⟩=4λ⟨x4⟩=4⟨V⟩2\langle T \rangle =\langle xV' \rangle =4\lambda\langle x^4 \rangle =4\langle V \rangle、よって ⟨T⟩=2⟨V⟩\langle T \rangle =2\langle V \rangle。E0=⟨T⟩+⟨V⟩=3⟨V⟩E_0=\langle T \rangle +\langle V \rangle =3\langle V \rangle から ⟨V⟩=0.22266\langle V \rangle =0.22266、⟨T⟩=0.44533\langle T \rangle =0.44533。数値の対角化で求めた基底状態の ⟨x4⟩\langle x^4 \rangle がこの値に一致することを検算で確かめた。古典の周期運動でも時間平均で T‾=2V‾\overline T=2\overline V が成り立ち、ビリアル定理は2次以下でない HH でも括弧の形が古典と一致する例である(xpxp が2次なので)。

木下の勉強ステータス

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

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

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