tecktom.com

橋:期待値を各構成の言葉で書く

2026年10月17日15分で読めますAI生成教科書『虎の巻』第1部 第1章 この巻で扱う構成方法とその長所と短所
目次
  1. 期待値の書き方の一覧
  2. 状態ベクトルと密度行列:純粋状態の場合
  3. 密度行列と位相空間:Wigner 関数の期待値
  4. 学部の授業の形との対応
  5. 木下の層

AI 生成・木下が未勉強

量子力学の様々な構成方法 第1部 第1章 §1.1.6 ★★★ 全部の橋は「同じ期待値を与える」ことで示すので、期待値の書き方を構成ごとに並べ、純粋状態と Wigner 関数の2本をこの章で証明する

期待値の書き方の一覧

1.1.1節で決めたとおり、2つの構成が同値であるとは同じ期待値を与えることである。そこで、基準の期待値 (1.1.1) を各構成の言葉で書いた形を並べる。各行の記号は、それぞれの部で定義する。

表 1.1.5 同じ期待値 ⟨A^⟩\langle \hat A \rangle の構成ごとの書き方
構成期待値の形示す所
波動関数(学部の授業の形)∫ψ(x)∗ A(x,−iℏddx)ψ(x) dx\int\psi(x)^*\,A\bigl(x,-\mathrm{i}\hbar\tfrac{\mathrm{d}}{\mathrm{d} x}\bigr)\psi(x)\,\mathrm{d} x第1部第3章
状態ベクトル(基準)⟨ψ∣A^∣ψ⟩\langle \psi | \hat A | \psi \rangle第1部第2章
行列力学∑m,ncm∗Amn(t)cn\sum_{m,n}c_m^*A_{mn}(t)c_n(Amn(t)=⟨m∣A^H(t)∣n⟩A_{mn}(t)=\langle m | \hat A_H(t) | n \rangle)第2部第6章
密度行列Tr⁡(ρ^A^)\operatorname{Tr}(\hat\rho\hat A)この節(純粋状態)・第4部第2章
位相空間∫W(x,p) AW(x,p) dx dp\int W(x,p)\,A_W(x,p)\,\mathrm{d} x\,\mathrm{d} pこの節・第5部第2章
経路積分1Z∫Dx A(x(τ)) e−SE[x]/ℏ\frac1Z\int\mathcal Dx\,A(x(\tau))\,\mathrm{e}^{-S_{\mathrm E}[x]/\hbar}(熱平衡・虚時間)第3部第6章
第二量子化⟨Ψ∣A^∣Ψ⟩\langle \Psi | \hat{\mathcal A} | \Psi \rangle(A^\hat{\mathcal A} は a^†\hat a^\dagger, a^\hat a で書いた物理量)第6部第2章
流体∫ρ(x)[12m∣v∣2+Q+V]dx\int\rho(x)\bigl[\tfrac12m\lvert \bm{v}\rvert ^2+Q+V\bigr]\mathrm{d} x など、密度で重みをつけた積分第7部第1章
代数状態 ω\omega の値 ω(A)\omega(A)第8部第2章
再構成状態空間の点と効果(測定の結果)を組にした確率第9部第4章

流体の行はエネルギーの期待値の例で、QQ は第7部第1章で定義する量子ポテンシャルである。表の多くの行は各部で示すが、そのうち2本は短いので、この節で両向きを証明する。

状態ベクトルと密度行列:純粋状態の場合

命題 1.1.1(純粋状態の対応 ★★★)

