密度行列(みつどぎょうれつ、英語: density matrix)は、量子力学における混合状態を表現するために使われる行列である。そこで本項ではまず混合状態とは何かについて解説し、その後に密度行列について解説する。

背景 - 混合状態編集

本節では、密度行列の概念の背後にある、混合状態の概念を説明する。次節ではこれを踏まえて密度行列の概念を説明する。

概要編集

量子力学では、系の状態は状態ベクトルもしくは純粋状態と呼ばれるベクトル  で書き表され、その振幅は波動関数

 

によって書き表されるが、こうした方法による系の記述方法は実験者がψ(x)の値を完全に知っているというのが暗黙の前提である。さもなければψ(x)を数式で書き表す事ができないので、系を数学的に解析ができなくなる。

しかしこの暗黙の前提は、実験者が系に関する情報を不完全にしか知らない場合には成り立たない。特に量子統計力学で想定されるような、数モル 個もの粒子を扱う状況下において、全ての粒子の情報を実験者が完全に知っていると仮定するのは現実的ではない。

そこでこうした、系に対する情報の不足石坂 et.al. 12:p104がある状況下における量子力学を記述するため、混合状態と呼ばれる複数の純粋状態に確率を付加した状態を考える必要がある。これは例えば「半分の確率で純粋状態 であり、残り半分の確率で純粋状態 になる」といったものが混合状態である。

関連概念との違い編集

混合状態の正確な理解のために、関連概念との違いを述べる。

観測との違い編集

まず混合状態でいうところの「確率」は古典的な確率論に従う確率(ベイズ確率)のことであり、系を観測した場合に波束がどこに収束するかを示す量子力学的な確率ではない。したがって混合状態は系の観測とは無関係に定義できる概念である。

これを踏まえた上で混合状態の定義をより正確に述べると、系が混合状態であっても系の状態は実は通常の純粋状態のいずれかである。しかし実験者の知識不足により、系がどの純粋状態にあるのかを(ベイズ確率的な意味で)推定する必要があり、「系が実はこの純粋状態にあると推定される確率」の分布が混合状態である。

重ね合わせとの違い編集

また混合状態とは量子力学的な状態の重ね合わせではない。これを偏光の例で説明する。光子には右円偏光と左円偏光がある。以下、右偏光と左偏光をそれぞれ純粋状態  で表すことにする。

量子力学では状態の重ね合わせが可能なので、ある光子の  の状態とこれと同じ光子の の状態を1/2ずつ重ね合わせると、光子は

 

という状態になる。これを規格化すれば

 

である。この状態にある光子は垂直方向に偏光であり、垂直方向の偏光板は通過できる。

これに対し、無偏光の状態にある光は、上記のような重ね合わせでは表現できず、混合状態によって記述する必要がある。無偏光の光とは例えば、光に含まれる複数の光子のうち50%の光子が の状態にあり、これらとは別の光子である残り50%の光子が の状態にある場合である。このような状態にある光の中に重ね合わせ状態 の光子は存在しない。また上述の状態にある光は横向きの偏光板を完全に通過するが、縦向きの偏光板にはある程度吸収されるなど、物理的性質も とは異なる。

このような「50%が 、残り50%が 」という統計的な状態、すなわち個々の粒子が「確率1/2で 、確率1/2で 」となっている状態を記述するのが混合状態である。

単純な記述方法の欠点編集

混合状態を記述する単純な記述方法は、状態ベクトルとその生起確率を並べて書く、というものである。例えば「確率1/2で 、確率1/2で 」という混合状態であれば、

 

と状態を記述する。しかしこの記述方法は、実質的に同一(=観測によって区別できない)の状態が複数の異なる表記を持ってしまうという欠点を持つ石坂 et.al. 12:p104-105

このため密度行列という表記方法を採用する必要があるのだが、これについては次章で述べることとし、本節ではまず上述した欠点を具体例で示す。

具体例編集

  

と定義し、

 

という状態記述を考えると、s'は前述のsは見かけ上全く異なるにも関わらず、観測によって両者は区別できない石坂 et.al. 12:p104-105。すなわち、どのような物理量Aを持ってきても、sAを観測したときの観測値の確率分布とs'Aを観測したときの観測値の確率分布は同一となる石坂 et.al. 12:p104-105

証明編集

実際、任意の観測値aに対し、

 の観測値が の観測値が の観測値が 

であり、

 の観測値が の観測値が の観測値が 

が成立するので両者は等しい。ここでPaaの固有空間への射影である。

密度行列編集

密度行列は、混合状態を数学的に記述する為の道具立てであり、しかも上述した単純な記述方法のような欠点を持たない事である。

本章では、密度行列の定義とその性質を述べ、次章において単純な記述方法の欠点が解消されている事を見る。

数学的考察編集

