tecktom.com

この章のまとめ

2026年10月15日10分で読めますAI生成教科書『虎の巻』第4部 第6章 熱力学的不確定性関係と速度限界

AI 生成・木下が未勉強

相対エントロピーから情報と物理を組み立てる 第4部 第6章 §4.6.7 ★★★ 比べるのは経路の分布と流れを強めた経路の分布。その D=θ2⟨Σ⟩/4D=\theta^2\langle \Sigma \rangle /4 に Cramér–Rao を当てると TUR が出て、同じ補題から速度限界と熱機関の交換関係が出る

この章では、定常状態の経路の分布 PP と、分布を変えずに流れだけをそろって (1+θ)(1+\theta) 倍にした経路の分布 PθP_\theta を比べた。過減衰の Langevin 方程式ではその DD が厳密に θ2⟨Σ⟩/4\theta^2\langle \Sigma \rangle /4 で(命題 4.6.1)、経路の Fisher 情報は ⟨Σ⟩/2\langle \Sigma \rangle /2 である。Markov 跳躍過程では Fisher 情報は t∑jij2/bijt\sum j_{ij}^2/b_{ij} で、対数平均の不等式(補題 4.6.3)によって ⟨Σ⟩/2\langle \Sigma \rangle /2 以下に押さえられる。当てた DD の性質は、DD の2次の係数が Fisher 情報になり、Fisher 情報が推定の精度を Cramér–Rao の不等式で押さえる、という第1部第7章の結果である。

出てきたものは3つある。第一に TUR Var(J)/⟨J⟩2≥2/⟨Σ⟩\mathrm{Var}(\mathcal J)/\langle \mathcal J \rangle ^2\ge2/\langle \Sigma \rangle(定理 4.6.2)で、相対誤差 εJ\varepsilon_{\mathcal J} の流れを作るには 2kBT/εJ22 k_{\mathrm B} T/\varepsilon_{\mathcal J}^2 以上の散逸が要る。等号に近づくのは駆動の弱い線形応答の領域で、そこでは Einstein の関係(揺動散逸定理)のもとで等号が成り立つ。第二に速度限界で、分布を p0p_0 から pτp_\tau へ移すのに要るエントロピー生成は、移した距離の2乗に比例し時間に反比例する。距離を全変動で測ると Shiraishi–Funo–Saito の限界(定理 4.6.4)、Wasserstein 距離で測ると最適輸送の限界(定理 4.6.5)、Fisher 情報で測ると情報幾何の限界 (4.6.21) になり、最後のものは量子力学の Mandelstam–Tamm の限界と同じ形をしている。第三に、定常的に働く熱機関の出力は ΔP(ηC−η)/(2kBTcη)\Delta_P(\eta_{\mathrm C}-\eta)/(2 k_{\mathrm B} T_{\mathrm c}\eta) 以下で、Carnot 効率と有限の出力と小さな揺らぎは同時に得ることはできない。

3つとも、「エントロピー生成は、ある変化を経路の上で見分けられるようにするための代価である」という1つの見方から出ている。流れを正確に作ることも、分布を速く移すことも、経路の分布を変えることで、その変化の大きさを DD と Fisher 情報で測ると、エントロピー生成が下から押さえる。これらの関係が実験でどう確かめられてきたかは、次の 第4部第7章 で扱う。

