密度行列

1 密度行列の基礎

1.1 定義と基本性質

1.1.1 エルミート性

密度行列は量子系の統計的状態を表す行列（演算子）であり、計算上は一般に有限次元ヒルベルト空間上の正方行列として扱われる。密度行列 ρ は自己随伴（エルミート）であることが要求される。これは確率解釈の整合性を支える条件であり、観測量の期待値が実数になることを保証する役割を持つ。

1.1.2 正の半定義性

密度行列は正の半定義でなければならない。つまり ρ の固有値はすべて非負である。固有値が非負であることにより、後述する測定確率が負になる状況を排除できる。物理的には「いかなる測定結果に対しても確率が整合的に定義される」ことを保証する数学条件として理解できる。

1.1.3 トレース1規格化

密度行列はトレースが 1 に等しい。トレースはスペクトル分解における固有値の総和に一致するため、非負かつ総和が 1 という性質は確率分布としての総量を与える。これにより、測定の結果全体にわたる確率の和が 1 となることが確保される。

1.2 混合状態と純粋状態

1.2.1 ユニタリと純粋状態の表現

純粋状態は波動関数あるいは状態ベクトル	ψ⟩ により記述され、密度行列では ρ =	ψ⟩⟨ψ	と表される。このとき ρ は冪等性（ρ^2 = ρ）を満たし、ランクが 1 になる。ユニタリ変換 U によって状態が変わる場合、密度行列は ρ ↦ UρU† の形で更新され、純粋性が保たれる。

1.2.2 混合状態の古典的混合との関係

混合状態は、異なる純粋状態が確率 p_i で混ざっている状況として直感的に捉えられることが多い。密度行列の言い換えとして ρ = Σ_i p_i	ψ_i⟩⟨ψ_i	が用いられ、係数 p_i は確率として解釈される。重要な点は、同じ ρ を与える分解が一意ではないことであり、これは「古典的な確率混合だけで完全に情報が定まる」と限らないことを示す。量子の特徴は、取りうる状態の重ね合わせと相関の扱いに現れる。

1.3 期待値の計算規則

1.3.1 観測量の期待値

観測量をエルミート演算子 A とするとき、期待値は ⟨A⟩ = Tr(ρA) により一様に与えられる。純粋状態の場合は	ψ⟩⟨ψ	を代入すれば ⟨ψ	A	ψ⟩ に戻り、密度行列が状態の一般化記述になっていることが分かる。この枠組みでは、測定以前の「統計的平均」が ρ によって統制される。

1.3.2 測定確率との対応

一般化測定として射影測定 {Π_k}（あるいはポジティブ作用素による測定）を考えると、測定結果 k の確率は p(k) = Tr(ρΠ_k) と表される。正の半定義性とトレース1規格化がこの確率の非負性と総和1を保証する。つまり密度行列は、確率分布の生成器としても機能し、観測の統計を直接与える。

2 密度行列の表現と操作

2.1 行列表現

2.1.1 計算基底での書き方

有限次元では基底 {	e_i⟩} を選ぶと、密度行列は成分 ρ_{ij} = ⟨e_i	ρ	e_j⟩ をもつ行列として表される。基底の取り方に応じて成分は変化するが、トレースや固有値のような不変量は同じままである。計算ではこの自由度を利用して、問題の対称性に合わせた基底で表現することが多い。

2.1.2 固有分解（スペクトル分解）

エルミート性により密度行列は実数の固有値を持ち、固有ベクトルの直交基底に関して ρ = Σ_k λ_k	k⟩⟨k	と書ける。固有値 λ_k は非負で総和が 1 になる。固有分解は、状態の混合度や計算可能な量（エントロピーなど）を評価する際の基本的な手がかりとなる。

2.2 ブラ・ケットによる記述

2.2.1 背景となる状態空間

密度行列の演算子としての意味は、状態が属するヒルベルト空間の構造に依存する。有限次元では直交基底により容易に扱えるが、連続変数を含む系では汎関数解析的な拡張が必要になることがある。いずれの場合も、物理量の期待値が Tr(ρA) で与えられるという構造が中心概念である。

2.2.2 作用素としての密度行列

密度行列は「状態ベクトル」そのものではなく、「作用素」としての役割を持つ。したがって線形代数の演算（積や随伴）は適切な条件のもとで意味を持ち、一般に ρ は 2点間の情報（複素位相を含む）を反映する。ユニタリ変換や量子チャネルを作用素として適用できる点が、密度行列が操作論的に扱いやすい理由である。

2.3 部分系と縮約密度行列

2.3.1 部分トレース

合成系の状態が ρ_AB として与えられているとき、部分系 A の状態は縮約密度行列 ρ_A = Tr_B(ρ_AB) により得られる。部分トレースは、他方の部分系の自由度を「観測しない」ことに対応し、統計的には周辺化の役割を担う。非負性とトレース1は保存されるため、結果として密度行列として正しい形になる。

2.3.2 出現する代表的状況（混合・トレーシング）

合成系が純粋でも、片方の部分系を見れば混合になることがある。これはもつれに由来し、全体の情報が部分系だけでは完結しないために周辺化が「不確実性」を生む。縮約密度行列は、観測可能な部分に対する統計の全てを集約し、相関の有無は ρ_A の性質だけでは決まらず、合成系側の構造と合わせて評価される。

