搜索
File 2026.10.29物理

热力学与统计物理期中复习

Date / 日期
Updated / 更新
2026.10.29
Reading / 阅读
193 Min
Words / 字数
38429
Category / 分类
物理热力学

热力学与统计物理期中复习

· 更新于 2026年10月29日· 约 193 分钟· 38429 字

写在前面

这场期中考的范围是前两章,第一章热力学配第二章经典统计力学。第一章的 §1.10 有限时间热力学和第二章末尾关于负温度的讨论属于边角,主战场在 §1.2 到 §1.8 与 §2.2 到 §2.7。

这份材料的排法和课堂笔记不一样。课堂笔记顺着讲义节次走,用来第一遍学;这份按考试当天的动作排,先看全景把两章接起来,再用公式速查定位,然后靠必会推导和典型题型进入能落笔的状态,最后用易错点和自测清单收口。所以它不重述定义,只写考场上要调用的东西。

距考试还有六天,时间紧就按这个顺序走一遍。先扫一遍公式速查里的状态方程、熵方程、四个势、卡诺与微正则这几张表,把记号认熟;再啃必会推导里的卡诺效率与 S=klog⁡ΓS=k\log\Gamma 两条,这两条几乎每次都从里面出题;然后做典型题型的算熵变和循环效率,这两类是基础分;精力够的话把 Maxwell 关系套用和 van der Waals 那部分补上,那是拉分的地方。易错点放在考前最后一天看,那时候错题已经定型,看清单比看推导有用。

记号沿用前两章。自由能用 AA,玻尔兹曼常量写 kk,β=1/kT\beta=1/kT,Γ\Gamma 是 Γ\Gamma 空间里的壳层体积,Σ\Sigma 是 H<EH<E 那整块体积,ω=∂Σ/∂E\omega=\partial\Sigma/\partial E 是态密度。

全景

热力学那条线

起点是热平衡。两个系统热接触足够久之后变得一样冷热,这叫热平衡;第零定律说与第三个系统各自热平衡的两个系统彼此热平衡。这条定律的作用是把温度 TT 变成一个有良好定义的东西,否则不同的系统对会给出不同的温度,TT 就没法比较。有了 TT,再把它塞进完备集 (T,x1,…,xM)(T,x_1,\dots,x_M) 里,就得到状态方程 f(T,x1,…,xM)=0f(T,x_1,\dots,x_M)=0,理想气体的 PV=NkTPV=NkT 是最简单的一例。温度和状态方程这两样东西一到位,后面的推导才有落脚点。

接着是第一定律。功 dW=PdV\dbar W=P\dd V 依赖路径,所以写成 d\dbar 不写成 d\dd;热也一样。第一定律的真正内容是 ΔU=ΔQ−ΔW\Delta U=\Delta Q-\Delta W 对给定初末态的所有变换都取同一个值,于是 dU\dd U 是恰当微分,UU 是态函数。这一步把 dQ\dbar Q 换成了 dU+PdV\dd U+P\dd V,热就从”某种流进来的东西”变成了可以用状态量算的量。

第二定律管方向。Kelvin 表述堵住”从单一热库取热全部变功”,Clausius 表述堵住”热自发从冷到热”,两种表述可以互推,所以只取其一当公理。由它推出卡诺定理,两个给定温度之间任何热机的效率都不超过卡诺热机,而卡诺热机的效率是 1−T1/T21-T_1/T_2。这条定理顺手给出一个不依赖工作物质的绝对温标,把 TT 的定义从”选一支温度计”升级成”由效率比固定”。

再往下是熵。把任意循环用等温线和绝热线切成无穷多个小卡诺循环,每个小循环满足 Qj/Tj+Qj+1/Tj+1=0Q_j/T_j+Q_{j+1}/T_{j+1}=0,相邻小循环的内部分界线互相抵消,剩下的就是 ∮dQ/T≤0\oint\dbar Q/T\le0,可逆取等号。等号成立时 ∫ABdQ/T\int_A^B\dbar Q/T 与路径无关,于是 SS 存在,且 TdS=dU+PdVT\dd S=\dd U+P\dd V。这条式子把第一和第二定律装进一个微分式,剩下的全是它的改写。

改写走两条路。一条是 Legendre 变换,得到 A=U−TSA=U-TS、H=U+PVH=U+PV、G=H−TSG=H-TS,四个势的全微分各给一条 Maxwell 关系,并且给出等温等容下 AA 最小、等温等压下 GG 最小的平衡判据。另一条是相变,把 GG 最小用到两相共存上,得 g1=g2g_1=g_2,沿相界线微分一次就是 Clapeyron 方程 dP/dT=Δs/Δv\dd P/\dd T=\Delta s/\Delta v。van der Waals 方程从另一头补上这块,它用 V−bV-b 扣掉分子的硬球体积、用 P+a/V2P+a/V^2 扣掉器壁附近的吸引,等温线在低温下出现负压缩率的一段,Maxwell 等面积构造把那段换成一水平平台,平台的端点是两相共存。

最后是第三定律。熵的定义要接一个参考态,参考态选得不唯一,SS 就定到差一个常数;第三定律说绝对零度时熵是普适常数,可以取零,于是 SS 唯一了。代价是热容必须在 T→0T\to0 时归零,因为 S=∫0TC dT/TS=\int_0^T C\,\dd T/T 要收敛。

统计力学那条线

统计力学从另一个方向进来。等概率先验假设说孤立系统在平衡时,能量在 EE 到 E+ΔE+\Delta 之间的每个微观状态出现的概率相同,分布 ρ\rho 在这层壳上是常数。由此得到微正则系综,壳层体积记为 Γ(E)\Gamma(E)。Γ\Gamma 量的是 Γ\Gamma 空间里的体积,要数状态得先规定单位体积,这一点讲义特意提了一句。

熵的定义是 S=klog⁡Γ(E)S=k\log\Gamma(E)。它有两个必须验证的地方。广延性靠”把系统切成两个子系统后,能量分解 E=E1+E2E=E_1+E_2 的求和中只有一项占主导”来保证,一项占主导才有 log⁡Γ1Γ2=log⁡Γ1+log⁡Γ2\log\Gamma_1\Gamma_2=\log\Gamma_1+\log\Gamma_2;第二定律靠 Σ(E)\Sigma(E) 随 VV 不减来保证,因为孤立系统的 VV 只能增大。

温度压强从偏导里长出来。定义 (∂S/∂E)V=1/T(\partial S/\partial E)_V=1/T,再用 P=T(∂S/∂V)EP=T(\partial S/\partial V)_E,代进 dS\dd S 的两个偏导就回到 dE=TdS−PdV\dd E=T\dd S-P\dd V。这一步的意义是热力学第一定律可以从这里推出来。讲义特意提醒,第一定律本身在统计力学里仍然是一条假设。两个子系统能量交换到平衡时 (∂S1/∂E1)=(∂S2/∂E2)(\partial S_1/\partial E_1)=(\partial S_2/\partial E_2),正是第零定律的 T1=T2T_1=T_2,所以温度的定义和热力学是自洽的。

接着是能量均分。对任意 Hamiltonian,⟨xi∂H/∂xj⟩=δijkT\langle x_i\partial H/\partial x_j\rangle=\delta_{ij}kT,用的是把壳层积分写成差分、再分部积分、边界项为零这一套。若 HH 是二次型,每条二次型自由度分到 kT/2kT/2,于是 ⟨H⟩=fkT/2\langle H\rangle=fkT/2,单原子理想气体的 CV=3Nk/2C_V=3Nk/2 就是这么来的。

理想气体把整条链跑通一遍。Σ(E)=C3N[V(2mE)3/2]N\Sigma(E)=C_{3N}[V(2mE)^{3/2}]^N,S=klog⁡ΣS=k\log\Sigma 算出后反解出 E=E(S,V)E=E(S,V),微分得 T=2U/(3Nk)T=2U/(3Nk) 即 U=3NkT/2U=3NkT/2,以及 P=NkT/VP=NkT/V。这两个结果和第一章的实验结论对上,说明 S=klog⁡ΓS=k\log\Gamma 这个定义没选错。混合同种气体时熵会多出一项 Nklog⁡VNk\log V,破坏广延性,加 Gibbs 因子 1/N!1/N! 修掉。

最后落到具体分布。准独立粒子的情形下,把 μ\mu 空间切成小格,用占据数 {nl}\{n_l\} 描述,Γ\Gamma 空间里对应 {nl}\{n_l\} 的体积是 N!/∏nl!⋅∏glnlN!/\prod n_l!\cdot\prod g_l^{n_l}。这个量在给定 NN 与 EE 下取极大,用两个 Lagrange 乘子 α,β\alpha,\beta 得 nˉl=gle−α−βϵl\bar n_l=g_le^{-\alpha-\beta\epsilon_l},即 Maxwell-Boltzmann 分布。把求和写成单粒子配分函数 Z=∑le−βϵlZ=\sum_le^{-\beta\epsilon_l},热力学量变成对 β\beta 和 VV 的求导,E=−N∂log⁡Z/∂βE=-N\partial\log Z/\partial\beta,P=(N/β)∂log⁡Z/∂VP=(N/\beta)\partial\log Z/\partial V,温度由 β=1/kT\beta=1/kT 定。

两章怎么接起来

统计力学在这里干的事,是给热力学里的 TT、SS、PP 配一套微观算法。形式上两边的接口只有一条式子,(∂S/∂E)V=1/T(\partial S/\partial E)_V=1/T 和 P=T(∂S/∂V)EP=T(\partial S/\partial V)_E 合起来就是 dE=TdS−PdV\dd E=T\dd S-P\dd V,和第一章的热力学基本方程逐字相同,只是 SS 换成了 klog⁡Γk\log\Gamma。讲义说 (∂S/∂E)=1/T(\partial S/\partial E)=1/T 本身就是第一章的 Maxwell 关系之一,这话的意思是两边不是两套理论,是同一套关系的两种算法。

考试里这个接口常见的用法有两种。一种是给你一个微观模型,让你算 Γ\Gamma 或 ZZ,再问 UU、CVC_V、PP、SS,这时候走 S=klog⁡ΓS=k\log\Gamma 或者走 ZZ 都行,走 ZZ 快一些。另一种是给你一个宏观关系,问它在统计上对应什么,比如问能量均分为什么给出 U=3NkT/2U=3NkT/2,那就要把 Σ∝E3N/2\Sigma\propto E^{3N/2} 这一步补上。

公式速查

温度、状态方程与温标