ヒルベルト空間 H\mathcal{H} の単位ベクトル ∣ψ⟩| \psi \rangle に ρ^=∣ψ⟩⟨ψ∣\hat\rho=| \psi \rangle \langle \psi | を対応させると、すべての有界な物理量 A^\hat A について Tr⁡(ρ^A^)=⟨ψ∣A^∣ψ⟩\operatorname{Tr}(\hat\rho\hat A)=\langle \psi | \hat A | \psi \rangle が成り立ち、ρ^\hat\rho は ρ^≥0\hat\rho\ge0、Tr⁡ρ^=1\operatorname{Tr}\hat\rho=1、ρ^2=ρ^\hat\rho^2=\hat\rho を満たす。逆に、ρ^≥0\hat\rho\ge0、Tr⁡ρ^=1\operatorname{Tr}\hat\rho=1、ρ^2=ρ^\hat\rho^2=\hat\rho を満たす ρ^\hat\rho は、位相因子 eiθ\mathrm{e}^{\mathrm{i}\theta} を除いて一意な単位ベクトル ∣ψ⟩| \psi \rangle によって ρ^=∣ψ⟩⟨ψ∣\hat\rho=| \psi \rangle \langle \psi | と書ける。

証明

(基準→密度行列)正規直交基底 {∣k⟩}\{| k \rangle \} でトレースをとると Tr⁡(∣ψ⟩⟨ψ∣A^)=∑k⟨k∣ψ⟩⟨ψ∣A^∣k⟩=∑k⟨ψ∣A^∣k⟩⟨k∣ψ⟩\operatorname{Tr}(| \psi \rangle \langle \psi |\hat A)=\sum_k\langle k | \psi \rangle \langle \psi | \hat A | k \rangle =\sum_k\langle \psi | \hat A | k \rangle \langle k | \psi \rangle となり、∑k∣k⟩⟨k∣=1\sum_k| k \rangle \langle k |=1 を使えば ⟨ψ∣A^∣ψ⟩\langle \psi | \hat A | \psi \rangle に等しい。⟨ϕ∣ρ^∣ϕ⟩=∣⟨ψ∣ϕ⟩∣2≥0\langle \phi | \hat\rho | \phi \rangle =\lvert \langle \psi | \phi \rangle \rvert ^2\ge0、Tr⁡ρ^=⟨ψ∣ψ⟩=1\operatorname{Tr}\hat\rho=\langle \psi | \psi \rangle =1、ρ^2=∣ψ⟩⟨ψ∣ψ⟩⟨ψ∣=ρ^\hat\rho^2=| \psi \rangle \langle \psi | \psi \rangle \langle \psi |=\hat\rho である。

(密度行列→基準)ρ^2=ρ^\hat\rho^2=\hat\rho と ρ^=ρ^†\hat\rho=\hat\rho^\dagger から ρ^\hat\rho は直交射影で、固有値は 00 か 11 である。Tr⁡ρ^=1\operatorname{Tr}\hat\rho=1 なので固有値 11 はちょうど1つで、その固有空間は1次元である。その固有空間の単位ベクトルを ∣ψ⟩| \psi \rangle とすると ρ^=∣ψ⟩⟨ψ∣\hat\rho=| \psi \rangle \langle \psi |。固有空間の単位ベクトルは eiθ∣ψ⟩\mathrm{e}^{\mathrm{i}\theta}| \psi \rangle の形に限られ、∣ψ⟩| \psi \rangle を eiθ∣ψ⟩\mathrm{e}^{\mathrm{i}\theta}| \psi \rangle に替えても ∣ψ⟩⟨ψ∣| \psi \rangle \langle \psi | は変わらない。

命題 1.1.1 は、基準の状態ベクトルが「位相を除いて」密度行列と同じものであることを言っている。全体の位相 eiθ\mathrm{e}^{\mathrm{i}\theta} はどの期待値にも現れないので、状態ベクトルの構成が持つ情報のうち予言に関係する部分は、純粋状態の密度行列がすべて持っている。混合状態(ρ^2≠ρ^\hat\rho^2\neq\hat\rho)への拡張と、その逆向きの浄化は第4部第2章で扱う(相対エントロピーから情報と物理を組み立てる 第3部第1章 も密度行列と量子チャネルを扱う)。

密度行列と位相空間:Wigner 関数の期待値

第5部で使う Wigner 関数と Weyl 記号を、この巻の約束どおり次の形で書く。

