相対エントロピーから情報と物理を組み立てる 第4部 第2章 §4.2.3 ★★★ 前向きの仕事の分布と、逆回しの仕事の分布を裏返したものの比は eβ(W−ΔF)。2つの分布は W=ΔF で交わる
仕事の分布の比
前向きの操作を pλ0eq から始めたときの仕事の分布を PF(W)、逆回しの操作を pλNeq から始めたときの仕事の分布を PR(W) と書く。連続の値をとる場合は確率密度、離散の値の場合は確率とする。
定理 4.2.5(Crooks の定理 ★★★)
局所詳細つり合いのもとで、前向きと逆回しをそれぞれの平衡分布から始めると
PR(−W)PF(W)=eβ(W−ΔF).(4.2.17)
証明
選び方 (b) では、式 (4.2.16) により Σ=β(W−ΔF) で、Σ の値と W の値は1対1に対応する。逆回しの側では (4.2.7) から Σ~=−Σ=β(W~+ΔF) なので、Σ~=−s は W~=−W と同じことである。詳細ゆらぎの定理 (4.2.14) に s=β(W−ΔF) を入れればよい。定理 4.2.3 の仮定(逆回しで起こる経路は前向きでも起こる)は、局所詳細つり合いで遷移確率の正負が両向きでそろうことと、平衡分布がどの状態にも正の確率を与えることから満たされる。
Crooks(1999)、Phys. Rev. E 60、2721 がこの形で示した。式 (4.2.17) から読めることを3つ挙げる。
第一に、2つの分布 PF(W) と PR(−W) は、ちょうど W=ΔF で値が等しくなる。2つの分布を測って交わる点を読めば、操作がどれほど速く非可逆でも、平衡の量 ΔF が求まる。Collin ほか(2005)、Nature 437、231 は、RNA の分子を光ピンセットで引き伸ばす操作と縮める操作の仕事の分布を測り、この交点から折りたたみの自由エネルギーを求めた。
第二に、仕事の値だけから「いま見ているのは前向きか逆回しか」を当てる統計の問題として読める。W の値を見たときの対数尤度比は ln[PF(W)/PR(−W)]=β(W−ΔF) で、散逸仕事そのものである。これを前向きの分布で平均すると
D(PF∥P~R)=β(⟨W⟩F−ΔF),P~R(W):=PR(−W)(4.2.18)
を得る。平均の散逸仕事は、前向きの仕事の分布と、逆回しの仕事の分布を裏返したものの相対エントロピーに kBT を掛けたものである。経路の分布 P から仕事の分布 PF への移りは Γ↦W[Γ] という粗視化なので、データ処理不等式から D(PF∥P~R)≤D(P∥P~) である。選び方 (b) では経路の比が W だけで決まるので等号が成り立つ。比が W だけで決まらない場合(測定やフィードバックがある場合)にこの不等式がどう分かれるかは、第4部第3章 と 第4部第4章 で扱う。
第三に、PF が正規分布 N(μ,σ2) なら、式 (4.2.17) から PR(−W)∝e−(W−μ)2/(2σ2)−βW も正規分布で、平均は μ−βσ2、分散は σ2 になる。規格化の定数を比べると
ΔF=μ−2βσ2,⟨W⟩F−ΔF=2βσ2(4.2.19)
が出る。散逸仕事は仕事のばらつき(分散)だけで決まる。図 4.2.2 にこの場合を示す。

図 4.2.2 Crooks の定理(定理 4.2.5)の正規分布の場合。前向きの仕事の分布 PF(W) を平均 3kBT・分散 2(kBT)2 の正規分布とすると、逆回しの仕事の分布を裏返した PR(−W) は平均 1kBT の正規分布になり、2つは W=ΔF=2kBT で交わる(式 (4.2.19))。
例題 4.2.3 ばねの強さを急に変える ★★
オリジナル調和ポテンシャル U(x,k)=21kx2 の中の粒子を考える。ばね定数 k0 の平衡状態から、瞬間的にばね定数を k1=4k0 に変える(前向き)。逆回しは k1 の平衡状態から瞬間的に k0 に戻す。
- ΔF、前向きと逆回しの平均の仕事 ⟨W⟩F、⟨W⟩R を求めよ。
- PF(W) と PR(−W) を求め、式 (4.2.17) を確かめよ。
- D(PF∥P~R) を求めよ。
解答
(1) 1次元の調和振動子の分配関数は 2πkBT/k に比例するので ΔF=−kBTlnk0/k1=2kBTlnk0k1=kBTln2 である。瞬間的な変化のあいだ粒子は動かないので、仕事は W=21(k1−k0)x2 で、x は初めの平衡分布 ⟨x2⟩=kBT/k0 に従う。⟨W⟩F=2k0k1−k0kBT=1.5kBT。逆回しは ⟨W⟩R=2k1k0−k1kBT=−0.375kBT である。−⟨W⟩R=0.375≤ΔF/kBT=0.693≤⟨W⟩F/kBT=1.5 で、どちらの向きでも平均では第二法則が成り立つ。
(2) c:=21(k1−k0) と置く。x が平均 0・分散 s2 の正規分布のとき、W=cx2(W>0)の確率密度は 2πs2cW1e−W/(2cs2) である。前向きは s2=kBT/k0、逆回しの −WR=cx2 は s2=kBT/k1 なので
PR(−W)PF(W)=k1k0exp[−2cW(kBTk0−kBTk1)]=k1k0eβW=eβ(W−ΔF)
となる。指数の中で k0−k1=−2c を使った。
(3) 式 (4.2.18) から D=β(⟨W⟩F−ΔF)=1.5−ln2=0.807 nat である。W の値を1回見て前向きか逆回しかを当てる問題で、前向きの側にこれだけの手がかりがある。逆向きの組 D(PR∥P~F)=β(⟨W⟩R+ΔF)=0.318 nat は小さく、2つの向きで散逸は対称でない。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com