温度这一块考得不多,但定义链条容易被问。第零定律不只是一句常识,它的作用是保证”与 CC 热平衡”这个关系可以传递,从而 TT 唯一。温标的做法分三步,选一个系统加一个随温度单调变化的性质,假定这个性质与 TT 的关系(通常取线性),再定原点和单位。第三步依赖实验,所以水银温度计这种方案不普适,只在有限温区好用。

理想气体温标好在它不需要指定具体气体。固定 VV 时设 TV=100(P−Pi)/(Pf−Pi)T_V=100(P-P_i)/(P_f-P_i),固定 PP 时设 TP=100(V−Vi)/(Vf−Vi)T_P=100(V-V_i)/(V_f-V_i),由 Boyle 定律 PfVi=VfPiP_fV_i=V_fP_i 可以证明 TV=TPT_V=T_P,证明的办法是把 PV=aVTV+bVPV=a_VT_V+b_V 和 PV=aPTP+bPPV=a_PT_P+b_P 在两条等容、等压线上比较系数,得到 aV=aPa_V=a_P、bV=bPb_V=b_P。

关系式子适用条件
第零定律A∼CA\sim C 且 B∼CB\sim C 则 A∼BA\sim B热平衡的传递性,定义 TT 的前提
状态方程f(T,x1,…,xM)=0f(T,x_1,\dots,x_M)=0平衡态,{xi}\{x_i\} 是完备集,个数由系统固定
理想气体温标PV=NkTPV=NkT稀薄气体;TT 用理想气体温标
一般温标PV=F(T)PV=F(T)FF 单值单调,任何一种选择都构成温标
绝对温标T1/T2=1−ηT_1/T_2=1-\eta卡诺热机,不依赖工作物质
绝对温标单位Tn+1−Tn=1 KT_{n+1}-T_n=1\ \mathrm{K}固定每一级卡诺循环的功
强度与广延V→λVV\to\lambda V 时广延量乘 λ\lambda,强度量不变热力学极限,可分系统才可加

第一定律与熵

这张表是整场考试用得最多的一张。要留意两件事,一是 d\dbar 与 d\dd 的区别,功和热不是态函数,只有 UU、SS、HH、AA、GG 是;二是每个式子背后的独立变量,熵方程有两种写法,分别对应 (T,V)(T,V) 和 (T,P)(T,P) 两套独立变量,混用必错。

关系式子适用条件
第一定律ΔU=ΔQ−ΔW\Delta U=\Delta Q-\Delta W任何变换,不要求准静态
无穷小形式dU=dQ−dW\dd U=\dbar Q-\dbar W任何变换,d\dbar 表示非恰当
准静态功dW=PdV\dbar W=P\dd V准静态,且只有 PdVP\dd V 这一种功
热容CV=(∂U/∂T)VC_V=(\partial U/\partial T)_V,CP=(∂H/∂T)PC_P=(\partial H/\partial T)_P任何系统,下标是路径
理想气体热容差CP−CV=NkC_P-C_V=Nk理想气体
热力学基本方程TdS=dU+PdVT\dd S=\dd U+P\dd V可逆变换,或连接两个平衡态
熵方程TdS=CVdT+T(∂P/∂T)VdVT\dd S=C_V\dd T+T(\partial P/\partial T)_V\dd V可逆,独立变量取 (T,V)(T,V)
熵方程TdS=CPdT−T(∂V/∂T)PdPT\dd S=C_P\dd T-T(\partial V/\partial T)_P\dd P可逆,独立变量取 (T,P)(T,P)
内能对体积的偏导(∂U/∂V)T=T(∂P/∂T)V−P(\partial U/\partial V)_T=T(\partial P/\partial T)_V-P由 dS\dd S 恰当性得出,任何物质
可逆熵变S(B)−S(A)=∫ABdQ/TS(B)-S(A)=\int_A^B\dbar Q/T沿可逆路径积分
不可逆熵变S(B)−S(A)≥∫ABdQ/TS(B)-S(A)\ge\int_A^B\dbar Q/T任意变换,等号只在可逆时
Clausius 定理∮dQ/T≤0\oint\dbar Q/T\le0循环,可逆取等号
绝热系统ΔS≥0\Delta S\ge0绝热,可逆取等号
等温膨胀ΔS=Nkln⁡(V2/V1)\Delta S=Nk\ln(V_2/V_1)理想气体,TT 不变
任意两态ΔS=CVln⁡(T2/T1)+Nkln⁡(V2/V1)\Delta S=C_V\ln(T_2/T_1)+Nk\ln(V_2/V_1)理想气体,CVC_V 为常数
偏导恒等式(∂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)w(∂y/∂z)w=(∂x/∂z)w(\partial x/\partial y)_w(\partial y/\partial z)_w=(\partial x/\partial z)_w固定同一个 ww

热力学势与 Maxwell 关系

四个势的记忆法则是自变量成对换。UU 配 (S,V)(S,V),AA 配 (T,V)(T,V),HH 配 (S,P)(S,P),GG 配 (T,P)(T,P),负号出现在换掉的那个变量上。把势写成它自己那对变量的函数(特性函数)之后,其余所有量都能靠微分得到,这就是 Masius 定理。

关系式子适用条件
自由能A=U−TSA=U-TS,dA=−SdT−PdV\dd A=-S\dd T-P\dd VA=A(T,V)A=A(T,V),可逆
焓H=U+PVH=U+PV,dH=TdS+VdP\dd H=T\dd S+V\dd PH=H(S,P)H=H(S,P),可逆
Gibbs 势G=H−TSG=H-TS,dG=−SdT+VdP\dd G=-S\dd T+V\dd PG=G(T,P)G=G(T,P),可逆
Maxwell 关系(∂T/∂V)S=−(∂P/∂S)V(\partial T/\partial V)_S=-(\partial P/\partial S)_V由 dU\dd U 的交叉偏导相等
Maxwell 关系(∂S/∂V)T=(∂P/∂T)V(\partial S/\partial V)_T=(\partial P/\partial T)_V由 dA\dd A,用过最多的一条
Maxwell 关系(∂T/∂P)S=(∂V/∂S)P(\partial T/\partial P)_S=(\partial V/\partial S)_P由 dH\dd H
Maxwell 关系−(∂S/∂P)T=(∂V/∂T)P-(\partial S/\partial P)_T=(\partial V/\partial T)_P由 dG\dd G,负号在左边
特性函数S=−(∂A/∂T)VS=-(\partial A/\partial T)_V,P=−(∂A/∂V)TP=-(\partial A/\partial V)_TA=A(T,V)A=A(T,V) 已知
由 AA 求 CVC_VU=A−T(∂A/∂T)VU=A-T(\partial A/\partial T)_V,CV=−T(∂2A/∂T2)VC_V=-T(\partial^2A/\partial T^2)_VA=A(T,V)A=A(T,V) 已知
等温等容平衡AA 不增,平衡时最小dW=0\dbar W=0 且 TT 固定
等温等压平衡GG 不增,平衡时最小T,PT,P 固定
两子系统平衡(∂A1/∂V1)T=(∂A2/∂V2)T(\partial A_1/\partial V_1)_T=(\partial A_2/\partial V_2)_T 即 P1=P2P_1=P_2总体积固定,活塞可动

卡诺循环

关系式子适用条件
效率η=W/Q2=1−Q1/Q2\eta=W/Q_2=1-Q_1/Q_2任何热机,Q2Q_2 吸热 Q1Q_1 放热
卡诺效率η=1−T1/T2\eta=1-T_1/T_2可逆循环,两个热库,TT 为绝对温标
卡诺定理η≤1−T1/T2\eta\le1-T_1/T_2任何热机
可逆循环的热比Q2/T2=Q1/T1Q_2/T_2=Q_1/T_1可逆,两个热库
绝热方程VγP=constV^\gamma P=\mathrm{const},Vγ−1T=constV^{\gamma-1}T=\mathrm{const}理想气体,准静态绝热,CVC_V 常数
绝热指数γ=1+Nk/CV\gamma=1+Nk/C_V理想气体,CVC_V 常数
等温吸热Q2=NkT2ln⁡(Vb/Va)Q_2=NkT_2\ln(V_b/V_a)理想气体,等温可逆
卡诺推论同两热库间所有卡诺热机效率相同效率与 WW、与工作物质都无关

换到 TT-SS 平面上看同一个循环,图形会变成一个矩形,这件事本身就值一张图。两个绝热段是等熵的,所以在图上是两条竖线;两个等温段温度不变,是两条横线。

TikZ 源码
\begin{tikzpicture}[>=stealth, thick, font=\small, x=1.55cm, y=1.15cm]
\draw[->] (0,0) -- (5.1,0) node[below] {$S$};
\draw[->] (0,0) -- (0,5.4) node[left] {$T$};
\node[below left] at (0,0) {$O$};
\draw[dashed] (1.2,0) -- (1.2,4.6);
\draw[dashed] (3.8,0) -- (3.8,4.6);
\draw[dashed] (0,2.0) -- (3.8,2.0);
\draw[dashed] (0,4.6) -- (3.8,4.6);
\node[left, font=\scriptsize] at (0,2.0) {$T_1$};
\node[left, font=\scriptsize] at (0,4.6) {$T_2$};
\node[below, font=\scriptsize] at (1.2,0) {$S_1$};
\node[below, font=\scriptsize] at (3.8,0) {$S_2$};
\fill[gray!25] (1.2,2.0) rectangle (3.8,4.6);
\draw[thick] (1.2,2.0) -- (3.8,2.0) -- (3.8,4.6) -- (1.2,4.6) -- cycle;
\draw[->, thick] (2.1,4.6) -- (2.9,4.6);
\draw[->, thick] (3.8,3.6) -- (3.8,3.0);
\draw[->, thick] (2.9,2.0) -- (2.1,2.0);
\draw[->, thick] (1.2,3.0) -- (1.2,3.6);
\node[above left, font=\scriptsize] at (1.2,4.6) {$a$};
\node[above right, font=\scriptsize] at (3.8,4.6) {$b$};
\node[below right, font=\scriptsize] at (3.8,2.0) {$c$};
\node[below left, font=\scriptsize] at (1.2,2.0) {$d$};
\node[font=\scriptsize] at (2.5,4.95) {$Q_2=T_2(S_2-S_1)$};
\node[font=\scriptsize] at (2.5,1.5) {$Q_1=T_1(S_2-S_1)$};
\node[font=\scriptsize] at (2.5,3.3) {$W=(T_2-T_1)(S_2-S_1)$};
\end{tikzpicture}

