わかる温度第 4 回 / 第二法則は、「ほぼ」しか成り立たない

小さい系では、エントロピーが減る事象が実際に起きる ── しかも起きる確率が厳密に計算できる

第二法則は、「ほぼ」しか成り立たない 熱力学第二法則は、絶対の掟ではなく圧倒的に有利な賭けです。
破れる確率はきっちり計算できて、\(P(+\sigma)/P(-\sigma)=e^{\sigma/k_B}\)。
コロイド粒子や RNA の実験で確認済み ── 法則に見えるのは、ただ \(N\) が大きいから。

必要な道具:第2回のボルツマン因子、第3回の熱、指数関数、確率の気分 この回の芯:P(+σ)/P(−σ) = e^(σ/k_B)

「エントロピーは増える」── 熱力学第二法則は、物理でいちばん破れそうにない法則として教わります。ところが 1990 年代以降、それが破れる瞬間を実際に観測して数えることができるようになりました。破っているのは、コロイド粒子ひとつや、ほどけかけの RNA 一本といった小さい系です。しかも重要なのは、「たまに例外がある」という曖昧な話ではないこと ── 破れる確率と守る確率の比が、厳密な等式で与えられます。それがゆらぎの定理:\(P(+\sigma)/P(-\sigma)=e^{\sigma/k_B}\)。ここでも \(k_B\) が両替レートとして働いていて、生成エントロピーが \(k_B\) 何個ぶんかという無次元の比だけが効きます。第二法則が「法則」に見えるのは、日常の系ではこの比が \(10^{20}\) といった数になるから、それだけのことです。

01ブラウン運動の頃から、答えは見えていた

水に浮かぶ花粉の微粒子は、まわりの水分子に叩かれてふらふら動きます。このとき微粒子は、時々まわりから熱をもらって速くなり、時々まわりに熱を返して遅くなる。もらっている瞬間だけを切り取れば ── それは熱が冷たい側から熱い側へ流れているということです。

姉妹編「わかる繰り込み」第1回で見たとおり、こうした逆行が見えるのは \(N\) が小さいから。相対的なゆらぎは \(1/\sqrt N\) で、粒子が \(10^{23}\) 個あれば \(10^{-12}\) となって決して見えない。だからブラウン運動の逆行と、第二法則は矛盾しません

問題は、その先です。「たまに逆行する」を、どれくらいの頻度で、と言えるか。

02ゆらぎの定理 ── 逆行の確率が、等式で決まる

この回の芯
$$\frac{P(+\sigma)}{P(-\sigma)}=e^{\sigma/k_B}$$

ある実験で生成エントロピーが \(+\sigma\) になる確率と、\(-\sigma\) になる(=第二法則に逆らう)確率のが、これだけで決まります。
ここでも効くのは無次元の比 \(\sigma/k_B\) ── 「\(k_B\) 何個ぶんのエントロピーか」。第1回で「\(k_B\) は換算レート」と書いた意味が、ここでいちばんはっきりします。

この一行が言っていることは、二つあります。

  1. 逆行は起きる。 \(P(-\sigma)\) はゼロではない。第二法則は「絶対に起きない」ではなく「起きにくい」。
  2. どれくらい起きにくいかが、指数で決まる。 \(\sigma\) が \(k_B\) 数個ぶんなら日常茶飯事、\(k_B\) 100 個ぶんなら \(e^{-100}\)、\(10^{20}\) 個ぶんなら宇宙の年齢でも起きない。
桁で見る ── 第二法則はどれくらい「法則」か
典型的な σ/k_B逆行の相対確率 e^(−σ/k_B)実際に見えるか
光ピンセットで引くコロイド粒子(1 秒)数〜数十10⁻¹ 〜 10⁻⁵見える(実験で日常的に観測)
分子モーター 1 歩≈ 20≈ 2×10⁻⁹ときどき逆走する(観測済み)
コップの水が 1 mK 冷える≈ 10¹⁸10^(−4×10¹⁷)絶対に見えない

第二法則が「法則」なのは、\(\sigma/k_B\) が桁違いに大きい世界にだけ住んでいるからです。ナノスケールに降りると、法則は賭けに戻ります。

03動かしてみる ── 逆行を数える

下の図は、水中の微小な粒子を光ピンセットで一定速度で引きずるという、実験室で実際にやられている設定を、その場で数値計算したものです(過減衰ランジュバン方程式を 3000 本の軌道について解いています)。1 本の軌道ごとに、外から加えた仕事 \(W\) を測ります。

左が \(W\) のヒストグラム。平均は必ず正(第二法則)ですが、赤く塗った左側の裾 ── \(W<0\) の軌道が実在します。これは「粒子が水から熱をもらって、引き手を押し返した」回。第二法則が破れた個々の実験です。

右が検証で、各ビンについて \(\ln\bigl[P(W)/P(-W)\bigr]\) を \(W\) に対して打ったもの。ゆらぎの定理が正しければ、これは傾き 1 の直線に乗るはずです(\(k_BT\) を単位にとっているので)。乗ります。

