こちらのバージョンをご利用になると鮮明に印刷されます。 目次付フレーム版 詳細目次付フレーム版 NOフレーム版
熱力学に於ける絶対温度の意味を説明します。温度はとても難しい量でエントロピー(これも難しい)と表裏一体です。説明にはかなりの準備が必要ですが、なにとぞ御容赦下さい。
温度を定めるに当たって熱平衡という概念が導入される。これは熱学の基礎になる概念で温度(次で説明)の違う二物体を接触させると、温かいものは冷え、冷たい方は温まり、温度の変化が両方とも止まった状態になる。その場合を熱平衡にあると言う。
このとき物体Aが物体Bと物体Cのそれぞれと熱平衡にあるならば、BとCも熱平衡である。これは熱平衡の概念に附随する法則と言ってもよいし、経験的事実と言っても良い。
この事実があるからこそ、二物体BとCの温度をそれらが離れて存在しても比較することができる。つまり温度計Aを作り、まず温度計AとBと接触させて熱平衡状態を実現し、次にAとCを接触させて熱平衡状態を実現する。温度計Aを構成する物質の体積の観測によってBとCが熱平衡にあるか、あるいはそのいずれがより温かいか冷たいかを決定するのである。
温度とは何かと改めて聞かれると、これは本当に難しい量です。その当たりはおいおいに説明するが、当面は簡単に考える。
温度は、我々の涼しい・温かい、冷たい・熱いの感覚に関係している。それを量的に表すものとして温度計というものが作られた。それらは、容器に詰めた液体や気体の膨張を測ってその体積の増減でもって温度の高低を示す。通常の温度計は膨張量をより詳細に測定するために物質溜めに細管を繋いで膨張の有様を観測しやすくしている。つまり、温度の高低は物体の体積の変化量に比例していると考えられおり、それらの量の変化を引き起こすものが温度の実体である。
温度は、物質の示す二つの基準状態と熱平衡にある温度計について、それぞれの平衡時における温度計構成物体の体積を量り、二状態間の体積差を等分割する事によって目盛りを決めるような量である。つまり任意の物体の温度は、その物体と熱平衡状態にある温度計の構成物質の示す体積を量ることにより決定する。
しかし、よく考えてみると温度がどの様なメカニズムで気体や液体の膨張を引き起こすのかは何も解っていない。またその膨張によって温度を定義し、その値を知ることはできるかも知れないが、温度の実体については何も教えてくれない。
1.ファーレンハイトの温度目盛り t
(Fahrenheit、中国語で華倫海特、単位は°Fまたはファーレンハイト氏が作ったと言う意味で華氏〜度と言う)
1724年ファーレンハイトは前記の二つの基準状態として寒剤により達成される状態と、動物の体温を選んだ。つまり、当時人工的に造られる最低温度(氷と塩を混合した寒剤による)を0度とし人の体温を96度としてその間の温度計の体積差を12×8=96等分割して温度を定めた。これはほぼ、1気圧の元での氷りの融点を32°Fとし、沸点を212°Fとすることに相当する。
2.セルシウスの温度目盛り t
(Celsius、中国語で摂爾修斯、単位は℃またはセルシウス氏が作ったと言う意味で摂氏〜度と言う)
1742年セルシウスは二つの基準状態として、1気圧の元での水の凝固点(氷と水が共存する状態)と、1気圧の元での水が沸騰している状態を選んだ。つまり、1気圧下での水の凝固点を0℃、沸点を100℃として、温度計の体積変化を100に等分割して温度を定めた。
これ以外にも様々の温度目盛りが考案されたが、いずれも二つの物理的状況と熱平衡にある温度計の体積を量り、二状態間の体積差を等分割して温度目盛りを定めるものである。(ダンネマン「大自然科学史」第6巻(三省堂)等参照)
ここで問題になるのは、例えばセルシウスの目盛りを例にして説明すると、二種類の温度計物質の示す温度は0℃と100℃では一致させることはできても、その間では決して一致しないと言うことです。つまり一般に0℃と100℃の間の体積変化の様子は物質によりそれぞれ異なる。つまり同じ温度の熱浴に浸けて膨張率を測定しても、例えば30℃付近では物質Aの方の膨張率が物質Bより大きいが、70℃付近では逆にBの方がAより大きいと言うようなことが起こる。
そのとき、温度計の物質によらない絶対的な温度を定めることはできないのだろうか?これは”大いなる疑問?”ですが、その解決は可能です。熱や温度についての深い考察が必要ですが、その問題は1848年のケルビンKelvin(ウィリアム・トムソンWilliam
Thomson)による絶対温度の定義から始まって、最終的にクラウジウスやケルビンによって解決された。
この問題は、現実には存在しない理想気体の膨張率を基準にして上記のグラフを書くことができるかという事と同じです。つまり現実に存在する気体の熱膨張率から、現実には存在しない理想気体の熱膨張率を決めることができるかという問題です。
3.ケルビン(ウィリアム・トムソン)の絶対温度目盛り T
(単位はK(ケルビンと読む)とする。このとき”°”を付けないことに注意)
これはカルノーサイクル(ある種の熱機関)を用いて温度を定義するものです。この定義によると温度測定に用いる物質に依存しない温度目盛り(絶対温度目盛り)を打ち立てることができる。また、この定義によると必然的に温度がそれ以下に下がることがない基準点(それを0Kとする)が存在することが解る(6.(2)2.参照)。
そのため絶対温度については、温度目盛りを定める基準状態はただ一つ選べば十分で、その状態の値を定めさえすれば絶対温度目盛りを定めることができる。今日その基準点として水の三重点(水の固相、液相、気相が共存する状態、温度=0.01℃、圧力=611.73Pa=0.006気圧)をえらび、その状態の絶対温度を 273.16 [K] としている。なぜこの値を選んだかと言うと、目盛りの幅がそれまでに普及していた摂氏温度目盛りと一致するようにするためです。そのため絶対温度の1[K]はその1/273.16倍であり、 1[K]=1[℃]となる。そのため絶対温度の0Kは摂氏-273.15℃となる。以後の議論では温度の単位として実用温度には[℃]を、絶対温度には[K]を用いて説明する。
絶対温度目盛りの定義はかなり込み入っているので後で詳しく説明します。そのため、当面は1.や2.で述べたような現実の温度計の示す温度tで持って議論する。温度計構成物質の違いによる温度目盛りの違いは実際はごく僅かであるから当面の議論にはその温度定義で十分です。
温度と同様、熱も当初は訳の解らないものでした。熱は物を燃やしたり、物体を摩擦したりしたときに現れるもので最初は様々な解釈があった。熱はある種の不滅の物質のようなものカロリック(熱素)と考えられたり、またある種の運動のようなものであると考えられたりした。前者の立場では、熱素は物質を構成する原子と結びついたり離れたりして移動する事ができる。そのさい原子と結びつくとき潜んだ形(潜熱)にもなれるし解放されて温度計で検出可能にも成るというのである。また後者の立場では熱とは物体の目に見えない程小さな各部分が極めて激しく振動することであるとするもので、摩擦熱や、我々が感じる熱いと言う感覚を旨く説明した。
このように熱の実体は当初よく解らないものでしたが、燃焼による加熱にしろ、摩擦熱による加熱にしろとにかく何らかの方法で水を熱したとき水1gの温度を1度(1℃)だけ上げるのに必要な熱量Qを1calと定める。この定義は熱量の測定法を定めているだけで熱の実体が何であるかは何も説明していない。しかし、その量を定めることはできるので当面の議論はこれで十分です。ちなみに単位のcalはラテン語calor(熱)に由来する。
物体の融点や沸点に於いて生じる現象から明らかなように、熱量が温度の上昇に用いられるのではなくて、物体の集合状態の変化に用いられる場合がある。温度が変化しないのに熱が吸収されたり放出されたりするため、温度や熱量の理解を非常に難しくするが、熱量の出入り自体は厳密に測定することができる。そのとき物質1gを一つの集合状態から他の集合状態に移すのに必要な熱量の事を潜熱という。
潜熱にはそれぞれの状態変化に応じて、固体→液体の変化に伴う物を融解熱、その逆の変化に伴う物を凝固熱、液体→気体の変化に伴う物を蒸発熱(気化熱)、その逆の変化に伴う物を凝縮熱(液化熱)、固体と気体の間の変化に伴う物を昇華熱と呼ぶ。
集合状態の変化以外に、物体の混合過程、溶解過程や化学変化などにも熱の出入りが伴う。それらを混合熱、溶解熱、反応熱、解離熱等々と言う。この中で反応熱については別項「反応熱と熱化学方程式」で詳しく説明した。
このように温度変化を伴わない熱の移動は熱というものが何であるか考えるのをとても難しくする。しかしながら、そのことによって生じた熱量は比熱の測定に用いたのと同様な装置(例えばヘスHessが反応熱の測定に用いたような)で一定量の水の温度の上昇に置き換えることができ熱量を厳密に測定することができる。そのため潜熱や反応熱も上記の熱量の定義に含ませることができる。
物体の温度を単位温度(1℃)だけ上昇させるのに要する熱量のことをその物体の熱容量Cという。
定義から明らかなように物体の得た熱量Qと、そのために起こった温度上昇t’−t=Δtの間には
の関係が成り立つ。そのため熱容量の単位は[cal/℃]となる。
単位質量(1g)の物質の温度を単位温度(1℃)だけ上昇させるのに要する熱量をその物質の比熱cという。
比熱とは物質1g当たりの熱容量のことであるから、物質の質量をm[g]とすると前記の熱容量Cとは
の関係にある。そのため比熱の単位は[cal/g・℃]となる。
比熱の定義と 1.(2)1.の熱量の定義から直ちに水の比熱は1.0[cal/g・℃]であることが解る。
(1)、(2)式より重要な関係式
が得られる。これは熱と温度を関係づける最も重要な式です。
種々の物質の比熱は、具体的には以下の様にして測定します。
これは高校物理で習う事柄ですが、このようにして比熱を測定できることこそが熱量という量の定義を可能にする。
ここで述べた測定には、摩擦や圧縮・膨張といった機械的・力学的な仕事のやり取りはなく、熱の移動だけで実行できる。そのとき、同じ物体の比熱を、例えば上記の水を他の液体に変えたり、他の物体を仲介にして測定するなどして様々な異なった方法で測定してもいつも同一の値が得られることが経験で知られている。
そのことから熱の移動のみを伴う現象については、いわゆる熱量保存の法則がなりたち、熱というものは何も無いところから生まれたり、あるいは消えてしまったりせずいつも保存される量であることが理解できる。歴史的には、このことが上記の熱量、熱容量、比熱と言った量の有用性を保証し、それらを定義することの根拠を与えてきたのである。
しかしながら、この法則は、熱量について様々な誤った憶測(例えば熱素説のような)を呼び起こしてきた。そのため熱量について次に述べる事柄は特に肝に銘じておく必要がある。
今日熱量のみでは保存則は成り立たないことが解っており、熱は後に説明するエネルギーという量の一形態であるとされている。
一般の熱量測定は1.(2)1.における定義より明らかなように、外部から加えられた熱量、又は外部へ与えた熱量の値を与えるのみで、ある温度の物体中にいかほどの熱量が含まれているかという問いには何も答えることはできない。
つまり、ある標準状態Aにおける熱量を0と定めた後に、その状態から出発して別の状態Bに到達するまでに、その物体に供給された熱量で持って状態Bの熱量を定めてもその熱量は何の意味もないのである。
多くの経験事実から明らかな様に、状態Aと状態B定めてもAからBへ移動する方法・経路によって物体に出入りする熱量は様々に変化する。例えば0℃、1気圧の気体を100℃、10気圧にするのに、まず1気圧に保って100℃まで熱し次に温度を一定に保って10気圧に圧縮する場合と、最初に0℃の下で等温的に10気圧まで圧縮してその後定圧的に100℃まで加熱した場合で加えるべき熱量は全く異なる。
熱のこのような性質を「熱は状態量ではない」という言い方をする。
一般に比熱は温度によって変化する。同じ1gの水の温度を1℃だけ上げるときでも、最初の温度が異なれば必要な熱量は少しではあるが変化する。そのため1.(2)1.熱量の定義も正確には水1gを14.5℃から15.5℃に温めるときの熱量を1calと決めている。
熱力学的状態は温度の他に圧力とか体積の様な他の変数にも依存する。そのためそれらが規定されないと比熱も確定しない。固体や液体の場合の比熱は外圧が一定で行われるか、体積を一定で行われるかの違いによる差はその値の数%程度以下なので特にその条件を特記しないが、気体の場合はかなり異なる。
そのため気体の比熱については、測定を一定圧力のもとで行な場合を定圧比熱cp、一定体積のもとで行う場合を定積比熱cvと呼んで特に区別する。そのとき、後に述べる理想気体については cp−cv=R/M (R:気体定数、M:モル質量=分子量[g/mol])という重要な関係式(マイヤーの関係Mayer's relation)が得られる。[証明は4.(2)2.参照]
比熱の定義は普通物質1g当たりについて定めるが、物質1mol当たりに付いて定義する場合もある。それを特にモル比熱Cといい単位は[cal/mol・℃]となる。モル比熱にすると理想気体について成り立つマイヤーの関係式は Cp−Cv=R (R:気体定数〜2.0cal/mol・℃)となる。[証明は4.(2)2.参照]
なぜ、物質のモル当たり(つまり物質を構成する原子や分子の個数当たり)が重視されるかというと次に述べる興味深い事実があるからです。
気体については、定圧モル比熱Cpや定積モル比熱Cvが気体の種類によらずほぼ一定の値(単原子気体についてはCp〜5.0cal/mol・℃、Cv〜3.0cal/mol・℃)をしめす。気体のモル比熱が気体の種類によらないと言うことは、とても不思議な事で、それは何か一般的な自然法則に基づくことを推測させます。しかし、それをここで論じている熱力学によって説明することはできない。熱力学公理系の守備範囲外の事柄です。
固体の比熱についても、特に原子量の重い元素の固体については、特別の低温でなければ定積モル比熱Cvは元素の種類によらずほぼ一定の値(〜6.0cal/mol・℃)を示す(デュロン=プティの法則 Dulong-Petit law) ことが知られています。この事実も何らかの自然法則が関係していると推察されるが、これも熱力学公理系の領域外の事柄です。
今日、これらの不思議な事実は統計力学や量子力学によって説明さている。
以下で理想気体と言われるものの性質を論じます。これは熱力学の理論を展開する仮想的な物質例としてとても大切です。高校物理で習うので簡単に説明します。
気体の示す圧力を大気圧を例に説明する。大気の圧力はトリチェリー(伊)の指示で弟子のヴィヴィアーニが1643年行った有名な実験から計算できる。
ここで水銀の密度は13.6g/cm3だから、ガラス管の断面積をSm2とすると
となる。
この水銀柱の高さが高度により変化する事を確かめて大気圧にその原因を帰着させたのはパスカル(仏)である。パスカルの依頼で義兄のペリエが様々な高度でトリチェリーの実験を繰り返し、高度との関係を確かめた。
ボイル(英)は尖端を閉じた曲がったガラス管に水銀を注ぎ込み、尖端部に閉じこめられた気体の体積と水銀柱の高さとの関係から、気体に対して成り立つボイルの法則を見つけた。温度を一定に保たねばならない事にボイルは特に気を付けていなかったが、そのことを明示したのはマリオット(仏)で1680年頃ボイルとは独立にこの法則を発見した。
シャルル(仏)は1787年頃にこの法則を発見していたが報告していなかったので世に知られなかった。
詳しく研究したのがゲーリュサック(仏)である。ゲー・リュサックは、水銀上に立てたフラスコDを試験しようとする気体で満たして、装置全体を水槽に浸けた。そして100℃の沸騰水と0℃の氷り水の温度間で体積変化を一定圧力の下で測定した。多くの気体(空気、O2、H2、N2、NO2、NH3、HCl、SO2、CO2、エーテル等)で調べたがいずれも0℃から100℃の間に0℃の体積の0.3749〜0.3752倍膨張(膨張係数=1/266)した。気体の種類によらずほぼ一定の値を示すという驚くべき結論が得られた。
ゲー・リュサック以前に同様な測定を試みた人は多くいたが水上置換の水蒸気圧誤差の為に旨くいかなかった。ゲー・リュサックが成功した理由は、蒸気圧の小さい水銀上でおこない、しかもフラスコに詰める気体をあらかじめ塩化カルシウムを通して十分乾燥したことにある。
その後の精密な測定によると気体の膨張は、気体の種類によらず1℃当たり、0℃の体積V0の1/273倍ずつ変化するということがわかる。実はここのところは非常に微妙なところで、実用温度t(それは当然1.(1)3.に述べた何らかの物質の膨張の様子によって測る)は温度計物質に依存していてどの温度計を用いるかによって結果は微妙に変わります。だから0℃以下、100℃以上でtが1℃変化したとき体積が本当にV01/273倍ずつ変化するのか微妙な所で、この当たりが実用温度tのあやふやなところです。次の議論はその微妙な所は無視した上での話です。
以後はこの値1/273を用いて説明する。
ここで、273+tをTと置き換えると法則を表す式が特別簡単な形になることは特筆に値する。その当たりは膨張係数の意味をグラフの上で考えてみると良く解る。
式の特性から-273℃が特別な意味を持つことが解る。気体温度計を用いる限り-273℃を温度目盛りの零点とした温度(絶対温度)がより正統性を持つことになる。つまり-273℃を0に置き換えると、y=ax+bのグラフがy=axのグラフに成り、シャルル=ゲーリュサックの法則の”体積が温度に比例する”という言い方になる。このときの温度は273+t[℃]=T[K]と置き換えたいわゆる絶対温度Tについて成り立つものであることは特に注意すべきである。ここで導入された絶対温度Tはまことにあやふやなもの[実用温度tがあやふやなので]ですが、当面はそれを絶対温度としておく。
絶対零度や絶対温度の着想自体は意外と古くからある。絶対零度の概念に思い至った最初の人はおそらくフランス人のグイヨム・アモントン(Guillaume
Amontons 1663〜1705年)であろう。彼は1600年代後半の人でボイルやマリオットの晩年と同時代の人です。彼は下図の様な空気温度計を作った。
温度は容積の変化ではなくて圧力の変化で測定される。彼は気圧の効果も考慮した上で、圧力がゼロとなる様な空気の状態のときを最低温度とする考えに至った。それを基準にするとパリの最も寒いときと、暑いときの温度比は5:6であるとした。今日の摂氏温度でいくと約-240℃が最低温度と言うことになる。
空気温度計の開発当初から空気の体積がゼロになる点が-270℃付近にあり、それを絶対零度とする考え方はアンモントン以降もランベルト(728〜1777年)やドールトン(1766〜1844年)により1700年代から1800年代初頭に繰り返し取り上げられている。
ボイルの法則とシャルルの法則は以下の手順で統合できる。
このような統合が可能なのは下図のようにボイルの法則に沿った経路ととシャルルの法則に沿った経路をつなぎ合わせればP-V座標上の任意の点を繋ぐ事ができるからです。
この法則は一定質量の気体を逃がさないようにしておけば、それを温めようが冷やそうが、また圧縮しようが膨張させようが、どんなことをしても良い。そのとき最初の状態のP1,V1,T1と最後の状態のP3,V3,T3を上記の式の形で組み合わせれば等号が成り立つと言っているのである。これは気体の振る舞いに関して大きな制約を課す驚くべき法則です。このとき上記の形でないとだめです。またTは絶対温度でないとだめです。
ボイル=シャルルの法則は気体の種類によらず普遍的に成り立つ法則ですが、厳密な測定(絶対的な温度が定まっていないのでこの言い方には語弊がある)を行うと気体ごとに法則から微少なずれ(一定値からずれる)が存在することが解る。そこで、改めてボイル=シャルルの法則が厳密に成り立つような仮想的な気体を理想気体と言うことにし、その理想気体の膨張係数でもって当面の絶対温度Tを測ることにする。理想気体そのものが仮想的なものだから、その膨張係数で定める絶対温度Tも全く仮想的なものです。今の段階では実用温度tとどの様な関係にあるのかは知ることはできません。
ここの所は、現実には定めようがない絶対温度Tに関係する性質(つまりボイル=シャルルの法則)で理想気体を定義し、その理想気体で絶対温度Tを定義している。つまり世の中に流布している理想気体の定義は循環論法の矛盾に陥っていることを心にとめておいて下さい。
その理想気体に対して成り立つボイル=シャルルの法則にアボガドロの法則(同数の気体は同温度、同圧力で同体積を占める)を適用するとボイル・シャルルの法則の一定値を決めることができる。アボガドロの法則によると1mol(6.02×1023個)の気体は0℃(=273K)、1気圧で22.4リットルを占める事が解っているので
となる。そのときnmolの気体はそのn倍の体積を占めるので
となる。上記の式を気体の状態方程式というが、モル当たりにしてこの様に表したはファント・ホッフのアイディアです。このとき、高校化学で習うように右の二式は様々な気体状物質の分子量測定に利用できる。
気体は(P,V,T,n)or(P,V,T,w)の内の三つを決めれば他は自ずと決まる。普通はn or wを定めた一定質量の気体について考察するので、その場合は二変数を定めれば残りの変数は自ずと定まりP=f(V,T)、V=f’(P,T)やT=f”(P,V)の様な関数として表される。
気体定数Rの値はPやVの単位の取り方により変化する。このときPVの単位が仕事(エネルギー)の単位になることに注意。
このように定めた1mol当たりの気体定数Rは気体の種類によらずどの気体でも同じ値となる。
このHPではモル当たりではなくて単位質量についての状態方程式を用いるが、その場合には上記の値をモル質量M[g/mol]で割っておけばよい。その様にグラム単位で表した気体定数Rは気体の種類ごとに異なった値となる。
[熱力学第一法則]
仕事は熱に、また熱は仕事に変わることができるが、そのとき一方の量は他方の量に常に比例する。つまり両者の間に当量関係がある。
熱力学第一法則は、エネルギー保存則を熱の発生や吸収を伴う現象に拡張適用したものに他ならない。このときあらゆる原理・法則について言えることであるが、法則の正しさを保証するものは、あらゆる場合に絶えず新しくその正しいことが確認されてきた長年にわたる人類の経験した事実しかない。
一つの物体系のエネルギーとは、系がそのとき存在する物理的状態によって定まる一つの量を意味する。ある状態にある物体系のエネルギー値はその物体系のある標準状態に対する値としてのみ決定できる。つまり、ある任意の状態にある物体系のエネルギーとは、その系をその状態から何らかの方法によって標準状態に移すときに、系外に持ち出されるあらゆる作用の力学的当量の総和に等しい。つまりエネルギーとは、系が外部に対して作用をおよぼすことができる能力のことに他ならない。この様にエネルギーがある物体系の内部に存在する部分を意味するとき、これを今後その物体系の内部エネルギーと言うことにする。
上記の定義を完全にするには、系が外部へ及ぼす作用の力学的当量についての考察が必要です。
外部作用として最初に考えられるのは力学的仕事である。これは多くの経験事実から外部系へ働く力と移動距離の積となることが解っている。そのとき
の図からあきらかなように、気体の塊のような均一系においては、その仕事量は圧力×体積変化量で表せる。
次に大切なものは熱的作用である。1.(2)で定義説明した熱量の移動によるものであるが、多くの経験事実から熱量の単位と力学的単位との間には 1[cal]=4.1868[J] の関係があることが解っている。そのため出入りする熱量にこの当量関係を適応して仕事の単位に変換することができる。
これ以外にも必要ならば、電磁気的な作用、化学的な作用等々任意の作用についても当量関係を探し出してこの議論に取り入れればよい。
このような作用の当量関係の正当性は経験事実によってのみ検証される。つまり一つの物体系を様々な道筋を通って標準状態に移行させてみる。そのとき道筋毎に作用の現れ方は異なるであろうが、それらの道筋毎に現れるあらゆる作用を力学的当量に換算して総和を取ってみる。その総和が道筋の違いによらず常に等しくなることをもって検証されるのである。Jouleが行った様々な検証実験はまさにその様なものの一例です。これらの事実の積み重ねの全体がエネルギー保存則といわれるものです。
そのとき、エネルギー保存則は、無からエネルギーを生み出すことはできないことを意味する。そのためエネルギー保存則の別表現として力学的・熱学的・化学的・その他いかなる種類の装置を用いても、それによって絶えず仕事あるいは運動エネルギーを無から生じる機関を構成することは不可能であると言ってもよい。このような機関を第一種永久機関ということにすると、エネルギー保存則は第一種永久機関不可能の原理が成り立つことと等価である。
エネルギー保存則あるいは第一種永久機関不可能の原理を今後熱力学第一法則と呼ぶことにする。
物体系の持つエネルギーは熱のやり取りや摩擦・圧縮などの仕事のやり取り、またそれ以外の諸々の物理的作用によるやり取りによってはかられる。そのとき、物体系の二つの状態が確定されていれば、その二つの状態を推移するとき外界とやり取りされる作用の総和は常に等しいのであるから、物体の持つエネルギーは物体の状態が定まれば一意に定まる量だといえる。
これこそランフォード、ゲーリュサック、マイヤー、ヘルムホルツ、ジュール等の人々によって多くの経験事実から明らかにされたエネルギーという量のもっとも重要な性質で「エネルギーは状態量である」と言われる所である。
3.(2)2.で説明するように熱d’Qも仕事d’Wも完全微分ではない。それにもかかわらず、その和d’Q+d’W=dUは完全微分であると言っているのが熱力学第一法則の内容です。(完全微分の意味は5.(2)で説明)
熱力学理論で重要な状態変化は、無限にゆっくりと行われ、従って常に平衡状態にあると言われる過程で行われる。字義通りにとると、この言葉づかいは実に曖昧である。というのは、過程は必然的に変化であり、従って平衡の乱れを前提としているからである。しかしながら、変化の速さを問題とせず唯その結果のみを問題とするならば、この乱れは欲するままに幾らでも小さく採ることができる。
例えば、各瞬間において外圧を気体の圧力よりも極僅か大きくして、気体をいくらでも緩慢にその体積を任意の程度に圧縮することができる。またそのとき、一定の有限の圧縮に用いられた仕事を計算する場合に、外圧の大きさを取り扱う代わりに気体の圧力を採用しても非常に僅かの誤差を生じるに過ぎない。極限に持って行けば此の微少誤差も消失し、無限にゆっくりした圧縮に於いては、このようにして得られた結果が厳密に正しい。
ここで述べたことは、一定の圧力での圧縮・膨張にも圧力が変化する場合の圧縮・膨張についても同様に行われる。後の場合においては、小さな錘を添加したり除去したりして、外圧を常に所用の大きさに等しくしなければならない。錘を添加・除去するのに水平に押す等の無限小の力で作業できる様な配置を工夫して仕事を要しないで圧力の変化を生じせしめる様な装置を構築して作業するものとする。
熱の授受についても、外部からの仕事の場合と全く同様である。系が外部から受けるか外部へ与える値のみを取扱い、その際時間が何らの役目を勤めないならば、系が熱を受け取るか又は与えるかに従って、熱源の温度を系の温度より極めて僅かだけ高くするか、又は低くすれば十分である。この微少の差は過程の向きを決定するが過程によって最後に系に誘起される熱量の移動量には関係しない。その様にして気体の圧縮と同様に、同温度にある一物体から他物体への熱の移動を論じることができる。
これは等温過程のみならず温度の変化する過程も含まれている。そのときには、もちろんただ一つの一定温度の熱源だけでは間に合わない。例えば適当な圧縮や膨張によって自在に温めたり冷やしたりできる気体の様なものを用いたり、あるいは様々な温度を有する多数の熱源を準備して、各瞬間に於いてその瞬間の系の温度にできるだけ近い温度を有する熱源と次々に接触させる等のメカニズムを利用するものとする。
以上の様な変化を準静的過程と言う。この考え方の重要なことは、無限に緩慢な過程はすべて反対の向きにも進ませることができる点にある。これらの過程は無限小の偏差を無視すると平衡状態にあると言えるので、その過程を逆向きに進ませるためにはやはり無限小の変化を与えれば足りる。そのため準静的過程は必然的に可逆過程である。
物体系のもつ状態というのは具体的にどの様に表現されるのであろうか。熱力学に於ける均質系の状態は、一般に三つの物理量、即ち圧力P、体積V、温度Tによって定まり、これらの量は一般に f(V,P,T)=0 と言う形の状態方程式を満たす。均質とはその系を構成する物体の任意の微少部分に於いて完全に同質の場合を言う。例えば理想気体の場合、2.(5)で説明したように状態方程式はPV−RT=0 となる。
このとき、均質な系に於ける温度T、圧力P、体積V等は状態が定まれば一意に定まる状態量である。このことはそれらの量の性質を考えれば明らかである。
物体の状態を表すこの三つの量の間に一つの関係式(状態方程式)が存在するので、V-T座標平面に垂直な座標軸としてPを選んだ場合、状態方程式 f(V,P,T)=0 は P=f’(V,T)なる一つの曲面をあらわすことになる。その曲面上の点でもって物体系の状態を表すことができる。そして、均質系の状態変化は曲面上の点の連なり、すなわち曲線でもって表される。そして曲線の端点が最初の状態と最後の状態を表す。
実際には三つの内の二つ、例えばTとV、を決めれば残りのPは状態方程式により必然的に定まるので、三変数の内の二つで状態変化を表すことができる。状態が変化すれば、状態を決める点(V,T)はV-T平面上で一つの曲線を描く。それらの端点を1と2として、その端点1が最初の状態、端点2が最後の状態を表す。また曲線途中の点の連なりが状態変化を表す。その様な二変数座標平面上の曲線図形を状態変化の経路図(ダイアグラム)という。
そのとき、物体系が最終的に初めの状態に戻るときには経路図(ダイアグラム)は閉じた曲線になり、循環過程(サイクル)と言われる。
この当たりの例としては別稿「気体の断熱変化」の理想気体の場合の図などを参照されたし。
状態を表す変数として(V,T)の組み合わせ以外に(T,P)や(P,V)の組を選んでもよい。そのとき P=f(V,T)曲面上の経路図をT-P平面、P-V平面へ射影すれば、それぞれの平面上で一つの曲線が得られ、いずれも今考えている状態変化を表す経路図(ダイアグラム)となる。そのとき、状態変化を説明するのにどのダイアグラムを用いてもよい。
さらに内部エネルギーUや後で説明するエントロピーSも状態量なので(U,P)、(U,V)、(U,T)や(S,P)、(S,V)や(S,T)などの組み合わせの経路図を描くことも可能です。中でも縦軸をT、横軸をSにしたT−S経路図は極めて重要です。それは二つのS=S1とS=S1+dSの直線と経路曲線とS軸に囲まれた領域がちょうどTdS=d’Qとなり、P−V経路図の仕事PdV=d’Wと同じように、状態変化に伴って系が外界とやり取りする熱量を意味するからです。また熱機関の効率はまさに縦軸Tの作動温度に関係するからです。この当たりは後ほど6.で詳しく説明する。
物質系の状態1に於けるエネルギーをU1、状態2におけるエネルギーをU2とすると、熱力学第一法則はU2−U1が状態1から状態2へ推移するとき外界から物体系へ成された作用の総和に等しいことを言っている。つまり、系外から系内へ移動した熱量をQ'、系外から系へ加えられた仕事量をW'とすると
となる。[熱量や仕事量は、次で説明するように状態量ではないので、その事を明示するために’(ダッシュ)をつけます。]
前記の経路図の考え方を質量M以外に圧力Pと体積Vとによって決定される均質系に適用する。例として下図の様な状態1と状態2の内部エネルギー変化を考える。
すでに述べたように内部エネルギーUが圧力Pや体積Vの変化によりどの様に変化するかを調べるには、その系を何らかの方法で最終状態(圧力や体積が異なった)に移行させ、その移行に伴う外部作用(熱の出入りや仕事のやり取り)を測定すればよい。その時、状態1から状態2へ移行するとき系に対してなされる仕事は
と表されるが、その積分値は経路が異なると全く異なった値(上図の灰色部分の面積×(-1)、経路αと経路βでの違いに注意)になる。このことから明らかなように
仕事は状態量ではない。経路毎に全く異なった値となる。これが仕事量の微少変分d’W=PdVが完全微分でないと言われることの意味です。
このときU1とU2は状態量だから上式の左辺U2−U1は状態1から状態2へ移行する経路にかかわらず同じ値であるが、仕事W’が経路により変化するため、移行時に出入りする熱量Q’も経路により異なることになる。つまり熱力学第一法則より
の関係があるのだから
熱は状態量ではない。経路毎に全く異なった値となる。これは1.(4)2.ですでに注意しておいた事です。これが熱量の微少変分d’Qが完全微分でないと言われることの意味です。
仕事と熱が状態量ではないということは、どんなに強調しても強調しすぎることが無いくらいに重要な事柄です。熱力学が良く理解できないのはひとえにこの事実を忘れてしまうからです。
物質系がある状態から任意の変化を経て元の状態に戻るような循環過程(サイクル)をたどる場合には、上記の熱力学第一法則は
となる。つまり循環過程においては外部から系に及ぼされる作用の仕事当量の総和は零に等しい。下図のP-V図の場合循環経路で囲まれた部分の面積が1サイクルの間に系に対して成された仕事であり、それを打ち消すだけの熱が外部へ流出した事になる。
ここで、物質変化を特徴付ける幾つかの係数を定義しておきます。
定圧膨張係数(体積膨張率):α
定圧の条件下で1℃だけ温める際の体積増加を0℃における体積V0で割ったものです。すなわち
のことです。理想気体の場合は
となります。
定積圧力係数:αP
定積の条件下で1℃だけ温める際の圧力増加を0℃における圧力P0で割ったものです。すなわち
のことです。理想気体の場合は
となります。
等温体積弾性率:κ
定温の条件下で圧力の無限小増加によって生じる単位体積に生じる体積減少量の比です。すなわち
のことです。理想気体の場合には
となります。これ以外に断熱体積弾性率というのもあります。その場合はγPとなること(γ=cp/cv)は別項「音速の理論」で説明しました。
等温圧縮率:κ’
弾性係数の逆数、すなわち単位体積についての無限に小なる圧縮の、それに対応する圧力の増加に対する比です。すなわち
のことです。理想気体の場合には
となります。これ以外に断熱圧縮率というのもあります。
ここで述べた定圧、定積、等温の各変化に於ける物質の変化性状を示す三個の係数は互いに独立なものではなくて、任意の物質について一つの定まった関係式で結び付けられる。すでに述べたように、均質な物体系には、三つの状態量P、V、Tに対して状態方程式P=f(V,T)が存在する。f(V,T)を今後P(V,T)と書く事にすると、dPは完全微分なので
となる[完全微分の意味については5.(2)参照]。dVやdTを利用しても同様な関係式が得られる。
このことの意味を理想気体を例にして説明すると、下図のΔV,ΔT,ΔPについて
が成り立つ事である。これらの関係式は単に状態方程式曲面に対して成り立つ多変数関数の数学的な事実を表しているに過ぎない。
この関係式を用いると未知の係数を残りの二つから計算する事ができる。例えば
となる。すなわち0℃から1℃温める際に水銀の体積を不変に保つためには46atmの圧力増加を必要とする。
熱力学第一法則から多くの簡単ではあるが重要な結論が導かれる。ただし、絶対温度の厳密な導入・定義をまだしていないので、ここでの温度は1.(1)で説明した実用温度tか、又は単にt+273=Tと置いた温度だと思って下さい。このTは、2.(3)and(5)で注意したように、かなりいい加減な意味での絶対温度です。
ある断熱的な容器に詰めた一様な温度にある気体に真空の容器を接続して、その接続部の活栓を開けて気体を真空中へ自由膨張させる。このとき気体は外部に対して何らの力学的な作用を及ぼさないから仕事のやり取りはしない。なぜなら真空中に噴出した気体がその容器壁へ衝突したとしても、その硬い壁は動かないのでその衝突による作用は外部に何らの力学的な仕事をしないからである。また断熱的に行っているので熱のやり取りもない。そのため熱力学第一法則(エネルギー保存則)より気体が持つ全エネルギーは最初の状態と二つの容器に拡散した最後の状態では変化しないと考えなければならない。
このような実験はゲーリュサックやジュールなどにより繰り返し行われた。そのとき彼らは気体の真空中への膨張で気体の温度がどうなるか調べたのである。具体的な手順としては高圧の気体を詰めた容器と真空の容器を活栓で接続して水槽に浸けて、熱平衡が実現されたときの水槽の水温をまず測定しておく。その後活栓を開いて暫く経ったとき水槽の温度が変化するかどうかを確かめたのです。
彼らはそのとき気体の体積が変化(当然圧力も変わる)したにもかかわらず温度が少しも変化しないことを見つけた。厳密に言うと極僅か変化する場合もあるが、彼らは理想的な気体については変化しないとした。つまり上記の実験で温度が変化しないような仮想的な気体を理想気体と言うことにしたのである。
この経験的な事実から「理想気体の内部エネルギーは温度のみに関係し体積や圧力には無関係である。」という重要な結論が得られる。これは状態方程式からは決して得られない、理想気体について一歩踏み込んだ性質です。状態方程式はPとVとTが互いにどの様に関係して変化するかを示しているだけで、状態が決まったとき(つまりP、V、Tが決まったとき)の内部エネルギーがどういった値になるのかは何も教えてくれない。その何も解らない内部エネルギーについて一歩踏み込んだ性質を教えてくれています。ただし内部エネルギーが温度にどの様に関係するのかは何も教えてくれない。
高校物理の教科書ではこの結論がどこから出たのか曖昧で初学者は面食らうところですが、これは単なる経験事実(ゲーリュサック=ジュールの実験や次に述べるジュール=トムソンの実験)を理想化したものでしかない。熱力学ではそれ以外の根拠は無いことに注意されたし。
前項で述べたゲーリュサックやジュールの実験は気体の熱容量が器壁や水槽のそれと比較して非常に小さいので精度的には十分とは言えない。そのため精度を改良するために行われたのがW・トムソンが工夫しジュールの協力のもとで行われた、有名なジュール=トムソンの実験です。
この実験の改良点は、気体の流出を緩慢にすることにより直接に第二の平衡状態に移し第二の平衡状態の温度T2を気体の中で直接測定することです。そのため、気体の熱容量が小さいことは不利な事ではなくむしろ有利に働く。
この実験では、気体が高圧P1の場所から低圧P2の場所(大気圧)へ、途中の柘植(つげ)材の管(その中に綿花もしくはほぐした絹がつめられており、気体が非常に通過しにくくしてある)を通って定常的にゆっくりと移動させられる。
実際の測定は下図に示す装置を用いて、次のように行われる。一定温度の水槽中に置かれた長い螺旋パイプの末端部に綿、絹、その他の繊維をほぐしたものを詰めて細孔栓とする。パイプは細孔栓の近くのみ断熱されており、それ以前は熱伝導性の良い金属銅パイプで構成してあり細孔栓に流入する空気の温度が厳密に水槽の水温に一致するようにしてある。細孔栓の前後の断熱された部分の管内に温度計を挿入して定常的に流れる気体の細孔栓直後の温度を測定する。測定部より下流側は大気圧(1気圧)になるように調整されいる。また上流側にポンプが設置されてあり、大気圧よりも少し高い圧力で定常的に気体を送り込むようにしてある。温度測定はポンプを連続的に働かせて定常状態になったときに行われる。
測定の結果は、空気については非常に小さい温度変化があり、水素についてはそれが更に一層小さくほとんど測定にかからないほどであった。それゆえ理想気体に対しては温度変化が無いとして良いだろう。その様な気体を仮想的に理想気体と言うことにしたのである。
この測定は断熱的におこなわれ、しかも温度変化が少しも起こらないと見なして良いのだから理想気体に於けるボイルの法則を用いれば、細孔栓の手前での体積をV1、栓を通過した後の体積をV2とすると、P1V1=P2V2が成り立つ。
そのため、高圧側で外界から気体に対してなされた仕事P1V1は低圧側で気体が外界(大気)に対してなした仕事P2V2に常に等しくなる。それは気体に対して正味の仕事はなされず、気体の内部エネルギーが変化しないことを意味する。
これはかなり込み入った考え方です。つまり、普通の気体ならば温度変化が起こらないからと言ってP1V1=P2V2の保証などはないのだが、理想気体というものは温度が一定ならボイルの法則が成り立つようなものとしていたのだから理想気体というものを仮想的に考えるときにはP1V1=P2V2が成り立って仕事は零で、そのため内部エネルギーは変化していないと言っても良いだろうということです。
つまり、気体の体積は変わるのに、その温度は大体一定に保たれているのであるから理想気体の内部エネルギーは温度のみに関係して体積には無関係であると言える。すなわち
となる。
普通この実験はエンタルピーとの絡みで説明されるが、もともと理想気体というものの性質を見極めるために行われた実験であることを忘れないでほしい。この実験の一番の意義はそこにあるのだから。
理想気体ではない現実の気体に対しては温度変化を生じる。これは断熱的な変化であるにも関わらず内部エネルギーが変化する事を意味し、外部から気体になされる仕事と気体が外部に対してする仕事の間に差がある事を意味する。
例えばΔP=P1−P2=1気圧当たりで説明すると、T1=0℃のとき、空気に対してはΔT=-0.3℃、二酸化炭素ではΔT=-1.4℃であり、水素では逆に温度が上昇してΔT=+0.033℃となる。水素でも十分低温になるとΔTは負になる。温度を高くし、密度を小さくするとΔT→0となり理想気体に近づく。
このように実在気体に於いて温度が少し変化する事をジュール・トムソン効果という。これは気体分子の間に引力が働き、分子が大きさを持っているために生じるものであるが低温になるほど大きくなる。この効果は気体の液化や極低温をつくる装置に利用される。また実在気体を用いて絶対温度を定めるのに利用される[7.(2)参照]。
多くの教科書では、理想気体の性質として
が成り立つとしています。
そして、1.が成り立てば2.や3.が証明でき、逆に2.と3.が成り立てば1.が証明できるとしている(例えば戸田盛和著「いまさら熱力学?」丸善(1997年刊)p73〜82など)。
しかし、この当たりの説明は不十分なのでもう少し補足しておきます。2.(5)で述べたようにPV=RTを用いる理想気体の定義は絶対温度Tが厳密に定義されていないので循環論法の論理矛盾があり、もともとこの式は意味をなしていないのです。
だから正しくは”2.や3.は実験事実を理想化してその様な性質を持つ気体を理想気体と言うことにする”です。そのときの2.や3.で言う温度は1(1)で述べた実用温度tでかまいません。それでは2.や3.の事実と1.のPV=RTとの関係はどうなるのかと言うと、”状態方程式中のTを、熱力学第二法則を用いて定めた絶対温度(7.(1)、(2)で定めます)であると解釈するとPV=RTが2.や3.と論理矛盾が生じていないことが言える”、です。決して、”1.から2.や3.が導ける”、と言われるような関係ではありません。プランクの本を読むとこの当たりの事情がよく解る。詳しくは7.(1)2.をご覧下さい。上記の戸田盛和氏の本の説明も論理矛盾が無いことを確認しているのだと解釈されれば良いと思います。
ここでは物質量1モル(n=1mol)の均一な一成分系を考える。
熱力学第一法則より
であるが、物質系の内部エネルギーUは状態量で物質系の温度Tと体積Vの連続的な関数U(T,V)と考えることができるので
が成り立つ。そのため
となる。
定積の条件下で加熱される場合、仕事は無いので、系のエネルギー変化は流れ込んだ熱に等しい。そのため定積モル比熱Cvは
となる。
一方、定圧条件下で加熱される場合、定圧モル比熱Cpは
となる。上記の二式を一緒にすると
となる。この式は任意の物質について成り立つもので、右辺は定圧過程で体積膨張により消費されるエネルギーを補うために必要な余分の熱量を表す。
これは別稿「気体のモル比熱」で解りやすく説明しているので、ここでは簡単に記す。4.(1)で述べたように理想気体の内部エネルギーUは温度だけの関数です。1molの理想気体についてはPV=RTが成り立つのでモル比熱に対して
となる。これが1.(5)2.or3.で説明したマイヤーの関係(Mayer's relation)です。
理想気体では
において、Uは温度のみに関係し体積には無関係であるから、Cvは温度のみに関係し体積に無関係として良いであろう。またマイヤーの関係式を援用するとCpも同様に温度のみに関係し体積に無関係として良いであろう。
この結論はH.V.Regnault(アンリ・ヴィクトル・ルニョー1810-1878)が多くの気体の測定により確認した。そのとき彼はさらに、Cpは温度のかなり広い範囲にわたって一定である事も発見した。したがって再びマイヤーの関係式を適用するとCvもまた、その範囲内で一定であるとして良いであろう。
それ故に、Cp、Cvが温度および体積に無関係であるとする事を、理想気体の条件に加えても良いであろうし、また今後そうすることにする。
ここで次のことに注意してほしい。理想気体の比熱が温度、体積によらず一定値であることは、今まで述べた理想気体の性質から理論的に証明できるような事柄ではない。実在の気体の性状から推し量って、その様に置くことにしたのである。理想気体の内部エネルギーは温度のみに関係して体積には無関係であるというのが、実験結果を理想化したものであったように、この比熱についての性質も実験事実の理想化から、そのようにしたのだということを銘記すべし。世の一般の教科書ではその当たりが実に曖昧にしか記載されていない。
これは別稿「気体の断熱変化」で解りやすく説明したので、ここでは簡単に記す。
当然のことですが可逆過程の断熱変化とする。つまり外界と可逆的に仕事をやり取りしながら膨張・収縮する。断熱変化ではd’Q=0だから
音は圧力変化が伝播する現象ですが、変化速度が速いため音の伝播に伴う圧縮・膨張は断熱的に行われると考えて良い。さらに理想気体の状態方程式が成り立つとして良い。
力学理論によると音速は体積弾性率κ、気体密度ρと
の関係にあることが解っている[別項「音速の理論」参照]。この式に熱力学の成果を代入すると
となる。これは測定結果に良く一致する。
不均一系とは明確な境界面で分かたれて均質な物質が幾つか共存して成り立ってるような系や様々な化学物質が混在して化学平衡が成り立っている様な系を意味する。前者では液体がその蒸気と接触している場合や、固体がその融解した液体と共存している場合などが考えられる。また後者では一般の化学反応がすべて含まれる。さらに、固体がその溶媒に解ける溶解現象や、二種類の気体や液体が混合する現象などが含まれる。
これらの現象の多くは後に述べる熱力学第二法則が関係するので8.まで詳細の議論は置いておく。しかしながら、熱力学第一法則のみで説明できる事柄もある。中でも大切な応用は化学反応における反応熱と物質の性質についての議論です。化学反応熱の分野はエネルギー保存則(熱力学第一法則)だけからでも様々な重要な情報を得ることができる。ただし、その当たりは別稿「反応熱と熱化学方程式」で説明したのでここでは議論しない。
6.で説明する熱力学第二法則は極めて抽象的・数学的な法則です。使われている数学概念を理解しないと真の理解は得られないので、まずその準備から始めます。熱力学では偏微分を含む数式が沢山出てきますが、次に述べる数学概念が理解できていないと、その意味を理解するのは難しい。
直交座標空間の中に、定まった方向を持った曲線sが与えられているとする。(x1,y1)を曲線のはじめの点、(xn,yn)を終わりの点とし、曲線sの弧の長さをはじめの点から測るものとする。この空間の各点に連続な関数X(x,y)とY(x,y)が定義されているとする。
微少部分dsのx座標方向の成分をdx、y座標方向の成分をdyとして下記の様な(X(x,y),Y(x,y))との積を作り、その和を取る。曲線の分割の数nを限りなく大きくし、曲線要素dsを限りなく小さくした極限の値を関数(X(x,y),Y(x,y))の線積分という。
このときX(x,y)とY(x,y)を点(x,y)で定義されるベクトルのx成分、y成分と考えるとXdx+Ydyはベクトル(X,Y)とベクトル(dx,dy)の内積を意味するので、線積分とは下図のベクトル(X(x,y),Y(x,y))を経路の微少線分方向へ射影した成分にdsを乗じて経路長に沿って加え合わせて行くことを意味する。
微分方程式
を考える。
(1)の左辺がx,yのある関数の完全微分に等しいとき、すなわち
を満足する関数φ(x,y)が存在するとき、つまり
が成り立つとき、(1)は完全微分方程式という。
このような関数φ(x,y)に対して
は(1)の一般解である。なんとなれば、yをxの関数と考えて(4)の両辺xで微分すれば
となるからである。
次に、(1)が完全微分方程式であるための条件を求める。X(x,y)、Y(x,y)はx,yについて連続な偏導関数を持つと仮定すれば、完全微分方程式ならば(3)式が成り立つので
を得る。(5)は(1)が完全微分方程式であるための必要条件である。
次に、(5)は(1)が完全微分方程式であるための十分条件でもあることを証明する。
とおけば、
となる。よって
となるが、この式はXdx+Ydyが関数
の完全微分に等しいことを示している。よって(1)は完全微分方程式である。
以上をまとめると
が完全微分方程式[つまり左辺がある関数φ(x,y)の完全微分dφとなる]であるためには
が成立する事が必要十分条件である。
そのとき(1)の解は
で与えられる。
一般論だけでは何のことか解らないので幾つか例を説明する。
[例1]
X(x,y)がxだけの関数、Y(x,y)がyだけの関数である場合、常に条件式(5)が満足されるので
完全微分方程式となる。それゆえに一般解は
となる。
[例2]
は条件式(5)を満足するので完全微分方程式である。
よって一般解は
であるから、
で与えられる。この結論は下記のくくり直し法でも直ちに導ける。
[例3]
も条件式(5)を満足する。
よって完全微分方程式である。ゆえに
前節[例2]を利用して完全微分方程式の意味を説明する。完全微分方程式の
と言う条件は、一つのx,yの連続関数φ(x,y)=x2y+3x2−yが存在して
が成り立つから言えるのである。
この例の場合X(x,y)、Y(x,y)、φ(x,y)の値をxy平面に垂直にとれば、それぞれの関数は下図の曲面を表す。つまり(5)式が成り立てば関数曲面X(x,y)とY(x,y)を利用して常に一つの関数曲面φ(x,y)を定めることができる。このことは、下図に於いてそれぞれの関数の偏微分が関数曲面の傾きを表していることを読み取れば了解できる。
つまり完全微分方程式を解くと言うことは、ある座標点(x0,y0)のφ(x0,y0)を定めてからdxとdyを少しずつ変化させて座標点(x0+dx,y0+dy)へ移動するときの変分dφ(完全微分)を求め、それをφ(x0,y0)に積み上げながら関数曲面φ(x,y)を決めていくことです。
このとき φ(x,y)=C (ある定数)は 関数曲面φ(x,y) の等高線を意味する。
例えば x2y+3x2−y=100、x2y+3x2−y=0、x2y+3x2−y=−100のグラフは下図の様になる。
つまり完全微分方程式の (2xy+6x)dx+(x2−1)dy=0 の意味はdφが零になる方向(dx,dy)が存在し、その方向の軌跡を定めることです。解を求める操作とは等高線を少しずつ伸ばして等高線を表す曲線を求めることです。例えば上図の(x,y)=(7,7)の位置では
となるので、その点に於ける等高線のxy平面に於ける傾きは-2.92(下図参照)となる。
ここで、(X(x,y),Y(x,y))=(2xy+6x,x2-1)と(dx,dy)をそれぞれ2次元平面のベクトル成分だと考えるとベクトル(X,Y)とベクトル(dx,dy)の内積が零と言うことは二つのベクトル場が直交していることを意味する。つまり(dx,dy)ベクトルの方向を連ねたものが解φ(x,y)=C (ある定数)であるが、それはちょうど電場ベクトルに直交する等電位線のようなもので、どの(x,y)点に於いても上図の等高線ベクトル(dx,dy)は下図のベクトル場(X,Y)に直交していることが解る。
等高線の集合が関数曲面φ(x,y) に他ならないのだから、xy平面上で座標点(x,y)をあらゆる方向に移動させたときの変分 (2xy+6x)dx+(x2−1)dy=dφ を加え合わせながらプロットしていくと一つの曲面φ(x,y)ができる。そのとき、xy平面上の曲線に沿ってdφ加え合わせてゆくことは、5.(1)で説明した線積分を実施することに他ならない。
このとき最初の点と最後の点が同じならどの経路を通って線積分を実施してもそれらの値はすべて同じになる。なぜなら積分値はまさに点1と点2の標高差に他ならないのだから、経路によらないことは明らかです。
これが、完全微分方程式の左辺がある関数φ(x,y)の完全微分であると言うことの意味です。その完全微分成分dφをある経路に沿って線積分することは関数曲面φ(x,y)上を移動することを意味する。熱力学的な言い方をすれば、φ(x,y)は変数xとyの状態量であると言うことです。
これをベクトル場(X(x,y),Y(x,y))=(2xy+6x,x2-1)に付いての線積分と考えると、Xdx+Ydyはベクトル(X,Y)とベクトル(x,y)の内積を意味するので、下図のベクトル(X(x,y),Y(x,y))を経路の微少線分方向へ射影した成分を経路長に沿って加え合わせて行くことを意味する。
この結果が上図の三経路ですべて等しいのである。これはちょうど、保存力の場に於ける力のベクトル(ポテンシャルの勾配)と仕事の積算値との関係と同じです。
一般の微分方程式では
つまり完全微分方程式ではない。しかしながら、この両辺にx,yのある関数 λ(x,y) を乗じて得られる微分方程式
は完全微分方程式になることがある。このような関数 λ(x,y) を(1)の積分因子という。
例えば微分方程式
は、
であるから、完全微分方程式ではない。しかし両辺に 2x をかけて得られる微分方程式
は
であるから、完全微分方程式である。だからλ(x,y)=2xは積分因子である。このとき当然のことですが積分因子はX(x,y)=1+y2、Y(x,y)=xyの関数形に依存する。
一般に積分因子となるものは沢山あります。例えばCを任意定数とするとλ(x,y)=C・xやλ(x,y)=C・(1+y2)-1/2、λ(x,y)=x・(1+y2)-1/2、・・・・・等々のいずれも微分方程式 (1+y2)dx+xydy=0 の積分因子となります。積分因子であることは条件式に代入してみれば簡単に確かめることができます。
このとき(1+y2)dx+xydy=0が完全微分方程式ではないということは、(1+y2)dx+xydy=d’ψなる量をxy平面上の曲線に沿って計算しながらd’ψを積み上げていったとき、つまり線積分
を実施したとき、最初の点(x1,y1)と最後の点(x2,y2)は同じでも、その線積分値は経路が異なれば全く違った値となることを意味する。熱力学的な言い方をすれば、”(1+y2)dx+xydy=d’ψは(dx,dy)が決まったからと言って一意には決まらず、状態量ではない”、と言うことです。d’ψがある関数ψの完全微分ではないことを示すためにdではなく’(ダッシュ)を付けたd’で表していることに注意。今後完全微分では無い場合’(ダッシュ)を付けることにする。本によってはδψのようにδを用いて表す場合もある。
(1+y2)dx+xydy=0に於いてX(x,y)=1+y2、Y(x,y)=xyのグラフは下図のようになる。
このグラフの中の領域(2≦x≦3、2≦y≦3)に着目して、(x1,y1)=(2,2)から(x2,y2)=(3,3)までを二つの経路 (2,2)→(3,2)→(3,3) と (2,2)→(2,3)→(3,3) で移動する。そのときdx=1、dy=1を単位として移動して、各経路について (1+y2)dx+xydy=d’ψの積算値を計算してみる。
経路をもう少し拡大して (2,2)→(3,2)→(4,2)→(4,3)→(4,4) と (2,2)→(2,3)→(2,4)→(3,4)→(4,4) で積算値を比較すると
となる。図から明らかなようにΣd’ψは積算経路によって全く異なることが解る。
ところが積分因子2xを乗じた(2x+2xy2)dx+2x2ydyについては、この微少成分はある連続関数φ(x,y)の完全微分dφとして表せ、xy平面上の線積分は最初と最後の点が同じなら積分経路に関係なく同一の値になることが言える。熱力学的な言い方をすれば”(2x+2xy2)dx+2x2ydy=dφは(dx,dy)が決まれば一意に決まる状態量である”。
積分因子=2xを掛けて完全微分方程式にした(2x+2xy2)dx+2x2ydy=0に付いて同様な考察を行う。まずX(x,y)=2x+2xy2、Y(x,y)=2x2yのグラフは下図のようになる。
このグラフの中の領域(2≦x≦3、2≦y≦3)に着目して、(x1,y1)=(2,2)から(x2,y2)=(3,3)までを二つの経路 (2,2)→(3,2)→(3,3) と (2,2)→(2,3)→(3,3) で移動する。そのときdx=1、dy=1を単位として移動して、各経路について(2x+2xy2)dx+2x2ydy=dφの積算値を計算してみる。
経路をもう少し拡大して (2,2)→(3,2)→(4,2)→(4,3)→(4,4) と (2,2)→(2,3)→(2,4)→(3,4)→(4,4) で積算値を比較すると
となる。
dxもdyも微少量ではなく1にしているので、有限な大きさであることによる誤差のため二つの経路の積算値は完全には一致していない。しかし、移動単位をもっと小さくしてゆけば極限に於いてΣdφが厳密に一致するであろうことは予想できる。この一致が完全微分方程式の必要十分条件、つまり任意の(x,y)点に於いて
を満たすという意味です。
ちなみにくくり直し法を用いるとφ(x,y)=x2+x2y2であることは直ちに導け、そのグラフは下図のようになる。ただし縦軸のスケールは縮小してある。
このとき等高線のグラフは下図の様になる。
ここで、完全微分方程式ではない任意の微分方程式を完全微分方程式にする積分因子は常に存在するのかが問題となるが、二変数の場合必ず存在することが数学的に証明される。つまり
X(x,y)とY(x,y)がxとyに関して連続な関数であれば、任意の二変数の微分方程式X(x,y)dx+Y(x,y)dy=0に対して、この両辺にx,yのある関数 λ(x,y) を乗じて得られる微分方程式
が完全微分方程式になるような関数 λ(x,y)[積分因子]は必ず存在する。
ここで”二変数の”という但し書きは重要です。三変数の場合には常に存在するとは限らないからです。三変数の場合の積分因子の存在性については参考文献(犬井、スミルノフ等)を参照して下さい。
[証明の概略]
二変数の完全微分方程式ではない微分方程式X(x,y)dx+Y(x,y)dy=0のX(x,y)とY(x,y)を二次元ベクトルの(x,y)成分と考えると、この微分方程式の解はベクトル場(X(x,y),Y(x,y))に直交する曲線群を求めることに相当する。なぜなら微分方程式の意味はベクトル(X,Y)とベクトル(dx,dy)の内積が零であると言っており、これは二つのベクトル場(X(x,y),Y(x,y))と(dx,dy)が直交していることを意味するのだから。
例として、前項で取り上げた完全微分方程式ではない(1+y2)dx+xydy=0について説明すると、ベクトル場(X(x,y),Y(x,y))=(1+y2,xy)は下図の様(長さは1/50倍にしてある)になる。
このようなベクトル場に対して常に垂直に交わる曲線群は必ず存在する。これは電場に対する等電位線の関係に似ています。実際のところ存在することを数学的に証明すること[参考文献スミルノフ(第3巻§75と§51)等を参照してください]は結構面倒です。しかし、存在するであろうことは了解してもらえると思いますので、その関数をF(x,y)=C(Cは任意定数)とする。数学では、これを微分方程式X(x,y)dx+Y(x,y)dy=0の積分と言う。この関数の両辺をxで微分すると、右辺は定数だから零となり、これは
を満足する。ここでX(x,y)dx+Y(x,y)dy=0の要請から
が成り立たねばならない。この相等しい両辺の商を共通の関数λ(x,y)と置けば
となるが、これはまさしくλ(x,y)X(x,y)dx+λ(x,y)Y(x,y)dy=0が完全微分方程式であることを意味している。故に積分因子λ(x,y)は必ず存在して、それは
となる。
[証明終わり]
このとき次のことに注意されたし。このF(x,y)はλ(x,y)X(x,y)dx+λ(x,y)Y(x,y)dy=dφ(x,y)がある関数φ(x,y)の完全微分になるというときの関数φ(x,y)そのものです。つまり、前項の例で求めたφ(x,y)=x2+x2y2=F(x,y) です。なぜなら、積分因子λ(x,y)=2xを乗じた(2x+2xy2)dx+2x2ydy=0についてのベクトル場(X(x,y),Y(x,y))=(2x+2xy2,2x2y)のグラフは
の様になる(長さは1/500倍にしてある)のだが、この図を前出のベクトル場(λ(x,y)X(x,y),λ(x,y)Y(x,y))=(1+y2,xy)のグラフと比較してみと、両図の同じ(x,y)点のでベクトルは、向きが逆になる所もあり、長さの分布の様子も異なるが、その傾きは全く同一になる。
これはベクトル(λ(x,y)X(x,y),λ(x,y)Y(x,y))はベクトル(X(x,y),Y(x,y))を単にλ(x,y)倍したものにすぎないのだから当然です。ベクトルの傾きが全く同じなのだから、それらのベクトル場に垂直な曲線群が等しくなるのは当然です。
実際、積分 F(x,y)=x2+x2y2 と X(x,y)=1+y2、Y(x,y)=xyを
に代入すると、確かにλ(x,y)=2xが得られる。
このとき、微分方程式(1+y2)dx+xydy=0の一般解である積分と完全微分方程式(2x+2xy2)dx+2x2ydy=0の一般解である積分が同じφ(x,y)=x2+x2y2=C(Cは任意定数)だからといっても、ψ(x1,y1)から出発して任意の微少変分(dx,dy)方向への積算値を線積分した
の関数ψ(x2,y2)が任意の点(x2,y2)に関して一意に定まるわけではありません。何度も出てきたように
なのですから。
このことは5.(3)2.で説明した通りですが、今一度、5.(3)2.で述べたベクトル場について説明する。(x1,y1)=(2,2)から(x2,y2)=(4,4)までを二経路 (2,2)→(3,2)→(4,2)→(4,3)→(4,4) と (2,2)→(2,3)→(2,4)→(3,4)→(4,4) で線積分すると
(長さは1/20倍にしてある)となり、線積分値は二経路で一致しない。しかし
(長さは1/110倍にしてある)では、一致する。
両者の違いは、各点(x,y)に於けるベクトル(X,Y)の傾きは同じだが長さの分布が違うことから生じた。積分因子とは上記の線積分が経路に依存しなくなるように各点のベクトルの長さを調整する関数のことです。どの様に調整するのかというと各点に於けるベクトルの長さがある関数φ(x,y)のgrad(勾配)になるように調整する。つまり微分方程式X(x,y)dx+Y(x,y)dy=0の一般解F(x,y)=C(Cは任意定数)が等高線となるように調整するのです。
一般解の一つの曲線に沿って(つまり等高線に沿って)移動すると隣の解曲線との間隔が密になったりまばらになったりする。等高線が密の所は傾斜が急でまばらな所はなだらかであることを思い出して下さい。その密度の変化に対応するようにそれに垂直なベクトルの長さ(傾斜を意味する)を調整する。つまり間隔が狭いところでは長く、広い所では短くする。実際の所、5.(3)4.の一般解の集合図(等高線図)を上記のベクトル図に重ねてみると、前者ではベクトルの長さが等高線間隔と対応していないが、後者では対応していることが解る。
このとき、一般解の曲線の一つ一つに沿ってベクトルの長さを調整していくのだが、曲線に沿って出発するとき最初に取る基準値の大きさには任意性がある。この任意性から解るように積分因子となる関数は沢山ある。実際すでに述べたようにλ(x,y)=x、λ(x,y)=3x、・・・、λ(x,y)=C・(1+y2)-1/2、λ(x,y)=x・(1+y2)-1/2、・・等々はいずれも積分因子となる。
二変数の場合には積分因子λ(x,y)は必ず存在することは解っているのだが、それを求めるのは必ずしも簡単ではありません。積分因子λ(x,y)を決めるにはφ(x,y)=F(x,y)が解らねばなりません。完全微分方程式の場合、φ(x,y)は5.(2)1.の(4)式の様に簡単に求まりましたが、今は先ず完全微分方程式にするためのλ(x,y)を求めなければならないからです。
それは結局、λ(x,y)X(x,y)dx+λ(x,y)Y(x,y)dy=0が完全微分方程式であることの必要十分条件
から得られる
なる線形一階偏微分方程式を未知関数λ(x,y)について解くことと同じです。
当然のことですが解λ(x,y)はX(x,y)やY(x,y)の関数形に依存します。その関数の値を考慮した上であらゆる(x,y)点で上記の式を満足する関数λ(x,y)を求める。
偏微分方程式を解くのは一般に難しいのですが、ある特殊な仮定が成り立つ場合には簡単に解くことができます。そのような例を幾つかを説明します。
5.(3)で出てきた (1+y2)dx+xydy=0 を例に説明します。いまこの微分方程式の積分因子がxだけの関数λ(x)であることが解っているとします。その場合上記の偏微分方程式は
となります。
となり、唯の積分の式に簡単化できます。そのため両辺をそれぞれの変数で積分すれば
となる。これはまさに5.(3)1.で述べた積分因子です。この解き方が7.(1)、(2)で重要になります。
同様に (1+y2)dx+xydy=0 を例にして説明します。この場合も積分因子がyだけの関数λ(y)であることが解っているとします。その場合上記の偏微分方程式は
となるので、両辺を積分して
となる。これも積分因子です。実際に条件式に代入してみると
となることから了解できる。5.(4)3.と同じ範囲のベクトル図を描いてみると
(長さは1/5倍にしてある)となり、矢印の長さの分布状況は変わっているが、線積分値は二経路で同じであることが解る。つまり、この場合も矢印に垂直な曲線群の隣同士の間隔が矢印の長さに逆比例するように矢印の長さが調整されている。
このとき5.(4)3.で説明した積分因子の意味を考えれば明らかなように前記二つの解を積の形で組み合わせた
もすべて (1+y2)dx+xydy=0 の積分因子となることに注意してください。一般に積分因子は無限に存在します。ところが後で述べるように、熱力学第二法則の縛りを入れると一つに定まります。
以下の議論には積分因子λ(x,y)は必ず存在することと、λがxだけの関数λ(x)の場合の求め方が解れば十分ですので、一般的な解き方の説明は省略します。
前項までの結論を3.章で説明した事柄に当てはめてみます。
U2−U1=dUと置くと熱力学第一法則(エネルギー保存則)は
と表される。ここで内部エネルギーUが状態量であるからと言って、Q-V平面上にある連続曲面U(Q,V)が存在する訳ではありません。d’Qやd’Wは状態量ではないので物体系に出入りする熱の総量Qや仕事の総量Wを座標軸に取ることはできません。同じQやWの値(Vは共通)でもUは一つに定まらないので、QやWは状態方程式を表す変数にはなれないからです。
熱力学第一法則を
の形に書き直したとき、3.(2)2.で説明したように物質系の状態変化に伴って出入りする熱の総量は状態変化の経路ごとに異なるから熱量の変分d’Qは完全微分ではありません。
しかし、UとVは状態量なので状態方程式の変数となり得ます。このときd’Qが定まらない、つまりdUとdVの変化経路に依存すると言うことは5.(3)の議論で言うところのX(U,V)=1とY(U,V)=P(U,V)に関して完全微分の条件が成り立たないということです。つまり
と言うことです。もちろんそのとき、着目している物質系のP(U,V)がどの様な関数になるのか何の情報もありませんが、とにかくU-V平面上の経路に沿ってd’Qを加え合わせると始点と終点が同じでもQの線積分値は経路毎に違った値になることだけは確かです。
内部エネルギーUは直接測定することは難しいので、Uが幾ら状態量だからといっても状態変数としては適当ではありません(ただし、次で説明するエントロピーSは経路図の状態変数として重要です6.(4)2.参照)。普通は変数Uを温度Tと体積Vなどに変換しておきます。
ただし温度に関しては絶対温度Tの完全な定義はまだ行われていないので、1.(1)で説明した様な任意に選んだ寒暖計(例えば水銀寒暖計や気体温度計など)の目盛りから適宜検定される実用温度tを採用します。そして絶対温度Tと区別するために小文字の t で表す。このとき圧力PはtとVの関数P(t,V)となります。
いま着目している物質系の内部エネルギーUは状態量で物質系の温度tと体積Vの連続的な関数U(t,V)と考えることができるので
が成り立ち
となる。これらの量は今までの議論の
に相当する。
仕事は状態量ではないので
のd’Wは完全微分ではありません。だから’(ダッシュ)を付けています。ところが
と書いてみると明らかなようにVは状態量ですから、これは完全微分dVとなります。同じことがd’Qについても言えないだろうか。つまりd’Qに何かを掛けたら、それが完全微分になるような関数はないのだろうか?そのような関数はあります。
のd’Qは5.(2)で説明した意味で完全微分ではないが、二変数の場合5.(3)3.で証明したように、ある関数λ(t,V)を乗じることにより
がある関数φ(t,V)の完全微分になるようにする事ができます。つまり
が常に成り立つようにできる。ただし、熱力学では積分因子を普通
の様に書きN(t,V)を積分分母と呼んだりするが、単に逆数をその様に置くだけなので、今後はN(t,V)自身を積分因子と言うことにする。そのとき上式は
と表される。ここに出てくるd’Q/N=dSはある関数のφ(t,V)=S(t,V)の完全微分だから、それをt−V座標上の任意の経路図(ダイアフラム)に沿って線積分した値
は経路によらない。つまりS(t,V)は系の状態(t,V)が定まれば一意に定まる状態量です。このSのことを熱力学ではエントロピーと言う。クラウジウスが1865年の論文の中でその様に命名しました。これはギリシア語の”変化”を意味する言葉(τροπη)から取ったものだそうです[ギリシャ語の”中(en)”+”変転(trop)”]。
このとき積分因子N(t,V)は当然
の関数形に依存するので、U(t,V)やP(t,V)が変数(t,V)のどの様な関数か解らないとN(t,V)は決まらない。
1.(1)で述べたように、実用温度tは温度を測るための温度計を構成する物質の性質に依存する。そのため同じ熱力学的現象でも、用いる温度計によってX(t,V)やY(t,V)の関数形は異なる。その為当然温度計ごとに積分因子1/N(t,V)の関数形は異なる。確かに同じ状態変化経路(もちろん座標変数の実用温度tが別の実用温度t’に変わると経路図は変わる)をたどる熱力学的現象なら、どの様な温度計を用いようとd’Qは同じはずです。また、d’Q/N(t,V)が状態量になる事情は同じです。しかし、このように温度計毎に異なる積分因子1/N(t,V)に何か意味はあるのだろうか?
これは温度計の説明のとき述べた”大いなる疑問?”と同じことで、ここで議論が行き詰まってしまう。ところが、実は思わぬ所から解決の糸口が出てくるのです。
次に述べる熱力学第二法則を用いて定義される絶対温度をTと書くと、熱力学第二法則からd’Q/T=dSが状態量であることが保証される。そのためN(t,V)の関数形など解らなくても関数N(t,V)を絶対温度Tとすれば良いのです。次章6.で定義する絶対温度Tは、そこで説明するように物体の性質には依存しない[つまりX(t,V)やY(t,V)の関数形に依存しない]が、物体の状態がきまればその値が一意に決まる状態量です。
そして実用温度tやVで表した積分因子N(t,V)は、何度も強調するように物体の性質に依存している。しかしN(t,N)=Tと置くTは熱力学第二法則より物質の性質に依存しない。ならば、T=N(t,V)を逆に解いて実用温度tについての関数t=f(T,V)にしてみれば、この関数こそ物質の性質に依存する実用温度tを絶対温度Tで表す関数であることになる。この関数f(T,V)こそ、現実の温度tを絶対温度Tに関係づけるものです。
すでに述べたように積分因子となる関数は無限に存在します。しかし熱力学第二法則による制限を課すことにより一つに定まる事になる。
このとき異なる温度計物質を用いた温度計の温度目盛をt’で表せばX’(t’,V)やY’(t’,V)の関数形はX(t,V)やY(t,V)の関数形と異なるであろうから積分因子のN’(t’,V)もN(t,V)と異なるであろう。そのため、当然のことであるがt’=f’(T,V)の関数形はt=f(T,V)とは異なる。そのためこの関数形をそれぞれの場合に実験で求めれば絶対温度Tをそれぞれの実用温度t、t’、t”・・・・から決定できることになる。
そして、様々な物質、様々な熱力学的現象を利用して定めた実用温度計t、t’、t”、・・・・を絶対温度Tと対応付けておき、ある一つの熱浴の温度をそれらの様々な温度計で測ってみる。それらすべての温度計が同じ絶対温度を指し示せば、絶対温度の正当性、さらには熱力学第二法則の正当性が実証されることになるであろう。これは素晴らしい着想です!
実際にその関数形を求める方法は7.で説明します。
ここで注意してほしいのは、エントロピー自体は熱力学第二法則とは無関係に導入・定義できることです。積分因子N(t,V)が存在してd’Q/N=dSを完全微分にできるということは数学的な真理なのですから。これはクラウジウスが気付いていたことです。
熱力学第二法則が関係するところは、積分因子N(t,V)を絶対温度Tにすればよいというところにあります。そのとき絶対温度Tはどこから導入・定義されたのかと言えば、まさしく熱力学第二法則によって導入されたのです。それはケルビン(W.トムソン)が1848年に気付いた事柄です。クラウジウスは、トムソンの導入した絶対温度の逆数こそ積分因子そのものだと言うことに1850年に気付きます。実はケルビンもクラウジウスに僅か遅れて(1851年)独自にそのことに気付いています。
その当たりを明らかにするために、いよいよ熱力学第二法則の説明に入りましょう。
ここはたいていの熱力学の教科書に載っているので簡単に説明します。参考文献の朝永先生の本などを参照されて下さい。
熱力学第二法則は、様々な熱機関に於いて高温物体だけからは動力を得られないと言う経験事実に基づく。熱機関を働かせるには必ず低温の熱の捨て場が必要なのです。これは熱が絡む動力機関の非常に奇妙な性質です。カルノーやトムソンはこのことの理由が解らず、混乱し思い悩みます。その事実がエネルギー保存則とどういった関係なのかも悩ましい事柄でした。
そういった状況の中でクラウジウスはそれこそ熱の本性であると考えて、それをそのまま受け入れることにしました。つまり、それを一つの自然法則とすればよい。そして、彼はこの新法則はエネルギー保存則と何ら矛盾するものではなくてエネルギー保存則と互いに補い合うものであることを明らかにします。この法則を用いればカルノーの考察に内在していた疑問点を解決する事ができたのです。
熱力学第一法則のところで強調しましたが、法則を保証するものは、あらゆる場合に絶えず新しくその正しいことが確認されてきた長年にわたる人類の経験した事実しかないのですから、これを一つの法則の形に昇華させても何ら不都合はありません。その表現方法は種々ありますが、普通以下のいずれかの表現を用います。
1.何らかの他の変化を残さずに熱は低温の物体から高温の物体に移ることはできない。(Clausiusの原理)
2.循環的な過程によって、一つの物体から熱を取りだしそれを当量の仕事に変える様な機関はあり得ない。(Thomsonの原理)
Carnotサイクルは、準静的過程であって、二つの断熱過程(d’Q=0)と二つの等温過程(dt=0)とから構成されている循環過程である。準静的過程は可逆であるから、Carnotサイクルは可逆である。
Carnotサイクルを具体的に実現する可逆熱機関としては、実用温度t1の高温熱浴と実用温度t2の低温熱浴の間で働く下図の様な機構を考えればよい。ここで温度測定は1.(1)で述べたような実用温度計による温度で測られたものとする。またシリンダー内の気体は任意の実在気体でよい。
ただしシリンダーとピストン間や重りを載せるテーブルなどには摩擦は無く、ピストンやテープルの質量は無視できるように補助の滑車と重りの装置などが理想的に配置されているものとする。シリンダーは熱浴と接する部分以外は断熱的な壁でできているものとする。また重りや重りを載せる台の上下間隔は準静的過程が実現されるように無限小の大きさに分割されているものとする。
![]()
Carnotサイクルは次の手順で運行される。まず、状態の等温変化(等温線t1に沿って)によって着目している物体を状態A(V0,t1)からB(V1,t1)に移す。この途中熱量Q1が加えられたとする。次にその物体を断熱変化(断熱線に沿って)させて状態C(V2,t2)に移す。さらに続けて、等温曲線t2に沿って熱量Q2を取り去りながら状態D(V3,t2)に移す。この三番目の状態Dから、状態の断熱変化によって初めの状態A(V0,t1)にもどす。
この四つの過程の途中ではすべて仕事が行われる。熱力学第一法則によって
が成り立つ。
上図矢印の方向に変化させた場合(順Carnotサイクル)加えられた熱量Q1の方が取り去られた熱量Q2よりも多く、その分だけこの物体は外部に仕事を行う。したがってサイクルの途中で熱が仕事に変わり、それとともに熱量Q2が、t1からt2へ、つまり高温熱浴から低温熱浴へ運ばれることになる。
Carnotサイクルは可逆機関だから上図と逆方向に運転することもできる。その場合(逆Carnotサイクル)加えられた熱量Q2より取り去られた熱量Q1の方が多く、その差に相当する分だけ、外部から仕事がなされる。したがってサイクルの途中で仕事が熱に変えられ、それとともに熱量Q2が、t2からt1へ、つまり低温熱浴から高温熱浴へ運ばれたことになる。
ここで注意してほしいのはCarnotサイクルは具体的な熱機関を指すのではなくて、二つの温度が定まっている熱浴と温度差無しに熱量の交換を行う(準静的過程には必須の条件)ために、断熱膨張による冷却と断熱圧縮による加熱により二つの等温過程が結びつけられている可逆循環過程を指すのである。
それを実現する熱機関の機構やその作業物質は任意で良いことに注意。つまり、作業物質は実在気体で良いし、変化の途中に凝縮液化が起こっても良い。またこのときの温度は実在物質で作られた1.(1)で言う実用温度tで測ったもので全くかまわない。
このようなCarnotサイクルを行う任意の可逆熱機関に対して次のカルノーの原理が成り立つ。
二つの温度の熱浴の間で働く任意の可逆機関(Carnotサイクルはその一例)が一循環行程で吸収、排出する熱量の比は機関の構造や作業物質によらず、二つの熱浴の温度のみに関係する。つまり
の関係式が常に成立する。[カルノーの原理]
[証明]
もし、上記の結論が正しくないとすると、二個のCarnotサイクルを組み合わせて第二種の永久機関が実現できる。今物質Tを用いてCarnotサイクルを運転し、他方物質Uを用いて全く同一のCarnotサイクルを作り、Q高1=Q高2となるようにして正反対に運転する。この場合はTが温度t高の高温熱源からもらう熱量を、Uがその熱源にもどすことになっているから、温度t高の熱源からの熱の出入りはゼロである。
そのとき、もしQ低1<Q低2ならば、温度t低の熱源のからUが奪う熱量の方が、Tがその熱源へもどす熱量より多いということである。熱力学第一法則からW1−W2>0となり、正の仕事が外界に対して成されることになる。従って他に何の代償も残さずに、温度t低の熱源から熱を奪って仕事にかえたことになり、熱力学第二法則に反する。
このとき、もしQ低1>Q低2ならば、TとUの運転を逆向きにして、結局上と同じ矛盾に到達する。結局熱力学第二法則から次の等式が成立しなければならない。そのため
が成り立つ。つまり、高温熱源と低温熱源でやり取りされる熱量の比は、その熱源の温度にのみ依存して熱機関の構造やそれを構成する物質の性質によらない。
[証明終了]
ここで用いられている仮定は、今考えている機関が(準静的に働く)可逆機関であると言うことと熱浴の温度が定まっているということだけであることに注意。この仮定が成り立てば、カルノーの原理は熱力学第一法則と第二法則の当然の帰結だと言うことです。
熱伝導についてのClausiusの原理を否定すれば熱量Q低を低温熱浴t低から高温熱浴t高へ何の代償変化も無く移動する熱機関Uが存在することになる。次に別の順CarnotサイクルTで高温熱浴t高からQ高1(>Q低)の熱を受け取りW1=Q高1−Q低の仕事を外部に行えば、TとUの機関を合わせて結局高温熱浴t高から熱(Q高1−Q低)をうばってすべて仕事W1に代えたことになり、Thomsonの原理の否定が成立する。(下左図)
またThomsonの原理が否定されれば、高温熱浴t高から熱量Q高1を奪って他に何の変化も残さずにすべて仕事W1=Q高1に変える熱機関Tが作れる。その熱機関が行う仕事を、逆CarnotサイクルUに使用する。この機関TとUが1行程を終えれば、熱量Q高2−Q高1=Q低2が低温側熱浴t低から高温側熱浴t高へ移り、他に何の変化も生じないことになる。すなわちClausiusの原理の否定が成立する。(上右図)
以上から、第二法則の二つの表現が同等な内容であることが示される。
ケルビン(W・トムソン)は1848年にカルノーの可逆サイクルの真の意味や意義を研究しているとき、これを使えば、温度を測るのにどういう物体を使うかに関係なく、エネルギー保存則だけに基づいて、温度の普遍的で絶対的な尺度を作ることができることに気付いた。その考え方は次のようなものです。
カルノーサイクルを二つの温度t高とt低の間に動かします。ここで温度t高とt低は任意の熱力学的尺度で測った1.(1)で述べた実用温度です。熱力学第二法則によるとカルノー機関が高温熱浴及び低温熱浴とやり取りする熱量の比は高温熱浴と低温熱浴の温度のみに関係し、それ以外には依存しない[カルノーの原理]。そのため、カルノー熱機関の効率は
ここでf(t高,t低)がどの様な関数であるのかは熱力学第二法則からは解りません。解っているのは高温熱浴と低温熱浴の温度のみに関係することと、その温度差が大きければ効率が大きくなると言う経験事実だけです。
そのとき、この可逆機関は上記の熱浴間で働く機関としては最高の効率を実現するものです。もしそれよりも効率の良い機関(超能機関)があれば、その超能機関と上記の可逆機関を逆方向に運転させたものを同じ熱浴間で同時に働かせれば低温熱浴から熱を汲み出してそれをすべて仕事に代える第二種の永久機関ができることになり、熱力学第二法則が否定されるからです。
このf(t高,t低)の関数形は以下の様な三つの可逆機関を考えるともう少し限定する事ができます。機関AとCは、温度t1の高温熱浴から同じ大きさの熱量Q1を受け取る。機関Cは温度t3の低温熱浴にQ3の熱を放出する。一方、機関Bは温度t2で機関Aから放出された熱Q2を受け取り、温度t3の熱浴にQ3の熱を放出する。このとき機関A+Bは機関Cと同じ熱源の間で運転される一つの可逆熱機関と考えることができる。そのため組み合わされた熱機関A+Bの効率は単一機関Cと同じでなければならないので、両者が熱浴t3に排出する熱Q3は等しくなければならない。
三つすべてに対して前記の式を適応すると
となる。これを恒等式
に代入すると
となる。つまり右辺もt1とt2のみの関数とならねばならない。このようなことが満たされるのは関数fが以下の形をしているときです。
ゆえに、φ(t2)はf(t1,t2)とf(t2,t3)の積からは削除されて
となる。そのため温度t高とt低の二つの熱源の間で作用する可逆機関に対しては
のように書くことができる。
このときケルビン(W.Thomson)は次のように考えました。このφの関数形がどの様なものか、熱力学の法則からは決して定めることはできない。第一法則は温度に無関係だし、第二法則からは温度に関係していることが言えるだけですから。ならば
と置いて、逆に上記の関係を絶対的な温度Tを定める式だと解釈すればよい。熱力学第二法則から可逆熱機関は最高の効率を実現するものであり二つの熱源の温度が定まれば一意に定まるのだから、そこに熱力学第一法則を援用すると、それは更にQ低/Q高が一意に定まることを意味する。その値を用いれば、物質に依存しない[可逆機関の効率は作業物質に依存しない]熱力学的な絶対温度を定めることができる。つまり二つの熱浴の間に可逆機関を働かせて熱の出入りを測定してその熱量でもって二つの熱浴の絶対的な温度を定めればよい。つまり
と決めればよい。これは素晴らしい着想です!実際のケルビンの推論はかなり込み入っていますが、結論としてはこのように言って良いと思います。
W.Thomson,"On an absolute thermometric scale founded on Carnot's theory of the motive power of heat, and calculated from Regnault's observations"[Cambridge Philosophical Society Proceedings for June 5, 1848; and Phil.
Mag., Oct. 1848.]
まことに素晴らしい着想なのですが、残念ながらその様な理想的な可逆熱機関は現実には存在しません。そのため具体的な測定をおこなうことはできません。残念ながら、ここで行き詰まってしまうように見えます。
ところが、ケンビンの導入した絶対温度にはとても素晴らしい性質が備わっていることをクラウジウスは発見します。
前節で述べたように、可逆熱機関(カルノーサイクルも含まれる)を動かすとき、高温の熱源と低温の熱源の温度を一定にしておくと熱機関が受け取る熱と放出する熱との比は、常に一定の値を持つ。つまり、T高をTout、T低をTin、Q高をQout、Q低をQinと書き換えると
となる。クラウジウスはこの関係を次のように書き直した。
この関係式は1サイクルの間にカルノーの可逆熱機関に入る熱を、その時の絶対温度で割ったものは、同じサイクル中に放出する熱をその時の絶対温度で割ったものに等しいことを表している。その時、入るものを正、出るものを負と考えて、系に出入りする熱をその時の絶対温度で割ったものを集計すると、1サイクル終わるとゼロになることを示している。
クラウジウスは上記の関係式を更に一般的な可逆熱機関に適応します。任意の物体が準静的可逆的に任意の変化(図ではA点から出発してBを経過してもとの状態に戻る)をしたとする。このサイクルを多数の熱源間で働くカルノーサイクルを結合したもので置き換えることができる。下図の様に、それぞれの温度圧力の状況に対応するカルノーサイクルを適当な熱源を選択しながら働かせればよい。そのとき下図の反平行で逆行する部分は互いに打ち消し合うので積算経路の中らか削除できる。
各熱源T1、T2、・・・とやり取りする熱量をQ1、Q2、・・・とする。そのとき受け取るときを正、放出するときを負として以下の和を取る。各カルノーサイクルには前記の関係式が成り立つので、全体としても
が成り立つ。無数のカルノーサイクルに分けた極限では和Σは線積分になるので
と書くことができる。すでに何度も説明したようにQは状態量ではなくdQが完全微分ではないので、それを明示するために’(ダッシュ)を付けておく。このとき積分路を図のようにB点で二つに分けてそれぞれをT、Uとすると
と書ける。準静的変化Uを逆にたどれば熱の吸収、放出が逆になるので
となり、従って
となる。このとき積分路T、Uは任意に取ることができるので、これらの積分値は途中の道によらないで、状態Aと状態Bだけで決まることを意味する。そのため
と置くことができる。
つまり、”d’Q/Tを積算した値Sは状態が決まれば一意に決まる量[状態量]である”ことが証明できたことになる。
つまり、カルノーの原理とは「可逆機関が最高の熱効率を達成するもので、その効率は熱をやり取りした熱源の温度(実用温度でよい)だけの関数である」と言うものでしたが、そのことから絶対温度が定義[ケルビンによる]でき、それ絶対温度Tを用いて系がやり取り(もちろん可逆的に)する熱量d’Qをやり取りした熱源の絶対温度Tで割ったd’Q/Tを系の変化経路に沿って積算した値は、その系がたどった経路によらず系の最終的な状態のみによって決まることが言える[クラウジウスによる]わけです。
ただし、この段階では熱源の絶対温度Tをどうやって測るのかはまだ何も解っていません。それを可能にするのが次の考察です。
前節6.(3)で明らかになったことは、前々節6.(2)で定義した絶対温度Tを用いると準静的可逆変化で物質系に出入りする熱量の変分をd’Qとすると、物質系の準静的可逆変化の経路に沿ってd’Q/Tを積算した値は、その物質系の状態量となることです。
一方5.(6)で導入したエントロピーdS=d’Q/N(t,V)の積算値も物質系の状態量である。それならばこのN(t,V)こそ絶対温度Tそのものだと考えれば良いのではないか。そうすれば、その関数N(t,V)を使って絶対温度Tを現実の温度tで表すことができるだろう。
5.(6)3.で説明したようにまさに革命的な素晴らしい着想です。このことに最初に気付いたのはクラウジウスです。ケルビンもクラウジウスに僅か遅れましたが独自に気付きました。
ここで非常に重要な仮定をします。今まで積分因子N(t,V)はtとVあるいはtとPの関数として説明してきましたが、絶対温度Tはカルノーサイクルの定義からも明らかなように物質系の状態量です。もちろん1.(1)2.で説明した現実の温度計を用いて定めた温度tも物質系の状態量です。その両者は値が違っていても同じ物理量だと考えるべきです。
そうならば、積分因子Nを絶対温度Tとするのですから積分因子は温度tのみの関数のはずです。つまりTとtは一対一の対応が付くのだから、tが定まればその時の体積や圧力とは無関係にTは一意に定まる量のはずです。だからT=N(t,V)=N’(t,P)と置くのではなくてT=N(t)と置いても良いはずです。これはプランクが「熱力学」§129で特に注意していることです。そのため今後はT=N(t)と仮定しましょう。そうすると積分因子を求めるための偏微分方程式を解く操作が劇的に簡単になります。5.(5)2.のやり方で解けばよいのですから。
以上のような事情で今後は6.(3)で導入したd’Q/TもdSと書いてエントロピーの完全微分、またそのその積算値Sをエントロピーと呼ぶことにしましょう。
絶対温度の測定の話に入る前に、積算値S(エントロピー)について幾つかの補足説明をします。
可逆過程で出入りした微少熱量を絶対温度で割ったものはエントロピーの微少変化
であるが、これを書き直すと
となる。これは可逆過程において出入りした熱量が絶対温度にエントロピーの変化をかけたものに等しく、熱が仕事と同じように二つの[状態量である]物理量の積となることを表している。
そのため、TとSを状態変数として経路図(ダイアグラム)を書いてみると、P-V図で面積が仕事量を表したように、T-S図では面積が状態変化に伴って系に出入りする熱量を表すことになる。
T-S図上で起こっている現象を物理的に解釈したり、エントロピーが何を意味するのかを理解するのは難しいが、TdSの意味は極めて明瞭なので、エントロピーは状態を表す変数としてとても重要です。
物質系の状態を示す状態変数(当然状態量)としては、普通TとS、PとVの4つの内の二つを変数に選びます。熱力学の法則は[物体系の状態変化]を[熱の出入り]や[仕事の出入り]と結びつけるもので、熱の出入りに関係する変数がTとS、仕事の出入りに関係する変数がPとVだからです。
理想気体を例にして、状態方程式のグラフを示します。3.(2)1.でPVTを状態変数とした状態方程式のグラフを示しましたが、STV、STP、SVPに関する状態方程式のグラフを描くこともできます。
熱力学第一法則を熱力学第二法則によって書き直した
から出発します。理想気体の場合
各状態方程式をグラフにしてみる。簡単化の為にnCv=nCp=nR=1としてT、V、PはT0、V0、P0を単位としたグラフを示す。当然のことですが(T,V,P)=(T0,V0,P0)のときのSをS0=0としている。
S=S(T,V)
S=S(P,V)
上記の簡単化の仮定の基ではS(P,V)のグラフはS(T,V)のグラフのT軸をP軸に置き換えたものと同じになる。
S=S(T,P)
S(P,V)のグラフはS(T,V)のグラフと同じなので省略する。
次に、工学で良く用いられるP=P(S,T)のグラフを描いてみる。今度も簡単化してnCv=nCp=nR=1としてT、PはT0、P0を単位としたグラフを示す。当然のことですが(T,P)=(T0,P0)のときのSをS0=0としている。このグラフはT→0、S→+∞とともにP→0に近づく。
同じくV=V(S,T)のグラフを記す。nCv=nCp=nR=1とし、T、VはT0、V0を単位としたグラフを示す。当然のことですが(T,V)=(T0,V0)のときのSをS0=0としている。このグラフはT→+∞とともにV→0に近づく。
最後に述べた二つの形式は工学でとても重要です。[別稿「熱機関の効率(ガス動力サイクル)」等参照]
不可逆変化のときエントロピーがどの様に変化するのかを前項のP=P(S,T)図を用いて説明します。すでに述べたようにd’Q/Tの積分は、その積分が二つの状態の間で内的に可逆な経路に沿って実行されるときにのみ、エントロピー値の変化を与える。不可逆な経路に沿って実行されたd’Q/Tの積分値は状態量ではなく、可逆な経路に沿って積分が実行された値とは異なる。
不可逆過程である拡散の現象を、6.(1)2.で用いた機構を利用して説明します。ただし、今回はシリンダーとピストンは断熱材で作られており、状態変化の過程で外界と熱のやり取りはしないものとする。またピストンや重りを載せたテープルの質量はゼロになるように適当な滑車と重りの系が組み合わされており、それらに摩擦は無いものとする。簡単化のために実験は真空の部屋の中で行う。
P0=5atm、V0、T0の初期状態から、次に述べる3っの方法で変化させる。
状態0から状態1〜3への過程をT-S座標上に記入すると下図のようになる。
この図から、エントロピーが各過程の不可逆性の程度をはかる目安になっていることが解る。不可逆変化が起こっているときはいつでもエントロピーは増加し、同時に体系のエネルギー利用度は減少する。1の場合は棚に置いた重りをピストンのテーブルに移して圧縮して行けば元の状態0に帰れるが、2や3の状態から0の状態に帰る事はできない。つまり不可逆変化が起こると、その系から仕事を取り出すチャンスを永久に失ってしまうのである。
ここで注意してほしいことは、たとえ不可逆過程で到達した状態2や状態3に対しても状態量であるエントロピーは定まっている事です。それはもちろん前項3.で説明した準静的可逆過程を使って計算した値です。状態3は状態0から温度T0=T3が一定になるように温度T0の熱浴から準静的に熱を加えながら、すこしづつおもりを取り除いて無限にゆっくり膨張させていく。圧力が1atmになるまでに加えられた熱量をQ03とするとS3−S0=Q03/T0で計算できる。また状態2は状態0から準静的断熱変化で温度T2まで膨張させた後、温度T2の熱浴に接した状態で温度をT2に保った準静的等温変化で錘をゆっくり取り除きながら圧力1atmまで膨張させる。そのとき加えられた熱量をQ02とするとS2−S0=Q02/T2で計算できる。このとき準静的な仮定だから、図のピストンシリンダーの左側にある梯子状の錘置き場に錘のやり取りの過程が残っています。そのためその錘を再び無限にゆっくりとピストンに載せながら同じ温度の熱浴に準静的に熱を排出して圧縮して行くと元の状態0に可逆的に復帰できます。
前記のように不可逆仮定で状態3や2に推移させることもできますが、その場合には熱の流入はないのにエントロピーは増大します。また、錘のやり取りの過程が梯子状の錘置き場に残っていませんから、その分の仕事を取り出すチャンスは永久に失われているわけです。
温度差のある熱伝導も不可逆過程です。例えば下図の様な、等圧での相変化を考える。その際定圧(飽和蒸気圧が同じ)だから当然等温変化(20℃)である。このとき左側1は20℃の熱源から無限小の温度差の元で、右側2は80℃の熱源から温度差がある状況で熱をもらい定圧・等温の相変化を伴った膨張が生じる。そして、シリンダーとピストンを含めた状況変化が1と2で同じように起こって、両者とも500kJの熱量を得たとする。
このとき、シリンダー内の水に起こるエントロピー変化ΔSwaterはどちらの場合もシリンダーの内壁に接する部分の水の温度は20℃=293Kと考えて良いので
となる。一方熱浴に起こるエントロピー変化ΔSbathは熱浴の温度が異なるため異なった値となる。
そのため、水と熱浴を合わせた全系のエントロピーは、可逆変化のみで変化する左側1の場合は変化しないが、不可逆変化を含む右側2の場合には増大する。
このとき、2の場合のエントロピーの増大は何処で起こるのかと言うと温度差を生じるシリンダー壁内を熱が移動する過程で生じると考えねばならない。ただしシンリンダーの質量が無視できて、その熱容量が小さい場合や、今の場合のようにシリンダー内の温度勾配が定常的だと考えられる場合は、熱の通過にともなってエントロピーの増大を生み出すが、シリンダー内に含まれるエントロピー自身の変化は無視することができる。
もちろんシンリダー内や熱浴内の撹拌が不十分でシリンダー壁面近くの温度分布が均一で無ければ、エントロピーの生成はシリンダー壁面近くの水や、熱浴の中でも生じることになる。
ここで、”温度差がある場合の熱伝導は不可逆過程なのにエントロピーが計算できるのはおかいしのではないか?”と思われる方がおられるかもしれませんので、具体的な例でもう少し補足します。
下図の様に断熱壁で囲まれた容器を、熱を良く通す仕切りで二等分し、それぞれに温度T1とT3(T1<T3)の同種の理想気体を1モル量ずつ入れる。この状態をAとし、エントロピーをSAとする。次に仕切りを通して高温から低温に熱が移り、やがて二つの気体が熱平衡状態を実現して同じ温度T2になったとする。この状態をBとし、そのときのエントロピーをSBとする。
ここで、状態AからBへ変化したときのエントロピー変化を求めてみる。何度も強調したようにエントロピー変化を求めるためには、状態AからBへ準静的可逆過程で変化させて、上記の系に対して出入りする熱量を計算してdQ/Tを積算しなければならないのですが、6.(1)1.で述べたように、上記の過程は不可逆過程です。そのため可逆過程で同じ変化を実現する方法を考えます。中央の仕切りを断熱壁で置き換え、T1〜T3の範囲で少しずつ温度の違う熱源をたくさん準備する。そして、高温理想気体に接する熱源を次々に切り替えながら無限小の温度差で等温的に温度T3〜T2の熱源へ熱を移して自らの温度をT2とする。
次に、低温理想気体を温度T1〜T2の熱源に次々と繋ぎ替えて等温的に熱を少しずつもらって低温理想気体の温度をT2にする。このとき、外界の熱源に変化を残しながらではあるが、全ての段階の可逆性を保って着目する系の状態をAからBへ変化できたことになる。このとき外界の熱源に残っている熱量移動の痕跡をを逆向きに戻せば可逆的に元の状態に戻れるから確かに可逆過程である。
この可逆過程において系に出入りしたdQ/Tの合計値を計算するのは簡単です。その過程で体積は変化せず、理想気体の定積モル比熱Cvは一定と見なせるので
上記手順を用いれば状態Aから状態Bへ可逆的に変化させることができる。そのとき系のエントロピーは増大しているが、外界の熱源のエントロピーは同じ量だけ減少している。上記の過程をそっくり逆に行えば、外界の熱源に一時的に預けたり、それから借りたりしている熱量をすべて元に戻して最初の状態に復帰できる。つまり状態AからBへの変化が可逆的にできるのは、外界に変化の痕跡をを残しているからです。
もし状態AからBへの変化を外界に頼らず、直接真ん中の仕切り板を通して温度差がある状態で熱を移動させたならば、同じBの状態になるのだが、その場合には外界も含めて元の状態に戻すことは二度とできないのです。
この例を考察すれば、前記の水と熱浴の温度差がある場合の熱移動において、エントロピー変化がそこでの議論の様にして計算しても良いことが理解できる。仮想的に上記例の様な外部熱源を仲介にしていると考えれば良いのです。
異種物質の混合過程も不可逆変化です。下図の様に断熱壁で囲まれた容器を仕切り板で二つの部分にわけ、一方に理想気体Tを圧力P、温度T、体積V1でn1モル、他方に理想気体Uを同じ圧力PとT、体積V2でn2モル積める。この状態をAとし、エントロピーをSAとする。次に仕切り板を取り除き二種類の気体が拡散・混合できるようにする。完全に混合して熱平衡状態になった状態をBとし、そのときのエントロピーをSBとする。
ここで、状態AからBへ変化したときのエントロピー変化を求める。エントロピー変化を求めるには、状態AからBへ準静的可逆過程で変化させて、系に対して出入りする熱量を測定しd'Q/Tを積算しなければならない。可逆変化を実現するために次のようなメカニズムを考える。二つのピストンTとUを考える。ピストンTは気体Tは通さないが、気体Uは自由に通過できるフィルターが張ってあるとする。同様にピストンUは気体Uは通さないが、気体Tは自由に通過できるものとする。このようなフィルターは現実には存在しないが、あくまで思考の上では可能であるとする。
最初に次図の様にピストンUを固定してピストンTだけが動くようにし、温度Tの熱浴に接した状態で準静的等温変化により理想気体Tをゆっくりと膨張させる。このときピストンUに働く圧力は、ピストンUが気体Tを自由に通過させるので、最初の圧力のままです。一方、ピストンTに働く圧力は、気体Uを自由に通過させるため気体Uの中を移動してもならん影響を受けず、通過させない気体Tによる圧力だけが働く。その圧力には気体Tの状態方程式にしたがった等温膨張の圧力変化が生じます。
気体Tが体積V1から体積V1+V2まで膨張するとき、等温変化をさせるために系に加えなければならない熱量は理想気体Tが膨張に際してピストンTを押すことにより外部にした仕事量と同じです。理想気体の場合内部エネルギーは温度が変わらなければな同じなので、エネルギー保存則から明らかなように外にした仕事と等しいだけの熱量を加える必要があるからです。そのためこの変化によるエントロピー変化ΔSv1→v1+v2は
となる。
引き続いて、ピストンUのロックを外して自由に動くようにし、全く同様に温度Tの熱浴に接した状態で準静的等温変化により理想気体Uをゆっくりと膨張させる。このときも気体Tによる圧力は働きません。気体Uの膨張による圧力変化だけが生じます。
このとき気体Uが体積V2から体積V1+V2まで膨張するとき等温的に変化させるために系に加えなければならない熱量は理想気体Uが膨張に際してピストンUを押すことにより外部にする仕事と同じです。この場合も、理想気体の内部エネルギーは温度が変わらなければな同じなので、エネルギー保存則から明らかなように外にした仕事と等しいだけの熱量を加える必要があるからです。そのためこの変化によるエントロピー変化ΔSv2→v1+v2は、同様にして
となる。
両方の変化を加えると結局状態Aから状態Bへ変化したことに伴うエントロピー変化量は
となり、エントロピーは増大する。
最初に述べた中央の仕切り板を取り除くことによる混合では、外界と仕事や熱をやり取りすることはないので、内部エネルギーは変化しません。そのため当然のことですが温度も圧力も最初のTとPがそのまま成り立ちます。
この場合も状態Aから状態Bへの可逆的な変化において系のエントロピーは増大しているが、外界の熱源のエントロピーは同じ量だけ減少している。そのとき、仕事や熱のやり取りの過程は外部の錘の位置変化や熱浴の熱量の変化として外部に残っているので、上記の過程をそっくり逆に行えば、外界の熱浴から一時的借りた熱量や、梯子に錘を乗せることで預けた仕事の結果をすべて元に戻して最初の状態に復帰できる。つまり状態AからBへの変化が可逆的にできるのは、外界に変化の痕跡をを残しているからです。
もし状態AからBへの変化を外界に頼らず、直接真ん中の仕切り板を外して拡散混合させたならば、同じBの状態になるのだが、その場合には外界も含めて元の状態に戻すことは二度とできないのです。
不可逆な過程としては上記の[拡散]、[温度差のある熱伝導]、[異種物体の混合]以外に、[摩擦による熱の発生]、[電気抵抗によるジュール熱]・・・等々がありますが、いずれについても同様な事情が当てはまり、各状態に対して一意に決まるエントロピー値を定義にしたがって計算できます。その当たりは適当な教科書を参照されて下さい。
上記の例は教訓的です。状態変化を可逆過程で行うと、着目する系のエントロピーが変化しても、それを打ち消すだけの変化が外界のエントロピー変化として生じており、系と外界を含めた全体としてのエントロピー変化はありません。ところが着目する系の同様な状態変化が不可逆な過程によって起こると系のエントロピーは必ず増大していく。
エントロピーは系の状態が決まれば、その状態に応じて定まる状態量ですが、孤立系のエントロピーは一定かあるいは増大していくかで、決して減少することはありません。
これがエントロピーという量に関して熱力学理論が教えてくれる驚くべき結論です。
ここで、いよいよ絶対温度Tを現実の温度計で測定した実用温度tに関係付けます。
以後の議論はすべて、プランクに習って、絶対温度を測定する物体やそれと熱量をやり取りする熱機関の作業物質の単位質量当たりについて論じます。そのため示量性変数をすべて単位質量当たりに付いてのものとし、そのことを強調するために小文字で表します。また、特にそうする必要は無いのですが、間違いを防ぐために示強性変数も小文字で表すことにします。
以下の説明のなかで最も重要な事柄は、5.(5)2.で説明した積分因子を求めるための偏微分方程式
を解く操作であることに注目して下さい。
5.(6)3.の式から出発します。6.(4)で述べたように積分因子λ(t,V)中のN(t,V)はtだけの関数と見なせますのでN(t,V)=T(t)とおける。また単位質量当たりで考えているので
となる。ここでsは状態量ですから、dsが完全微分であることの条件式から
が導けます。式の変形をたどれば明らかなように偏微分が常微分に変換できる。そのため、この微分方程式は簡単に積分できて
となります。これは5.(5)2.の計算過程と全く同じですからそこの例と比較してみて下さい。
これによって物質系の絶対温度Tが温度tによって完全に定められる。右辺積分記号内の分子は物質の状態方程式から積分途中の各温度tにおける値が直接得られるのでそれを用いればよい。例えば実在気体を使う温度計の場合はファン・デア・ワールスの状態方程式のような実験式を用いればよい。また分母は、熱力学第一法則(エネルギー保存則)から得られる
を考慮すると、物質が等温可逆膨張するときに外部から流入する熱量、あるいは等温可逆圧縮されるときに外部へ流出する熱量だから、各温度に於ける体積変化にともなう発熱or吸熱量の測定値を用いればよい。
ここで温度計物質として理想気体を用いる場合を説明します。現実には理想気体は存在しないのでその様な温度計はできませんが仮に存在するとしたらの話です。その温度計が示す温度tは理想気体の体積で測るとします。その温度目盛りは次のようにして定められます。単位質量の理想気体(仮に存在するとして)を1気圧の下で沸騰している水の中に浸けてそのときの体積を量ります。それをv100とする。次に同じく単位質量の理想気体を1気圧の下で氷と平衡にある水の中に浸けてその時の体積を量ります。それをv0とします。そうして単位質量の理想気体に定圧の下で(v100−v0)/100の体積変化が起これば1度の変化が生じたとします。
から明らかなようにp=一定の下でvが零に近づけばRは定数ですのでtも零に近づきます。体積は零以下はありませんから体積零の点が温度の零度となります。そして0℃の理想気体温度目盛りは、その時の体積v0を1度に付いての体積変化量(v100−v0)/100で割ったものに他なりません。
ところで理想気体の膨張係数αは1℃=1Kだけ温める際の体積増加分を0℃=T0Kの体積v0で割ったものです。つまり
ですから、0℃の理想気体温度計による温度t0は膨張係数を用いると
となります。
ところで、理想気体の状態方程式から
が成り立ち、また理想気体の内部エネルギーは温度のみの関数で体積に関係しませんから
です。これらを前節の式に代入すると
となりますから、まさに理想気体温度計が示す温度tは絶対温度Tそのものです。そのため理想気体の膨張係数αは0℃の絶対温度T0の逆数
に他なりません。ただし、残念ながら理想気体は現実には存在しませんからその膨張係数を直接測定することはできません。だからこういった意味での絶対温度Tは2.(5)で注意したように定めようがありません。この当たりは7.(2)5.でもう一度詳しく説明します。
このとき次のことに注意してください。上の証明をみれば明らかなように理想気体の状態方程式p=Rt/vのtが絶対温度であると見なせるためには、“理想気体の内部エネルギーは温度のみの関数で体積には関係しない”という事が必要です。
4.(1)3.で注意したように、元々理想気体の定義には未だ定義できていない絶対温度Tを用いて状態方程式を定義し、その状態方程式でもって絶対温度Tを定義するという循環論法の論理矛盾がある。そのため、PV=RTが成り立てば理想気体の内部エネルギーUが体積によらないということが証明できるというような関係ではないのです。
正しくは、“理想気体の内部エネルギーが温度のみの関数で体積に関係しないのなら、p=Rt/vのtをケルビンが定義した意味での絶対温度であると解釈できる。”ということです。またその逆に、そのように解釈できる為には、上で見たように“理想気体の内部エネルギーが温度のみの関数で体積には関係しない。”でなければならないのです。
測定精度を上げるには4.(1)2.で説明した実在気体のJoule-Thomson効果を用いればよい。
Joule-Thomson効果を適用するために、5.(6)1.の式の状態を表す変数を(u,v)から(t,p)に変換しておく。今は単位質量あたりで考えているので
となります。uもvも状態量ですから、duもdvも完全微分となり上記の様な変換ができます。
ここでsは状態量ですから、dsが完全微分であることの条件式から
が得られる。
4.(1)2.で説明したJoule-Thomsonの実験は、もともと理想気体の性質の決定に利用されたのですが、それよりも更に立ち入った意味を与えることができます。
すなわち、気体に、熱の授受なしに(Q=0)体積膨張V2−V1を(単位質量についてはv1−v2を)与える。そのとき気体の単位体積に及ぼされる外部仕事はp1v1−p2v2=wで表される。理想気体では、そのとき温度が一定に保たれるので内部エネルギーは変化せずw=0となる。ただし、実際の気体に対しては、熱力学第一法則(エネルギー保存則)からq=0だから
となる。つまりJoule-Thomsonの実験では状態変化の前後(つまり細孔栓の前後)でエネルギーuが一定に保たれるのではなくて、h=u+pvなる状態量(エンタルピーと呼ぶ)が単位質量について一定に保たれる。このとき実在気体ではp1v1とp2v2は等しくないので当然u1とu2は等しくない。そのため温度が変わるのです。
ここでは、簡単のために細孔栓の両側の圧力差は極めて小さくなるようにして気体を流す。両側に於ける状態を表す量の差をΔを付けて表すと
となる。これがJoule-Thomsonの実験に於いて圧力差Δpで細孔栓を押し出される気体に起こる温度変化Δtの式である。公式から解るようにその時の気体の絶対温度Tと、その温度での実用温度tと絶対温度Tの変化率dt/dTに関係している。そのためこれをdt/dTの微分方程式と見なして積分すれば
となる。右辺積分記号内の分子は物質の状態方程式から積分途中の各圧力pにおける値が直接得られるのでそれを用いればよい。また分母のvは流れている気体のその温度に於ける単位質量の体積、c’pは温度tにおける気体の定圧比熱(実在気体の定圧比熱は温度に依存する)、ΔtはJoule-Thomsonの実験の温度差、ΔpはJoule-Thomsonの実験における圧力差であるからすべて実測可能です。つまり、この実験により気体の絶対温度Tが温度tによって完全に定められる。
右辺の積分はtのみで決まり他の変数には無関係なので、最も測定がやりやすい状況で行えばよい。例えば圧力を一定(1気圧)に保った状態で温度tを変えながら微少な圧力差Δpのもとで気体を細孔栓に流し込み温度変化Δtの測定を行うのである。もちろん各温度でのvやc’pはあらかじめ測定しておく。
そのとき細孔栓の前後の気体温度をその管の中に挿入した水銀温度計で測ったとすると、tは水銀柱の体積膨張によって測られ定義されている実用温度tのことで、その温度計を用いて絶対温度Tが測定されたと言うことになる。
絶対温度の目盛りを決めるために、1気圧の下で沸騰しつつある水の絶対温度をT100、1気圧の下で凍りつつある水の絶対温度をT0としてその差を
と置く。これは絶対温度の目盛り幅を摂氏温度と同じにするための措置です。さらに
とおくと
となり、実用温度tが絶対温度Tに関係付けられる。もちろんJtもJ100も任意の温度計を用いて測定した実用温度tによる、任意の気体に付いてのJoule-Thomson効果の測定値についての積分値です。
具体的な計算は以下のようにやります。すでに述べたようにJの積分はtのみで決まり他の変数には関係しないので、最も測定がやりやすい状況で行えばよい。そのため以下では一定圧力(例えば1気圧、以下の数値例では1kgw/cm2)に限定し、実用温度tの温度計物質として、細孔栓実験に用いたのと同じ気体(たとえば空気)を採用する。
このとき実用温度t0(=0℃とする)とt100(=100℃とする)は絶対温度T0とT100の熱浴[つまり1気圧の下で氷と平衡にある水と、1気圧の下で沸騰している水]と平衡にあるときの単位質量(1g)の温度計気体(空気)の体積から定める。t0のときの体積をv0、t100のときの体積をv100とし、その実在気体(空気)の膨張係数をα’(実在気体のものと言う意味で’を付けている)とする。
気体の膨張係数α’は気体の温度を1℃だけ温める際の体積増加分をt0(=0℃)の体積v0で割ったもの[3.(2)4.参照]ですから実用温度t℃のときの気体(空気)の体積(単位質量)vtは
となる。これが実用温度tと実用温度目盛りの指標となる空気体積vtとの関係を表す式です。このとき
となるから、前述の式は
となる。それでは実際のデータを用いて0℃が絶対温度で何Kとなるか求めてみよう。
圧力pmmHg、温度t℃における空気の密度ρtは下記の関係式で良く近似される。それを用いると0℃に於ける単位質量の体積v0と膨張係数α’は
となる。
また空気の定圧比熱c’Pは圧力1[kgw/cm2]の下で0℃から100℃までの間、ほぼ一定でc’P=0.240cal/(g・℃)です。理想気体の場合は温度に依存しない[4.(2)3.参照]としましたが、実在気体の定圧比熱は温度と共に少しずつ変化します。例えば「化学便覧」(丸善)によると空気は1気圧のもとで100℃では0.24568cal/g・℃、200℃では0.24045cal/g・℃、300℃では0.24043cal/g・℃、500℃では0.24608cal/g・℃です。しかし積分の中で関与する項が他の項に比べて小さい(下表参照)ので、あまり変化しない場合には定数として計算しても差し支えない。
また空気のJoule-Thomson係数Δt/Δpは実験によると下表第二欄の様になる。これはp=1[kgw/cm2]の圧力の空気を、微少圧力差Δpの下で細孔栓を押し出したときの微少温度変化Δtを微少圧力差Δpで割ったものです。単位は[℃/(kgw/cm2)]で表してあります。工学でよく使われる力の単位[kgw]を使っているのでcgs単位系に変換したものを同じ表の第三欄に付け加えておきます。[具体的な測定法は4.(1)2.のJoule=Thomsonの実験を参照されたし。]
であるから下表の様に計算される。
これらの値α’、v0、c’p、Δt/Δpを用いてJ100の数値積分を行う。数値積分は、台形公式では誤差が大きすぎるので、シンプソン公式を用いると
となる。このJ100の値を前項で求めた0℃の絶対温度T0を与える式に代入すると
となる。実験誤差から小数点以下に0.1K程度の誤差はあるが、ここで利用したデータは、まさに-273.25℃が絶対零度0Kであることを示している。つまり、これより低い温度は無いのである。そして、もしこの温度の熱浴が準備できれば理想的な熱機関の効率を1にできる訳です。
計算に使ったExcel元ファイルはここをマウスで右クリックでダウンロード可。htmlファイルに変換していますが、もともとExcelファイルですから、Excelで読み込み再編集可。xls形式で再保存するとExcelファイルにもどります。
任意の摂氏温度t℃を絶対温度TKに直すにはすでに述べたJtを同じように測定・計算して
に代入すればよい。その時、当然の事ながらJt=J0=0を代入するとT=T0となり、Jt=J100を代入するとT=T100=T0+100となる。
ThomsonとJouleが1852年以降の一連の実験から得た値(マッハの本からの引用)を下記に示します。彼らの測定では0℃は絶対温度でT0=273.7Kとなります。これらは新しい目盛りでの温度(T−273.7)と0℃で760mmHgの圧力を示す定積空気温度計による温度t℃との対応関係を示している。
この表において、当然の事ながら0℃と100℃においてT-273.7とtは完全に一致しているが、それ以外の温度ではT-273.7とtは異なる事に注意。ただしその差はごく僅かです。それは、実在気体(空気)が凝縮(液化・凝固)するような温度範囲以外では、実在気体と理想気体の性質の違いはごく僅かだからです。(J.
P. Joule, and W. Thomson, "On the Thermal Effects of Fluids in Motion.Part
II.", Phil. Trans. 1854, vol. cxliv. p. .321-364.)もとデータはこちら
マッハが強調しているように、彼らが得た結果の精度がどこまで精密であるかは、さして重要なことではない。新たに定義された目盛り[絶対温度]へ、その他の個別的な各目盛り[実用温度]を引き直すための原理はどんなものか、熱力学が明快な理論的関係を獲得するに至るための原理がどのようなものなのかがここに示されたのである。
上記の計算に用いた数値は久保亮五偏 「大学演習 熱学・統計力学」裳華房(1969年刊)P59第2章問題[C]の問[45][46]から引用しました。この問題集は若い頃かなり気合いを入れて集中して勉強しました。上記問[45][46]の解答も熟読したのですが、思考力の無い私には、なぜこの問の解答に書かれている方法で絶対温度が求まるのか幾ら考えても良く解らなかったのです。問[45][46]の解答として記載されている内容はプランクの本のそのままなのですが、私には全く意味不明でした。しかし、このHPを書いた今なら、一点の曇りもなく理解できます。理解して初めて言えるのですが、いきなりこの本を読んだのでは絶対に理解できないでしょうね。
1968年の国際度量衡総会 (CGPM) で、水の三重点 (0.01℃) を 273.16 K と定義したため、絶対温度とセルシウス温度の目盛の幅(1度の温度差)は等しく、相互の値には[絶対温度(K)
= セルシウス温度(℃) + 273.15]の関係が成立する。
しかし次の事柄に注意してほしい。1.(1)3.3.で注意したように絶対温度目盛りの基準点は、二つの基準点が必要な摂氏温度目盛りと違って、一点だけで良いのです。そのため、水の三重点を273.16と定めた時点で目盛りの幅も同時に決定されていることになります。もともとは、セルシウス温度目盛りの幅と等しくなるように定めたのですが、今日では、そのこととは無関係に定められていると言って良いと思います。
前節の計算では理想気体との違いが見えてこないので上記の計算を別なやり方でやってみます。まず
の右辺の分母に於いて
が成り立つので、Joule-Thomson係数Δt/Δpが掛かっている項がほぼゼロであると見なせる場合[理想気体ではΔt=0だから完全にゼロ]には
となり、0℃の絶対温度はほぼ膨張係数α’の逆数に近くなる[理想気体の膨張係数αは測定することはできないが、もし求めることができたなら、その逆数に完全に一致する]。これは2.(3)ですでに注目していた事柄です。このような事情から
と置いて計算してみる。下段の式は
と置くことと同じです。そうすると前記の積分は
となるが、ここで
だから左辺を展開してxの1次の項だけを取ると
が得られる。前項の数表のデータを用いて数値積分を行うと
となり、前項と全く同様な結論が得られる。
この解き方から、気体が理想気体に近づけばx→0へ近づき、α’→α(理想気体の膨張係数)に近づくことが解る。
プランクも強調しているように、”絶対温度の設定に対して、熱力学第二法則が原理上からみてはなはだ重要であり、また極端な温度においては欠くべからざるものであることは明らかである。しかし、実際のところ、理想気体の膨張係数αの最も精密な数値は、熱力学第二法則には無関係に、実在気体の膨張係数α’の測定から得られている。”[プランク「熱力学」§164]
その具体的な方法は下記のようなものです。
理想気体ではその法則から明らかなように、一定圧力の下での体積V、及び一定体積のときの圧力Pは、ともに絶対温度Tに比例する。例えば100℃と0℃に於けるそれぞれの値に添え字を付けて表すと
となる。
そこで、下図のような定積気体温度計や、定圧気体温度計を準備する。まず、定積気体温度計に二酸化炭素などの気体を入れ最初の圧力がPであったとする。この温度計をまず沸騰する水の中に入れてその圧力がP100に変わり、次に融解しつつある氷りの中に入れたら圧力はP0になったとする。これから(P100/P0)vの比が解り、グラフの点(1)がとれる。
次にはじめの圧力がPと同じ二酸化炭素の定圧気体温度計を使い、同じことをすれば(V100/V0)pの比がわかり、グラフの点(2)が決まる。
グラフの点(1)と(2)が違うのは、CO2が理想気体ではないからです。初めの圧力を変えて実験すれば、図のような曲線を描くのに十分な点を取ることができる。さらに、気体をCO2から水素やヘリウムなどに代えて同様な実験を行えば、圧力が減少するのに伴って互いに近づいていく曲線群が得られる。すべての実在気体は圧力がゼロに近づくにつれて、理想気体の性質に近くなることが解っているので、様々な種類の気体を、様々な圧力で測定して得られる曲線群を圧力P→0まで伸ばしたところが求める点である。それは1.3661へ収束しているように見える。この値から
であることが解るというわけです。
この方法に原理的な裏付けがあるわけではないので何となく納得できない所があるのは偽らざる心境です。測定できない圧力の部分のグラフが直線的に伸びている保証などないわけですから。これを教えられたとき、ほんとにこのやり方しか無いのか疑問に思われたと思う。その疑問に対する答えが上に述べた熱力学第二法則を用いる方法です。それによって求めたT0の逆数こそ現実には測定できない理想気体の膨張係数αなのです。
クラウジウス=クライペイロンの式は、最初ClapeyronによってCarnotの定理から導かれ、後にClausiusによって厳密に証明されたもので、気相と液相の熱平衡に対して成り立つ関係式です。プランクは彼の「熱力学」§177で、この関係式を用いても絶対温度Tが実用温度tの関数として決定できることを説明しています。
クラウジウス=クラペイロン過程は、一定圧力pに調整され自由に動くピストンで蓋をした円筒状容器に封入された物質の二相平衡系に生じる変化です。系に熱を加えると圧力一定の下で相1(液体)が蒸発して相2(飽和蒸気)に変化して体積が変化する現象です。
液相(相1)にあるときの物質1gが持つ内部エネルギーをu液、気相(相2)にある物質1gが持つ内部エネルギーをu気とすると、この相変化は一定圧力、一定温度の下での準静的可逆過程だから、実用温度で測った転移温度をtとすると
となります。これは
と変形してみれば明らかなように、相転移前と相転移後の1gの物質のギブスの自由エネルギー(状態量の一種)が等くなるように生じる変化です。
このクラウジウス=クラペイロン過程を実用温度dtだけ異なった温度(t+dt)で行わせる事を考える。そのときの圧力は当然pからdpだけ異なった圧力(p+dp)となります。このとき(dp/dt)がp−t座標面における飽和蒸気圧曲線の傾きになります。この傾きに対して成り立つ関係式がクラウジウス=クラペイロンの式です。
この式変形中に含まれている微少変分の意味は非常にわかりにくいので是非、別稿「ファン・デル・ワールスの状態方程式」4.(3)で確認されて下さい。このとき式中のsやvは全て実用温度tの関数ですが、変分dTだけは絶対温度Tで無ければいけないことに注意して下さい。ここで、飽和蒸気圧曲線に沿った積分を実用温度tに関して実施すると実用温度計の温度目盛 t を絶対温度 T に結びつける関係式T(t)が得られます。
つまり、実用温度で測った各転移温度tに於ける転移熱q転と相変化に伴う体積変化(v気−v液)を用いることにより、各実用温度tが絶対温度Tで何度になるかが決定できるのです。(∂p/∂t)は転移温度tにおける飽和蒸気圧曲線の傾きであり、積分は飽和蒸気圧曲線に沿って実用温度t0からtまでおこなうと、そのときの温度差(t−t0)が絶対温度で(T−T0)に相当するということです。具体的な計算法は7.(2)のジュール=トムソン効果を用いる場合と全く同じですので、そこの議論を御覧下さい。
ここまでのプロセスをたどれば解るように原則的には第二法則から導かれる方程式は何れもこれを絶対温度Tの測定に利用することができる。そして、このことは5.(6)3.の□囲み枠内に書いた様に熱力学第二法則の検証に利用できる。これらはプランクが彼の著書「熱力学」の中で何度も強調していることです。
このページを読まれた方はお解りのように絶対温度とエントロピーは極めて抽象的・数学的概念です。熱力学が解らなくなる最大の原因は、ここのプロセスで出てきた式中の実用温度tを最初から絶対温度Tで置き換えてしまうところにある。本来積分因子として導入された絶対温度Tは積分因子の位置のT(t)だけで、そこに限って使っても良いものです。それ以外は絶対温度を使うべきではないのです。
ところがここで見てきたように目盛りの幅を実用温度(例えば摂氏温度)と同じにしたら(つまり100℃を373.15K、0℃を273.15Kにしたら)、その範囲外の温度でも普通の温度計による測定値と絶対温度との差は(T=t+273.15の違いを除いて)極僅かです[7.(2)4.のThomsonとJouleの測定結果参照]。
そのため多くの教科書では最初からtを絶対温度Tに置き換えて論じており、そのときdT/dt=1とした式を公式として利用している。そうしていきなり、エントロピーが状態量であるという性質を最大限利用して様々な物理量と温度Tを結びつける熱力学の有用な公式を導き出し、利用している。
絶対温度の最も重要な事柄はその値に下限があり、それは知ることができるということと、その下限から測られた絶対温度の高さが、熱機関の効率に深く関わっていると言うことです。それと絶対温度を導入する事でエントロピーという状態量を定義でき、熱の出入りが、熱の代わりにエントロピーを用いることによって、とても解析しやすくなることです。熱力学の法則から導かれる重要な公式・法則の大半は、エントロピーが状態量であるという事実を最大限活用して得られたものです。それにも拘わらず最終的に、それらの公式・法則にはエントロピーは含まれていない。
このとき、前述したように絶対温度は極僅かの誤差を無視すれば普通の温度計で簡単に測定[もちろん-273.15℃が基準ですが]することができる。ならば、現実の温度計で測定した温度tを絶対温度Tと同じと見なして、dT/dt=1として差し支えないだろう[7.(1)〜(3)の式でt=T、dT/dt=1としたものが普通の教科書に出てくる式だということを確かめて下さい]。そして、後は絶対温度とエントロピーが導入された過程やその理由を忘れてしまってもかまわない。最終的な公式には現れないのだから、エントロピーという量の意味はよく解らなくてもかまわない。温度という量が測定できさえすればよい。エントロピーが状態量であるという事実を最大限活用して現実の問題に適応してゆけばよい。
熱が絡んだ現象で最も重要な概念が温度ですが、これは極めて難しい量です。熱力学第一法則のdU=dQ−PdVだけでは、その式の中の何処にも温度が入ってきていない。d’Q=TdSとしてd’QをTとdSに置き換えたとき、初めて温度が出てくるのですが、このとき状態量でないd’QをTで割ったd’Q/T=dSとして導入されるSこそ状態量エントロピーです。熱力学第二法則はまさに温度とエントロピーを導入することと等価です。つまり熱力学第二法則は熱をエネルギーの他の形態から区別して特徴付けるものであり、熱の発生、移動の不可逆性を温度という量との関係で説明するものです。Carnotサイクルの議論を読めばその当たりの事情は了解できる。
熱力学の最大の成果はエントロピーの微分が完全微分である[エントロピーが状態量である]という事実
を最大限利用して、熱が絡む現象の間に成り立つ関係を導き出すことです。
そのときエントロピーの導入に伴って現れた温度という物理量が、様々な物質定数、物理量を結びつける関係式に必然的に介在してきて重要な役割を果たす。
そのとき、Tとtの差は無視できるのだから、それを同等と見なせばよい。実際に知りたいのは温度という量[実用温度である摂氏温度に273.15を足せば直ちにわかるもの]と様々な物理量との関係なのですから。熱力学とはまさに熱が絡む[当然温度が絡む]現象が示す様々な法則を研究する学問なのだから。
事実世の中の大半の教科書はその立場で議論が展開されている。だからそれらの教科書を幾ら読んでも絶対温度とは何か、エントロピーとは何かがよく解らないのですが、利用する立場から行くと差し支えないわけです。その当たりは割り切って考えるべきかもしれません。しかし、そのやりかたが熱力学を難しくしており、温度やエントロピーという量の理解を困難にしている。本来ここで述べた議論をしっかりやってから上記の応用に入るべきなのだと肝に銘じたほうが良い。
最近の高校教科書ではcal(カロリー)をやめていきなりJ(ジュール)という単位で熱量や比熱を説明しています。しかし、私はそのやり方には反対です。1.で述べたような従来の導入方法できちんと教えた方が良い。そのやり方では温度とか熱量の測定法しか定義していなくて、それらの量が何を意味するのか全く明らかではありませんが、その導入法こそ温度と熱量の本質をよく表していると思います。それであるからこそ後に明らかになった熱の仕事当量や、後に導入された絶対温度とエントロピーの概念が重要になるのですから。
この章は未完成です。前章7.(4)「まとめ」で説明した t=T and dT/dt=1とする普通の教科書のやり方で、熱力学第二法則から導かれる成果の幾つかを説明するもので、このHPのテーマとは関係ありません。公式集として利用するために、今後少しずつ付け加えてゆこうと思います。
以下のように、色々な観点で見ると熱力学第二法則の特徴が浮かび上がってくる。
理想気体の状態方程式はその気体の状態を一意に定めます。その時の変数温度T、圧力P、体積Vはその状態が決まれば一意に決まる状態量です。
このときすでに出てきた内部エネルギー: U やエンタルピー: H=U+PV も、その状態の物質系に附随する状態量でそのときの状態が決まれば一意に決まる量です。さらに付け加えると
エントロピー: S、
ヘルムホルツの自由エネルギー: F=U−TS、
ギブスの自由エネルギー: G=H−TS=U+PV−TS
も物質系の状態が定まれば一意に定まる状態量です。
これらの量が状態量であることはPVやTSなどの量も状態が定まれば一意に定まる状態量であることから明らかです。このときdS=d’Q/T=(dU−d’W)/T=(dU−PdV)/Tとおけるからと言って、決して Q’=ST 、 W’=PV 、 S=(U−W’)/T=Q’/T 、 F=U−TS=W’ 、 G=PV+W’ と言った置き方はできません。仕事W’や熱量Q’は状態変化の経路によって異なり状態量ではないのですから。つまり
であってそれは積分経路に依存し、決してW’=−(P2V2−P1V1)やQ’=(T2S2−T1S1)ではありません。ところがP1V1やP2V2、T1S1やT2S2は状態量であって(P1,V1)や(P2,V2)、(T1,S1)や(T2、S2)が決まれば一意に決まる量です。
PVやTSが状態量であるということは、P-V座標平面に垂直にPV座標を取ると一つの連続な曲面になるということです。つまり
ということです。TSについても同様です。
その為
が言えますので
が成り立ちます。そのため、U、H、F、GはいずれもP、V、T、Sの中の二つの変数でそれぞれの微少変化量をあらわすことができます。そのとき変数の微少変分の前に掛かっているT、P、V、Sはその二つの変数の関数とみなせる。
このときUやH、F、Gも状態を記述する状態変数と見なすこともできますが、6.(4)3.で注意したように普通はTとS、PとVの4つの内の二つを変数に選びます。熱力学の法則は[物体系の状態変化]を[熱の出入り]や[仕事の出入り]と結びつけるものです。そのとき熱の出入りに関係する変数がTとSで、仕事の出入りに関係する変数がPとVですから。
この中でS以外は直接測定することができます。dSも可逆変化の場合はd’Q/Tのように熱量の出入りd’QとTを測定すれば直接測定できますから、これらの量を状態変数とするのは妥当な選択です。
ここで導入したエンタルピーH、ヘルムホルツの自由エネルギーF、ギブスの自由エネルギーGは、単に状態量である内部エネルギーUに状態量であるPVやTSを足したり引いたりしたものでしかありません。それだけのものですが、物質系の状態が決まれば一意に決まる状態量であることは確かです。
これらはある特殊な状態変化に伴う法則を解析するとき便利だから導入されました。エンタルピーHを導入する意義については別稿「反応熱と熱化学方程式」1.(3)を、Gを導入す意義は7.(3)などを御覧下さい。Fを導入する意義については適当な教科書を参照して下さい。
U(S,V)、H(S,P)、F(T,V)、G(T,P)に対して、U=f(S,V)、H=f’(S,P)、F=f”(T,V)、G=f’”(T,P)の様な関数が存在します。つまりそれぞれの量に対する状態方程式が存在すると言うことです。これはちょうど理想気体のとき説明したP=f(T,V)の様なものです。理想気体のときそうしたように、それらの状態方程式を最大限利用して現実の問題を解いてゆけば良い。その具体的な例を別項「ファン・デル・ワールスの状態方程式」4.(2)で紹介していますので御覧下さい。
Maxwellの関係式とは、dU、dH、dF、dGの変分がそれぞれの独立変数についての完全微分であることから成り立つ、完全微分の必要十分条件[5.(2)1.参照]そのものです。元々はエントロピーSが状態量であり、そのためdSが完全微分になることから派生してきたものです。
これらの関係式は、とても便利なものですが、安易に利用するとエントロピー導入の意義を見失ってしまい、わけがわからなくなります。熱力学が難しいのもその当たりにあるように思いますので、充分注意して利用しなければなりません。
関係式そのものは、前記dU、dH、dF、dGの定義式から、その完全微分必要十分条件式[5.(2)1.参照]を書き記せば直ちに求まります。詳細は適当な参考書を御覧下さい。
この題材は高校物理の範囲を逸脱していますが、いつか解りやすく説明するHPを作りたいと思っていたものです。
準静的可逆機関を用いて状態変化を追跡したとき出入りする熱量を測定することによってしか定義できない絶対温度をどうして知るのか?現実にはその様な理想的可逆機関などありはしないのですから、現実に存在する熱機関(熱膨張する温度計もその一種)から得られる仮の温度tから絶対温度Tを導くしかない。そのとき、どの様にして現実の温度計が示す仮の温度tからその物体の絶対温度Tを知るのだろうか?この疑問に答えてくれないと絶対温度の概念は意味を成さないし、決して理解できない。
いままで熱力学の本は沢山(専門書は10冊以上、入門書・解説書のたぐいを含めれば30冊以上)読みましたが、上記の疑問に真面目に答えてくれているのはプランクの本だけでした。[マッハの本にも書いてあるがプランクの本より更に解りにくい。プランクの本を翻訳された芝亀吉氏が書かれた「熱力学」(岩波全書1967年刊)の§16、§19、パウリの物理学講座第3巻(講談社1976年)第2章、岩波講座現代物理学の基礎「古典物理学U」第W部(碓井恒丸著)(岩波書店1973年)等にもそれらしいことは書かれていますが、極めて難解で凡人にはとうてい理解しがたい]。しかし、そのプランクの本も理解するのがとても難しい。良く解らない所を推察しながら解釈したのでずが、本当にここに書いたような理解でよいのか疑問が残ります。いずれにしても第二法則を理解する鍵は積分因子の数学的な考察を十分に行うことではないでしょうか。このHPを作っていて改めて感じました。
世の中に熱力学の教科書は沢山ありますが、私の見た限りでは積分因子と偏微分方程式の積分を説明しているのは上記の本だけのようです。ただし説明はしてあるのですが、普通に読んだのではそうであるとはなかなか解りません。かなり推測を働かせるとなるほどそうなのだと解るような書き方です。いずれにしてもそこの所が説明されていないと絶対温度を真に理解するのは難しいと思います。
絶対温度の導入過程についてはどの教科書にも説明してあり、下記の本なども秀逸で解りやすい。しかし積分因子そのものを説明していないので、これらの本から絶対温度を真に理解するのは難しいと思います。
高校生が熱力学を勉強されるには下記の本を薦めます。これらはいずれも工学系の本ですが、とても丁寧に説明してあり解りやすい。そして、これらを読めば熱力学が実用的にいかに重要な学問であるかがよく解る。
<数学参考書>完全微分方程式と積分因子はたいていの本で説明されていますが、なかなか解りづらいのでできるだけ沢山の本を読み比べて見られることを薦めます。私が参考にしたものの幾つかを記しておきます。
参考文献に上げた本が古い本ばかりで申し訳ありません。いずれも若い頃買って読んだ本です。読むには読んだのですが理解したという状況にはほど遠く、また教職についてからは教師の仕事が忙しくて振り返る暇もなく、それらの内容はほとんど忘れてしまいました。しかし、プランクの「熱力学」§160〜§164に絶対温度の本質的な説明があるということは記憶の片隅にありました。記憶はしていましたが正直なところその内容は良く理解できていなかったのです。
退職して時間の余裕ができましたので、2ヶ月かけて昔に読んだ熱力学関係の本を集中的に読み直してみました。読み直してもなかなか本質が見えなかったのですが、5.(6)3.の□囲み枠内に書いたことが見えてきて、プランクの教科書の§131、§153〜§164に展開されている説明は、表現こそ違うが5.(5)2.で説明した積分因子λ(x)の偏微分方程式を解く手順そのものだと解ったときすべてが理解できました。長年の胸のつかえが下りた感じです。恥ずかしながら、やっとこの歳になって絶対温度とエントロピーのなんたるかが理解できました。
プランクの「熱力学」は今日絶版で手に入りませんが、機会がありましたら是非読んでみられることを薦めます。