热量在 TT-SS 图上是线下的面积,因为 Q=∫TdSQ=\int T\dd S。顶边下方面积是吸热 Q2Q_2,底边下方面积是放热 Q1Q_1,两者相减就是矩形自己那块灰面积,也就是净功。于是效率一行就出来,η=1−Q1/Q2=1−T1/T2\eta=1-Q_1/Q_2=1-T_1/T_2,连绝热方程都不用碰。这也是 TT-SS 图在复习时比 pp-VV 图顺手的原因,熵是态函数,绝热线在图上必须是竖线,循环的面积关系一目了然。

相变与 van der Waals

相变这一块,两相共存的条件是 Gibbs 势最小,而 G=m1g1+m2g2G=m_1g_1+m_2g_2,对质量分配求变分得 g1=g2g_1=g_2。相变在等温等压下进行,对应的判据是 GG 最小,等温等容才轮到 AA。gg 是单位质量的 Gibbs 势,也叫化学势。

关系式子适用条件
相平衡g1(P,T)=g2(P,T)g_1(P,T)=g_2(P,T)一阶相变,等温等压
单位质量微分dg=−sdT+vdP\dd g=-s\dd T+v\dd P单相内部
Clapeyron 方程dPdT=ΔsΔv=lTΔv\dfrac{\dd P}{\dd T}=\dfrac{\Delta s}{\Delta v}=\dfrac{l}{T\Delta v}一阶相变,沿相界线
潜热l=TΔsl=T\Delta s一阶相变
一级判据Δv≠0\Delta v\ne0 或 Δs≠0\Delta s\ne0,gg 的一阶导不连续一级相变
二级相变Δs=Δv=0\Delta s=\Delta v=0,Clapeyron 失效临界点及以上
三相点三条 P(T)P(T) 相界线的交点三相共存,P,TP,T 都定死
临界点相界线的终点 (Tc,Pc)(T_c,P_c)两相差别消失,一级转二级
van der Waals 方程(P+a/V2)(V−b)=NkT(P+a/V^2)(V-b)=NkT实际气体的简单模型
vdW 临界量NkTc=8a/27bNkT_c=8a/27b,Pc=a/27b2P_c=a/27b^2,Vc=3bV_c=3b由 (V−Vc)3=0(V-V_c)^3=0 三重根定出
对应态方程(Vˉ−1/3)(Pˉ+3/Vˉ2)=8Tˉ/3(\bar V-1/3)(\bar P+3/\bar V^2)=8\bar T/3Tˉ=T/Tc\bar T=T/T_c 等,不含物质常数
内压(∂U/∂V)T=a/V2(\partial U/\partial V)_T=a/V^2vdW 气体,推导见后面「Maxwell 关系套用」

van der Waals 的两处修正各有物理来源。硬球排斥让每个分子可用的体积从 VV 减到 V−bV-b,bb 是物质常数;吸引让器壁附近的分子被往回拉,器壁感受到的压强比动理学压强小,减少的量正比于器壁附近分子层与内部的对数,约 N2/V2N^2/V^2,写成 a/V2a/V^2。把修正后的假设 PkineticVeff=NkTP_{\mathrm{kinetic}}V_{\mathrm{eff}}=NkT 写出来就是上面那条方程。

TikZ 源码
\begin{tikzpicture}[>=stealth, thick, font=\small, x=1.0cm, y=1.05cm]
\fill[gray!25] (1.8,2.35) .. controls (2.0,1.75) and (2.15,1.55) .. (2.35,1.55)
.. controls (2.6,1.6) and (2.75,2.0) .. (2.95,2.35) -- cycle;
\fill[gray!25] (2.95,2.35) .. controls (3.2,2.7) and (3.6,2.85) .. (3.9,2.85)
.. controls (4.3,2.75) and (4.8,2.5) .. (5.2,2.35) -- cycle;
\draw[->] (0,0) -- (7.9,0) node[below] {$V$};
\draw[->] (0,0) -- (0,7.7) node[left] {$P$};
\node[below left] at (0,0) {$O$};
\draw[gray, dashed] (0.84,7.55)
.. controls (1.26,5.1) and (1.56,3.6) .. (1.95,3.30)
.. controls (2.6,3.15) and (3.2,3.05) .. (3.95,3.00)
.. controls (4.7,2.92) and (5.5,2.5) .. (6.2,2.10)
.. controls (6.8,1.8) and (7.05,1.62) .. (7.35,1.50);
\node[gray, font=\scriptsize, anchor=west] at (7.5,1.45) {$T>T_c$};
\draw[blue] (0.9,7.2) .. controls (1.25,4.2) and (1.6,2.7) .. (1.8,2.35)
.. controls (2.0,1.75) and (2.15,1.55) .. (2.35,1.55)
.. controls (2.6,1.6) and (2.75,2.0) .. (2.95,2.35)
.. controls (3.2,2.7) and (3.6,2.85) .. (3.9,2.85)
.. controls (4.3,2.75) and (4.8,2.5) .. (5.2,2.35)
.. controls (5.8,2.0) and (6.5,1.4) .. (7.2,0.95);
\draw[orange, line width=1.1pt] (2.35,1.55)
.. controls (2.6,1.6) and (2.75,2.0) .. (2.95,2.35)
.. controls (3.2,2.7) and (3.6,2.85) .. (3.9,2.85);
\draw[red] (1.8,2.35) -- (5.2,2.35);
\fill (1.8,2.35) circle (2pt);
\fill (5.2,2.35) circle (2pt);
\node[below left, font=\scriptsize] at (1.76,2.3) {$1$};
\node[above right, font=\scriptsize] at (5.24,2.4) {$2$};
\node[font=\scriptsize] at (4.5,2.48) {$A$};
\node[font=\scriptsize] at (2.45,1.88) {$B$};
\node[font=\scriptsize, anchor=north] at (3.65,1.2) {橙段 $\partial P/\partial V>0$,不物理};
\end{tikzpicture}

图中蓝线是 T<TcT<T_c 的一条 van der Waals 等温线,在点 1 与点 2 之间它有一段压强随体积上升,压缩率为负,物理上不允许。红线的画法是让两块阴影面积相等,即 Maxwell 等面积构造,用它代替蓝色那一段。红线左端是液相、右端是气相,中间是两相共存,PP 与 TT 都保持不变。蓝色曲线上,从点 1 走到极小值这一段的压强低于共存压强,是亚稳的过热液体;从极大值走到点 2 这一段的压强高于共存压强,是亚稳的过冷气体;极小到极大之间那一段就是图里的橙色,压缩率为负,属于不稳定区,既不算液体也不算气体。两个亚稳段的自由能都比共存态高,一有凝结核或者气泡就会跳到红线上。灰虚线是 T>TcT>T_c 的等温线,单调下降,没有平台,两相也就分不出来了。

相图里的三条线各有来历。固液线、液气线、固气线在一点交汇,那点是三相点,三个相同时共存,PP 和 TT 都被定死,所以三相点可以拿来做温标的固定点。液气线的终点是临界点,越过它之后液体和气体的差别消失,转变从一级变成二级。

TikZ 源码
\begin{tikzpicture}[>=stealth, thick, font=\small, x=1.15cm, y=1.0cm]
\draw[->] (0,0) -- (7.4,0) node[below] {$T$};
\draw[->] (0,0) -- (0,7.2) node[left] {$P$};
\node[below left] at (0,0) {$O$};
\coordinate (TP) at (2.3,2.3);
\coordinate (CP) at (5.7,5.0);
\draw[thick] (TP) -- (2.0,6.8);
\draw[thick] (TP) .. controls (3.4,3.4) and (4.8,4.6) .. (CP);
\draw[thick] (TP) .. controls (1.85,1.4) and (1.5,0.9) .. (1.1,0.3);
\fill (TP) circle (2.2pt);
\fill (CP) circle (2.2pt);
\node[below left, font=\scriptsize] at (TP) {三相点};
\node[above right, font=\scriptsize] at (CP) {临界点 $(T_c,P_c)$};
\node at (1.25,4.9) {固};
\node at (3.6,4.3) {液};
\node at (6.5,1.9) {气};
\node[font=\scriptsize] at (1.15,6.2) {固液线};
\node[font=\scriptsize] at (4.9,5.35) {液气线};
\end{tikzpicture}

三条线的斜率都由 Clapeyron 方程决定,符号看 Δv\Delta v。从三相点向左下延伸的那条是固气线,也叫升华线。固液线的 Δv\Delta v 很小,所以它几乎竖直,也就意味着改变压强对熔点的影响远小于对沸点的影响。

微正则系综

Γ\Gamma 与 Σ\Sigma 的区别在考卷上很爱出。Σ(E)\Sigma(E) 是 H<EH<E 那整块体积,Γ(E)\Gamma(E) 是 EE 到 E+ΔE+\Delta 那层壳,ω=∂Σ/∂E\omega=\partial\Sigma/\partial E 是单位能量的体积。三者在 S=klog⁡S=k\log 里可以互换,代价只是差一个 O(log⁡N)O(\log N) 的加性项,所以算熵时挑好算的那个。

TikZ 源码
\begin{tikzpicture}[>=stealth, thick, font=\small, x=1.0cm, y=1.0cm]
\fill[gray!30] (0,0) circle (2.05);
\fill[white] (0,0) circle (1.75);
\draw[thick] (0,0) circle (1.75);
\draw[thick] (0,0) circle (2.05);
\draw[->] (-3.3,0) -- (3.3,0) node[below] {$q$};
\draw[->] (0,-3.3) -- (0,3.3) node[left] {$p$};
\node[font=\scriptsize, below left] at (0,0) {$O$};
\draw[->, thin] (0,0) -- (50:1.75);
\node[font=\scriptsize, above left] at (0.95,1.25) {$\sqrt{E}$};
\draw[->, thin] (0,0) -- (50:2.05);
\node[font=\scriptsize, anchor=north west] at (1.25,1.65) {$\sqrt{E+\Delta}$};
\node[font=\scriptsize] at (-0.55,-1.05) {$\Sigma(E)$};
\draw[->, thin] (2.1,-1.15) -- (1.62,-1.1);
\node[font=\scriptsize, anchor=west] at (2.18,-1.15) {$\Gamma(E)=\omega(E)\Delta$};
\end{tikzpicture}

以 H=p2+q2H=p^2+q^2 为例,H<EH<E 是一张圆盘,半径 E\sqrt{E};壳层就是圆环,宽度对应 Δ\Delta。Γ\Gamma 是圆环的面积,Σ\Sigma 是整张圆盘的面积,ω\omega 是周长。这个例子里圆环面积正比于 Δ\Delta、与 EE 无关,圆盘面积正比于 EE,取对数只差一个常数;推广到 NN 个自由度,差别是 O(log⁡N)O(\log N) 量级,不影响热力学量,这就是三者可以互换的来源。

