tecktom.com

出てくるもの (1):Crooks の定理

2026年10月15日8分で読めますAI生成教科書『虎の巻』第4部 第2章 ゆらぎの定理
目次
  1. 仕事の分布の比
  2. 木下の層

AI 生成・木下が未勉強

相対エントロピーから情報と物理を組み立てる 第4部 第2章 §4.2.3 ★★★ 前向きの仕事の分布と、逆回しの仕事の分布を裏返したものの比は eβ(W−ΔF)\mathrm{e}^{\beta(W-\Delta F)}。2つの分布は W=ΔFW=\Delta F で交わる

仕事の分布の比

前向きの操作を pλ0eqp^{\mathrm{eq}}_{\lambda_0} から始めたときの仕事の分布を PF(W)P_{\mathrm F}(W)、逆回しの操作を pλNeqp^{\mathrm{eq}}_{\lambda_N} から始めたときの仕事の分布を PR(W)P_{\mathrm R}(W) と書く。連続の値をとる場合は確率密度、離散の値の場合は確率とする。

定理 4.2.5(Crooks の定理 ★★★)

局所詳細つり合いのもとで、前向きと逆回しをそれぞれの平衡分布から始めると

PF(W)PR(−W)=eβ(W−ΔF).(4.2.17)\frac{P_{\mathrm F}(W)}{P_{\mathrm R}(-W)}=\mathrm{e}^{\beta(W-\Delta F)}.\tag{4.2.17}

証明

選び方 (b) では、式 (4.2.16) により Σ=β(W−ΔF)\Sigma=\beta(W-\Delta F) で、Σ\Sigma の値と WW の値は1対1に対応する。逆回しの側では (4.2.7) から Σ~=−Σ=β(W~+ΔF)\tilde\Sigma=-\Sigma=\beta(\tilde W+\Delta F) なので、Σ~=−s\tilde\Sigma=-s は W~=−W\tilde W=-W と同じことである。詳細ゆらぎの定理 (4.2.14) に s=β(W−ΔF)s=\beta(W-\Delta F) を入れればよい。定理 4.2.3 の仮定(逆回しで起こる経路は前向きでも起こる)は、局所詳細つり合いで遷移確率の正負が両向きでそろうことと、平衡分布がどの状態にも正の確率を与えることから満たされる。

Crooks(1999)、Phys. Rev. E 60、2721 がこの形で示した。式 (4.2.17) から読めることを3つ挙げる。

第一に、2つの分布 PF(W)P_{\mathrm F}(W) と PR(−W)P_{\mathrm R}(-W) は、ちょうど W=ΔFW=\Delta F で値が等しくなる。2つの分布を測って交わる点を読めば、操作がどれほど速く非可逆でも、平衡の量 ΔF\Delta F が求まる。Collin ほか(2005)、Nature 437、231 は、RNA の分子を光ピンセットで引き伸ばす操作と縮める操作の仕事の分布を測り、この交点から折りたたみの自由エネルギーを求めた。

第二に、仕事の値だけから「いま見ているのは前向きか逆回しか」を当てる統計の問題として読める。WW の値を見たときの対数尤度比は ln⁡[PF(W)/PR(−W)]=β(W−ΔF)\ln[P_{\mathrm F}(W)/P_{\mathrm R}(-W)]=\beta(W-\Delta F) で、散逸仕事そのものである。これを前向きの分布で平均すると

D(PF ∥ P~R)=β(⟨W⟩F−ΔF),P~R(W)≔PR(−W)(4.2.18)D\bigl(P_{\mathrm F}\,\|\,\tilde P_{\mathrm R}\bigr)=\beta\bigl(\langle W \rangle _{\mathrm F}-\Delta F\bigr),\qquad \tilde P_{\mathrm R}(W)\coloneqq P_{\mathrm R}(-W)\tag{4.2.18}

を得る。平均の散逸仕事は、前向きの仕事の分布と、逆回しの仕事の分布を裏返したものの相対エントロピーに kBTk_{\mathrm B} T を掛けたものである。経路の分布 PP から仕事の分布 PFP_{\mathrm F} への移りは Γ↦W[Γ]\Gamma\mapsto W[\Gamma] という粗視化なので、データ処理不等式から D(PF ∥ P~R)≤D(P ∥ P~)D(P_{\mathrm F}\,\|\,\tilde P_{\mathrm R})\le D(P\,\|\,\tilde P) である。選び方 (b) では経路の比が WW だけで決まるので等号が成り立つ。比が WW だけで決まらない場合(測定やフィードバックがある場合)にこの不等式がどう分かれるかは、第4部第3章 と 第4部第4章 で扱う。

第三に、PFP_{\mathrm F} が正規分布 N(μ,σ2)\mathcal N(\mu,\sigma^2) なら、式 (4.2.17) から PR(−W)∝e−(W−μ)2/(2σ2)−βWP_{\mathrm R}(-W)\propto\mathrm{e}^{-(W-\mu)^2/(2\sigma^2)-\beta W} も正規分布で、平均は μ−βσ2\mu-\beta\sigma^2、分散は σ2\sigma^2 になる。規格化の定数を比べると

ΔF=μ−βσ22,⟨W⟩F−ΔF=βσ22(4.2.19)\Delta F=\mu-\frac{\beta\sigma^2}{2},\qquad \langle W \rangle _{\mathrm F}-\Delta F=\frac{\beta\sigma^2}{2}\tag{4.2.19}

