量子力学の様々な構成方法 第5部 第4章 §5.4.1 ★★★ 負の値の場所と量。直交性から負の値の存在が従うことと、負の量を測る負の体積
第5部の構成は、密度行列の代わりに位相空間の関数 W ( x , p ) W(x,p) W ( x , p ) を状態とみなす。第5部第1章 では、この関数に古典の分布と同じ性質を求めると、実数値・規格化・線形性・周辺分布・平行移動との整合の5つは両立するが、「至る所で非負」だけは両立しないことを示した。第5部第2章 では、W W W と密度行列 ρ ^ \hat\rho ρ ^ が1対1に対応し、期待値がすべて ∫ W A W d x d p \int W A_{\mathrm W}\,\mathrm{d} x\,\mathrm{d} p ∫ W A W d x d p の形で書けることを証明した。すると、W W W が負になる点を含む状態でも、計算した期待値や確率はすべて基準(第1部第2章 )と一致している。負の値を持つ関数が確率のように使えて、しかも矛盾を起こさないのはなぜか。この章はこの問いに答える。
部の型(出発点の要請 → 基本の対象 → 基準との同値 → 共通例題 → 得意な問題 → 限界)のうち、この章は「基準との同値」の後に、この構成だけが持つ特徴(負の値)を調べる位置にある。結論は次の4点にまとめられる。第一に、負の値は干渉と同じ起源を持ち、相空間の面積 2 π ℏ 2\pi\hbar 2 π ℏ より細かい所にしか現れない。どんな状態の W W W でも、別の状態の W W W でならすと非負になる。第二に、W W W の負の値はどの測定の結果の確率でもない。「点 ( x , p ) (x,p) ( x , p ) のまわりの小さな領域に粒子がいる」という事象には、それを測る量子力学の測定が存在しない。第三に、純粋状態で W W W が至る所で非負になるのは Gauss 波束だけである(Hudson の定理)。この章で証明する。第四に、確率論の言葉では、「すべての1次結合 a x + b p ax+bp a x + b p の分布を同時に再現する確率分布」は存在すれば W W W そのものに限られる(Cramér–Wold の定理)。W W W が負なら、そのような同時分布は存在しない。
この章では1次元・1粒子に話を限る。調和振動子の計算では無次元の量 u = 2 H / ( ℏ ω ) u=2H/(\hbar\omega) u = 2 H / ( ℏ ω ) (H = p 2 / ( 2 m ) + m ω 2 x 2 / 2 H=p^2/(2m)+m\omega^2x^2/2 H = p 2 / ( 2 m ) + m ω 2 x 2 /2 )を使う。u u u だけの関数を積分するときは、相空間の面積要素が d x d p = π ℏ d u \mathrm{d} x\,\mathrm{d} p=\pi\hbar\,\mathrm{d} u d x d p = π ℏ d u になる。
1点の値はずらしたパリティの期待値
第5部第1章 で示したとおり、W W W の1点 ( x 0 , p 0 ) (x_0,p_0) ( x 0 , p 0 ) での値は、状態をずらしてその点を原点に持ってきたときのパリティの期待値に比例する。位置を x 0 x_0 x 0 、運動量を p 0 p_0 p 0 だけずらすユニタリを D ^ ( x 0 , p 0 ) = exp ( i ( p 0 x ^ − x 0 p ^ ) / ℏ ) \hat D(x_0,p_0)=\exp\bigl(\mathrm{i}(p_0\hat x-x_0\hat p)/\hbar\bigr) D ^ ( x 0 , p 0 ) = exp ( i ( p 0 x ^ − x 0 p ^ ) /ℏ ) 、パリティを ( Π ^ ψ ) ( x ) = ψ ( − x ) (\hat\Pi\psi)(x)=\psi(-x) ( Π ^ ψ ) ( x ) = ψ ( − x ) として、ずらしたパリティ を
Π ^ x 0 , p 0 = D ^ ( x 0 , p 0 ) Π ^ D ^ ( x 0 , p 0 ) † (5.4.1) \hat\Pi_{x_0,p_0}=\hat D(x_0,p_0)\,\hat\Pi\,\hat D(x_0,p_0)^\dagger\tag{5.4.1} Π ^ x 0 , p 0 = D ^ ( x 0 , p 0 ) Π ^ D ^ ( x 0 , p 0 ) † ( 5.4.1 )
と置くと、混合状態まで含めて
W ( x 0 , p 0 ) = 1 π ℏ Tr ( ρ ^ Π ^ x 0 , p 0 ) (5.4.2) W(x_0,p_0)=\frac{1}{\pi\hbar}\operatorname{Tr}\bigl(\hat\rho\,\hat\Pi_{x_0,p_0}\bigr)\tag{5.4.2} W ( x 0 , p 0 ) = π ℏ 1 Tr ( ρ ^ Π ^ x 0 , p 0 ) ( 5.4.2 )
が成り立つ(純粋状態の式を (W2) の線形性で混合状態へ広げた)。Π ^ x 0 , p 0 \hat\Pi_{x_0,p_0} Π ^ x 0 , p 0 は点 ( x 0 , p 0 ) (x_0,p_0) ( x 0 , p 0 ) を中心とする鏡映で、固有値は ± 1 \pm1 ± 1 だけである。したがって W ( x 0 , p 0 ) < 0 W(x_0,p_0)<0 W ( x 0 , p 0 ) < 0 は、「点 ( x 0 , p 0 ) (x_0,p_0) ( x 0 , p 0 ) のまわりで鏡映したとき、奇の成分の重みが偶の成分の重みより大きい」ことを意味する。負の値は「確率が負」なのでなく、固有値 ± 1 \pm1 ± 1 の物理量の期待値が負であることを表している。この読み方は 2節と6節でもう一度使う。
直交性から負の値の存在が従う
第5部第2章 のトレースの公式から、2つの状態の重なりは
Tr ( ρ ^ 1 ρ ^ 2 ) = 2 π ℏ ∫ W 1 W 2 d x d p (5.4.3) \operatorname{Tr}(\hat\rho_1\hat\rho_2)=2\pi\hbar\int W_1W_2\,\mathrm{d} x\,\mathrm{d} p\tag{5.4.3} Tr ( ρ ^ 1 ρ ^ 2 ) = 2 π ℏ ∫ W 1 W 2 d x d p ( 5.4.3 )
で与えられる。直交する2つの純粋状態では左辺が 0 0 0 になる。もし W 1 W_1 W 1 と W 2 W_2 W 2 がどちらも至る所で正なら、右辺の被積分関数は正で、積分も正になる。したがって次の命題が成り立つ。
命題 5.4.1(至る所で正の W W W を持つ状態に直交する状態 ★★★)
状態 ρ ^ 1 \hat\rho_1 ρ ^ 1 の Wigner 関数 W 1 W_1 W 1 が至る所で正(W 1 > 0 W_1>0 W 1 > 0 )で、ρ ^ 2 \hat\rho_2 ρ ^ 2 が Tr ( ρ ^ 1 ρ ^ 2 ) = 0 \operatorname{Tr}(\hat\rho_1\hat\rho_2)=0 Tr ( ρ ^ 1 ρ ^ 2 ) = 0 を満たすなら、W 2 W_2 W 2 はどこかで負になる。
証明 W 2 ≥ 0 W_2\ge0 W 2 ≥ 0 と仮定する。W 2 W_2 W 2 は連続で ∫ W 2 = 1 \int W_2=1 ∫ W 2 = 1 なので、ある開集合の上で正である。そこで W 1 W 2 > 0 W_1W_2>0 W 1 W 2 > 0 だから2 π ℏ ∫ W 1 W 2 > 0 2\pi\hbar\int W_1W_2>0 2 π ℏ ∫ W 1 W 2 > 0 となり、(5.4.3) の左辺が 0 0 0 であることに矛盾する。
調和振動子の基底状態の W 0 = 1 π ℏ e − u W_0=\frac{1}{\pi\hbar}\mathrm{e}^{-u} W 0 = π ℏ 1 e − u は至る所で正なので、励起状態 ∣ n ⟩ | n \rangle ∣ n ⟩ (n ≥ 1 n\ge1 n ≥ 1 )の W n W_n W n はすべてどこかで負になる。共通例題①の形(第5部第5章 )
W n = ( − 1 ) n π ℏ e − u L n ( 2 u ) , u = 2 H ℏ ω (5.4.4) W_n=\frac{(-1)^n}{\pi\hbar}\mathrm{e}^{-u}L_n(2u),\qquad u=\frac{2H}{\hbar\omega}\tag{5.4.4} W n = π ℏ ( − 1 ) n e − u L n ( 2 u ) , u = ℏ ω 2 H ( 5.4.4 )
でこれを確かめると、∫ W 0 W n d x d p = ( − 1 ) n π ℏ ∫ 0 ∞ e − 2 u L n ( 2 u ) d u = 0 \int W_0W_n\,\mathrm{d} x\,\mathrm{d} p=\frac{(-1)^n}{\pi\hbar}\int_0^\infty\mathrm{e}^{-2u}L_n(2u)\,\mathrm{d} u=0 ∫ W 0 W n d x d p = π ℏ ( − 1 ) n ∫ 0 ∞ e − 2 u L n ( 2 u ) d u = 0 (n ≥ 1 n\ge1 n ≥ 1 )となり、被積分関数が符号を変えることがわかる。命題 5.4.1 の議論は Gauss 波束1つと直交するどんな状態にも当てはまる。奇の状態はすべて原点の Gauss 波束に直交するので、どれも負の値を持つ。
負の体積
負の値の大きさを1つの数で測るには、∣ W ∣ \lvert W\rvert ∣ W ∣ の積分が1をどれだけ超えるかを見ればよい。
δ = ∫ ∣ W ∣ d x d p − 1 = 2 ∫ W < 0 ∣ W ∣ d x d p (5.4.5) \delta=\int\lvert W\rvert \,\mathrm{d} x\,\mathrm{d} p-1=2\int_{W<0}\lvert W\rvert \,\mathrm{d} x\,\mathrm{d} p\tag{5.4.5} δ = ∫ ∣ W ∣ d x d p − 1 = 2 ∫ W < 0 ∣ W ∣ d x d p ( 5.4.5 )
を負の体積 と呼ぶ。2つ目の等号は ∫ W = 1 \int W=1 ∫ W = 1 から出る。W ≥ 0 W\ge0 W ≥ 0 なら δ = 0 \delta=0 δ = 0 で、負の部分が大きいほど δ \delta δ は大きい。この量は A. Kenfack と K. Życzkowski が非古典性の指標として導入した(A. Kenfack and K. Życzkowski (2004)、J. Opt. B 6 、396)。
例題 5.4.1 第1励起状態の負の円板 ★★★
オリジナル 調和振動子の第1励起状態 ∣ 1 ⟩ | 1 \rangle ∣1 ⟩ について次を求めよ。
(1) W 1 W_1 W 1 が負になる領域と、その領域の面積。
(2) 負の領域での W 1 W_1 W 1 の積分と、負の体積 δ 1 \delta_1 δ 1 。
(3) (1) の面積を、純粋状態1つが占める面積 2 π ℏ 2\pi\hbar 2 π ℏ (第5部第1章 の命題)と比べよ。
解答 (1) L 1 ( y ) = 1 − y L_1(y)=1-y L 1 ( y ) = 1 − y なので (5.4.4) は W 1 = − 1 π ℏ e − u ( 1 − 2 u ) W_1=-\frac{1}{\pi\hbar}\mathrm{e}^{-u}(1-2u) W 1 = − π ℏ 1 e − u ( 1 − 2 u ) 。負になるのは u < 1 / 2 u<1/2 u < 1/2 、つまり H < ℏ ω / 4 H<\hbar\omega/4 H < ℏ ω /4 の楕円の内側である。d x d p = π ℏ d u \mathrm{d} x\,\mathrm{d} p=\pi\hbar\,\mathrm{d} u d x d p = π ℏ d u なので、面積は π ℏ ⋅ 1 2 = π ℏ 2 \pi\hbar\cdot\frac12=\frac{\pi\hbar}{2} π ℏ ⋅ 2 1 = 2 π ℏ 。
(2) d d u ( e − u ( 1 + 2 u ) ) = e − u ( 1 − 2 u ) \frac{\mathrm{d}}{\mathrm{d} u}\bigl(\mathrm{e}^{-u}(1+2u)\bigr)=\mathrm{e}^{-u}(1-2u) d u d ( e − u ( 1 + 2 u ) ) = e − u ( 1 − 2 u ) を使うと、
∫ u < 1 / 2 W 1 d x d p = − ∫ 0 1 / 2 e − u ( 1 − 2 u ) d u = − [ e − u ( 1 + 2 u ) ] 0 1 / 2 = − ( 2 e − 1 ) = − 0.2131. \int_{u<1/2}W_1\,\mathrm{d} x\,\mathrm{d} p=-\int_0^{1/2}\mathrm{e}^{-u}(1-2u)\,\mathrm{d} u=-\bigl[\mathrm{e}^{-u}(1+2u)\bigr]_0^{1/2}=-\Bigl(\frac{2}{\sqrt{\mathrm{e}}}-1\Bigr)=-0.2131 . ∫ u < 1/2 W 1 d x d p = − ∫ 0 1/2 e − u ( 1 − 2 u ) d u = − [ e − u ( 1 + 2 u ) ] 0 1/2 = − ( e 2 − 1 ) = − 0.2131.
(5.4.5) から δ 1 = 2 × 0.2131 = 0.4261 \delta_1=2\times0.2131=0.4261 δ 1 = 2 × 0.2131 = 0.4261 。
(3) 負の領域の面積 π ℏ / 2 \pi\hbar/2 π ℏ/2 は 2 π ℏ 2\pi\hbar 2 π ℏ の 1 / 4 1/4 1/4 である。W 1 W_1 W 1 の負の部分は、状態1つの大きさより小さい楕円の中に収まっている。原点の値 − 1 / ( π ℏ ) -1/(\pi\hbar) − 1/ ( π ℏ ) は 第5部第1章 の上限 ∣ W ∣ ≤ 1 / ( π ℏ ) \lvert W\rvert \le1/(\pi\hbar) ∣ W ∣ ≤ 1/ ( π ℏ ) の下側にちょうど達している。
同じ計算を n = 2 , … , 5 n=2,\dots,5 n = 2 , … , 5 で数値的に行うと、δ n = 0.729 , 0.977 , 1.191 , 1.383 \delta_n=0.729,\ 0.977,\ 1.191,\ 1.383 δ n = 0.729 , 0.977 , 1.191 , 1.383 で、n n n とともに増える。W n W_n W n は L n ( 2 u ) L_n(2u) L n ( 2 u ) の n n n 個の零点で符号を変える同心の輪になり、輪の幅は相空間の面積で測っておよそ 2 π ℏ 2\pi\hbar 2 π ℏ 程度である。もう1つの典型は、離れた2つの Gauss 波束の重ね合わせ(猫状態)で、第5部第2章 の総合問題で見たとおり、2つの山の中間に cos ( 2 a p / ℏ ) \cos(2ap/\hbar) cos ( 2 a p /ℏ ) の縞(2 a 2a 2 a は山の間隔)が立ち、縞の谷で W W W が負になる。縞の運動量方向の周期は π ℏ / a \pi\hbar/a π ℏ/ a で、山の間隔と縞の周期の積はつねに 2 π ℏ 2\pi\hbar 2 π ℏ である。どちらの例でも、負の値は相空間の面積 2 π ℏ 2\pi\hbar 2 π ℏ より細かい構造の中にある。
図 5.4.1 第1励起状態の W 1 W_1 W 1 (点線)と、それを基底状態の W 0 W_0 W 0 でならした分布 W ˉ 1 \bar W_1 W ˉ 1 (実線。2節の命題 5.4.2 )の p = 0 p=0 p = 0 での断面。縦軸は ℏ \hbar ℏ を掛けて無次元にした。W 1 W_1 W 1 は ∣ x ∣ < ℏ / ( 2 m ω ) \lvert x\rvert <\sqrt{\hbar/(2m\omega)} ∣ x ∣ < ℏ/ ( 2 mω ) で負になり、原点で − 1 / π -1/\pi − 1/ π に達する。ならした分布は原点で 0 0 0 に下がるが負にはならない。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com