关系式子适用条件
壳层体积Γ(E)=∫E<H<E+Δdp dq\Gamma(E)=\int_{E<H<E+\Delta}\dd p\,\dd q微正则,Δ≪E\Delta\ll E 但 Δ\Delta 是宏观量
累计体积Σ(E)=∫H<Edp dq\Sigma(E)=\int_{H<E}\dd p\,\dd q微正则
态密度ω(E)=∂Σ(E)/∂E\omega(E)=\partial\Sigma(E)/\partial E,Γ=ωΔ\Gamma=\omega\Delta微正则
熵S(E,V)=klog⁡Γ(E)S(E,V)=k\log\Gamma(E)微正则,Γ\Gamma 是体积需先定单位体积
三式等价S=klog⁡Γ≃klog⁡ω≃klog⁡ΣS=k\log\Gamma\simeq k\log\omega\simeq k\log\Sigma差 O(log⁡N)O(\log N) 的加性常数
温度(∂S/∂E)V=1/T(\partial S/\partial E)_V=1/T微正则
压强P=T(∂S/∂V)EP=T(\partial S/\partial V)_E微正则
等价形式P=−(∂E/∂V)SP=-(\partial E/\partial V)_S由 dE=TdS−PdV\dd E=T\dd S-P\dd V 直接读
基本方程dE=TdS−PdV\dd E=T\dd S-P\dd V只有两个独立参数
广延性S=S1(Eˉ1)+S2(Eˉ2)S=S_1(\bar E_1)+S_2(\bar E_2),Eˉ1+Eˉ2=E\bar E_1+\bar E_2=E求和只有一项占主导
平衡条件(∂S1/∂E1)=(∂S2/∂E2)(\partial S_1/\partial E_1)=(\partial S_2/\partial E_2)即 T1=T2T_1=T_2,第零定律

能量均分与理想气体

关系式子适用条件
均分定理⟨xi∂H∂xj⟩=δijkT\left\langle x_i\dfrac{\partial H}{\partial x_j}\right\rangle=\delta_{ij}kTxix_i 是 pp 或 qq,边界项可忽略
二次型结果⟨H⟩=fkT/2\langle H\rangle=fkT/2H=∑AiPi2+∑BiQi2H=\sum A_iP_i^2+\sum B_iQ_i^2
单原子气体CV=3Nk/2C_V=3Nk/2,γ=5/3\gamma=5/3只有平动,f=3f=3
谐振子⟨H⟩=kT\langle H\rangle=kT,CV=NkC_V=NkH=p2/2m+ω2q2/2H=p^2/2m+\omega^2q^2/2,f=2f=2
态数目Σ(E)=C3N[V(2mE)3/2]N\Sigma(E)=C_{3N}[V(2mE)^{3/2}]^N理想气体,Ωn(R)=CnRn\Omega_n(R)=C_nR^n
熵(未加修正)S=Nklog⁡[V(4πmE/3N)3/2]+3Nk/2S=Nk\log[V(4\pi mE/3N)^{3/2}]+3Nk/2无 Gibbs 因子
熵(加修正)S=Nklog⁡[(V/N)(4πmE/3N)3/2]+5Nk/2S=Nk\log[(V/N)(4\pi mE/3N)^{3/2}]+5Nk/2同种粒子需要 1/N!1/N!
内能T=(∂U/∂S)V=2U/(3Nk)T=(\partial U/\partial S)_V=2U/(3Nk),U=3NkT/2U=3NkT/2理想气体
压强P=−(∂U/∂V)S=NkT/VP=-(\partial U/\partial V)_S=NkT/V理想气体

配分函数与 MB 分布

关系式子适用条件
最可几分布nˉl=gle−α−βϵl\bar n_l=g_le^{-\alpha-\beta\epsilon_l}准独立粒子,NN 很大
归一化N=e−αZN=e^{-\alpha}Z,即 α=log⁡(Z/N)\alpha=\log(Z/N)由 ∑nˉl=N\sum\bar n_l=N
单粒子配分函数Z(β,y)=∑le−βϵl=∫dp dq e−βϵZ(\beta,y)=\sum_le^{-\beta\epsilon_l}=\int\dd p\,\dd q\,e^{-\beta\epsilon}准独立,ϵ≤E\epsilon\le E 的边界可忽略
内能E=−N∂log⁡Z∂βE=-N\dfrac{\partial\log Z}{\partial\beta}准独立
广义力Yk=Nβ∂log⁡Z∂ykY_k=\dfrac{N}{\beta}\dfrac{\partial\log Z}{\partial y_k}yky_k 为外参量
压强P=Nβ∂log⁡Z∂VP=\dfrac{N}{\beta}\dfrac{\partial\log Z}{\partial V}y=Vy=V,要把积分限随 VV 的依赖算进去
温度β=1/kT\beta=1/kT由 dQ=TdS\dbar Q=T\dd S 定出
熵S=Nk(log⁡Z−β∂log⁡Z∂β)+constS=Nk\left(\log Z-\beta\dfrac{\partial\log Z}{\partial\beta}\right)+\mathrm{const}与微正则差 Gibbs 因子
热量dQ=dE+∑kYkdyk\dbar Q=\dd E+\sum_kY_k\dd y_k第一定律,yy 固定时 dQ=∑ϵldnˉl\dbar Q=\sum\epsilon_l\dd\bar n_l
自由粒子Z=V(2πmkT)3/2Z=V(2\pi mkT)^{3/2}ϵ=p2/2m\epsilon=p^2/2m,无外势
理想气体E=3NkT/2E=3NkT/2,P=NkT/VP=NkT/V,CV=3Nk/2C_V=3Nk/2由上面那条 ZZ 直接求导

涨落

关系式子适用条件
相对涨落⟨f2⟩−⟨f⟩2⟨f⟩2≪1\dfrac{\langle f^2\rangle-\langle f\rangle^2}{\langle f\rangle^2}\ll1大系统;热力学极限下趋于零
最可几值⟨f⟩≃f0\langle f\rangle\simeq f_0上式成立;长程相互作用或连续相变附近失效
能量分配⟨(δE1)2⟩=kT2CV\langle(\delta E_1)^2\rangle=kT^2C_V两子系统交换能量,峰宽二阶展开
相对量级δE/E∼1/N\delta E/E\sim1/\sqrt{N}广延系统

必会推导

卡诺定理与卡诺效率

推导分两段,先算理想气体卡诺热机的效率,再证任何热机都超不过它。两段的结论合起来才是卡诺定理。

理想气体卡诺循环由两条等温线和两条绝热线围成,按 a→b→c→d→aa\to b\to c\to d\to a 走,T2T_2 是高温,T1T_1 是低温。

TikZ 源码
\begin{tikzpicture}[>=stealth, thick, font=\small, x=1.15cm, y=0.72cm]
\draw[->] (0,0) -- (7.5,0) node[below] {$V$};
\draw[->] (0,0) -- (0,10.3) node[left] {$P$};
\node[below left] at (0,0) {$O$};
\draw[red] plot[domain=1.2:3.4, samples=60] (\x, {10.8/\x});
\draw[blue] plot[domain=3.4:6.8, samples=60] (\x, {17.6204/pow(\x,1.4)});
\draw[red] plot[domain=2.4:6.8, samples=60] (\x, {8.1849/\x});
\draw[blue] plot[domain=1.2:2.4, samples=60] (\x, {11.6171/pow(\x,1.4)});
\draw[->, red] (2.0,5.400) -- (2.2,4.909);
\draw[->, blue] (4.0,2.530) -- (4.2,2.363);
\draw[->, red] (5.0,1.637) -- (4.7,1.741);
\draw[->, blue] (2.0,4.402) -- (1.85,4.910);
\fill (1.2,9.0) circle (1.8pt) node[above left, font=\scriptsize] {$a$};
\fill (3.4,3.176) circle (1.8pt) node[above right, font=\scriptsize] {$b$};
\fill (6.8,1.204) circle (1.8pt) node[below right, font=\scriptsize] {$c$};
\fill (2.4,3.410) circle (1.8pt) node[below left, font=\scriptsize] {$d$};
\node[red, font=\scriptsize, anchor=west] at (4.30,6.35) {等温 $T_2$,吸热 $Q_2$};
\draw[->, red, thin] (4.20,6.05) -- (2.95,3.80);
\node[red, font=\scriptsize] at (4.55,0.72) {等温 $T_1$,放热 $Q_1$};
\node[blue, font=\scriptsize] at (5.45,3.05) {绝热 $Q=0$};
\draw[->, blue, thin] (5.32,2.82) -- (5.12,1.92);
\node[font=\scriptsize, anchor=west] at (4.55,9.6) {红 等温线};
\node[font=\scriptsize, anchor=west] at (4.55,8.7) {蓝 绝热线};
\end{tikzpicture}

图中 a→ba\to b 与 c→dc\to d 是两条等温线,b→cb\to c 与 d→ad\to a 是两条绝热线,四条线围出的就是卡诺循环。同一体积处两条线的陡峭程度不一样,等温线的斜率是 −P/V-P/V,绝热线的是 −γP/V-\gamma P/V,比值正好是 γ\gamma。这个差别在图上表现为绝热线离原点更近,也是判断一段曲线属于等温还是绝热的依据。

先推绝热方程。设 PV=NkTPV=NkT、U=CVT+constU=C_VT+\mathrm{const} 且 CVC_V 为常数,绝热时 dQ=0\dbar Q=0,所以

CVdT=−PdV\begin{align} C_V\dd T=-P\dd V \end{align}

对状态方程微分得 PdV+VdP=NkdTP\dd V+V\dd P=Nk\dd T,消去 dT\dd T,

(1+NkCV)PdV+VdP=0⟹VγP=const,γ=1+NkCV\begin{align} \left(1+\frac{Nk}{C_V}\right)P\dd V+V\dd P=0 \quad\Longrightarrow\quad V^\gamma P=\mathrm{const},\qquad \gamma=1+\frac{Nk}{C_V} \end{align}

把 P=NkT/VP=NkT/V 代进去还可以写成 Vγ−1T=constV^{\gamma-1}T=\mathrm{const}。这一步只能对理想气体做,因为消 dT\dd T 用了 U=U(T)U=U(T)。

