密度行列

英

density matrix

Wikipedia preview

出典(authority):フリー百科事典『ウィキペディア（Wikipedia）』「2015/04/16 15:40:41」(JST)

wiki ja

[Wiki ja表示]

密度行列（みつどぎょうれつ、英語: density matrix）は、量子力学における混合状態を表現するために使われる行列である。

定義

ある系（ここでは時間依存は考えない）の状態が、|Ψ₁〉, |Ψ₂〉, ..., |Ψ_k〉, ... という状態ベクトルが古典的に混ざったものとして記述できるとする。つまり、ある系のどんなオブザーバブルの測定値の確率分布も、|Ψ₁〉, |Ψ₂〉, ..., |Ψ_k〉, ... での測定値の確率分布に重みをつけて平均したものとして表せるとする。

k 番目の状態が重み p_k で混合されているとき、

で定義される演算子 を密度演算子 (density operator ) と言う。密度行列 (density matrix ) は、密度演算子を行列表示したものである。

ここで、|Ψ_k〉 は互いに直交している必要はない。たとえば、スピンの z 成分の +1 の固有状態と、y 成分の −1 の固有状態を混ぜることも可能である。

このような複数の状態が混ざってできた状態を混合状態 (mixed state ) といい、そうでない状態を純粋状態 (pure state ) という。

性質

一般に密度演算子 はエルミート演算子で、ヒルベルト空間上の任意の状態ベクトル に対し、

を満たす。ここで Tr{ρ} は密度行列 ρ のトレース (trace) である。

一般には、密度演算子の二乗は、

$\begin{alignat}{4} \rho^2 \, & = \bigl\{ \sum_k p_k | \Psi_k \rangle \langle \Psi_k | \bigr\} \bigl\{ \sum_{k'} p_{k'} | \Psi_{k'} \rangle \langle \Psi_{k'} | \bigr\} \\ & = \sum_{k} \sum_{k'} p_k p_{k'} | \Psi_k \rangle \langle \Psi_k | \Psi_{k'} \rangle \langle \Psi_{k'} | \\ & = \sum_k p_k^2 | \Psi_k \rangle \langle \Psi_k | \\ & \ne \rho \end{alignat}$

となる。ここで、簡単のために状態 |Ψ_k〉 と状態 |Ψ_k'〉 は規格直交性 〈Ψ_k|Ψ_k'〉 = δ_k,k' を持つとした（密度行列を対角化することで、このような表示は常に可能である）。

特別な場合として、複数の状態が混ざっていない純粋状態 |Ψ〉 に対する密度行列 ρ = |Ψ〉〈Ψ| では、状態ベクトルの規格化条件 〈Ψ|Ψ〉 = 1 より、

となる。よって、純粋状態の場合は ρⁿ = ρ（n は正の整数）であることが分かる。また、純粋状態に対する密度行列は射影演算子の形となっている。

物理量の期待値

量子力学において、純粋状態 |Ψ_k〉 におけるオブザーバブル の期待値は、

となる。この純粋状態が混ざってできた混合状態の期待値 は、重み p_k を使って、

$\left\langle A \right\rangle = \sum_k p_k \langle \Psi_k | \hat{A} | \Psi_k \rangle = \operatorname{Tr} \{\rho A\}$

つまり密度行列 ρ と、オブザーバブル の行列表示 A との積の、トレースで表される。

時間発展

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

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

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

$\begin{align} i \hbar {\partial |\Psi_k \rangle \over {\partial t}} &= \hat{H} | \Psi_k \rangle,\\ - i \hbar {\partial \langle \Psi_k| \over {\partial t}} &= \langle \Psi_k | \hat{H}, \end{align}$

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

統計力学への応用

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

密度行列 ρ は、たとえばカノニカル分布では、

$\mathbf{\rho} = \frac{\mathrm{e}^{-\beta H}}{\operatorname{Tr} (\mathrm{e}^{-\beta H}) }$

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

$\rho = \frac{\mathrm{e}^{-\beta H_\mathrm{G}}}{\operatorname{Tr} (\mathrm{e}^{-\beta H_\mathrm{G}}) } = \mathrm{e}^{\beta(\Omega - H_\mathrm{G})}$

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

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

$\langle A \rangle = \operatorname{Tr}\{\rho A\} = \frac{\operatorname{Tr}\{\mathrm{e}^{-\beta H}A\}}{\operatorname{Tr}\{\mathrm{e}^{-\beta H}\}}$

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

