跳到主要内容

東京大学 工学系研究科 2011年8月実施 物理学 第3問

Author​

祭音Myyura (co-authored with GPT 6 Astra)

Description​

理想気体の状態方程式は PV=nRTPV=nRT、気体定数は R=kNAR=kN_A である。P,V,n,T,k,NAP,V,n,T,k,N_A はそれぞれ圧力、体積、モル数、温度、ボルツマン定数、アボガドロ定数を表す。圧縮因子 Z=PV/(nRT)Z=PV/(nRT) は、不完全気体について

Z=1+nVB(T)+(nV)2C(T)+(nV)3D(T)+⋯Z=1+\frac nV B(T)+\left(\frac nV\right)^2C(T)+\left(\frac nV\right)^3D(T)+\cdots

とビリアル展開できる。B,C,DB,C,D をそれぞれ第 2、第 3、第 4 ビリアル係数という。

I​

第 2 ビリアル係数と分子間ポテンシャル ϕ(r)\phi(r) の関係は

B(T)=2πNA∫0∞[1−exp⁡(−ϕ(r)kT)]r2 drB(T)=2\pi N_A\int_0^\infty\left[1-\exp\left(-\frac{\phi(r)}{kT}\right)\right]r^2\,dr

である。正の定数 σ,ε0\sigma,\varepsilon_0 に対して

ϕ(r)={+∞,r<σ,−ε0(σ/r)6,r≥σ\phi(r)=\begin{cases} +\infty,&r<\sigma,\\ -\varepsilon_0(\sigma/r)^6,&r\ge\sigma \end{cases}

とするとき、B(T)B(T) を求めよ。気相が維持される条件を利用してよい。

II​

van der Waals の状態方程式

(P+an2V2)(Vn−b)=RT\left(P+\frac{an^2}{V^2}\right)\left(\frac Vn-b\right)=RT

をビリアル展開し、B(T),C(T)B(T),C(T) を求めよ。a,ba,b は気体固有の正の定数である。

III​

I、II の第 2 ビリアル係数を関係付けて、a,ba,b を導出せよ。

IV​

a,b,σ,ε0a,b,\sigma,\varepsilon_0 それぞれの物理的意味を説明せよ。

题目描述​

理想气体满足 PV=nRTPV=nRT 且 R=kNAR=kN_A,其中 P,V,n,TP,V,n,T 分别为压强、体积、摩尔数和温度,kk 为玻尔兹曼常数,NAN_A 为阿伏伽德罗常数。真实气体的压缩因子为

Z=PVnRT=1+nVB(T)+(nV)2C(T)+(nV)3D(T)+⋯ .Z=\frac{PV}{nRT}=1+\frac nVB(T)+\left(\frac nV\right)^2C(T)+\left(\frac nV\right)^3D(T)+\cdots.

B,C,DB,C,D 分别称为第二、第三、第四维里系数。

I. 已知

B(T)=2πNA∫0∞[1−e−ϕ(r)/(kT)]r2 dr,ϕ(r)={+∞,r<σ,−ε0(σ/r)6,r≥σ,B(T)=2\pi N_A\int_0^\infty[1-e^{-\phi(r)/(kT)}]r^2\,dr, \qquad \phi(r)=\begin{cases}+\infty,&r<\sigma,\\-\varepsilon_0(\sigma/r)^6,&r\ge\sigma,\end{cases}

其中 σ,ε0>0\sigma,\varepsilon_0>0。求 B(T)B(T),可利用保持气相的条件进行近似。

II. 对范德瓦耳斯方程

(P+an2V2)(Vn−b)=RT\left(P+\frac{an^2}{V^2}\right)\left(\frac Vn-b\right)=RT

作维里展开,求 B(T),C(T)B(T),C(T),a,ba,b 为气体特有的正常数。

III. 比较 I、II 的第二维里系数,导出 a,ba,b。

IV. 解释 a,b,σ,ε0a,b,\sigma,\varepsilon_0 的物理意义。

Kai​

I​

u=ε0/(kT)u=\varepsilon_0/(kT) と置き、積分を r=σr=\sigma で分ける。引力部分の指数関数を展開すると

B(T)=2πNAσ33−2πNA∑j=1∞ujσ6jj!∫σ∞r2−6j dr=2πNAσ33[1−∑j=1∞ujj!(2j−1)].\begin{aligned} B(T) &=\frac{2\pi N_A\sigma^3}{3} -2\pi N_A\sum_{j=1}^\infty\frac{u^j\sigma^{6j}}{j!} \int_\sigma^\infty r^{2-6j}\,dr\\ &=\frac{2\pi N_A\sigma^3}{3} \left[1-\sum_{j=1}^\infty\frac{u^j}{j!(2j-1)}\right]. \end{aligned}

高温の気相で u≪1u\ll1 とする一次近似では、

B(T)≃2πNAσ33(1−ε0kT).\boxed{B(T)\simeq\frac{2\pi N_A\sigma^3}{3} \left(1-\frac{\varepsilon_0}{kT}\right)}.

II​

c=n/Vc=n/V と置くと、

P=cRT1−bc−ac2,Z=11−bc−aRTc=1+(b−aRT)c+b2c2+⋯ .P=\frac{cRT}{1-bc}-ac^2, \qquad Z=\frac1{1-bc}-\frac{a}{RT}c =1+\left(b-\frac a{RT}\right)c+b^2c^2+\cdots.

したがって

B(T)=b−aRT,C(T)=b2.\boxed{B(T)=b-\frac a{RT}}, \qquad \boxed{C(T)=b^2}.

III​

I の一次近似と係数を比較し、R=kNAR=kN_A を用いると

b=2π3NAσ3,a=2π3NA2ε0σ3.\boxed{b=\frac{2\pi}{3}N_A\sigma^3}, \qquad \boxed{a=\frac{2\pi}{3}N_A^2\varepsilon_0\sigma^3}.

IV​

  • σ\sigma:分子中心同士が近づける最小距離、すなわち剛体球としての分子直径。
  • ε0\varepsilon_0:接触距離 r=σr=\sigma における引力ポテンシャルの深さ。
  • bb:分子の有限な大きさによる 1 モル当たりの排除体積補正。剛体球の実体積の 4 倍に相当する。
  • aa:分子間引力の強さを表す圧力補正係数。引力による圧力低下は a(n/V)2a(n/V)^2 である。