tecktom.com

密度演算子による構成の要請

2026年10月17日10分で読めますAI生成教科書『虎の巻』第4部 第1章 なぜ状態ベクトルでは足りないか
目次
  1. 5つの要請
  2. von Neumann 方程式
  3. 状態の集合は凸集合で、純粋状態はその端点
  4. 木下の層

AI 生成・木下が未勉強

量子力学の様々な構成方法 第4部 第1章 §4.1.3 ★★★ 第4部の公理:状態=密度演算子・確率=Tr⁡(ρ^P^a)\operatorname{Tr}(\hat\rho\hat P_a)・時間発展=U^ρ^U^†\hat U\hat\rho\hat U^\dagger(von Neumann 方程式)・測定後=P^aρ^P^a\hat P_a\hat\rho\hat P_a・合成系=テンソル積と部分トレース

5つの要請

第4.1.2節で、期待値と確率が ρ^\hat\rho だけで書けることを見た。第4部では、これを出発点の要請として最初に置く。基準の公理(第1部第2章)の各項目を ρ^\hat\rho で書き直すと次のようになる。

定義 4.1.3(密度演算子による構成の要請 ★★★)

  1. (D1) 状態:系の状態は、ヒルベルト空間 H\mathcal{H} 上の密度演算子 ρ^\hat\rho(定義 4.1.1)で表す。
  2. (D2) 測定の確率:物理量 A^=∑aaP^a\hat A=\sum_aa\hat P_a を測って値 aa を得る確率は Pr⁡(a)=Tr⁡(ρ^P^a)\Pr(a)=\operatorname{Tr}(\hat\rho\hat P_a)。したがって期待値は Tr⁡(ρ^A^)\operatorname{Tr}(\hat\rho\hat A)。
  3. (D3) 時間発展:孤立した系では ρ^(t)=U^(t)ρ^(0)U^(t)†\hat\rho(t)=\hat U(t)\hat\rho(0)\hat U(t)^\dagger、U^(t)=e−iH^t/ℏ\hat U(t)=\mathrm{e}^{-\mathrm{i}\hat Ht/\hbar}(H^\hat H が時刻によらないとき)。
  4. (D4) 測定後の状態:値 aa を得た後の状態は P^aρ^P^a/Tr⁡(ρ^P^a)\hat P_a\hat\rho\hat P_a/\operatorname{Tr}(\hat\rho\hat P_a)(Lüders の規則)。結果を読まずに平均すると ∑aP^aρ^P^a\sum_a\hat P_a\hat\rho\hat P_a。
  5. (D5) 合成系:系 A と B を合わせた系の状態は HA⊗HB\mathcal{H}_{\mathrm A}\otimes\mathcal{H}_{\mathrm B} 上の密度演算子 ρ^AB\hat\rho_{\mathrm{AB}}。A だけを見るときの状態は部分トレース ρ^A=Tr⁡Bρ^AB\hat\rho_{\mathrm A}=\operatorname{Tr}_{\mathrm B}\hat\rho_{\mathrm{AB}}(第4.1.5節)。

(D1)〜(D5) はどれも、アンサンブル {pi,∣ψi⟩}\{p_i,| \psi_i \rangle \} の各成分に基準の公理を当てて平均したものと一致するように選んである。(D2) は (4.1.5) そのものである。(D3) は各 ∣ψi⟩| \psi_i \rangle を U^∣ψi⟩\hat U| \psi_i \rangle に送れば∑ipiU^∣ψi⟩⟨ψi∣U^†=U^ρ^U^†\sum_ip_i\hat U| \psi_i \rangle \langle \psi_i |\hat U^\dagger=\hat U\hat\rho\hat U^\dagger となることから従う。(D4) で結果を読まない場合は、値 aa が確率 Tr⁡(ρ^P^a)\operatorname{Tr}(\hat\rho\hat P_a) で出て、その後の状態が P^aρ^P^a/Tr⁡(ρ^P^a)\hat P_a\hat\rho\hat P_a/\operatorname{Tr}(\hat\rho\hat P_a) になるアンサンブルを平均したものである。

von Neumann 方程式