密度行列の概念を導入する前準備として、簡単な数学的考察を行う。

系が純粋状態 にあるとき、物理量 を観測すると、観測値の期待値は

 

となる。これを変形すれば以下のようになる:

 

ここで  、…は完全正規直交系であり、  は行列 トレース (trace) である。

したがってより一般に「p1の確率で p2の確率で 、…」という混合状態を観測すればその期待値は

 

である。

そこでこの混合状態の密度演算子を対角行列

 

によって定義し、そのその行列表示密度行列と呼ぶことにすると、混合状態にある際の の観測値の期待値は

 

という簡単な形で書き表す事ができる。以上のことから、密度行列は混合状態にある系の観測値の期待値を計算するのに便利である。

定義編集

以上を踏まえた上で、密度行列とその関連概念を以下のように定義する。

密度行列の厳密な定義編集

状態空間上の完全正規直交系  、…に対し、状態空間における 方向の射影作用素Pkとするとき、

    ....(M1)

という形で表記できる演算子  密度演算子 (density operator )もしくは密度行列 (density matrix ) という。(上式右辺の収束はトレースノルムに関するものである新井08:p81)。

なお、射影作用素Pkブラ-ケット記法では

 

と書けるので、上述の定義は前節で述べた定義と実質的に一致する。しかしブラ-ケット記法は文脈により数学的な定式化方法が異なるので、本節では定義を厳密に記述する為、射影作用素Pkを用いて密度行列を定義した。

また上の定義では、  、…が正規直交系をなしている事を仮定したが、必ずしもこれは必須ではない。しかし正規直交ではない  、…に対して同様に密度行列を定義したとしても、必ず完全正規直交基底の表現に書き換えられる事が知られている(次節の別定義との同値性から従う)。

別定義編集

 が上述したように書ける必要十分条件は、以下の3つを満たす事が知られている新井08:p81

  •  有界自己共役作用素
  •  は非負の作用素である。すなわち が状態空間上の任意の状態ベクトル に対して成立する。
  •  

よってこの3条件を満たす事を密度行列の定義としても良い。

なお、 は状態空間上の完全正規直交系  、…を用いて

 

により定義されるH13:p421。この値は完全正規直交系の取り方に依存しない為、well-definedであるH13:p421

本節の方法で定義した密度行列を前節の(M1)式の形で表す事を、密度行列のシャッテン分解という新井08:p81

純粋状態の定義編集

状態ベクトル に対し、状態空間における 方向の射影作用素Pψとする。

密度行列が何らかの純粋状態の状態ベクトル を用いて

 

と書ける時、 純粋状態にあるというH13:p426。前節で述べたように、ブラ-ケット記法では上式は

 

を意味する。

観測編集

観測の期待値編集

密度演算子 と有界な自己共役作用素  に対し、   が定義可能で

 

が成立することが知られているH13:p423

既に述べたように、密度行列 で表現される混合状態において を観測した際の観測値の期待値はこの値になる。

波束の収縮編集

密度行列 で記述される混合状態に対して物理量 を観測した結果、 の固有値λを得たとすると、波束の収縮が起こり、密度行列は

 

になるH13:p428

密度行列集合の凸性編集

密度行列の全体の集合は凸集合である事が知られている。すなわち、  を密度行列とし、u0≦u≦1満たす実数とする時、  重ね合わせ

 

も密度行列であるH13:p426

また、この凸集合の「端っこ」にあるのは純粋状態である。すなわち が純粋状態である必要十分条件は、

 

を満たす密度行列  、および実数0<u<1が存在しない事であるH13:p426

注意編集

上で定義した重ね合わせの概念は、状態ベクトルの重ね合わせとは異なる概念である。実際、一般には

 

であるH13:p426

両者を区別するため、状態ベクトルの重ね合わせをコヒーレントな重ね合わせ、密度行列の重ね合わせをインコヒーレントな重ね合わせというH13:p427

密度行列の公理的特徴づけ編集

本章では密度行列を全く別の角度から公理的に特徴づける。そしてこの特徴づけができる事の結果として、前述した単純な表記方法の持つ欠点が密度行列では解消されている事を見る。

公理的特徴付け編集

Aを物理量、すなわち状態空間上の自己共役作用素であるとする。

今何らかの量子力学的な系が与えられていたとし、この系でAを観測した観測値の期待値を

 

と書くことにする。なお系の具体的な状態は問わない。したがって系が純粋状態であっても混合状態であってもよい。


 は自己共役作用素Aに実数を対応させる関数

 

とみなす事ができるが、物理的に考えると、この関数は次の2性質を満たさねばならないはずである。なお以下でIは単位行列である。さらにA非負であるとは、任意の状態ベクトル に対し が成立する事を言う:

(1)  
(2) Aが非負なら、 