W(x,p)=1πℏ∫dy ⟨x+y∣ρ^∣x−y⟩ e−2ipy/ℏ,AW(x,p)=∫dy ⟨x−y2∣A^∣x+y2⟩ eipy/ℏ.(1.1.9)W(x,p)=\frac{1}{\pi\hbar}\int\mathrm{d} y\,\langle x+y | \hat\rho | x-y \rangle \,\mathrm{e}^{-2\mathrm{i} py/\hbar},\qquad A_W(x,p)=\int\mathrm{d} y\,\langle x-\tfrac y2 | \hat A | x+\tfrac y2 \rangle \,\mathrm{e}^{\mathrm{i} py/\hbar}.\tag{1.1.9}

命題 1.1.2(Wigner 関数による期待値 ★★★)

式 (1.1.9) の WW と AWA_W について(積分の順序を交換してよい ρ^\hat\rho, A^\hat A に対して)、∫dx dp  W(x,p) AW(x,p)=Tr⁡(ρ^A^)\displaystyle\int\mathrm{d} x\,\mathrm{d} p\;W(x,p)\,A_W(x,p)=\operatorname{Tr}(\hat\rho\hat A) が成り立つ。

証明

左辺に式 (1.1.9) を代入すると、yy と y′y' の2つの積分変数を持つ式になる。

∫dx dp W ⁣AW=1πℏ∫dx dy dy′ ⟨x+y∣ρ^∣x−y⟩⟨x−y′2∣A^∣x+y′2⟩∫dp eip(y′−2y)/ℏ.\int\mathrm{d} x\,\mathrm{d} p\,W\!A_W=\frac{1}{\pi\hbar}\int\mathrm{d} x\,\mathrm{d} y\,\mathrm{d} y'\,\langle x+y | \hat\rho | x-y \rangle \langle x-\tfrac{y'}2 | \hat A | x+\tfrac{y'}2 \rangle \int\mathrm{d} p\,\mathrm{e}^{\mathrm{i} p(y'-2y)/\hbar}.