(D3) を時間で微分すると、dU^dt=−iℏH^U^\frac{\mathrm{d} \hat U}{\mathrm{d} t}=-\tfrac{\mathrm{i}}{\hbar}\hat H\hat U、dU^†dt=iℏU^†H^\frac{\mathrm{d} \hat U^\dagger}{\mathrm{d} t}=\tfrac{\mathrm{i}}{\hbar}\hat U^\dagger\hat H から

iℏdρ^dt=H^ρ^−ρ^H^=[H^, ρ^](4.1.7)\mathrm{i}\hbar\frac{\mathrm{d} \hat\rho}{\mathrm{d} t}=\hat H\hat\rho-\hat\rho\hat H=\left[\hat H,\,\hat\rho\right]\tag{4.1.7}

を得る。これを von Neumann 方程式と呼ぶ。H^\hat H が時刻による場合も、各 ∣ψi(t)⟩| \psi_i(t) \rangle が Schrödinger 方程式に従うことから同じ式が出る。Tr⁡[H^, ρ^]=0\operatorname{Tr}\left[\hat H,\,\hat\rho\right]=0 なので Tr⁡ρ^\operatorname{Tr}\hat\rho は保たれ、Tr⁡ρ^2\operatorname{Tr}\hat\rho^2 も ddtTr⁡ρ^2=2iℏTr⁡(ρ^[H^, ρ^])=0\frac{\mathrm{d} }{\mathrm{d} t}\operatorname{Tr}\hat\rho^2=\tfrac{2}{\mathrm{i}\hbar}\operatorname{Tr}(\hat\rho\left[\hat H,\,\hat\rho\right])=0 で保たれる。孤立した系の時間発展では、純粋状態は純粋状態のまま、混合状態は同じ純粋度の混合状態のままである。純粋度は、(D4) の測定をした場合か、系が環境と相互作用して (D5) の部分トレースを取る場合にだけ変わる。後者は第4部第5章(第4部第5章)で Lindblad 方程式として扱う。

符号に注意しておく。Heisenberg 描像の演算子は iℏdA^Hdt=[A^H, H^]\mathrm{i}\hbar\frac{\mathrm{d} \hat A_{\mathrm H}}{\mathrm{d} t}=\left[\hat A_{\mathrm H},\,\hat H\right] に従う(第2部第3章)。(4.1.7) は交換子の順序が逆である。ρ^\hat\rho は Schrödinger 描像の状態で、物理量ではないので、Heisenberg の運動方程式に従う理由はない。両者が同じ期待値を与えることは第4.1.6節で確かめる。

状態の集合は凸集合で、純粋状態はその端点

2つの密度演算子 ρ^1,ρ^2\hat\rho_1,\hat\rho_2 と 0≤p≤10\le p\le1 に対して、pρ^1+(1−p)ρ^2p\hat\rho_1+(1-p)\hat\rho_2 もまた (4.1.6) を満たす。これは「確率 pp で ρ^1\hat\rho_1 を、確率 1−p1-p で ρ^2\hat\rho_2 を準備する」という操作に当たる。状態ベクトルでは p∣ψ1⟩+1−p∣ψ2⟩\sqrt p| \psi_1 \rangle +\sqrt{1-p}| \psi_2 \rangle のような重ね合わせは作れても、確率的に混ぜる操作を状態の中で書けなかった。(D1) の状態の集合は、確率的な混合で閉じている(凸集合である)。

命題 4.1.4(純粋状態は混ぜて作れない)

∣ψ⟩⟨ψ∣=pρ^1+(1−p)ρ^2| \psi \rangle \langle \psi |=p\hat\rho_1+(1-p)\hat\rho_2(0<p<10<p<1)なら ρ^1=ρ^2=∣ψ⟩⟨ψ∣\hat\rho_1=\hat\rho_2=| \psi \rangle \langle \psi | である。つまり純粋状態は状態の集合の端点である。

証明