两条绝热线给出 Vbγ−1T2=Vcγ−1T1V_b^{\gamma-1}T_2=V_c^{\gamma-1}T_1 和 Vdγ−1T1=Vaγ−1T2V_d^{\gamma-1}T_1=V_a^{\gamma-1}T_2,两式相除得 Vb/Va=Vc/VdV_b/V_a=V_c/V_d,这个比值后面要用。

两个等温段的热量,等温时 ΔU=0\Delta U=0,所以 Q=∫PdVQ=\int P\dd V,

Q2=NkT2ln⁡VbVa,Q1=NkT1ln⁡VcVd\begin{align} Q_2=NkT_2\ln\frac{V_b}{V_a},\qquad Q_1=NkT_1\ln\frac{V_c}{V_d} \end{align}

代入效率定义,对数部分的比值相同,一起约掉,

η=1−Q1Q2=1−T1ln⁡(Vc/Vd)T2ln⁡(Vb/Va)=1−T1T2\begin{align} \eta=1-\frac{Q_1}{Q_2} =1-\frac{T_1\ln(V_c/V_d)}{T_2\ln(V_b/V_a)} =1-\frac{T_1}{T_2} \end{align}

第二段证一般热机超不过卡诺热机。取一台卡诺热机 CC 和一台任意热机 XX,都工作在 T2>T1T_2>T_1 之间,记 W=Q2−Q1W=Q_2-Q_1、W′=Q2′−Q1′W'=Q_2'-Q_1'。选正整数 N,N′N,N' 使 Q2/Q2′=N′/NQ_2/Q_2'=N'/N,整数取得足够大就能逼近任何比值。

让 CC 反向运行 NN 个循环,XX 正向运行 N′N' 个循环。反向运行的意思是 W,Q1,Q2W,Q_1,Q_2 的符号同时改变,这一步是卡诺热机可逆才做得到的。净效果是

Wtotal=N′W′−NW(Q2)total=N′Q2′−NQ2=0(Q1)total=N′Q1′−NQ1\begin{align} W_{\mathrm{total}}&=N'W'-NW\\ (Q_2)_{\mathrm{total}}&=N'Q_2'-NQ_2=0\\ (Q_1)_{\mathrm{total}}&=N'Q_1'-NQ_1 \end{align}

高温热库的净吸热为零,于是整个装置只跟低温热库交换热。一圈下来 ΔU=0\Delta U=0,由第一定律 Wtotal=−(Q1)totalW_{\mathrm{total}}=-(Q_1)_{\mathrm{total}}。如果 Wtotal>0W_{\mathrm{total}}>0,就是从单一热库取热全部变功,违反 Kelvin 表述,所以 Wtotal≤0W_{\mathrm{total}}\le0,即 (Q1)total≥0(Q_1)_{\mathrm{total}}\ge0,也就是 N′Q1′≥NQ1N'Q_1'\ge NQ_1。两边除以 N′Q2′=NQ2N'Q_2'=NQ_2,

Q1Q2≤Q1′Q2′⟹1−Q1Q2≥1−Q1′Q2′⟹η≥η′\begin{align} \frac{Q_1}{Q_2}\le\frac{Q_1'}{Q_2'} \quad\Longrightarrow\quad 1-\frac{Q_1}{Q_2}\ge1-\frac{Q_1'}{Q_2'} \quad\Longrightarrow\quad \eta\ge\eta' \end{align}

让 XX 也取卡诺热机,把两者的角色对调再证一遍,会得到反向的不等式,两条合起来取等号,于是所有卡诺热机效率相同,与工作物质和功的大小都无关。这个”与工作物质无关”正是绝对温标的立足点。定义 T1/T2≡1−ηT_1/T_2\equiv1-\eta,效率只由两个温度的比值决定;再固定每一级卡诺循环的功为 1 K1\ \mathrm{K} 对应的量,单位也定死了。这样定义的 TT 与理想气体温标 PV=NkTPV=NkT 给出的 TT 一致。

熵是态函数

要证的是 ∫ABdQ/T\int_A^B\dbar Q/T 在可逆路径上与路径无关。出发点是一条任意循环的 Clausius 不等式

∮dQT≤0\begin{align} \oint\frac{\dbar Q}{T}\le0 \end{align}

可逆取等号。做法是用等温线和绝热线把循环切成很多小卡诺循环。每个小卡诺循环满足 Qj/Tj+Qj+1/Tj+1=0Q_j/T_j+Q_{j+1}/T_{j+1}=0,把所有小循环加起来,落在内部的分界线被相邻两个小循环各走一次、方向相反,互相抵消,只剩外边界,于是 ∮revdQ/T=0\oint_{\mathrm{rev}}\dbar Q/T=0。如果循环不可逆,同一步骤给出严格小于零,因为不可逆过程是向平衡态弛豫,热失去了做功的能力。

有了 ∮revdQ/T=0\oint_{\mathrm{rev}}\dbar Q/T=0,取两条从 AA 到 BB 的可逆路径 C1C_1、C2C_2,让 C1C_1 接上 C2C_2 的反向,构成一个可逆循环,于是

∫C1dQT−∫C2dQT=0⟹∫C1dQT=∫C2dQT\begin{align} \int_{C_1}\frac{\dbar Q}{T}-\int_{C_2}\frac{\dbar Q}{T}=0 \quad\Longrightarrow\quad \int_{C_1}\frac{\dbar Q}{T}=\int_{C_2}\frac{\dbar Q}{T} \end{align}

积分与路径无关,就可以拿它定义态函数 SS,S(B)−S(A)=∫ABdQ/TS(B)-S(A)=\int_A^B\dbar Q/T,等价说法是 dS=dQ/T\dd S=\dbar Q/T 恰当,1/T1/T 是积分因子。

两个直接推论。对任意变换,构造一条从 BB 回到 AA 的可逆路径,与原来的任意路径拼成循环,Clausius 不等式给出 S(B)−S(A)≥∫ABdQ/TS(B)-S(A)\ge\int_A^B\dbar Q/T;绝热时右边为零,于是 ΔS≥0\Delta S\ge0,等号只在可逆时成立。

四个 Maxwell 关系

来源只有一句话,四个势的全微分都是恰当的,所以交叉偏导相等。以 dA=−SdT−PdV\dd A=-S\dd T-P\dd V 为例,把 −S-S 看成 TT 方向的系数、−P-P 看成 VV 方向的系数,

(∂S∂V)T=(∂P∂T)V\begin{align} \left(\frac{\partial S}{\partial V}\right)_T=\left(\frac{\partial P}{\partial T}\right)_V \end{align}

另外三个照同样的动作做,只是换势、换变量对,

(∂T∂V)S=−(∂P∂S)V由 dU=TdS−PdV(∂S∂V)T=(∂P∂T)V由 dA=−SdT−PdV(∂T∂P)S=(∂V∂S)P由 dH=TdS+VdP−(∂S∂P)T=(∂V∂T)P由 dG=−SdT+VdP\begin{align} \left(\frac{\partial T}{\partial V}\right)_S&=-\left(\frac{\partial P}{\partial S}\right)_V &&\text{由 } \dd U=T\dd S-P\dd V\\ \left(\frac{\partial S}{\partial V}\right)_T&=\left(\frac{\partial P}{\partial T}\right)_V &&\text{由 } \dd A=-S\dd T-P\dd V\\ \left(\frac{\partial T}{\partial P}\right)_S&=\left(\frac{\partial V}{\partial S}\right)_P &&\text{由 } \dd H=T\dd S+V\dd P\\ -\left(\frac{\partial S}{\partial P}\right)_T&=\left(\frac{\partial V}{\partial T}\right)_P &&\text{由 } \dd G=-S\dd T+V\dd P \end{align}

负号的来源是换变量时符号会翻。dU\dd U 里 −PdV-P\dd V 带着负号,交叉偏导就多一个负号;dG\dd G 里 −SdT-S\dd T 带着负号,负号跑到左边。所以背这四条不如每次现场推,两行就能写完,还不会记错符号。

从 S=klog⁡ΓS=k\log\Gamma 到温度与压强

第一步证 SS 广延。把系统切成两个固定 N1,V1,N2,V2N_1,V_1,N_2,V_2 的子系统,H=H1+H2H=H_1+H_2,忽略相互作用但保留交换能量的能力。整个系统的壳层体积是

Γ(E)=∫E1<H1<E1+Δ, E2<H2<E2+Δdp1dq1dp2dq2=Γ1(E1)Γ2(E2)\begin{align} \Gamma(E)=\int_{E_1<H_1<E_1+\Delta,\,E_2<H_2<E_2+\Delta}\dd p_1\dd q_1\dd p_2\dd q_2 =\Gamma_1(E_1)\Gamma_2(E_2) \end{align}

麻烦在于 E=E1+E2E=E_1+E_2 的分解不唯一,E1E_1 可以从 00 变到 EE。把 EE 按 Δ\Delta 切成 E/ΔE/\Delta 段,

Γ(E)=∑i=0E/ΔΓ1(Ei)Γ2(E−Ei)\begin{align} \Gamma(E)=\sum_{i=0}^{E/\Delta}\Gamma_1(E_i)\Gamma_2(E-E_i) \end{align}

关键是求和里只有一项占主导。设最大的那一项对应 Eˉ1+Eˉ2=E\bar E_1+\bar E_2=E,则 Γ1(Eˉ1)Γ2(Eˉ2)≤Γ(E)≤Γ1(Eˉ1)Γ2(Eˉ2)⋅E/Δ\Gamma_1(\bar E_1)\Gamma_2(\bar E_2)\le\Gamma(E)\le\Gamma_1(\bar E_1)\Gamma_2(\bar E_2)\cdot E/\Delta,取对数后上下界只差 klog⁡(E/Δ)k\log(E/\Delta),而 S∝NS\propto N,这一项是 O(log⁡N)O(\log N),可以丢掉。于是

S(E,V)=S1(Eˉ1,V1)+S2(Eˉ2,V2)+O(log⁡N)\begin{align} S(E,V)=S_1(\bar E_1,V_1)+S_2(\bar E_2,V_2)+O(\log N) \end{align}

主导项的位置由 δlog⁡[Γ1Γ2]=0\delta\log[\Gamma_1\Gamma_2]=0 配上 δE1+δE2=0\delta E_1+\delta E_2=0 定出,即

∂S1∂E1∣Eˉ1=∂S2∂E2∣Eˉ2\begin{align} \left.\frac{\partial S_1}{\partial E_1}\right|_{\bar E_1}=\left.\frac{\partial S_2}{\partial E_2}\right|_{\bar E_2} \end{align}