要点

  • 仮想的な駆動 F→F+θγνF\to F+\theta\gamma\nu(ν=js/ps\nu=j_{\mathrm s}/p_{\mathrm s})で流れだけが (1+θ)(1+\theta) 倍。D(P ∥ Pθ)=θ2⟨Σ⟩/4D(P\,\|\,P_\theta)=\theta^2\langle \Sigma \rangle /4、Jpath=⟨Σ⟩/2J_{\mathrm{path}}=\langle \Sigma \rangle /2。
  • TUR:Var(J)/⟨J⟩2≥2/⟨Σ⟩\mathrm{Var}(\mathcal J)/\langle \mathcal J \rangle ^2\ge2/\langle \Sigma \rangle。跳躍過程では ≥1/Jpath≥2/⟨Σ⟩\ge1/J_{\mathrm{path}}\ge2/\langle \Sigma \rangle((x−y)ln⁡(x/y)≥2(x−y)2/(x+y)(x-y)\ln(x/y)\ge2(x-y)^2/(x+y))。
  • 一様な環 R=(Ac/N)coth⁡(Ac/2N)≥2\mathcal R=(\mathcal A_{\mathrm c}/N)\coth(\mathcal A_{\mathrm c}/2N)\ge2。弱い駆動・細かい段で等号に近づく。
  • 速度限界:τ≥∥pτ−p0∥12/(2⟨Σ⟩⟨A⟩τ)\tau\ge\lVert p_\tau-p_0\rVert _1^2/(2\langle \Sigma \rangle \langle \mathcal A \rangle _\tau)、kBT⟨Σ⟩≥γW22/τk_{\mathrm B} T\langle \Sigma \rangle \ge\gamma\mathcal W_2^2/\tau、∫Jt dt≥2arccos⁡∑p0pτ\int\sqrt{J_{\mathrm t}}\,\mathrm{d} t\ge2\arccos\sum\sqrt{p_0p_\tau}。
  • 熱機関:P≤ΔP(ηC−η)/(2kBTcη)P\le\Delta_P(\eta_{\mathrm C}-\eta)/(2 k_{\mathrm B} T_{\mathrm c}\eta)。量子:ΔE τ≥ℏarccos⁡∣⟨ψ(0)∣ψ(τ)⟩∣\Delta E\,\tau\ge\hbar\arccos\lvert \langle \psi(0) | \psi(\tau) \rangle \rvert。

総合問題 3状態の環で回る分子モーター ★★★

オリジナル3つの状態 1→2→3→11\to2\to3\to1 を環状に回るモーターがある。どの状態からも前へ k+=2 s−1k_+=2\,\mathrm s^{-1}、後ろへ k−=1 s−1k_-=1\,\mathrm s^{-1} の速さで移り、1周するごとに周回の数 J\mathcal J が1増える(後ろへ1周すると1減る)。

  1. (1) 定常分布、各辺の流れ jj、周回の数の平均の速さ ⟨J⟩/t\langle \mathcal J \rangle /t、エントロピー生成の速さ Σ˙\dot\Sigma を求めよ。(使う道具:4.6.2 節、第4部第1章・第4部第3章 のエントロピー生成)
  2. (2) 長い時間での Var(J)/t\mathrm{Var}(\mathcal J)/t を求め、TUR の比 R\mathcal R を計算して R≥2\mathcal R\ge2 を確かめよ。(使う道具:4.6.3 節)
  3. (3) 仮想的な駆動の経路の Fisher 情報 Jpath(0)/tJ_{\mathrm{path}}(0)/t を求め、Σ˙/2\dot\Sigma/2 と比べよ。式 (4.6.13) の2つの不等号のどちらが等号かを述べよ。(使う道具:4.6.2 節、第1部第7章 の Cramér–Rao)
  4. (4) 周回の数で時間を測り、相対誤差 1% を得るのに要る時間と ⟨Σ⟩\langle \Sigma \rangle を求め、TUR の下限と比べよ。(使う道具:4.6.3 節)
  5. (5) 時刻0に全確率を状態1に置いて放す。τ=0.2 s\tau=0.2\,\mathrm s での p(τ)p(\tau)、∥pτ−p0∥1\lVert p_\tau-p_0\rVert _1、その間のエントロピー生成 ⟨Σ⟩\langle \Sigma \rangle(数値)を使い、速度限界 (4.6.16) を確かめよ。(使う道具:4.6.4 節、第4部第1章 の緩和)
解答

(1) どの状態も同じ速さで出入りするので、定常分布は一様 pi=13p_i=\frac13 である。各辺の流れは j=13(k+−k−)=13 s−1j=\frac13(k_+-k_-)=\frac13\,\mathrm s^{-1} で、周回の数は辺を1本通るごとに 13\frac13 進むから⟨J⟩/t=j=13 s−1\langle \mathcal J \rangle /t=j=\frac13\,\mathrm s^{-1} である。式 (4.6.11) で各辺の対数は ln⁡(k+/k−)=ln⁡2\ln(k_+/k_-)=\ln2 なので Σ˙=3×13ln⁡2=ln⁡2=0.693 s−1\dot\Sigma=3\times\frac13\ln2=\ln2=0.693\,\mathrm s^{-1} である。1周あたりの親和力は Ac=3ln⁡2\mathcal A_{\mathrm c}=3\ln2 で、1周するごとに 3ln⁡23\ln2 のエントロピーが生成される。

