相対エントロピーから情報と物理を組み立てる 第4部 第8章 §4.8.5 ★★ 熱的操作で自由に使えるのは熱状態とエネルギーを保つユニタリ。1回きりの仕事の上限は D 0 D_0 D 0 、作るコストは D ∞ D_\infty D ∞ 、多数回の平均で D D D に揃う。コヒーレンスの分は単独では仕事にならない
熱的操作
第2節の式 (4.8.7) で、熱浴 B の大きさと U ^ \hat U U ^ (全体のエネルギーを保つもの)を自由に選んでよいとした操作を熱的操作と呼ぶ。熱力学の資源理論は、熱的操作を「ただで行える操作」、熱状態を「ただで手に入る状態」とみなし、熱状態でない状態を資源として、どの状態からどの状態へ移れるか、仕事をどれだけ取り出せるかを問う1 1 M. Horodecki, J. Oppenheim (2013), Nature Commun. 4 , 2059; F. G. S. L. Brandão, M. Horodecki, J. Oppenheim, J. M. Renes, R. W. Spekkens (2013), Phys. Rev. Lett. 111 , 250404. 。熱的操作には2つの性質がある。
第一に、熱状態を動かさない(第2節)。したがって、データ処理不等式を満たすすべての量は熱的操作で増えない。D ( ρ ^ ∥ ρ ^ β ) D(\hat\rho\,\|\,\hat\rho_\beta) D ( ρ ^ ∥ ρ ^ β ) はその1つで、他に Rényi の相対エントロピー D α ( p ∥ q ) ≔ 1 α − 1 ln ∑ i p i α q i 1 − α D_\alpha(p\,\|\,q)\coloneqq\frac{1}{\alpha-1}\ln\sum_ip_i^\alpha q_i^{1-\alpha} D α ( p ∥ q ) : = α − 1 1 ln ∑ i p i α q i 1 − α (エネルギーの基底で対角な状態のとき)の族がある。α → 1 \alpha\to1 α → 1 で D D D に戻り、α \alpha α について増える関数で、両端は
D 0 ( p ∥ q ) = − ln ∑ i : p i > 0 q i , D ∞ ( p ∥ q ) = ln max i p i q i (4.8.19) D_0(p\,\|\,q)=-\ln\sum_{i:\,p_i>0}q_i,\qquad D_\infty(p\,\|\,q)=\ln\max_i\frac{p_i}{q_i}\tag{4.8.19} D 0 ( p ∥ q ) = − ln i : p i > 0 ∑ q i , D ∞ ( p ∥ q ) = ln i max q i p i ( 4.8.19 )
である。Brandão らは、この族のすべてが増えないことが、熱的操作での状態の移り変わりの必要条件になることを示した2 2 F. G. S. L. Brandão, M. Horodecki, N. H. Y. Ng, J. Oppenheim, S. Wehner (2015), Proc. Natl. Acad. Sci. USA 112 , 3275. 。マクロな熱力学では自由エネルギー1つが移り変わりを決めたが、1個の量子の系では「第二法則」が無数に並ぶ。
第二に、熱的操作は時間の並進と可換である。V ^ t = e − i H ^ S t / ℏ \hat V_t=\mathrm{e}^{-\mathrm{i}\hat H_{\mathrm S}t/\hbar} V ^ t = e − i H ^ S t /ℏ とすると、U ^ \hat U U ^ が全体の Hamiltonian と可換で ρ ^ B β \hat\rho_{\mathrm B}^\beta ρ ^ B β が H ^ B \hat H_{\mathrm B} H ^ B と可換なことからN t h ( V ^ t ρ ^ V ^ t † ) = V ^ t N t h ( ρ ^ ) V ^ t † \mathcal N_{\mathrm{th}}(\hat V_t\hat\rho\hat V_t^\dagger)=\hat V_t\mathcal N_{\mathrm{th}}(\hat\rho)\hat V_t^\dagger N th ( V ^ t ρ ^ V ^ t † ) = V ^ t N th ( ρ ^ ) V ^ t † が成り立つ。これを t t t について平均すると、エネルギーの縮退が無いとき N t h ∘ Δ B = Δ B ∘ N t h \mathcal N_{\mathrm{th}}\circ\Delta_{\mathcal B}=\Delta_{\mathcal B}\circ\mathcal N_{\mathrm{th}} N th ∘ Δ B = Δ B ∘ N th を得る。対角の状態からは対角の状態しか出てこないので、熱的操作はエネルギーの重ね合わせを作れない。
1回きりの仕事の上限は D 0 D_0 D 0
仕事を蓄える装置として、間隔 w w w の量子ビット W(バッテリー)を系に添える。熱的操作で ρ ^ S ⊗ ∣ 0 ⟩ ⟨ 0 ∣ W → σ ^ S ⊗ ∣ w ⟩ ⟨ w ∣ W \hat\rho_{\mathrm S}\otimes| 0 \rangle \langle 0 |_{\mathrm W}\to\hat\sigma_{\mathrm S}\otimes| w \rangle \langle w |_{\mathrm W} ρ ^ S ⊗ ∣0 ⟩ ⟨ 0 ∣ W → σ ^ S ⊗ ∣ w ⟩ ⟨ w ∣ W ができれば、仕事 w w w を確実に(確率1で)取り出したことになる。占有数 p p p の対角の状態から取り出せる w w w には次の上限がある。
命題 4.8.3(1回きりの仕事の上限 ★★)
熱的操作で p ⊗ ∣ 0 ⟩ ⟨ 0 ∣ W → σ ⊗ ∣ w ⟩ ⟨ w ∣ W p\otimes| 0 \rangle \langle 0 |_{\mathrm W}\to\sigma\otimes| w \rangle \langle w |_{\mathrm W} p ⊗ ∣0 ⟩ ⟨ 0 ∣ W → σ ⊗ ∣ w ⟩ ⟨ w ∣ W ができるなら、w ≤ k B T D 0 ( p ∥ γ ) w\le k_{\mathrm B} T\,D_0(p\,\|\,\gamma) w ≤ k B T D 0 ( p ∥ γ ) である。ここで γ \gamma γ は系の熱状態の占有数(Boltzmann 分布)。
証明 バッテリーの熱状態を γ W = ( 1 , e − β w ) / Z W \gamma_{\mathrm W}=(1,\mathrm{e}^{-\beta w})/Z_{\mathrm W} γ W = ( 1 , e − β w ) / Z W とする。D 0 D_0 D 0 も熱状態を動かさない操作で増えないので、D 0 ( p ⊗ ∣ 0 ⟩ ⟨ 0 ∣ ∥ γ ⊗ γ W ) ≥ D 0 ( σ ⊗ ∣ w ⟩ ⟨ w ∣ ∥ γ ⊗ γ W ) D_0(p\otimes| 0 \rangle \langle 0 |\,\|\,\gamma\otimes\gamma_{\mathrm W})\ge D_0(\sigma\otimes| w \rangle \langle w |\,\|\,\gamma\otimes\gamma_{\mathrm W}) D 0 ( p ⊗ ∣0 ⟩ ⟨ 0∣ ∥ γ ⊗ γ W ) ≥ D 0 ( σ ⊗ ∣ w ⟩ ⟨ w ∣ ∥ γ ⊗ γ W ) が成り立つ。D 0 D_0 D 0 は積の状態について和に分かれ、D 0 ( ∣ 0 ⟩ ⟨ 0 ∣ ∥ γ W ) = ln Z W D_0(| 0 \rangle \langle 0 |\,\|\,\gamma_{\mathrm W})=\ln Z_{\mathrm W} D 0 ( ∣0 ⟩ ⟨ 0∣ ∥ γ W ) = ln Z W 、D 0 ( ∣ w ⟩ ⟨ w ∣ ∥ γ W ) = β w + ln Z W D_0(| w \rangle \langle w |\,\|\,\gamma_{\mathrm W})=\beta w+\ln Z_{\mathrm W} D 0 ( ∣ w ⟩ ⟨ w ∣ ∥ γ W ) = β w + ln Z W なので、D 0 ( p ∥ γ ) ≥ D 0 ( σ ∥ γ ) + β w ≥ β w D_0(p\,\|\,\gamma)\ge D_0(\sigma\,\|\,\gamma)+\beta w\ge\beta w D 0 ( p ∥ γ ) ≥ D 0 ( σ ∥ γ ) + β w ≥ β w を得る。
Horodecki と Oppenheim は、対角の状態ではこの上限に届くこと、逆に熱状態から状態 p p p を作るのに要る仕事が k B T D ∞ ( p ∥ γ ) k_{\mathrm B} T\,D_\infty(p\,\|\,\gamma) k B T D ∞ ( p ∥ γ ) であることを示した。D 0 ≤ D ≤ D ∞ D_0\le D\le D_\infty D 0 ≤ D ≤ D ∞ なので、1個の系では取り出せる仕事と作るのに要る仕事が食い違い、移り変わりは可逆でない。同じ状態を n n n 個並べると、1個あたりの取り出せる仕事と作るコストはどちらも k B T D ( p ∥ γ ) k_{\mathrm B} T\,D(p\,\|\,\gamma) k B T D ( p ∥ γ ) に近づく(Brandão ほか 2013)。これは 第1部第8章 で見た大偏差の考え方と同じ構造で、多数の複製では典型的な系列だけが残り、D D D だけが意味を持つ。
対角の状態どうしの移り変わりには、図で判定できる必要十分条件がある。p i / γ i p_i/\gamma_i p i / γ i の大きい順に準位を並べ、点 ( ∑ j ≤ i γ j , ∑ j ≤ i p j ) (\sum_{j\le i}\gamma_j,\sum_{j\le i}p_j) ( ∑ j ≤ i γ j , ∑ j ≤ i p j ) を結んだ折れ線を p p p の熱的 Lorenz 曲線と呼ぶ。熱的操作で p → q p\to q p → q ができるのは、p p p の曲線がどこでも q q q の曲線より上(または同じ)にあるときに限る(熱的多数化。Horodecki–Oppenheim 2013)。熱状態の曲線は対角線で、どの状態の曲線もその上にある。D 0 D_0 D 0 は曲線が高さ 1 1 1 に届く横座標の − ln -\ln − ln 、D ∞ D_\infty D ∞ は最初の線分の傾きの ln \ln ln として図から読める。
図 4.8.3 β ℏ ω = 1 \beta\hbar\omega=1 β ℏ ω = 1 の量子ビットの熱的 Lorenz 曲線。熱状態(対角線)・最大混合 ( 1 2 , 1 2 ) (\frac12,\frac12) ( 2 1 , 2 1 ) (実線)・基底状態(破線)。基底状態の曲線は最大混合の曲線より上にあるので、基底状態から最大混合へは熱的操作で移れるが、逆はできない。最大混合の曲線は横座標 1 1 1 でやっと高さ 1 1 1 に届くので D 0 = 0 D_0=0 D 0 = 0
コヒーレンスの分は単独では仕事にならない
第1節の式 (4.8.4) で、D D D は占有数の分とコヒーレンスの分に分かれた。熱的操作では、コヒーレンスの分を仕事に変えることができない。実際、ρ ^ ⊗ ∣ 0 ⟩ ⟨ 0 ∣ W \hat\rho\otimes| 0 \rangle \langle 0 |_{\mathrm W} ρ ^ ⊗ ∣0 ⟩ ⟨ 0 ∣ W から対角の目標 σ ^ ⊗ ∣ w ⟩ ⟨ w ∣ W \hat\sigma\otimes| w \rangle \langle w |_{\mathrm W} σ ^ ⊗ ∣ w ⟩ ⟨ w ∣ W が作れたとすると、熱的操作は Δ B \Delta_{\mathcal B} Δ B と可換なので、同じ操作を Δ B ( ρ ^ ) ⊗ ∣ 0 ⟩ ⟨ 0 ∣ W \Delta_{\mathcal B}(\hat\rho)\otimes| 0 \rangle \langle 0 |_{\mathrm W} Δ B ( ρ ^ ) ⊗ ∣0 ⟩ ⟨ 0 ∣ W に当てると Δ B ( σ ^ ⊗ ∣ w ⟩ ⟨ w ∣ ) = σ ^ ⊗ ∣ w ⟩ ⟨ w ∣ \Delta_{\mathcal B}(\hat\sigma\otimes| w \rangle \langle w |)=\hat\sigma\otimes| w \rangle \langle w | Δ B ( σ ^ ⊗ ∣ w ⟩ ⟨ w ∣ ) = σ ^ ⊗ ∣ w ⟩ ⟨ w ∣ が出る。したがって ρ ^ \hat\rho ρ ^ から確実に取り出せる仕事は、非対角を落とした Δ B ( ρ ^ ) \Delta_{\mathcal B}(\hat\rho) Δ B ( ρ ^ ) から取り出せる仕事を超えない。コヒーレンスの分を使うには、外から位相の基準となるコヒーレントな駆動(レーザーや古典的な場)を持ち込む必要がある。Brandão らの 2013 年の結果でも、n n n 個の複製から k B T D k_{\mathrm B} T\,D k B T D を取り出すには、n n n に比べて少ない量のエネルギーの重ね合わせを外から借りることが前提になっている。
例題 4.8.5 1回きりの仕事と作るコスト ★★
オリジナル β ℏ ω = 1 \beta\hbar\omega=1 β ℏ ω = 1 の量子ビットについて考える(Boltzmann 分布 γ = ( 0.7311 , 0.2689 ) \gamma=(0.7311,0.2689) γ = ( 0.7311 , 0.2689 ) )。
(1) 最大混合 p = ( 1 2 , 1 2 ) p=(\frac12,\frac12) p = ( 2 1 , 2 1 ) の D 0 D_0 D 0 、D D D 、D ∞ D_\infty D ∞ を求め、1回きりに取り出せる仕事・平均で取り出せる仕事・作るのに要る仕事を比べよ。
(2) 基底状態 p = ( 1 , 0 ) p=(1,0) p = ( 1 , 0 ) について同じことをせよ。
(3) 例題 4.8.1 の ∣ + ⟩ | + \rangle ∣ + ⟩ から熱的操作だけで確実に取り出せる仕事はいくらか。
解答 (1) 両方の準位に確率があるので D 0 = − ln ( 0.7311 + 0.2689 ) = 0 D_0=-\ln(0.7311+0.2689)=0 D 0 = − ln ( 0.7311 + 0.2689 ) = 0 である。D = 0.1201 D=0.1201 D = 0.1201 、D ∞ = ln max ( 0.5 / 0.7311 , 0.5 / 0.2689 ) = ln 1.859 = 0.620 D_\infty=\ln\max(0.5/0.7311,\,0.5/0.2689)=\ln1.859=0.620 D ∞ = ln max ( 0.5/0.7311 , 0.5/0.2689 ) = ln 1.859 = 0.620 である。1回きりには仕事を1つも確実に取り出せないが、多数の複製なら1個あたり 0.120 k B T 0.120\, k_{\mathrm B} T 0.120 k B T まで取り出せる。逆にこの状態を熱状態から1個作るには 0.620 k B T 0.620\, k_{\mathrm B} T 0.620 k B T が要る。図 4.8.3 の実線は横座標 1 1 1 で高さ 1 1 1 に届き(D 0 = 0 D_0=0 D 0 = 0 )、最初の線分の傾きは 0.5 / 0.2689 = 1.859 0.5/0.2689=1.859 0.5/0.2689 = 1.859 である。
(2) 準位が1つだけ占有されているので D 0 = − ln 0.7311 = 0.3133 D_0=-\ln0.7311=0.3133 D 0 = − ln 0.7311 = 0.3133 、D = 1 ⋅ ln ( 1 / 0.7311 ) = 0.3133 D=1\cdot\ln(1/0.7311)=0.3133 D = 1 ⋅ ln ( 1/0.7311 ) = 0.3133 、D ∞ = ln ( 1 / 0.7311 ) = 0.3133 D_\infty=\ln(1/0.7311)=0.3133 D ∞ = ln ( 1/0.7311 ) = 0.3133 で、3つが一致する。純粋なエネルギー固有状態では、取り出す仕事と作るコストが同じで、1個でも可逆である。
(3) 熱的操作は Δ B \Delta_{\mathcal B} Δ B と可換なので、∣ + ⟩ | + \rangle ∣ + ⟩ から確実に取り出せる仕事は Δ B ( ∣ + ⟩ ⟨ + ∣ ) = 1 / 2 \Delta_{\mathcal B}(| + \rangle \langle + |)=\mathbf 1/2 Δ B ( ∣ + ⟩ ⟨ + ∣ ) = 1 /2 から取り出せる仕事を超えず、(1) により 0 0 0 である。∣ + ⟩ | + \rangle ∣ + ⟩ の D = 0.8133 D=0.8133 D = 0.8133 のうちコヒーレンスの分 ln 2 \ln2 ln 2 は、熱浴だけでは使えない。ergotropy 0.5 k B T 0.5\, k_{\mathrm B} T 0.5 k B T を取り出すには、∣ + ⟩ | + \rangle ∣ + ⟩ を ∣ 0 ⟩ | 0 \rangle ∣0 ⟩ へ回すコヒーレントな駆動が要る。
1 M. Horodecki, J. Oppenheim (2013), Nature Commun. 4 , 2059; F. G. S. L. Brandão, M. Horodecki, J. Oppenheim, J. M. Renes, R. W. Spekkens (2013), Phys. Rev. Lett. 111 , 250404.
2 F. G. S. L. Brandão, M. Horodecki, N. H. Y. Ng, J. Oppenheim, S. Wehner (2015), Proc. Natl. Acad. Sci. USA 112 , 3275.
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com