热物理学

组成物质的分子或者粒子都在做永不停息的无规则运动,称为热运动。大量分子的热运动的集体效应在宏观上表现为物体的热现象和热性质。

宏观理论——热力学 以观察和实验为基础,通过归纳和推理得出有关热现象的基本规律,因而其结论普遍而且可靠。

微观理论——分子动理论 从分子结构和分子运动出发,应用力学规律和统计方法,研究大量分子热运动的集体效应,从微观本质上解释热现象和热性质。

平衡态

热力学系统—热力学研究的对象,包括极大量的分子、原子

外界:热力学系统外的物体。

宏观量与微观量 宏观量:在整体上描述系统的状态量,如 M,V,E 可以叠加,称为广延量。P,T,不可累加,称为强度量。

微观量 分子的质量,直径,位置,速度,动量,能量等等。

热力学第零定律

热平衡的概念:在与外界影响隔绝的条件下,使两物体(热力学系统)接触,让他们之间能发生传热(这种接触叫做热接触),之后它们的宏观性质不再发生改变,称为热平衡。

当A、B与C同时达到热平衡时,A与B也必然处于热平衡。Thermal equilibrium state 热平衡有 可传递性。

达到热平衡的系统具有共同的内部属性—温度 热平衡具有可传递性:原因具有 共同属性

确定温度数值的表示方法:温标

  1. 选定测温物体

  2. 选定与温度单调变化的属性

  3. 假定测温属性与温度成线性关系 t=aX+b, 其中 t 为温度,X 为水银体积。

  4. 选定温度标准点,将温度计分度

由定压气体温度计和定容气体温度计得到,所建立的温标与温度计的测温物质和测温属性有关。有没有与测温物质无关的统一的温标?

当气体非常稀薄的时候,用实际气体所建立的温标趋于一个共同的极限,称为理想气体温标。

image-20230511085027937

理想气体温标和状态方程

对于一个简单系统,存在以下函数:

F(p,V,T)=0

实验发现:

得到:pV=DT.

当 T0=273.15 K,p0=一个大气压时,得到

D=p0V0T0=ν⏟物质的量p0V0 mol⏞22.4T0≡νR

因此算出 R=8.31 J/(mol⋅K−1).

设气体质量为 m,摩尔质量为 M,气体的分子数为 N,阿伏伽德罗常数 NA,最后得到理想气体的状态方程为:

pV=νRT=mMRT=NNART

引入玻尔兹曼常数:

kB=RNA=1.38×10−23(J⋅K−1)
⇒pV=NNART=kBNT

设气体分子数密度(单位体积内的分子数)为 n,则 p=nkBT.

理想气体微观模型 压强和温度的统计意义

气体分子热运动的特征:

理想气体的微观模型

分子间间距

(VmolNA)13

统计假设:气体处在平衡态时,分子在容器中的空间分布平均来说是均匀的。

n=dNdV=NV

dV 无限大的小,宏观上是小的,微观上是大的。

气体在平衡态时,具有相同速率的分子向各个方向运动的平均分子数是相同的,定义

vx―=∑ivxiNvx2―=∑ivxi2N

由于各向同性,得到 vx―=vy―=vz―=0 和 vx2―=vy2―=vz2―。

所以 vx2―=vy2―=vz2―=13v2―。

设第 i 组分子的速度在 vi―∼vi―+dvi― 区间内,以 ni 表示单位体积内第 i 组分子的数量。m:分子质量,n=∑ni

一次碰撞分子动量变化 2mvix,完全弹性碰撞,在 dt 时间内与 dA 碰撞的分子数 nivixdt⏟能碰到容器壁对应的距离dS⏟面元。

他们给容器壁的总冲量:

2mnivix2dtdS

考虑不同速度分子,一半朝左一半朝右,朝左碰右壁,朝右碰左壁:

dI=∑i,vix>0朝左2mnivix2dtdS=∑imnivix2dtdS

得到

P=dFdS=dI/dtdS=∑imnivix2=nm∑inivxi2∑ini⏟n=nmvx2―=13nmv2―

代入平均分子动能 εt―=12mv2―,得到 P=23nεt―。

代入 P=nkT,结论 εt―=32kT, 温度标志着物体内部分子热运动的剧烈程度,它是大量分子热运动的平均平动动能 εt― 的统计平均值的量度。

均方根速率:εt―=12mv2―=32kT,得到

v2―=3kTm=×NA3RTM

道尔顿分压定律 混合理想气体是由若干种化学纯的气体构成的,设 n1∼nm 为气体的数密度,则总数密度为 n=∑ni. 各组气体的平动动能相等,得到

