tecktom.com

出てくるもの (3):熱機関の出力・効率・揺らぎ

2026年10月15日7分で読めますAI生成教科書『虎の巻』第4部 第6章 熱力学的不確定性関係と速度限界
目次
  1. 出力の揺らぎと効率
  2. 木下の層

AI 生成・木下が未勉強

相対エントロピーから情報と物理を組み立てる 第4部 第6章 §4.6.5 ★★ 定常的に働く熱機関では P≤ΔP(ηC−η)/(2kBTcη)P\le\Delta_P(\eta_{\mathrm C}-\eta)/(2 k_{\mathrm B} T_{\mathrm c}\eta)。Carnot 効率に近づくには、出力を0にするか出力の揺らぎを無限にするしかない

出力の揺らぎと効率

温度 ThT_{\mathrm h} の高温の熱浴と TcT_{\mathrm c} の低温の熱浴の間で、時間によらない駆動のもとで働き続ける小さな熱機関(熱電素子や、温度差で回る分子モーター)を考える。系が高温の熱浴から受け取る熱の速さを Q˙h\dot Q_{\mathrm h}、低温の熱浴から受け取る熱の速さを Q˙c\dot Q_{\mathrm c}(熱機関なら負)、外へ出す仕事の速さ(出力)を PP とする。定常状態のエネルギーの保存は P=Q˙h+Q˙cP=\dot Q_{\mathrm h}+\dot Q_{\mathrm c}、効率は η≔P/Q˙h\eta\coloneqq P/\dot Q_{\mathrm h}、Carnot 効率は ηC≔1−Tc/Th\eta_{\mathrm C}\coloneqq1-T_{\mathrm c}/T_{\mathrm h} である。系の分布は変わらないので、エントロピー生成は2つの熱浴のエントロピーの増えの和で

Σ˙=−Q˙hkBTh−Q˙ckBTc=−Q˙hkBTh+Q˙h−PkBTc=PkBTc⋅ηC−ηη(4.6.22)\dot\Sigma=-\frac{\dot Q_{\mathrm h}}{ k_{\mathrm B} T_{\mathrm h}}-\frac{\dot Q_{\mathrm c}}{ k_{\mathrm B} T_{\mathrm c}}=-\frac{\dot Q_{\mathrm h}}{ k_{\mathrm B} T_{\mathrm h}}+\frac{\dot Q_{\mathrm h}-P}{ k_{\mathrm B} T_{\mathrm c}}=\frac{P}{ k_{\mathrm B} T_{\mathrm c}}\cdot\frac{\eta_{\mathrm C}-\eta}{\eta}\tag{4.6.22}

になる(Q˙h=P/η\dot Q_{\mathrm h}=P/\eta と 1/Tc−1/Th=ηC/Tc1/T_{\mathrm c}-1/T_{\mathrm h}=\eta_{\mathrm C}/T_{\mathrm c} を使った)。Σ˙≥0\dot\Sigma\ge0 から η≤ηC\eta\le\eta_{\mathrm C} である。

時間 tt のあいだに出した仕事 WoutW_{\mathrm{out}} は、逆回しの経路で符号が変わる流れである。その平均は PtPt、分散は長い時間で ΔPt\Delta_Pt と書ける(ΔP\Delta_P は出力の揺らぎの強さ)。TUR (4.6.9) に式 (4.6.22) を入れると ΔPtP2t2≥2kBTcηPt(ηC−η)\frac{\Delta_Pt}{P^2t^2}\ge\frac{2 k_{\mathrm B} T_{\mathrm c}\eta}{Pt(\eta_{\mathrm C}-\eta)}、整理して

P ≤ ΔP ηC−η2kBTc η(4.6.23)P\ \le\ \Delta_P\,\frac{\eta_{\mathrm C}-\eta}{2 k_{\mathrm B} T_{\mathrm c}\,\eta}\tag{4.6.23}

を得る(P. Pietzonka and U. Seifert (2018), Phys. Rev. Lett. 120, 190602)。

式 (4.6.23) は、大きな出力・Carnot に近い効率・小さな揺らぎの3つを同時には満たせないことを表す。η→ηC\eta\to\eta_{\mathrm C} で右辺は0に近づくので、揺らぎ ΔP\Delta_P が有限のままなら出力は0にならなければならない。逆に有限の出力で Carnot 効率に近づくには、出力の揺らぎが発散しなければならない。マクロな熱力学の教科書では、Carnot 効率は準静的な(したがって出力0の)過程で達成されると習う。式 (4.6.23) は、その「出力0」の理由が揺らぎと散逸の帳尻にあることを、不等式として示している。

図