なぜこれらの条件が要請されるかというと、単位行列Iの固有値は全て1なので、Aを観測した結果は常に1でなければならない(=条件(1))。またAが非負になるにはその固有値(や連続スペクトル)が全て非負になる場合だけなので、 が成立しなければならない(=条件(2))。


さらに関数 が以下の連続性を満たしている事を要請する:

(3)  のとき となる任意の自己共役作用素の列 に対し、 

ここで は実数としての収束であり、 はL2ノルムに関するweak-*収束である。


このとき次が成立する事が知られている:

定理 ― 有界な自己共役作用素Aに実数を対応させる線形汎関数 が(1)、(2)、(3)をすべて満たす必要十分条件は、

 

を満たす密度行列ρが存在する事であるH13:p423-424。しかもそのような密度行列は一意に定まるH13:p423-424

欠点が解消されている事編集

前の章で述べたように、混合状態を単純な方法で記述した場合、見かけ上の記述が異なるにも関わらず、実質的に同一の量子状態を表している(=観測結果では両者を区別できない)、という事が起こりうる。

しかし密度行列を用いて混合状態を記述した場合にはこのような問題は生じない。

実際、2つの量子状態が実質的に同一である(=観測結果では両者を区別できない)という事は、この2つの量子状態に対する関数 が同一であるということを意味し、 が同一であるという事は対応する密度行列ρが同一だという事を意味するからである。したがって実質的に同一の量子状態が相異なる2つの密度行列で表示できる事はありえない。

純粋状態を記述する上での数学的利点編集

本項では密度行列を、混合状態を記述する上での便利な道具立てとして導入した。しかし純粋状態を記述する際にも密度行列は有効に働く[1]

これは状態ベクトル表記もやはり、(位相以外にも差がある)全く別の状態ベクトルが、同一の純粋状態を表す場合があるからである。前述のように密度行列であればこうした問題は生じない。

状態ベクトルに対してこの問題が生じるのは、合成系の場合である。2つの系を合成した場合、合成系の状態ベクトルは、各々の系の状態ベクトル  のテンソル積である:

 

位相にしか差がない2つの状態ベクトルは同一の物理状態を表すので、

 

 

は同一の物理状態を表す。しかしθが0でない限り、

 

を満たすαは存在しない。

すなわち  は(位相以外にも差がある)全く別の状態ベクトルであるにもかかわらず、同一の量子状態を表す。

密度行列による問題解決編集

これに対し、密度行列を利用した場合は、上述の問題は生じない。そもそも、上述の問題が生じたのは、状態ベクトルに位相分の自由度

 

が存在したからである。しかし密度行列で記述した場合、純粋状態は

 

という形式なので、位相分の自由度は消え去る:

 

よって前述の問題はそもそも生じない。

フォン・ノイマンエントロピー編集

密度行列は何らかの混合状態を表し、混合状態とは純粋状態の集合に何らかの確率分布を付与したものである。よってこの確率分布に対して情報理論におけるシャノンエントロピー(情報量)を定義することができ、これにボルツマン定数をかけたものを密度行列のフォン・ノイマンエントロピーという。本項ではまず、シャノンエントロピーの概念を復習し、これをベースにフォン・ノイマンエントロピーの概念を定義する。

シャノンエントロピーとは編集

情報理論では、確率1/2で表がでるコインを単位として、事象の確率がコイン何枚分に相当するかを考える。例えば確率1/8=(1/2)3で起こる事象があったとき、この確率はコイン3枚全てが表になる確率に相当するので、この事象の「自己情報量」は

 

であると定義する。より一般に、確率pで起こる事象があった場合、この事象のaに対する自己情報量

 

により定義する。コインを単位にする場合は、底のaは2である。

また値1、2、3、…を取る確率変数Xがあった時、X=jであるという事象の自己情報量は

 

であるので、Lj期待値

 

を定義でき、この値をXaに対する情報量、もしくはaに対するシャノンエントロピーという。

ただし である項に関しては、

 

であるので、

 

とみなす。

本項で重要なのは底a自然対数eの場合なので、底eに対するシャノンエントロピーを単にシャノンエントロピーと呼び、

 

と略記する。

フォン・ノイマンエントロピーの定義編集

定義編集

  、…を完全正規直交系とする。密度行列

 

に対し、 フォン・ノイマンエントロピー 

 

と呼ぶ。ここでkBボルツマン定数である。

なおシャノンエントロピーの場合と同様、上述の定義で とみなす。

別定義編集

密度行列 に対し、作用素解析の手法により

 

を定義する事ができ、

 

によりフォン・ノイマンエントロピーを定義する事ができる新井08:p190-191。この定義は前述した定義と一致する新井08:p190-191

性質編集

任意の密度行列 に対し、フォン・ノイマンエントロピーは

 

を満たす。