pp の積分は ∫dp eip(y′−2y)/ℏ=2πℏ δ(y′−2y)\int\mathrm{d} p\,\mathrm{e}^{\mathrm{i} p(y'-2y)/\hbar}=2\pi\hbar\,\delta(y'-2y) なので、y′y' の積分で y′=2yy'=2y となり、

∫dx dp W ⁣AW=2∫dx dy ⟨x+y∣ρ^∣x−y⟩⟨x−y∣A^∣x+y⟩.\int\mathrm{d} x\,\mathrm{d} p\,W\!A_W=2\int\mathrm{d} x\,\mathrm{d} y\,\langle x+y | \hat\rho | x-y \rangle \langle x-y | \hat A | x+y \rangle .

変数を u=x+yu=x+y、v=x−yv=x-y に変えると dx dy=12du dv\mathrm{d} x\,\mathrm{d} y=\frac12\mathrm{d} u\,\mathrm{d} v なので、右辺は ∫du dv ⟨u∣ρ^∣v⟩⟨v∣A^∣u⟩=Tr⁡(ρ^A^)\int\mathrm{d} u\,\mathrm{d} v\,\langle u | \hat\rho | v \rangle \langle v | \hat A | u \rangle =\operatorname{Tr}(\hat\rho\hat A) になる。最後の等号は位置の完全系 ∫dv∣v⟩⟨v∣=1\int\mathrm{d} v| v \rangle \langle v |=1 を使った。

逆向きは、WW から ρ^\hat\rho の行列要素を取り戻す式で示す。式 (1.1.9) の WW を pp について Fourier 逆変換すると、

⟨x+y∣ρ^∣x−y⟩=∫dp W(x,p) e2ipy/ℏ\langle x+y | \hat\rho | x-y \rangle =\int\mathrm{d} p\,W(x,p)\,\mathrm{e}^{2\mathrm{i} py/\hbar}

となる(∫dp e2ip(y−y′′)/ℏ=πℏ δ(y−y′′)\int\mathrm{d} p\,\mathrm{e}^{2\mathrm{i} p(y-y'')/\hbar}=\pi\hbar\,\delta(y-y'') を使う)。左辺の xx と yy はすべての ρ^\hat\rho の行列要素 ⟨u∣ρ^∣v⟩\langle u | \hat\rho | v \rangle を尽くすので、WW と ρ^\hat\rho は互いに一意に決まる。これで密度行列から位相空間へ、位相空間から密度行列への両向きの対応がそろった。⋆\star積と Moyal 方程式による動力学の対応は第5部第3章で示す。

学部の授業の形との対応

学部の量子力学の授業では、期待値を波動関数で ⟨A⟩=∫ψ∗Aψ dx\langle A \rangle =\int\psi^*A\psi\,\mathrm{d} x と書き、運動量は −iℏ d/dx-\mathrm{i}\hbar\,\mathrm{d}/\mathrm{d} x に置き換える。これは表 1.1.5 の1行目で、基準の ⟨ψ∣A^∣ψ⟩\langle \psi | \hat A | \psi \rangle に位置の完全系を2回挿入したものにあたる(第1部第3章)。時間発展を Schrödinger 方程式で計算するのは Schrödinger 描像、エネルギー固有状態の基底で行列要素を計算するのは行列力学の書き方、統計力学の授業で Z=∑ne−βEnZ=\sum_n\mathrm{e}^{-\beta E_n} と書くのは密度行列 e−βH^\mathrm{e}^{-\beta\hat H} のトレースである。学部の授業は、同じ1つの問題の中でこれらの構成を断らずに使い分けている。この巻はその使い分けを、どこで同じになるかを示したうえで行う。

例題 1.1.6 基底状態と熱平衡の位置のゆらぎを4つの構成で出す ★★★

オリジナル共通例題①の調和振動子について、⟨x^2⟩\langle \hat x^2 \rangle を求める。

  1. (1) 波動関数 ψ0(x)=(mω/πℏ)1/4e−mωx2/(2ℏ)\psi_0(x)=(m\omega/\pi\hbar)^{1/4}\mathrm{e}^{-m\omega x^2/(2\hbar)} で ∫ψ02x2 dx\int\psi_0^2x^2\,\mathrm{d} x を計算せよ。
  2. (2) 行列力学の x^=ℏ/(2mω)(a^+a^†)\hat x=\sqrt{\hbar/(2m\omega)}(\hat a+\hat a^\dagger) を使って ⟨0∣x^2∣0⟩\langle 0 | \hat x^2 | 0 \rangle を計算せよ。
  3. (3) 基底状態の Wigner 関数 W0=1πℏe−2H/(ℏω)W_0=\frac{1}{\pi\hbar}\mathrm{e}^{-2H/(\hbar\omega)}(H=p2/(2m)+mω2x2/2H=p^2/(2m)+m\omega^2x^2/2)で ∫W0x2 dx dp\int W_0x^2\,\mathrm{d} x\,\mathrm{d} p を計算せよ(x2x^2 の Weyl 記号は x2x^2)。
  4. (4) 熱平衡の Wigner 関数 Wβ=tanh⁡(βℏω/2)πℏexp⁡(−2tanh⁡(βℏω/2)Hℏω)W_\beta=\frac{\tanh(\beta\hbar\omega/2)}{\pi\hbar}\exp\bigl(-\frac{2\tanh(\beta\hbar\omega/2)H}{\hbar\omega}\bigr) で ⟨x^2⟩\langle \hat x^2 \rangle を求め、密度行列で Tr⁡(ρ^x^2)=∑npn⟨n∣x^2∣n⟩\operatorname{Tr}(\hat\rho\hat x^2)=\sum_np_n\langle n | \hat x^2 | n \rangle(pn=e−βEn/Zp_n=\mathrm{e}^{-\beta E_n}/Z)と計算したものと一致することを示せ。
解答

(1) α=mω/ℏ\alpha=m\omega/\hbar とすると、Gauss 積分 ∫x2e−αx2dx=π2α3/2\int x^2\mathrm{e}^{-\alpha x^2}\mathrm{d} x=\frac{\sqrt\pi}{2\alpha^{3/2}} から∫ψ02x2 dx=απ⋅π2α3/2=12α=ℏ2mω\int\psi_0^2x^2\,\mathrm{d} x=\sqrt{\frac{\alpha}{\pi}}\cdot\frac{\sqrt\pi}{2\alpha^{3/2}}=\frac{1}{2\alpha}=\frac{\hbar}{2m\omega}。

(2) (a^+a^†)2=a^2+a^a^†+a^†a^+a^†2(\hat a+\hat a^\dagger)^2=\hat a^2+\hat a\hat a^\dagger+\hat a^\dagger\hat a+\hat a^{\dagger2} の ∣0⟩| 0 \rangle での期待値は ⟨0∣a^a^†∣0⟩=1\langle 0 | \hat a\hat a^\dagger | 0 \rangle =1 だけが残るので、⟨0∣x^2∣0⟩=ℏ2mω\langle 0 | \hat x^2 | 0 \rangle =\frac{\hbar}{2m\omega}。

(3) −2Hℏω=−mωℏx2−p2mωℏ-\frac{2H}{\hbar\omega}=-\frac{m\omega}{\hbar}x^2-\frac{p^2}{m\omega\hbar} なので W0W_0 は xx と pp の独立な Gauss 関数の積である。pp の積分は πmωℏ\sqrt{\pi m\omega\hbar}、xx の積分は π2(ℏ/mω)3/2\frac{\sqrt\pi}{2}(\hbar/m\omega)^{3/2} で、

∫W0x2 dx dp=1πℏπmωℏ⋅π2(ℏmω)3/2=ℏ2mω.\int W_0x^2\,\mathrm{d} x\,\mathrm{d} p=\frac{1}{\pi\hbar}\sqrt{\pi m\omega\hbar}\cdot\frac{\sqrt\pi}{2}\Bigl(\frac{\hbar}{m\omega}\Bigr)^{3/2}=\frac{\hbar}{2m\omega}.

同じ計算で ∫W0 dx dp=1\int W_0\,\mathrm{d} x\,\mathrm{d} p=1 も確かめられる。

(4) τ=tanh⁡(βℏω/2)\tau=\tanh(\beta\hbar\omega/2) とすると、WβW_\beta の xx についての指数は −τmωx2/ℏ-\tau m\omega x^2/\hbar である。WβW_\beta は xx と pp の独立な Gauss 関数の積で規格化されているので、xx の分散は ℏ2τmω=ℏ2mωcoth⁡βℏω2\frac{\hbar}{2\tau m\omega}=\frac{\hbar}{2m\omega}\coth\frac{\beta\hbar\omega}{2}。密度行列では ⟨n∣x^2∣n⟩=ℏ2mω(2n+1)\langle n | \hat x^2 | n \rangle =\frac{\hbar}{2m\omega}(2n+1) なので

Tr⁡(ρ^x^2)=ℏ2mω(2⟨n⟩+1)=ℏ2mω(2eβℏω−1+1)=ℏ2mωcoth⁡βℏω2\operatorname{Tr}(\hat\rho\hat x^2)=\frac{\hbar}{2m\omega}\bigl(2\langle n \rangle +1\bigr)=\frac{\hbar}{2m\omega}\Bigl(\frac{2}{\mathrm{e}^{\beta\hbar\omega}-1}+1\Bigr)=\frac{\hbar}{2m\omega}\coth\frac{\beta\hbar\omega}{2}

(⟨n⟩\langle n \rangle は例題 1.1.3)。2つの結果は一致し、命題 1.1.2 の具体例になっている。T→0T\to0 では coth⁡→1\coth\to1 で(1)〜(3)の基底状態の値に、高温では coth⁡(βℏω/2)≈2kBT/(ℏω)\coth(\beta\hbar\omega/2)\approx2 k_{\mathrm B} T/(\hbar\omega) で古典の等分配 mω2⟨x2⟩=kBTm\omega^2\langle x^2 \rangle = k_{\mathrm B} T に戻る。

木下の勉強ステータス

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

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

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