物理を解説 ♪
ツイッター用のシェアボタン フェイスブック用のシェアボタン はてなブックマーク用のシェアボタン ライン用のシェアボタン
ツイッター用のシェアボタン フェイスブック用のシェアボタン はてなブックマーク用のシェアボタン ライン用のシェアボタン

状態方程式の微分形

最小限必要な偏微分の知識。
作成:2004/7/15
更新:2023/11/24

全微分形式

理想気体の圧力,体積,温度を結びつける式については前に 数式 であるとした.つまりp,V,Tの内の 2 つの量が決まれば,残りの 1 つは自動的に決まってしまうということだ.そこで,体積Vは温度Tと圧力pの関数V(T,p)であると考えて,次のような式を作ることが出来る. 数式 なぜこのような表現が出来るかについてきちんと説明しておこう.軽々しく納得していいところではない.

温度と圧力がわずかに変化することで,体積がわずかに変化したとする.その変化は次のように書ける. 数式 これを変形してやれば, 数式 と書けるが,ΔTΔpについて無限小の極限を考えれば,分数で表した部分は微分の定義式そのものである.ただし高校で習う 1 変数のみの微分とは少しだけ違っていて,初めの項の中の分数の部分はpの値を固定したままでのTによる微分を表しており,2 項目の分数の部分はTの値を固定したままでのpによる微分になっている.

このように他の変数を固定して行う微分を,通常の微分と区別して「偏微分」と呼び,記号もdの代わりに∂を使って区別する.計算自体は全く難しくない.考えている以外の変数を定数と見なして普通に微分すればいいだけのことだ.

偏微分を書き表すのに,例えば, 数式 と書いた場合,これは変数x,yを固定してzで微分することを意味するのだが,もっと簡略化して, 数式 のように表すことがある.固定した変数を右下に書いておくわけだ.変数が何であるかが分かっている場合には,右下に変数名を書くのを省略するし,そうする方が普通なのだが,熱力学ではどの変数を一定に保ったまま状態を変化させるかというところが重要なので,忘れないようにメモ代わりに書いておく習慣になっている.

それで初めに書いたような表現が出来ることになるわけだ.もう一度書いておこう. 数式 ここで使っているdVdTdpというのは,ΔVΔTΔpとは違って「無限小変化」であるというニュアンスを表している.つまりこの関係式は無限小変化の極限で成り立つものであるということを意味している.そんな繊細な関係式は現実世界での大きな変化には当てはめられないではないかと思うかもしれないが,そういう場合には積分を使って変化の合計を計算してやって,現実世界の大きな変化と比較できる式に変えるのである.そういう考えはこの記事でも少し後の方で出てくる.

この (1) 式のような表現を「完全微分」あるいは「全微分」と呼ぶ.もちろん圧力pや温度Tについてもそれぞれp(V,T)T(V,p)という関数であると考えることが出来て,全く同じ操作を経て次のような全微分形式で書くことが出来る. 数式 こうして 3 通りの表現が出てきたが,どれを使っても本質的には同じ式である.その時々に応じて一番便利だと思うものを使うことになる.


相関係数

色々な気体の性質の違いを比較するために,実験でp,V,Tの間の関係を調べ,幾つかの相関係数として表すことが行われる.(1) 式が (2),(3) 式に比べて便利な点は,式の中で使われている偏微分係数が,よく使われる相関係数の定義に近い形になっていることである.

例えば (1) 式の第 1 項の偏微分は,圧力を一定に保ったまま温度を変化させた時の体積変化を表しているが,これを体積で割ったものは「定圧膨張率」あるいは「熱膨張率」として良く使われる量である. 数式 また第 2 項の偏微分を体積で割ったものにマイナスをつけたものは「等温圧縮率」として良く使われる量である. 数式 マイナスが付くだけで難しく見えてしまうが,このマイナスは大した理由ではない.圧力が増せば体積は減るのでこの偏微分の値は常にマイナスになってしまう.係数の値が常にマイナスになるのはかっこ悪いのでそれを防ぐために付けてあるだけだ.

これらの係数を使って (1) 式を次のように表現し直しても良い. 数式 もちろん (2),(3) 式も全く同等な式なのでうまく変形すればこれと全く同じ式に行き着くこともできるわけだが,そのためにはそれぞれの偏微分係数の間にある関係を知らなくてはならない.それは少し後で説明しよう.

