相対エントロピーから情報と物理を組み立てる 第4部 第4章 §4.4.3 ★★★ 測定を含む経路の確率と、結果ごとに操作を逆回しにした経路の確率の比の対数は β(W−ΔF)+I(x;y)
前の2節は、準静的な操作と平均の仕事に限っていた。この節では、有限の速さの操作も含め、1本1本の経路について同じことを言う。道具は前の章で示した経路の確率の比である。
測定の無い過程の復習
系の状態の時間列(経路)を Γ=(x0,x1,…,xN) と書く。Hamiltonian Eλ(x) のパラメータ λ を、決めておいた手順 Λ=(λ0,λ1,…,λN) で動かし、各段で系は温度 T の熱浴に接して遷移する。1段の遷移は、そのときの λ の Gibbs 分布に対して詳細つり合いを満たすとする(局所詳細つり合い)。系が受け取る仕事は、パラメータを変えた瞬間のエネルギーの変化の和 W[Γ]=∑n=0N−1[Eλn+1(xn)−Eλn(xn)] である。
初めの分布を λ0 の Gibbs 分布とした前向きの経路の確率を P[Γ∣Λ]、初めの分布を終わりのパラメータ λN の Gibbs 分布とし、手順を逆に並べた Λ†=(λN,…,λ0) で動かしたときに、逆回しの経路 Γ~=(xN,…,x0) をたどる確率を P~[Γ~∣Λ†] と書く。前の章で示した Crooks の関係は
P~[Γ~∣Λ†]P[Γ∣Λ]=eβ(W[Γ]−ΔF),ΔF=FλN−Fλ0(4.4.10)
である(第4部第2章)。両辺を P で平均すると、左辺の対数の平均は経路の相対エントロピー D(P∥P~) で、
D(P∥P~)=β(⟨W⟩−ΔF)≥0(4.4.11)
が出る。これが前の章(第4部第3章)の「第二法則は経路の D の非負性」の、初めに平衡にある場合の形である。
測定とフィードバックを入れる
時刻 0 で系の状態 x0 を測り、結果 y を確率 p(y∣x0) で得る。その後のパラメータの手順を結果に合わせて Λy と選ぶ。手順の終わりのパラメータも y によってよく、自由エネルギーの差を ΔFy と書く。初めの分布 peq(x0) は λ0 の Gibbs 分布で、すべての y に共通である。経路と結果の同時確率は
P[Γ,y]=p(y∣x0)P[Γ∣Λy](4.4.12)
である。P[Γ∣Λy] には初めの分布 peq(x0) が含まれている。
比べる相手として、次の「逆回しの過程」を考える。まず結果 y を、前向きの過程と同じ確率 p(y) でくじ引きで決める。次に、終わりのパラメータの Gibbs 分布から出発して、手順 Λy を逆に並べた Λy† で系を動かす。測定はしない。この過程で (Γ~,y) が出る確率は
P~[Γ~,y]=p(y)P~[Γ~∣Λy†](4.4.13)
で、∑y∑Γ~P~[Γ~,y]=∑yp(y)=1 と規格化されている。
定理 4.4.3(フィードバックつきの経路の比 ★★★)
式 (4.4.12)・(4.4.13) の2つの経路の確率について
lnP~[Γ~,y]P[Γ,y]=β(W[Γ]−ΔFy)+I(x0;y),I(x0;y)=lnp(y)p(y∣x0).(4.4.14)
証明
式 (4.4.12) を式 (4.4.13) で割ると
P~[Γ~,y]P[Γ,y]=p(y)p(y∣x0)⋅P~[Γ~∣Λy†]P[Γ∣Λy]
である。y を固めると Λy は決めておいた手順なので、第2の因子に Crooks の関係 (4.4.10) を当てることができ、eβ(W−ΔFy) になる。対数をとれば式 (4.4.14) である。
式 (4.4.14) の右辺の第2項 I(x0;y) が、測定の無い式 (4.4.10) からの違いのすべてである。この項は、測定の結果 y が系の状態 x0 についてどれだけのことを教えたかを、経路1本ごとに測ったものである。y が x0 について何も教えないなら(p(y∣x0)=p(y))、I(x0;y)=0 で式 (4.4.10) に戻る。
測定の時刻が途中にある場合
測定を途中の時刻 tm に行い、それまでの手順はすべての y に共通、それ以降を y によって変える場合も、同じ式が成り立つ。前向きの確率は P[Γ,y]=p(y∣xtm)P[Γ∣Λy](Λy は tm の前で共通の手順と、後で y によって変わる手順をつないだもの)と書けるので、証明の x0 を xtm に替えればよい。式 (4.4.14) の I は、測った時刻の状態と結果の I(xtm;y) になる。測定を何回も行う場合は、結果の列とそれぞれの時刻の状態のあいだの I の和に置き換わる。
例題 4.4.3 測定せずに仕切りを動かす ★★
オリジナルSzilard エンジンで、測定をせずに、いつも仕切りを左の壁から測って体積の割合 0.9 の位置まで準静的に動かしてから取り去る。
- 分子が左にいた回と右にいた回の仕事 W を求め、平均 β⟨W⟩ を計算せよ。
- 逆回しの過程(仕切りを割合 0.9 の位置に入れ、真ん中まで準静的に動かす)で分子が左右にいる確率を求め、経路の比の対数が βW に等しいことを確かめよ。
- 式 (4.4.11) の D(P∥P~) が D((1/2,1/2)∥(0.9,0.1)) に等しいことを示せ。
解答
(1) 左にいた回は体積が V/2 から 0.9V に広がるので βW=−ln1.8、右にいた回は V/2 から 0.1V に縮むので βW=−ln0.2=ln5 である。平均は
β⟨W⟩=−21ln1.8+21ln5=21ln95+21ln5=−0.2939+0.8047=0.5108
で、正である。1回の操作で系は元に戻るので ΔF=0 で、平均で仕事を払う。測定が無いので第二法則 ⟨W⟩≥ΔF が成り立つ。
(2) 逆回しの過程では、全体の一様分布から仕切りを割合 0.9 の位置に入れるので、分子は確率 0.9 で左、0.1 で右にいる。その後は準静的に動かすだけで、分子はその側にとどまる。前向きの過程で分子が左右にいる確率は 1/2 ずつなので、経路の比の対数は、左なら ln(0.5/0.9)=−ln1.8、右なら ln(0.5/0.1)=ln5 で、どちらも (1) の βW に等しい。
(3) 準静的な部分では、経路は分子のいる側 x だけで決まり、細かい動きの確率は前向きと逆向きで等しく打ち消し合う。したがってD(P∥P~)=∑x21lnvx1/2=D((1/2,1/2)∥(0.9,0.1))=0.5108 で、(1) の β⟨W⟩ に一致する。式 (4.4.9) で I=0、qy=(0.9,0.1) とおいた場合でもある。
近似の裏付け
型 B(極限) | 小さいもの 熱浴の相関時間と系の遷移の時間の比 →0(弱結合の van Hove 極限) | 証明 条件付き(熱浴の相関が速く減衰し、系は有限次元。量子系で Davies 1974)。局所詳細つり合いはこの上の仮定 | 行き先 近似の数理 第14章
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com