§1.1 基本概念§1.2 第零定律与温度§1.3 第一定律§1.4 第二定律与熵§1.5 平衡条件与稳定性条件§1.6 热力学势与勒让德变换§1.7 麦克斯韦关系§1.8 响应函数与偏导数技巧§1.9 广延性:欧拉关系与吉布斯–杜安关系§1.10 第三定律§1.11 理想气体与化学平衡速查§1.12 本章小结自测题 本章目的
只 回顾统计热力学中真正会用到的热力学内容。热力学定律是经验总结(公设),不加证明;由定律推出的关系则逐一推导。
§1.1 基本概念
系统的类型 :孤立系统(与外界既无能量交换也无物质交换);封闭系统(有能量交换、无物质交换);开放系统(两者都有)。
平衡态 :在不变的外界条件下,宏观性质不随时间变化、且内部没有宏观流动(热流、粒子流)的状态。
状态函数 :只由平衡态本身决定、与如何到达该状态无关的量,如 U , V , p , T , S U, V, p, T, S U , V , p , T , S 。
广延量与强度量 :把系统复制 λ \lambda λ 份并合并,广延量(U , V , N , S , F , G U,V,N,S,F,G U , V , N , S , F , G 等)变为原来的 λ \lambda λ 倍;强度量(T , p , μ T,p,\mu T , p , μ )不变。
准静态过程 :进行得足够缓慢、每一时刻都可视为处于平衡态的过程。没有耗散的准静态过程是可逆过程。
物态方程 :平衡态下状态变量之间的函数关系。最重要的例子是理想气体:p V = N k B T = n R T pV = Nk_{\mathrm B}T = nRT p V = N k B T = n RT (n n n 为物质的量)。
§1.2 第零定律与温度
若 A 与 C 处于热平衡,B 与 C 也处于热平衡,则 A 与 B 处于热平衡。这使我们可以定义一个"所有相互处于热平衡的系统都相同"的状态量——温度。第3章 将给出温度的微观定义 (3.5) 。
§1.3 第一定律
d U = δ Q + δ W (1.1) dU = \delta Q + \delta W \tag{1.1} d U = δ Q + δ W ( 1.1 )
U U U 是内能(状态函数);δ Q \delta Q δ Q 是系统吸收的热;δ W \delta W δ W 是外界对系统做的功,对准静态过程 δ W = − p d V \delta W = -p\,dV δ W = − p d V 。记号 δ \delta δ 强调 Q Q Q 、W W W 不是状态函数:δ Q \delta Q δ Q 、δ W \delta W δ W 只是过程中传递的小量,其累积值依赖于路径。
全微分与恰当微分
函数 f ( x , y ) f(x,y) f ( x , y ) 的全微分是 d f = ( ∂ f / ∂ x ) y d x + ( ∂ f / ∂ y ) x d y df = (\partial f/\partial x)_y\,dx + (\partial f/\partial y)_x\,dy df = ( ∂ f / ∂ x ) y d x + ( ∂ f / ∂ y ) x d y 。微分式 M ( x , y ) d x + L ( x , y ) d y M(x,y)\,dx + L(x,y)\,dy M ( x , y ) d x + L ( x , y ) d y 能写成某个函数的全微分(称为恰当微分 )的充要条件(在单连通区域上)是 ∂ M / ∂ y = ∂ L / ∂ x \partial M/\partial y = \partial L/\partial x ∂ M / ∂ y = ∂ L / ∂ x ;必要性来自"混合二阶偏导数与求导次序无关":∂ 2 f / ∂ y ∂ x = ∂ 2 f / ∂ x ∂ y \partial^2 f/\partial y\partial x = \partial^2 f/\partial x\partial y ∂ 2 f / ∂ y ∂ x = ∂ 2 f / ∂ x ∂ y 。恰当微分的线积分只依赖于端点,与路径无关。(详见附录 A2 ,特别是 §A2.2 –§A2.3 。)
例 :以 ( T , V ) (T,V) ( T , V ) 为自变量,理想气体的 δ W = − p d V = 0 ⋅ d T + ( − N k B T / V ) d V \delta W = -p\,dV = 0\cdot dT + \left(-Nk_{\mathrm B}T/V\right)dV δ W = − p d V = 0 ⋅ d T + ( − N k B T / V ) d V 。这里 M = 0 M = 0 M = 0 ,L = − N k B T / V L = -Nk_{\mathrm B}T/V L = − N k B T / V ,而 ∂ M / ∂ V = 0 ≠ ∂ L / ∂ T = − N k B / V \partial M/\partial V = 0 \ne \partial L/\partial T = -Nk_{\mathrm B}/V ∂ M / ∂ V = 0 = ∂ L / ∂ T = − N k B / V ,所以 δ W \delta W δ W 不是恰当微分——功与路径有关。
热容 。定容热容 C V ≡ ( δ Q / d T ) V C_V \equiv (\delta Q/dT)_V C V ≡ ( δ Q / d T ) V 。体积不变时 δ W = 0 \delta W = 0 δ W = 0 ,由 (1.1) 得 δ Q = d U \delta Q = dU δ Q = d U ,所以
C V = ( ∂ U ∂ T ) V (1.2) C_V = \left(\frac{\partial U}{\partial T}\right)_V \tag{1.2} C V = ( ∂ T ∂ U ) V ( 1.2 )
定义焓 H ≡ U + p V H \equiv U + pV H ≡ U + p V 。准静态过程中 d H = d U + p d V + V d p = ( δ Q − p d V ) + p d V + V d p = δ Q + V d p dH = dU + p\,dV + V\,dp = (\delta Q - p\,dV) + p\,dV + V\,dp = \delta Q + V\,dp d H = d U + p d V + V d p = ( δ Q − p d V ) + p d V + V d p = δ Q + V d p ;压强不变时 δ Q = d H \delta Q = dH δ Q = d H ,所以
C p = ( ∂ H ∂ T ) p (1.3) C_p = \left(\frac{\partial H}{\partial T}\right)_p \tag{1.3} C p = ( ∂ T ∂ H ) p ( 1.3 )
§1.4 第二定律与熵
第二定律(克劳修斯表述):热量不能自发地从低温物体传到高温物体。在热力学中,由它(借助卡诺循环的论证)得到以下两条结论。本讲义把它们作为第二定律的内容直接接受;第3 、4章 将从微观上重新得到它们。
(1) 熵的存在 :对可逆过程,δ Q r e v / T \delta Q_{\mathrm{rev}}/T δ Q rev / T 是某个状态函数——熵 S S S ——的全微分:
d S = δ Q r e v T (1.4) dS = \frac{\delta Q_{\mathrm{rev}}}{T} \tag{1.4} d S = T δ Q rev ( 1.4 )
用提示 A2 的语言:δ Q r e v \delta Q_{\mathrm{rev}} δ Q rev 本身不是恰当微分,但乘上"积分因子" 1 / T 1/T 1/ T 后成为恰当微分。
(2) 熵增原理 :对任意过程有克劳修斯不等式 d S ≥ δ Q / T dS \ge \delta Q/T d S ≥ δ Q / T (等号对应可逆过程)。对孤立系统 δ Q = 0 \delta Q = 0 δ Q = 0 ,故
d S ≥ 0 ( 孤立系统 ) (1.5) dS \ge 0 \qquad(\text{孤立系统}) \tag{1.5} d S ≥ 0 ( 孤立系统 ) ( 1.5 )
孤立系统的熵永不减少;在给定约束下,平衡态是熵取极大值的状态 。第3章 将看到它的微观含义:平衡态是包含微观态数目最多的宏观态。
热力学基本方程 。对封闭系统的可逆过程,δ Q = T d S \delta Q = T\,dS δ Q = T d S ,δ W = − p d V \delta W = -p\,dV δ W = − p d V ,代入 (1.1) 得 d U = T d S − p d V dU = T\,dS - p\,dV d U = T d S − p d V 。由于此式只含状态函数,它对连接两个邻近平衡态的任何过程都成立。若粒子数也可变,则
d U = T d S − p d V + μ d N (1.6) dU = T\,dS - p\,dV + \mu\,dN \tag{1.6} d U = T d S − p d V + μ d N ( 1.6 )
其中 μ ≡ ( ∂ U / ∂ N ) S , V \mu \equiv (\partial U/\partial N)_{S,V} μ ≡ ( ∂ U / ∂ N ) S , V 称为化学势 。(1.6) 说明 U U U 的"自然变量"是 ( S , V , N ) (S,V,N) ( S , V , N ) 。把 (1.6) 解出 d S dS d S ,得到熵表象 的基本方程:
d S = 1 T d U + p T d V − μ T d N (1.7) dS = \frac{1}{T}\,dU + \frac{p}{T}\,dV - \frac{\mu}{T}\,dN \tag{1.7} d S = T 1 d U + T p d V − T μ d N ( 1.7 )
偏导数与全微分
( ∂ S / ∂ U ) V , N (\partial S/\partial U)_{V,N} ( ∂ S / ∂ U ) V , N 表示固定 V , N V,N V , N 、只让 U U U 变化时 S S S 的变化率。热力学中下标注明"哪些变量保持不变",这一点至关重要:例如固定 V V V 与固定 p p p 时的 ∂ S / ∂ T \partial S/\partial T ∂ S / ∂ T 是两个不同的量。若 S = S ( U , V , N ) S = S(U,V,N) S = S ( U , V , N ) ,则 d S = ( ∂ S / ∂ U ) V , N d U + ( ∂ S / ∂ V ) U , N d V + ( ∂ S / ∂ N ) U , V d N dS = (\partial S/\partial U)_{V,N}dU + (\partial S/\partial V)_{U,N}dV + (\partial S/\partial N)_{U,V}dN d S = ( ∂ S / ∂ U ) V , N d U + ( ∂ S / ∂ V ) U , N d V + ( ∂ S / ∂ N ) U , V d N 。(详见附录 A1 。)
把这个全微分与 (1.7) 逐项比较:
( ∂ S ∂ U ) V , N = 1 T , ( ∂ S ∂ V ) U , N = p T , ( ∂ S ∂ N ) U , V = − μ T (1.8) \left(\frac{\partial S}{\partial U}\right)_{V,N} = \frac{1}{T},\qquad \left(\frac{\partial S}{\partial V}\right)_{U,N} = \frac{p}{T},\qquad \left(\frac{\partial S}{\partial N}\right)_{U,V} = -\frac{\mu}{T} \tag{1.8} ( ∂ U ∂ S ) V , N = T 1 , ( ∂ V ∂ S ) U , N = T p , ( ∂ N ∂ S ) U , V = − T μ ( 1.8 )
要点
(1.8) 是统计力学与热力学"接头"的地方:微正则系综(第3章 )先从微观上算出 S ( U , V , N ) S(U,V,N) S ( U , V , N ) ,再由 (1.8) 读出 T T T 、p p p 、μ \mu μ 。
§1.5 平衡条件与稳定性条件
考虑一个孤立系统,由子系统 1、2 组成,二者可以交换能量、体积和粒子:U 1 + U 2 = U U_1+U_2 = U U 1 + U 2 = U ,V 1 + V 2 = V V_1+V_2 = V V 1 + V 2 = V ,N 1 + N 2 = N N_1+N_2 = N N 1 + N 2 = N 均固定。总熵
S = S 1 ( U 1 , V 1 , N 1 ) + S 2 ( U − U 1 , V − V 1 , N − N 1 ) S = S_1(U_1,V_1,N_1) + S_2(U-U_1,\ V-V_1,\ N-N_1) S = S 1 ( U 1 , V 1 , N 1 ) + S 2 ( U − U 1 , V − V 1 , N − N 1 )
利用 (1.8) 和链式法则(注意 d U 2 = − d U 1 dU_2 = -dU_1 d U 2 = − d U 1 等):
d S = ( 1 T 1 − 1 T 2 ) d U 1 + ( p 1 T 1 − p 2 T 2 ) d V 1 − ( μ 1 T 1 − μ 2 T 2 ) d N 1 dS = \left(\frac{1}{T_1}-\frac{1}{T_2}\right)dU_1 + \left(\frac{p_1}{T_1}-\frac{p_2}{T_2}\right)dV_1 - \left(\frac{\mu_1}{T_1}-\frac{\mu_2}{T_2}\right)dN_1 d S = ( T 1 1 − T 2 1 ) d U 1 + ( T 1 p 1 − T 2 p 2 ) d V 1 − ( T 1 μ 1 − T 2 μ 2 ) d N 1
平衡时 S S S 取极大,对任意的 d U 1 , d V 1 , d N 1 dU_1, dV_1, dN_1 d U 1 , d V 1 , d N 1 都有 d S = 0 dS = 0 d S = 0 ,因此三个括号分别为零:
T 1 = T 2 , p 1 = p 2 , μ 1 = μ 2 (1.9) T_1 = T_2,\qquad p_1 = p_2,\qquad \mu_1 = \mu_2 \tag{1.9} T 1 = T 2 , p 1 = p 2 , μ 1 = μ 2 ( 1.9 )
分别称为热平衡、力学平衡与化学平衡条件。
趋向平衡的方向 。若只交换能量且 T 1 < T 2 T_1<T_2 T 1 < T 2 ,熵增要求 d S = ( 1 / T 1 − 1 / T 2 ) d U 1 > 0 dS = (1/T_1 - 1/T_2)\,dU_1 > 0 d S = ( 1/ T 1 − 1/ T 2 ) d U 1 > 0 ;括号为正,故 d U 1 > 0 dU_1>0 d U 1 > 0 :能量(热)从高温流向低温。若温度相同、只交换粒子且 μ 1 > μ 2 \mu_1 > \mu_2 μ 1 > μ 2 ,则 d S = − ( μ 1 − μ 2 ) d N 1 / T > 0 dS = -(\mu_1-\mu_2)\,dN_1/T > 0 d S = − ( μ 1 − μ 2 ) d N 1 / T > 0 要求 d N 1 < 0 dN_1 < 0 d N 1 < 0 :粒子从化学势高处流向化学势低处。
稳定性 。极大值还要求二阶变化为负。只考虑能量交换,并取两个完全相同的子系统,令 U 1 = U / 2 + Δ U_1 = U/2 + \Delta U 1 = U /2 + Δ ,U 2 = U / 2 − Δ U_2 = U/2 - \Delta U 2 = U /2 − Δ 。
泰勒展开
f ( x 0 + Δ ) = f ( x 0 ) + f ′ ( x 0 ) Δ + 1 2 f ′ ′ ( x 0 ) Δ 2 + ⋯ f(x_0+\Delta) = f(x_0) + f'(x_0)\Delta + \tfrac12 f''(x_0)\Delta^2 + \cdots f ( x 0 + Δ ) = f ( x 0 ) + f ′ ( x 0 ) Δ + 2 1 f ′′ ( x 0 ) Δ 2 + ⋯ ;多元函数类似,二阶项含各二阶偏导数。(详见附录 A3 ,特别是 §A3.1 、§A3.4 。)
S ( U / 2 + Δ ) + S ( U / 2 − Δ ) = 2 S ( U / 2 ) + ∂ 2 S ∂ U 2 Δ 2 + ⋯ S(U/2+\Delta) + S(U/2-\Delta) = 2S(U/2) + \frac{\partial^2 S}{\partial U^2}\Delta^2 + \cdots S ( U /2 + Δ ) + S ( U /2 − Δ ) = 2 S ( U /2 ) + ∂ U 2 ∂ 2 S Δ 2 + ⋯
(一阶项相互抵消。)Δ = 0 \Delta = 0 Δ = 0 为极大要求 ∂ 2 S / ∂ U 2 < 0 \partial^2 S/\partial U^2 < 0 ∂ 2 S / ∂ U 2 < 0 。而由 (1.8) 、(1.2) :
∂ 2 S ∂ U 2 = ∂ ∂ U ( 1 T ) = − 1 T 2 ( ∂ T ∂ U ) V = − 1 T 2 C V , 故 C V > 0 (1.10) \frac{\partial^2 S}{\partial U^2} = \frac{\partial}{\partial U}\left(\frac1T\right) = -\frac{1}{T^2}\left(\frac{\partial T}{\partial U}\right)_V = -\frac{1}{T^2 C_V}, \qquad\text{故}\qquad C_V > 0 \tag{1.10} ∂ U 2 ∂ 2 S = ∂ U ∂ ( T 1 ) = − T 2 1 ( ∂ U ∂ T ) V = − T 2 C V 1 , 故 C V > 0 ( 1.10 )
§4.4 将从统计力学给出这一结论的另一证明:C V C_V C V 正比于能量涨落的方差。
§1.6 热力学势与勒让德变换
实验上容易控制的通常是温度而不是熵。我们希望把自变量 S S S 换成 T = ( ∂ U / ∂ S ) V , N T = (\partial U/\partial S)_{V,N} T = ( ∂ U / ∂ S ) V , N ,同时不丢失任何信息。这正是勒让德变换所做的事。
勒让德变换
设 f ( x ) f(x) f ( x ) 的斜率 y = f ′ ( x ) y = f'(x) y = f ′ ( x ) 随 x x x 单调变化,因而可反解出 x = x ( y ) x = x(y) x = x ( y ) 。定义 g ( y ) ≡ f ( x ) − x y g(y) \equiv f(x) - xy g ( y ) ≡ f ( x ) − x y (右边的 x x x 用 x ( y ) x(y) x ( y ) 代入)。则 d g = f ′ ( x ) d x − x d y − y d x = − x d y dg = f'(x)\,dx - x\,dy - y\,dx = -x\,dy d g = f ′ ( x ) d x − x d y − y d x = − x d y ,即 g ′ ( y ) = − x g'(y) = -x g ′ ( y ) = − x 。g ( y ) g(y) g ( y ) 与 f ( x ) f(x) f ( x ) 包含相同的信息,但以斜率 y y y 为自变量。多元函数可对需要替换的变量逐个进行。(详见附录 A5 。)
亥姆霍兹自由能 F ≡ U − T S F \equiv U - TS F ≡ U − TS :
d F = d U − T d S − S d T = − S d T − p d V + μ d N (1.11) dF = dU - T\,dS - S\,dT = -S\,dT - p\,dV + \mu\,dN \tag{1.11} d F = d U − T d S − S d T = − S d T − p d V + μ d N ( 1.11 )
(第二个等号代入了 (1.6) 。)同理可得:
H ≡ U + p V : d H = T d S + V d p + μ d N (1.12) H \equiv U + pV:\qquad dH = T\,dS + V\,dp + \mu\,dN \tag{1.12} H ≡ U + p V : d H = T d S + V d p + μ d N ( 1.12 )
G ≡ U − T S + p V : d G = − S d T + V d p + μ d N (1.13) G \equiv U - TS + pV:\qquad dG = -S\,dT + V\,dp + \mu\,dN \tag{1.13} G ≡ U − TS + p V : d G = − S d T + V d p + μ d N ( 1.13 )
J ≡ U − T S − μ N : d J = − S d T − p d V − N d μ (1.14) J \equiv U - TS - \mu N:\qquad dJ = -S\,dT - p\,dV - N\,d\mu \tag{1.14} J ≡ U − TS − μ N : dJ = − S d T − p d V − N d μ ( 1.14 )
从这些全微分可以直接读出偏导数,例如
S = − ( ∂ F ∂ T ) V , N , p = − ( ∂ F ∂ V ) T , N , μ = ( ∂ F ∂ N ) T , V , N = − ( ∂ J ∂ μ ) T , V S = -\left(\frac{\partial F}{\partial T}\right)_{V,N},\quad p = -\left(\frac{\partial F}{\partial V}\right)_{T,N},\quad \mu = \left(\frac{\partial F}{\partial N}\right)_{T,V},\quad N = -\left(\frac{\partial J}{\partial \mu}\right)_{T,V} S = − ( ∂ T ∂ F ) V , N , p = − ( ∂ V ∂ F ) T , N , μ = ( ∂ N ∂ F ) T , V , N = − ( ∂ μ ∂ J ) T , V
吉布斯–亥姆霍兹关系 (第4章 要用)。由 U = F + T S = F − T ( ∂ F / ∂ T ) V U = F + TS = F - T(\partial F/\partial T)_V U = F + TS = F − T ( ∂ F / ∂ T ) V ,以及
∂ ∂ T ( F T ) = 1 T ∂ F ∂ T − F T 2 = − T S + F T 2 = − U T 2 \frac{\partial}{\partial T}\left(\frac FT\right) = \frac1T\frac{\partial F}{\partial T} - \frac{F}{T^2} = -\frac{TS + F}{T^2} = -\frac{U}{T^2} ∂ T ∂ ( T F ) = T 1 ∂ T ∂ F − T 2 F = − T 2 TS + F = − T 2 U
得 U = − T 2 ∂ ( F / T ) / ∂ T U = -T^2\,\partial(F/T)/\partial T U = − T 2 ∂ ( F / T ) / ∂ T 。改用 β = 1 / ( k B T ) \beta = 1/(k_{\mathrm B}T) β = 1/ ( k B T ) 为变量:∂ / ∂ β = ( d T / d β ) ∂ / ∂ T = − k B T 2 ∂ / ∂ T \partial/\partial\beta = (dT/d\beta)\,\partial/\partial T = -k_{\mathrm B}T^2\,\partial/\partial T ∂ / ∂ β = ( d T / d β ) ∂ / ∂ T = − k B T 2 ∂ / ∂ T ,而 β F = F / ( k B T ) \beta F = F/(k_{\mathrm B}T) βF = F / ( k B T ) ,于是
U = ( ∂ ( β F ) ∂ β ) V , N (1.15) U = \left(\frac{\partial(\beta F)}{\partial\beta}\right)_{V,N} \tag{1.15} U = ( ∂ β ∂ ( βF ) ) V , N ( 1.15 )
极小原理 。设系统与温度为 T T T 的大热源接触,V V V 、N N N 固定。系统 + 热源构成孤立系统,由 (1.5) ,d S + d S R ≥ 0 dS + dS_{\mathrm R} \ge 0 d S + d S R ≥ 0 。热源温度恒定且内部过程可逆,它吸收的热为 − δ Q -\delta Q − δ Q (δ Q \delta Q δ Q 是系统吸收的热),故 d S R = − δ Q / T dS_{\mathrm R} = -\delta Q/T d S R = − δ Q / T ;体积不变时 δ Q = d U \delta Q = dU δ Q = d U 。于是 d S − d U / T ≥ 0 dS - dU/T \ge 0 d S − d U / T ≥ 0 ,即(T T T 不变)
d F = d ( U − T S ) ≤ 0 (1.16) dF = d(U - TS) \le 0 \tag{1.16} d F = d ( U − TS ) ≤ 0 ( 1.16 )
F F F 只减不增,在 T , V , N T,V,N T , V , N 固定时平衡态使 F F F 极小 。同理:若系统还与压强为 p p p 的外界交换体积,d S R = − ( d U + p d V ) / T dS_{\mathrm R} = -(dU + p\,dV)/T d S R = − ( d U + p d V ) / T ,得 d ( U − T S + p V ) ≤ 0 d(U - TS + pV) \le 0 d ( U − TS + p V ) ≤ 0 ,即 T , p , N T,p,N T , p , N 固定时 G G G 极小;若与源交换能量和粒子(化学势 μ \mu μ ),d S R = − ( d U − μ d N ) / T dS_{\mathrm R} = -(dU - \mu\,dN)/T d S R = − ( d U − μ d N ) / T ,得 T , V , μ T,V,\mu T , V , μ 固定时 J J J 极小。
§1.7 麦克斯韦关系
由 (1.11) ,F ( T , V ) F(T,V) F ( T , V ) 的两个一阶偏导数是 − S -S − S 与 − p -p − p 。由混合偏导数相等(提示 A2 ),∂ ∂ V ∂ F ∂ T = ∂ ∂ T ∂ F ∂ V \dfrac{\partial}{\partial V}\dfrac{\partial F}{\partial T} = \dfrac{\partial}{\partial T}\dfrac{\partial F}{\partial V} ∂ V ∂ ∂ T ∂ F = ∂ T ∂ ∂ V ∂ F ,即 ( ∂ S / ∂ V ) T = ( ∂ p / ∂ T ) V (\partial S/\partial V)_T = (\partial p/\partial T)_V ( ∂ S / ∂ V ) T = ( ∂ p / ∂ T ) V 。对 G G G 、U U U 、H H H 做同样的事,得到四个麦克斯韦关系(N N N 固定,下标中省略):
( ∂ S ∂ V ) T = ( ∂ p ∂ T ) V , ( ∂ S ∂ p ) T = − ( ∂ V ∂ T ) p , ( ∂ T ∂ V ) S = − ( ∂ p ∂ S ) V , ( ∂ T ∂ p ) S = ( ∂ V ∂ S ) p (1.17) \left(\frac{\partial S}{\partial V}\right)_T = \left(\frac{\partial p}{\partial T}\right)_V,\quad \left(\frac{\partial S}{\partial p}\right)_T = -\left(\frac{\partial V}{\partial T}\right)_p,\quad \left(\frac{\partial T}{\partial V}\right)_S = -\left(\frac{\partial p}{\partial S}\right)_V,\quad \left(\frac{\partial T}{\partial p}\right)_S = \left(\frac{\partial V}{\partial S}\right)_p \tag{1.17} ( ∂ V ∂ S ) T = ( ∂ T ∂ p ) V , ( ∂ p ∂ S ) T = − ( ∂ T ∂ V ) p , ( ∂ V ∂ T ) S = − ( ∂ S ∂ p ) V , ( ∂ p ∂ T ) S = ( ∂ S ∂ V ) p ( 1.17 )
它们的用处是把难以直接测量的熵的导数换成可测的 p p p –V V V –T T T 关系。
§1.8 响应函数与偏导数技巧
由于 δ Q r e v = T d S \delta Q_{\mathrm{rev}} = T\,dS δ Q rev = T d S ,热容也可写成
C V = T ( ∂ S ∂ T ) V , C p = T ( ∂ S ∂ T ) p (1.18) C_V = T\left(\frac{\partial S}{\partial T}\right)_V,\qquad C_p = T\left(\frac{\partial S}{\partial T}\right)_p \tag{1.18} C V = T ( ∂ T ∂ S ) V , C p = T ( ∂ T ∂ S ) p ( 1.18 )
另定义体膨胀系数和等温压缩系数:
α ≡ 1 V ( ∂ V ∂ T ) p , κ T ≡ − 1 V ( ∂ V ∂ p ) T (1.19) \alpha \equiv \frac1V\left(\frac{\partial V}{\partial T}\right)_p,\qquad \kappa_T \equiv -\frac1V\left(\frac{\partial V}{\partial p}\right)_T \tag{1.19} α ≡ V 1 ( ∂ T ∂ V ) p , κ T ≡ − V 1 ( ∂ p ∂ V ) T ( 1.19 )
偏导数恒等式
若 x , y , z x,y,z x , y , z 三个变量之间有一个关系(只有两个独立),则:(i) 倒数关系 ( ∂ x / ∂ y ) z = 1 / ( ∂ y / ∂ x ) z (\partial x/\partial y)_z = 1/(\partial y/\partial x)_z ( ∂ x / ∂ y ) z = 1/ ( ∂ y / ∂ x ) z ;(ii) 循环关系 ( ∂ x / ∂ y ) z ( ∂ y / ∂ z ) x ( ∂ z / ∂ x ) y = − 1 (\partial x/\partial y)_z(\partial y/\partial z)_x(\partial z/\partial x)_y = -1 ( ∂ x / ∂ y ) z ( ∂ y / ∂ z ) x ( ∂ z / ∂ x ) y = − 1 ;(iii) 换固定变量:若 f = f ( x , y ) f = f(x,y) f = f ( x , y ) 而 y y y 又可看作 ( x , z ) (x,z) ( x , z ) 的函数,则 ( ∂ f / ∂ x ) z = ( ∂ f / ∂ x ) y + ( ∂ f / ∂ y ) x ( ∂ y / ∂ x ) z (\partial f/\partial x)_z = (\partial f/\partial x)_y + (\partial f/\partial y)_x(\partial y/\partial x)_z ( ∂ f / ∂ x ) z = ( ∂ f / ∂ x ) y + ( ∂ f / ∂ y ) x ( ∂ y / ∂ x ) z 。(详见附录 A4 。)
例:C p − C V C_p - C_V C p − C V 。把 S S S 看成 ( T , V ) (T,V) ( T , V ) 的函数,再把 V V V 看成 ( T , p ) (T,p) ( T , p ) 的函数,由 (iii):
( ∂ S ∂ T ) p = ( ∂ S ∂ T ) V + ( ∂ S ∂ V ) T ( ∂ V ∂ T ) p \left(\frac{\partial S}{\partial T}\right)_p = \left(\frac{\partial S}{\partial T}\right)_V + \left(\frac{\partial S}{\partial V}\right)_T\left(\frac{\partial V}{\partial T}\right)_p ( ∂ T ∂ S ) p = ( ∂ T ∂ S ) V + ( ∂ V ∂ S ) T ( ∂ T ∂ V ) p
两边乘以 T T T ,用 (1.18) 及麦克斯韦关系 ( ∂ S / ∂ V ) T = ( ∂ p / ∂ T ) V (\partial S/\partial V)_T = (\partial p/\partial T)_V ( ∂ S / ∂ V ) T = ( ∂ p / ∂ T ) V :C p = C V + T ( ∂ p / ∂ T ) V ( ∂ V / ∂ T ) p C_p = C_V + T(\partial p/\partial T)_V(\partial V/\partial T)_p C p = C V + T ( ∂ p / ∂ T ) V ( ∂ V / ∂ T ) p 。再由循环关系与倒数关系,
( ∂ p ∂ T ) V = − ( ∂ V / ∂ T ) p ( ∂ V / ∂ p ) T = V α V κ T = α κ T \left(\frac{\partial p}{\partial T}\right)_V = -\frac{(\partial V/\partial T)_p}{(\partial V/\partial p)_T} = \frac{V\alpha}{V\kappa_T} = \frac{\alpha}{\kappa_T} ( ∂ T ∂ p ) V = − ( ∂ V / ∂ p ) T ( ∂ V / ∂ T ) p = V κ T V α = κ T α
于是
C p − C V = T V α 2 κ T (1.20) C_p - C_V = \frac{TV\alpha^2}{\kappa_T} \tag{1.20} C p − C V = κ T T V α 2 ( 1.20 )
对理想气体 V = N k B T / p V = Nk_{\mathrm B}T/p V = N k B T / p ,α = 1 / T \alpha = 1/T α = 1/ T ,κ T = 1 / p \kappa_T = 1/p κ T = 1/ p ,得 C p − C V = p V / T = N k B C_p - C_V = pV/T = Nk_{\mathrm B} C p − C V = p V / T = N k B (迈耶关系)。
能态方程 。在 (1.6) 中固定 T T T 、N N N ,两边除以 d V dV d V ,并用麦克斯韦关系:
( ∂ U ∂ V ) T = T ( ∂ S ∂ V ) T − p = T ( ∂ p ∂ T ) V − p (1.21) \left(\frac{\partial U}{\partial V}\right)_T = T\left(\frac{\partial S}{\partial V}\right)_T - p = T\left(\frac{\partial p}{\partial T}\right)_V - p \tag{1.21} ( ∂ V ∂ U ) T = T ( ∂ V ∂ S ) T − p = T ( ∂ T ∂ p ) V − p ( 1.21 )
对理想气体右边 = T ⋅ N k B / V − p = 0 = T\cdot Nk_{\mathrm B}/V - p = 0 = T ⋅ N k B / V − p = 0 :理想气体的内能只是温度的函数。
κ T > 0 \kappa_T > 0 κ T > 0 。在 T , N T,N T , N 固定下把系统等分为两半,允许它们交换体积:F 总 = F ( V / 2 + Δ ) + F ( V / 2 − Δ ) = 2 F ( V / 2 ) + ( ∂ 2 F / ∂ V 2 ) Δ 2 + ⋯ F_{\text{总}} = F(V/2+\Delta) + F(V/2-\Delta) = 2F(V/2) + (\partial^2F/\partial V^2)\Delta^2 + \cdots F 总 = F ( V /2 + Δ ) + F ( V /2 − Δ ) = 2 F ( V /2 ) + ( ∂ 2 F / ∂ V 2 ) Δ 2 + ⋯ 。由 (1.16) ,平衡要求 F F F 极小,故 ∂ 2 F / ∂ V 2 > 0 \partial^2F/\partial V^2 > 0 ∂ 2 F / ∂ V 2 > 0 。而 ∂ 2 F / ∂ V 2 = − ( ∂ p / ∂ V ) T = 1 / ( V κ T ) \partial^2 F/\partial V^2 = -(\partial p/\partial V)_T = 1/(V\kappa_T) ∂ 2 F / ∂ V 2 = − ( ∂ p / ∂ V ) T = 1/ ( V κ T ) ,所以
κ T > 0 (1.22) \kappa_T > 0 \tag{1.22} κ T > 0 ( 1.22 )
§1.9 广延性:欧拉关系与吉布斯–杜安关系
U U U 是广延量,其自然变量 S , V , N S,V,N S , V , N 也都是广延量,因此 U ( λ S , λ V , λ N ) = λ U ( S , V , N ) U(\lambda S,\lambda V,\lambda N) = \lambda U(S,V,N) U ( λ S , λV , λ N ) = λ U ( S , V , N ) 。
齐次函数的欧拉定理
若 f ( λ x 1 , … , λ x n ) = λ k f ( x 1 , … , x n ) f(\lambda x_1,\dots,\lambda x_n) = \lambda^k f(x_1,\dots,x_n) f ( λ x 1 , … , λ x n ) = λ k f ( x 1 , … , x n ) 对一切 λ > 0 \lambda>0 λ > 0 成立,两边对 λ \lambda λ 求导(左边用链式法则)后令 λ = 1 \lambda = 1 λ = 1 ,得 ∑ i x i ∂ f / ∂ x i = k f \sum_i x_i\,\partial f/\partial x_i = k f ∑ i x i ∂ f / ∂ x i = k f 。(详见附录 A9 ,特别是 §A9.2 –§A9.4 。)
取 k = 1 k=1 k = 1 ,并用 ∂ U / ∂ S = T \partial U/\partial S = T ∂ U / ∂ S = T 等:
U = T S − p V + μ N (1.23) U = TS - pV + \mu N \tag{1.23} U = TS − p V + μ N ( 1.23 )
于是
G = U − T S + p V = μ N , J = U − T S − μ N = − p V (1.24) G = U - TS + pV = \mu N,\qquad J = U - TS - \mu N = -pV \tag{1.24} G = U − TS + p V = μ N , J = U − TS − μ N = − p V ( 1.24 )
即化学势是每个粒子的吉布斯自由能,巨势等于 − p V -pV − p V (第8章 要用)。对 (1.23) 求全微分并减去 (1.6) :
S d T − V d p + N d μ = 0 (1.25) S\,dT - V\,dp + N\,d\mu = 0 \tag{1.25} S d T − V d p + N d μ = 0 ( 1.25 )
这是吉布斯–杜安关系 :T T T 、p p p 、μ \mu μ 三个强度量不是独立的。写成每个粒子的量(s = S / N s = S/N s = S / N ,v = V / N v = V/N v = V / N ):d μ = − s d T + v d p d\mu = -s\,dT + v\,dp d μ = − s d T + v d p 。
§1.10 第三定律
能斯特–普朗克表述:当 T → 0 T\to 0 T → 0 时,系统的熵趋于一个常数 S 0 S_0 S 0 (对完美晶体可取 S 0 = 0 S_0 = 0 S 0 = 0 )。一个直接推论:由 (1.18) ,定容下 d S = C V d T / T dS = C_V\,dT/T d S = C V d T / T ,所以
S ( T ) − S ( 0 ) = ∫ 0 T C V ( T ′ ) T ′ d T ′ S(T) - S(0) = \int_0^T \frac{C_V(T')}{T'}\,dT' S ( T ) − S ( 0 ) = ∫ 0 T T ′ C V ( T ′ ) d T ′
若 T → 0 T\to 0 T → 0 时 C V C_V C V 趋于非零常数 c c c ,积分在下限附近表现为 c ln T ′ c\ln T' c ln T ′ 而发散,与 S ( 0 ) S(0) S ( 0 ) 有限矛盾。因此
C V → 0 ( T → 0 ) (1.26) C_V \to 0 \qquad (T\to 0) \tag{1.26} C V → 0 ( T → 0 ) ( 1.26 )
(C p C_p C p 同理。)§3.3 将给出第三定律的统计解释(基态简并度)。
§1.11 理想气体与化学平衡速查
(a) 理想气体的熵 。以 ( T , V ) (T,V) ( T , V ) 为变量,d S = ( ∂ S / ∂ T ) V d T + ( ∂ S / ∂ V ) T d V = ( C V / T ) d T + ( ∂ p / ∂ T ) V d V = ( C V / T ) d T + ( N k B / V ) d V dS = (\partial S/\partial T)_V dT + (\partial S/\partial V)_T dV = (C_V/T)\,dT + (\partial p/\partial T)_V\,dV = (C_V/T)\,dT + (Nk_{\mathrm B}/V)\,dV d S = ( ∂ S / ∂ T ) V d T + ( ∂ S / ∂ V ) T d V = ( C V / T ) d T + ( ∂ p / ∂ T ) V d V = ( C V / T ) d T + ( N k B / V ) d V 。设 C V = N c v C_V = Nc_v C V = N c v 为常数,积分得 S = N c v ln T + N k B ln V + f ( N ) S = Nc_v\ln T + Nk_{\mathrm B}\ln V + f(N) S = N c v ln T + N k B ln V + f ( N ) ,其中 f ( N ) f(N) f ( N ) 是积分"常数"(与 T T T 、V V V 无关,但可以依赖于 N N N )。要求 S S S 是广延量:S ( T , λ V , λ N ) = λ S ( T , V , N ) S(T,\lambda V,\lambda N) = \lambda S(T,V,N) S ( T , λV , λ N ) = λ S ( T , V , N ) ,比较两边得 f ( λ N ) = λ f ( N ) − λ N k B ln λ f(\lambda N) = \lambda f(N) - \lambda N k_{\mathrm B}\ln\lambda f ( λ N ) = λ f ( N ) − λ N k B ln λ 。令 N = 1 N=1 N = 1 、再把 λ \lambda λ 改记为 N N N :f ( N ) = N f ( 1 ) − N k B ln N f(N) = Nf(1) - Nk_{\mathrm B}\ln N f ( N ) = N f ( 1 ) − N k B ln N 。记 s 0 = f ( 1 ) s_0 = f(1) s 0 = f ( 1 ) :
S = N [ c v ln T + k B ln V N + s 0 ] (1.27) S = N\left[c_v\ln T + k_{\mathrm B}\ln\frac{V}{N} + s_0\right] \tag{1.27} S = N [ c v ln T + k B ln N V + s 0 ] ( 1.27 )
热力学无法 确定常数 s 0 s_0 s 0 ;统计力学将把它算出来(萨克尔–特多德公式,§3.6 )。
(b) 理想气体的化学势 。由吉布斯–杜安关系,等温时 d μ = v d p = ( k B T / p ) d p d\mu = v\,dp = (k_{\mathrm B}T/p)\,dp d μ = v d p = ( k B T / p ) d p ,从标准压强 p ∘ p^\circ p ∘ (通常取 1 bar)积分到 p p p :
μ ( T , p ) = μ ∘ ( T ) + k B T ln p p ∘ (1.28) \mu(T,p) = \mu^\circ(T) + k_{\mathrm B}T\ln\frac{p}{p^\circ} \tag{1.28} μ ( T , p ) = μ ∘ ( T ) + k B T ln p ∘ p ( 1.28 )
对理想气体混合物,每个组分 i i i 的化学势具有同样的形式,只是 p p p 换成分压 p i p_i p i :μ i = μ i ∘ ( T ) + k B T ln ( p i / p ∘ ) \mu_i = \mu_i^\circ(T) + k_{\mathrm B}T\ln(p_i/p^\circ) μ i = μ i ∘ ( T ) + k B T ln ( p i / p ∘ ) 。这一点将在 §7.8 由统计力学证明。
(c) 化学平衡 。把化学反应写成 ∑ i ν i A i = 0 \sum_i \nu_i A_i = 0 ∑ i ν i A i = 0 ,其中 ν i \nu_i ν i 为化学计量数(生成物取正、反应物取负)。若反应进度改变 d ξ d\xi d ξ ,则 d N i = ν i d ξ dN_i = \nu_i\,d\xi d N i = ν i d ξ 。在 T T T 、p p p 固定时,由 (1.13) 的多组分形式 d G = ∑ i μ i d N i = ( ∑ i ν i μ i ) d ξ dG = \sum_i\mu_i\,dN_i = \left(\sum_i\nu_i\mu_i\right)d\xi d G = ∑ i μ i d N i = ( ∑ i ν i μ i ) d ξ 。G G G 极小要求
∑ i ν i μ i = 0 (1.29) \sum_i \nu_i\mu_i = 0 \tag{1.29} i ∑ ν i μ i = 0 ( 1.29 )
对理想气体代入 (1.28) :∑ i ν i μ i ∘ + k B T ∑ i ν i ln ( p i / p ∘ ) = 0 \sum_i\nu_i\mu_i^\circ + k_{\mathrm B}T\sum_i \nu_i\ln(p_i/p^\circ) = 0 ∑ i ν i μ i ∘ + k B T ∑ i ν i ln ( p i / p ∘ ) = 0 。利用 ∑ i ν i ln x i = ln ∏ i x i ν i \sum_i \nu_i \ln x_i = \ln\prod_i x_i^{\nu_i} ∑ i ν i ln x i = ln ∏ i x i ν i :
K ∘ ( T ) ≡ ∏ i ( p i p ∘ ) ν i = exp ( − ∑ i ν i μ i ∘ k B T ) (1.30) K^\circ(T) \equiv \prod_i\left(\frac{p_i}{p^\circ}\right)^{\nu_i} = \exp\left(-\frac{\sum_i\nu_i\mu_i^\circ}{k_{\mathrm B}T}\right) \tag{1.30} K ∘ ( T ) ≡ i ∏ ( p ∘ p i ) ν i = exp ( − k B T ∑ i ν i μ i ∘ ) ( 1.30 )
按摩尔计时把 k B k_{\mathrm B} k B 换成 R R R ,指数中即为 − Δ r G ∘ / R T -\Delta_{\mathrm r}G^\circ/RT − Δ r G ∘ / RT 。§7.8 将从分子配分函数计算 K ∘ K^\circ K ∘ 。
§1.12 本章小结
热力学势 定义 全微分 自然变量 平衡判据 对应的系综 S S S — d S = d U / T + p d V / T − μ d N / T dS = dU/T + p\,dV/T - \mu\,dN/T d S = d U / T + p d V / T − μ d N / T U , V , N U,V,N U , V , N 孤立系统中极大 微正则(第3章 ) F F F U − T S U-TS U − TS d F = − S d T − p d V + μ d N dF = -S\,dT - p\,dV + \mu\,dN d F = − S d T − p d V + μ d N T , V , N T,V,N T , V , N 极小 正则(第4章 ) G G G U − T S + p V = μ N U-TS+pV = \mu N U − TS + p V = μ N d G = − S d T + V d p + μ d N dG = -S\,dT + V\,dp + \mu\,dN d G = − S d T + V d p + μ d N T , p , N T,p,N T , p , N 极小 等温等压系综(本讲义不展开) J J J U − T S − μ N = − p V U-TS-\mu N = -pV U − TS − μ N = − p V d J = − S d T − p d V − N d μ dJ = -S\,dT - p\,dV - N\,d\mu dJ = − S d T − p d V − N d μ T , V , μ T,V,\mu T , V , μ 极小 巨正则(第8章 )
自测题
证明 ( ∂ 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 ,并由此说明理想气体的 C V C_V C V 与体积无关。[提示:C V = T ( ∂ S / ∂ T ) V C_V = T(\partial S/\partial T)_V C V = T ( ∂ S / ∂ T ) V ,交换求导次序后用麦克斯韦关系。]
范德瓦尔斯气体满足 ( p + a N 2 / V 2 ) ( V − N b ) = N k B T (p + aN^2/V^2)(V - Nb) = Nk_{\mathrm B}T ( p + a N 2 / V 2 ) ( V − N b ) = N k B T 。用 (1.21) 求 ( ∂ U / ∂ V ) T (\partial U/\partial V)_T ( ∂ U / ∂ V ) T 。[答:a N 2 / V 2 aN^2/V^2 a N 2 / V 2 ]
从 G = μ N G = \mu N G = μ N 与 (1.13) 出发,重新推导吉布斯–杜安关系 (1.25) 。
这一篇已记为读完。标为未读