相関係数は他にもあるが,直接的には意味をつかみにくいかも知れない.例えば先ほど出てきたκの逆数である「体積弾性率」と呼ばれる量がある. 数式 この定義をいきなり見せられても意味がつかみにくいだろう.これは結局は圧縮率の逆数に過ぎなくて圧力変化と体積変化の比であるので,同じ事ではあるが 数式 と理解した方が分かりやすい.これを (2) 式の第 1 項の偏微分に代入すれば先ほど (1) 式でやったような書き換えができそうだが,話の流れに関係ないのでここではやらないでおこう.

(3) 式の第 2 項にある係数∂p/∂T,温度によって圧力がどう変わるかという意味で「圧力係数」と呼ばれているが,あまり馴染みがないだろう.これについてはかなり後の方で少しだけ使うことになるかも知れないし,私の説明ではその範囲まで手を出さないかも知れない.

まだ大して深くない内容のところであるにも関わらず,すでに似たような関係式が沢山出て来てしまっている.熱力学とはこういうものだ.この先,意味も考えずに片っ端から丸暗記しようとすると非常に効率が悪いことになることは予想が付くだろう.だから覚えるべきは式を作り出す手順だけでいい.


理想気体の場合

この微分形式は高校で習うような状態方程式とは比べ物にならないほど便利な表現であって,理想気体以外についても使えるように出来ている.

なぜなら,偏微分係数の値が一定でなければならないという制限がないので,そこは現実に合わせて自由に変化する部分だと考えてやればいいからである.ずるいやり方だと思うかも知れないが,もっと良く言えば,より柔軟性のある表現になっているわけだ.

すると理想気体の状態方程式というのは (1) 式に含まれる状況の一つに過ぎないということか.では理想気体のみを表す微分形の状態方程式というのはどんな形をしているのだろう確かめておこう.

V=nRT/ppで偏微分することにより, 数式 となる.またVTで偏微分して, 数式 となる.これらを (1) 式に代入して形を整えてやれば,理想気体の状態方程式を微分形で表したものが出来上がる. 数式 なかなかすっきりしていて気持ちいい.しかし,本当にこれでいいのかモル数を表すnも気体定数Rも消えてしまっているではないか.確かにこの表現は元の式とは完全に等価というわけではない.しかし初期条件さえ与えてやれば,この式を積分することで元の式が復元できるのである.やってみよう.

積分するためには積分範囲を考える必要がある.ある初期状態(p_{0},V_{0},T_{0})から始めて,任意の状態(p,V,T)にまで到達することを考えよう.すると, 数式 と計算できる.ここで初期状態(p_{0},V_{0},T_{0})において初期条件p_{0}V_{0}=nRT_{0}が成り立っているとすると,任意の状態でpV=nRTが成り立つことになり,元の式が復元できたことになる.


全微分で表せる条件

上の積分計算で元の式が復元できたのはこれが「全微分形式」だからである.ところでまだ全微分の本当の意味を話していなかった.全微分形式で書けるためには条件があって,上の例ではその条件が自動的に満たされていたのだ.

例えば,3 つの変数(x,y,z)があって,それらの微小量の間の関係が次のような形で表されているとする. 数式 これを積分して関数z(x,y)を求めようとした時,先ほどのように(x_{0},y_{0})から出発して(x,y)まで進むことを考えるだろう.もし高校で学ぶような 1 変数の関数の積分ならば計算の方法は一通りしかないので問題はない.しかしこの場合は 2 変数関数になっているために,その計算の道筋は無数に考えられるわけだ.

どういうことかと言えば,例えば,yy_{0}に固定したままxだけ移動して積分を行い,その後でyだけ移動して目標地点に到達することも出来るし,初めにxを固定しておいてyを先に移動して積分することも出来る.しかしどの道筋を通ってもzの値が同じになるなんて保証は何もないのだ.普通は選んだ経路によって違った値になってしまう.

しかし次の条件を満たしていれば,どんなコースで積分してもその結果が一通りに定まることが分かっている. 数式 これは数学の問題なので,ここではその理由を説明はしない.詳しくは数学の教科書を調べてもらいたい.ここまで分かっていれば数学も楽しく学べるだろう.

(1) 式はこの条件を満たしていることが分かる.結局,関数をxyで偏微分するとき,全ての点でどちらを先にしようとも結果が同じになるようなら問題ないということだ.

熱力学において全微分が大切な意味を持つ理由が分かるだろうかもし温度や圧力の変化のさせ方の違いによって体積が変わってしまうなら,状態方程式が成り立たないことになるではないか.変化のさせ方によらずに,ある状態なら必ず同じ値を取るような量を「状態量」と呼ぶ.全微分の形に書けなければ「状態量」の資格がないわけだ.

全微分の条件を満たさないものは「不完全微分」と呼ばれる.熱力学にはそういう量もあり,近いうちに出てくることになるだろう.