図 4.6.3 例題 4.6.5 の1段で回る熱機関(Th=400 KT_{\mathrm h}=400\,\mathrm K、Tc=300 KT_{\mathrm c}=300\,\mathrm K、ηC=0.25\eta_{\mathrm C}=0.25、1歩で高温から 10 kBTc10\, k_{\mathrm B} T_{\mathrm c} を受け取る)。黒の点線が上限 (4.6.23) の右辺、赤の実線が模型の値(どちらも厳密な式)。Carnot 効率に近づくほど両者は近づき、ともに0へ向かう。

例題 4.6.5 1段で回る熱機関 ★★

オリジナル1歩前へ進むたびに、高温の熱浴から熱 qhq_{\mathrm h} を受け取り、仕事 ww を外へ出し、残りの qh−wq_{\mathrm h}-w を低温の熱浴へ捨てる機械がある。後ろへの1歩はこの逆である。前後の速さを k±k_\pm とし、局所詳細つり合い k+/k−=eak_+/k_-=\mathrm{e}^{a}(aa は1歩あたりのエントロピー生成)が成り立つとする。

  1. (1) aa を qhq_{\mathrm h}、ww、ThT_{\mathrm h}、TcT_{\mathrm c} で表し、効率 η\eta を求めよ。
  2. (2) PP と ΔP\Delta_P を求め、式 (4.6.23) が tanh⁡(a/2)≤a/2\tanh(a/2)\le a/2 と同値であることを示せ。
  3. (3) Th=400 KT_{\mathrm h}=400\,\mathrm K、Tc=300 KT_{\mathrm c}=300\,\mathrm K、qh=10 kBTcq_{\mathrm h}=10\, k_{\mathrm B} T_{\mathrm c}、w=2 kBTcw=2\, k_{\mathrm B} T_{\mathrm c}、k−=1 s−1k_-=1\,\mathrm s^{-1} で、PP と上限 (4.6.23) の右辺を比べよ。
解答

(1) 前への1歩で高温の熱浴のエントロピーは −qh/Th-q_{\mathrm h}/T_{\mathrm h}、低温の熱浴のエントロピーは +(qh−w)/Tc+(q_{\mathrm h}-w)/T_{\mathrm c} 変わるのでa=qh−wkBTc−qhkBTha=\frac{q_{\mathrm h}-w}{ k_{\mathrm B} T_{\mathrm c}}-\frac{q_{\mathrm h}}{ k_{\mathrm B} T_{\mathrm h}} である。1歩ごとに受け取る熱と出す仕事の比が決まっているので、効率は歩数によらず η=w/qh\eta=w/q_{\mathrm h} である。w=ηqhw=\eta q_{\mathrm h} を入れると a=qhkBTc(ηC−η)a=\frac{q_{\mathrm h}}{ k_{\mathrm B} T_{\mathrm c}}(\eta_{\mathrm C}-\eta) で、a≥0a\ge0 と η≤ηC\eta\le\eta_{\mathrm C} が同値になる。

(2) 正味の歩数は例題 4.6.2 と同じ統計に従うので P=w(k+−k−)P=w(k_+-k_-)、ΔP=w2(k++k−)\Delta_P=w^2(k_++k_-) である。(1) から ηC−ηkBTcη=aw\frac{\eta_{\mathrm C}-\eta}{ k_{\mathrm B} T_{\mathrm c}\eta}=\frac{a}{w} なので、式 (4.6.23) の右辺は w2(k++k−)a2w=wa2(k++k−)w^2(k_++k_-)\frac{a}{2w}=\frac{wa}{2}(k_++k_-) である。不等式は k+−k−≤a2(k++k−)k_+-k_-\le\frac a2(k_++k_-)、すなわち tanh⁡(a/2)≤a/2\tanh(a/2)\le a/2 と同値で、これは a≥0a\ge0 で常に成り立つ。

(3) η=0.2\eta=0.2、ηC=0.25\eta_{\mathrm C}=0.25、a=8−10×300400=0.5a=8-10\times\frac{300}{400}=0.5 で、k+=e0.5=1.6487 s−1k_+=\mathrm{e}^{0.5}=1.6487\,\mathrm s^{-1} である。P=2×0.6487=1.297 kBTc s−1P=2\times0.6487=1.297\, k_{\mathrm B} T_{\mathrm c}\,\mathrm s^{-1}、ΔP=4×2.6487=10.59 (kBTc)2 s−1\Delta_P=4\times2.6487=10.59\,( k_{\mathrm B} T_{\mathrm c})^2\,\mathrm s^{-1}、右辺は 10.59×0.05/(2×0.2)=1.324 kBTc s−110.59\times0.05/(2\times0.2)=1.324\, k_{\mathrm B} T_{\mathrm c}\,\mathrm s^{-1} である。出力は上限の 98% に達している。aa が小さい(Carnot に近い)ので、tanh⁡(a/2)\tanh(a/2) と a/2a/2 の差が小さいからである(図 4.6.3)。

木下の勉強ステータス

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

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

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