p=23(∑ini)ε―=∑i23niε―=∑ipi

image-20230511223125559

分子热运动的速度和速率统计分布规律

速率分布函数

设速率分布在 v→v+Δv 间隔内的分子数为 ΔN,则 ΔN 正比于 Δv。比值 ΔN/Δv 代表在速率 v 附近 单位速率区间内分布的分子数。比值 ΔN/NΔv 代表在速率 v 附近 单位速率区间内分布的分子数占总分子数的比例。

f(v)=limΔv→0ΔNNΔv=dNNdv

当速率间隔 Δv 很小的时候,即 Δv=|v1−v2|≪v1,可以用小矩形替代积分,分布在 v→v+Δv 或者 v−Δv/2→v+Δv/2 区间内的分子数可以近似为 dN=f(v)Ndv.

归一化条件:

∫0∞f(v)dv=∫0∞dNN=1N∫0∞dN=1

归一化条件说明 f(v)−v 曲线和 x 轴包围的面积为 1。

平均速率:

v―=∫0∞vf(v)dv

分子速率平方的平均值:

v2―=∫0∞v2f(v)dv

分子动能的平均值为:

εt―=∫0∞12mv2f(v)dv=12m∫0∞v2f(v)dv=12mv2―

f(v) 可以看成权重。广义上来讲:

g―=∫0∞gf(v)dv

速度分布函数

令 dN 表示速度分布在速度范围 (vx→vx+dvx,vy→vy+dvy,vz→vz+dvz) 之间的分子数,简单记为 v→v+d3v,其中 d3v 称为速度体积元。

如果把整个容器内分子的运动速度拿出来,画在三维坐标系里面,就可以发现应该形成一个分布不均的球,这样分布仅仅和速率有关,形成 number density。因为速度分布应该是各向同性的。(类比三维坐标系中的位矢)

速度分布函数可以表示为:

F(vx,vy,vz)=dNNdvxdvydvzF(v)=dNNd3v

其中

dN=NF(vx,vy,vz)dvxdvydvz

归一化条件:

∫−∞∞∫−∞∞∫−∞∞F(v)d3v=1
v2―=∫−∞∞∫−∞∞∫−∞∞v2F(vx,vy,vz)dvxdvydvz
dNvx=Ndvx∫−∞∞∫−∞∞F(vx,vy,vz)dvydvz

定义:

g(vx)=dNvxNdvx=∫−∞∞∫−∞∞F(vx,vy,vz)dvydvz

处于 vx→vx+dvx 的概率为 g(vx)。

vx―=∫−∞∞vxdNvxN=0

下面定量计算处于平衡态的理想气体对容器壁的压强。假设气体分子的数密度为 n,每个分子的质量为 m,分布在 v→v+d3v 的分子具有相同的速度 v。

研究 dA 为底,vxdt 为高的斜柱体中的满足分布在 v→v+d3v 分子数:

dN=[nF(v)dvxdvydvz]⏟v→v+d3v分子数vxdtdA⏟体积

每个分子都会产生 2mvx 的冲量,因此冲量之和:

dI′=2mvxdN

碰撞时要求 vx≥0,因此所有能碰到面元上的分子在 dt 时间给面元 dA 的冲量为积分:

dI=∫0∞∫−∞∞∫−∞∞dI′=∫0∞dvx∫−∞∞dvy∫−∞∞dvz⋅2mvx⋅nF(v)vxdtdA

由 dI=Fdt=pdAdt,得到

p=dIdtdA=2mn∫0∞dvx∫−∞∞dvy∫−∞∞dvzF(v)vx2

速度分布函数应该具有以下特点:

  1. 各向同性,速度分布函数与速度方向无关,也就是:

    F(vx,vy,vz)=F(v2)=F(vx2+vy2+vz2)
  2. 分布对等

    vx2―=vy2―=vz2―=13v2―
  3. 分子对三个速度分量的分布应该是相互独立的,同时处于三个间隔的概率。

    F(vx2+vy2+vz2)=F(v2)=g(vx2)g(vy2)g(vz2)

由于速度分布函数是偶函数,得到

p=mnvx2―=13nmv2―=23nεt―

对比 p=nkBT,得到

εt―=32kBT=12mv2―
v2―=3kBTm=×NA3RTM

道尔顿分压定理:

n=∑inip=∑23niεi―=23∑iniε―=kBT∑ini=∑ipi

混合气体的压强等于各气体的分压之和。

麦克斯韦速度分布律与速率分布律

速度分布律

由 F(vx2+vy2+vz2)=F(v2)=g(vx2)g(vy2)g(vz2) 可得应该取 g(x) 为指数函数,得到