$\langle \operatorname{Id} \rangle = \frac{\operatorname{Tr}\{\mathrm{e}^{-\beta H}\operatorname{Id}\}}{\operatorname{Tr}\{\mathrm{e}^{-\beta H}\}} = \frac{\operatorname{Tr}\{\mathrm{e}^{-\beta H}\}}{\operatorname{Tr}\{\mathrm{e}^{-\beta H}\}} = 1$

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

$\langle H \rangle = \frac{\operatorname{Tr}\{\mathrm{e}^{-\beta H}H\}}{\operatorname{Tr}\{\mathrm{e}^{-\beta H}\}} = \frac{\sum_i E_i\mathrm{e}^{-\beta E_i}}{\sum_i \mathrm{e}^{-\beta E_i}}$

と書き換えられる。

関連記事

• 量子状態
• 統計力学
  • 量子統計力学
• リウヴィルの定理 (物理学)

UpToDate Contents

全文を閲覧するには購読必要です。 To read the full text you will need to subscribe.

• 1. リポ蛋白の分類、代謝、およびアテローム性動脈硬化症における役割 lipoprotein classification metabolism and role in atherosclerosis
• 2. 閉経後の女性における骨粗鬆症のマネージメントの概要 overview of the management of osteoporosis in postmenopausal women
• 3. 肝線維症の病因 pathogenesis of hepatic fibrosis
• 4. 不安定プラークの病理および病因 pathology and pathogenesis of the vulnerable plaque
• 5. 大腸癌の病理および予後決定因子 pathology and prognostic determinants of colorectal cancer

Japanese Journal

• 並列FMOプログラムOpenFMOの性能最適化

• 稲富 雄一,眞木 淳,高見 利也,本田 宏明,小林 泰三,南里 豪志,青柳 睦,南 一生
• 情報処理学会研究報告. [ハイパフォーマンスコンピューティング] 2012-HPC-133(35), 1-8, 2012-03-19
• … 指して,並列フラグメント分子軌道プログラム OpenFMO の性能最適化を行った.その際に,FMO 計算のカーネルコードである分子積分計算の負荷均等化と,中間データであるモノマー密度行列データの保存,アクセス方法に注目して最適化を行った.その結果,分子積分計算の並列化効率向上と,モノマー密度行列へのアクセス性能向上を達成することができ,数万並列実行時での効率的な FMO 計算が可能になった. …
• NAID 110008803070

• 並列密度行列繰りこみ群法の開発と量子物理学への応用(<小特集>計算物理)

• 町田 昌彦,奥村 雅彦,山田 進
• シミュレーション 31(1), 35-44, 2012-03-15
• NAID 110009437352

• Duality of Linearity and Nonlinearity to Reveal Irreversibility in Transport Phenomena (量子科学における双対性とスケール)

• 鈴木 増雄
• 素粒子論研究 119(4A), 145-157, 2012-02-20
• … 前回は,電気伝導という典型的な輸送現象におけるエントロピー生成が密度行列の対称成分ρ_s(t)(=ρ_0+ρ_2(t)+…)によって記述されることを見つけ,この事実に基づいてフォン・ノイマン方程式を対称成分ρ_s(t)と反対称成分ρ_a(t))=ρ_1(t)+ρ_3(t)+…とに分離し,対称成分の時間発展方程式に糸の内部で発生した熱を外部に取り出す緩和項をつけ加える定式化を行い,定常状態を記述することに成功した。 …
• NAID 110009418987

Related Links

• 密度行列とは - goo Wikipedia (ウィキペディア)

つまり密度行列 ρ と、オブザーバブル の行列表示 A との積の、トレースで表される。 エントロピー [編集] 混合状態のフォン・ノイマンエントロピー は、 の固有値の項、あるいは、密度演算子 のトレースと 行列の対数 （英語版 ...

• 第4章 密度行列入門 - インターネット回線 | au

4. 1 密度行列の時間変化 43 小さなH1(t)からなっているときには，以下のようにρ，H1を変換する． *( ) exp(0) ( )exp( ) iH t iH t ρt=ρt == −0 (4.1.17) * 0 11( ) exp( ) exp( ) iH t iH t Ht=H == −0 (4.1.18) このような表現を相互作用表示という ...

Related Pictures

---

★リンクテーブル★

関連記事

「行列」「列」

「行列」

　　[★]

英

matrix

• (数学)matrix、matrices

関

基質、マトリックス、マトリクス

「列」

　　[★]

英

row、tier、train

関

作条、訓練、層、連、段