ホーム/レッスン/レッスン 4

熱力学第一法則

仕事、熱、内部エネルギー、そして第一法則。

内部エネルギー熱力学第一法則圧力による仕事熱移動状態関数完全微分準静的過程断熱過程示量性

前の講義では熱力学系、平衡状態、状態関数を導入した。ここでは中心的な問いに答えよう。系のエネルギーをこの巨視的な記述にどう組み込むか?

いくつかの講義では、この問いがあまりに手早く片づけられている。第一法則は単なるエネルギー保存として提示され、そのまますぐにエネルギー収支の式が、それ以上の正当化なしに与えられる。これらの式は第 3 節と第 4.3 節で示すので、急ぐ読者は実際そこへ直接進んでもよい(第 7 節の短いまとめも参照)。

しかし、これらの式に至る構成は実のところ微妙である。まず第一法則は、微視的な観点からは(今日では)自明に見えるエネルギー保存よりも多くのことを主張している。この法則が実際に埋めているのは、およそ 6N6N 個の座標をもつ空間の上で力学的に定義されたエネルギー UmicroU_{\mathrm{micro}} の存在と、わずかな巨視的変数だけに依存する状態関数としてのエネルギー UthU_{\mathrm{th}} の存在との間にある論理的な隙間である。これを示していこう。

次に、エネルギーの相加性が、孤立系のエネルギー保存から閉じた系の境界におけるエネルギー収支へ移るために必要な追加の仮定であることを示す。この最後の段階によってこそ、熱の認識論的・論理的な身分をようやく明確にできる。すなわち熱とは、閉じた系になされた巨視的な仕事を考慮したうえで、この収支に残るあらゆるエネルギー移動として定義される:Q  =def  ΔUWQ \equiv \Delta U - W

この講義ではまず閉じた系を考える。すなわち物質がその境界を横切らない系である。理由は単純で、開いた系では物質そのものが境界を越えてエネルギーを運び、エネルギー収支を複雑にするからである。開いた系の場合は本講義の後の方で扱う。

1. 微視的内部エネルギー

慣性系における古典的で、量子的でも相対論的でもない系を考える1。質量 mim_i、位置 ri\vec r_i、速度 vi\vec v_i をもつ NN 個の質点を考えよう。任意の時刻において、その微視的状態は次によって完全に定まる。

Note 1 : これらの場合には相応の扱いが必要であり、ここでは範囲外とする。
(r1,...,rN,p1,...,pN)Γ,ただしpi=mivi,\boxed{ (\vec r_1,...,\vec r_N, \vec p_1,...,\vec p_N) \in\Gamma, \, \text{ただし} \,\,\, \vec p_i=m_i\vec v_i, }

ここで Γ\Gamma は系の古典的な位相空間を表し、この単純なモデルでは次元は dimΓ=6N\dim\Gamma=6N である。分子間の引力や斥力など、すべての内力は保存力であると仮定する。これらには全ポテンシャルエネルギー Epint(r1,...,rN),ただしFiint=riEpintE_p^{\mathrm{int}}(\vec r_1,...,\vec r_N), \, \text{ただし} \,\,\, \vec F_i^{\,\mathrm{int}} =-\nabla_{\vec r_i}E_p^{\mathrm{int}} を対応させることができる。系はさらに、重力のような外部の保存力を受けていてもよい。同様に、これらには外部ポテンシャルエネルギー Epext(r1,...,rN)E_p^{\mathrm{ext}}(\vec r_1,...,\vec r_N) を対応させる。

系は選んだ基準系において必ずしも静止していない。任意の時刻にその質量中心と、場合により時間に依存する速度 V\vec V を定義できる。系の全質量と質量中心に対する相対速度

M=i=1Nmi,V=1Mi=1Nmivi,ui=viV,M=\sum_{i=1}^{N} m_i, \qquad \vec V=\frac{1}{M}\sum_{i=1}^{N} m_i\vec v_i, \qquad \vec u_i=\vec v_i-\vec V,

を導入すると、全運動エネルギーは次のように書ける。

Ec=12i=1Nmivi2=12i=1Nmi(V+ui)2.E_c=\frac12 \sum_{i=1}^{N}m_i \vec {v_i}^2 = \frac12 \sum_{i=1}^{N}m_i (\vec V + \vec u_i)^2.

平方を展開すると、imiui=0\sum_i m_i\vec u_i=\vec 0 であるから交差項が消えることがわかる。この結果は古典力学においてケーニヒの定理として知られている。

Ec=12MV2+12i=1Nmiui2+Vi=1NmiuiimiviMV=0E_c =\frac12MV^2 +\frac12 \sum_{i=1}^{N}m_i \vec{u_i}^2 +\vec V\mathbin{\cdot} \underbrace{\sum_{i=1}^{N} m_i\vec u_i}_{\sum_i m_i\vec v_i-M\vec V=\vec0}

したがって:

Ec=Ecmacro+Ecmicro\boxed{E_c = E_c^{\mathrm{macro}} + E_c^{\mathrm{micro}}}

(1)

全運動エネルギーは、このように常に2、質量中心の巨視的な運動エネルギーと、「無秩序な」残余の相対運動のエネルギーとに分解される。後者が熱運動と呼ばれるものを構成する。