(2) 状態によらずいつも前へ k+k_+、後ろへ k−k_- で跳ぶので、前後の跳躍の回数は独立な Poisson 過程で、正味の跳躍数の分散は (k++k−)t=3t(k_++k_-)t=3t である。周回の数はその 13\frac13 なのでVar(J)/t=3/9=13 s−1\mathrm{Var}(\mathcal J)/t=3/9=\frac13\,\mathrm s^{-1} である。R=Σ˙⋅1/3(1/3)2=3ln⁡2=2.079≥2\mathcal R=\dot\Sigma\cdot\frac{1/3}{(1/3)^2}=3\ln2=2.079\ge2 で、式 (4.6.15) の N=3N=3 の値 ln⁡2coth⁡(ln⁡22)=ln⁡2×3\ln2\coth(\frac{\ln2}{2})=\ln2\times3 と一致する。

(3) 各辺で b=13(k++k−)=1 s−1b=\frac13(k_++k_-)=1\,\mathrm s^{-1} なので、式 (4.6.12) から Jpath(0)/t=3×j2/b=3×19=13 s−1J_{\mathrm{path}}(0)/t=3\times j^2/b=3\times\frac19=\frac13\,\mathrm s^{-1} で、Σ˙/2=0.347 s−1\dot\Sigma/2=0.347\,\mathrm s^{-1} より小さい。1段目の下限は 1/Jpath(0)=3/t1/J_{\mathrm{path}}(0)=3/t で、実際の Var(J)/⟨J⟩2=(t/3)/(t/3)2=3/t\mathrm{Var}(\mathcal J)/\langle \mathcal J \rangle ^2=(t/3)/(t/3)^2=3/t に一致する。したがって Cramér–Rao の不等号は等号で、TUR が開いている理由は対数平均の不等号(各辺の往きと帰りの比が2で、平衡から離れていること)だけである。例題 4.6.2 の格子の上のモーターと同じ構造である。

(4) εJ2=3/t\varepsilon_{\mathcal J}^2=3/t なので、εJ=0.01\varepsilon_{\mathcal J}=0.01 には t=3/10−4=3.00×104 st=3/10^{-4}=3.00\times10^4\,\mathrm s(約8.3時間)が要る。そのあいだの ⟨Σ⟩=ln⁡2×3×104=2.08×104\langle \Sigma \rangle =\ln2\times3\times10^4=2.08\times10^4 で、TUR の下限 2/εJ2=2.00×1042/\varepsilon_{\mathcal J}^2=2.00\times10^4 の 1.04 倍(=R/2=\mathcal R/2)である。散逸に換算すると 2.08×104 kBT2.08\times10^4\, k_{\mathrm B} T で、TUR の下限からの上乗せは4%にとどまる。

(5) 遷移の行列は巡回行列なので、p(t)p(t) は離散 Fourier 変換で対角化することができる。ω=e2πi/3\omega=\mathrm{e}^{2\pi\mathrm{i}/3} として固有値は k+(ωm−1)+k−(ω−m−1)k_+(\omega^m-1)+k_-(\omega^{-m}-1)(m=0,1,2m=0,1,2)で、m=1m=1 は −92+32i-\frac92+\frac{\sqrt3}{2}\mathrm{i} である。数値で解くと p(0.2)=(0.6003, 0.2403, 0.1594)p(0.2)=(0.6003,\,0.2403,\,0.1594)、∥pτ−p0∥1=2(1−0.6003)=0.7994\lVert p_\tau-p_0\rVert _1=2(1-0.6003)=0.7994 である。活動度は分布によらず A=k++k−=3 s−1\mathcal A=k_++k_-=3\,\mathrm s^{-1} で、Σ˙\dot\Sigma を式 (4.6.11)(その時刻の pp で)から時間で積分すると ⟨Σ⟩=1.080\langle \Sigma \rangle =1.080 である(t=0t=0 では p2=p3=0p_2=p_3=0 で被積分関数は対数的に発散するが、積分は有限である)。右辺は 0.79942/(2×1.080×3)=0.0986 s0.7994^2/(2\times1.080\times3)=0.0986\,\mathrm s で、実際の τ=0.2 s\tau=0.2\,\mathrm s を下回り、速度限界が確かに成り立つ。初めの分布は定常分布から遠く、緩和のエントロピー生成と、環を回し続けるためのエントロピー生成の両方が ⟨Σ⟩\langle \Sigma \rangle に入っている。

木下の勉強ステータス

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

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

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