が出る。散逸仕事は仕事のばらつき(分散)だけで決まる。図 4.2.2 にこの場合を示す。

図

図 4.2.2 Crooks の定理(定理 4.2.5)の正規分布の場合。前向きの仕事の分布 PF(W)P_{\mathrm F}(W) を平均 3kBT3 k_{\mathrm B} T・分散 2(kBT)22( k_{\mathrm B} T)^2 の正規分布とすると、逆回しの仕事の分布を裏返した PR(−W)P_{\mathrm R}(-W) は平均 1kBT1 k_{\mathrm B} T の正規分布になり、2つは W=ΔF=2kBTW=\Delta F=2 k_{\mathrm B} T で交わる(式 (4.2.19))。

例題 4.2.3 ばねの強さを急に変える ★★

オリジナル調和ポテンシャル U(x,k)=12kx2U(x,k)=\frac12kx^2 の中の粒子を考える。ばね定数 k0k_0 の平衡状態から、瞬間的にばね定数を k1=4k0k_1=4k_0 に変える(前向き)。逆回しは k1k_1 の平衡状態から瞬間的に k0k_0 に戻す。

  1. ΔF\Delta F、前向きと逆回しの平均の仕事 ⟨W⟩F\langle W \rangle _{\mathrm F}、⟨W⟩R\langle W \rangle _{\mathrm R} を求めよ。
  2. PF(W)P_{\mathrm F}(W) と PR(−W)P_{\mathrm R}(-W) を求め、式 (4.2.17) を確かめよ。
  3. D(PF ∥ P~R)D(P_{\mathrm F}\,\|\,\tilde P_{\mathrm R}) を求めよ。
解答

(1) 1次元の調和振動子の分配関数は 2πkBT/k\sqrt{2\pi k_{\mathrm B} T/k} に比例するので ΔF=−kBTln⁡k0/k1=kBT2ln⁡k1k0=kBTln⁡2\Delta F=- k_{\mathrm B} T\ln\sqrt{k_0/k_1}=\frac{ k_{\mathrm B} T}{2}\ln\frac{k_1}{k_0}= k_{\mathrm B} T\ln2 である。瞬間的な変化のあいだ粒子は動かないので、仕事は W=12(k1−k0)x2W=\frac12(k_1-k_0)x^2 で、xx は初めの平衡分布 ⟨x2⟩=kBT/k0\langle x^2 \rangle = k_{\mathrm B} T/k_0 に従う。⟨W⟩F=k1−k02k0kBT=1.5 kBT\langle W \rangle _{\mathrm F}=\frac{k_1-k_0}{2k_0} k_{\mathrm B} T=1.5\, k_{\mathrm B} T。逆回しは ⟨W⟩R=k0−k12k1kBT=−0.375 kBT\langle W \rangle _{\mathrm R}=\frac{k_0-k_1}{2k_1} k_{\mathrm B} T=-0.375\, k_{\mathrm B} T である。−⟨W⟩R=0.375≤ΔF/kBT=0.693≤⟨W⟩F/kBT=1.5-\langle W \rangle _{\mathrm R}=0.375\le\Delta F/ k_{\mathrm B} T=0.693\le\langle W \rangle _{\mathrm F}/ k_{\mathrm B} T=1.5 で、どちらの向きでも平均では第二法則が成り立つ。

(2) c≔12(k1−k0)c\coloneqq\frac12(k_1-k_0) と置く。xx が平均 00・分散 s2s^2 の正規分布のとき、W=cx2W=cx^2(W>0W>0)の確率密度は 12πs2cWe−W/(2cs2)\frac{1}{\sqrt{2\pi s^2cW}}\mathrm{e}^{-W/(2cs^2)} である。前向きは s2=kBT/k0s^2= k_{\mathrm B} T/k_0、逆回しの −WR=cx2-W_{\mathrm R}=cx^2 は s2=kBT/k1s^2= k_{\mathrm B} T/k_1 なので

PF(W)PR(−W)=k0k1exp⁡[−W2c(k0kBT−k1kBT)]=k0k1 eβW=eβ(W−ΔF)\frac{P_{\mathrm F}(W)}{P_{\mathrm R}(-W)}=\sqrt{\frac{k_0}{k_1}}\exp\Bigl[-\frac{W}{2c}\Bigl(\frac{k_0}{ k_{\mathrm B} T}-\frac{k_1}{ k_{\mathrm B} T}\Bigr)\Bigr]=\sqrt{\frac{k_0}{k_1}}\,\mathrm{e}^{\beta W}=\mathrm{e}^{\beta(W-\Delta F)}

となる。指数の中で k0−k1=−2ck_0-k_1=-2c を使った。

(3) 式 (4.2.18) から D=β(⟨W⟩F−ΔF)=1.5−ln⁡2=0.807D=\beta(\langle W \rangle _{\mathrm F}-\Delta F)=1.5-\ln2=0.807 nat である。WW の値を1回見て前向きか逆回しかを当てる問題で、前向きの側にこれだけの手がかりがある。逆向きの組 D(PR ∥ P~F)=β(⟨W⟩R+ΔF)=0.318D(P_{\mathrm R}\,\|\,\tilde P_{\mathrm F})=\beta(\langle W \rangle _{\mathrm R}+\Delta F)=0.318 nat は小さく、2つの向きで散逸は対称でない。

木下の勉強ステータス

木下はまだこの記事を勉強していません。

AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com

← AI生成教科書『虎の巻』の記事一覧へ