実在気体の近似式

現実の気体は理想気体の状態方程式に厳密には従わず,少しズレがあるということを前に話したが,次のような「ファンデルワールスの式」は現実の気体についての実験値とかなり良く一致することが知られている. 数式 abとはごく小さな値の定数であって,気体の種類によって異なる.もしa=b=0ならばこの式は理想気体の式と一致する.この式の妥当性について細かい議論をするつもりはないのだが,簡単な説明だけしておこう.

理想気体と実在気体の差は何によって生じるのだろうか.一つは,実在気体の分子の一つ一つには大きさがあることが原因だ.そのために気体分子は体積Vの広がりよりも少し狭い範囲しか運動することが出来ない.それでその分のbを引いておいてやるとうまくいくのである.

もう一つは,気体分子同士が接近するときに「ファンデルワールス力」と呼ばれるごく小さな引力が働くことによる.静電力によって原子核と電子の位置がわずかにずれることで電荷が外に漏れて生じる力であるが,aはその分だけあたかも大き目の圧力がかかっているかのように凝集していることを表している.

この近似式についても,先ほどと同じように微分形で表すことが出来るはずである.今からちょっとチャレンジしてみたい.


偏微分の計算

ところが先ほどのようにVpTで偏微分したくとも,この式は複雑すぎてVについて解く事ができない.∂V/∂T∂V/∂pを計算するための何かいいテクニックはないだろうか

例えば,(1) (2) (3) 式にdp=0という条件を課してやると, 数式 という 3 つの式が出来る.ここで (1') 式のdVを (2') 式に代入して整理すれば, 数式 が得られる.偏微分をあたかも微小量同士の割り算であるかのように計算しても問題ないと思わせるような結果だが,これはたまたまなので気をつけよう.pの値が固定されているせいで残された変数がVTの 2 つのみになっており,偏微分がそれらの変数の変化の比を表している状況だからこの関係が言えるのである.3 つ以上の変数が残っていればこのような単純な関係は成り立たない.

このような関係式を使えばVについて解けないと言っていた先ほどの困難が回避できて目的が果たせそうだ.pTについてなら簡単に解くことができるからである.

ところで (1') 式のdVを (3') 式に代入するとどうなるだろうか. 数式 この結果はしっかりと心に留めておく必要がある.「微分なんて単なる微小量同士の割り算だ」などいう考えが間違っていることを思い知らされるのはこの瞬間である.割り算っぽいことにはなっているが,マイナスが余計に付くではないかある変数がV(T,p)のように残りの変数の関数になっているとき,あるいは同じ事だが,変数同士がF(V,T,p)=0のような陰関数による条件で結ばれている時にはこのようなことが起こるのである.この辺りは気を付けておかないと変形の途中で符号が逆になって,とんでもない結論を導いてしまう危険があるところだ.

これらの関係を駆使すれば,先ほどの「圧力係数」が実は,β/κとして表すことが出来るということや,(1) (2) (3) 式が同じ中身のものであることなどを簡単に確かめることが出来るだろう.

さらに (5) 式の左辺を,先ほどの (4) 式のような「割り算に似ている」やり方で変形して少し整理してやれば, 数式 となり,対称的で美しい関係が導かれる.この公式は「マクスウェルの規則」と呼ばれている.このマクスウェルはあの電磁気学で有名なマクスウェルのことである.彼は熱力学や統計力学でも天才的な大活躍をしており,あちこちに名前を残している.マクスウェルの名前を冠した法則はこれからも色々と出てくるので混乱しないように気をつけよう.

ファンデルワールスの方程式についての計算をするつもりだったが説明が長くなり過ぎた.もう終わろう.計算できそうなところまで説明したので,あとは読者の研究に任せることにする.問題なく計算できるかどうかについての責任も持たない.

この近似式ではVp,Tに対して一意に定まらないので完全微分の要件を満たさず,よってこの形に表すべきではないのかも知れないが,そのような領域はごく一部であってそこを避けて考える限りは問題はないのではないかと思う.

とにかくこのようなわけで,微分形で表す限り,これからは理想気体にこだわらないで話を続けてゆけることになる.これまで科学は堅物だと思っていた人に誤解が解けたかどうか聞いてみたいところだ.科学はこのように結構柔軟に現実に対処できるのである.


修正履歴

(2023/11/24)
 変数の無限小変化を表す記号の説明と積分との関連の説明を追加.

(2023/11/24)
 偏微分が割り算のように計算できる特別な場合について,その根拠の説明を追加.



趣味の物理学書店の案内バナー

EMAN物理note出張所の案内バナー