二酸化炭素の比熱比(物性)
今日は,理想気体の断熱過程での気体法則であるポアソン(Poisson)の公式PVγ=一定,または TVγ-1=一定で使用される比熱比 γ= Cp/Cvの値について,考察します。
統計力学によれば,比熱比は対象とする気体1分子を構成する原子の個数,つまり気体分子が単原子分子,2原子分子,3原子分子etcのいずれであるかによって決まります。
ここで, Cv は定積比熱,Cpは定圧比熱です。
(理想)気体に対する定積比熱,と定圧比熱の間にはマイヤー(Mayer)の法則というルールがあり,nモルの気体に対してはCp=Cv+nR (1モルなら Cp=Cv+R )が成り立ちます。
ただし,Rは気体定数と呼ばれる定数で,R≒8.31J/(mol・K)です。
そして,気体の定積比熱 Cvは絶対温度をT,内部エネルギーをUとすると Cv=dU/dTで与えられます。
理想気体ではUは温度だけの関数なので,T=0 での零点エネルギーを無視すると,気体の内部エネルギーはU =CvT と書けます。
古典統計力学によると,物体の常温での内部エネルギーUは,1粒子の運動する自由度1つごとに kBT/2 だけの値を割り当てられます。
ここで kB はボルツマン(Boltzmann)定数と呼ばれる気体分子1個当たりの気体定数です。
kBは気体1分子当たりの気体定数ですから,R=N0kB,またはN0=R/kBとすると気体1モルというのはN0個の分子の集合体を意味することがわかります。
N0はAvogadro数と呼ばれる物理定数で6.02×1023 なる値です。
nモルの気体を構成する分子数はnN0個ですから,それの1自由度あたりの内部エネルギーはnN0kBT/2=nRT/2 です。
以上の事実はエネルギー等分配の法則といわれますす。
単原子分子気体では分子1個の自由度は並進運動の自由度3だけなのでnモルの気体の内部エネルギーはU=3nRT/2 です。そこでCv=3nR/2, Cp=Cv+nR =5nR/2です。
また,2原子分子気体は回転の自由度2 が加わるので,分子1個の自由度は並進運動(重心運動)の自由度3と合わせて5となります。そこでnモルの気体の内部エネルギーはU=5nRT/2となります。Cv=5nR/2, Cp=7nR/2です。
3原子分子以上では重心の周りの回転の自由度が最大の3になるため,これを並進運動(重心運動)の自由度3と合わせると分子1個の自由度は6となります。
それ故,nモルの内部エネルギーはU=3nRT で,Cv=3nR, Cp=4nRです。
そこで,比熱比γ=Cp/Cvは単原子分子気体なら1.67で2原子分子気体なら 1.4,そして3原子分子以上なら特別な対称性がない限り1.33になるはずです。
そこで本当にそうなっているのかどうか?を理科年表で確かめてみると,He 1.66, Ar 1.67, H2 1.40, N2 1.40, H2O 1.31, NH3 1.33 とありました。
これを見ると,必ずしも近似的に理想気体と見なせる希薄気体ではないような実在気体でも,かなり良く適合値を示しているようです。
ここで,ニフティ「物理フォーラム」でのある方からの質問を呈示してみます。
"3原子分子であっても,二酸化炭素 CO2が典型例であるように,一直線に並ぶ3原子分子の場合にはどうなるのだろうか?
もし厳密に一直線なら回転の自由度は2なので2原子分子と同じγ,つまり 1.4になるはずですが,理科年表によると二酸化炭素 CO2のγは1.30でしたから,これは普通の3原子分子に近い値です。"
上記が質問の内容です。
そこで,これに対する答えを見出すために,これまで考えてきた並進や回転の自由度だけではなく,振動の自由度も考慮するとどうなるかを考えてみます。
重心の並進運動や回転の運動とは異なり,振動の自由度なら1方向の調和振動に対しては,位置エネルギーと運動エネルギーの2つの自由度があるので,1方向についての平均エネルギーは 1分子当たりkBTになります。
たとえば静止した固体は3方向に熱振動しているので,常温では1モルにつき,比熱は気体定数をRとして固体の種類によらず3Rとなります。(デュロン・プティ(Dulog-Petit)の法則)
つまり,1次元調和振動子のエネルギーは E=p2/(2m)+(1/2)kx2であり,"Maxwell-Boltzmann分布=古典確率分布"によれば,振動子の座標が(x,p)である確率密度はGibbs因子exp{-E/(kBT)}に比例します。
そこで,エネルギー Eを表わす式の中の1つの変数の2乗を与える変数自由度について,それぞれ平均をとると kBT/2 となりますが, E=p2/(2m)+(1/2)kx2の右辺にはp2と x2 の2つの2乗項があるので振動のエネルギーを考えた場合には,平均エネルギーへの寄与は 1分子当たり一つの方向(1次元)について kBTとなります。
これに対して,重心の自由な並進運動とか,回転運動では位置エネルギーの項はなくて運動エネルギーの項しかない,つまり p2の項しかないので,平均エネルギーへの寄与は1分子当たり1次元について kBT/2 となるのですね。
とにかく,古典統計力学ではax2 exp {- ax2/(kBT)} なる式をx で積分したものを,exp{-ax2/(kBT)}をx で積分したもので割ると,必ずkBT/2 になるということを直接計算で確かめることができます。
これは自由度が1つでもあればそうで,係数aの大きさには無関係です。
ところで常温での固体では,格子を構成する原子のイオンの熱振動がメインになる(電子振動は無視される)のに対して,気体では、原子の重心運動と回転運動のみがメインとなり,熱振動の自由度や電子の運動の自由度が何故無視されるのかという問題があります。
これは量子論ではエネルギーが量子化され,統計分布がPlanck定数hに関係した量子確率分布で与えられるためです。
こうしたことの理由を簡単に言うなら,物質内部のエネルギーを E としその構成粒子の主要な振動数をνとすると,Eは量子論では大体においてhνの倍数で与えられ,量子統計分布では,先のGobbs因子exp{(-E/(kBT)}がexp{-nhν/(kBT)}という形で現われるからです。
常温のTでは固体の電子の振動や気体での原子振動の振動数や電子の自由度に関わる周期運動の振動数νに対しては,一般にhν>>kBTが成立するので,exp{-nhν/(kBT)} ~ 0 となるためこれらは内部エネルギーにはほとんど寄与しないのです。
ところが,問題の二酸化炭素:CO2について「甘泉法師さん」から得た情報によると,"二酸化炭素分子の振動データは,次の振動モードのそれぞれについて,全対称伸縮は実測=1333/cm,計算=1373/cm(12CO2),逆対称伸縮は実測=2349/cm,計算=2420/cm(12CO2),変角振動は 実測=667/cm,計算=669/cm(12CO2)となっているそうです。
一番エネルギーの小さい変角振動について,そのエネルギーを温度に換算すると赤外線温度 1.4387752・667 = 953Kで,常温(300K)の約3倍"なので,振動を無視できないそうです。
実際,変角の振動モードに対して,例えば摂氏(Celsius)16度:T=289Kで x = E/(kBT)=3.32を用いて量子論でのモル比熱を求める式(固体のEinsyeinモデルと同じ式)であるCvib=R x2 exp ( x2 )/[exp ( x2 )-1]2 に代入すると,Cvib=0.43Rとなります。
変角振動は横波なので縦振動を除き自由度が2 であるため,結局Cvib=0.43R ×2=0.86Rであり,比熱比はγ=1+ R/ (5/2R+Cvib)=1.30となって,めでたく理科年表の値と一致します。
ただし,こうして正しい値が得られたのは,振動を除く自由度としては原子が1直線状であることを考慮して2原子分子と同様,定積モル比熱がCv=5/2Rの場合に対応する自由度を想定して計算した結果ですから,やはりCO2では回転の自由度は2である,と考えるのが正解のようです。
http://fphys.nifty.com/(ニフティ「物理フォーラム」サブマネージャー)
http://blog.with2.net/link.php?269343(ブログ・ランキングの投票)
| 固定リンク
「109. 物性物理」カテゴリの記事
- 再掲記事'5):ボーズ・アインシュタイン凝縮 関連(2019.08.01)
- 電気伝導まとめ(2)(2019.07.26)
- 電気伝導まとめ(2)(2019.07.26)
- 電気伝導まとめ(1)(2019.07.25)
- 水蒸気の比熱(2009.02.09)
「104. 熱力学・統計力学」カテゴリの記事
- ボルツマン方程式と緩和時間(2019.08.03)
- -熱交換器・ヒートポンプなど(2019.07.12)
- 統計力学の基礎(6)(量子統計力学4)(2014.12.12)
- 統計力学の基礎(8)(統計力学の基礎付け2)(2014.11.27)
- 統計力学の基礎(6)(量子統計力学3)(2014.11.24)
コメント