熱力学量に関する方程式
熱力学 において 、 基本的な熱力学関係 は、4つの重要な熱力学的量が、実験的に制御および測定できる変数にどのように依存するかを示す4つの基本方程式です。したがって、それらは本質的に状態方程式であり、基本方程式を使用すると、実験データを使用して、 G ( ギブズ自由エネルギー )や H ( エンタルピー )などの必要な量を決定することができます。 [1] この関係は、通常、 エントロピー の微視的変化と、熱平衡にある 閉鎖系の 体積 の観点から、 内部エネルギー の微視的変化として 次のように表現されます。
d
U
=
T
d
S
−
P
d
V
{\displaystyle \mathrm {d} U=T\,\mathrm {d} S-P\,\mathrm {d} V\,}
ここで、 U は 内部エネルギー 、 T は 絶対温度 、 S は エントロピー 、 P は 圧力 、 Vは 体積 です 。
これは基本的な熱力学関係の1つの表現に過ぎません。これは、異なる変数(例えば、熱力学的ポテンシャル )を使用して、他の方法で表現することもできます。たとえば、基本的な関係は、 エンタルピー H で 次のよう
に表現できます。
d
H
=
T
d
S
+
V
d
P
{\displaystyle \mathrm {d} H=T\,\mathrm {d} S+V\,\mathrm {d} P\,}
ヘルムホルツ自由エネルギー F に関して 、
d
F
=
−
S
d
T
−
P
d
V
{\displaystyle \mathrm {d} F=-S\,\mathrm {d} T-P\,\mathrm {d} V\,}
ギブスの自由エネルギー G に関しては 、
d
G
=
−
S
d
T
+
V
d
P
{\displaystyle \mathrm {d} G=-S\,\mathrm {d} T+V\,\mathrm {d} P\,}
。
熱力学の第一法則と第二法則
熱力学の第一法則は 次 のように述べています。
d
U
=
δ
Q
−
δ
W
{\displaystyle \mathrm {d} U=\delta Q-\delta W\,}
ここで 、およびは それぞれ、システムの周囲からシステムに供給される微量の熱と、システムが周囲に対して行う仕事です。
δ
Q
{\displaystyle \delta Q}
δ
W
{\displaystyle \delta W}
熱力学の第二法則 によれば、 可逆プロセスは次のようになります。
d
S
=
δ
Q
T
{\displaystyle \mathrm {d} S={\frac {\delta Q}{T}}\,}
したがって:
δ
Q
=
T
d
S
{\displaystyle \delta Q=T\,\mathrm {d} S\,}
これを第一法則に代入すると次のようになります。
d
U
=
T
d
S
−
δ
W
{\displaystyle \mathrm {d} U=T\,\mathrm {d} S-\delta W\,}
システムが周囲に対して行う
可逆的な 圧力-体積仕事 とすると、
δ
W
{\displaystyle \delta W}
δ
W
=
P
d
V
{\displaystyle \delta W\ =P\mathrm {d} V\,}
我々は持っています:
d
U
=
T
d
S
−
P
d
V
{\displaystyle \mathrm {d} U=T\,\mathrm {d} S-P\,\mathrm {d} V\,}
この式は可逆的な変化の場合に導出されたものです。しかし、 U 、 S 、 Vは 熱力学的プロセスの初期状態と最終状態にのみ依存する熱力学的 状態関数 であるため、上記の関係は非可逆的な変化にも当てはまります。 均一な温度と圧力のシステムにおける組成、つまり化学成分の量も、例えば化学反応により変化する可能性がある場合、基本的な熱力学関係は次のように一般化されます。
n
i
{\displaystyle n_{i}}
d
U
=
T
d
S
−
P
d
V
+
∑
i
μ
i
d
n
i
{\displaystyle \mathrm {d} U=T\,\mathrm {d} S-P\,\mathrm {d} V\ +\sum _{i}\mu _{i}\,\mathrm {d} n_{i}\,}
は 、 型の粒子に対応する 化学ポテンシャル です 。
μ
i
{\displaystyle \mu _{i}}
i
{\displaystyle i}
システムが体積だけでなく変化できる外部パラメータを持つ場合、基本的な熱力学関係は次のように一般化される。
d
U
=
T
d
S
+
∑
j
X
j
d
x
j
+
∑
i
μ
i
d
n
i
{\displaystyle \mathrm {d} U=T\,\mathrm {d} S+\sum _{j}X_{j}\,\mathrm {d} x_{j}+\sum _{i}\mu _{i}\,\mathrm {d} n_{i}\,}
ここで、 は 外部パラメータに対応する 一般化された力 です 。(圧力で使用される負の符号は珍しく、圧力は体積を減少させる傾向がある圧縮応力を表すために生じます。その他の一般化された力は共役変位を増加させる傾向があります。)
X
j
{\displaystyle X_{j}}
x
j
{\displaystyle x_{j}}
統計力学との関係
基本的な熱力学関係と統計力学原理は相互に導き出すことができます。
統計力学原理からの導出
上記の導出では、熱力学の第一法則と第二法則を使用しています。熱力学の第一法則は本質的に 熱 の定義であり、つまり熱はシステムの外部パラメータの変化によって引き起こされないシステムの内部エネルギーの変化です。
しかし、熱力学の第二法則はエントロピーを定義する関係ではありません。ある量のエネルギーを含む孤立したシステムのエントロピーの基本的な定義は次の とおりです。
E
{\displaystyle E}
S
=
k
log
[
Ω
(
E
)
]
{\displaystyle S=k\log \left[\Omega \left(E\right)\right]\,}
ここで、 は と の 間の小さな間隔にある量子状態の数です 。ここでは、 固定されたマクロ的に小さなエネルギー間隔があります。厳密に言えば、これはエントロピーが の選択に依存することを意味します 。ただし、熱力学的極限 (つまり、無限に大きいシステム サイズの極限) では、比エントロピー (単位体積または単位質量あたりのエントロピー) は に依存しません 。したがって、システムのエネルギーがサイズ の何らかの間隔にあることがわかっている場合、エントロピーはシステムが正確にどの量子状態にあるかについての不確実性の尺度です 。
Ω
(
E
)
{\displaystyle \Omega \left(E\right)}
E
{\displaystyle E}
E
+
δ
E
{\displaystyle E+\delta E}
δ
E
{\displaystyle \delta E}
δ
E
{\displaystyle \delta E}
δ
E
{\displaystyle \delta E}
δ
E
{\displaystyle \delta E}
したがって、第一原理から基本的な熱力学関係を導くことは、上記のエントロピーの定義が可逆プロセスに対して次のことを意味することを証明することに相当する。
d
S
=
δ
Q
T
{\displaystyle dS={\frac {\delta Q}{T}}}
統計力学 の基本的な仮定は、特定のエネルギーにおける すべての状態が等しく起こり得るというものです 。これにより、関心のあるすべての熱力学的量を抽出できます。温度は次のように定義されます。
Ω
(
E
)
{\displaystyle \Omega \left(E\right)}
1
k
T
≡
β
≡
d
log
[
Ω
(
E
)
]
d
E
{\displaystyle {\frac {1}{kT}}\equiv \beta \equiv {\frac {d\log \left[\Omega \left(E\right)\right]}{dE}}\,}
この定義は、ミクロカノニカル集団 から導き出すことができます。ミクロカノニカル集団 は、粒子の数が一定で体積が一定であり、環境とエネルギーを交換しないシステムです。システムには、変更可能な外部パラメータ x があるとします。一般に、 システムの エネルギー 固有状態は xに依存します。量子力学の 断熱定理 によれば 、システムのハミルトニアンの変化が無限に遅い限界では、システムは同じエネルギー固有状態にとどまり、現在のエネルギー固有状態のエネルギー変化に応じてエネルギーが変化します。
外部パラメータ x に対応する一般化された力 X は、 x が量 dx だけ増加した場合にシステムが実行する仕事となる ように定義されます 。たとえば、 x が体積の場合、 X は圧力です。エネルギー固有状態にあることがわかっているシステムの一般化された力は、次 のように表されます。
X
d
x
{\displaystyle Xdx}
E
r
{\displaystyle E_{r}}
X
=
−
d
E
r
d
x
{\displaystyle X=-{\frac {dE_{r}}{dx}}}
システムは の範囲内の任意のエネルギー固有状態をとることができるため 、システムの一般化された力を上記の式の期待値として定義します。
δ
E
{\displaystyle \delta E}
X
=
−
⟨
d
E
r
d
x
⟩
{\displaystyle X=-\left\langle {\frac {dE_{r}}{dx}}\right\rangle \,}
平均を評価するには、 から までの範囲内で の値を持つ固有状態の数を数えることによってエネルギー固有状態を分割します 。この数を と呼ぶと 、次の式が得られます。
Ω
(
E
)
{\displaystyle \Omega (E)}
d
E
r
d
x
{\displaystyle {\frac {dE_{r}}{dx}}}
Y
{\displaystyle Y}
Y
+
δ
Y
{\displaystyle Y+\delta Y}
Ω
Y
(
E
)
{\displaystyle \Omega _{Y}\left(E\right)}
Ω
(
E
)
=
∑
Y
Ω
Y
(
E
)
{\displaystyle \Omega (E)=\sum _{Y}\Omega _{Y}(E)\,}
一般化された力を定義する平均は次のように記述できます。
X
=
−
1
Ω
(
E
)
∑
Y
Y
Ω
Y
(
E
)
{\displaystyle X=-{\frac {1}{\Omega (E)}}\sum _{Y}Y\Omega _{Y}(E)\,}
これを、一定のエネルギー E での x に関するエントロピーの微分と次のように関連付けることができます。 x を x + dx に変更するとします 。すると は 変化します。これは、エネルギー固有状態が x に依存し、エネルギー固有状態が と の間の範囲に出入りするからです。 が と の間の範囲内にある エネルギー固有状態に再び注目してみましょう。これらのエネルギー固有状態はエネルギーが Y dx 増加するため、 E − Y dx から E の 範囲にあるすべてのエネルギー固有状態は、 E の下から E の上に移動 し ます 。
Ω
(
E
)
{\displaystyle \Omega \left(E\right)}
E
{\displaystyle E}
E
+
δ
E
{\displaystyle E+\delta E}
d
E
r
d
x
{\displaystyle {\frac {dE_{r}}{dx}}}
Y
{\displaystyle Y}
Y
+
δ
Y
{\displaystyle Y+\delta Y}
N
Y
(
E
)
=
Ω
Y
(
E
)
δ
E
Y
d
x
{\displaystyle N_{Y}(E)={\frac {\Omega _{Y}(E)}{\delta E}}Y\,dx}
そのようなエネルギー固有状態。 の場合、これらすべてのエネルギー固有状態は と の間の範囲に移動し 、 の増加に寄与します。下から 上に 移動するエネルギー固有状態の数は 、もちろん で与えられます 。その差
Y
d
x
≤
δ
E
{\displaystyle Ydx\leq \delta E}
E
{\displaystyle E}
E
+
δ
E
{\displaystyle E+\delta E}
Ω
{\displaystyle \Omega }
E
+
δ
E
{\displaystyle E+\delta E}
E
+
δ
E
{\displaystyle E+\delta E}
N
Y
(
E
+
δ
E
)
{\displaystyle N_{Y}\left(E+\delta E\right)}
N
Y
(
E
)
−
N
Y
(
E
+
δ
E
)
{\displaystyle N_{Y}(E)-N_{Y}(E+\delta E)\,}
したがって、 の増加に対する純粋な寄与となります 。Y dx が より大きい場合、の下から の 上へ 移動するエネルギー固有状態が存在することに注意してください 。これらは と の両方でカウントされる ため、その場合でも上記の式は有効です。
Ω
{\displaystyle \Omega }
δ
E
{\displaystyle \delta E}
E
{\displaystyle E}
E
+
δ
E
{\displaystyle E+\delta E}
N
Y
(
E
)
{\displaystyle N_{Y}(E)}
N
Y
(
E
+
δ
E
)
{\displaystyle N_{Y}(E+\delta E)}
上記の式を E に関する微分として表し、Y について合計すると、次の式が得られます。
(
∂
Ω
∂
x
)
E
=
−
∑
Y
Y
(
∂
Ω
Y
∂
E
)
x
=
(
∂
(
Ω
X
)
∂
E
)
x
{\displaystyle \left({\frac {\partial \Omega }{\partial x}}\right)_{E}=-\sum _{Y}Y\left({\frac {\partial \Omega _{Y}}{\partial E}}\right)_{x}=\left({\frac {\partial (\Omega X)}{\partial E}}\right)_{x}\,}
したがって、 のx に関する 対数微分は 次のように表されます。
Ω
{\displaystyle \Omega }
(
∂
log
(
Ω
)
∂
x
)
E
=
β
X
+
(
∂
X
∂
E
)
x
{\displaystyle \left({\frac {\partial \log \left(\Omega \right)}{\partial x}}\right)_{E}=\beta X+\left({\frac {\partial X}{\partial E}}\right)_{x}\,}
最初の項は強度的であり、システムのサイズに応じて変化しません。対照的に、最後の項はシステムサイズの逆数に応じて変化し、熱力学的限界では消えます。したがって、次のことが分かりました。
(
∂
S
∂
x
)
E
=
X
T
{\displaystyle \left({\frac {\partial S}{\partial x}}\right)_{E}={\frac {X}{T}}\,}
これと
(
∂
S
∂
E
)
x
=
1
T
{\displaystyle \left({\frac {\partial S}{\partial E}}\right)_{x}={\frac {1}{T}}\,}
与えるもの:
d
S
=
(
∂
S
∂
E
)
x
d
E
+
(
∂
S
∂
x
)
E
d
x
=
d
E
T
+
X
T
d
x
{\displaystyle dS=\left({\frac {\partial S}{\partial E}}\right)_{x}\,dE+\left({\frac {\partial S}{\partial x}}\right)_{E}\,dx={\frac {dE}{T}}+{\frac {X}{T}}\,dx\,}
これは次のように書くことができます。
d
E
=
T
d
S
−
X
d
x
{\displaystyle dE=T\,dS-X\,dx}
基本的な熱力学関係からの統計力学原理の導出
基本的な熱力学関係は、以下の3つの公理と合わせて示されている [2]
確率密度関数は、アンサンブル パラメーターとランダム変数の何らかの関数に比例します。 熱力学的状態関数は、ランダム変数の集合平均によって記述されます。 ギブスのエントロピー公式 によって定義されるエントロピーは、 古典熱力学 で定義されるエントロピーと一致します 。
等先験的確率の公理なしに統計力学の理論を構築するには十分です。
例えば、 ボルツマン分布を導くためには、ミクロ状態 i の確率密度が を満たすと仮定する 。したがって、正規化係数(分配関数)は
Pr
(
i
)
∝
f
(
E
i
,
T
)
{\textstyle \Pr(i)\propto f(E_{i},T)}
Z
=
∑
i
f
(
E
i
,
T
)
.
{\displaystyle Z=\sum _{i}f(E_{i},T).}
したがって、エントロピーは次のように表される。
S
=
k
B
∑
i
f
(
E
i
,
T
)
Z
log
(
f
(
E
i
,
T
)
Z
)
.
{\displaystyle S=k_{B}\sum _{i}{\frac {f(E_{i},T)}{Z}}\log \left({\frac {f(E_{i},T)}{Z}}\right).}
系の体積を一定に保ちながら
温度 Tを dT だけ変化させると、エントロピーの変化は
d
S
=
(
∂
S
∂
T
)
V
d
T
{\displaystyle dS=\left({\frac {\partial S}{\partial T}}\right)_{V}dT}
どこ
(
∂
S
∂
T
)
V
=
−
k
B
∑
i
Z
⋅
∂
f
(
E
i
,
T
)
∂
T
⋅
log
f
(
E
i
,
T
)
−
∂
Z
∂
T
⋅
f
(
E
i
,
T
)
⋅
log
f
(
E
i
,
T
)
Z
2
{\displaystyle \left({\frac {\partial S}{\partial T}}\right)_{V}=-k_{B}\sum _{i}{\frac {Z\cdot {\frac {\partial f(E_{i},T)}{\partial T}}\cdot \log f(E_{i},T)-{\frac {\partial Z}{\partial T}}\cdot f(E_{i},T)\cdot \log f(E_{i},T)}{Z^{2}}}}
=
−
k
B
∑
i
∂
∂
T
(
f
(
E
i
,
T
)
Z
)
⋅
log
f
(
E
i
,
T
)
{\displaystyle =-k_{B}\sum _{i}{\frac {\partial }{\partial T}}\left({\frac {f(E_{i},T)}{Z}}\right)\cdot \log f(E_{i},T)}
それを考慮すると
⟨
E
⟩
=
∑
i
f
(
E
i
,
T
)
Z
⋅
E
i
{\displaystyle \left\langle E\right\rangle =\sum _{i}{\frac {f(E_{i},T)}{Z}}\cdot E_{i}}
我々は持っています
d
⟨
E
⟩
=
∑
i
∂
∂
T
(
f
(
E
i
,
T
)
Z
)
⋅
E
i
⋅
d
T
{\displaystyle d\left\langle E\right\rangle =\sum _{i}{\frac {\partial }{\partial T}}\left({\frac {f(E_{i},T)}{Z}}\right)\cdot E_{i}\cdot dT}
基本的な熱力学関係から、
−
d
S
k
B
+
d
⟨
E
⟩
k
B
T
+
P
k
B
T
d
V
=
0
{\displaystyle -{\frac {dS}{k_{B}}}+{\frac {d\left\langle E\right\rangle }{k_{B}T}}+{\frac {P}{k_{B}T}}dV=0}
Tを摂動させる際に Vを 一定に 保ったので 、次の式が得られます 。上記の式を組み合わせると、次の式が得られます。
d
V
=
0
{\textstyle dV=0}
∑
i
∂
∂
T
(
f
(
E
i
,
T
)
Z
)
⋅
[
log
f
(
E
i
,
T
)
+
E
i
k
B
T
]
⋅
d
T
=
0
{\displaystyle \sum _{i}{\frac {\partial }{\partial T}}\left({\frac {f(E_{i},T)}{Z}}\right)\cdot \left[\log f(E_{i},T)+{\frac {E_{i}}{k_{B}T}}\right]\cdot dT=0}
物理法則は普遍的であるべきであり、つまり、上記の式は任意のシステムに対して成立するはずであり、これが起こる唯一の方法は
log
f
(
E
i
,
T
)
+
E
i
k
B
T
=
0
{\displaystyle \log f(E_{i},T)+{\frac {E_{i}}{k_{B}T}}=0}
つまり
f
(
E
i
,
T
)
=
exp
(
−
E
i
k
B
T
)
.
{\displaystyle f(E_{i},T)=\exp \left(-{\frac {E_{i}}{k_{B}T}}\right).}
上記の形式論における第3公準は次のように置き換えることができることが示されている: [3]
無限温度では、すべてのミクロ状態は同じ確率を持ちます。
しかし、数学的な導出ははるかに複雑になります。
参考文献
^ 「基本方程式の微分形式」。Chemistry LibreTexts 。2013年10月2日。
^ Gao, Xiang; Gallicchio, Emilio; Roitberg, Adrian (2019). 「一般化ボルツマン分布は、ギブス・シャノンエントロピーが熱力学的エントロピーに等しい唯一の分布である」。The Journal of Chemical Physics . 151 (3): 034113. arXiv : 1903.02121 . Bibcode :2019JChPh.151c4113G. doi :10.1063/1.5111333. PMID 31325924. S2CID 118981017.
^ Gao, Xiang (2022年3月). 「アンサンブル理論の数学」. Results in Physics . 34 : 105230. arXiv : 2006.00485 . Bibcode :2022ResPh..3405230G. doi : 10.1016/j.rinp.2022.105230 . S2CID 221978379.
外部リンク