量子力学の様々な構成方法 第5部 第4章 §5.4.2 ★★★ なぜ観測と矛盾しないか。確率はすべて非負の作用素との重なりで、相空間の小さな領域を測る測定は存在しない
確率はすべて非負の作用素との重なり
基準の公理では、測定の結果 k k k が出る確率は、結果 k k k に対応する作用素 E ^ k \hat E_k E ^ k を使って Pr ( k ) = Tr ( ρ ^ E ^ k ) \Pr(k)=\operatorname{Tr}(\hat\rho\hat E_k) Pr ( k ) = Tr ( ρ ^ E ^ k ) で与えられる。射影測定なら E ^ k \hat E_k E ^ k はスペクトル分解の射影で(第1部第2章 の Q3)、一般の測定では0 ≤ E ^ k ≤ 1 0\le\hat E_k\le1 0 ≤ E ^ k ≤ 1 、∑ k E ^ k = 1 \sum_k\hat E_k=1 ∑ k E ^ k = 1 を満たす作用素の組になる(第4部第5章 の Kraus 表示)。E ^ k \hat E_k E ^ k の Weyl 記号を ( E k ) W (E_k)_{\mathrm W} ( E k ) W と書くと、第5部第2章 のトレースの公式から
Pr ( k ) = Tr ( ρ ^ E ^ k ) = ∫ W ( x , p ) ( E k ) W ( x , p ) d x d p (5.4.6) \Pr(k)=\operatorname{Tr}(\hat\rho\hat E_k)=\int W(x,p)\,(E_k)_{\mathrm W}(x,p)\,\mathrm{d} x\,\mathrm{d} p\tag{5.4.6} Pr ( k ) = Tr ( ρ ^ E ^ k ) = ∫ W ( x , p ) ( E k ) W ( x , p ) d x d p ( 5.4.6 )
となる。左辺は ρ ^ ≥ 0 \hat\rho\ge0 ρ ^ ≥ 0 と E ^ k ≥ 0 \hat E_k\ge0 E ^ k ≥ 0 から必ず 0 0 0 以上である。右辺の積分の中で W W W が負になっていても、重みの ( E k ) W (E_k)_{\mathrm W} ( E k ) W が「非負の作用素の記号」である限り、積分の結果は負にならない。観測にかかる量はすべてこの形をしているので、W W W の負の値が観測された確率として表に出ることはない。
古典の確率論では、事象 D D D (相空間の領域)の確率は ∫ D ρ c l \int_D\rho_{\mathrm{cl}} ∫ D ρ cl 、つまり重みが指示関数 χ D \chi_D χ D (D D D の中で1、外で0)のときの積分である。量子でも ∫ D W \int_DW ∫ D W が何かの確率になるには、χ D \chi_D χ D が非負の作用素の Weyl 記号でなければならない。次の例題が示すように、領域が小さいとそうならない。
例題 5.4.2 円板の指示関数を量子化した作用素 ★★★
オリジナル D a D_a D a を相空間の楕円の内側 { u < a } \{u<a\} { u < a } (u = 2 H / ( ℏ ω ) u=2H/(\hbar\omega) u = 2 H / ( ℏ ω ) 、面積 π ℏ a \pi\hbar a π ℏ a )とし、その指示関数 χ D a \chi_{D_a} χ D a を Weyl 量子化した作用素χ ^ a = Q W [ χ D a ] \hat\chi_a=\mathcal Q_{\mathrm W}[\chi_{D_a}] χ ^ a = Q W [ χ D a ] を考える(第5部第2章 )。
(1) χ ^ a \hat\chi_a χ ^ a は数状態 ∣ n ⟩ | n \rangle ∣ n ⟩ で対角になる。固有値が λ n ( a ) = ( − 1 ) n ∫ 0 a e − u L n ( 2 u ) d u \lambda_n(a)=(-1)^n\int_0^a\mathrm{e}^{-u}L_n(2u)\,\mathrm{d} u λ n ( a ) = ( − 1 ) n ∫ 0 a e − u L n ( 2 u ) d u であることを示せ。
(2) λ 0 ( a ) \lambda_0(a) λ 0 ( a ) と λ 1 ( a ) \lambda_1(a) λ 1 ( a ) を求め、a = 1 / 2 a=1/2 a = 1/2 での値を出せ。
(3) χ ^ a \hat\chi_a χ ^ a が非負の作用素でなくなる a a a の範囲を、λ 1 \lambda_1 λ 1 について求めよ。
解答 (1) χ D a \chi_{D_a} χ D a は H H H だけの関数なので、Weyl 量子化した作用素は相空間の回転(調和振動子の時間発展)で不変で、n ^ \hat n n ^ と交換する。n ^ \hat n n ^ の固有値は縮退していないので χ ^ a \hat\chi_a χ ^ a は ∣ n ⟩ | n \rangle ∣ n ⟩ で対角である。固有値は (5.4.6) と同じトレースの公式で
λ n ( a ) = ⟨ n ∣ χ ^ a ∣ n ⟩ = ∫ χ D a W n d x d p = π ℏ ∫ 0 a ( − 1 ) n π ℏ e − u L n ( 2 u ) d u = ( − 1 ) n ∫ 0 a e − u L n ( 2 u ) d u . \lambda_n(a)=\langle n | \hat\chi_a | n \rangle =\int\chi_{D_a}W_n\,\mathrm{d} x\,\mathrm{d} p=\pi\hbar\int_0^a\frac{(-1)^n}{\pi\hbar}\mathrm{e}^{-u}L_n(2u)\,\mathrm{d} u
=(-1)^n\int_0^a\mathrm{e}^{-u}L_n(2u)\,\mathrm{d} u . λ n ( a ) = ⟨ n ∣ χ ^ a ∣ n ⟩ = ∫ χ D a W n d x d p = π ℏ ∫ 0 a π ℏ ( − 1 ) n e − u L n ( 2 u ) d u = ( − 1 ) n ∫ 0 a e − u L n ( 2 u ) d u .
(2) L 0 = 1 L_0=1 L 0 = 1 から λ 0 ( a ) = 1 − e − a \lambda_0(a)=1-\mathrm{e}^{-a} λ 0 ( a ) = 1 − e − a 。L 1 ( 2 u ) = 1 − 2 u L_1(2u)=1-2u L 1 ( 2 u ) = 1 − 2 u と例題 5.4.1 の原始関数から λ 1 ( a ) = 1 − e − a ( 1 + 2 a ) \lambda_1(a)=1-\mathrm{e}^{-a}(1+2a) λ 1 ( a ) = 1 − e − a ( 1 + 2 a ) 。a = 1 / 2 a=1/2 a = 1/2 では λ 0 = 1 − e − 1 / 2 = 0.3935 \lambda_0=1-\mathrm{e}^{-1/2}=0.3935 λ 0 = 1 − e − 1/2 = 0.3935 、λ 1 = 1 − 2 / e = − 0.2131 \lambda_1=1-2/\sqrt\mathrm{e}=-0.2131 λ 1 = 1 − 2/ e = − 0.2131 。
(3) λ 1 ( a ) < 0 \lambda_1(a)<0 λ 1 ( a ) < 0 は e a < 1 + 2 a \mathrm{e}^{a}<1+2a e a < 1 + 2 a と同じで、e a = 1 + 2 a \mathrm{e}^a=1+2a e a = 1 + 2 a の正の解 a ∗ = 1.2564 a_*=1.2564 a ∗ = 1.2564 より小さい a a a で成り立つ。面積でいえば π ℏ a ∗ = 0.628 × 2 π ℏ \pi\hbar a_*=0.628\times2\pi\hbar π ℏ a ∗ = 0.628 × 2 π ℏ より小さい楕円では、χ ^ a \hat\chi_a χ ^ a は負の固有値を持つ。このとき χ ^ a \hat\chi_a χ ^ a は測定の作用素の条件 0 ≤ E ^ ≤ 1 0\le\hat E\le1 0 ≤ E ^ ≤ 1 を満たさず、「粒子が楕円 D a D_a D a の中にいる」という結果を持つ測定は存在しない。∫ D a W \int_{D_a}W ∫ D a W は χ ^ a \hat\chi_a χ ^ a という作用素の期待値ではあるが、何かの確率ではない。
固有値を a a a の関数として描くと図 5.4.2 になる。λ 2 \lambda_2 λ 2 は a = 1.71 a=1.71 a = 1.71 で − 0.239 -0.239 − 0.239 、λ 3 \lambda_3 λ 3 は a = 3.15 a=3.15 a = 3.15 で − 0.248 -0.248 − 0.248 まで下がる。どの n n n でも a → ∞ a\to\infty a → ∞ で λ n → 1 \lambda_n\to1 λ n → 1 になる(∫ 0 ∞ e − u L n ( 2 u ) d u = ( 1 − 2 ) n = ( − 1 ) n \int_0^\infty\mathrm{e}^{-u}L_n(2u)\,\mathrm{d} u=(1-2)^n=(-1)^n ∫ 0 ∞ e − u L n ( 2 u ) d u = ( 1 − 2 ) n = ( − 1 ) n から)。楕円が十分大きければ χ ^ a \hat\chi_a χ ^ a はほぼ「楕円の中にいる」という射影のようにふるまい、∫ D a W \int_{D_a}W ∫ D a W はほぼ確率として読める。面積が 2 π ℏ 2\pi\hbar 2 π ℏ の数倍以下の領域では、この読み方は成り立たない。
図 5.4.2 楕円 u < a u<a u < a の指示関数を Weyl 量子化した作用素の固有値 λ n ( a ) \lambda_n(a) λ n ( a ) (n = 0 , 1 , 2 n=0,1,2 n = 0 , 1 , 2 、厳密な値なので点線)。n = 1 n=1 n = 1 は a < 1.2564 a<1.2564 a < 1.2564 、n = 2 n=2 n = 2 は 0.74 < a < 2.86 0.74<a<2.86 0.74 < a < 2.86 で負になる。どの n n n でも a → ∞ a\to\infty a → ∞ で1に近づく。
ならせば負は消える
W W W の負の部分が面積 2 π ℏ 2\pi\hbar 2 π ℏ より細かい所にしか無いなら、その程度の幅でならせば消えるはずである。これは正確な定理になる。
命題 5.4.2(別の状態の W W W でならすと非負 ★★★)
ρ ^ \hat\rho ρ ^ と σ ^ \hat\sigma σ ^ を状態とし、σ ^ \hat\sigma σ ^ を ( x 0 , p 0 ) (x_0,p_0) ( x 0 , p 0 ) だけずらした状態を σ ^ x 0 , p 0 = D ^ σ ^ D ^ † \hat\sigma_{x_0,p_0}=\hat D\hat\sigma\hat D^\dagger σ ^ x 0 , p 0 = D ^ σ ^ D ^ † と書く。
W ˉ ( x 0 , p 0 ) = ∫ W ρ ^ ( x , p ) W σ ^ ( x − x 0 , p − p 0 ) d x d p = 1 2 π ℏ Tr ( ρ ^ σ ^ x 0 , p 0 ) ≥ 0 (5.4.7) \bar W(x_0,p_0)=\int W_{\hat\rho}(x,p)\,W_{\hat\sigma}(x-x_0,p-p_0)\,\mathrm{d} x\,\mathrm{d} p=\frac{1}{2\pi\hbar}\operatorname{Tr}\bigl(\hat\rho\,\hat\sigma_{x_0,p_0}\bigr)\ge0\tag{5.4.7} W ˉ ( x 0 , p 0 ) = ∫ W ρ ^ ( x , p ) W σ ^ ( x − x 0 , p − p 0 ) d x d p = 2 π ℏ 1 Tr ( ρ ^ σ ^ x 0 , p 0 ) ≥ 0 ( 5.4.7 )
で、∫ W ˉ d x 0 d p 0 = 1 \int\bar W\,\mathrm{d} x_0\,\mathrm{d} p_0=1 ∫ W ˉ d x 0 d p 0 = 1 である。
証明 (W4)(第5部第1章 )から W σ ^ ( x − x 0 , p − p 0 ) W_{\hat\sigma}(x-x_0,p-p_0) W σ ^ ( x − x 0 , p − p 0 ) は σ ^ x 0 , p 0 \hat\sigma_{x_0,p_0} σ ^ x 0 , p 0 の Wigner 関数である。(5.4.3) を当てると等号が出て、ρ ^ ≥ 0 \hat\rho\ge0 ρ ^ ≥ 0 、σ ^ x 0 , p 0 ≥ 0 \hat\sigma_{x_0,p_0}\ge0 σ ^ x 0 , p 0 ≥ 0 から2つの非負の作用素の積のトレースは 0 0 0 以上になる。規格化は、∫ W σ ^ ( x − x 0 , p − p 0 ) d x 0 d p 0 = 1 \int W_{\hat\sigma}(x-x_0,p-p_0)\,\mathrm{d} x_0\,\mathrm{d} p_0=1 ∫ W σ ^ ( x − x 0 , p − p 0 ) d x 0 d p 0 = 1 を先に積分すれば ∫ W ρ ^ = 1 \int W_{\hat\rho}=1 ∫ W ρ ^ = 1 に帰着する。
σ ^ \hat\sigma σ ^ を基底状態 ∣ 0 ⟩ | 0 \rangle ∣0 ⟩ (W 0 = 1 π ℏ e − u W_0=\frac{1}{\pi\hbar}\mathrm{e}^{-u} W 0 = π ℏ 1 e − u 、幅 σ x σ p = ℏ / 2 \sigma_x\sigma_p=\hbar/2 σ x σ p = ℏ/2 の Gauss)に取ると、σ ^ x 0 , p 0 \hat\sigma_{x_0,p_0} σ ^ x 0 , p 0 はコヒーレント状態 ∣ α ⟩ ⟨ α ∣ | \alpha \rangle \langle \alpha | ∣ α ⟩ ⟨ α ∣ (α = ( m ω x 0 + i p 0 / m ω ) / 2 ℏ \alpha=(\sqrt{m\omega}\,x_0+\mathrm{i} p_0/\sqrt{m\omega})/\sqrt{2\hbar} α = ( mω x 0 + i p 0 / mω ) / 2ℏ )になり、W ˉ = ∣ ⟨ α ∣ ψ ⟩ ∣ 2 / ( 2 π ℏ ) \bar W=\lvert \langle \alpha | \psi \rangle \rvert ^2/(2\pi\hbar) W ˉ = ∣⟨ α ∣ ψ ⟩ ∣ 2 / ( 2 π ℏ ) となる。測度を d 2 α = d x 0 d p 0 / ( 2 ℏ ) \mathrm{d}^2\alpha=\mathrm{d} x_0\,\mathrm{d} p_0/(2\hbar) d 2 α = d x 0 d p 0 / ( 2ℏ ) に替えたものが Husimi の Q Q Q 関数 Q ( α ) = ⟨ α ∣ ρ ^ ∣ α ⟩ / π Q(\alpha)=\langle \alpha | \hat\rho | \alpha \rangle /\pi Q ( α ) = ⟨ α ∣ ρ ^ ∣ α ⟩ / π で、第5部第6章 で光の状態の地図に使う。図 5.4.1 の実線は ∣ 1 ⟩ | 1 \rangle ∣1 ⟩ の場合で、⟨ α ∣ 1 ⟩ = e − ∣ α ∣ 2 / 2 α ∗ \langle \alpha | 1 \rangle =\mathrm{e}^{-\lvert \alpha\rvert ^2/2}\alpha^* ⟨ α ∣1 ⟩ = e − ∣ α ∣ 2 /2 α ∗ から W ˉ 1 = e − ∣ α ∣ 2 ∣ α ∣ 2 / ( 2 π ℏ ) \bar W_1=\mathrm{e}^{-\lvert \alpha\rvert ^2}\lvert \alpha\rvert ^2/(2\pi\hbar) W ˉ 1 = e − ∣ α ∣ 2 ∣ α ∣ 2 / ( 2 π ℏ ) である。幅の積が ℏ / 2 \hbar/2 ℏ/2 以上の中心0の Gauss 関数は、どれもある状態(一般にはスクイーズした熱平衡状態)の Wigner 関数になっている(第5部第1章 の総合問題で、熱平衡の W W W の幅の積が ℏ 2 coth β ℏ ω 2 \frac\hbar2\coth\frac{\beta\hbar\omega}{2} 2 ℏ coth 2 β ℏ ω だった)。したがって、面積 π ℏ \pi\hbar π ℏ 程度以上の Gauss の窓でならすと、どんな W W W も非負の確率密度になる。
この「ならした分布」は、実際の測定で得られる確率分布でもある。E. Arthurs と J. L. Kelly は、位置と運動量を別々の測定器に同時に結合し、両方の読みを1回の測定で得る模型を作った(E. Arthurs and J. L. Kelly (1965)、Bell Syst. Tech. J. 44 、725–729)。その読みの同時分布は、測定器の初期状態から決まる Gauss の窓で W W W をならした (5.4.7) の形になり、位置と運動量の両方の分散に雑音が上乗せされる。位置と運動量を同時に読む測定は存在するが、その結果の分布は W W W でなく W ˉ \bar W W ˉ で、W ˉ \bar W W ˉ は必ず非負である。
例題 5.4.3 ならした分布の原点の値 ★★
オリジナル 状態 ∣ 1 ⟩ | 1 \rangle ∣1 ⟩ を、基底状態の W 0 W_0 W 0 でならした分布 W ˉ 1 \bar W_1 W ˉ 1 と、幅の積が ℏ / 2 \hbar/2 ℏ/2 の別の Gauss の窓(σ x 2 = ℏ / ( m ω ) \sigma_x^2=\hbar/(m\omega) σ x 2 = ℏ/ ( mω ) 、σ p 2 = m ℏ ω / 4 \sigma_p^2=m\hbar\omega/4 σ p 2 = m ℏ ω /4 のスクイーズした窓)でならした分布について、原点 ( 0 , 0 ) (0,0) ( 0 , 0 ) での値を求めよ。
解答 (5.4.7) の右辺の形を使う。基底状態の窓では σ ^ 0 , 0 = ∣ 0 ⟩ ⟨ 0 ∣ \hat\sigma_{0,0}=| 0 \rangle \langle 0 | σ ^ 0 , 0 = ∣0 ⟩ ⟨ 0∣ で、∣ ⟨ 0 ∣ 1 ⟩ ∣ 2 = 0 \lvert \langle 0 | 1 \rangle \rvert ^2=0 ∣⟨ 0∣1 ⟩ ∣ 2 = 0 だから W ˉ 1 ( 0 , 0 ) = 0 \bar W_1(0,0)=0 W ˉ 1 ( 0 , 0 ) = 0 。スクイーズした窓は、∣ 0 ⟩ | 0 \rangle ∣0 ⟩ の位置の幅を 2 \sqrt2 2 倍、運動量の幅を 1 / 2 1/\sqrt2 1/ 2 倍にした純粋状態 ∣ ϕ ⟩ | \phi \rangle ∣ ϕ ⟩ (Gauss 波束で、実数の偶関数)である。∣ 1 ⟩ | 1 \rangle ∣1 ⟩ の波動関数は奇関数なので、偶関数との重なりは 0 0 0 で、やはり W ˉ ( 0 , 0 ) = 0 \bar W(0,0)=0 W ˉ ( 0 , 0 ) = 0 になる。どちらの窓でも、原点を中心とする窓は偶の状態で、奇の状態 ∣ 1 ⟩ | 1 \rangle ∣1 ⟩ と直交する。ならした分布は原点で 0 0 0 まで下がるが、負にはならない。W 1 ( 0 , 0 ) = − 1 / ( π ℏ ) W_1(0,0)=-1/(\pi\hbar) W 1 ( 0 , 0 ) = − 1/ ( π ℏ ) の負の値は、窓の中で周りの正の値と打ち消し合って消える。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com