また、S(ρ)=0となる必要十分条件は は純粋状態にある事であるH13:p426。したがってフォン・ノイマンエントロピーは純粋状態からの「ズレ」を表す量だと解釈できる。

フォン・ノイマンエントロピーは通常の観測(射影観測)を行った場合には、増加するかも知れないが、減少する事はない。しかしより一般的な観測をした場合には減少する場合がある[2][3]

量子相互作用を混合系の中で消去することにより、観測は「情報を減少させる」。—量子もつれ, einselection, や量子デコヒーレンスを参照。すなわち孤立していない系のフォン・ノイマンエントロピーを減少させる事はできるが、これは系の外部のフォン・ノイマンエントロピーを上昇させている場合のみであり、系の内外のフォン・ノイマンエントロピーは減少しない。。熱力学の第二法則熱力学と情報理論のエントロピー英語版(Entropy in thermodynamics and information theory)を参照。

時間発展編集

密度演算子の時間発展は、次のフォン・ノイマン方程式 (von Neumann equation ) で記述される。フォン・ノイマン方程式は古典論におけるリウヴィル方程式 (Liouville equation ) に対応するので、リウヴィル=フォン・ノイマン方程式 (Liuville–von Neumann equation)、あるいは単に(量子)リウヴィル方程式とも呼ばれる。

 

ここで ħ = h/2π換算プランク定数hプランク定数)、 ハミルトニアン、括弧 [·,·]交換子である。

フォン・ノイマンの式は、純粋状態(状態ベクトル)の時間発展を記述するシュレーディンガー方程式

 

密度演算子の定義式だけを用いて導出できる。ここでブラ・ベクトル 〈Ψ| はケット・ベクトル |Ψ〉双対であること 〈Ψ| = |Ψ〉 に注意。

統計力学への応用編集

統計力学においては、状態のアンサンブルを混合状態と考えることができる。量子統計力学では、あるハミルトニアンの各エネルギー固有状態が混合していると考えて密度行列を表現することがよくある。

密度行列 ρ は、たとえば混合の比率がカノニカル分布で表せるとすると、

 

グランドカノニカル分布では、

 

で表される。ここで β = 1/kBT逆温度kBボルツマン定数ΩグランドポテンシャルHG はグランドカノニカル分布でのハミルトニアンである。

このときオブザーバブルの期待値 A は、

 

と書くことができる。特に A恒等演算子 A = Id の場合、

 

を満たす。また、A がハミルトニアン A = H の場合、ハミルトニアンの固有値を {Ei} とすれば、

 

と書き換えられる。

量子リウヴィル方程式、モーヤル方程式編集

密度行列演算子は相空間の中でも実現される。ウィグナー函数の下では、等価なウィグナー函数への密度行列変換は、

 

である。このウィグナー函数の時間発展の方程式は、上記のフォン・ノイマン函数のウィグナー変換である。

  .

ここに H(q,p) はハミルトニンであり、{ { •,• } } はモーヤルの括弧英語版(Moyal bracket)、量子交換子の変換関係である。

ウィグナー函数の発展方程式は、古典極限の発展方程式、古典物理学リウヴィル方程式の類似である。プランク定数 ħ が 0 となる極限では、W(q,p,t)相空間の古典リウヴィル確率分布函数へと還元される。

古典リウヴィル方程式は、偏微分方程式の特性曲線法を使い解くことができ、特性曲線はハミルトン方程式である。同じように、量子力学でのモーヤル方程式は、量子特性曲線法英語版(quantum characteristic)を用いて解、すなわち、相空間のモーヤル積英語版(Moyal−product)を求めることができる。実践的には、解を求める方法は異る方法を用いる。

脚注編集

注釈編集

出典編集

  1. ^ 本節はH13の19.1節を参考にした。
  2. ^ Nielsen, Michael; Chuang, Isaac (2000), Quantum Computation and Quantum Information, Cambridge University Press, ISBN 978-0-521-63503-5 . Chapter 11: Entropy and information, Theorem 11.9, "Projective measurements cannot decrease entropy"
  3. ^ Everett, Hugh (1973), “The Theory of the Universal Wavefunction (1956) Appendix I. "Monotone decrease of information for stochastic processes"”, The Many-Worlds Interpretation of Quantum Mechanics, Princeton Series in Physics, Princeton University Press, pp. 128–129, ISBN 978-0-691-08131-1 

参考文献編集

  • [H13] Brian C.Hall (2013/7/1). Quantum Theory for Mathematicians. Graduate Texts in Mathematics 267. Springer 
  • [新井08] 新井朝雄 (2008/7/10). 量子統計力学の数理. 共立出版. ISBN 978-4320018655 
  • [石坂 et.al. 12] 石坂智 、小川朋宏、河内亮周、木村元、林正人 (2012/6/8). 量子情報科学入門. 共立出版. ISBN 978-4320122994 

関連記事編集