tecktom.com

出てくるもの (1):精度にはエントロピー生成が要る

2026年10月15日8分で読めますAI生成教科書『虎の巻』第4部 第6章 熱力学的不確定性関係と速度限界
目次
  1. 精度の代価
  2. 環状の模型と等号の条件
  3. TUR の成り立つ範囲
  4. 木下の層

AI 生成・木下が未勉強

相対エントロピーから情報と物理を組み立てる 第4部 第6章 §4.6.3 ★★★ 相対誤差 εJ\varepsilon_{\mathcal J} を得るには ⟨Σ⟩≥2/εJ2\langle \Sigma \rangle \ge2/\varepsilon_{\mathcal J}^2。等号に近づくのは駆動が弱い極限で、そのかわり時間がかかる

精度の代価

式 (4.6.9) を、相対誤差 εJ\varepsilon_{\mathcal J} を達成するために要るエントロピー生成の下限として読むと

kBT⟨Σ⟩ ≥ 2kBTεJ2(4.6.14)k_{\mathrm B} T\langle \Sigma \rangle \ \ge\ \frac{2 k_{\mathrm B} T}{\varepsilon_{\mathcal J}^2}\tag{4.6.14}

である。左辺は、時間 tt のあいだに系が熱浴へ捨てたエネルギーのうち、熱浴が平衡のまま吸い込んだ分(散逸)である。右辺には時間 tt も系の速さも模型の細部も入っていない。相対誤差1% の流れを作るには、どんな仕組みを使っても 2×104 kBT2\times10^4\, k_{\mathrm B} T 以上を散逸しなければならない。Barato と Seifert が強調したのはこの読み方で、細胞の中の時計や分子モーターの正確さには、熱力学が決める最低の値段がついている。

tt が入らないのは、定常状態では Var(J)\mathrm{Var}(\mathcal J) も ⟨J⟩\langle \mathcal J \rangle も ⟨Σ⟩\langle \Sigma \rangle も(長い時間では)tt に比例し、Var/⟨J⟩2∝1/t\mathrm{Var}/\langle \mathcal J \rangle ^2\propto1/t と ⟨Σ⟩∝t\langle \Sigma \rangle \propto t の積で tt が消えるからである。長く測れば相対誤差は 1/t1/\sqrt t で下がるが、散逸も tt に比例して増える。その交換の比率が、系によらず2以上に決まっている。

環状の模型と等号の条件

TUR がどこまで厳しいかを、NN 個の状態を環状につないだ一様な模型で見る。各状態から前へ k+k_+、後ろへ k−k_- の速さで移り、1周すると流れが1つ進むとする(N=1N=1 は前の例題の格子の上のモーター)。1周あたりのエントロピー生成(親和力)を Ac≔Nln⁡(k+/k−)\mathcal A_{\mathrm c}\coloneqq N\ln(k_+/k_-) とおく。定常分布は一様で、周回の数 J\mathcal J について長い時間で⟨J⟩=(k+−k−)t/N\langle \mathcal J \rangle =(k_+-k_-)t/N、Var(J)=(k++k−)t/N2\mathrm{Var}(\mathcal J)=(k_++k_-)t/N^2、⟨Σ⟩=(k+−k−)tln⁡(k+/k−)\langle \Sigma \rangle =(k_+-k_-)t\ln(k_+/k_-) だから、TUR の比は

R≔⟨Σ⟩ Var(J)⟨J⟩2=AcNcoth⁡Ac2N(4.6.15)\mathcal R\coloneqq\langle \Sigma \rangle \,\frac{\mathrm{Var}(\mathcal J)}{\langle \mathcal J \rangle ^2}=\frac{\mathcal A_{\mathrm c}}{N}\coth\frac{\mathcal A_{\mathrm c}}{2N}\tag{4.6.15}

になる((k++k−)/(k+−k−)=coth⁡(ln⁡(k+/k−)/2)(k_++k_-)/(k_+-k_-)=\coth(\ln(k_+/k_-)/2) を使った)。R≥2\mathcal R\ge2 が TUR で、xcoth⁡(x/2)≥2x\coth(x/2)\ge2 は補題 4.6.3 と同じ不等式である。

図

図 4.6.1 NN 個の状態を環状につないだ一様な模型の TUR の比 R\mathcal R(式 (4.6.15)、厳密な式から描いた)。黒の実線が N=1N=1、青の破線が N=3N=3、赤の一点鎖線が N=10N=10、灰色の実線が下限 2。親和力 Ac→0\mathcal A_{\mathrm c}\to0 ではどれも 2 に近づき、同じ親和力なら状態を多く分けるほど 2 に近い。

図 4.6.1 から2つのことを読み取ることができる。第一に、Ac→0\mathcal A_{\mathrm c}\to0(平衡に近い弱い駆動)で R→2\mathcal R\to2 になり、TUR は等号に近づく。これは第6節で見る線形応答の領域である。第二に、同じ親和力でも1周を多くの段に分けるほど R\mathcal R は 2 に近い。N→∞N\to\infty では R=2+Ac2/(6N2)+⋯→2\mathcal R=2+\mathcal A_{\mathrm c}^2/(6N^2)+\cdots\to2 である。段の数を増やすことは、1段あたりの駆動を弱くすることなので、これも「各辺が平衡に近いほど等号に近い」という補題 4.6.3 の等号条件の表れである。第4部第1章で、ばねを NN 段に分けて動かすと余分な仕事が 1/N1/N で下がることを見た。ここでも、1回の大きな跳躍を小さな段に分けることが散逸を無駄にしない方法になっている。