そしてスライダーで引く時間(=系の「大きさ」)を伸ばすと ── 逆行の割合が指数関数的に消えていきます。第二法則が、目の前で「法則」になっていく

図:光ピンセットで引かれる微粒子。3000 本の軌道を数値計算して仕事 W を測る。左=W のヒストグラム(赤い部分が「第二法則を破った」軌道)、右=ln[P(W)/P(−W)] が傾き 1 の直線に乗ることの検証。単位は k_BT。スライダーで引く時間を伸ばすと逆行が消える
W > 0(第二法則どおり) W < 0(破れた軌道) 理論:傾き 1 の直線

04ヤジンスキー等式 ── 破れる側から、平衡量を取り出す

もうひとつ、同じ家族の等式があります。こちらは実用的です。

ヤジンスキー等式(1997)
$$\bigl\langle e^{-W/k_BT}\bigr\rangle=e^{-\Delta F/k_BT}$$

左辺はどんなに速く(非平衡に)引いてもよい実験の平均。右辺は平衡の量(自由エネルギー差)。乱暴な実験を何度も繰り返して指数平均を取ると、平衡の値がぴったり出るという驚くべき等式です。

しかも、ここから熱力学第二法則が系として導けます ── イェンセンの不等式より \(\langle e^{-W/k_BT}\rangle\ge e^{-\langle W\rangle/k_BT}\) なので、\(\langle W\rangle\ge\Delta F\)。第二法則は等式の系であって、独立な公理ではなかった。

この等式は実験でも使われていて、RNA のヘアピンを光ピンセットで無理やり引きほどく実験から、ほどくのに必要な自由エネルギー(平衡量)が実測されています。ゆらぎに支配された乱暴なデータから、きれいな平衡量が出てくる ── 指数平均が、まれな逆行イベントを拾い上げているからです。

「わかる繰り込み」第1回の 1/√N が、ここで牙をむく 姉妹編第1回では「\(N\) が大きいから相対ゆらぎが \(1/\sqrt N\) で消え、マクロ量は鋭くなる」と書きました。この回はその裏側です ── \(N\) が小さければ、ゆらぎは消えずに法則を破る
そして重要なのは、破れ方がでたらめではないこと。\(P(+\sigma)/P(-\sigma)=e^{\sigma/k_B}\) という、時間反転対称性から出る厳密な構造を持っています。ミクロの法則が時間反転対称であること(姉妹編 番外編①のロシュミット)が、マクロの非対称性の破れ方まで規定している ── これがゆらぎの定理のいちばん深いところです。

05では、第二法則とは何だったのか

言い換え

従来:「孤立系のエントロピーは減少しない」

正確には:「エントロピーが \(\sigma\) だけ減る確率は、同じだけ増える確率の \(e^{-\sigma/k_B}\) 倍である
── そして \(\sigma/k_B\) が巨大な日常の系では、前者が後者の \(10^{-10^{18}}\) 倍になるので、実質的に「減らない」と言ってよい。

第1回から一貫している構図が、ここでも出ています ── 効くのは無次元の比だけ。第1回は \(E/k_BT\)、第3回は \(\hbar\omega/k_BT\)、そしてこの回は \(\sigma/k_B\)。単位のある量そのものではなく、\(k_B\) で割った「何個ぶんか」だけが物理を決めています。

◇ ◇ ◇
正直な線 ── 「破れる」の正確な意味

確立していること:ゆらぎの定理(Evans–Cohen–Morriss 1993、Gallavotti–Cohen 1995、Crooks 1999 ほか)と、その形 \(P(+\sigma)/P(-\sigma)=e^{\sigma/k_B}\);ヤジンスキー等式 \(\langle e^{-W/k_BT}\rangle=e^{-\Delta F/k_BT}\)(1997)とそこから第二法則が不等式として従うこと;コロイド粒子の光ピンセット実験(Wang ら 2002 ほか)や RNA ヘアピンのほどき実験(Liphardt ら 2002、Collin ら 2005 ほか)による実験的検証;過減衰ランジュバン系で仕事分布がガウスになりゆらぎの定理を満たすこと ── いずれも確立した物理です。

誤解しないこと:第二法則は破れていません。破れるのは「個々の試行」であって、平均は必ず \(\langle\sigma\rangle\ge0\)。むしろゆらぎの定理は第二法則を精密化したもので、否定するものではありません。② 永久機関は作れません ── 逆行イベントを選んで取り出すには、どちらが逆行かを測って記憶する必要があり、その情報処理のコスト(ランダウアー、姉妹編「わかる繰り込み」番外編②)が必ず上回ります(マクスウェルの悪魔の現代的決着)。③ 定理が成り立つには、初期状態が平衡である・熱浴が正準分布であるなど、いくつかの条件が要ります。形も設定によって少しずつ違う版(過渡的/定常状態/Crooks/Jarzynski)があり、本文は代表形です。④ 図は過減衰ランジュバン方程式(\(k=\gamma=k_BT=1\)、トラップを一定速度で並進)の数値解で、実験装置の再現ではありません。この設定では \(\Delta F=0\) なので、ゆらぎの定理が \(P(W)/P(-W)=e^{W/k_BT}\) の形になります。