∣ψ⟩| \psi \rangle に直交する任意の ∣χ⟩| \chi \rangle で挟むと 0=p⟨χ∣ρ^1∣χ⟩+(1−p)⟨χ∣ρ^2∣χ⟩0=p\langle \chi | \hat\rho_1 | \chi \rangle +(1-p)\langle \chi | \hat\rho_2 | \chi \rangle。両項は非負なので ⟨χ∣ρ^j∣χ⟩=0\langle \chi | \hat\rho_j | \chi \rangle =0。半正定値の演算子で ⟨χ∣ρ^j∣χ⟩=0\langle \chi | \hat\rho_j | \chi \rangle =0 なら ρ^j1/2∣χ⟩=0\hat\rho_j^{1/2}| \chi \rangle =0、したがって ρ^j∣χ⟩=0\hat\rho_j| \chi \rangle =0 である。ρ^j\hat\rho_j は ∣ψ⟩| \psi \rangle の直交補空間を 00 に送るエルミート演算子なので ρ^j=cj∣ψ⟩⟨ψ∣\hat\rho_j=c_j| \psi \rangle \langle \psi | の形で、トレース1から cj=1c_j=1。

混合状態は端点ではないので、2つ以上の異なる状態の混合として書ける。書き方が1通りでないことが、第4部第3章の主題につながる。

例題 4.1.3 結果を読まない測定は純粋状態を混合状態にする ★★★

オリジナル量子ビットを ∣+⟩| + \rangle に準備し、σz\sigma_z を測る。P^↑=∣↑⟩⟨↑∣\hat P_\uparrow=| \uparrow \rangle \langle \uparrow |、P^↓=∣↓⟩⟨↓∣\hat P_\downarrow=| \downarrow \rangle \langle \downarrow | とする。

  1. (1) 結果が +1+1 である確率と、そのときの測定後の状態を (D2)・(D4) で求めよ。
  2. (2) 結果を読まずに平均した測定後の状態を求め、純粋度が 11 から 12\tfrac12 に下がることを確かめよ。
  3. (3) 測定の前後で ⟨σx⟩\langle \sigma_x \rangle はどう変わるか。状態ベクトルの公理だけでは (2) の状態が書けないことと合わせて述べよ。
解答

(1) ρ^=∣+⟩⟨+∣=12(1111)\hat\rho=| + \rangle \langle + |=\tfrac12\begin{pmatrix}1&1\\1&1\end{pmatrix}。Pr⁡(+1)=Tr⁡(ρ^P^↑)=⟨↑∣ρ^∣↑⟩=12\Pr(+1)=\operatorname{Tr}(\hat\rho\hat P_\uparrow)=\langle \uparrow | \hat\rho | \uparrow \rangle =\tfrac12。P^↑ρ^P^↑=12∣↑⟩⟨↑∣\hat P_\uparrow\hat\rho\hat P_\uparrow=\tfrac12| \uparrow \rangle \langle \uparrow | なので、測定後の状態は ∣↑⟩⟨↑∣| \uparrow \rangle \langle \uparrow |。−1-1 も確率 12\tfrac12 で、その後は ∣↓⟩⟨↓∣| \downarrow \rangle \langle \downarrow |。

(2) P^↑ρ^P^↑+P^↓ρ^P^↓=12(1000)+12(0001)=121\hat P_\uparrow\hat\rho\hat P_\uparrow+\hat P_\downarrow\hat\rho\hat P_\downarrow=\tfrac12\begin{pmatrix}1&0\\0&0\end{pmatrix}+\tfrac12\begin{pmatrix}0&0\\0&1\end{pmatrix}=\tfrac12\mathbf 1。Tr⁡ρ^2\operatorname{Tr}\hat\rho^2 は測定前が 11、測定後が Tr⁡(141)=12\operatorname{Tr}(\tfrac14\mathbf 1)=\tfrac12。非対角成分(∣↑⟩| \uparrow \rangle と ∣↓⟩| \downarrow \rangle の位相の関係)が 00 になり、対角成分だけが残った。

(3) 測定前は Tr⁡(ρ^σx)=1\operatorname{Tr}(\hat\rho\sigma_x)=1、測定後は Tr⁡(12σx)=0\operatorname{Tr}(\tfrac12\sigma_x)=0。σz\sigma_z の測定は σx\sigma_x の値を壊す。測定装置の結果を読まない人にとって、測定後の系は ∣↑⟩| \uparrow \rangle か ∣↓⟩| \downarrow \rangle のどちらかが半々で入ったアンサンブルで、例題 4.1.1 (2) と同じ理由でどの状態ベクトルでも書けない。結果を読まない測定は、場面1に並ぶ4つ目の場面である。

木下の勉強ステータス

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

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

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