3 時間発展と量子ダイナミクス

3.1 ユニタリ進化

3.1.1 リヤプノフ的な安定性の見方（概念整理）

密度行列の時間発展を考える際、数学的には写像（ダイナミクス）が状態空間上でどのような性質を持つかが論点になる。特に閉鎖系のユニタリ進化では、状態の距離やスペクトルといった性質が保たれるため、「安定性」を議論する際はリヤプノフ関数のような量が一定あるいは単調に変化するかを調べる、という観点がしばしば整理の枠組みとして使われる。ただし量子力学では観測による測定バックアクションや開放性が絡むと別の考え方が必要になる。

3.1.2 状態と密度行列の同時進化

閉鎖系でハミルトニアン H が与えられると、状態ベクトルは Schrödinger 方程式に従う。密度行列ではフォン・ノイマン方程式 iħ dρ/dt = [H,ρ] が対応し、解は ρ(t)=U(t)ρ(0)U†(t) と書ける。ユニタリが保つ性質により、純粋性の有無や固有値分布は変化しない。

3.2 非ユニタリ過程（開放系）

3.2.1 量子チャネルと完全正性

環境との相互作用がある場合、一般には ρ ↦ Φ(ρ) のような非ユニタリ写像が現れる。物理的妥当性のためには、写像 Φ が確率解釈を壊さないことが重要で、特に完全正性（任意の補助系を添えても正性が保たれる性質）が要件になる。これにより縮約や合成の操作と矛盾しない更新が保証される。トレース保存であれば、総確率が保存される。

3.2.2 マルコフ近似と動力学方程式の位置づけ

環境の相関時間が十分短いという仮定のもとでは、時間発展が時間局所的に記述される場合がある。マルコフ近似では、密度行列の時間微分が ρ に対する線形な作用素で与えられ、典型的にはリンドブラッド形の方程式が用いられる。ここでは「メモリ（履歴依存）」を無視したモデルとして位置づけられ、長時間挙動の記述や緩和過程の特徴づけに使われる。

3.3 代表的なノイズモデルの形式

3.3.1 ディフェージング型

ディフェージングは、位相情報が環境との相互作用により平均化されることで生じるノイズとして理解される。密度行列では基底選択に依存して非対角成分が減衰し、人口（対角成分）は必ずしも大きくは変わらない。結果として干渉は弱まり、状態は古典的な混合に近づく。

3.3.2 緩和型

緩和型ノイズはエネルギー準位の間で遷移が起こる形で現れ、平衡状態へ向かう緩和が特徴になる。一般には対角成分が時間とともに再配分され、最終的に環境温度や摂動条件に対応した定常状態へ収束する。緩和は振幅減衰や熱化として観測されることがある。

4 密度行列が担う情報量

4.1 エントロピー

4.1.1 ヴォン・ノイマンエントロピー

ヴォン・ノイマンエントロピーは S(ρ) = −Tr(ρ log ρ) により定義される。基底に依らず、固有値 λ_k に対して S(ρ) = −Σ_k λ_k log λ_k として計算できる。純粋状態では ρ の固有値が 1 と 0 になり、エントロピーは 0 になる。一方で混合が強まるほど値は増加し、状態の不確実性の尺度として用いられる。

4.1.2 純粋状態と混合状態の違い

純粋状態は位相を含む完全な情報が状態ベクトルに集約されているのに対し、混合状態は統計的な不確実性が残る。エントロピーはこの差を定量化し、同じ平均値でも混合度が異なるとエントロピーが変わることがある。さらに、ある測定の結果分布とエントロピーとの関係は、測定種別によっても姿を変える。

4.2 量子相関と指標

4.2.1 もつれの観点（縮約状態からの読み取り）

合成系で片側を縮約した状態が混合であることは、もつれの存在を示唆するが、それ単独では十分条件とはならない場合がある。直感的には、全体が純粋でも部分系が混合になる状況がもつれの典型である。もつれを判定・定量化するには、縮約密度行列だけでなく合成系側の情報（構造）を反映した指標が用いられる。

4.2.2 相関量の一般的な考え方

相関の測度は、単なる「一致度」だけでなく、古典的な相関と量子的な相関を分離して扱う必要がある。密度行列の枠組みでは、同一の周辺状態を持つ別の合成状態が存在しうるため、相関は合成系の演算子構造に依存する。したがって相関量の定義は、どの操作を許しどの性質を“量子的”と見なすかにより形が変わる。

4.3 状態の同定とトモグラフィ

4.3.1 測定からの復元の概要

量子状態トモグラフィは、複数の測定結果を集めて密度行列を再構成する手続きである。測定設定を十分に多様化し、期待値や確率のデータから線形方程式（あるいは推定問題）として ρ の成分を求める。データの有限性や誤差を考慮し、物理的制約（エルミート性、正性、トレース1）を満たす形で推定する方法が用いられる。

4.3.2 実験設計の基本方針

実験設計では、どの観測量（測定基底）を選ぶか、サンプル数をどう配分するかが精度を左右する。特に再構成したいサブスペースに対して十分な情報が得られるよう、測定点の配置を決めることが重要である。またノイズモデルを想定した推定（例えば不確かさの伝播やベイズ的更新）を組み合わせることで、現実の測定条件に合わせた密度行列の推定が可能になる。