F(vx,vy,vz)=Ae−α(vx2+vy2+vz2)=Ae−αv2

由归一化条件:

∫−∞∞∫−∞∞∫−∞∞F(vx,vy,vz)dvxdvydvz=1

可以得到

A(∫−∞∞e−αvx2dvx)3=1
A=(απ)3/2

接下来确定常数 α,由

vx2―=∫−∞∞vx2g(vx2)dvx∫−∞∞g(vx2)dvx=12α

由温度公式就可以得到

vx2―=13v2―=12α=kTm

得到

α=m2kT

因此得到麦克斯韦速度分布函数

F(vx,vy,vz)=(m2πkT)3/2e−m2kT(vx2+vy2+vz2)

速率分布律

f(v)=4πv2(m2πkBT)3/2e−m2kBTv2

其中 m 是一个分子的质量。

f(v)dv=F(v)4πv2dv

最概然速率

ddvf(v)=0⇒vp=2RTM=1.414kBTm

改造公式:

f(v)=4π⋅π−3/2vp−2v2e−v2/vp2
f(v)=4π−1/2(vvp)2e−v2/vp2/vp

证明速度分布在 vp∼2vp/π 分子比例与气体种类、温度无关。(也就是最概然速率和平均速率之间)

∫vpv―f(v)dv=∫vpv―4π−1/2(vvp)2e−v2/vp2dv/vp=∫vp2vp/π4π−1/2x2e−x2dx=∫12π4π−1/2x2e−x2dx=C

平均速率

∫0∞vf(v)dv=v―=2πvp=8kBTπm=8πkBTm=1.60kBTm

均方根速率

∫0∞v2f(v)dv=v2―=3RTM=1.73kBTm

分子的平均动能(平动)

ε―=12m3RTM=32kBT

结论相同。

image-20230516113344874

没有一个线恰好把曲线平分为面积相同的两半……


image-20230518081715160

注:温度必须一样,否则速率分布不一样。

image-20230518082057258

N(S2−S1)=N(1−A).

分子的速率分布律实验验证 泄流

思考,在容器壁上开一个小孔,从小孔里面泄流出来的气体分子的速率分布能不能代表容器内部的气体分子速率分布。答案是不能,因为气体分子速率越大,从小孔泄出的概率越高,通过下面的推导,我们能够验证泄流速率和速度分布有以下关系:fe(v)∝vf(v).

碰壁频率:单位时间内撞到单位面积上的分子数,以 Γ 表示。

计算得到碰壁频率为:

Γ=n14(8kBTπm)1/2=14nv―

如果从数密度 n 的容器向真空泄流,写成微分方程的形式,其中 S 为小孔的面积:

dNdt=−14nv―S

用麦克斯韦速率分布律可以将碰壁频率表示为:

Γ=14n∫0∞vf(v)dv⏟v―≡∫0∞γ(v)dv

其中 γ(v)=14nvf(v),物理意义为单位时间,从小孔单位面积上射出的在速率 v 附近单位速率间隔内的分子数。

泄流速率分布函数是指分子束中速率在 v 附近单位速率间隔内的分子数占总泄流分子数的比率,也就是代表了泄流气体分子速率的分布函数,由归一化,得到:

fe(v)=γ(v)Γ=14nvf(v)14nv―=vv―f(v)=12(mkBT)2e−m2kBTv2v3

可见:fe(v)∝vf(v).


理想气体,温度 T,求气体分子按平动动能的分布律 f(ε) (ε=12μv2) 并求最可几速率。

需要把 f(v) 变成 f(ε),满足在 v 附近,ε=12mv2 附近:

f(v)dv=f(ε)dε=dNN

代入 v2=2ε/m,dv=dε2με.

f(ε)=4π(μ2πkT)3/2e−μv2/(2kT)v212με=2π(1kT)3/2ε1/2e−ε/(kT)

玻尔兹曼能量分布

问题:在重力场的作用下,如何求 h→h+Δh 范围内的分子数。

玻尔兹曼能量分布:

如果分子的势能相同,由麦克斯韦速度分布函数:

dN=NF(vx,vy,vz)d3v=N(m2πkBT)3/2e−m2kT(vx2+vy2+vz2)d3v∝e−εk/(kBT)

如果气体系统处于外场(如重力场中),每个分子的能量 ε=εk+εp.

可以推出同时处在六个间隔 v→v+d3v,r→r+d3r 的分子几率和指数函数成正比,也就是

dN∝e−εk+εpkT=e−εkT

(体现了能量越小,几率越大)