既然两边相等,就把它定义成 1/T1/T,

(∂S∂E)V≡1T\begin{align} \left(\frac{\partial S}{\partial E}\right)_V\equiv\frac{1}{T} \end{align}

这一步同时给出了第零定律 T1=T2T_1=T_2。第二步定压强。有两个独立参数时,

dS=(∂S∂E)VdE+(∂S∂V)EdV=1TdE+(∂S∂V)EdV\begin{align} \dd S=\left(\frac{\partial S}{\partial E}\right)_V\dd E+\left(\frac{\partial S}{\partial V}\right)_E\dd V =\frac{1}{T}\dd E+\left(\frac{\partial S}{\partial V}\right)_E\dd V \end{align}

定义 P≡T(∂S/∂V)EP\equiv T(\partial S/\partial V)_E,上式就变成 dE=TdS−PdV\dd E=T\dd S-P\dd V。要检查这个 PP 是不是真的压强,看两条。它得是强度量,SS 和 VV 都广延,偏导是强度量,乘 TT 还是强度量。它还得和第一章的压强对得上,在 SS 固定时 dE=−PdV\dd E=-P\dd V,所以 P=−(∂E/∂V)SP=-(\partial E/\partial V)_S,正是第一章那个压强。

还要说明 SS 满足第二定律。用 Σ(E)\Sigma(E) 的写法,孤立系统 NN 与 EE 不变、VV 只增不减,而 Σ(E)\Sigma(E) 是 VV 的非减函数(体积变大,H<EH<E 的区域只会变大),所以 SS 不减。

能量均分定理

要算的是 ⟨xi∂H/∂xj⟩\langle x_i\partial H/\partial x_j\rangle 在微正则系综下的值,xix_i 取 pip_i 或 qiq_i。按定义

⟨xi∂H∂xj⟩=1Γ(E)∫E<H<E+Δdp dq xi∂H∂xj\begin{align} \left\langle x_i\frac{\partial H}{\partial x_j}\right\rangle =\frac{1}{\Gamma(E)}\int_{E<H<E+\Delta}\dd p\,\dd q\,x_i\frac{\partial H}{\partial x_j} \end{align}

把壳层积分写成两个累计积分之差,再写成对 EE 的导数,

⟨xi∂H∂xj⟩=ΔΓ(E)∂∂E∫H<Edp dq xi∂H∂xj\begin{align} \left\langle x_i\frac{\partial H}{\partial x_j}\right\rangle =\frac{\Delta}{\Gamma(E)}\frac{\partial}{\partial E}\int_{H<E}\dd p\,\dd q\,x_i\frac{\partial H}{\partial x_j} \end{align}

里面的积分做分部积分。注意 ∂H/∂xj=∂(H−E)/∂xj\partial H/\partial x_j=\partial(H-E)/\partial x_j,

∫H<Edp dq xi∂(H−E)∂xj=∫H<Edp dq ∂∂xj[xi(H−E)]−δij∫H<Edp dq (H−E)\begin{align} \int_{H<E}\dd p\,\dd q\,x_i\frac{\partial(H-E)}{\partial x_j} =\int_{H<E}\dd p\,\dd q\,\frac{\partial}{\partial x_j}\big[x_i(H-E)\big] -\delta_{ij}\int_{H<E}\dd p\,\dd q\,(H-E) \end{align}

第一项是全导数,化成 H=EH=E 曲面上的面积分,而在那个曲面上 H−E=0H-E=0,所以整项为零。剩下

⟨xi∂H∂xj⟩=δijω(E)∂∂E∫H<Edp dq (E−H)=δijω(E)Σ(E)\begin{align} \left\langle x_i\frac{\partial H}{\partial x_j}\right\rangle =\frac{\delta_{ij}}{\omega(E)}\frac{\partial}{\partial E}\int_{H<E}\dd p\,\dd q\,(E-H) =\frac{\delta_{ij}}{\omega(E)}\Sigma(E) \end{align}

这里用了 ∂(E−H)/∂E=1\partial(E-H)/\partial E=1。最后把 Σ/ω\Sigma/\omega 换成温度,因为 ω=∂Σ/∂E\omega=\partial\Sigma/\partial E,

Σω=(∂log⁡Σ∂E)−1=k(∂S∂E)−1=kT\begin{align} \frac{\Sigma}{\omega}=\left(\frac{\partial\log\Sigma}{\partial E}\right)^{-1} =k\left(\frac{\partial S}{\partial E}\right)^{-1}=kT \end{align}

于是 ⟨xi∂H/∂xj⟩=δijkT\langle x_i\partial H/\partial x_j\rangle=\delta_{ij}kT。这一条对任何 Hamiltonian 都成立,只要边界项能丢掉。

要得到 ⟨H⟩=fkT/2\langle H\rangle=fkT/2 还得加一个条件,HH 必须是二次型。设 H=∑iAiPi2+∑iBiQi2H=\sum_iA_iP_i^2+\sum_iB_iQ_i^2,由齐次函数的 Euler 定理

∑i(Pi∂H∂Pi+Qi∂H∂Qi)=2H\begin{align} \sum_i\left(P_i\frac{\partial H}{\partial P_i}+Q_i\frac{\partial H}{\partial Q_i}\right)=2H \end{align}

两边取系综平均,左边每个自由度贡献 kTkT,共 ff 个,于是 fkT=2⟨H⟩fkT=2\langle H\rangle,即 ⟨H⟩=fkT/2\langle H\rangle=fkT/2。单原子理想气体 f=3f=3,U=3NkT/2U=3NkT/2;谐振子 H=p2/2m+ω2q2/2H=p^2/2m+\omega^2q^2/2 的两个变量都是二次的,f=2f=2,⟨H⟩=kT\langle H\rangle=kT。

从微正则到 MB 分布

把 μ\mu 空间切成 KK 个小格,第 ll 格的体积记为 glg_l(最后可以取 11),粒子数 nln_l,约束是 ∑nl=N\sum n_l=N 与 ∑ϵlnl=E\sum\epsilon_ln_l=E。给定 {nl}\{n_l\} 之后,Γ\Gamma 空间里对应的体积正比于把 NN 个分子分到这 KK 组里的分法数,

Ω({nl})∝N!∏lnl!∏lglnl\begin{align} \Omega(\{n_l\})\propto\frac{N!}{\prod_ln_l!}\prod_lg_l^{n_l} \end{align}

N!N! 来自分子的排法,nl!n_l! 除掉同一格内互换(同一格里分子对 f(p,q)f(p,q) 的贡献相同,互换不算新状态)。分子数很大,用 Stirling 公式 log⁡n!≃n(log⁡n−1)\log n!\simeq n(\log n-1),

log⁡Ω({nl})=N(log⁡N−1)−∑lnl(log⁡nl−1)+∑lnllog⁡gl\begin{align} \log\Omega(\{n_l\})=N(\log N-1)-\sum_ln_l(\log n_l-1)+\sum_ln_l\log g_l \end{align}

在两条约束下求极值,引入乘子 α\alpha、β\beta,

δ[log⁡Ω−α∑lnl−β∑lϵlnl]=0\begin{align} \delta\left[\log\Omega-\alpha\sum_ln_l-\beta\sum_l\epsilon_ln_l\right]=0 \end{align}

把 δnl\delta n_l 看成独立变分,每个系数都得为零,

−log⁡nˉl+log⁡gl−α−βϵl=0⟹nˉl=gle−α−βϵl\begin{align} -\log\bar n_l+\log g_l-\alpha-\beta\epsilon_l=0 \quad\Longrightarrow\quad \bar n_l=g_le^{-\alpha-\beta\epsilon_l} \end{align}

二阶变分是 −∑l(δnl)2/nˉl<0-\sum_l(\delta n_l)^2/\bar n_l<0,所以这个驻点确实是极大。再用 ∑nˉl=N\sum\bar n_l=N 和 ∑ϵlnˉl=E\sum\epsilon_l\bar n_l=E 定 α\alpha 与 β\beta。写成配分函数就是 N=e−αZN=e^{-\alpha}Z,α=log⁡(Z/N)\alpha=\log(Z/N)。

β\beta 的身份靠熵来定。把 nˉl\bar n_l 代回 log⁡Ω\log\Omega,用 N=e−αZN=e^{-\alpha}Z 化掉 α\alpha,得到

S=Nk(log⁡Z−β∂log⁡Z∂β)+const\begin{align} S=Nk\left(\log Z-\beta\frac{\partial\log Z}{\partial\beta}\right)+\mathrm{const} \end{align}

固定外参量 yy 对它微分,log⁡Z\log Z 的两项正好抵消,只剩

dS=Nk d ⁣(log⁡Z−β∂log⁡Z∂β)=kβ dE\begin{align} \dd S=Nk\,\dd\!\left(\log Z-\beta\frac{\partial\log Z}{\partial\beta}\right) =k\beta\,\dd E \end{align}

要求 dQ=TdS\dbar Q=T\dd S,也就是 dS=dE/T\dd S=\dd E/T(yy 固定时不做功),比较即得 β=1/kT\beta=1/kT。到这一步热力学第二定律在统计力学里被推出来了,前提是第一定律仍然是假设。

典型题型

算熵变

题面给初末态和一条路径,问 ΔS\Delta S,有时顺带问环境的熵变和”损失了多少可用能”。判断第一步是这个过程可不可逆。可逆就直接沿路径积 ∫dQ/T\int\dbar Q/T;不可逆就不能用路径上的 dQ\dbar Q,要另设计一条连接同样初末态的可逆路径,因为 SS 是态函数,走哪条路都行。常用的是先等温后等容,或者直接用 dS=(CV/T)dT+(∂P/∂T)VdV\dd S=(C_V/T)\dd T+(\partial P/\partial T)_V\dd V 沿一条方便的折线积分。

例,一摩尔理想气体在 300 K300\ \mathrm{K} 下自由膨胀到两倍体积,求熵变。自由膨胀 ΔW=0\Delta W=0,绝热 ΔQ=0\Delta Q=0,所以由第一定律 ΔU=0\Delta U=0,理想气体 U=U(T)U=U(T) 给出 ΔT=0\Delta T=0。实际的路径上 dQ≡0\dbar Q\equiv0,直接积会得零,那是错的。换一条可逆等温膨胀路径,温度不变,ΔU=0\Delta U=0,气体对外做的功全部来自吸热,

ΔS=∫dQT=1T∫PdV=Nkln⁡V2V1=Nkln⁡2\begin{align} \Delta S=\int\frac{\dbar Q}{T}=\frac{1}{T}\int P\dd V=Nk\ln\frac{V_2}{V_1}=Nk\ln2 \end{align}