練習問題(この回の考え方で解けます)
  1. \(\sigma=10\,k_B\) のとき、逆行が起きる相対確率は。
    答えを見る
    \(e^{-10}\approx4.5\times10^{-5}\)。2 万回に 1 回くらい。ナノスケールの実験では十分観測できる頻度。分子モーターが「ときどき逆走する」のはこの領域。
  2. ゆらぎの定理は第二法則を否定するのか。
    答えを見る
    否定しない、精密化する。破れるのは個々の試行だけで、平均は必ず \(\langle\sigma\rangle\ge0\)。むしろヤジンスキー等式+イェンセンの不等式から第二法則が導けるので、第二法則は独立な公理ではなく系だった、ということになる。
  3. 逆行イベントだけを選んで取り出せば、永久機関になるか。
    答えを見る
    ならない。どれが逆行かを測って記憶する必要があり、その記憶を消すコスト \(k_BT\ln2\)(ランダウアー)が取り出せるエネルギーを必ず上回る。マクスウェルの悪魔の現代的決着そのもの。
  4. 図で引く時間 τ を伸ばすと逆行が消えるのはなぜか。
    答えを見る
    平均の生成エントロピー \(\langle\sigma\rangle\) が τ に比例して増えるから。ゆらぎの幅は \(\sqrt{\langle\sigma\rangle}\) でしか増えないので、比 \(\langle\sigma\rangle/\sqrt{\langle\sigma\rangle}=\sqrt{\langle\sigma\rangle}\) が大きくなり、負側の裾が指数的に消える。姉妹編「わかる繰り込み」第1回の \(1/\sqrt N\) と同じ勘定

第4回まとめ法則ではなく、圧倒的に有利な賭けだった

小さい系・短い時間では、エントロピーが減る事象が実際に起きます。しかもその頻度は厳密な等式で決まる ── \(P(+\sigma)/P(-\sigma)=e^{\sigma/k_B}\)(ゆらぎの定理)。効くのは無次元の比 \(\sigma/k_B\)=「\(k_B\) 何個ぶんのエントロピーか」だけで、ここでも \(k_B\) は両替レートとして働いています。

コロイド粒子や RNA の実験で確認済み。さらにヤジンスキー等式 \(\langle e^{-W/k_BT}\rangle=e^{-\Delta F/k_BT}\) を使うと、乱暴な非平衡実験から平衡の自由エネルギーが取り出せ、しかもそこから第二法則が不等式として導けます ── 第二法則は独立な公理ではなく、等式の系だった

ただし破れるのは個々の試行だけで、平均は必ず守られます。そして逆行だけを拾おうとすれば、測って記憶するコスト(ランダウアー)が必ず上回る ── 永久機関にはなりません。第二法則が「法則」に見えるのは、私たちが \(\sigma/k_B\sim10^{18}\) の世界に住んでいるから、それだけです。

この文書は「わかる温度」シリーズ第4回、物理好きの高校生・大学生向け読み物です。ゆらぎの定理(Evans–Cohen–Morriss 1993、Gallavotti–Cohen 1995、Crooks 1999 ほか)とその形 \(P(+\sigma)/P(-\sigma)=e^{\sigma/k_B}\)、ヤジンスキー等式(1997)とそこからイェンセンの不等式を通じて熱力学第二法則が従うこと、コロイド粒子の光ピンセット実験および RNA ヘアピンのほどき実験による実験的検証、ならびに過減衰ランジュバン系での仕事分布がガウスとなりゆらぎの定理を満たすことは、いずれも確立した物理です。第二法則が破れるのは個々の試行についてであって平均については破れないこと、逆行イベントの選別には測定と記憶のコスト(ランダウアー原理)が伴い永久機関にならないこと、定理には初期平衡性など条件があり設定ごとに複数の版が存在すること、および図が過減衰ランジュバン方程式(\(k=\gamma=k_BT=1\)、トラップの等速並進、したがって \(\Delta F=0\))の数値解であって実験装置の再現ではないことは、本文「正直な線」に明記しました。 ── 印刷する場合はブラウザの「印刷」から「PDF に保存」を(印刷版ではスライダーと解答は静止・非表示になります)。前後の回:第3回 熱とは、追跡をやめた運動第5回 負の温度は、無限大より熱い目次/姉妹編 わかる繰り込み

印刷 / PDF 化:⌘+P(Windows は Ctrl+P)。画面ではスライダーで引く時間を変えると、第二法則を破った軌道の割合が変わります。右のグラフが傾き 1 の直線に乗るのがゆらぎの定理の検証です。「答えを見る」で解答が開きます。