ただし、R→2\mathcal R\to2 に近づける方法はどれも時間を犠牲にする。駆動を弱くすると流れが遅くなり、同じ相対誤差に達するまでの時間が延びる(例題 4.6.3)。TUR は散逸の下限を時間によらずに与えるが、その下限に近づくためには遅く動かなければならない。

TUR の成り立つ範囲

定理 4.6.2 の証明で使ったのは、(i) 系が時間によらない駆動のもとの定常状態にあること、(ii) 経路を見分ける問題が過減衰の Langevin 方程式か Markov 跳躍過程で書けること、の2つである。(i) が崩れて、駆動が時間に周期的に依存する場合や、初めの分布が定常分布でない有限の時間の過程では、式 (4.6.9) がそのままの形では成り立たない例がある。その理由は、仮想的な駆動で分布を変えずに流れだけを (1+θ)(1+\theta) 倍にする操作ができなくなることにある。そうした場合には、PθP_\theta の分布の変化の分も Fisher 情報に足した一般化された形が使われる。(ii) が崩れて、速度のように時間反転で符号を変える変数を持つ慣性のある系では、逆回しの経路で速度の符号も反転させるため、エントロピー生成と流れの関係が変わり、やはり元の形の TUR が破れる例がある。どちらの場合も「DD の2次の展開と Cramér–Rao」という骨組みは同じで、比べる相手の PθP_\theta の作り方と、その DD がエントロピー生成でどう押さえられるかが変わる。

例題 4.6.3 分子の時計の精度と値段 ★★

オリジナル前の例題のモーターを時計に使い、正味の歩数 J\mathcal J で時間を測る。後ろへの速さを k−=1 s−1k_-=1\,\mathrm s^{-1} に固定し、1歩あたりの親和力 a=ln⁡(k+/k−)a=\ln(k_+/k_-) を選べるとする。相対誤差 εJ=0.01\varepsilon_{\mathcal J}=0.01 を目標にする。

  1. (1) 目標の精度に達するのに要るエントロピー生成 ⟨Σ⟩\langle \Sigma \rangle を aa の関数として求めよ。
  2. (2) a=0.1a=0.1、22、2020 のそれぞれで ⟨Σ⟩\langle \Sigma \rangle と、要る時間 tt を求めよ。
  3. (3) 結果を TUR の下限と比べ、駆動の強さの選び方について何が言えるかを述べよ。
解答

(1) 例題 4.6.2 から εJ2=(k++k−)/((k+−k−)2t)\varepsilon_{\mathcal J}^2=(k_++k_-)/((k_+-k_-)^2t)、⟨Σ⟩=(k+−k−)ta\langle \Sigma \rangle =(k_+-k_-)ta なので、2つを掛けて tt を消すとεJ2⟨Σ⟩=acoth⁡(a/2)\varepsilon_{\mathcal J}^2\langle \Sigma \rangle =a\coth(a/2)、つまり ⟨Σ⟩=acoth⁡(a/2)/εJ2\langle \Sigma \rangle =a\coth(a/2)/\varepsilon_{\mathcal J}^2 である(式 (4.6.15) の N=1N=1)。

(2) 時間は t=(k++k−)/((k+−k−)2εJ2)=(ea+1)/((ea−1)2εJ2)t=(k_++k_-)/((k_+-k_-)^2\varepsilon_{\mathcal J}^2)=(\mathrm{e}^a+1)/((\mathrm{e}^a-1)^2\varepsilon_{\mathcal J}^2) 秒である。a=0.1a=0.1:acoth⁡(a/2)=2.0008a\coth(a/2)=2.0008 で ⟨Σ⟩=2.0008×104\langle \Sigma \rangle =2.0008\times10^4、t=1.90×106 st=1.90\times10^6\,\mathrm s(約22日)。a=2a=2:2coth⁡1=2.6262\coth1=2.626 で ⟨Σ⟩=2.63×104\langle \Sigma \rangle =2.63\times10^4、t=2.06×103 st=2.06\times10^3\,\mathrm s。a=20a=20:20coth⁡10=20.020\coth10=20.0 で ⟨Σ⟩=2.00×105\langle \Sigma \rangle =2.00\times10^5、t=2.06×10−5 st=2.06\times10^{-5}\,\mathrm s である。

(3) TUR の下限は 2/εJ2=2×1042/\varepsilon_{\mathcal J}^2=2\times10^4 である。弱い駆動(a=0.1a=0.1)は下限のほぼ等号を達成するが、22日かかる。強い駆動(a=20a=20)は一瞬で目標に達するが、下限の10倍の散逸を払う。a=2a=2 は下限の1.3倍の散逸で、時間も現実的である。散逸を下限に近づけることと、速く測ることは両立しない。これが次の節の速度限界の主題でもある。

木下の勉強ステータス

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

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

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