取 Nk=R≃8.3 J/KNk=R\simeq8.3\ \mathrm{J/K},得 ΔS≃5.8 J/K\Delta S\simeq5.8\ \mathrm{J/K}。环境没变,ΔSenv=0\Delta S_{\mathrm{env}}=0,总熵增加了 5.8 J/K5.8\ \mathrm{J/K},符合不可逆过程的 ΔS>∫dQ/T=0\Delta S>\int\dbar Q/T=0。这 TΔST\Delta S 就是那句”损失了可用的能量”,同一份熵增如果在可逆过程里,还能换成 5.8×300≈1.7 kJ5.8\times300\approx1.7\ \mathrm{kJ} 的功。

循环效率

题面给 pp-VV 图上的一圈,问效率、功或某一支的吸放热。一圈下来 ΔU=0\Delta U=0,所以 W=QnetW=Q_{\mathrm{net}},而 WW 就是回路围出的面积(顺时针为正)。效率的分母只是吸热的那几段之和,不要把放热段混进去,这是最常见的失分点。若回路由等温线和绝热线组成,先认卡诺,直接用 1−T1/T21-T_1/T_2;若不是卡诺,就逐段算 QQ。

例,工作在 400 K400\ \mathrm{K} 与 300 K300\ \mathrm{K} 之间的卡诺热机,每循环从高温热库吸热 800 J800\ \mathrm{J}。效率 η=1−300/400=0.25\eta=1-300/400=0.25,做功 W=ηQ2=200 JW=\eta Q_2=200\ \mathrm{J},向低温热库放热 Q1=Q2−W=600 JQ_1=Q_2-W=600\ \mathrm{J}。核对 Q2/T2=2 J/KQ_2/T_2=2\ \mathrm{J/K}、Q1/T1=2 J/KQ_1/T_1=2\ \mathrm{J/K},比值为一,可逆循环的判据成立。同样两个热库之间的任何其他热机,效率都不超过 0.250.25;想提高效率只能抬高 T2T_2 或压低 T1T_1,这是卡诺定理的实用含义。

Maxwell 关系套用

题面要一个不好测的偏导,比如 (∂S/∂P)T(\partial S/\partial P)_T、(∂U/∂V)T(\partial U/\partial V)_T、(∂H/∂P)T(\partial H/\partial P)_T。做法是先认它出现在哪个势的全微分里,用 Maxwell 关系换成可测量((∂V/∂T)P(\partial V/\partial T)_P、(∂P/∂T)V(\partial P/\partial T)_V、CPC_P 这类),再把状态方程代进去。不要凭记忆写关系式,现场从 dA\dd A 或 dG\dd G 推,两行的事。

例,求 van der Waals 气体的 (∂U/∂V)T(\partial U/\partial V)_T。这一条正好把第一章的公式和第二章要用到的状态方程接上。由能量方程

(∂U∂V)T=T(∂P∂T)V−P\begin{align} \left(\frac{\partial U}{\partial V}\right)_T=T\left(\frac{\partial P}{\partial T}\right)_V-P \end{align}

对 (P+a/V2)(V−b)=NkT(P+a/V^2)(V-b)=NkT 求偏导,aa、bb 与 TT 无关,

(∂P∂T)V=NkV−b\begin{align} \left(\frac{\partial P}{\partial T}\right)_V=\frac{Nk}{V-b} \end{align}

代回去,

(∂U∂V)T=NkTV−b−(NkTV−b−aV2)=aV2\begin{align} \left(\frac{\partial U}{\partial V}\right)_T=\frac{NkT}{V-b}-\left(\frac{NkT}{V-b}-\frac{a}{V^2}\right)=\frac{a}{V^2} \end{align}

结果是正的,说明等温膨胀时内能增加,因为要克服分子间的吸引做功。a=0a=0 时回到理想气体的 (∂U/∂V)T=0(\partial U/\partial V)_T=0,也就回到 U=U(T)U=U(T)。这道题还顺带给出一个结论,vdW 气体的内能可以写成 U=CVT−a/V+constU=C_VT-a/V+\mathrm{const},把它对 TT 求偏导得 CVC_V,与体积无关。

相变与潜热

题面给潜热、比容差和温度,问相界线的斜率,或者反过来问沸点随压强怎么变。核心只有一条,g1=g2g_1=g_2。沿相界线两边都满足这个等式,把 dg=−sdT+vdP\dd g=-s\dd T+v\dd P 对两相各写一遍再相减,沿同一条相界线走时两边的 dT\dd T 与 dP\dd P 相同,于是 ΔsdT=ΔvdP\Delta s\dd T=\Delta v\dd P,整理即得 Clapeyron 方程。用 l=TΔsl=T\Delta s 还可以写成 l/(TΔv)l/(T\Delta v)。

例,取潜热的量级 l∼2×106 J/kgl\sim2\times10^{6}\ \mathrm{J/kg}、比容差的量级 Δv∼1 m3/kg\Delta v\sim1\ \mathrm{m^3/kg}、T=373 KT=373\ \mathrm{K},则

dPdT=lTΔv∼2×106373×1≈5×103 Pa/K\begin{align} \frac{\dd P}{\dd T}=\frac{l}{T\Delta v}\sim\frac{2\times10^{6}}{373\times1}\approx5\times10^{3}\ \mathrm{Pa/K} \end{align}

数量级是每开尔文几 kPa,即每升高 1 K1\ \mathrm{K} 饱和蒸气压上升几千帕。这个量级说明为什么高压锅里沸点能显著抬高。符号也值得记,气相的熵与比容都比液相大,Δs>0\Delta s>0、Δv>0\Delta v>0,所以 dP/dT>0\dd P/\dd T>0,沸点随压强上升。若某处 Δv\Delta v 为负(凝固时体积变大),相界线的斜率就反过来,固液线几乎竖直就是 Δv\Delta v 很小的结果。

理想气体配分函数

题面给 ϵ(p)\epsilon(p) 或 HH,要求 EE、PP、CVC_V、SS。路线固定,先算 ZZ,再按表里的公式求导。动量积分是高斯积分,∫d3p e−βp2/2m=(2πm/β)3/2\int\dd^3p\,e^{-\beta p^2/2m}=(2\pi m/\beta)^{3/2}。

例,ϵ=p2/2m\epsilon=p^2/2m,无外势。位置积分直接给出体积 VV,

Z=∫d3p d3q e−βp2/2m=V(2πmβ)3/2=V(2πmkT)3/2\begin{align} Z=\int\dd^3p\,\dd^3q\,e^{-\beta p^2/2m}=V\left(\frac{2\pi m}{\beta}\right)^{3/2}=V(2\pi mkT)^{3/2} \end{align}

取对数 log⁡Z=log⁡V+32log⁡(2πmkT)\log Z=\log V+\frac32\log(2\pi mkT),再求导,

E=−N∂log⁡Z∂β=3N2β=32NkTP=Nβ∂log⁡Z∂V=NkTV\begin{align} E&=-N\frac{\partial\log Z}{\partial\beta}=\frac{3N}{2\beta}=\frac32NkT\\ P&=\frac{N}{\beta}\frac{\partial\log Z}{\partial V}=\frac{NkT}{V} \end{align}

CV=(∂E/∂T)V=3Nk/2C_V=(\partial E/\partial T)_V=3Nk/2。这三个结果和微正则路线算出来的一致,可以互相验算。注意求 PP 时 ZZ 里的 VV 只出现在一个因子上,如果 ϵ\epsilon 本身依赖 VV(比如把 qq 的积分限放大),就得把那个依赖也算进去,讲义在这一处专门提了一句。

能量均分应用

题面给一个 HH,问 ⟨H⟩\langle H\rangle、CVC_V 或 γ\gamma。做法是数二次型自由度的个数 ff,然后 ⟨H⟩=fkT/2\langle H\rangle=fkT/2、CV=fk/2C_V=fk/2、γ=1+2/f\gamma=1+2/f。数的时候只数二次项,一次项、四次项、1/r1/r 项都不算。

例,谐振子 H=p2/2m+ω2q2/2H=p^2/2m+\omega^2q^2/2。两个变量 pp、qq 都以二次方出现,f=2f=2,所以 ⟨H⟩=2⋅kT/2=kT\langle H\rangle=2\cdot kT/2=kT,CV=NkC_V=Nk,γ=1+2/2=2\gamma=1+2/2=2。这个结果也可以直接验证,pp 的积分和 qq 的积分各是高斯型,各自贡献 kT/2kT/2。单原子理想气体只有三个平动动量,f=3f=3,CV=3Nk/2C_V=3Nk/2,γ=5/3\gamma=5/3。反过来,如果某个坐标只以 q4q^4 的形式出现,那一项不给 kT/2kT/2,这个自由度不能算。

涨落估算

题面问”能不能把两个子系统的能量当成确定值相加”,或者问某个量的相对误差有多大。要算的是峰宽。把 Γ1(E1)Γ2(E−E1)\Gamma_1(E_1)\Gamma_2(E-E_1) 的对数在极值附近展开到二阶,一阶项给 T1=T2T_1=T_2,二阶项给宽度。

log⁡[Γ1(E1)Γ2(E−E1)]=log⁡Γ1Γ2∣Eˉ1+12k(S1′′+S2′′)(δE1)2+⋯\begin{align} \log\big[\Gamma_1(E_1)\Gamma_2(E-E_1)\big] =\log\Gamma_{1}\Gamma_{2}\Big|_{\bar E_1} +\frac{1}{2k}\left(S_1''+S_2''\right)(\delta E_1)^2+\cdots \end{align}

其中 S′′=∂2S/∂E2S''=\partial^2S/\partial E^2。由 (∂S/∂E)V=1/T(\partial S/\partial E)_V=1/T 再求一次导,

∂2S∂E2=∂∂E(1T)=−1T2(∂T∂E)V=−1T2CV\begin{align} \frac{\partial^2S}{\partial E^2} =\frac{\partial}{\partial E}\left(\frac{1}{T}\right) =-\frac{1}{T^2}\left(\frac{\partial T}{\partial E}\right)_V =-\frac{1}{T^2C_V} \end{align}

所以 S1′′+S2′′=−(1/T2)(1/CV1+1/CV2)=−1/(T2CV)S_1''+S_2''=-(1/T^2)(1/C_{V1}+1/C_{V2})=-1/(T^2C_V),用 CV=CV1+CV2C_V=C_{V1}+C_{V2}。分布近似高斯,方差为 −k/(S1′′+S2′′)-k/(S_1''+S_2''),

