設問(a) 式(2)の導出
使う関係と成立条件
熱力学第一法則と第二法則を合わせた d U = T d S − p d V \mathrm{d}U=T\,\mathrm{d}S-p\,\mathrm{d}V d U = T d S − p d V は、組成が変わらない閉じた系で、体積変化以外の仕事がなく、平衡状態を結ぶ可逆変化について成り立つ。U , S , V U,S,V U , S , V はいずれも状態量なので、この式は経路によらず、状態 ( T , V ) (T,V) ( T , V ) の近傍での変数間の関係を表す。変数を ( S , V ) (S,V) ( S , V ) から ( T , V ) (T,V) ( T , V ) に取り替えるには、T T T を独立変数にするルジャンドル変換 A = U − T S A=U-TS A = U − T S を使う。独立変数を ( T , V ) (T,V) ( T , V ) とすると、d A \mathrm{d}A d A の係数が − S -S − S と − p -p − p になる。
完全微分の係数の間には、f ( x , y ) f(x,y) f ( x , y ) について ∂ 2 f / ∂ x ∂ y = ∂ 2 f / ∂ y ∂ x \partial^2 f/\partial x\partial y=\partial^2 f/\partial y\partial x ∂ 2 f / ∂ x ∂ y = ∂ 2 f / ∂ y ∂ x が成り立つ(2階偏導関数が連続であれば微分の順序は入れ替えられる)。この性質を A ( T , V ) A(T,V) A ( T , V ) に使うと、直接測定しにくい ( ∂ S / ∂ V ) T (\partial S/\partial V)_T ( ∂ S / ∂ V ) T が、状態方程式から計算できる ( ∂ p / ∂ T ) V (\partial p/\partial T)_V ( ∂ p / ∂ T ) V に置き換わる。これが Maxwell の関係式である。
問題のどこが合図か
式(2)の左辺は「温度一定で体積を変えたときの内部エネルギーの変化」であり、状態方程式 p ( T , V ) p(T,V) p ( T , V ) だけから内部エネルギーの体積依存性を決める式になっている。問題文が A A A を定義していることは、d A = − S d T − p d V \mathrm{d}A=-S\,\mathrm{d}T-p\,\mathrm{d}V d A = − S d T − p d V から出発せよという合図と読める。
Joule の法則との整合と符号の誤り
理想気体 p = n R T / V p=nRT/V p = n R T / V では T ( ∂ p / ∂ T ) V = n R T / V = p T(\partial p/\partial T)_V=nRT/V=p T ( ∂ p / ∂ T ) V = n R T / V = p となり、右辺は − p + p = 0 -p+p=0 − p + p = 0 になる。これは Joule の法則(理想気体の内部エネルギーは体積によらない)に一致する。
Maxwell の関係式の符号を取り違え、( ∂ S / ∂ V ) T = − ( ∂ p / ∂ T ) V (\partial S/\partial V)_T=-(\partial p/\partial T)_V ( ∂ S / ∂ V ) T = − ( ∂ p / ∂ T ) V とすると、理想気体でも ( ∂ U / ∂ V ) T = − 2 p (\partial U/\partial V)_T=-2p ( ∂ U / ∂ V ) T = − 2 p となり Joule の法則に反する。− ( ∂ V / ∂ T ) p -(\partial V/\partial T)_p − ( ∂ V / ∂ T ) p が現れる関係式 ( ∂ S / ∂ p ) T = − ( ∂ V / ∂ T ) p (\partial S/\partial p)_T=-(\partial V/\partial T)_p ( ∂ S / ∂ p ) T = − ( ∂ V / ∂ T ) p はギブズエネルギー G G G から得られるもので、体積で微分する本問では使えない。独立変数が ( T , V ) (T,V) ( T , V ) の対応は A A A 、( T , p ) (T,p) ( T , p ) の対応は G G G と整理しておくと取り違えない。また、第1項だけを残して ( ∂ U / ∂ V ) T = − p (\partial U/\partial V)_T=-p ( ∂ U / ∂ V ) T = − p とするのは、熱の出入りがない断熱過程での関係 d U = − p d V \mathrm{d}U=-p\,\mathrm{d}V d U = − p d V と混同した誤りである。
設問(b) van der Waals気体での式(3)
偏微分で固定する量
( ∂ p / ∂ T ) V (\partial p/\partial T)_V ( ∂ p / ∂ T ) V では V V V と n n n を固定し、式(1)の T T T を含む項だけを微分する。a a a と b b b は問題で「p , T , V p,T,V p , T , V に依存しない定数」と指定されているので、− a n 2 / V 2 -an^2/V^2 − a n 2 / V 2 は T T T について定数として消える。この指定がなければ、たとえば a a a が温度の関数のときは − T d a / d T ⋅ n 2 / V 2 -T\,\mathrm{d}a/\mathrm{d}T\cdot n^2/V^2 − T d a / d T ⋅ n 2 / V 2 のような項が残る。
次元と極限
次元を見ると、a n 2 / V 2 an^2/V^2 a n 2 / V 2 は圧力の次元であり、a a a の単位は Pa\,m6 ^6 6 \,mol− 2 ^{-2} − 2 である。( ∂ U / ∂ V ) T (\partial U/\partial V)_T ( ∂ U / ∂ V ) T の単位 J\,m− 3 ^{-3} − 3 は Pa に等しいので整合する。a → 0 a\to0 a → 0 とすると ( ∂ U / ∂ V ) T → 0 (\partial U/\partial V)_T\to0 ( ∂ U / ∂ V ) T → 0 となり、排除体積 b b b だけをもつ気体では内部エネルギーが体積によらないことが分かる。つまり体積依存性は引力の項 a a a だけが担っている。b → 0 b\to0 b → 0 でも結果は変わらない。
数値で確かめると、n = 1.0 n=1.0 n = 1.0 mol、a = 0.14 a=0.14 a = 0.14 Pa\,m6 ^6 6 \,mol− 2 ^{-2} − 2 、b = 3.2 × 10 − 5 b=3.2\times10^{-5} b = 3.2 × 1 0 − 5 m3 ^3 3 \,mol− 1 ^{-1} − 1 、T = 300 T=300 T = 300 K、V = 1.0 × 10 − 4 V=1.0\times10^{-4} V = 1.0 × 1 0 − 4 m3 ^3 3 で、p = 2.27 × 10 7 p=2.27\times10^{7} p = 2.27 × 1 0 7 Pa、T ( ∂ p / ∂ T ) V = n R T / ( V − n b ) = 3.67 × 10 7 T(\partial p/\partial T)_V=nRT/(V-nb)=3.67\times10^{7} T ( ∂ p / ∂ T ) V = n R T / ( V − nb ) = 3.67 × 1 0 7 Pa なので、− p + T ( ∂ p / ∂ T ) V = 1.40 × 10 7 -p+T(\partial p/\partial T)_V=1.40\times10^{7} − p + T ( ∂ p / ∂ T ) V = 1.40 × 1 0 7 Pa となる。これは a n 2 / V 2 = 0.14 / ( 1.0 × 10 − 4 ) 2 = 1.40 × 10 7 an^2/V^2=0.14/(1.0\times10^{-4})^2=1.40\times10^{7} a n 2 / V 2 = 0.14/ ( 1.0 × 1 0 − 4 ) 2 = 1.40 × 1 0 7 Pa に一致する。
T ( ∂ p / ∂ T ) V T(\partial p/\partial T)_V T ( ∂ p / ∂ T ) V を取り違えたとき
T ( ∂ p / ∂ T ) V T(\partial p/\partial T)_V T ( ∂ p / ∂ T ) V を n R T / V nRT/V n R T / V と書き、理想気体の T ( ∂ p / ∂ T ) V = p T(\partial p/\partial T)_V=p T ( ∂ p / ∂ T ) V = p をそのまま使うと
( ∂ U ∂ V ) T = − n R T V − n b + a n 2 V 2 + n R T V = a n 2 V 2 − n 2 b R T V ( V − n b )
\left(\frac{\partial U}{\partial V}\right)_T
=-\frac{nRT}{V-nb}+\frac{an^2}{V^2}+\frac{nRT}{V}
=\frac{an^2}{V^2}-\frac{n^2bRT}{V(V-nb)}
( ∂ V ∂ U ) T = − V − nb n R T + V 2 a n 2 + V n R T = V 2 a n 2 − V ( V − nb ) n 2 b R T
となり、右辺が T T T に依存する誤った式が出る。n b nb nb を引くべき分母を忘れた結果である。理想気体の p p p との対応は n R T / ( V − n b ) nRT/(V-nb) n R T / ( V − nb ) のまま保つ。
設問(c) 体積増加で内部エネルギーが増える理由
式(3)の物理的な意味
式(3)を V V V で積分すると、C V C_V C V が一定の仮定のもとで
U ( T , V ) = C V T − a n 2 V + U 0
U(T,V)=C_VT-a\frac{n^2}{V}+U_0
U ( T , V ) = C V T − a V n 2 + U 0
の形になる。C V T C_VT C V T は並進などの運動に由来する項、− a n 2 / V -an^2/V − a n 2 / V は引力による位置エネルギーである。分子数密度 n / V n/V n / V を平均場として各分子が感じる引力ポテンシャルは n / V n/V n / V に比例し、それを全分子について足すと n ⋅ ( n / V ) n\cdot(n/V) n ⋅ ( n / V ) に比例する。この平均場の見積りから − a n 2 / V -an^2/V − a n 2 / V という形が出る。V → ∞ V\to\infty V → ∞ では位置エネルギーの項が 0 0 0 に近づき、理想気体の内部エネルギーに戻る。
「引力」を使った説明の組み立て
問題が求めているのは、a > 0 a>0 a > 0 という引力の強さを表す定数が ( ∂ U / ∂ V ) T > 0 (\partial U/\partial V)_T>0 ( ∂ U / ∂ V ) T > 0 を与えることの説明である。「体積が増える」→「平均距離が増える」→「引力に逆らう仕事で位置エネルギーが増える」→「温度一定なら運動エネルギーは同じ」という順に書くと筋が通る。a n 2 / V 2 an^2/V^2 a n 2 / V 2 は内部圧力(引力による凝集圧)と呼ばれ、気体を膨張させる際の分子間の引力の強さを圧力の単位で測ったものである。
斥力や圧力で説明したとき
「体積が増えると圧力が下がるので内部エネルギーが増える」という説明は、圧力は状態量 p ( T , V ) p(T,V) p ( T , V ) であって内部エネルギーの原因ではないので誤りである。また、斥力の寄与を使って「分子が近づくと斥力で位置エネルギーが上がる」と書くと、( ∂ U / ∂ V ) T (\partial U/\partial V)_T ( ∂ U / ∂ V ) T の符号が逆になる。式(3)に b b b が現れないことが示すように、体積の増加で内部エネルギーが増えることに関わるのは引力であり、排除体積ではない。
設問(d) 等温変化での Δ U \Delta U Δ U と Δ S \Delta S Δ S
積分できる理由
U ( T , V ) U(T,V) U ( T , V ) について d U = C V d T + ( ∂ U / ∂ V ) T d V \mathrm{d}U=C_V\,\mathrm{d}T+(\partial U/\partial V)_T\,\mathrm{d}V d U = C V d T + ( ∂ U / ∂ V ) T d V が成り立つ。等温では第1項が消え、残る被積分関数 ( ∂ U / ∂ V ) T = a n 2 / V 2 (\partial U/\partial V)_T=an^2/V^2 ( ∂ U / ∂ V ) T = a n 2 / V 2 が T T T を含まないので、積分は V V V だけの積分になる。エントロピーも同様に d S = ( C V / T ) d T + ( ∂ p / ∂ T ) V d V \mathrm{d}S=(C_V/T)\,\mathrm{d}T+(\partial p/\partial T)_V\,\mathrm{d}V d S = ( C V / T ) d T + ( ∂ p / ∂ T ) V d V で等温ならば第2項だけが残り、( ∂ p / ∂ T ) V = n R / ( V − n b ) (\partial p/\partial T)_V=nR/(V-nb) ( ∂ p / ∂ T ) V = n R / ( V − nb ) は T T T を含まない。このため結果に温度が現れない。問題の「必要なものを用いて」という指定は、各式に実際に現れる量を答えることを求めている。
数値例
仮に n = 1.0 n=1.0 n = 1.0 mol、a = 0.14 a=0.14 a = 0.14 Pa\,m6 ^6 6 \,mol− 2 ^{-2} − 2 、b = 3.2 × 10 − 5 b=3.2\times10^{-5} b = 3.2 × 1 0 − 5 m3 ^3 3 \,mol− 1 ^{-1} − 1 (アルゴンに近い値)、V i = 1.0 × 10 − 4 V_{\mathrm i}=1.0\times10^{-4} V i = 1.0 × 1 0 − 4 m3 ^3 3 、V f = 2.0 × 10 − 4 V_{\mathrm f}=2.0\times10^{-4} V f = 2.0 × 1 0 − 4 m3 ^3 3 とすると、
Δ U = 0.14 × ( 10 4 − 5 × 10 3 ) = 7.0 × 10 2 J , Δ S = 8.314 ln 168 68 = 7.52 J K − 1
\Delta U=0.14\times\left(10^{4}-5\times10^{3}\right)=7.0\times10^{2}\ \mathrm{J},\qquad
\Delta S=8.314\ln\frac{168}{68}=7.52\ \mathrm{J\,K^{-1}}
Δ U = 0.14 × ( 1 0 4 − 5 × 1 0 3 ) = 7.0 × 1 0 2 J , Δ S = 8.314 ln 68 168 = 7.52 J K − 1
である(R = 8.314 R=8.314 R = 8.314 J\,K− 1 ^{-1} − 1 \,mol− 1 ^{-1} − 1 )。極限では、a → 0 a\to0 a → 0 で Δ U → 0 \Delta U\to0 Δ U → 0 、b → 0 b\to0 b → 0 で Δ S → n R ln ( V f / V i ) \Delta S\to nR\ln(V_{\mathrm f}/V_{\mathrm i}) Δ S → n R ln ( V f / V i ) となり、後者は理想気体の等温膨張のエントロピー変化に一致する。V f = V i V_{\mathrm f}=V_{\mathrm i} V f = V i では両方が 0 0 0 になる。
第一法則による確かめもできる。可逆等温変化で系が吸収する熱は Q = T Δ S = n R T ln V f − n b V i − n b Q=T\Delta S=nRT\ln\dfrac{V_{\mathrm f}-nb}{V_{\mathrm i}-nb} Q = T Δ S = n R T ln V i − nb V f − nb 、系がされる仕事は W = − ∫ p d V = − n R T ln V f − n b V i − n b + a n 2 ( 1 V i − 1 V f ) W=-\displaystyle\int p\,\mathrm{d}V=-nRT\ln\dfrac{V_{\mathrm f}-nb}{V_{\mathrm i}-nb}+an^2\left(\dfrac1{V_{\mathrm i}}-\dfrac1{V_{\mathrm f}}\right) W = − ∫ p d V = − n R T ln V i − nb V f − nb + a n 2 ( V i 1 − V f 1 ) である。Q + W = a n 2 ( 1 / V i − 1 / V f ) Q+W=an^2(1/V_{\mathrm i}-1/V_{\mathrm f}) Q + W = a n 2 ( 1/ V i − 1/ V f ) となり、Δ U \Delta U Δ U と一致する。数値例では T = 300 T=300 T = 300 K として Q = 2.26 × 10 3 Q=2.26\times10^{3} Q = 2.26 × 1 0 3 J、W = − 1.56 × 10 3 W=-1.56\times10^{3} W = − 1.56 × 1 0 3 J、Q + W = 7.0 × 10 2 Q+W=7.0\times10^{2} Q + W = 7.0 × 1 0 2 J である。
積分の向きと b b b の扱い
積分の上限と下限を逆にして Δ U = a n 2 ( 1 / V f − 1 / V i ) \Delta U=an^2(1/V_{\mathrm f}-1/V_{\mathrm i}) Δ U = a n 2 ( 1/ V f − 1/ V i ) とすると、上の数値例で − 7.0 × 10 2 -7.0\times10^{2} − 7.0 × 1 0 2 J となり、(c)の結論(膨張で内部エネルギーが増える)に反する。また Δ S = n R ln ( V f / V i ) \Delta S=nR\ln(V_{\mathrm f}/V_{\mathrm i}) Δ S = n R ln ( V f / V i ) と書くと 5.76 5.76 5.76 J\,K− 1 ^{-1} − 1 となり、b b b を無視した理想気体の値である。n b nb nb の n n n を落として V − b V-b V − b と書くと、n ≠ 1 n\ne1 n = 1 の気体で次元が合わない(b b b は物質量あたりの体積)。
設問(e) 等温膨張でのエントロピー変化の比較
理由の述べ方
エントロピーの体積依存性は ( ∂ S / ∂ V ) T = n R / ( V − n b ) (\partial S/\partial V)_T=nR/(V-nb) ( ∂ S / ∂ V ) T = n R / ( V − nb ) で、分子が実際に動ける体積 V − n b V-nb V − nb (自由体積)で決まる。排除体積があると、同じ絶対量 V f − V i V_{\mathrm f}-V_{\mathrm i} V f − V i の膨張が、自由体積の相対的な増加としては大きくなる。実際、( V f − n b ) / ( V i − n b ) (V_{\mathrm f}-nb)/(V_{\mathrm i}-nb) ( V f − nb ) / ( V i − nb ) は V f / V i V_{\mathrm f}/V_{\mathrm i} V f / V i より大きい。統計力学では S = k B ln W S=k_{\mathrm B}\ln W S = k B ln W の W W W が自由体積の N N N 乗に比例するので、自由体積の比が効く、と読める。
引力の項は効かない
( ∂ S / ∂ V ) T = ( ∂ p / ∂ T ) V (\partial S/\partial V)_T=(\partial p/\partial T)_V ( ∂ S / ∂ V ) T = ( ∂ p / ∂ T ) V に引力の項 − a n 2 / V 2 -an^2/V^2 − a n 2 / V 2 は現れない(温度に依存しない)ので、Δ S \Delta S Δ S は a a a に依存しない(∂ Δ S / ∂ a = 0 \partial\Delta S/\partial a=0 ∂ Δ S / ∂ a = 0 )。「分子間引力で動きが制限され、エントロピーの増加が小さくなる」という説明は誤りである。引力は内部エネルギーの体積依存性((c))に入るが、等温変化のエントロピーには入らない。
数値と極限
(d)と同じ数値例では Δ S v d W = 7.52 \Delta S_{\mathrm{vdW}}=7.52 Δ S vdW = 7.52 J\,K− 1 ^{-1} − 1 、Δ S i d = 8.314 ln 2 = 5.76 \Delta S_{\mathrm{id}}=8.314\ln2=5.76 Δ S id = 8.314 ln 2 = 5.76 J\,K− 1 ^{-1} − 1 で、前者のほうが約 1.76 1.76 1.76 J\,K− 1 ^{-1} − 1 大きい。n b → 0 nb\to0 nb → 0 では差が 0 0 0 に近づき、V i → n b V_{\mathrm i}\to nb V i → nb では ln [ ( V f − n b ) / ( V i − n b ) ] → ∞ \ln[(V_{\mathrm f}-nb)/(V_{\mathrm i}-nb)]\to\infty ln [( V f − nb ) / ( V i − nb )] → ∞ となって差が大きくなる。分子がほとんど動けないほど詰まった状態から膨張させると、エントロピー増加が大きいことを表している。
設問(f) 断熱可逆変化で ( V − n b ) T c = (V-nb)T^{c}= ( V − nb ) T c = 一定
使う式と条件
断熱可逆変化は、系が熱を出入りさせず、かつ平衡を保って進む変化である。ここで δ Q = T d S = 0 \delta Q=T\,\mathrm{d}S=0 δ Q = T d S = 0 から、エントロピーは一定に保たれる。この変化では熱の項が消えるので、d U = δ W = − p d V \mathrm{d}U=\delta W=-p\,\mathrm{d}V d U = δ W = − p d V と書ける。内部エネルギーの全微分 d U = C V d T + ( ∂ U / ∂ V ) T d V \mathrm{d}U=C_V\,\mathrm{d}T+(\partial U/\partial V)_T\,\mathrm{d}V d U = C V d T + ( ∂ U / ∂ V ) T d V は、U U U が ( T , V ) (T,V) ( T , V ) の状態関数であることだけで成り立つ。
本問の C V C_V C V は一定という仮定がある。この仮定は van der Waals 気体で矛盾しない。式(2)から ( ∂ C V / ∂ V ) T = T ( ∂ 2 p / ∂ T 2 ) V (\partial C_V/\partial V)_T=T(\partial^2p/\partial T^2)_V ( ∂ C V / ∂ V ) T = T ( ∂ 2 p / ∂ T 2 ) V が導かれ、式(1)の p p p は T T T の1次式なので、右辺は 0 0 0 になる。
エントロピーによる別解
d S = ( C V / T ) d T + ( ∂ p / ∂ T ) V d V = ( C V / T ) d T + n R d V / ( V − n b ) \mathrm{d}S=(C_V/T)\,\mathrm{d}T+(\partial p/\partial T)_V\,\mathrm{d}V=(C_V/T)\,\mathrm{d}T+nR\,\mathrm{d}V/(V-nb) d S = ( C V / T ) d T + ( ∂ p / ∂ T ) V d V = ( C V / T ) d T + n R d V / ( V − nb ) で d S = 0 \mathrm{d}S=0 d S = 0 とおいても、同じ微分方程式が得られる。引力の項 a a a はどちらの方法でも消えるので、断熱可逆変化の不変量は a a a によらない。
ポアソンの式との対応と、d U \mathrm{d}U d U の取り違え
b → 0 b\to0 b → 0 とすると V T c = VT^{c}= V T c = 一定となり、c = 1 / ( γ − 1 ) c=1/(\gamma-1) c = 1/ ( γ − 1 ) より T V γ − 1 = TV^{\gamma-1}= T V γ − 1 = 一定という理想気体のポアソンの式に戻る。単原子分子では c = 3 / 2 c=3/2 c = 3/2 である。
誤りとして多いのは、C V d T = − p d V C_V\,\mathrm{d}T=-p\,\mathrm{d}V C V d T = − p d V と書き、d U = C V d T \mathrm{d}U=C_V\,\mathrm{d}T d U = C V d T とみなしてしまうことである。この式では p p p の中の − a n 2 / V 2 -an^2/V^2 − a n 2 / V 2 が残り、T T T と V V V が分離できず、a a a を含む複雑な式になってしまう。d U \mathrm{d}U d U に ( ∂ U / ∂ V ) T d V (\partial U/\partial V)_T\,\mathrm{d}V ( ∂ U / ∂ V ) T d V の項があることを忘れない。
設問(g) 断熱可逆膨張後の温度
使う関係と条件
断熱可逆膨張は等エントロピー変化であり、(f)の不変量は Δ S = 0 \Delta S=0 Δ S = 0 の結果である。T i , V i , V f , C V T_{\mathrm i},V_{\mathrm i},V_{\mathrm f},C_V T i , V i , V f , C V が両気体で同じという指定により、2つの終端温度の違いはすべて因子 [ ( V i − n b ) / ( V f − n b ) ] 1 / c [(V_{\mathrm i}-nb)/(V_{\mathrm f}-nb)]^{1/c} [( V i − nb ) / ( V f − nb ) ] 1/ c と ( V i / V f ) 1 / c (V_{\mathrm i}/V_{\mathrm f})^{1/c} ( V i / V f ) 1/ c の大小に帰着する。V i > n b V_{\mathrm i}>nb V i > nb であれば両辺は正で、1 / c > 0 1/c>0 1/ c > 0 なので不等式の向きは保たれる。
数値例:断熱膨張後の温度
(d)の数値例に T i = 300 T_{\mathrm i}=300 T i = 300 K、単原子分子 c = 3 / 2 c=3/2 c = 3/2 、V f = 2 V i V_{\mathrm f}=2V_{\mathrm i} V f = 2 V i を加えると
T f v d W = 300 ( 68 168 ) 2 / 3 = 164 K , T f i d = 300 ( 1 2 ) 2 / 3 = 189 K
T_{\mathrm f}^{\mathrm{vdW}}=300\left(\frac{68}{168}\right)^{2/3}=164\ \mathrm K,\qquad
T_{\mathrm f}^{\mathrm{id}}=300\left(\frac12\right)^{2/3}=189\ \mathrm K
T f vdW = 300 ( 168 68 ) 2/3 = 164 K , T f id = 300 ( 2 1 ) 2/3 = 189 K
である。温度が約 25 25 25 K 異なる。上の Δ S = 0 \Delta S=0 Δ S = 0 の式に代入すると、温度項は C V ln ( T f / T i ) = 1.5 n R ln ( 0.547 ) = − 0.905 n R C_V\ln(T_{\mathrm f}/T_{\mathrm i})=1.5\,nR\ln(0.547)=-0.905\,nR C V ln ( T f / T i ) = 1.5 n R ln ( 0.547 ) = − 0.905 n R 、体積項は n R ln ( 168 / 68 ) = + 0.905 n R nR\ln(168/68)=+0.905\,nR n R ln ( 168/68 ) = + 0.905 n R となり、確かに打ち消し合う。
T f T_{\mathrm f} T f が a a a によらない点
(f)の式に a a a は現れないので、断熱可逆変化の終端温度は引力の強さによらない。引力に逆らう仕事のため内部エネルギーが増えること((c))と、引力で圧力が下がり外にする仕事が減ることが、T T T の変化に対して打ち消し合うからである。したがって冷え方に差を生むのは排除体積 b b b だけであり、これは(e)で Δ S \Delta S Δ S の差の原因が b b b だけだったことに対応する。
引力が冷却に効くのは、外部に仕事をしない断熱自由膨張(Joule 膨張)である。この場合は Δ U = 0 \Delta U=0 Δ U = 0 で、C V Δ T = − a n 2 ( 1 / V i − 1 / V f ) C_V\Delta T=-an^2(1/V_{\mathrm i}-1/V_{\mathrm f}) C V Δ T = − a n 2 ( 1/ V i − 1/ V f ) から、数値例では Δ T = − 700 / ( 1.5 × 8.314 ) = − 56 \Delta T=-700/(1.5\times8.314)=-56 Δ T = − 700/ ( 1.5 × 8.314 ) = − 56 K となり、冷却の原因は a a a だけで、理想気体では Δ T = 0 \Delta T=0 Δ T = 0 である。同じ V i → V f V_{\mathrm i}\to V_{\mathrm f} V i → V f でも、可逆断熱膨張(− 136 -136 − 136 K、理想気体は − 111 -111 − 111 K)と自由膨張(− 56 -56 − 56 K)で、温度変化の大きさも原因も異なる。
*[3pt]
自作図:始状態から V V V を増やしたときの断熱可逆変化の温度。n b / V i = 0.32 nb/V_{\mathrm i}=0.32 nb / V i = 0.32 、c = 3 / 2 c=3/2 c = 3/2 の場合
n b nb nb を落とした式と a a a を含めた式
( V − n b ) T c (V-nb)T^{c} ( V − nb ) T c の n b nb nb を落とし、van der Waals 気体にも理想気体の式 T f = T i ( V i / V f ) 1 / c T_{\mathrm f}=T_{\mathrm i}(V_{\mathrm i}/V_{\mathrm f})^{1/c} T f = T i ( V i / V f ) 1/ c を使うと、両者が同じ 189 189 189 K となり差が出ない。反対に、引力が冷却を大きくすると考えて a a a を含む式を立てても、(f)で a a a が消えているので正しくない。
解答
設問(a) 式(2)の導出
閉じた1成分系で、仕事が体積変化の仕事だけの可逆変化を考える。
d U = T d S − p d V
\mathrm{d}U=T\,\mathrm{d}S-p\,\mathrm{d}V
d U = T d S − p d V
を T T T 一定で V V V について偏微分すると
( ∂ U ∂ V ) T = T ( ∂ S ∂ V ) T − p .
\left(\frac{\partial U}{\partial V}\right)_T=T\left(\frac{\partial S}{\partial V}\right)_T-p.
( ∂ V ∂ U ) T = T ( ∂ V ∂ S ) T − p .
ヘルムホルツエネルギー A = U − T S A=U-TS A = U − T S の全微分は
d A = − S d T − p d V
\mathrm{d}A=-S\,\mathrm{d}T-p\,\mathrm{d}V
d A = − S d T − p d V
であり、A A A は状態関数なので d A \mathrm{d}A d A は完全微分である。したがって
( ∂ A ∂ T ) V = − S , ( ∂ A ∂ V ) T = − p
\left(\frac{\partial A}{\partial T}\right)_V=-S,\qquad
\left(\frac{\partial A}{\partial V}\right)_T=-p
( ∂ T ∂ A ) V = − S , ( ∂ V ∂ A ) T = − p
であり、混合2階偏導関数が等しいことから
− ( ∂ S ∂ V ) T = − ( ∂ p ∂ T ) V ⟹ ( ∂ S ∂ V ) T = ( ∂ p ∂ T ) V .
-\left(\frac{\partial S}{\partial V}\right)_T=-\left(\frac{\partial p}{\partial T}\right)_V
\quad\Longrightarrow\quad
\left(\frac{\partial S}{\partial V}\right)_T=\left(\frac{\partial p}{\partial T}\right)_V.
− ( ∂ V ∂ S ) T = − ( ∂ T ∂ p ) V ⟹ ( ∂ V ∂ S ) T = ( ∂ T ∂ p ) V .
これを最初の式に代入して
( ∂ U ∂ V ) T = − p + T ( ∂ p ∂ T ) V
\left(\frac{\partial U}{\partial V}\right)_T=-p+T\left(\frac{\partial p}{\partial T}\right)_V
( ∂ V ∂ U ) T = − p + T ( ∂ T ∂ p ) V
を得る。
設問(b) van der Waals気体での式(3)
式(1)を V V V と n n n を一定にして T T T で偏微分する。a , b a,b a , b は T T T によらないので、第2項は消える。
( ∂ p ∂ T ) V = n R V − n b .
\left(\frac{\partial p}{\partial T}\right)_V=\frac{nR}{V-nb}.
( ∂ T ∂ p ) V = V − nb n R .
したがって
T ( ∂ p ∂ T ) V = n R T V − n b = p + a n 2 V 2 .
T\left(\frac{\partial p}{\partial T}\right)_V=\frac{nRT}{V-nb}=p+a\frac{n^2}{V^2}.
T ( ∂ T ∂ p ) V = V − nb n R T = p + a V 2 n 2 .
最後の等号は、式(1)を n R T / ( V − n b ) = p + a n 2 / V 2 nRT/(V-nb)=p+an^2/V^2 n R T / ( V − nb ) = p + a n 2 / V 2 と書き換えたものである。式(2)に代入すると
( ∂ U ∂ V ) T = − p + ( p + a n 2 V 2 ) = a n 2 V 2 .
\left(\frac{\partial U}{\partial V}\right)_T=-p+\left(p+a\frac{n^2}{V^2}\right)=a\frac{n^2}{V^2}.
( ∂ V ∂ U ) T = − p + ( p + a V 2 n 2 ) = a V 2 n 2 .
設問(c) 体積増加で内部エネルギーが増える理由
分子間に引力が働くと、分子同士が近いほど相互作用の位置エネルギーは低い(負で大きい)。等温で体積を増やすと分子間の平均距離が大きくなり、分子を引力に逆らって引き離すことになるので、分子間の位置エネルギーは増える。温度一定なら分子の運動エネルギーの総和は変わらない。内部エネルギーは運動エネルギーと分子間の位置エネルギーの和なので、体積の増加によって増える。
設問(d) 等温変化での Δ U \Delta U Δ U と Δ S \Delta S Δ S
U U U と S S S は状態量なので、T T T 一定のまま V i V_{\mathrm i} V i から V f V_{\mathrm f} V f まで積分した値が変化量になる。式(3)から
Δ U = ∫ V i V f a n 2 V 2 d V = a n 2 ( 1 V i − 1 V f ) .
\Delta U=\int_{V_{\mathrm i}}^{V_{\mathrm f}}a\frac{n^2}{V^2}\,\mathrm{d}V
=an^2\left(\frac1{V_{\mathrm i}}-\frac1{V_{\mathrm f}}\right).
Δ U = ∫ V i V f a V 2 n 2 d V = a n 2 ( V i 1 − V f 1 ) .
(a)の関係 ( ∂ S / ∂ V ) T = ( ∂ p / ∂ T ) V = n R / ( V − n b ) (\partial S/\partial V)_T=(\partial p/\partial T)_V=nR/(V-nb) ( ∂ S / ∂ V ) T = ( ∂ p / ∂ T ) V = n R / ( V − nb ) を使うと
Δ S = ∫ V i V f n R V − n b d V = n R ln V f − n b V i − n b .
\Delta S=\int_{V_{\mathrm i}}^{V_{\mathrm f}}\frac{nR}{V-nb}\,\mathrm{d}V
=nR\ln\frac{V_{\mathrm f}-nb}{V_{\mathrm i}-nb}.
Δ S = ∫ V i V f V − nb n R d V = n R ln V i − nb V f − nb .
設問(e) 等温膨張でのエントロピー変化の比較
van der Waals 気体と理想気体の Δ S \Delta S Δ S の差は、(d)の結果から
Δ S v d W − Δ S i d = n R ln V f − n b V i − n b − n R ln V f V i = n R ln ( V f − n b ) V i ( V i − n b ) V f .
\Delta S_{\mathrm{vdW}}-\Delta S_{\mathrm{id}}
=nR\ln\frac{V_{\mathrm f}-nb}{V_{\mathrm i}-nb}-nR\ln\frac{V_{\mathrm f}}{V_{\mathrm i}}
=nR\ln\frac{(V_{\mathrm f}-nb)V_{\mathrm i}}{(V_{\mathrm i}-nb)V_{\mathrm f}}.
Δ S vdW − Δ S id = n R ln V i − nb V f − nb − n R ln V i V f = n R ln ( V i − nb ) V f ( V f − nb ) V i .
対数の中の分子と分母の差は
( V f − n b ) V i − ( V i − n b ) V f = n b ( V f − V i ) > 0
(V_{\mathrm f}-nb)V_{\mathrm i}-(V_{\mathrm i}-nb)V_{\mathrm f}=nb\,(V_{\mathrm f}-V_{\mathrm i})>0
( V f − nb ) V i − ( V i − nb ) V f = nb ( V f − V i ) > 0
である(b > 0 b>0 b > 0 、V f > V i V_{\mathrm f}>V_{\mathrm i} V f > V i )。分母は正なので対数の中は1より大きく、差は正になる。
設問(f) 断熱可逆変化で ( V − n b ) T c = (V-nb)T^{c}= ( V − nb ) T c = 一定
断熱可逆変化では δ Q = 0 \delta Q=0 δ Q = 0 なので d U = − p d V \mathrm{d}U=-p\,\mathrm{d}V d U = − p d V である。一方、van der Waals 気体では(b)の結果と C V C_V C V の仮定から
d U = C V d T + ( ∂ U ∂ V ) T d V = C V d T + a n 2 V 2 d V .
\mathrm{d}U=C_V\,\mathrm{d}T+\left(\frac{\partial U}{\partial V}\right)_T\mathrm{d}V
=C_V\,\mathrm{d}T+a\frac{n^2}{V^2}\,\mathrm{d}V.
d U = C V d T + ( ∂ V ∂ U ) T d V = C V d T + a V 2 n 2 d V .
両者を等しいとおき、式(1)の p p p を代入すると
C V d T + a n 2 V 2 d V = − ( n R T V − n b − a n 2 V 2 ) d V .
C_V\,\mathrm{d}T+a\frac{n^2}{V^2}\,\mathrm{d}V=-\left(\frac{nRT}{V-nb}-a\frac{n^2}{V^2}\right)\mathrm{d}V.
C V d T + a V 2 n 2 d V = − ( V − nb n R T − a V 2 n 2 ) d V .
a n 2 / V 2 an^2/V^2 a n 2 / V 2 の項が両辺で打ち消され、
C V d T = − n R T V − n b d V ⟹ c d T T = − d V V − n b ( c = C V / n R ) .
C_V\,\mathrm{d}T=-\frac{nRT}{V-nb}\,\mathrm{d}V
\quad\Longrightarrow\quad
c\,\frac{\mathrm{d}T}{T}=-\frac{\mathrm{d}V}{V-nb}\qquad(c=C_V/nR).
C V d T = − V − nb n R T d V ⟹ c T d T = − V − nb d V ( c = C V / n R ) .
積分すると c ln T + ln ( V − n b ) = c\ln T+\ln(V-nb)= c ln T + ln ( V − nb ) = 一定、すなわち
( V − n b ) T c = 一定 .
(V-nb)\,T^{c}=\text{一定}.
( V − nb ) T c = 一定 .
設問(g) 断熱可逆膨張後の温度
(f)の不変量を始状態 ( T i , V i ) (T_{\mathrm i},V_{\mathrm i}) ( T i , V i ) と終状態 ( T f , V f ) (T_{\mathrm f},V_{\mathrm f}) ( T f , V f ) に使う。
( V i − n b ) T i c = ( V f − n b ) T f c ⟹ T f v d W = T i ( V i − n b V f − n b ) 1 / c .
(V_{\mathrm i}-nb)T_{\mathrm i}^{c}=(V_{\mathrm f}-nb)T_{\mathrm f}^{c}
\quad\Longrightarrow\quad
T_{\mathrm f}^{\mathrm{vdW}}=T_{\mathrm i}\left(\frac{V_{\mathrm i}-nb}{V_{\mathrm f}-nb}\right)^{1/c}.
( V i − nb ) T i c = ( V f − nb ) T f c ⟹ T f vdW = T i ( V f − nb V i − nb ) 1/ c .
理想気体では b = 0 b=0 b = 0 なので T f i d = T i ( V i / V f ) 1 / c T_{\mathrm f}^{\mathrm{id}}=T_{\mathrm i}(V_{\mathrm i}/V_{\mathrm f})^{1/c} T f id = T i ( V i / V f ) 1/ c である。(e)の結果 ( V f − n b ) / ( V i − n b ) > V f / V i (V_{\mathrm f}-nb)/(V_{\mathrm i}-nb)>V_{\mathrm f}/V_{\mathrm i} ( V f − nb ) / ( V i − nb ) > V f / V i から
V i − n b V f − n b < V i V f ⟹ T f v d W < T f i d .
\frac{V_{\mathrm i}-nb}{V_{\mathrm f}-nb}<\frac{V_{\mathrm i}}{V_{\mathrm f}}
\quad\Longrightarrow\quad
T_{\mathrm f}^{\mathrm{vdW}}<T_{\mathrm f}^{\mathrm{id}}.
V f − nb V i − nb < V f V i ⟹ T f vdW < T f id .
理由はエントロピー保存で説明できる。断熱可逆変化では
Δ S = C V ln T f T i + n R ln V f − n b V i − n b = 0
\Delta S=C_V\ln\frac{T_{\mathrm f}}{T_{\mathrm i}}+nR\ln\frac{V_{\mathrm f}-nb}{V_{\mathrm i}-nb}=0
Δ S = C V ln T i T f + n R ln V i − nb V f − nb = 0
である。(e)より、同じ V i → V f V_{\mathrm i}\to V_{\mathrm f} V i → V f の膨張でも van der Waals 気体の体積項のエントロピー増加が大きいので、これを打ち消すため温度項の減少が大きくなる必要があり、T f T_{\mathrm f} T f が低くなる。
最終答
設問(a) 式(2)の導出
式(2)が成り立つ。途中で ( ∂ S / ∂ V ) T = ( ∂ p / ∂ T ) V (\partial S/\partial V)_T=(\partial p/\partial T)_V ( ∂ S / ∂ V ) T = ( ∂ p / ∂ T ) V を使った。
設問(b) van der Waals気体での式(3)
( ∂ U ∂ V ) T = a n 2 V 2 \displaystyle\left(\frac{\partial U}{\partial V}\right)_T=a\frac{n^2}{V^2} ( ∂ V ∂ U ) T = a V 2 n 2 が成り立つ。
設問(c) 体積増加で内部エネルギーが増える理由
体積が増えると、引力に逆らって分子が引き離され位置エネルギーが上がる。温度が同じなら運動エネルギーは変わらないので、内部エネルギーが増える。
設問(d) 等温変化での Δ U \Delta U Δ U と Δ S \Delta S Δ S
Δ U = a n 2 ( 1 V i − 1 V f ) \displaystyle\Delta U=an^2\left(\frac1{V_{\mathrm i}}-\frac1{V_{\mathrm f}}\right) Δ U = a n 2 ( V i 1 − V f 1 ) 、
Δ S = n R ln V f − n b V i − n b \displaystyle\Delta S=nR\ln\frac{V_{\mathrm f}-nb}{V_{\mathrm i}-nb} Δ S = n R ln V i − nb V f − nb 。Δ U \Delta U Δ U に T , R , b T,R,b T , R , b は現れず、Δ S \Delta S Δ S に T , a T,a T , a は現れない。
設問(e) 等温膨張でのエントロピー変化の比較
Δ S v d W > Δ S i d \Delta S_{\mathrm{vdW}}>\Delta S_{\mathrm{id}} Δ S vdW > Δ S id で、van der Waals 気体のほうが大きい。差は n R ln ( V f − n b ) V i ( V i − n b ) V f nR\ln\dfrac{(V_{\mathrm f}-nb)V_{\mathrm i}}{(V_{\mathrm i}-nb)V_{\mathrm f}} n R ln ( V i − nb ) V f ( V f − nb ) V i で、原因は排除体積 b b b である。
設問(f) 断熱可逆変化で ( V − n b ) T c = (V-nb)T^{c}= ( V − nb ) T c = 一定
断熱可逆過程では ( V − n b ) T c = (V-nb)T^{c}= ( V − nb ) T c = 一定(c = C V / n R c=C_V/nR c = C V / n R )が成り立つ。
設問(g) 断熱可逆膨張後の温度
T f v d W < T f i d T_{\mathrm f}^{\mathrm{vdW}}<T_{\mathrm f}^{\mathrm{id}} T f vdW < T f id 。van der Waals 気体のほうが低温になる。