相対エントロピーから情報と物理を組み立てる 第3部 第2章 §3.2.6 ★★ 量子統計力学の S=−kBTrρ^lnρ^、量子情報の講義の底2の書き方、統計力学の Gibbs–Bogoliubov の不等式を、この章の式と並べる
量子統計力学のエントロピー
量子統計力学の授業では、エントロピーを
S=−kBTr(ρ^lnρ^)
と kB を付けて定義し、カノニカル分布 ρ^=e−βH^/Z に代入して S=(U−F)/T、F=−kBTlnZ を導くことが多い。この章の von Neumann エントロピーは無次元(nat)で、熱力学のエントロピーとは
S熱=kBS(ρ^β)(3.2.29)
の関係にある。kB は nat を J/K に換える係数で、情報の量としての中身は S(ρ^β) が持っている。授業の式 TS熱=U−F は、式 (3.2.24) で ρ^=ρ^β としたものである。
量子情報の講義の書き方
量子情報の講義や教科書では、エントロピーを底2の対数で S(ρ)=−Trρlog2ρ(単位ビット)と書き、量子相対エントロピーを S(ρ∥σ)=Trρlog2ρ−Trρlog2σ と書くことが多い。この巻では、古典と量子で同じ記号 D を使い、自然対数で書く。換算はビットの値 = nat の値 /ln2 である。記号 S(ρ∥σ) を使わない理由は、S がすでに作用とエントロピーに使われていることと、古典の D(p∥q) と量子の D(ρ^∥σ^) が命題 3.2.2 で直接つながることを記号の上でも見せるためである。また講義では Klein の不等式を「ln が作用素単調である」ことや Lagrange の未定乗数法で示すことがあるが、この章の証明(二重確率行列と Jensen)は古典の Gibbs の不等式の証明をそのまま使っている。
統計力学の Gibbs–Bogoliubov の不等式
統計力学の変分法では、解ける Hamiltonian H^0 を選び、その平衡状態 ρ^0=e−βH^0/Z0 を使って
F≤F0+⟨H^−H^0⟩0(3.2.30)
を示す(Gibbs–Bogoliubov の不等式、Feynman の変分原理とも呼ぶ)。ここで F0=−kBTlnZ0、⟨⋅⟩0=Tr(ρ^0⋅) である。これは定理 3.2.12 で試行状態を ρ^0 に選んだものである。実際、S(ρ^0)=β(⟨H^0⟩0−F0) なので
F(ρ^0)=⟨H^⟩0−kBTS(ρ^0)=F0+⟨H^−H^0⟩0,
で、F(ρ^0)−F=kBTD(ρ^0∥ρ^β) が上界と真の値の差になる。授業で平均場近似を「1体の有効 Hamiltonian H^0=−∑iheffZi を選んでheff について右辺を最小にする」と習った場合、これは 3.2.5節の積状態の族の中での最小化と同じ問題である。
例題 3.2.7 2準位系の集まりのエントロピー ★★
オリジナル例題 3.2.5 の2準位系(βE=ln3)が 1mol あり、互いに独立とする。
- 1つの系の S(ρ^β) を求めよ。
- 統計力学の授業の式 S熱=NkB[ln(1+e−βE)+1+e−βEβEe−βE] と式 (3.2.29) が一致することを確かめ、値を JK−1mol−1 で求めよ(気体定数 R=NAkB=8.314JK−1mol−1)。
- 高温の極限の値を求め、ビットで言い直せ。
解答
(1) 固有値は (3/4,1/4) なので S(ρ^β)=hb(1/4)=−43ln43−41ln41=0.5623 nat。
(2) e−βE=1/3 を入れると括弧の中は ln34+4/3ln3⋅(1/3)=ln34+41ln3。一方 hb(1/4)=43ln34+41ln4=ln34+41ln3(41ln4−41ln34=41ln3)で一致する。独立な系のエントロピーは加法的(命題 3.2.5(b))なのでS熱=NAkBhb(1/4)=8.314×0.5623=4.675JK−1mol−1。
(3) βE→0 で ρ^β→1/2、S→ln2 なので S熱→Rln2=5.763JK−1mol−1。1つの系あたり ln2 nat =1 ビットである。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com