⟨(δE1)2⟩=kT2CV\begin{align} \langle(\delta E_1)^2\rangle=kT^2C_V \end{align}

若子系统 2 远大于 1,CV≃CV1C_V\simeq C_{V1},退化成 ⟨(δE1)2⟩≃kT2CV1\langle(\delta E_1)^2\rangle\simeq kT^2C_{V1},也就是小系统的能量涨落。相对涨落随 NN 变小,取理想气体 CV=3Nk/2C_V=3Nk/2、E=3NkT/2E=3NkT/2,

δEE=kT2⋅3Nk/23NkT/2=23N\begin{align} \frac{\delta E}{E}=\frac{\sqrt{kT^2\cdot3Nk/2}}{3NkT/2} =\sqrt{\frac{2}{3N}} \end{align}

N∼1023N\sim10^{23} 时这个比值是 10−1210^{-12} 量级。结论是宏观系统的能量分配几乎确定,所以系综平均可以用最可几值代替,这是整章统计力学能成立的前提。

题型题面特征第一动作
算熵变给初末态,问 ΔS\Delta S 或总熵判可逆,不可逆就另设一条可逆路径
循环效率一圈闭合曲线,问 η\eta 或 WW圈上 ΔU=0\Delta U=0,WW 取面积,Q2Q_2 只加吸热段
Maxwell 套用问一个难测的偏导找它出现在哪个势的全微分里
相变潜热给 ll、Δv\Delta v、TT写 g1=g2g_1=g_2,沿相界线微分
配分函数给 ϵ(p)\epsilon(p),问热力学量先算 ZZ,再对 β\beta 和 VV 求导
能量均分给 HH,问 CVC_V 或 γ\gamma数二次型自由度 ff
涨落问平均值能不能当确定值把 log⁡[Γ1Γ2]\log[\Gamma_1\Gamma_2] 展开到二阶

易错点

先把最容易失分的三块写成散文,剩下的列在表里。

不可逆过程的熵变是第一个大坑。dS=dQ/T\dd S=\dbar Q/T 只在可逆时成立,写成 ΔS=∫dQ/T\Delta S=\int\dbar Q/T 对不可逆过程是错的,因为不可逆时这个积分小于 ΔS\Delta S。自由膨胀就是判据,dQ≡0\dbar Q\equiv0 然而 ΔS=Nkln⁡2\Delta S=Nk\ln2。办法是换路径,只认初末态。另一个常见漏洞是只算气体,忘了环境,其实要判断的是孤立系(气体加热库)的总熵,自由膨胀时热库没变,总熵才等于气体的熵增。

偏导的固定量是第二个大坑。(∂U/∂T)V(\partial U/\partial T)_V 与 (∂U/∂T)P(\partial U/\partial T)_P 是两个不同的量,前者是 CVC_V,后者对理想气体等于 CVC_V 但一般不等于。同理 (∂S/∂V)T=(∂P/∂T)V(\partial S/\partial V)_T=(\partial P/\partial T)_V 里的下标不能省,省了以后别人没法判断你指的是哪条路径。写偏导时养成每次都写上下标的习惯,考试时省下标扣的分比想不起来公式还多。

CVC_V 与 CPC_P 何时能互换是第三个。CP−CV=NkC_P-C_V=Nk 只对理想气体成立,把它当普遍结论套到 van der Waals 气体上会错。用熵方程时也一样,TdS=CVdT+T(∂P/∂T)VdVT\dd S=C_V\dd T+T(\partial P/\partial T)_V\dd V 的独立变量是 (T,V)(T,V),另一条 TdS=CPdT−T(∂V/∂T)PdPT\dd S=C_P\dd T-T(\partial V/\partial T)_P\dd P 的独立变量是 (T,P)(T,P),两条不能混着用,也不能把一条里的 CVC_V 直接换成 CPC_P。

坑为什么错正确做法
不可逆过程用 ΔS=∫dQ/T\Delta S=\int\dbar Q/T等号只对可逆成立,不可逆时积分偏小换一条同初末态的可逆路径
只算系统的熵熵不减是对孤立系说的系统与热库一起算
偏导不写固定量不同固定量给出不同的值每次写全下标
把 CP−CV=NkC_P-C_V=Nk 当普遍式它只在理想气体成立先看状态方程是不是 PV=NkTPV=NkT
Maxwell 关系的负号丢一半负号来自换变量时的符号翻转从势的全微分现场推
Γ\Gamma、Σ\Sigma、ω\omega 混用Γ\Gamma 是壳层体积,Σ\Sigma 是整块体积,ω\omega 是态密度算熵时三者可换,算体积时不能换
S=klog⁡ΓS=k\log\Gamma 直接取对数Γ\Gamma 带 (dp dq)N(\dd p\,\dd q)^N 的量纲先规定单位体积,或明确 SS 只定到差一常数
用 Σ∝VNE3N/2\Sigma\propto V^NE^{3N/2} 直接得熵得到 S∼Nklog⁡VS\sim Nk\log V,混同种气体时破坏广延求和中只有一项占主导这一步不能省;同种粒子加 1/N!1/N!
能量均分乱套⟨xi∂H/∂xj⟩=δijkT\langle x_i\partial H/\partial x_j\rangle=\delta_{ij}kT 是一般的,⟨H⟩=fkT/2\langle H\rangle=fkT/2 只对二次型只数二次型自由度
涨落公式里放 CPC_P微正则下 VV 固定,S′′=−1/(T2CV)S''=-1/(T^2C_V)只能用 CVC_V
判断违反 Kelvin 时漏看其他变化等温膨胀从单一热库吸热全变功,但体积变了先检查”唯一效果”是否成立
相变用 AA 最小相变在等温等压下进行用 GG 最小,得 g1=g2g_1=g_2
把 vdW 等温线的上升段当相变平台那段 ∂P/∂V>0\partial P/\partial V>0,压缩率为负用 Maxwell 等面积构造画水平线
E=−N∂log⁡Z/∂βE=-N\partial\log Z/\partial\beta 少负号ZZ 随 β\beta 减小,log⁡Z\log Z 对 β\beta 的导数为负记住 EE 必须为正
把 Δ\Delta 当微观小量Δ≪E\Delta\ll E 但 Δ\Delta 是宏观量不要把它和 N→∞N\to\infty 混为一谈
温标当成温度PV=F(T)PV=F(T) 对任何单值单调 FF 都成立温度由第零定律定,温标还要定单位与原点

自测清单

能独立答出下面这些,这两章就算过了一遍。答不出的先记下来,回去翻对应的小节。

  • 能说出第零定律的内容,并解释没有它温度为什么不是唯一确定的。
  • 能写出状态方程的一般形式,并说清完备集里参数个数由什么决定。
  • 能判断一个参数是广延还是强度,并说出这样分类在讨论子系统时的用处。
  • 能说清水银温度计这类方案为什么不普适,理想气体温标为什么好一些。
  • 能证明 TV=TPT_V=T_P,并指出哪一步用了 Boyle 定律。
  • 能独立推出 VγP=constV^\gamma P=\mathrm{const} 与 Vγ−1T=constV^{\gamma-1}T=\mathrm{const},并说明哪一步用了理想气体。
  • 能证卡诺效率是 1−T1/T21-T_1/T_2,并把它和绝对温标的定义对上。
  • 能用反证法证卡诺定理,说清取整数 N,N′N,N' 的作用和”唯一效果”在哪里用到。
  • 能默写 Kelvin 与 Clausius 两种表述,并能互相反推。
  • 能用无穷多小卡诺循环证 ∮revdQ/T=0\oint_{\mathrm{rev}}\dbar Q/T=0,并由此说明 SS 是态函数。
  • 能说清可逆与不可逆过程算 ΔS\Delta S 的差别,并举出 dQ=0\dbar Q=0 而 ΔS>0\Delta S>0 的例子。
  • 能把 TdS=dU+PdVT\dd S=\dd U+P\dd V 改写成两条熵方程,并说明各自的独立变量。
  • 能从四个势的全微分现场推出四条 Maxwell 关系,包括负号的位置。
  • 能用 AA 最小与 GG 最小分别处理等温等容和等温等压的平衡。
  • 能写出 Clapeyron 方程,并说明它为什么只对一级相变成立。
  • 能写出 van der Waals 方程,并说清 aa 与 bb 各自的物理来源。
  • 能由 (V−Vc)3=0(V-V_c)^3=0 三重根定出 TcT_c、PcP_c、VcV_c,并写出对应态方程。
  • 能用 Maxwell 等面积构造解释 vdW 等温线上的平台,分清亚稳段与不物理段。
  • 能写出第三定律,并说明 T→0T\to0 时 C→0C\to0 是怎么从 S=∫0TC dT/TS=\int_0^TC\,\dd T/T 得到的。
  • 能写出等概率先验假设,并说清 Δ\Delta 为什么必须引入。
  • 能说清 Γ\Gamma、Σ\Sigma、ω\omega 的定义与差别,以及 SS 用哪一个都行的理由。
  • 能推出 SS 是广延量,并指出”求和只有一项占主导”这一步为什么不可省。
  • 能用 Γ(E)\Gamma(E) 定义温度与压强,并证明 P=−(∂E/∂V)SP=-(\partial E/\partial V)_S。
  • 能独立走一遍能量均分的积分骨架,包括分部积分和边界项为零。
  • 能说明 ⟨H⟩=fkT/2\langle H\rangle=fkT/2 需要二次型,而 ⟨xi∂H/∂xj⟩=δijkT\langle x_i\partial H/\partial x_j\rangle=\delta_{ij}kT 不需要。
  • 能算出理想气体的 Σ(E)\Sigma(E),并得到 U=3NkT/2U=3NkT/2 与 PV=NkTPV=NkT。
  • 能说清 Gibbs 佯谬,说明 1/N!1/N! 修正的作用,以及 N!N! 与 NNN^N 两种取法的区别。
  • 能用 Lagrange 乘子推出 nˉl=gle−α−βϵl\bar n_l=g_le^{-\alpha-\beta\epsilon_l},并用二阶变分说明它是极大。
  • 能由 β=1/kT\beta=1/kT 把 ZZ 与 EE、PP、SS 连起来,并写出自由粒子的 Z=V(2πmkT)3/2Z=V(2\pi mkT)^{3/2}。
  • 能估算两个子系统能量分配的涨落,并说明它为什么正比于 kT2CV\sqrt{kT^2C_V}。

Comment评论

添加新评论