電子レンジで食品を温めたとき、そのエネルギーはどこから来たのでしょうか。電磁波のエネルギーが食品中の水分子の運動エネルギーに変わった――つまり、エネルギーの「形」が変わっただけで、エネルギーの総量そのものは増えも減りもしていません。では、自動車のエンジンではどうでしょう。ガソリンの化学エネルギーが燃焼によって熱に変わり、その熱がピストンを押して力学的な仕事に変わります。ここでもやはり、エネルギーは形を変えているだけで、全体の帳尻は必ず合っています。
この「エネルギーは生まれも消えもしない。形を変えるだけだ」という自然界の根本原理を、熱と仕事と内部エネルギーの関係として数式に落とし込んだものが熱力学第1法則です。
熱力学第1法則を理解すると、次のような幅広い問題が解けるようになります。
- 熱機関の効率計算: エンジンやタービンが燃料の熱エネルギーのうち何割を仕事に変換できるかを定量的に評価できます
- 冷凍・空調サイクルの設計: 冷蔵庫やエアコンが消費する電力と移動させる熱量の関係を明確にできます
- 化学反応のエネルギー収支: 反応熱や生成エンタルピーの計算を通じて、化学プラントの設計に役立ちます
- 宇宙機の熱設計: 人工衛星が太陽光から受ける熱入力と放射冷却のバランスを解析できます
本記事では、第1法則の数学的な表現から出発し、各種の熱力学過程への適用を一つずつ丁寧に解説します。最後にPythonによるp-V図の可視化を通じて、数式で得た結論を視覚的に確認します。
本記事の内容
- 系・外界・境界の定義と、状態量・経路量の本質的な違い
- 熱力学第1法則の数学的定式化と符号規約
- 各種過程（定容・定圧・等温・断熱）での第1法則の適用
- 断熱過程のポアソン関係式の完全導出
- エンタルピーの導入動機と定圧過程での有用性
- Pythonによるp-V図の可視化とエネルギー収支の比較
前提知識
この記事を読む前に、以下の概念を理解しておくと理解が深まります。
- 理想気体の状態方程式 — $pV = nRT$ の意味と使い方
- 熱力学の基本概念 — 熱・温度・平衡の基礎
- 温度と熱平衡 — 温度の定義と第0法則
系と外界 — 「どこに注目するか」を決める
熱力学の問題を考えるとき、まず「何について考えるか」を明確に切り分ける必要があります。料理に例えると、鍋の中のスープについて議論したいのか、コンロの炎も含めて議論したいのかで、話の枠組みがまったく変わります。この「注目する対象の範囲」を決めることが、熱力学の出発点です。
系（System）
注目する対象そのものです。たとえば「シリンダー内のガス」「ビーカーに入った水溶液」「タービンを通過する蒸気」など、問題に応じて自由に選べます。系の選び方によって同じ現象でも立式が変わるため、問題文をよく読んで適切に系を定義することが大切です。
外界（Surroundings）
系以外の全てです。シリンダー内のガスを系に選んだなら、シリンダーの壁もピストンも、その外側の大気も全て外界に含まれます。
境界（Boundary）
系と外界を分ける面です。この境界を通じて、熱や仕事のやり取りが行われます。境界が熱を通すか通さないか（断熱壁か透熱壁か）、境界が動くか動かないか（ピストンのような可動壁か、固定壁か）によって、系の振る舞いが決まります。
系と外界を定義したら、次に考えるべきは「系の状態を記述する量」です。ここで、熱力学には他の物理分野にはない独特の区別――状態量と経路量――が登場します。
状態量と経路量 — 山登りのアナロジー
山の高さと歩いた距離
ある山の頂上を目指すことを考えてみましょう。頂上の標高は、北側の急斜面から登っても、南側の緩やかな登山道を行っても、同じ1,000mです。標高という量は経路に依存せず、「麓」と「頂上」という始点と終点だけで決まります。
一方、実際に歩いた距離はどうでしょうか。急斜面を直登すれば2km、ジグザグの登山道なら5km、大きく迂回すれば10kmかもしれません。歩いた距離は経路に大きく依存します。
熱力学でも、まったく同じ区別が重要になります。
状態量（State Function）
系の状態だけで値が決まり、途中の経路に一切依存しない量を状態量と呼びます。温度 $T$、圧力 $p$、体積 $V$、内部エネルギー $U$ が代表例です。系が状態Aから状態Bに変化したとき、$\Delta U = U_B – U_A$ は、途中でどのような過程を経たかに関係なく、始状態と終状態だけで決まります。これは山の標高差が経路によらないのと同じです。
微小変化を表す記号には、完全微分であることを示す $dU$、$dT$、$dp$、$dV$ を使います。
経路量（Path Function）
途中の過程（経路）に依存する量を経路量と呼びます。熱 $Q$ と仕事 $W$ がその代表です。同じ始状態Aから同じ終状態Bへ変化する場合でも、途中の過程が違えば、系がやり取りする熱や仕事の量は変わります。これは山頂への歩行距離が経路によって変わるのと同じです。
微小変化を表す記号には、完全微分ではないことを強調して $\delta Q$、$\delta W$ を使います。$d$ ではなく $\delta$ を使うことで、「この量は状態の関数ではなく、経路に依存するのだ」と明示しているわけです。
なぜこの区別が重要か
この区別を理解していないと、「$Q$ や $W$ は経路が変われば値が変わるのに、$\Delta U = Q – W$ の結果は常に同じになる」という第1法則の本質をつかめません。お金に例えると、銀行口座の残高（状態量）は一意に決まりますが、その残高に至る入金と出金の内訳（経路量）はいくらでも異なりうる、ということです。収入だけが大きくても支出が同じだけ大きければ、結局残高の変化は同じになります。
状態量と経路量の区別を押さえたところで、いよいよ熱力学第1法則の数式を見ていきましょう。第1法則は、まさにこの「状態量（内部エネルギーの変化）」と「経路量（熱と仕事）」を結びつける式です。
熱力学第1法則の定式化
エネルギーの出納帳
第1法則の数式に入る前に、家計簿のアナロジーで全体像を掴んでおきましょう。
銀行口座の残高変化を考えてください。「残高の変化 = 収入 – 支出」という関係は、どんな複雑なお金のやり取りがあっても必ず成り立ちます。ここで、残高の変化は始まりと終わりの残高だけで決まる量（状態量）であり、具体的な収入パターンや支出パターン（経路量）がどうであれ、差し引きの結果は一意です。
熱力学第1法則はこれとまったく同じ構造です。
- 残高の変化 $\leftrightarrow$ 内部エネルギーの変化 $\Delta U$
- 収入 $\leftrightarrow$ 系に加えられた熱 $Q$
- 支出 $\leftrightarrow$ 系がした仕事 $W$
微分形
熱力学第1法則を微小変化で書くと、次のようになります。
$$ \begin{equation} dU = \delta Q – \delta W \end{equation} $$
この式は「内部エネルギーの微小変化は、系に加えられた微小な熱から、系が外部にした微小な仕事を引いたもの」と読みます。左辺の $dU$ は完全微分（状態量）ですが、右辺の $\delta Q$ と $\delta W$ はそれぞれ不完全微分（経路量）です。経路量同士の差が、きれいに状態量の変化になる――これが第1法則の美しさであり、本質です。
積分形（有限変化）
有限の過程（始状態1から終状態2への変化）では、上式を積分して次のように書けます。
$$ \begin{equation} \Delta U = Q – W \end{equation} $$
ここで各量の意味を改めて整理します。
- $\Delta U = U_2 – U_1$: 内部エネルギーの変化。始状態と終状態だけで決まります
- $Q$: 過程全体で系に加えられた正味の熱量
- $W$: 過程全体で系が外界にした正味の仕事
符号の約束
物理量の正負は以下の規約に従います。
| 量 | 正の場合 | 負の場合 |
|---|---|---|
| $Q$ | 系が熱を受け取る（加熱） | 系が熱を放出する（冷却） |
| $W$ | 系が外界に仕事をする（膨張） | 外界が系に仕事をする（圧縮） |
この符号規約は「系を中心に考える」という視点に立っています。系にエネルギーが入ってくるとき（熱を受け取るとき）$Q > 0$、系からエネルギーが出ていくとき（仕事をするとき）$W > 0$ です。家計簿で言えば、収入がプラス、支出がプラスに対応します。
なお、教科書によっては $dU = \delta Q + \delta W$ と書き、$W$ を「系にされた仕事」として定義する場合もあります。物理学系の教科書は $dU = \delta Q – \delta W$、化学系や一部の工学系の教科書は $dU = \delta Q + \delta W$ を採用する傾向があります。本記事では物理学系の $dU = \delta Q – \delta W$ の規約を用います。
準静的過程における仕事の表現
では、仕事 $W$ を具体的にどう計算するのでしょうか。ピストン付きシリンダーの中にガスが入っている状況を想像してください。ガスが膨張してピストンをゆっくり押す場合、ガスがする仕事は「力 $\times$ 距離」です。ピストンの面積を $A$、ガスの圧力を $p$ とすると、ピストンにかかる力は $F = pA$ であり、ピストンが微小距離 $dx$ だけ動くと、仕事は $\delta W = pA \, dx = p \, dV$ となります（$A \, dx = dV$ は体積変化）。
$$ \begin{equation} \delta W = p \, dV \end{equation} $$
有限の過程では、これを積分して次のようになります。
$$ W = \int_{V_1}^{V_2} p \, dV $$
この積分はp-V図（圧力-体積図）における曲線の下の面積に等しくなります。p-V図はエンジニアにとって非常に重要なツールで、仕事の大きさを視覚的に比較できます。等温過程と断熱過程で同じ体積変化をしても、p-V曲線の形が異なるため、曲線の下の面積（= 仕事）が変わる、ということがグラフから一目でわかります。
ただし、この $\delta W = p \, dV$ の関係は準静的過程（系が常に平衡状態を保ちながら、無限にゆっくり変化する理想的な過程）でのみ成り立つことに注意してください。急激な膨張や爆発のような非平衡過程では、系内の圧力が一様ではないため、$p$ を一つの値で代表できず、この式は使えません。
第1法則の枠組みが整いました。次に、この法則を理想気体に適用するために、内部エネルギー $U$ の温度依存性を見ていきましょう。
理想気体の内部エネルギー
なぜ温度だけの関数なのか
理想気体とは、分子間の相互作用が一切なく、分子自体の体積も無視できるほど小さい、という理想化されたモデルです。この仮定のもとでは、分子が持つエネルギーは運動エネルギー（並進、回転、振動）だけであり、分子間のポテンシャルエネルギーはゼロです。
分子の運動エネルギーは温度だけで決まります（統計力学的には、各自由度に $\frac{1}{2}k_BT$ のエネルギーが等分配されます）。したがって、理想気体の内部エネルギーは体積や圧力に依存せず、温度のみの関数となります。
$$ U = U(T) $$
この性質は理想気体の状態方程式からの帰結ではなく、分子間相互作用がないという理想気体の定義そのものに由来しています。実在の気体では、分子間引力があるため、体積が変わると分子間距離が変わり、ポテンシャルエネルギーも変化します。そのため、実在気体の内部エネルギーは一般に $U = U(T, V)$ です。
内部エネルギーの微分形
内部エネルギーが温度のみの関数であることから、その微小変化は次のように書けます。
$$ \begin{equation} dU = nC_v \, dT \end{equation} $$
ここで $C_v$ は定容モル比熱（1モルの気体を体積一定のまま1K上昇させるのに必要な熱量）です。なぜ $C_v$（定容比熱）がここに現れるのかは、後ほど定容過程の議論で明らかになります。結論を先に言えば、定容過程では $W = 0$ なので $dU = \delta Q_v = nC_v \, dT$ となり、$C_v$ こそが内部エネルギーの温度変化率を表す量だからです。
マイヤーの関係式
理想気体では、定圧モル比熱 $C_p$ と定容モル比熱 $C_v$ の間に次の関係が成り立ちます。
$$ \begin{equation} C_p – C_v = R \end{equation} $$
ここで $R = 8.314 \, \text{J/(mol} \cdot \text{K)}$ は気体定数です。この式は「定圧で気体を加熱すると、温度上昇に使われるエネルギー（$nC_v \Delta T$）に加えて、膨張仕事（$nR\Delta T$）の分だけ余分に熱を供給しなければならない」ことを表しています。つまり、$C_p$ が $C_v$ より常に大きいのは、定圧では膨張仕事の分だけ追加の熱が必要だからです。マイヤーの関係式の詳しい導出は比熱と比熱比で解説しています。
比熱比
$C_p$ と $C_v$ の比を比熱比（断熱指数）$\gamma$ と定義します。
$$ \gamma = \frac{C_p}{C_v} $$
分子の種類によって自由度が異なるため、$\gamma$ の値も変わります。
| 分子の種類 | 自由度 | $C_v$ | $C_p$ | $\gamma$ |
|---|---|---|---|---|
| 単原子（He, Ar） | 3（並進のみ） | $\frac{3}{2}R$ | $\frac{5}{2}R$ | $\frac{5}{3} \approx 1.67$ |
| 2原子（N$_2$, O$_2$） | 5（並進 + 回転） | $\frac{5}{2}R$ | $\frac{7}{2}R$ | $\frac{7}{5} = 1.4$ |
| 多原子（CO$_2$, H$_2$O） | 6以上 | $3R$ 以上 | $4R$ 以上 | $\frac{4}{3} \approx 1.33$ |
$\gamma$ は断熱過程の解析で中心的な役割を果たします。
内部エネルギーの温度依存性と比熱の関係を整理したところで、いよいよ第1法則を具体的な熱力学過程に適用していきましょう。「何が固定されて、何が変化するのか」を一つずつ確認しながら進めます。
各種熱力学過程への適用
熱力学では、系がある状態から別の状態へ変化する際に、特定の条件（体積一定、圧力一定など）を課した過程を考えます。それぞれの過程で第1法則 $\Delta U = Q – W$ がどう簡略化されるかを見ていきましょう。各過程の詳細は熱力学の各種過程でも扱っています。
1. 定容過程（Isochoric Process, $V = \text{const}$）
物理的な状況
密閉された剛体容器（体積が変わらない容器）に入った気体を加熱する場面を想像してください。容器の壁は動かないので、気体は膨張できません。圧力鍋のような状況です（実際の圧力鍋には安全弁がありますが、それを無視した理想的な状況です）。
何が固定で、何が変化するか
- 固定: 体積 $V$（容器の壁が動かないため）
- 変化: 温度 $T$ と圧力 $p$（加熱によって上昇）
第1法則の適用
体積が変化しないので $dV = 0$ であり、仕事は $W = \int p \, dV = 0$ です。第1法則 $\Delta U = Q – W$ に代入すると、次のようになります。
$$ \begin{align} \Delta U &= Q_v – 0 \\ Q_v &= \Delta U = nC_v \Delta T \end{align} $$
この結果は明快です。加えた熱が全て内部エネルギーの増加に使われます。仕事に逃げるエネルギーがないので、投入した熱量がそのまま温度上昇に反映されます。だからこそ、この過程での比熱を「定容比熱 $C_v$」と呼ぶのです。
家計簿のアナロジーで言えば、「支出（仕事）がゼロなので、収入（熱）がそのまま残高増加（内部エネルギー増加）になる」という状況です。
定容過程では体積が変わらないため、p-V図上では垂直な直線として表されます。曲線の下の面積がゼロなので、仕事もゼロであることが図からも確認できます。
次に、体積ではなく圧力を一定に保つとどうなるでしょうか。
2. 定圧過程（Isobaric Process, $p = \text{const}$）
物理的な状況
大気圧のもとで自由に動けるピストン付きシリンダーに入った気体を加熱する場面を考えます。ピストンの上には一定の重りが載っており、気体にかかる圧力は常に一定です。気体を加熱すると、気体は膨張してピストンを押し上げます。
何が固定で、何が変化するか
- 固定: 圧力 $p$（ピストン上の重りで一定に保たれる）
- 変化: 体積 $V$ と温度 $T$（加熱によって増大）
第1法則の適用
圧力が一定なので、仕事は簡単に計算できます。
$$ W = \int_{V_1}^{V_2} p \, dV = p(V_2 – V_1) = p\Delta V $$
理想気体の状態方程式 $pV = nRT$ を使うと、$p\Delta V = nR\Delta T$ と書き換えられます。
第1法則 $\Delta U = Q_p – W$ に代入します。
$$ Q_p = \Delta U + W = nC_v\Delta T + nR\Delta T $$
ここで、マイヤーの関係式 $C_p = C_v + R$ を思い出すと、右辺を次のようにまとめられます。
$$ \begin{equation} Q_p = n(C_v + R)\Delta T = nC_p\Delta T \end{equation} $$
定圧過程では、加えた熱の一部が温度上昇（内部エネルギー増加: $nC_v\Delta T$）に、残りが膨張仕事（$nR\Delta T$）に使われます。そのため、同じ温度変化を起こすのに、定容過程（$Q_v = nC_v\Delta T$）よりも多くの熱が必要です。定圧比熱 $C_p$ が定容比熱 $C_v$ よりも常に大きい理由が、ここからも理解できます。
家計簿のアナロジーでは、「収入（熱）の一部が支出（膨張仕事）に回るので、残高増加（内部エネルギー増加）は収入より少ない」という状況です。
p-V図上では、定圧過程は水平な直線として表されます。仕事 $W = p\Delta V$ は、この水平線と体積軸に囲まれた長方形の面積として読み取れます。
定圧過程で現れた $Q_p = nC_p\Delta T$ という関係は非常に便利ですが、これをさらに見通しよく表現する量として「エンタルピー」が導入されます。エンタルピーについては、各種過程を一通り見た後に改めて議論しましょう。
次に、温度を一定に保つ過程を見ていきます。
3. 等温過程（Isothermal Process, $T = \text{const}$）
物理的な状況
気体をゆっくりと膨張（または圧縮）させながら、外部の熱浴（温度が変わらないほど十分に大きな熱源）と常に熱的に接触させて、系の温度を一定に保つ過程です。たとえば、大きな恒温槽に浸したシリンダー内のガスを、ゆっくりピストンで押す（あるいは引く）状況に対応します。
何が固定で、何が変化するか
- 固定: 温度 $T$（熱浴との熱交換で一定に保たれる）
- 変化: 体積 $V$ と圧力 $p$（$pV = nRT = \text{const}$ の関係を保ちながら変化）
第1法則の適用
理想気体の内部エネルギーは温度のみの関数なので、$\Delta T = 0$ ならば $\Delta U = nC_v \Delta T = 0$ です。第1法則 $\Delta U = Q – W$ より、次のことがわかります。
$$ Q = W $$
つまり、等温過程では加えた熱が全て仕事に変換されます（膨張の場合）。内部エネルギーが変化しないので、投入したエネルギーは全て外界への仕事として出ていきます。逆に、等温圧縮の場合は、外界からされた仕事が全て熱として放出されます。
仕事を具体的に計算しましょう。理想気体の状態方程式 $p = nRT/V$ を $\delta W = p \, dV$ に代入して積分します。
$$ W = \int_{V_1}^{V_2} p \, dV = \int_{V_1}^{V_2} \frac{nRT}{V} dV $$
$T$ は定数なので積分の外に出せて、次の結果を得ます。
$$ \begin{equation} W = nRT \ln\frac{V_2}{V_1} \end{equation} $$
膨張（$V_2 > V_1$）のとき $W > 0$（系が仕事をする）、圧縮（$V_2 < V_1$）のとき $W < 0$（系が仕事をされる）となり、符号も物理的な直感と合います。
p-V図上で等温過程は双曲線（$pV = \text{const}$）として描かれます。この双曲線は「等温線」と呼ばれ、温度が高いほどp-V図上で上方に位置します。
ここまでの3つの過程では、いずれも外界との熱のやり取りがありました。では、系が外界と一切熱をやり取りしない場合にはどうなるのでしょうか。
4. 断熱過程（Adiabatic Process, $Q = 0$）
物理的な状況
系の周囲が完全な断熱壁で囲まれていて、外界との熱のやり取りが一切ない過程です。実際の現象では、音波の振動のように変化が非常に速い場合（熱が伝わる時間がないほど速い変化）に、良い近似で断熱過程とみなせます。ディーゼルエンジンの圧縮行程も、ピストンの動きが速いため、概ね断熱過程です。
何が固定で、何が変化するか
- 固定: 熱の出入り $Q = 0$（断熱壁で遮断）
- 変化: 温度 $T$、圧力 $p$、体積 $V$ の全てが変化
断熱過程では固定される熱力学的な状態量はなく、代わりに $Q = 0$ という過程の制約が課されます。この制約のもとで $T$、$p$、$V$ が互いにどのような関係で変化するかが、次に導くポアソンの関係式で記述されます。
第1法則の適用
$Q = 0$ を第1法則に代入すると、次のようになります。
$$ \begin{equation} dU = -\delta W \quad \Longrightarrow \quad nC_v \, dT = -p \, dV \end{equation} $$
この式は「断熱過程では、系がした仕事の分だけ内部エネルギーが減少する（温度が下がる）」ことを表しています。家計簿のアナロジーでは、「収入（熱）がゼロなので、支出（仕事）はすべて貯金（内部エネルギー）の取り崩しでまかなう」という状況です。
断熱膨張すると気体の温度が下がり、断熱圧縮すると温度が上がります。スプレー缶を使い続けるとガスが断熱膨張に近い形で噴出するため、缶が冷たくなるのはこの原理です。
ここまでで断熱過程での第1法則の基本的な帰結がわかりました。次に、$p$、$V$、$T$ の間の具体的な関係式（ポアソンの関係式）を導出しましょう。
ポアソンの関係式の導出
断熱過程で $p$、$V$、$T$ がどのように連動するかを、第1法則と理想気体の状態方程式から導きます。ゴールは、有名な $pV^\gamma = \text{const}$ という関係式を得ることです。
出発点
断熱過程での第1法則は $nC_v \, dT = -p \, dV$ です。ここに理想気体の状態方程式 $p = nRT/V$ を代入します。
$$ nC_v \, dT = -\frac{nRT}{V} \, dV $$
変数分離
両辺を $nT$ で割って、$T$ と $V$ を左辺と右辺に分離します。
$$ \frac{dT}{T} = -\frac{R}{C_v} \frac{dV}{V} $$
ここで、マイヤーの関係式 $R = C_p – C_v$ と比熱比 $\gamma = C_p/C_v$ を使うと、$R/C_v = (C_p – C_v)/C_v = \gamma – 1$ と書けます。したがって、
$$ \frac{dT}{T} = -(\gamma – 1) \frac{dV}{V} $$
積分
両辺を積分します。左辺は $T$ について、右辺は $V$ について積分します。
$$ \int \frac{dT}{T} = -(\gamma – 1) \int \frac{dV}{V} $$
$$ \ln T = -(\gamma – 1) \ln V + C_1 $$
ここで $C_1$ は積分定数です。
整理
右辺の対数をまとめます。$-(\gamma – 1)\ln V = \ln V^{-(\gamma-1)}$ なので、
$$ \ln T + (\gamma – 1)\ln V = C_1 $$
$$ \ln\left(T V^{\gamma – 1}\right) = C_1 $$
両辺の指数をとると、次の関係式が得られます。
$$ \begin{equation} TV^{\gamma – 1} = \text{const} \end{equation} $$
これがポアソンの関係式の温度-体積版です。
圧力-体積版への変換
理想気体の状態方程式 $T = pV/(nR)$ を $TV^{\gamma-1} = \text{const}$ に代入します。
$$ \frac{pV}{nR} \cdot V^{\gamma-1} = \text{const} $$
$nR$ は定数なので左辺にまとめると、
$$ \begin{equation} pV^\gamma = \text{const} \end{equation} $$
これが最もよく知られた形のポアソンの関係式です。断熱過程のp-V曲線は $p \propto V^{-\gamma}$ で表され、$\gamma > 1$ なので等温線（$p \propto V^{-1}$）よりも急な傾きを持ちます。これは物理的にも理解できます。断熱膨張では熱の補給がないので温度が下がり、それによって圧力がさらに下がるため、等温過程より急速に圧力が減少するのです。
断熱過程の仕事
断熱過程での仕事を求めましょう。$Q = 0$ なので $W = -\Delta U$ です。
$$ W = -\Delta U = -nC_v(T_2 – T_1) = nC_v(T_1 – T_2) $$
理想気体の状態方程式を使って温度を消去すると、次のようにも書けます。
$$ \begin{equation} W = \frac{p_1V_1 – p_2V_2}{\gamma – 1} \end{equation} $$
この式の導出を示しておきましょう。$nC_v = nR/(\gamma-1)$（マイヤーの関係式 $R = C_p – C_v$ と $\gamma = C_p/C_v$ から導けます）を使い、さらに $nRT_1 = p_1V_1$、$nRT_2 = p_2V_2$ を代入すると、
$$ W = \frac{nR(T_1 – T_2)}{\gamma – 1} = \frac{p_1V_1 – p_2V_2}{\gamma – 1} $$
が得られます。
ポアソンの関係式によって、断熱過程での$p$-$V$-$T$の連動が完全に記述されました。次に、定圧過程で自然に現れた量「エンタルピー」について、なぜこの量が便利なのかを掘り下げましょう。
各過程の比較まとめ
ここまでに学んだ4つの過程を一覧表にまとめます。
| 過程 | 条件 | $Q$ | $W$ | $\Delta U$ |
|---|---|---|---|---|
| 定容 | $V = \text{const}$ | $nC_v\Delta T$ | $0$ | $nC_v\Delta T$ |
| 定圧 | $p = \text{const}$ | $nC_p\Delta T$ | $p\Delta V = nR\Delta T$ | $nC_v\Delta T$ |
| 等温 | $T = \text{const}$ | $nRT\ln\frac{V_2}{V_1}$ | $nRT\ln\frac{V_2}{V_1}$ | $0$ |
| 断熱 | $Q = 0$ | $0$ | $nC_v(T_1 – T_2)$ | $nC_v\Delta T$ |
この表から、いくつかの重要な特徴が読み取れます。
まず、$\Delta U$ は定容・定圧・断熱のいずれでも $nC_v\Delta T$ です（等温過程では $\Delta T = 0$ なので $\Delta U = 0$）。これは内部エネルギーが状態量であり、同じ温度変化 $\Delta T$ なら経路によらず同じ変化をすることの反映です。
次に、同じ温度上昇 $\Delta T$ を実現するために必要な熱量は、定圧過程のほうが定容過程より多い（$C_p > C_v$）。定圧過程では膨張仕事の分だけ余計にエネルギーを供給する必要があるからです。
さらに、等温過程では $Q = W$、つまりエネルギーの全てが熱と仕事の間で行き来するだけで、内部エネルギーは変わりません。
この表を頭に入れた上で、定圧過程の取り扱いをさらに便利にする「エンタルピー」という量を導入しましょう。
エンタルピーの導入 — なぜこの量が便利なのか
動機: 定圧過程の計算を楽にしたい
工学や化学の実験の多くは大気圧下、すなわち定圧条件のもとで行われます。ビーカー内の化学反応、ボイラーでの蒸気生成、大気中での燃焼――これらは全て定圧過程です。
定圧過程で第1法則を適用すると、次のように $\Delta U$ と $p\Delta V$ を毎回分けて考える必要があります。
$$ Q_p = \Delta U + p\Delta V $$
もし、$U + pV$ という組み合わせをひとまとまりの量として定義しておけば、定圧過程での熱量が一発で書けて便利ではないでしょうか。これがエンタルピー $H$ の導入動機です。
エンタルピーの定義
$$ \begin{equation} H = U + pV \end{equation} $$
$H$ は内部エネルギー $U$ に圧力と体積の積 $pV$ を加えたもので、状態量です。$U$、$p$、$V$ がすべて状態量なので、それらの組み合わせである $H$ も状態量です。
定圧過程でのエンタルピーの威力
定圧過程（$p = \text{const}$）でのエンタルピー変化を考えます。
$$ \Delta H = \Delta U + \Delta(pV) = \Delta U + p\Delta V $$
右辺はまさに $Q_p$（定圧過程での熱量）そのものです。したがって、
$$ \begin{equation} Q_p = \Delta H \end{equation} $$
この式の美しさは、定圧過程で系に加えた熱量が、エンタルピーの変化そのものであることです。$\Delta U$ と $p\Delta V$ を分けて計算する必要がなくなり、エンタルピーだけ追えばよくなります。
理想気体の場合、$H$ も温度のみの関数となり（$H = U + pV = U + nRT$ で、$U$ も $nRT$ も温度のみの関数）、次のように書けます。
$$ \Delta H = nC_p \Delta T $$
化学の分野では「反応エンタルピー」「生成エンタルピー」「蒸発エンタルピー」など、定圧条件での熱量をエンタルピーの変化として記述するのが標準的です。エンタルピーの詳しい性質や応用は内部エネルギーとエンタルピーで解説しています。
エンタルピーという便利な量を手に入れたところで、ここまでの理論をPythonで可視化して、数式で得た結論を目で確認しましょう。
Pythonでの実装
p-V図上の各種過程の可視化
まず、4つの熱力学過程がp-V図上でどのような曲線を描くかを可視化します。それぞれの過程が同じ初期状態から出発したとき、圧力と体積がどのように変化するかを比較することが目的です。
import numpy as np
import matplotlib.pyplot as plt
# 理想気体のパラメータ
n = 1.0 # モル数
R = 8.314 # 気体定数 [J/(mol·K)]
gamma = 1.4 # 2原子分子の比熱比（空気を想定）
# 初期状態
T1 = 300 # 初期温度 [K]
p1 = 1e5 # 初期圧力 [Pa]（≒ 1 atm）
V1 = n * R * T1 / p1 # 初期体積 [m^3]
# 体積の範囲
V = np.linspace(V1 * 0.4, V1 * 3.5, 500)
fig, ax = plt.subplots(figsize=(10, 7))
# 等温過程: pV = nRT = const → p = nRT/V
p_isothermal = n * R * T1 / V
ax.plot(V * 1e3, p_isothermal / 1e3, 'b-', linewidth=2.5,
label=f'Isothermal (T = {T1} K)')
# 断熱過程: pV^γ = const → p = p1 * (V1/V)^γ
p_adiabatic = p1 * (V1 / V)**gamma
ax.plot(V * 1e3, p_adiabatic / 1e3, 'r-', linewidth=2.5,
label=f'Adiabatic (γ = {gamma})')
# 定圧過程: p = p1 = const
ax.axhline(y=p1 / 1e3, color='g', linewidth=2.5, linestyle='--',
label=f'Isobaric (p = {p1/1e3:.0f} kPa)')
# 定容過程: V = V1 = const
ax.axvline(x=V1 * 1e3, color='m', linewidth=2.5, linestyle='-.',
label=f'Isochoric (V = {V1*1e3:.2f} L)')
# 初期状態の点
ax.plot(V1 * 1e3, p1 / 1e3, 'ko', markersize=12, zorder=5,
label='Initial state')
ax.set_xlabel('Volume V [L]', fontsize=13)
ax.set_ylabel('Pressure p [kPa]', fontsize=13)
ax.set_title('Thermodynamic Processes on p-V Diagram', fontsize=15)
ax.legend(fontsize=11, loc='upper right')
ax.grid(True, alpha=0.3)
ax.set_xlim(V[0] * 1e3, V[-1] * 1e3)
ax.set_ylim(0, 300)
plt.tight_layout()
plt.show()
print(f"初期状態: T = {T1} K, p = {p1/1e3:.0f} kPa, V = {V1*1e3:.2f} L")
このp-V図から、4つの過程の特徴が視覚的に読み取れます。
-
等温線（青の実線）は双曲線 $pV = \text{const}$ の形をしており、体積が増えると圧力は緩やかに減少します。外部の熱浴から熱が供給されるため、温度（したがって内部エネルギー）は一定に保たれます。
-
断熱線（赤の実線）は等温線より急勾配です。$\gamma = 1.4 > 1$ なので、$p \propto V^{-1.4}$ は $p \propto V^{-1}$（等温線）より速く減衰します。断熱膨張では熱の補給がなく温度が下がるため、圧力もそれだけ余計に低下するのです。
-
定圧線（緑の破線）は水平で、この線と体積軸に挟まれた長方形の面積が仕事 $W = p\Delta V$ に対応します。
-
定容線（紫の一点鎖線）は垂直で、体積が変化しないので曲線下の面積（= 仕事）はゼロです。
等温膨張と断熱膨張の仕事の比較
次に、同じ体積変化（体積を2倍に膨張）をしたとき、等温過程と断熱過程で仕事がどう違うかを定量的に比較します。p-V図の曲線下の面積を塗りつぶすことで、仕事の大きさを視覚的に比べます。
import numpy as np
import matplotlib.pyplot as plt
# パラメータ
n = 1.0
R = 8.314
gamma = 1.4
Cv = R / (gamma - 1) # 定容モル比熱
T1 = 300 # 初期温度 [K]
p1 = 1e5 # 初期圧力 [Pa]
V1 = n * R * T1 / p1 # 初期体積
V2 = V1 * 2 # 最終体積（2倍に膨張）
# 体積の配列
V_range = np.linspace(V1, V2, 300)
# 等温膨張
p_iso = n * R * T1 / V_range
W_isothermal = n * R * T1 * np.log(V2 / V1)
# 断熱膨張
p_adi = p1 * (V1 / V_range)**gamma
T2_adi = T1 * (V1 / V2)**(gamma - 1)
W_adiabatic = n * Cv * (T1 - T2_adi)
fig, ax = plt.subplots(figsize=(10, 7))
# 等温過程の仕事（面積）
ax.fill_between(V_range * 1e3, 0, p_iso / 1e3, alpha=0.25, color='blue')
ax.plot(V_range * 1e3, p_iso / 1e3, 'b-', linewidth=2.5,
label=f'Isothermal: W = {W_isothermal:.1f} J')
# 断熱過程の仕事（面積）
ax.fill_between(V_range * 1e3, 0, p_adi / 1e3, alpha=0.25, color='red')
ax.plot(V_range * 1e3, p_adi / 1e3, 'r-', linewidth=2.5,
label=f'Adiabatic: W = {W_adiabatic:.1f} J')
# 始点と終点
ax.plot(V1 * 1e3, p1 / 1e3, 'ko', markersize=10)
ax.annotate('Initial state', xy=(V1 * 1e3, p1 / 1e3),
xytext=(V1 * 1e3 + 3, p1 / 1e3 + 15),
fontsize=10, arrowprops=dict(arrowstyle='->', color='black'))
ax.set_xlabel('Volume V [L]', fontsize=13)
ax.set_ylabel('Pressure p [kPa]', fontsize=13)
ax.set_title('Work Comparison: Isothermal vs Adiabatic Expansion (V → 2V)',
fontsize=14)
ax.legend(fontsize=12)
ax.grid(True, alpha=0.3)
plt.tight_layout()
plt.show()
# 数値結果
print("=" * 55)
print("膨張仕事の比較（体積を2倍に膨張）")
print("=" * 55)
print(f"等温膨張: W = {W_isothermal:.1f} J, Q = {W_isothermal:.1f} J, ΔU = 0 J")
print(f"断熱膨張: W = {W_adiabatic:.1f} J, Q = 0 J, ΔU = {-W_adiabatic:.1f} J")
print(f"断熱膨張後の温度: T₂ = {T2_adi:.1f} K (初期 {T1} K)")
print(f"仕事の比: W_iso / W_adi = {W_isothermal / W_adiabatic:.3f}")
このグラフと数値結果から、等温膨張と断熱膨張の物理的な違いが明確に読み取れます。
まず、等温膨張の仕事（青い面積）は断熱膨張の仕事（赤い面積）よりも大きいことがわかります。等温過程では外部の熱浴から熱が供給され続けるため、内部エネルギー（温度）が一定に保たれ、圧力の低下が穏やかです。その結果、p-V曲線が断熱線よりも上を通り、曲線下の面積（= 仕事）が大きくなります。
一方、断熱膨張では熱の供給がないため、仕事をした分だけ内部エネルギーが減少し、温度が下がります。温度が下がると圧力もより大きく下がるので、断熱線は等温線よりも急に下降します。その結果、曲線下の面積（= 仕事）は等温過程より小さくなります。
断熱膨張後の温度が初期温度（300 K）から約228 Kまで下がっていることも確認できます。この温度低下こそが、スプレー缶を使い続けると冷たくなる理由です。
エネルギー収支の棒グラフ比較
最後に、各過程でのエネルギー収支（$Q$、$W$、$\Delta U$）を棒グラフで並べて、第1法則 $\Delta U = Q – W$ が全ての過程で成り立っていることを視覚的に確認します。ここでは、同じ温度変化（$\Delta T = 100$ K の上昇）を実現する場合と、等温過程（体積を2倍に膨張）のケースを並べます。
import numpy as np
import matplotlib.pyplot as plt
# パラメータ
n = 1.0
R = 8.314
Cv = 5 / 2 * R # 2原子分子の定容比熱
Cp = 7 / 2 * R # 2原子分子の定圧比熱
gamma = Cp / Cv
T1 = 300 # 初期温度 [K]
dT = 100 # 温度変化 [K]（定容・定圧・断熱で使用）
# 各過程でのエネルギー収支 [J]
# 定容: V=const, W=0, Q=ΔU=nCvΔT
Q_isochoric = n * Cv * dT
W_isochoric = 0
dU_isochoric = n * Cv * dT
# 定圧: p=const, Q=nCpΔT, W=nRΔT, ΔU=nCvΔT
Q_isobaric = n * Cp * dT
W_isobaric = n * R * dT
dU_isobaric = n * Cv * dT
# 等温: T=const, ΔU=0, Q=W=nRT*ln(2)
Q_isothermal = n * R * T1 * np.log(2)
W_isothermal = n * R * T1 * np.log(2)
dU_isothermal = 0
# 断熱: Q=0, W=-ΔU, ΔU=nCvΔT（圧縮: 温度上昇）
Q_adiabatic = 0
W_adiabatic = -n * Cv * dT # 圧縮（外界が仕事をする）
dU_adiabatic = n * Cv * dT
processes = [
'Isochoric\n(V = const)',
'Isobaric\n(p = const)',
'Isothermal\n(T = const)',
'Adiabatic\n(Q = 0)'
]
Q_vals = [Q_isochoric, Q_isobaric, Q_isothermal, Q_adiabatic]
W_vals = [W_isochoric, W_isobaric, W_isothermal, W_adiabatic]
dU_vals = [dU_isochoric, dU_isobaric, dU_isothermal, dU_adiabatic]
x = np.arange(len(processes))
width = 0.25
fig, ax = plt.subplots(figsize=(12, 7))
bars1 = ax.bar(x - width, Q_vals, width, label='Q (Heat)', color='red', alpha=0.7)
bars2 = ax.bar(x, W_vals, width, label='W (Work)', color='blue', alpha=0.7)
bars3 = ax.bar(x + width, dU_vals, width, label='ΔU (Internal energy change)',
color='green', alpha=0.7)
ax.set_xlabel('Thermodynamic Process', fontsize=12)
ax.set_ylabel('Energy [J]', fontsize=12)
ax.set_title('First Law of Thermodynamics: ΔU = Q − W', fontsize=14)
ax.set_xticks(x)
ax.set_xticklabels(processes, fontsize=11)
ax.legend(fontsize=11, loc='upper left')
ax.grid(True, alpha=0.3, axis='y')
ax.axhline(y=0, color='k', linewidth=0.8)
plt.tight_layout()
plt.show()
# 数値の確認
print("=" * 60)
print("各過程のエネルギー収支 [J]")
print("=" * 60)
print(f"{'過程':<12} {'Q':>10} {'W':>10} {'ΔU':>10} {'Q-W':>10}")
print("-" * 60)
for i, proc in enumerate(['定容', '定圧', '等温', '断熱']):
check = Q_vals[i] - W_vals[i]
print(f"{proc:<12} {Q_vals[i]:>10.1f} {W_vals[i]:>10.1f} "
f"{dU_vals[i]:>10.1f} {check:>10.1f}")
この棒グラフから、第1法則の構造が一目でわかります。
定容過程（左端）では、仕事（青）がゼロで、熱（赤）と内部エネルギー変化（緑）が完全に一致しています。加えた熱が全て内部エネルギーに変わる、という定容過程の性質がグラフに表れています。
定圧過程では、熱（赤）が最も大きく、その一部が仕事（青）に、残りが内部エネルギー変化（緑）に分配されています。$Q = \Delta U + W$ の関係が視覚的に確認でき、$C_p > C_v$（定圧比熱が定容比熱より大きい）理由もここから読み取れます。
等温過程では、内部エネルギー変化（緑）がゼロで、熱（赤）と仕事（青）が完全に一致しています。加えた熱が全て仕事に変わる（温度を保つための膨張仕事にまわる）ことが見て取れます。
断熱過程（断熱圧縮のケース）では、熱（赤）がゼロで、外界からされた仕事（青の負の値）がそのまま内部エネルギー増加（緑）になっています。出力テーブルの最右列 $Q – W$ が全ての過程で $\Delta U$ と一致していることも確認でき、第1法則 $\Delta U = Q – W$ の正しさが数値的に裏付けられています。
断熱過程と等温過程のp-V曲線比較（比熱比の影響）
最後に、比熱比 $\gamma$ の値によって断熱線の形がどう変わるかを可視化します。$\gamma$ は分子の自由度（= 分子の種類）によって決まるので、気体の種類が違うと断熱過程の振る舞いも異なります。
import numpy as np
import matplotlib.pyplot as plt
# パラメータ
n = 1.0
R = 8.314
T1 = 300
p1 = 1e5
V1 = n * R * T1 / p1
V = np.linspace(V1 * 0.5, V1 * 3, 400)
fig, ax = plt.subplots(figsize=(10, 7))
# 等温線（基準）
p_iso = n * R * T1 / V
ax.plot(V * 1e3, p_iso / 1e3, 'k--', linewidth=2, label='Isothermal (T = 300 K)')
# 異なるγでの断熱線
gammas = [5/3, 7/5, 4/3]
colors = ['red', 'blue', 'green']
labels = ['Monoatomic (γ = 5/3)', 'Diatomic (γ = 7/5)', 'Polyatomic (γ = 4/3)']
for g, c, lab in zip(gammas, colors, labels):
p_adi = p1 * (V1 / V)**g
ax.plot(V * 1e3, p_adi / 1e3, color=c, linewidth=2, label=lab)
ax.plot(V1 * 1e3, p1 / 1e3, 'ko', markersize=10, zorder=5)
ax.set_xlabel('Volume V [L]', fontsize=13)
ax.set_ylabel('Pressure p [kPa]', fontsize=13)
ax.set_title('Effect of Heat Capacity Ratio γ on Adiabatic Curves', fontsize=14)
ax.legend(fontsize=11)
ax.grid(True, alpha=0.3)
ax.set_xlim(V[0] * 1e3, V[-1] * 1e3)
ax.set_ylim(0, 350)
plt.tight_layout()
plt.show()
# 体積2倍に膨張した場合の最終温度の比較
print("=" * 55)
print("体積を2倍に断熱膨張した場合の最終温度")
print("=" * 55)
for g, lab in zip(gammas, ['単原子', '2原子', '多原子']):
T2 = T1 * (V1 / (V1 * 2))**(g - 1)
print(f"{lab} (γ = {g:.4f}): T₂ = {T2:.1f} K (ΔT = {T2 - T1:.1f} K)")
このグラフから、比熱比 $\gamma$ が断熱過程に与える影響が明確にわかります。
$\gamma$ が大きいほど、断熱線はより急勾配になります。単原子気体（$\gamma = 5/3$、赤線）の断熱線が最も急で、多原子気体（$\gamma = 4/3$、緑線）の断熱線が最も緩やかです。全ての断熱線は等温線（黒の破線）よりも急であることも確認できます。
これは物理的に理解できます。$\gamma$ が大きいということは $C_v$ が小さい（自由度が少ない）ことを意味します。$C_v$ が小さいと、同じ仕事をしても温度がより大きく変化する（$\Delta T = -W/(nC_v)$）ため、圧力の変化もより大きくなるのです。数値結果を見ると、体積を2倍に断熱膨張した場合、単原子気体では約113 Kの温度低下が起こるのに対し、多原子気体では約62 Kの温度低下にとどまります。
等温線（黒の破線）は全ての断熱線よりも上にあることにも注目してください。等温膨張では外部から熱が供給され続けるので温度（= 圧力）が維持され、曲線が上に留まります。断熱膨張では熱の供給がないので温度が下がり、曲線がより下に沈みます。この差が、先に見た「等温膨張の仕事 > 断熱膨張の仕事」という結果の視覚的な根拠です。
第1法則の物理的意義 — エネルギー保存の普遍性
ここまでの議論を一歩引いて俯瞰してみましょう。熱力学第1法則 $\Delta U = Q – W$ は、突き詰めればエネルギー保存則の一つの表現です。力学では運動エネルギーとポテンシャルエネルギーの和が保存される保存則を学びますが、現実の世界では摩擦や空気抵抗によって力学的エネルギーは減少するように見えます。
しかし、摩擦で「失われた」エネルギーは実は熱に変わっています。ジュールの実験（1840年代）は、力学的仕事と熱が本質的に同じ量であることを定量的に示しました。第1法則は、この知見を定式化したものであり、エネルギーの形態（力学的エネルギー、熱、化学エネルギー、電磁気エネルギーなど）がどう変わっても、総量は変わらないという普遍的な原理を表しています。
ただし、第1法則だけでは答えられない問いがあります。「エネルギーの総量は保存されるとして、あらゆるエネルギー変換は実際に起こりうるのか？」という問いです。たとえば、第1法則は「熱を100%仕事に変える」ことを禁じていません。$Q = W$、$\Delta U = 0$ とすれば帳尻は合います。しかし、現実にはそのような都合の良い変換は不可能です。この「エネルギー変換の方向性」の問題に答えるのが熱力学第2法則であり、そこで中心的な役割を果たすのがエントロピーです。
まとめ
本記事では、熱力学第1法則の数学的定式化から、各種過程への適用、ポアソンの関係式の導出、エンタルピーの導入動機までを一貫して解説しました。
-
熱力学第1法則 $\Delta U = Q – W$ は、エネルギー保存則の熱力学的表現です。内部エネルギーの変化（状態量）が、熱と仕事（いずれも経路量）の差として書けるという構造を持ちます
-
状態量と経路量の区別は第1法則の理解に不可欠です。山の標高（経路によらない）と歩行距離（経路による）のアナロジーで直感的に捉えられます
-
定容過程（$V = \text{const}$）では $W = 0$ となり、加えた熱が全て内部エネルギーに変わります（$Q_v = nC_v\Delta T$）
-
定圧過程（$p = \text{const}$）では加えた熱の一部が膨張仕事に使われます（$Q_p = nC_p\Delta T$）。エンタルピー $H = U + pV$ を導入すると $Q_p = \Delta H$ と簡潔に書けます
-
等温過程（$T = \text{const}$）では内部エネルギーが変化せず、$Q = W = nRT\ln(V_2/V_1)$ となります
-
断熱過程（$Q = 0$）では仕事の分だけ内部エネルギーが変化し、ポアソンの関係式 $pV^\gamma = \text{const}$ が導かれます
-
p-V図を使うと、仕事の大きさを曲線下の面積として視覚的に比較でき、等温膨張の仕事が断熱膨張より大きい理由も直感的に理解できます
第1法則はエネルギーの「収支」を教えてくれますが、エネルギー変換の「方向」については何も語りません。次のステップとして、エネルギー変換の方向性を支配する第2法則と、その中心概念であるエントロピーを学びましょう。
次のステップとして、以下の記事も参考にしてください。
- カルノーサイクル — 各過程を組み合わせた理想的な熱機関。第1法則と第2法則の両方が活きる応用例です
- エントロピー — エネルギー変換の方向性を支配する第2法則の中心概念です
- 自由エネルギー — エンタルピーとエントロピーを組み合わせた、化学反応の自発性を判定する量です
- 熱伝導（フーリエの法則） — 第1法則における「熱」が具体的にどのように伝わるかを記述する法則です