设归一化因子 C,则处于这个区间内的分子数:

dN=Ce−εkTd3vd3r

写成:

dN=n0(m2πkBT)3/2e−εk+εpkTd3vd3r=n0e−εpkTF(v)d3vd3r

其中 n0 是待定常数。

则处在空间体积元 dxdydz 内的分子数为:

dN′=n0e−εpkTd3r∭F(v)d3v=n0e−εpkTd3r

因此在 x,y,z 附近分子的数密度为:

n=dN′d3r=n0e−εp/kT

其中 n0 代表 εp=0 处单位体积内具有各种速度的分子总数。称为玻尔兹曼分布律。

重力场中 εp=mgz,则

n=n0e−mgz/kT
P=nkT=n0kTe−mgz/kT=P0e−mgz/kT
z=kTmgln⁡P0/P

知道了地面的气体压强和待测高度的气体压强,可以测出高度。


稳定大气温度 T,以地面为重力势能零点,试证明大气分子 Ep―=kT.

取截面积为 dS,高度为 dz 的小体积。

Ep―=∬mgz⋅n0e−mgz/kTdSdz∬n0e−mgz/kTdSdz

分子原子在能级上的分布也遵守玻尔兹曼定律:Ni∝e−EikT

一般假定 Ni=Ce−EikT,然后根据 ∑Ni=N 求出 C.


系统有4000个粒子,能级0,E,2E,初始,三能粒子数分别为2000,1700,300,问是否是平衡态? 平衡态下应如何分布?

  1. 系统分子数 C(1+e−EkT+e−2EkT)=4000

  2. 系统能量守恒 Ce−0kT⋅0+Ce−EkT⋅E+Ce−2EkT⋅2E=2000⋅0+1700⋅E+300⋅2E=2300E.

能量均分定理

自由度:决定某物体在空间的位置所需要的独立坐标数。

能量均分定理(Equipartition Theorem)

在温度为 T 的平衡态下,物质分子的每一个自由度都具有相同的平均动能,其大小都等于 kBT/2,对每一振动自由度还有 kBT/2 的平均势能。

设平动自由度 t,转动自由度 r,振动自由度 s,则 一个分子 的平均能量:

E―=(t+r+2s)12kT=i2kT

理想气体内能

E=Ek(对应t+r+s)+Epa(对应s)+Epm(理想气体无分子间势能)

32kT 是一个分子的平均动能。

对于 ν mol 气体,得到

E=N(i2kT)=νNA(i2kT)=ν(i2RT)=i2PV

因此:

E=i2PV

理想气体的内能只能是温度的函数。


绝热容器被绝热板分隔为 A,B 两部分,两部分的体积、压强相等,均为 P0,V0,A 内储存 1 mol 单原子理想气体,B 内储存 2 mol 双原子理想气体。

求

  1. 内能 EA,EB。

    EA=1⋅32RTA=32P0V0

    EB=2⋅52RTB=52P0V0.

  2. 抽出绝热板,两种气体混合后处于平衡时的温度。

    混合气体能量守恒,也就是:

    EA+EB=1⋅32RT+2⋅52RT=132RT
    T=8P0V013R

容器内贮有质量为 m 摩尔质量为 M 的理想气体, 设容器以速度 u 作定向运动,今使容器突然停止,问:

  1. 定向运动机械能转化为什么形式的能量?转化为无规则热运动的动能。

  2. 分子速度平方的平均值增加多少?(单原子、双原子)

    对于单原子分子,容器动能全部转化为分子的平动动能,因此

    12mu2=Δ(mMNA32kBT)=Δ(mMNAεt―)=Δ(mMNA12MNAv2―)

    u2=Δv2―.

    对于双原子分子,还有一部分能量转化为分子的转动动能。

    12mu2=Δ(mMNA52kBT)=Δ(mMNA53εt―)=Δ(mMNA56MNAv2―)

    因此 Δv2―=35u2

习题

19-9

image-20230611133801013

最概然速率:vp=2kBTm

平均速率:v―=8πkBTm=4πvp.

将 f(v) 用 vp 代换:

f(v)=4π(1πvp2)3/2e−v2/vp2v2=4π−1/2e−v2/vp2v2/vp2⋅1/vp

速率在 vp 到 v― 范围内的分子数占总分子数的百分比可以由以下式子计算:

∫vpv―4π−1/2e−v2/vp2v2/vp2⋅1/vpdv=换元u=v/vp∫14π4π−1/2e−u2u2du=Const

如果是 vp∼v2―? 也是不变的。考虑将速度参量化:v=kTmx,因此只要 x 从常数到常数,占比不变。