系统、边界与热功符号
本章先研究单组分、简单可压缩的封闭系统。系统边界允许热量和体积功交换,但不允许物质穿过;粒子数 N 固定。扩展到开放系统时再把 N 作为变量。取热量 δQ>0 表示能量以热的形式进入系统,功 δW>0 表示系统对环境做功。第一定律写成
dU=δQ−δW.
对仅有准静态 pV 功的可逆过程,δW=pdV,且 δQrev=TdS,于是基本关系为
dU=TdS−pdV+μdN.
U 的单位为焦耳(J),T 为开尔文(K),S 为 JK−1,p 为帕斯卡(Pa),V 为 m3,化学势 μ 为每粒子 J 或每摩尔 Jmol−1。pdV 的单位为 Pam3=mathrmJ。
这条基本关系是相邻平衡态的状态函数微分,不等于宣称任意实际过程都可逆。不可逆过程的边界功一般是 pextdV,热量也不能直接写成 TdS;第二定律需要另加熵产生。势函数的平衡判据必须连同环境约束一起使用。
为什么要更换自然变量
内能的自然变量是 (S,V,N):
T=(∂S∂U)V,N,p=−(∂V∂U)S,N,μ=(∂N∂U)S,V.
实验中更常控制温度和压强,而不是熵。Legendre 变换用共轭变量替换难控制的自然变量,同时保留完整状态信息。以 T=(∂U/∂S)V,N 替换 S,要求在所考察单相稳定区域内 S↦T 至少局部可逆;通常对应正热容和适当凸性。若映射不单值、响应函数发散或系统位于一阶相变共存区,单一光滑分支会失效,需要用稳定包络和相共存条件处理。Legendre 变换不是把字母机械换掉。
Hamilton 力学也使用 Legendre 变换,但那里常在速度与动量之间变换;热力学这里更换的是状态空间中的控制变量,并由稳定性和相平衡限制允许分支。
四个势函数及其自然变量
内能 (U(S,V,N))
dU=TdS−pdV+μdN.
它适合绝热刚性边界并固定组成的系统。固定 S,V,N 时,稳定平衡态使 U 取局部最小值;更根本的孤立系统表述则是固定 U,V,N 时熵取最大值。
焓 (H(S,p,N))
定义 H=U+pV。全微分为
dH=TdS+Vdp+μdN.
它把体积 V 换成压强 p。在仅有 pV 功且外压恒定的封闭过程中,ΔH 常等于进入系统的热量,但还需要初末态平衡以及动能、势能变化可忽略;不能对任意过程把热量定义成焓变。
Helmholtz 自由能 (F(T,V,N))
定义 F=U−TS:
dF=−SdT−pdV+μdN.
所以
S=−(∂T∂F)V,N,p=−(∂V∂F)T,N.
与温度为 T 的热库接触、体积和组成固定时,自发过程满足 ΔF≤0,平衡态使 F 最小。该结论包含系统加热库的总熵不减;它不是孤立地观察系统熵得出的。定温定容且无其他功输入时,可逆过程给出的最大对外非体积功为 −ΔF。
Gibbs 自由能 (G(T,p,N))
定义 G=U+pV−TS=H−TS:
dG=−SdT+Vdp+μdN.
在恒温、恒压、固定组成并只与相应热压环境交换时,自发变化满足 ΔG≤0。可逆极限下,−ΔG 是最大非膨胀功。若外界温度或压强并不恒定,直接比较系统 G 不能替代完整熵与能量账本。
例 1:理想气体等温膨胀与 Helmholtz 自由能
1.00mol 理想气体与 300K 热库接触,从 V1=10.0L 可逆等温膨胀到 V2=20.0L。取功为系统对外做功。理想气体内能只依赖温度,故 ΔU=0。熵变为
ΔS=nRlnV1V2=8.314ln2=5.76JK−1. 于是
ΔF=ΔU−TΔS=−1.73×103J. 可逆膨胀功
W=∫V1V2VnRTdV=nRTln2=1.73×103J=−ΔF. 热库向系统提供同样大小的热量,使 ΔU=Q−W=0。若膨胀不可逆,系统仍有相同初末态 ΔF,但实际功小于该可逆上限,差额与熵产生有关。
Maxwell 关系:混合偏导的物理用途
若势函数在所研究单相区域具有连续二阶偏导,则混合偏导次序可交换。从 U(S,V,N) 得
(∂V∂T)S,N=−(∂S∂p)V,N.
从 H(S,p,N) 得
(∂p∂T)S,N=(∂S∂V)p,N.
从 F(T,V,N) 得
(∂V∂S)T,N=(∂T∂p)V,N.
从 G(T,p,N) 得
(∂p∂S)T,N=−(∂T∂V)p,N.
右侧常是可由状态方程或膨胀实验测得的量,左侧则涉及较难直接测量的熵变化。符号来自势函数全微分,不能仅凭记忆排列。
例 2:由状态方程求等温熵差
2.00mol 理想气体在固定 T 下从 10.0L 变到 15.0L。由 p=nRT/V,
(∂T∂p)V,N=VnR. 利用 F 导出的 Maxwell 关系,
ΔS=∫V1V2(∂V∂S)T,NdV=nRlnV1V2=6.74JK−1. 结果只依赖初末平衡态。实际路径即使不可逆,系统熵差仍相同;但用 δQ/T 计算时只能沿连接同一端点的可逆路径积分。
把熵和能量微分改写成可测量形式
Maxwell 关系的价值在于把含熵的导数换成状态方程导数。以 (T,V,N) 为独立变量展开熵:
dS=(∂T∂S)V,NdT+(∂V∂S)T,NdV=TCVdT+(∂T∂p)V,NdV.
代入 dU=TdS−pdV,得到
dU=CVdT+[T(∂T∂p)V,N−p]dV.
第二项说明真实物质的内能一般既依赖温度也依赖体积。理想气体有 p=nRT/V,所以方括号恰为零,才得到 U=U(T)。这不是所有气体的定义,而是理想气体状态方程与 Maxwell 关系共同推出的特殊结果。
类似地,以 (T,p,N) 展开焓可得
dH=CpdT+[V−T(∂T∂V)p,N]dp=CpdT+V(1−Tα)dp.
稳态节流装置在绝热、无轴功且进出口宏观动能与势能变化可忽略时近似满足单位质量或单位摩尔焓不变。由 dH=0 得 Joule–Thomson 系数
μJT=(∂p∂T)H=CpV(Tα−1),
单位为 KPa−1。μJT>0 时降压导致降温,μJT<0 时降压导致升温。理想气体 α=1/T,所以 μJT=0。这个结论只适用于节流的等焓边界,不能用于活塞绝热膨胀或自由膨胀而不重新做能量账本。
响应函数与稳定性
定义定容和定压热容、体膨胀系数、等温压缩率:
CV=T(∂T∂S)V,N,Cp=T(∂T∂S)p,N,
α=V1(∂T∂V)p,N,κT=−V1(∂p∂V)T,N.
CV,Cp 的单位为 JK−1,α 为 K−1,κT 为 Pa−1。对简单可压缩稳定系统,利用 G 的 Maxwell 关系可得
Cp−CV=κTTVα2≥0.
这里四个量必须都取总量或都取一致的摩尔、质量归一化形式。正的 CV 与 κT 对应通常的局部热稳定和机械稳定;接近临界点时响应函数可能很大,普通线性展开的范围随之缩小。
例 3:用膨胀与压缩数据求热容差
某均匀样品在 T=300K、V=2.00×10−3m3 附近有
α=3.00×10−4K−1、
κT=5.00×10−10Pa−1。则
Cp−CV=5.00×10−10(300)(2.00×10−3)(3.00×10−4)2=108JK−1. 量纲核验为
Km3K−2/Pa−1=mathrmPam3K−1=mathrmJK−1。若给的是摩尔体积,结果就是摩尔热容差,单位应改为 Jmol−1K−1。
曲率、稳定分支与相共存
一阶导数给出 T,p,S,V 等共轭量,二阶导数控制小扰动的稳定性。固定 V,N 时
(∂T2∂2F)V,N=−TCV≤0,
所以稳定系统的 F 对温度通常是凹的;固定 T,N 时
(∂V2∂2F)T,N=−(∂V∂p)T,N=VκT1≥0,
所以 F 对体积通常是凸的。若候选状态给出 CV<0 或 κT<0,它在普通短程宏观系统的相应约束下不是稳定均匀平衡分支。实际系统可能分裂成两相,使总自由能落在更低的凸包上。
一阶相变共存区内,Gibbs 自由能本身连续,但熵 S=−(∂G/∂T)p 或体积 V=(∂G/∂p)T 可出现跳跃;潜热和体积差正来自这些一阶导数的差。临界点处两相差异消失,而热容、压缩率等二阶响应可能显著增大。因而在相变附近用单一低阶 Taylor 展开时,必须写出所属相和有效邻域。
化学势与相平衡
允许粒子交换时,化学势是加入粒子的边际自由能:
μ=(∂N∂G)T,p.
对均匀单组分广延系统,Euler 齐次性给出 G=μN。两个相 α、β 在平衡时必须共享温度和压强,并满足
μα(T,p)=μβ(T,p).
若化学势不同,物质会向降低总 Gibbs 自由能的方向转移。沿相共存曲线,两相化学势的微分相等。对每摩尔量 s、v,dμ=−sdT+vdp,所以
dTdp=ΔvΔs=TΔvL,
这就是 Clapeyron 方程。L 为每摩尔潜热,单位 Jmol−1;Δv 为 m3mol−1,故斜率单位为 PaK−1。公式适用于两相平衡曲线,并假设所用潜热和体积差对应同一 T,p。
例 4:由潜热估算相界斜率
某物质在 T=350K 的液—气平衡处,摩尔潜热为
L=40.0kJmol−1,气液摩尔体积差为
Δv=2.80×10−2m3mol−1。相界局部斜率为
dTdp=(350)(2.80×10−2)4.00×104=4.08×103PaK−1. 正斜率来自 L>0、Δv>0。若固液转变的体积差为负而熔化潜热为正,斜率会为负。用这一局部斜率外推很大温区会忽略 L 和 Δv 随温压的变化。
Euler 关系、Gibbs–Duhem 约束与混合物
对均匀单相的广延系统,把 S,V,N 同时放大因子 λ,内能满足
U(λS,λV,λN)=λU(S,V,N)。Euler 齐次定理结合基本关系给出
U=TS−pV+μN,G=U+pV−TS=μN.
对 Euler 关系求微分,再与 dU=TdS−pdV+μdN 比较,得到 Gibbs–Duhem 关系
SdT−Vdp+Ndμ=0.
因此单组分体系的 T,p,μ 不是三个独立强度变量。固定温度时,化学势随压强的变化满足 dμ=Vmdp,其中 Vm=V/N;固定压强时,dμ=−SmdT。这些关系也是 Clapeyron 推导的局部基础。
多组分系统把基本关系扩展为
dU=TdS−pdV+i∑μidNi,G=i∑μiNi.
相平衡不仅要求两相温度和压强相等,还要求每个可交换组分分别满足 μiα=μiβ。若发生化学反应,粒子数变化还受化学计量约束;不能把总粒子数当成唯一独立组成变量。这里的单组分 Clapeyron 方程是这一更一般条件在两相、单组分情况下的结果。
总量、摩尔量与参考零点
U,H,F,G,S,V,N 是广延量:把两个互不作用的同类系统合并,它们相加。T,p 和单一纯相的 μ 是强度量,不随系统整体等比例放大。摩尔量用小写或下标 m 表示,例如 g=G/n 的单位为 Jmol−1;比量按质量归一化,例如比焓单位为 Jkg−1。在 Cp−CV=TVα2/κT 或 Clapeyron 方程中,总量、摩尔量和比量必须成套使用,混合归一化会留下错误的 mol 或 kg 量纲。
内能和自由能的绝对参考零点可按约定平移,许多宏观预测依赖差值或导数。平移所有候选态相同的常数不会改变压强、熵或平衡位置;若常数随粒子数或相的组成不同,则会改变化学势,不能随意删除。实际计算应声明参考态,并优先核对 ΔF、ΔG、势的一阶导数和两相比值,而不是把任意参考下的绝对数值解释成可直接测量能量。
常见错误
自由能下降适用于任何封闭系统
F 的下降判据要求恒温、恒容、固定组成及合适热库;G 的下降判据要求恒温、恒压、固定组成及合适环境。约束改变时,应换势函数或直接计算总熵。
Legendre 变换只是把 S 写成 T
变换同时减去共轭乘积并改变自然变量。还需状态关系在目标单相区域允许局部反演;相共存或不稳定分支要由稳定包络处理。
不可逆过程不能计算状态函数差
不可逆路径仍有确定的 ΔU、ΔS、ΔF、ΔG,只要初末态是确定平衡态。不能沿实际不可逆路径把 δQ/T 当成系统熵微分。
练习
练习 1:从定义重建四个全微分
- 所属知识
- Legendre 变换
- 难度
- 3/5
不查表推导 H、F、G 的全微分,并列出 U,H,F,G 的自然变量。
查看提示
从
dU=TdS−pdV+μdN 出发,分别对 pV 与 TS 使用乘积法则。
查看解答
dH=d(U+pV)=TdS+Vdp+μdN;
dF=d(U−TS)=−SdT−pdV+μdN;
dG=d(U+pV−TS)=−SdT+Vdp+μdN。相消项决定自然变量,不能漏掉 V dp 或
−SdT。
练习 2:选择正确的平衡势
- 所属知识
- 自发判据
- 难度
- 3/5
分别为恒温定容、恒温定压、绝热刚性和完全孤立四种边界选择平衡判据,并写明固定变量。
查看提示
先写清系统与热库、压强库的交换方式,再匹配固定变量。
查看解答
恒 T、V、N 选 Helmholtz 自由能 F;恒 T、p、N 选 Gibbs 自由能 G;固定 S、V、N 可用 U 的最小判据;孤立且固定 U、V、N 用 S 最大判据。若温度或压强环境随过程变化,不能直接套单一势的单调性。
练习 3:由 Maxwell 关系求熵的压强导数
- 所属知识
- Maxwell 关系
- 难度
- 3/5
推导 (∂S/∂p)T=−Vα,并解释普通正热膨胀材料在等温增压时熵的变化方向。
查看提示
从
dG=−SdT+Vdp 出发,交换 G 对 T、p 的混合偏导。
查看解答
由
S=−(∂G/∂T)p、
V=(∂G/∂p)T,得到
(∂S/∂p)T=−(∂V/∂T)p=−Vα。若
α>0,等温增压使熵降低。右侧单位为
m3⋅K−1,等于
J⋅K−1⋅Pa−1。
练习 4:检验热容差的单位与符号
- 所属知识
- 响应函数
- 难度
- 3/5
不代数值,证明稳定简单可压缩系统通常满足 Cp≥CV,并核对总量和摩尔量两种单位。
查看提示
使用
Cp−CV=TVα2/κT;稳定系统通常有
κT>0。
查看解答
α2 非负,T、V、
κT 对普通稳定系统为正,因此
Cp−CV≥0。单位为
K⋅m3⋅K−2/Pa−1=Pa⋅m3⋅K−1=J⋅K−1。若使用摩尔体积,结果为
J⋅mol−1⋅K−1。负压缩率提示所取分支不具普通机械稳定性。
练习 5:两相物质转移方向
- 所属知识
- 化学势
- 难度
- 4/5
证明两相在恒温恒压下会从高化学势相向低化学势相转移,并给出平衡条件。
查看提示
在恒 T、p 下写
dG=μαdNα+μβdNβ,并用
dNβ=−dNα。
查看解答
令
dNα>0 表示物质从
β 转入
α,则
dG=(μα−μβ)dNα。自发方向要求 dG<0,所以当
μα<μβ 时物质从
β 流向
α。平衡时
μα=μβ。该结论要求温度、压强和总粒子数约束明确。
练习 6:Clapeyron 符号判断
- 所属知识
- 相平衡
- 难度
- 4/5
对固—液共存线,分别讨论液相摩尔体积大于和小于固相时 dp/dT 的符号。
查看提示
熔化潜热为正,斜率符号由
Δv=vliquid−vsolid 决定。
查看解答
dp/dT=L/(TΔv)。L>0、T>0,因此若液相比固相体积大,
Δv>0,熔化线斜率为正;若液相体积更小,
Δv<0,斜率为负。这个判断只给局部共存线方向,不提供远离该点的完整相图。
关系与后续
已登记课程资源
课程 · 2008Thermodynamics & Kinetics
Keith A. Nelson, Moungi Bawendi
用于核对 P05 的符号约定、循环效率、熵平衡、热力学势、Maxwell 关系和相平衡计算。
打开官方来源
MIT OpenCourseWare 5.60 覆盖自由能、化学势和相平衡,可用于继续核对势函数和相图习题。本章中的符号约定、自然变量、适用条件和数值算例均已在正文中展开,资源卡不替代推导。