相対エントロピーから情報と物理を組み立てる 第4部 第7章 §4.7.5 ★★ 回転するコロイド粒子で、なされた仕事より多い自由エネルギーを粒子に与えた。一般化 Jarzynski 等式のもう1つの形 ⟨eβ(ΔF−W)⟩=γ を確かめ、γ は逆向きの手続きで元の軌道が出る確率の和である
らせん階段のポテンシャル
情報を自由エネルギーに変える最初の実験は S. Toyabe, T. Sagawa, M. Ueda, E. Muneyuki and M. Sano (2010), Nature Phys. 6, 988 である(数値は arXiv:1009.5287 の版による)。直径 287nm のポリスチレンの粒子を2つつなげた二量体を、ガラス面に1点で留める。粒子はその点のまわりに回転の Brown 運動をする。4つの電極に 1MHz の電圧をかけて、回転の角度について周期 180∘ の正弦波のポテンシャルと、一定のトルクを同時に加える。ポテンシャルは傾いた周期の山で、1段の高さは kBT の程度、平均すると粒子は傾きにそって下る(らせん階段を下りる)。
フィードバックは次のとおりである。tc=44ms ごとに粒子の角度を測り、粒子が階段の「上がり口」にあたる領域 S にいれば、遅れ td の後にポテンシャルの山と谷を入れかえる(位相を半周期ずらす)。入れかえると粒子は谷の上り側の斜面に置かれ、上の段へ転がる。S にいなければ何もしない。S は、入れかえの前のポテンシャルが入れかえの後より必ず高い領域に選んであるので、入れかえの仕事はいつも負か0である。遅れ td が短いと粒子は階段を登り、長いと下る。周期 tc=44ms は井戸の中の緩和(10ms ほど)より長く、隣の井戸への跳び移り(1s ほど)より短い。
最も短い遅れ td=1.1ms で、粒子が S にいる確率は p=0.059、1周あたりの自由エネルギーの増えと仕事の差の平均は ⟨ΔF−W⟩=0.062kBT だった。測定に誤りが無いので、得た情報は hb(0.059)=0.22 nat で、そのうち 0.062/0.22=28% が自由エネルギーに変わった。第二法則 ⟨ΔF−W⟩≤0 は破れているが、式 (4.7.14) の ⟨ΔF−W⟩≤kBTI は守られている。
一般化 Jarzynski 等式の γ の形
沙川と上田は、式 (4.7.14) と並べてもう1つの形
⟨eβ(ΔF−W)⟩=γ,γ:=m∑P~m(m)(4.7.15)
を示した(第4部第4章)。P~m(m) は、測定結果 m のときに行ったフィードバックの手続きを時間について逆にたどり、最後に測定したとき結果 m が出る確率である。γ は情報をどれだけ使えたかを表す量で、測定結果を全く使わなければ γ=1 で普通の Jarzynski の等式に戻り、K 通りの結果を完全に使えば γ=K まで大きくなる。Jensen の不等式から ⟨ΔF−W⟩≤kBTlnγ が出る。Toyabe らは逆の手続きを別に行って γ を測り、式 (4.7.15) が遅れ td の広い範囲で成り立つこと、普通の Jarzynski の等式(右辺が1)はフィードバックがあると破れることを示した。td を長くすると γ は1に近づいた。
4.7.4 節の理想化した単電子箱では γ を書き下せる。測定結果 m=0 の手続きの逆は、縮退点から ng を準静的に 1/2−Δng まで動かし(平衡のまま)、速く 1/2 に戻して測る、というものである。動かした先の平衡で状態 0 にいる確率は 1/(1+x) で、誤り率 ε の測定で結果 0 が出る確率は 1+x1−ε+1+xεx である。m=1 も対称なので
γ=1+x2(1−ε+εx)(4.7.16)
である。一方、式 (4.7.12) から ⟨e−βW⟩=(1−ε)1+x2+ε1+x2x で、確かに γ に等しい。
例題 4.7.5 2つの上限を比べる ★★
オリジナル(a) Toyabe らの p=0.059 から hb(p) を求め、情報から自由エネルギーへの変換の割合を確かめよ。p が 1/2 なら上限はいくらか。(b) 理想化した単電子箱で ε=0.05、βΔE=4.30 のとき、γ と lnγ、I(X;M)、実際の −β⟨W⟩ を比べよ。2つの上限のうちどちらが厳しいか。
解答
(a) hb(0.059)=0.059×2.830+0.941×0.0608=0.167+0.057=0.224 nat で、0.062/0.224=0.28 である(論文は I=0.22 と丸めて 28%)。p=1/2 なら情報は ln2=0.693 nat が上限で、実験の 0.062 はその9%にあたる。S を狭く選ぶと、S にいる確率 p が小さく情報も小さいが、S にいたときに確実に登らせることができる。
(b) x=e−4.30=0.0136 で、γ=2(0.95+0.05×0.0136)/1.0136=1.876、lnγ=0.629 である。I=ln2−hb(0.05)=0.495 で、実際の値は例題 4.7.4 の 0.465 である。0.465≤0.495≤0.629 で、この場合は相互情報量の上限のほうが厳しい。γ はフィードバックの手続きを決めてから定まる量で、I は測定だけで決まる量である。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com