Note 2 : 質量中心の並進が巨視的な運動エネルギーを尽くすとは限らない。系は全体としての回転という秩序ある運動をもちうる。その場合には次の分解を書くことができる。 vi=V+ω×(riR)+ui,\vec v_i = \vec V +\vec\omega\times(\vec r_i-\vec R) +\vec u_i', ここで R\vec R は質量中心の位置であり、ui\vec u_i' は全体の並進と回転を差し引いた残りの運動を表す。この残りの運動が質量中心に関して全角運動量を持たないように ω\vec\omega を選べば、 Ec=12MV2+12ωICMω+12i=1Nmiui2,E_c = \frac12 MV^2 + \frac12\vec\omega\cdot \mathbf I_{\mathrm{CM}}\vec\omega + \frac12\sum_{i=1}^{N}m_i {u_i'}^2, が得られる。ここで ICM\mathbf I_{\mathrm{CM}} は質量中心に関する慣性テンソルである。したがって巨視的な運動エネルギーと微視的な運動エネルギーの分離は依然として可能である。

各項をまとめると、全力学的エネルギーは最終的に次のように書ける。

Etot=Ecmacro+Epext+[Ecmicro+Epint(r1,...,rN)+Eintautres]内部エネルギー.E_{\mathrm{tot}} = E_c^{\mathrm{macro}} +E_p^{\mathrm{ext}} +\underbrace{\left[ E_c^{\mathrm{micro}} +E_p^{\mathrm{int}}(\vec r_1,...,\vec r_N) +E_{\mathrm{int}}^{\mathrm{autres}} \right]}_{\text{内部エネルギー}}.

注意:項 EintautresE_{\mathrm{int}}^{\mathrm{autres}} は、この最小限の質点モデルが記述しない内部自由度、たとえば分子の回転自由度や振動自由度の寄与をまとめたものである。それらを明示的にモデル化するならば、位相空間 Γ\Gamma をそれに応じて拡張しなければならない。

定義 1 (微視的内部エネルギー)
微視的内部エネルギーを Umicro=Ecmicro+Epint+EintautresU_{\mathrm{micro}}= E_c^{\mathrm{micro}} +E_p^{\mathrm{int}} +E_{\mathrm{int}}^{\mathrm{autres}} と書くことにする。これは位相空間の上に定義された関数 Umicro:ΓR,U_{\mathrm{micro}}:\Gamma\longrightarrow\mathbb R,

であり、一般に膨大な数の変数に依存し、次を満たす。

Etot=Ecmacro+Epext+Umicro.E_{\mathrm{tot}} =E_c^{\mathrm{macro}}+E_p^{\mathrm{ext}}+U_{\mathrm{micro}}.

(2)

以下では、巨視的に静止し、全体としての回転をもたず、外部ポテンシャルエネルギーが変化しない系を考える。したがって原点の選び方を除いて、次のように書ける。

Etot=Umicro.E_{\mathrm{tot}}=U_{\mathrm{micro}}.

これから、暫定的に UthU_{\mathrm{th}} と記す熱力学的エネルギーという状態関数を導入し、それをア・プリオリUmicroU_{\mathrm{micro}} と混同してはならない理由を理解しよう。

2. 熱力学的内部エネルギー

電熱プレートの上に置かれた水の入った鍋を系として考えよう。上で述べた単純な枠組みのままである。鍋は巨視的に静止しており、その外部ポテンシャルエネルギーは変化しない。したがって定数を除いて Etot=UmicroE_{\mathrm{tot}}=U_{\mathrm{micro}} と書ける。

水を加熱しても巨視的な運動エネルギーは現れず、鍋の位置も変わらない。あらゆる形態におけるエネルギーの厳密な保存を認めるならば、鍋に与えられたエネルギーは必然的にその微視的内部エネルギーの中に見出される。

ΔEtot=ΔUmicro.\Delta E_{\mathrm{tot}}=\Delta U_{\mathrm{micro}}.

熱力学はまさに同じ変化を記述しようとする。そこで、次を満たす量 UthU_{\mathrm{th}} を構成したい。

ΔUth=ΔUmicro.\Delta U_{\mathrm{th}}=\Delta U_{\mathrm{micro}}.

解決は自明に見える。熱力学的内部エネルギーを微視的内部エネルギーと同一視すれば十分ではないだろうか?

実はこの同一視は数学的に正しくない。というのも、この二つの関数は定義域が同じではないからである。UmicroU_{\mathrm{micro}} は位相空間の上に定義された関数であり、少なくとも 6N6N 個の変数に依存するのに対し、UthU_{\mathrm{th}} は平衡状態 X=(x1,...,xd)EX=(x^1,...,x^d)\in\mathcal E を記述するわずか dd 個の独立変数にしか依存しない(xix^i が温度、圧力、体積といった巨視的な状態変数であることを思い出そう)。

それでもなお、これら二つの関数は、Uth(X)U_{\mathrm{th}}(X) が巨視的なスケールで平衡状態 XX に対応する内部エネルギーを表すように結びついていなければならない。この次元の縮約がどのように働くのかを説明する必要がある。この問いは深く、この入門講義の枠を超えている。二つのスケールの間のこの移行を導出したいのであれば、微視的な記述が必要である。第 6.2 小節(範囲外)では、統計物理がこの問題をどう解決するかの大筋を述べる。

しかし厳密に熱力学的な観点からは、この微視的な同一視を正当化することはできず、わずかな巨視的変数だけの関数としての UthU_{\mathrm{th}} の存在は仮定されなければならない。それが熱力学第一法則の役割である。したがって第一法則はエネルギーの保存だけを主張するのではなく、それ以上に、このエネルギーが平衡において巨視的な空間 E\mathcal E の上に定義された状態関数によって表されうることを主張するのである。

dd 個の独立な状態変数 xix^i、たとえば理想気体に対する (TT, VV, NN) をひとたび選べば、関数 U(T,V,N)U(T,V,N) は状態空間の上のスカラー場を定める。座標 U,T,V,NU,T,V,N を張る拡張された空間において、そのグラフは超曲面であり、平衡面と呼ばれる。この講義の以後においてその役割は決定的である。図 1 を参照。

平衡状態の空間の上に定義されたエネルギー超曲面。
図 1. 平衡状態の空間の上に定義されたエネルギー超曲面。

最後に、UthU_{\mathrm{th}}UmicroU_{\mathrm{micro}} の違いをよく示す直接の帰結を指摘しておこう。ここで展開する平衡熱力学において、熱力学的内部エネルギーは一般に平衡状態の外では定義されない。たとえば気体の急激な変化の際には、もはや唯一の温度や唯一の圧力を割り当てられないことがありうるし、実際よくある。この場合、熱力学的内部エネルギー UthU_{\mathrm{th}} を直接結びつけることもできない。それは平衡状態の空間 E\mathcal E の上でしか定義されていないからである。問題はエネルギーが存在しなくなることではない。微視的なエネルギーは各瞬間に完全に定義されたままである。問題は、系がもはや点 XEX\in\mathcal E によって表されないということである。

以上を述べたうえで、次の節からは力学的な関数 UmicroU_{\mathrm{micro}} はもはや直接には現れない。熱力学的内部エネルギー UthU_{\mathrm{th}} を単に UU と書き、平衡状態の空間の上でのみ定義されていることを暗黙の了解とする。

注記 1 (平衡近傍への拡張)
平衡に十分近い系については、状態空間を大きくすることで熱力学的な記述を拡張できることを指摘しておく。たとえば、すでに見た二つの温度源の間に置かれた金属棒の例では、唯一の温度ではなく定常な局所温度場 T(x)T(x) を導入し、内部エネルギーの局所密度を定義できる。全エネルギーはそのときこれらの場の汎関数となる。ここではこの平衡近傍熱力学と呼ばれるものを展開しない。

3. 熱力学第一法則

前節は、いま述べる内容のための下地を十分に整えてくれた。

閉じた系に対する熱力学第一法則
すべての閉じた熱力学系に対して、状態関数 U:ER,U:\mathcal E\longrightarrow\mathbb R,

が存在する。これを内部エネルギーと呼び、加法定数を除いてのみ定まり、少なくとも二回微分可能であって、任意の平衡状態 XEX\in\mathcal E に対して次を満たす。

Etot=Ecmacro+Epext+U(X).E_{\mathrm{tot}} = E_c^{\mathrm{macro}} +E_p^{\mathrm{ext}} +U(X).
(3)

第一法則はエネルギーの保存も主張する。すなわち孤立系に対して、

ΔEtot=0.\Delta E_{\mathrm{tot}}=0.

関数 UU の正則性もまた自明でない仮定である。実際上は UUC\mathcal C^\infty と見なしてよい。等式 3 に関わるポテンシャルエネルギーのために UU は定数を除いてしか定まらないので、物理的な意味をもつのはその変化だけである。簡単のため、以下では巨視的な運動エネルギーとポテンシャルエネルギーが変化しない過程に限る。そのとき ΔEtot=ΔU\Delta E_{\mathrm{tot}}=\Delta U となる。

ここで、環境とエネルギーは交換できるが物質は交換できない閉じた系 SS を考えよう。全体系

S=S+ExtS'=S+\mathrm{Ext}

は常に孤立系の中に含めることができる。そこで第一法則は次を要求する。

ΔES=0.\Delta E_{S'}=0.

エネルギーが相加的である、すなわち両者の相互作用エネルギーが無視できるという追加の仮定のもとで、ES=ES+EExtE_{S'}=E_S+E_{\mathrm{Ext}} が成り立ち、したがって

ΔES=ΔEExt.\Delta E_S=-\Delta E_{\mathrm{Ext}}.

相加性の仮定は、こうしてエネルギーの保存を部分系の間の収支という形に翻訳する。系が得たエネルギーはすべてその環境が失い、逆もまた同様である。あとは、このエネルギーが系の境界をどのように横切って移動しうるのかを明確にすればよい。二つの移動の様式を区別しよう。すなわち、力学に由来する独立な定義をもつ仕事 WW と、 QQ である。

そのときエネルギー収支は次の形をとる。

ΔES=ΔU=Q+W\boxed{ \Delta E_S = \Delta U = Q+W }
(4)

ここで WW は系が受け取った仕事を、QQ は系が受け取った熱を表す。

エネルギーの相加性の仮定が自明でないことに注意しよう。初等的な水準では、通常の熱力学系についてはそれを単に認めてよい。しかしそれは常に成り立つわけではない。とりわけ長距離相互作用は、エネルギーの相加性と示量性を同時に破りうる。第 6.4 節を参照。

あとはこの二つの移動の様式を定義すればよい。しかしその前に、本講義の以後ずっと用いる、ごく標準的な符号の約束を述べておかねばならない。

要点 (銀行家の符号の約束)
本書全体を通じて、エネルギーの移動は系が受け取るときを正、系が与えるときを負とする。したがって、 Q>0, W>0Q>0,\ W>0

は系が受け取ったエネルギーに対応し、

Q<0, W<0Q<0,\ W<0

は外へ与えられたエネルギーに対応する。

4. 仕事、熱、エネルギー収支

4.1. 仕事

仕事の利点の一つは、熱力学とは独立な定義を力学から受け継いでいることである。外力 Fext\vec F_{\mathrm{ext}} が系の境界上のある点に作用し、その点が drd\vec r だけ変位するとき、この力がなし、したがって系が受け取る微小仕事は

δW=Fextdr.\delta W = \vec F_{\mathrm{ext}}\cdot d\vec r.

広がりをもつ系については、その境界の上で仕事をするすべての外力の寄与を足し合わせなければならない。

実際上もっともよく用いる場合は、可動な壁をもつ容器に入れられた流体の場合であり、これは圧力の仕事を評価することにつながる。断面積 S\mathcal S の平らなピストンによって円筒に閉じ込められた気体を考えよう。図 2 を参照。軸 xx は気体の外側に向けてとる。ピストンが dxdx だけ変位すると、体積は

dV=Sdx.dV=\mathcal S\,dx.

だけ変化する。

気体の境界に及ぼされる外部圧力を PextP_{\mathrm{ext}} と記そう。気体が受ける力は内側を向いており、

Fext=PextSex.\vec F_{\mathrm{ext}} = -P_{\mathrm{ext}}\mathcal S\,\vec e_x.

に等しい。

したがって気体が受け取る微小仕事は

δWpression=Fext(dxex)=PextSdx,\delta W_{\mathrm{pression}} = \vec F_{\mathrm{ext}}\cdot(dx\,\vec e_x) = -P_{\mathrm{ext}}\mathcal S\,dx,

すなわち

δWpression=PextdV\boxed{ \delta W_{\mathrm{pression}}=-P_{\mathrm{ext}} dV }
(5)

銀行家の規則によって符号をただちに確かめられる。圧縮のときは dV<0dV<0 ゆえ δW>0\delta W>0:気体は仕事を受け取る。膨張のときは dV>0dV>0 ゆえ δW<0\delta W<0:気体は環境に仕事を与える。

この式に現れるのが系にとって外部の圧力であり、一般に気体自身の圧力に等しい理由はないことに注意しよう。とりわけ、内部圧力が一般には定義さえされない急激な非平衡過程の場合には、この点を心に留めておくことが重要である。

断面積 S のピストンが変位するときに気体が受け取る仕事。膨張のときは dV>0 である一方、外部圧力が及ぼす力は変位と逆向きである。したがって受け取る仕事は負である。
図 2. 断面積 S\mathcal S のピストンが変位するときに気体が受け取る仕事。膨張のときは dV>0dV>0 である一方、外部圧力が及ぼす力は変位と逆向きである。したがって受け取る仕事は負である。

積分すれば、外部圧力が既知である有限の変化 ABA \to B について次の式が得られる。

Wpression[AB]=ABPextdV\boxed{ W_{\mathrm{pression}}[A \to B] = -\int_{A}^B P_{\mathrm{ext}} dV }
(6)

この積分は、全仕事を計算するには変化の全体にわたって PextP_{\mathrm{ext}} の値を知らねばならないことを示している。ここから重要な帰結が導かれる。仕事は一般に、初状態と終状態 AABB だけからは決定できない。その表式はたどった過程に依存する

ここで指摘しておくべき特別な場合が三つある。

  1. PextP_{\mathrm{ext}} が一定ならば、 Wpression[AB]=PextΔVW_{\mathrm{pression}}[A\to B] = -P_{\mathrm{ext}}\Delta V
  2. 定積変化に沿っては dV=0dV=0 であり、圧力の仕事はゼロである:δWpression=0\delta W_{\mathrm{pression}}=0 かつ Wpression[AB]=0W_{\mathrm{pression}}[A\to B]=0
  3. 真空への膨張では Pext=0P_{\mathrm{ext}}=0 ゆえ、圧力の仕事も同様にゼロである。

注記 2 (その他の形の仕事)
圧力の仕事は数ある例の一つにすぎない。引き伸ばされる糸について、LL をその長さ、FextF_{\mathrm{ext}} を受け取る引張力とすれば、 δW=FextdL.\delta W=F_{\mathrm{ext}}\,dL.

面積 AA を増加させる界面については、表面張力 γ\gamma が通常の条件下で

δW=γdA.\delta W=\gamma\,dA.

という形の仕事を与える。

目に見える機械的な変位なしにエネルギーを移すこともできる。たとえば電荷 dqdq を電位差を横切って移動させると、

δWeˊlec=Vextdq,\delta W_{\mathrm{élec}}=V_{\mathrm{ext}}\,dq,

という形の電気的な仕事が生じうる。ここで VextV_{\mathrm{ext}} は外部の電位であり、この文脈では体積と混同しないようにしばしば Φext\Phi_{\mathrm{ext}} と記される。ここでも dqdq の符号は銀行家の約束に従って選ばれる。

第 6.3 小節では外部の仕事の一般形を明確にする。

4.2. 熱

仕事がエネルギー移動の唯一の様式ではありえない。というのも、剛体の容器に入れた気体を加熱すればその温度とエネルギーは増加するが、外部からの仕事はまったくなされていないからである。したがってエネルギーは、仕事の形で移動することなく系の境界を横切ったのである。この第二のエネルギー移動の様式を、あるいは熱移動と呼ぶ。

ここで採用する構成において(同等の別の構成については第 6.1 小節を参照)、二つの平衡状態 AABB を結ぶ物理的な変化 ABA \to B について、系が受け取ったすべての仕事を同定したのち、受け取った熱は次のように定義される。

Q[AB]  =def  U(B)U(A)W[AB].\boxed{ Q[A \to B] \equiv U(B)-U(A)-W[A \to B]. }
(7)

こうして先に述べた関係

ΔU=Q+W\boxed{ \Delta U=Q+W }

が改めて得られる。仕事が変化 ABA \to B の道筋に依存し、エネルギーの変化 ΔU=U(B)U(A)\Delta U = U(B)-U(A) がそれに依存しない以上、やり取りされる熱量も必然的にそれに依存する

熱力学の言葉では、内部エネルギー UU は(仮定により)状態関数であるが、仕事 WW と熱 QQ はそうではない、と言う。これが意味するところを改めて述べよう。WWQQ がある時刻における系の状態の関数でありえないということは、系が「ある量の仕事や熱」をもつのではなく、むしろ WWQQ が任意の変化 ABA \to B の間にその境界で生じるエネルギー移動である、ということを意味する。これはまさに、第2講で詳しく述べたさまざまな実験に照らして、我々が得たかったものである。

循環過程 AAA \to A においては系が初めの状態に戻ることに注意しよう。したがって

ΔUcycle=0.\Delta U_{\mathrm{cycle}}=0.

そこで第一法則は次を要求する。

Qcycle+Wcycle=0.\boxed{Q_{\mathrm{cycle}}+W_{\mathrm{cycle}}=0.}

したがって循環機関は、環境から等量のエネルギーを受け取ることなしに無限に仕事を供給することはできない。これが歴史的に第一種永久機関の不可能性と呼ばれるものである。

4.3. 第一法則の微分形式

関係式

ΔU=Q+W\Delta U=Q+W

は二つの平衡状態 AABB を結ぶ。これは常に成り立つ。というのも、途中のすべての状態が平衡状態の空間 E\mathcal E の点として表されうるとは仮定していないからである。変化が急激であれば、U(A)U(A)U(B)U(B) が完全に定義されたままであっても、系は E\mathcal E の上に定義された平衡面を一時的に離れることがある。

ここで特別な場合を詳しく述べねばならない。変化が準静的であると仮定しよう。そのとき系は定義により各瞬間に平衡状態にある。この変化は、無限に近い平衡状態からなり平衡面を離れない

γE\gamma\subset\mathcal E

によって表せる。UUE\mathcal E の上で微分可能な状態関数であるから、無限に近い二つの状態の間のその変化は完全微分であり、dUdU と記される。

対応する熱と仕事の微小な移動をそれぞれ δQ\delta QδW\delta W と記す。この場合、第一法則は次の形をとる。

dU=δQ+δW\boxed{ dU=\delta Q+\delta W }
(8)

dUd UδW\delta W あるいは δQ\delta Q の記法の違いは単に見かけ上のものではない。それは前者が完全微分であり、後の二つがそうでないことを示している。これは、エネルギーの変化がたどった道に依存しないのに対し QQWW はそれに依存する、という事実の数学的な翻訳である。次節ではこの決定的な点をよく理解するために必要な数学を、初学の読者のために復習する。

仕事が圧力の仕事だけであり、PextP_{\mathrm{ext}} を系の圧力 PP と同一視できる特別な場合には、この関係は次のようになる。

dU=δQPdV.dU=\delta Q-P\,dV.
要点 (第一法則の二つの定式化)
その積分形 ΔU=Q+W\Delta U=Q+W

は二つの平衡状態を結び、途中の過程が準静的でない場合にも使える。その微分形

dU=δQ+δWdU=\delta Q+\delta W

は、平衡状態の空間の中の道によって表せる準静的な変化に対してのみ成り立つ。

5. 完全微分形式と不完全微分形式

本節は、一方の dUdU と他方の δQ\delta QδW\delta W の区別に正確な意味を与えるために必要な最小限の数学をまとめたものである。すでに知っている読者は飛ばしてよい。

5.1. 関数の微分

変数 x1,...,xdx^1,...,x^d の関数 ff が微分可能であるとしよう。その微分とは表式

df=i=1dfxidxi,df=\sum_{i=1}^{d}\frac{\partial f}{\partial x^i}\,dx^i,

であり、点 XX から近傍の点 X+dXX+dX へ移るときの ff の変化を測る。重要なのは次の点である。点 AA から点 BB へ至る道 γ\gamma をたどるとき、これらの微小変化をすべて足し合わせれば

γdf=f(B)f(A).\int_\gamma df=f(B)-f(A).

が得られる。とくに閉じた道の上では

df=0.\oint df=0.

5.2. 微分形式

いま同じ型の表式を考えよう。ただし、これを関数の微分と混同することは決して許されない。任意の関数 Ai(x1,...,xd)A_i(x^1,...,x^d) を与え、ω\omega を次で定義する。

ω=i=1dAidxi.\omega=\sum_{i=1}^{d}A_i\,dx^i.

このような対象を微分形式と呼ぶ。関数 AiA_i の表式が与えられていれば、道 γ\gamma に沿ってこれを問題なく積分できる。しかし、ω\omega がその微分となるような関数 ff、すなわちすべての ii について Ai=f/xiA_i=\partial f/\partial x^i を満たす ff が存在する保証はどこにもない。

  • そのような関数 ff が存在するとき、その形式は完全であると言い、ω=df\omega = d f と記す。 そのとき積分 γω=γdf=f(B)f(A)\int_\gamma\omega = \int_\gamma df = f(B)-f(A) は端点だけに依存する。
  • そうでないとき、その形式は不完全であると言い、いかなる関数の微分でもないことを思い起こさせるために、ω\omega ではなく δω\delta\omega と記す。(数学者なら単にこの形式を ω\omega と書くだろう。)

5.3. シュワルツの判定条件

微分形式が完全であるかどうかを判定する簡単な条件が、シュワルツの判定条件である。簡単のため二変数の場合に限ろう(一般化はただちにできる)。

ω=A(x,y)dx+B(x,y)dy.\omega=A(x,y)\,dx+B(x,y)\,dy.

とする。ω\omega が完全であれば A=f/xA=\partial f/\partial x かつ B=f/yB=\partial f/\partial y となるはずであり、したがって

Ay=2fyx=2fxy=Bx,\frac{\partial A}{\partial y} =\frac{\partial^2 f}{\partial y\,\partial x} =\frac{\partial^2 f}{\partial x\,\partial y} =\frac{\partial B}{\partial x},

が成り立つ。二回連続微分可能な関数については微分の順序が問題にならないからである。こうして便利な検査法が得られる。

AyBx    ω は完全でない。\boxed{ \frac{\partial A}{\partial y}\neq\frac{\partial B}{\partial x} \;\Longrightarrow\; \omega\ \text{は完全でない。} }

熱力学では起こらない細かな事情を除けば、逆も成り立つ。

5.4. 圧力の仕事の例

微小な移動量が準静的な道に沿って状態変数によって表せるとき、δQ\delta QδW\delta WE\mathcal E の上の微分形式となり、それらが完全かどうかを確かめられる。

例として、Pext=PP_{\mathrm{ext}} = P のもとで準静的な変化をする理想気体 nn モルを考えよう。このとき δW=PdV\delta W=-P\,dV である。ところが P=nRT/VP=nRT/V であるから、微小仕事は

δW=PdV=nRTVdV=0×dTnRTVdV,すなわちA=0,B=nRTV,\delta W=-P dV = -\frac{nRT}{V}\,dV = 0\times dT-\frac{nRT}{V}\,dV, \qquad\text{すなわち}\qquad A=0, \quad B=-\frac{nRT}{V},

と書ける。ここで x=Tx=Ty=Vy=V と対応させた。これはたしかに状態空間の上に定義された微分形式である。シュワルツの判定条件は

AV=0,BT=nRV0.\frac{\partial A}{\partial V}=0, \qquad \frac{\partial B}{\partial T}=-\frac{nR}{V}\neq0.

を与える。二つの交差微分は等しくない。微小仕事はたしかに完全微分ではない。したがってこの場合、δW\delta W をその変化とするような状態関数 W(T,V)W(T,V) は存在せず、二つの状態の間に受け取る仕事はまさしくたどった道に依存する。同じ論法が δQ=dUδW\delta Q=dU-\delta W にも当てはまる。dUdU は完全であり δW\delta W はそうでないから、その差もまた完全ではありえない。

要点 (なぜ dd ではなく δ\delta か)
dUdU は完全微分である。その変化は初状態と終状態だけに依存し、循環に沿った積分はゼロである。δQ\delta QδW\delta W は不完全な形式である。それらはたどった道に依存し、δW=dW\delta W = d W かつ δQ=dQ\delta Q = dQ となるような系の状態の関数 QQWW も存在しない。

6. さらに進んで

本節では、これまでに予告したより進んだ考察をまとめる。最初に読むときは飛ばしてよい。

6.1. 第一法則の二つの可能な構成

ここでは次の論理的な順序を選んだ。状態関数 UU の存在は第一法則によって仮定され、仕事は力学によって独立に定義され、次いで熱が収支

Q=ΔUW.Q=\Delta U-W.

によって定義される。ここから注目すべき事実が導かれる。変化が断熱的であれば ΔU=W\Delta U = W であり、これは断熱仕事 WadW_{\mathrm{ad}} の方が常にたどった道に依存しないことを要求する。

このことは、内部エネルギーの別の可能な構成を指し示している。それは操作的であり、したがっていかなる微視的な考察にも依存せず、この場合は歴史的な歩みにより近い。まず、推論の循環を避けるために、量 QQ を先に導入することなく断熱変化を特徴づけるところから始める。

その着想は巧妙である。移動そのものではなく装置を定義するのである。ある壁が囲む系の状態が、外部の力学的な座標、すなわちピストン、撹拌器、あるいは電流を動かすことによってしか変えられないとき、その壁は断熱的であると言う。そこで検査は直接的である。これらの座標を固定したまま、たとえば容器を氷浴に浸したり炎に近づけたりして、外部を任意に変えてみる。系のどの状態変数も変化しないなら、その壁は断熱的である。

この特徴づけは状態と力学的な変位だけを用い、エネルギーの移動を決して用いない。したがってそれは熱というあらゆる概念に先立つ。実際上は、よい断熱によって、あるいは熱緩和時間に比べて速く操作することによって、これに近づける。

一方、受け取る仕事は純粋に力学的または電気的な仕方で測定される。これこそがジュールの実験の経験的な内容のすべてである。これらの実験は1843年から1850年にかけて行われ、以後ますます高い精度で繰り返されてきた。高さ hh から落ちる質量 mm は、いかなる熱移動も介さずに mghmgh を供給する。彼が考案し上述の意味での断熱容器の中で実現したすべての装置、すなわち羽根つき撹拌器、加熱抵抗、気体の圧縮は、いずれも同じ結論に彼を導いた。同じ仕事を与えれば同じ状態変化が生じるのであり、これが彼にとって「熱の仕事当量」を構成した。第2講で説明したとおりである。こうしてジュールは、断熱仕事がたどった道に依存しないことを(その事実の実験的な示唆として)最初に示した人物となった。

この実験的事実を公準に高めることを提案したのが、1909年のカラテオドリ [2]、そしてとりわけ1921年のボルンである。これを認めるならば、

U(B)U(A)=Wad(AB),U(B)-U(A)=W_{\mathrm{ad}}(A\to B),

と置くことによって状態関数 UU構成し、次いで我々がしたように、一般の変化について差によって熱を定義することが容易になる。

この構成は多くの著作、とりわけピパード [3] とキャレン [4] の著作で採用されている。それは我々の構成と同等である。UU を仮定してから断熱仕事の道への非依存性を導くことは、この仕事の非依存性を仮定してから UU を構成することと同等である。ジュールからカラテオドリとボルンに至るこの概念的な発展の詳しい物語は、ローゼンバーグ [5] に見出される。

6.2. 微視状態と巨視状態

上に述べた操作的な道をとらないのであれば、第 2 節は根本的な問いを明るみに出した。6N6N 個の座標の空間の上に定義されたエネルギーが、どうしてわずか数個の巨視的変数の関数に要約されうるのか? 統計物理は二つのスケールの間のこの移行を構築する。その大筋は以下のとおりである。

統計物理では、平衡状態の空間 E\mathcal E の点 XX巨視状態と呼ぶ。同一の巨視状態には、一般に、同じ巨視的な拘束と両立する厖大な数の微視状態 γΓ\gamma\in\Gamma が対応する。

そして、これらの微視状態のエネルギーが、NN\to\infty という熱力学的極限においてその平均値のまわりに極度に集中して分布することが示される。言い換えれば、エネルギーの相対的な揺らぎがゼロに向かう。

σ(Umicro)Umicro0.\frac{\sigma(U_{\mathrm{micro}})} {\langle U_{\mathrm{micro}}\rangle}\longrightarrow0.

そこで Uth(X)=UmicroXU_{\mathrm{th}}(X)= \langle U_{\mathrm{micro}} \rangle_X と置くことができる。平均はその巨視状態と両立するすべての微視状態にわたる。上で述べた集中こそが、この定義に意味を与えるものである。系は各瞬間に唯一の微視状態を占めるのであって平均を占めるのではない。ほとんどすべての微視状態が同じエネルギーをもつからこそ、巨視状態エネルギーについて語ることができ、同一の準備が同じ測定値を与えるのである。

イメージをつかむために、一つの巨視状態には、微視的には異なるが巨視的には区別できない微視状態の非常に大きな集合が対応する、と想像してもよい。この描像は位相空間の上の同値類に似ている。もっとも、正確な構成はいくらか異なり、位相空間の上の確率分布を用いることになる。

次の図はこの手続きをまとめたものである。

同一の巨視状態 X は厖大な数の微視状態と両立する。それらのエネルギーは、我々のスケールでは、同じ一つの巨視的な値のまわりに極度に集中している。これによって U_ th(X) を定義できる。
図 3. 同一の巨視状態 XX は厖大な数の微視状態と両立する。それらのエネルギーは、我々のスケールでは、同じ一つの巨視的な値のまわりに極度に集中している。これによって Uth(X)U_{\mathrm{th}}(X) を定義できる。

6.3. 一般化された仕事

上で出会ったさまざまな形の仕事は共通の構造をもつ。圧力の仕事について

δW=PextdV\delta W=-P_{\mathrm{ext}}\,dV

糸の伸長について

δW=FextdL\delta W=F_{\mathrm{ext}}\,dL

さらに表面の仕事について

δW=γextdA\delta W=\gamma_{\mathrm{ext}}\,dA

という形を見た。いずれの場合も微分される量が示量的であり、その前因子が示強的であることに注意しよう。一般化して、力学的な仕事の項はすべて次の形に書くことにする。

δW=iνiextdxi,\boxed{ \delta W=\sum_i \nu_i^{\mathrm{ext}}\,dx^i, }

ここで νiext\nu_i^{\mathrm{ext}} は座標 xix^i(一般に示量的)に共役な外部の一般化された力(一般に示強的)を表し、符号は我々の約束に対応する。

外部の一般化された力が系の対応する熱力学的な力に一致する場合(例:Pext=PP_{\mathrm{ext}}=P)、第一法則は次の形に書ける。

dU=δQ+iνidxi.\boxed{ dU=\delta Q+\sum_i\nu_i\,dx^i.}

6.4. 内部エネルギーの示量性

前の講義では、内部エネルギーを示量的な量として提示した。示強変数を固定したまま均一な系の大きさを因子 λ\lambda 倍すれば、そのエネルギーも λ\lambda 倍になる。

UλU.U\longrightarrow\lambda U.

この性質は第一法則によって仮定されるものではない。これはしばしば置くことになる追加の仮定であるが、常に成り立つわけではない。実際、系の構成要素の間に長距離力(とりわけ重力)が存在すれば、その示量性が損なわれることが予想される。

これを詳しく見よう。次元 DD の空間の中に一定の密度で分布した NN 個の構成要素を考え、それらの相互作用ポテンシャルエネルギーが

Ep(r)1rα.E_p(r)\sim\frac{1}{r^\alpha}.

のようにふるまうと仮定する。密度 ρ=N/V\rho=N/V を固定すると、系の典型的な線寸法は

LN1/D.L\sim N^{1/D}.

のように増大する。そこで全相互作用エネルギーを三段階で評価しよう。

隣人を数える。

一つの構成要素を固定し、距離 rrr+drr+dr の間に他の構成要素がいくつあるかを問う。それは密度に対応する殻の体積を掛けたものである。

dn(r)=ρSDrD1dr,dn(r)=\rho\,S_D\,r^{D-1}\,dr,

ここで SDS_D は次元 DD における単位球面の面積を表し、三次元では 4π4\pi である。

距離について足し合わせる。

これらの隣人はそれぞれ Ep(r)rαE_p(r)\sim r^{-\alpha} だけ寄与する。したがって一つの構成要素と他のすべての構成要素との相互作用エネルギーは

uρaLrD11rαdr=ρaLrD1αdr.u\sim\rho\int_a^L r^{D-1}\,\frac{1}{r^{\alpha}}\,dr =\rho\int_a^L r^{D-1-\alpha}\,dr .

に等しい。下限 aa は最小接近距離であり、それより近くでは 1/rα1/r^\alpha の法則が成り立たなくなり、また r0r\to0 で積分が発散するのを防ぐ。上限 LL は系の大きさである。それより遠くに隣人はいない。

構成要素について足し合わせる。

構成要素の数を掛ける(各対を二度数えないように二で割る)。

UintN2uNρaLrD1αdr,U_{\mathrm{int}}\sim\frac{N}{2}\,u \sim N\rho \int_a^L r^{D-1-\alpha}\,dr,

ここで幾何学的な因子 SD/2S_D/2 は桁の見積もりの中に吸収した。そこで積分の値を見よう。

aLrD1αdr=[rDαDα]aL(αD).\int_a^L r^{D-1-\alpha}\,dr =\left[\frac{r^{D-\alpha}}{D-\alpha}\right]_a^L \qquad(\alpha\neq D).

α>D\alpha>D ならば指数は負であり、積分は LL\to\infty で収束し、下限によって支配される。したがって定数 aDα/(αD)a^{D-\alpha}/(\alpha-D) の程度である。α<D\alpha<D ならば逆に上限によって支配され、LDα=N1α/DL^{D-\alpha}=N^{1-\alpha/D} の程度になる。最後に α=D\alpha=D の場合、原始関数は対数であり積分は ln(L/a)=1DlnN\ln(L/a)=\tfrac1D\ln N に等しい。これらを代入して次を得る。

Uintcste×{N,α>D,NlnN,α=D,N2α/D,α<D.\begin{aligned} U_{\mathrm{int}} \sim \text{cste} \times \begin{cases} N, & \alpha>D,\\[1mm] N\ln N, & \alpha=D,\\[1mm] N^{\,2-\alpha/D}, & \alpha<D. \end{cases} \end{aligned}

物理的な解釈は明快である。α>D\alpha>D のときは短距離力である。EpE_p の減衰が隣人の数の増加に勝つ。各構成要素はすぐ近傍しか感じないので、そのエネルギーは系の大きさに依存せず、全体は NN に比例する。すなわち示量的である。

αD\alpha\leq D のときは逆である。はるかに数の多い遠くの隣人が勝ち、各構成要素は系全体を感じ、それ自身のエネルギーは系の大きさとともに増大する。全体はそのとき NN より速く増大し、示量性は失われる。したがって十分に短距離の相互作用は自然に示量的なエネルギーへ導くのに対し、長距離相互作用はこの性質を破壊しうる。

たとえば三次元のニュートン重力では

Ep(r)1r,D=3,α=1,E_p(r)\sim-\frac1r, \qquad D=3, \qquad \alpha=1,

であり、密度を固定すれば先の見積もりは

UgravN5/3.U_{\mathrm{grav}}\sim-N^{5/3}.

を与える。桁の議論だけでは得られない負の符号は、重力が引力であることを反映している。

エネルギーの示量性は本書の以後で大いに用いられる。重力はもちろん普遍的だが、非常に弱い力である。接触させた二つの気体の体積を調べるのであれば、実際上これを無視できる。他方、気体内部のファン・デル・ワールス型の力は非常に速く(1/r61/r^6 のように)減衰するので、この場合の示量性はほぼ厳密である。

これに対し、いわゆる自己重力系(恒星、銀河など)を記述するには重力がもちろん本質的な要素であり、エネルギーの示量性が必然的に失われる以上、熱力学的な解析全体を最初からやり直さねばならない。そのため、これらの系の熱力学はほとんど独立した一つの学問であり、思いがけないふるまいを示す。たとえば、エネルギーを放射する恒星は冷えるのではなく温まる。これについては本書の発展的な部分で改めて扱う。

7. まとめ

本講義が据えたこと、そして以後ずっと役立つことをまとめておこう。

  • 第一法則は状態関数 U:ERU:\mathcal E\longrightarrow\mathbb R、すなわち内部エネルギーの存在を仮定し、孤立系のエネルギーの保存を主張する。
  • 本書の以後では、それが示量的でもあり、すべての変数について十分に連続かつ微分可能であると仮定する。
  • エネルギーの相加性によって補われることで、第一法則は保存を、閉じた系とその環境との間のエネルギー収支として定式化することを可能にする。仕事は力学によって独立に定義されるので、熱はそのとき残りの移動として定義され、有用な式 ΔU=Q+W\Delta U=Q+W が導かれる。
  • 準静的な変化については、このエネルギー収支は微分形 dU=δQ+δWdU=\delta Q+\delta W をとり、そのうち完全微分であるのは dUdU だけである。

開いた系の解析は本講義の別の場所で扱う。

8. 参考文献

第一法則の概念的な発展について、ジュールからカラテオドリとボルンまでは、ローゼンバーグ [5] を見よ。古典的な提示については、ピパード [3] とキャレン [4] を見よ。

  1. J. P. Joule, “On the Mechanical Equivalent of Heat,” Phil. Trans. R. Soc. Lond. 140, 61—82 (1850)
  2. C. Carathéodory, “Untersuchungen über die Grundlagen der Thermodynamik,” Math. Ann. 67, 355—386 (1909)
  3. A. B. Pippard, Elements of Classical Thermodynamics for Advanced Students of Physics, Cambridge University Press (1957)
  4. H. B. Callen, Thermodynamics and an Introduction to Thermostatistics, 2nd ed., Wiley (1985)
  5. R. M. Rosenberg, “From Joule to Caratheodory and Born: A Conceptual Evolution of the First Law of Thermodynamics,” J. Chem. Educ. 87, 691—693 (2010)
  6. A. M. Steane, “First Law, internal energy,” chap. 7 in Thermodynamics: A Complete Undergraduate Course, Oxford University Press (2017)
  7. E. A. Gislason and N. C. Craig, “Cementing the foundations of thermodynamics: Comparison of system-based and surroundings-based definitions of work and heat,” J. Chem. Thermodynamics 37, 954—966 (2005)