相対エントロピーから情報と物理を組み立てる 第4部 第8章 §4.8.3 ★★★ 消去で熱浴へ流れる熱の下限 − k B T Δ S S - k_{\mathrm B} T\Delta S_{\mathrm S} − k B T Δ S S は前の節の等式の書き直しで、余分は相関と熱浴のずれ。段に分けた消去の余分は D D D の和
等式の形の Landauer 原理
定理 4.8.2 を、熱浴へ流れる熱 Q B ≔ − Q Q_{\mathrm B}\coloneqq-Q Q B : = − Q について解くと
β Q B = − Δ S S + I ( S ; B ) ρ ′ + D ( ρ ^ B ′ ∥ ρ ^ B β ) ≥ − Δ S S (4.8.11) \beta Q_{\mathrm B}=-\Delta S_{\mathrm S}+I(\mathrm S;\mathrm B)_{\rho'}+D(\hat\rho'_{\mathrm B}\,\|\,\hat\rho^\beta_{\mathrm B})\ \ge\ -\Delta S_{\mathrm S}\tag{4.8.11} β Q B = − Δ S S + I ( S ; B ) ρ ′ + D ( ρ ^ B ′ ∥ ρ ^ B β ) ≥ − Δ S S ( 4.8.11 )
となる。系のエントロピーを減らす操作(消去)では、熱浴へ少なくとも − k B T Δ S S - k_{\mathrm B} T\Delta S_{\mathrm S} − k B T Δ S S の熱が流れる。1ビットの記憶(最大混合の量子ビット)を基底状態へ消すとΔ S S = − ln 2 \Delta S_{\mathrm S}=-\ln2 Δ S S = − ln 2 なので Q B ≥ k B T ln 2 Q_{\mathrm B}\ge k_{\mathrm B} T\ln2 Q B ≥ k B T ln 2 である。これが Landauer 原理1 1 R. Landauer (1961), IBM J. Res. Dev. 5 , 183. で、式 (4.8.11) はそれを等式の形で書いたものになっている。Reeb と Wolf はこの等式から出発し、熱浴が有限の次元を持つときは下限に等号で届かず、熱浴の大きさで決まる正の補正が必ず残ることを示した2 2 D. Reeb, M. M. Wolf (2014), New J. Phys. 16 , 103011. 。等号には I = 0 I=0 I = 0 と D ( ρ ^ B ′ ∥ ρ ^ B β ) = 0 D(\hat\rho'_{\mathrm B}\,\|\,\hat\rho^\beta_{\mathrm B})=0 D ( ρ ^ B ′ ∥ ρ ^ B β ) = 0 の両方、つまり熱浴が系と相関を持たず、しかも平衡からずれないことが要る。
消去の操作には仕事も要る。系の Hamiltonian が始めと終わりで同じなら、系がされる仕事は W = Δ E S + Q B W=\Delta E_{\mathrm S}+Q_{\mathrm B} W = Δ E S + Q B (系と熱浴のエネルギーの増えの和)で、式 (4.8.11) を使うと W − Δ F = k B T [ I + D ( ρ ^ B ′ ∥ ρ ^ B β ) ] W-\Delta F= k_{\mathrm B} T[I+D(\hat\rho'_{\mathrm B}\,\|\,\hat\rho^\beta_{\mathrm B})] W − Δ F = k B T [ I + D ( ρ ^ B ′ ∥ ρ ^ B β )] となる。ここで Δ F = Δ E S − k B T Δ S S \Delta F=\Delta E_{\mathrm S}- k_{\mathrm B} T\Delta S_{\mathrm S} Δ F = Δ E S − k B T Δ S S は非平衡の自由エネルギーの変化である。消去と測定のどちらにこのコストを割り振るかは 第4部第5章 の主題で、Maxwell の悪魔の帳尻は 第4部第4章 で扱う。
熱浴の量子ビット1個で消す
最も簡単な消去は、熱浴の量子ビットと状態を入れ替えることである。系は縮退した量子ビット(H ^ S = 0 \hat H_{\mathrm S}=0 H ^ S = 0 )で最大混合状態 1 / 2 \mathbf 1/2 1 /2 にあり、熱浴は間隔 ℏ ω B \hbar\omega_{\mathrm B} ℏ ω B の量子ビットで熱状態にあるとする。x ≔ β ℏ ω B x\coloneqq\beta\hbar\omega_{\mathrm B} x : = β ℏ ω B 、p 1 ( x ) ≔ 1 / ( 1 + e x ) p_1(x)\coloneqq1/(1+\mathrm{e}^{x}) p 1 ( x ) : = 1/ ( 1 + e x ) (熱状態の励起確率)と書く。SWAP の後、系は熱浴の熱状態(励起確率 p 1 ( x ) p_1(x) p 1 ( x ) )に、熱浴は最大混合状態になる。熱浴のエネルギーの増えは Q B = ℏ ω B ( 1 2 − p 1 ( x ) ) Q_{\mathrm B}=\hbar\omega_{\mathrm B}(\frac12-p_1(x)) Q B = ℏ ω B ( 2 1 − p 1 ( x )) 、系のエントロピーの変化は Δ S S = h b ( p 1 ( x ) ) − ln 2 \Delta S_{\mathrm S}=h_{\mathrm b}(p_1(x))-\ln2 Δ S S = h b ( p 1 ( x )) − ln 2 、相関は残らないので I = 0 I=0 I = 0 である。式 (4.8.11) の余りは
β Q B + Δ S S = D ( ( 1 2 , 1 2 ) ∥ ( 1 − p 1 ( x ) , p 1 ( x ) ) ) (4.8.12) \beta Q_{\mathrm B}+\Delta S_{\mathrm S}=D\Bigl(\bigl(\tfrac12,\tfrac12\bigr)\,\Big\|\,\bigl(1-p_1(x),\,p_1(x)\bigr)\Bigr)\tag{4.8.12} β Q B + Δ S S = D ( ( 2 1 , 2 1 ) ( 1 − p 1 ( x ) , p 1 ( x ) ) ) ( 4.8.12 )
である。x x x を大きくすると系はほぼ基底状態へ消えるが、β Q B ≈ x / 2 \beta Q_{\mathrm B}\approx x/2 β Q B ≈ x /2 が際限なく増え、余りも x / 2 − ln 2 x/2-\ln2 x /2 − ln 2 のように増える。完全な消去を1回の入れ替えで行うと、熱浴へ流す熱はいくらでも大きくなる。
段に分けた消去
下限に近づくには、操作を段に分ける。系を間隔 ℏ ω \hbar\omega ℏ ω が動かせる量子ビットにし、x = β ℏ ω x=\beta\hbar\omega x = β ℏ ω を 0 = x 0 < x 1 < ⋯ < x N = X 0=x_0<x_1<\cdots<x_N=X 0 = x 0 < x 1 < ⋯ < x N = X と N N N 段で上げる。各段では、まず間隔を x k − 1 x_{k-1} x k − 1 から x k x_k x k へ一瞬で上げ(状態は変わらないので、式 (4.8.1) の仕事は k B T p 1 ( x k − 1 ) ( x k − x k − 1 ) k_{\mathrm B} T\,p_1(x_{k-1})(x_k-x_{k-1}) k B T p 1 ( x k − 1 ) ( x k − x k − 1 ) )、次に同じ間隔 x k x_k x k の新しい熱浴の量子ビットと SWAP する(エネルギーを保つので仕事は 0 0 0 )。SWAP の後、系は熱状態 π x k \pi_{x_k} π x k (励起確率 p 1 ( x k ) p_1(x_k) p 1 ( x k ) )になる。最後に間隔を一瞬で 0 0 0 に戻す(仕事 − k B T p 1 ( X ) X - k_{\mathrm B} T\,p_1(X)X − k B T p 1 ( X ) X )と、系は縮退した量子ビットで励起確率 p 1 ( X ) p_1(X) p 1 ( X ) の状態にあり、消去の誤りは p 1 ( X ) p_1(X) p 1 ( X ) である。
一瞬で間隔を変える段では状態が変わらないので、非平衡の自由エネルギーの変化と仕事は等しい。SWAP の段では仕事が 0 0 0 で、自由エネルギーは式 (4.8.9) により k B T D ( π x k − 1 ∥ π x k ) k_{\mathrm B} T\,D(\pi_{x_{k-1}}\,\|\,\pi_{x_k}) k B T D ( π x k − 1 ∥ π x k ) だけ減る。これを足し合わせると、全体の仕事と自由エネルギーの変化 Δ F = k B T [ ln 2 − h b ( p 1 ( X ) ) ] \Delta F= k_{\mathrm B} T[\ln2-h_{\mathrm b}(p_1(X))] Δ F = k B T [ ln 2 − h b ( p 1 ( X ))] の差は
W N − Δ F = k B T ∑ k = 1 N D ( π x k − 1 ∥ π x k ) (4.8.13) W_N-\Delta F= k_{\mathrm B} T\sum_{k=1}^ND(\pi_{x_{k-1}}\,\|\,\pi_{x_k})\tag{4.8.13} W N − Δ F = k B T k = 1 ∑ N D ( π x k − 1 ∥ π x k ) ( 4.8.13 )
となる。これは 第4部第1章 で古典のばねについて出した「段に分けた操作の余分は D D D の和」と同じ形である。段を細かくすると各段の D D D は2次の小ささになり、D ( π x ∥ π x + δ ) ≈ 1 2 J ( x ) δ 2 D(\pi_x\,\|\,\pi_{x+\delta})\approx\frac12J(x)\delta^2 D ( π x ∥ π x + δ ) ≈ 2 1 J ( x ) δ 2 と Fisher 情報 J ( x ) = p 1 ( x ) ( 1 − p 1 ( x ) ) = − p 1 ′ ( x ) J(x)=p_1(x)(1-p_1(x))=-p_1'(x) J ( x ) = p 1 ( x ) ( 1 − p 1 ( x )) = − p 1 ′ ( x ) (第1部第7章 )で書ける。等間隔の N N N 段では δ = X / N \delta=X/N δ = X / N で
W N − Δ F ≈ k B T δ 2 ∫ 0 X J ( x ) d x = k B T X 2 N ( 1 2 − p 1 ( X ) ) (4.8.14) W_N-\Delta F\approx k_{\mathrm B} T\frac{\delta}{2}\int_0^XJ(x)\,\mathrm{d} x= k_{\mathrm B} T\,\frac{X}{2N}\Bigl(\frac12-p_1(X)\Bigr)\tag{4.8.14} W N − Δ F ≈ k B T 2 δ ∫ 0 X J ( x ) d x = k B T 2 N X ( 2 1 − p 1 ( X ) ) ( 4.8.14 )
となり、N → ∞ N\to\infty N → ∞ で余分は 0 0 0 になる。準静的な極限では W ∞ = k B T [ ∫ 0 X p 1 d x − X p 1 ( X ) ] = k B T [ ln 2 − h b ( p 1 ( X ) ) ] = Δ F W_\infty= k_{\mathrm B} T[\int_0^Xp_1\,\mathrm{d} x-Xp_1(X)]= k_{\mathrm B} T[\ln2-h_{\mathrm b}(p_1(X))]=\Delta F W ∞ = k B T [ ∫ 0 X p 1 d x − X p 1 ( X )] = k B T [ ln 2 − h b ( p 1 ( X ))] = Δ F で、X → ∞ X\to\infty X → ∞ で k B T ln 2 k_{\mathrm B} T\ln2 k B T ln 2 に届く。有限の段数と有限の X X X では、消去の誤りと余分な仕事の両方が残る。
図 4.8.1 量子ビットの段に分けた消去(X = 5 X=5 X = 5 )で余分な仕事 W N − Δ F W_N-\Delta F W N − Δ F を段の数 N N N に対して描いた。点線は式 (4.8.13) の厳密な値、実線は Fisher 情報による近似 (4.8.14) 。N = 1 N=1 N = 1 は熱浴の量子ビット1個との SWAP と同じ値 1.814 1.814 1.814 で、段を増やすと 1 / N 1/N 1/ N で減る
例題 4.8.3 段に分けた量子ビットの消去 ★★★
オリジナル 最大混合状態の縮退した量子ビットを、間隔を x = 0 x=0 x = 0 から X = 5 X=5 X = 5 まで上げる段の消去で消す(温度 T T T )。
(1) N = 1 N=1 N = 1 のとき、仕事 W 1 W_1 W 1 、Δ F \Delta F Δ F 、余分を k B T k_{\mathrm B} T k B T 単位で求め、余分が D ( π 0 ∥ π 5 ) D(\pi_0\,\|\,\pi_5) D ( π 0 ∥ π 5 ) に等しいことを確かめよ。
(2) N = 2 N=2 N = 2 と N = 10 N=10 N = 10 の余分を求め、式 (4.8.14) の近似と比べよ。
解答 (1) p 1 ( 0 ) = 1 2 p_1(0)=\frac12 p 1 ( 0 ) = 2 1 、p 1 ( 5 ) = 1 / ( 1 + e 5 ) = 0.006693 p_1(5)=1/(1+\mathrm{e}^5)=0.006693 p 1 ( 5 ) = 1/ ( 1 + e 5 ) = 0.006693 である。仕事は間隔を上げる段の 1 2 × 5 = 2.5 \frac12\times5=2.5 2 1 × 5 = 2.5 と、戻す段の − 0.006693 × 5 = − 0.0335 -0.006693\times5=-0.0335 − 0.006693 × 5 = − 0.0335 の和で W 1 = 2.4665 k B T W_1=2.4665\, k_{\mathrm B} T W 1 = 2.4665 k B T である。h b ( 0.006693 ) = 5 × 0.006693 + ln ( 1 + e − 5 ) = 0.0402 h_{\mathrm b}(0.006693)=5\times0.006693+\ln(1+\mathrm{e}^{-5})=0.0402 h b ( 0.006693 ) = 5 × 0.006693 + ln ( 1 + e − 5 ) = 0.0402 なので Δ F = ( ln 2 − 0.0402 ) k B T = 0.6529 k B T \Delta F=(\ln2-0.0402) k_{\mathrm B} T=0.6529\, k_{\mathrm B} T Δ F = ( ln 2 − 0.0402 ) k B T = 0.6529 k B T 、余分は 1.8136 k B T 1.8136\, k_{\mathrm B} T 1.8136 k B T である。一方 D ( π 0 ∥ π 5 ) = 1 2 ln 0.5 0.99331 + 1 2 ln 0.5 0.006693 = 1.8136 D(\pi_0\,\|\,\pi_5)=\frac12\ln\frac{0.5}{0.99331}+\frac12\ln\frac{0.5}{0.006693}=1.8136 D ( π 0 ∥ π 5 ) = 2 1 ln 0.99331 0.5 + 2 1 ln 0.006693 0.5 = 1.8136 で一致する。W 1 W_1 W 1 は、間隔 x = 5 x=5 x = 5 の熱浴の量子ビット1個との SWAP で熱浴へ流れる熱 β Q B = 5 ( 1 2 − 0.006693 ) = 2.4665 \beta Q_{\mathrm B}=5(\frac12-0.006693)=2.4665 β Q B = 5 ( 2 1 − 0.006693 ) = 2.4665 に等しく、余分も式 (4.8.12) の D D D と同じである。2つの操作は、仕事を系の間隔の上げ下げに入れるか、SWAP のエネルギーの差に入れるかが違うだけである。
(2) N = 2 N=2 N = 2 では x 1 = 2.5 x_1=2.5 x 1 = 2.5 で、p 1 ( 2.5 ) = 0.07586 p_1(2.5)=0.07586 p 1 ( 2.5 ) = 0.07586 。余分は D ( π 0 ∥ π 2.5 ) + D ( π 2.5 ∥ π 5 ) = 0.6357 + 0.1175 = 0.7532 k B T D(\pi_0\,\|\,\pi_{2.5})+D(\pi_{2.5}\,\|\,\pi_5)=0.6357+0.1175=0.7532\, k_{\mathrm B} T D ( π 0 ∥ π 2.5 ) + D ( π 2.5 ∥ π 5 ) = 0.6357 + 0.1175 = 0.7532 k B T である。近似は 1.2333 / 2 = 0.6166 1.2333/2=0.6166 1.2333/2 = 0.6166 で、2割ほど小さい。N = 10 N=10 N = 10 では厳密な和が 0.1284 0.1284 0.1284 、近似が 0.1233 0.1233 0.1233 で、差は4% に縮む。段を細かくするほど各段の D D D が2次の展開で書けるようになり、近似が正確になる。
1 R. Landauer (1961), IBM J. Res. Dev. 5 , 183.
2 D. Reeb, M. M. Wolf (2014), New J. Phys. 16 , 103011.
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com