写在前面
这场期中考的范围是前两章,第一章热力学配第二章经典统计力学。第一章的 §1.10 有限时间热力学和第二章末尾关于负温度的讨论属于边角,主战场在 §1.2 到 §1.8 与 §2.2 到 §2.7。
这份材料的排法和课堂笔记不一样。课堂笔记顺着讲义节次走,用来第一遍学;这份按考试当天的动作排,先看全景把两章接起来,再用公式速查定位,然后靠必会推导和典型题型进入能落笔的状态,最后用易错点和自测清单收口。所以它不重述定义,只写考场上要调用的东西。
距考试还有六天,时间紧就按这个顺序走一遍。先扫一遍公式速查里的状态方程、熵方程、四个势、卡诺与微正则这几张表,把记号认熟;再啃必会推导里的卡诺效率与 S = k log Γ S=k\log\Gamma S = k log Γ 两条,这两条几乎每次都从里面出题;然后做典型题型的算熵变和循环效率,这两类是基础分;精力够的话把 Maxwell 关系套用和 van der Waals 那部分补上,那是拉分的地方。易错点放在考前最后一天看,那时候错题已经定型,看清单比看推导有用。
记号沿用前两章。自由能用 A A A ,玻尔兹曼常量写 k k k ,β = 1 / k T \beta=1/kT β = 1/ k T ,Γ \Gamma Γ 是 Γ \Gamma Γ 空间里的壳层体积,Σ \Sigma Σ 是 H < E H<E H < E 那整块体积,ω = ∂ Σ / ∂ E \omega=\partial\Sigma/\partial E ω = ∂ Σ/ ∂ E 是态密度。
全景
热力学那条线
起点是热平衡。两个系统热接触足够久之后变得一样冷热,这叫热平衡;第零定律说与第三个系统各自热平衡的两个系统彼此热平衡。这条定律的作用是把温度 T T T 变成一个有良好定义的东西,否则不同的系统对会给出不同的温度,T T T 就没法比较。有了 T T T ,再把它塞进完备集 ( T , x 1 , … , x M ) (T,x_1,\dots,x_M) ( T , x 1 , … , x M ) 里,就得到状态方程 f ( T , x 1 , … , x M ) = 0 f(T,x_1,\dots,x_M)=0 f ( T , x 1 , … , x M ) = 0 ,理想气体的 P V = N k T PV=NkT P V = N k T 是最简单的一例。温度和状态方程这两样东西一到位,后面的推导才有落脚点。
接着是第一定律。功 d W = P d V \dbar W=P\dd V d W = P d V 依赖路径,所以写成 d \dbar d 不写成 d \dd d ;热也一样。第一定律的真正内容是 Δ U = Δ Q − Δ W \Delta U=\Delta Q-\Delta W Δ U = Δ Q − Δ W 对给定初末态的所有变换都取同一个值,于是 d U \dd U d U 是恰当微分,U U U 是态函数。这一步把 d Q \dbar Q d Q 换成了 d U + P d V \dd U+P\dd V d U + P d V ,热就从”某种流进来的东西”变成了可以用状态量算的量。
第二定律管方向。Kelvin 表述堵住”从单一热库取热全部变功”,Clausius 表述堵住”热自发从冷到热”,两种表述可以互推,所以只取其一当公理。由它推出卡诺定理,两个给定温度之间任何热机的效率都不超过卡诺热机,而卡诺热机的效率是 1 − T 1 / T 2 1-T_1/T_2 1 − T 1 / T 2 。这条定理顺手给出一个不依赖工作物质的绝对温标,把 T T T 的定义从”选一支温度计”升级成”由效率比固定”。
再往下是熵。把任意循环用等温线和绝热线切成无穷多个小卡诺循环,每个小循环满足 Q j / T j + Q j + 1 / T j + 1 = 0 Q_j/T_j+Q_{j+1}/T_{j+1}=0 Q j / T j + Q j + 1 / T j + 1 = 0 ,相邻小循环的内部分界线互相抵消,剩下的就是 ∮ d Q / T ≤ 0 \oint\dbar Q/T\le0 ∮ d Q / T ≤ 0 ,可逆取等号。等号成立时 ∫ A B d Q / T \int_A^B\dbar Q/T ∫ A B d Q / T 与路径无关,于是 S S S 存在,且 T d S = d U + P d V T\dd S=\dd U+P\dd V T d S = d U + P d V 。这条式子把第一和第二定律装进一个微分式,剩下的全是它的改写。
改写走两条路。一条是 Legendre 变换,得到 A = U − T S A=U-TS A = U − T S 、H = U + P V H=U+PV H = U + P V 、G = H − T S G=H-TS G = H − T S ,四个势的全微分各给一条 Maxwell 关系,并且给出等温等容下 A A A 最小、等温等压下 G G G 最小的平衡判据。另一条是相变,把 G G G 最小用到两相共存上,得 g 1 = g 2 g_1=g_2 g 1 = g 2 ,沿相界线微分一次就是 Clapeyron 方程 d P / d T = Δ s / Δ v \dd P/\dd T=\Delta s/\Delta v d P / d T = Δ s /Δ v 。van der Waals 方程从另一头补上这块,它用 V − b V-b V − b 扣掉分子的硬球体积、用 P + a / V 2 P+a/V^2 P + a / V 2 扣掉器壁附近的吸引,等温线在低温下出现负压缩率的一段,Maxwell 等面积构造把那段换成一水平平台,平台的端点是两相共存。
最后是第三定律。熵的定义要接一个参考态,参考态选得不唯一,S S S 就定到差一个常数;第三定律说绝对零度时熵是普适常数,可以取零,于是 S S S 唯一了。代价是热容必须在 T → 0 T\to0 T → 0 时归零,因为 S = ∫ 0 T C d T / T S=\int_0^T C\,\dd T/T S = ∫ 0 T C d T / T 要收敛。
统计力学那条线
统计力学从另一个方向进来。等概率先验假设说孤立系统在平衡时,能量在 E E E 到 E + Δ E+\Delta E + Δ 之间的每个微观状态出现的概率相同,分布 ρ \rho ρ 在这层壳上是常数。由此得到微正则系综,壳层体积记为 Γ ( E ) \Gamma(E) Γ ( E ) 。Γ \Gamma Γ 量的是 Γ \Gamma Γ 空间里的体积,要数状态得先规定单位体积,这一点讲义特意提了一句。
熵的定义是 S = k log Γ ( E ) S=k\log\Gamma(E) S = k log Γ ( E ) 。它有两个必须验证的地方。广延性靠”把系统切成两个子系统后,能量分解 E = E 1 + E 2 E=E_1+E_2 E = E 1 + E 2 的求和中只有一项占主导”来保证,一项占主导才有 log Γ 1 Γ 2 = log Γ 1 + log Γ 2 \log\Gamma_1\Gamma_2=\log\Gamma_1+\log\Gamma_2 log Γ 1 Γ 2 = log Γ 1 + log Γ 2 ;第二定律靠 Σ ( E ) \Sigma(E) Σ ( E ) 随 V V V 不减来保证,因为孤立系统的 V V V 只能增大。
温度压强从偏导里长出来。定义 ( ∂ S / ∂ E ) V = 1 / T (\partial S/\partial E)_V=1/T ( ∂ S / ∂ E ) V = 1/ T ,再用 P = T ( ∂ S / ∂ V ) E P=T(\partial S/\partial V)_E P = T ( ∂ S / ∂ V ) E ,代进 d S \dd S d S 的两个偏导就回到 d E = T d S − P d V \dd E=T\dd S-P\dd V d E = T d S − P d V 。这一步的意义是热力学第一定律可以从这里推出来。讲义特意提醒,第一定律本身在统计力学里仍然是一条假设。两个子系统能量交换到平衡时 ( ∂ S 1 / ∂ E 1 ) = ( ∂ S 2 / ∂ E 2 ) (\partial S_1/\partial E_1)=(\partial S_2/\partial E_2) ( ∂ S 1 / ∂ E 1 ) = ( ∂ S 2 / ∂ E 2 ) ,正是第零定律的 T 1 = T 2 T_1=T_2 T 1 = T 2 ,所以温度的定义和热力学是自洽的。
接着是能量均分。对任意 Hamiltonian,⟨ x i ∂ H / ∂ x j ⟩ = δ i j k T \langle x_i\partial H/\partial x_j\rangle=\delta_{ij}kT ⟨ x i ∂ H / ∂ x j ⟩ = δ ij k T ,用的是把壳层积分写成差分、再分部积分、边界项为零这一套。若 H H H 是二次型,每条二次型自由度分到 k T / 2 kT/2 k T /2 ,于是 ⟨ H ⟩ = f k T / 2 \langle H\rangle=fkT/2 ⟨ H ⟩ = f k T /2 ,单原子理想气体的 C V = 3 N k / 2 C_V=3Nk/2 C V = 3 N k /2 就是这么来的。
理想气体把整条链跑通一遍。Σ ( E ) = C 3 N [ V ( 2 m E ) 3 / 2 ] N \Sigma(E)=C_{3N}[V(2mE)^{3/2}]^N Σ ( E ) = C 3 N [ V ( 2 m E ) 3/2 ] N ,S = k log Σ S=k\log\Sigma S = k log Σ 算出后反解出 E = E ( S , V ) E=E(S,V) E = E ( S , V ) ,微分得 T = 2 U / ( 3 N k ) T=2U/(3Nk) T = 2 U / ( 3 N k ) 即 U = 3 N k T / 2 U=3NkT/2 U = 3 N k T /2 ,以及 P = N k T / V P=NkT/V P = N k T / V 。这两个结果和第一章的实验结论对上,说明 S = k log Γ S=k\log\Gamma S = k log Γ 这个定义没选错。混合同种气体时熵会多出一项 N k log V Nk\log V N k log V ,破坏广延性,加 Gibbs 因子 1 / N ! 1/N! 1/ N ! 修掉。
最后落到具体分布。准独立粒子的情形下,把 μ \mu μ 空间切成小格,用占据数 { n l } \{n_l\} { n l } 描述,Γ \Gamma Γ 空间里对应 { n l } \{n_l\} { n l } 的体积是 N ! / ∏ n l ! ⋅ ∏ g l n l N!/\prod n_l!\cdot\prod g_l^{n_l} N ! / ∏ n l ! ⋅ ∏ g l n l 。这个量在给定 N N N 与 E E E 下取极大,用两个 Lagrange 乘子 α , β \alpha,\beta α , β 得 n ˉ l = g l e − α − β ϵ l \bar n_l=g_le^{-\alpha-\beta\epsilon_l} n ˉ l = g l e − α − β ϵ l ,即 Maxwell-Boltzmann 分布。把求和写成单粒子配分函数 Z = ∑ l e − β ϵ l Z=\sum_le^{-\beta\epsilon_l} Z = ∑ l e − β ϵ l ,热力学量变成对 β \beta β 和 V V V 的求导,E = − N ∂ log Z / ∂ β E=-N\partial\log Z/\partial\beta E = − N ∂ log Z / ∂ β ,P = ( N / β ) ∂ log Z / ∂ V P=(N/\beta)\partial\log Z/\partial V P = ( N / β ) ∂ log Z / ∂ V ,温度由 β = 1 / k T \beta=1/kT β = 1/ k T 定。
两章怎么接起来
统计力学在这里干的事,是给热力学里的 T T T 、S S S 、P P P 配一套微观算法。形式上两边的接口只有一条式子,( ∂ S / ∂ E ) V = 1 / T (\partial S/\partial E)_V=1/T ( ∂ S / ∂ E ) V = 1/ T 和 P = T ( ∂ S / ∂ V ) E P=T(\partial S/\partial V)_E P = T ( ∂ S / ∂ V ) E 合起来就是 d E = T d S − P d V \dd E=T\dd S-P\dd V d E = T d S − P d V ,和第一章的热力学基本方程逐字相同,只是 S S S 换成了 k log Γ k\log\Gamma k log Γ 。讲义说 ( ∂ S / ∂ E ) = 1 / T (\partial S/\partial E)=1/T ( ∂ S / ∂ E ) = 1/ T 本身就是第一章的 Maxwell 关系之一,这话的意思是两边不是两套理论,是同一套关系的两种算法。
考试里这个接口常见的用法有两种。一种是给你一个微观模型,让你算 Γ \Gamma Γ 或 Z Z Z ,再问 U U U 、C V C_V C V 、P P P 、S S S ,这时候走 S = k log Γ S=k\log\Gamma S = k log Γ 或者走 Z Z Z 都行,走 Z Z Z 快一些。另一种是给你一个宏观关系,问它在统计上对应什么,比如问能量均分为什么给出 U = 3 N k T / 2 U=3NkT/2 U = 3 N k T /2 ,那就要把 Σ ∝ E 3 N / 2 \Sigma\propto E^{3N/2} Σ ∝ E 3 N /2 这一步补上。
公式速查
温度、状态方程与温标
温度这一块考得不多,但定义链条容易被问。第零定律不只是一句常识,它的作用是保证”与 C C C 热平衡”这个关系可以传递,从而 T T T 唯一。温标的做法分三步,选一个系统加一个随温度单调变化的性质,假定这个性质与 T T T 的关系(通常取线性),再定原点和单位。第三步依赖实验,所以水银温度计这种方案不普适,只在有限温区好用。
理想气体温标好在它不需要指定具体气体。固定 V V V 时设 T V = 100 ( P − P i ) / ( P f − P i ) T_V=100(P-P_i)/(P_f-P_i) T V = 100 ( P − P i ) / ( P f − P i ) ,固定 P P P 时设 T P = 100 ( V − V i ) / ( V f − V i ) T_P=100(V-V_i)/(V_f-V_i) T P = 100 ( V − V i ) / ( V f − V i ) ,由 Boyle 定律 P f V i = V f P i P_fV_i=V_fP_i P f V i = V f P i 可以证明 T V = T P T_V=T_P T V = T P ,证明的办法是把 P V = a V T V + b V PV=a_VT_V+b_V P V = a V T V + b V 和 P V = a P T P + b P PV=a_PT_P+b_P P V = a P T P + b P 在两条等容、等压线上比较系数,得到 a V = a P a_V=a_P a V = a P 、b V = b P b_V=b_P b V = b P 。
关系 式子 适用条件 第零定律 A ∼ C A\sim C A ∼ C 且 B ∼ C B\sim C B ∼ C 则 A ∼ B A\sim B A ∼ B 热平衡的传递性,定义 T T T 的前提 状态方程 f ( T , x 1 , … , x M ) = 0 f(T,x_1,\dots,x_M)=0 f ( T , x 1 , … , x M ) = 0 平衡态,{ x i } \{x_i\} { x i } 是完备集,个数由系统固定 理想气体温标 P V = N k T PV=NkT P V = N k T 稀薄气体;T T T 用理想气体温标 一般温标 P V = F ( T ) PV=F(T) P V = F ( T ) F F F 单值单调,任何一种选择都构成温标绝对温标 T 1 / T 2 = 1 − η T_1/T_2=1-\eta T 1 / T 2 = 1 − η 卡诺热机,不依赖工作物质 绝对温标单位 T n + 1 − T n = 1 K T_{n+1}-T_n=1\ \mathrm{K} T n + 1 − T n = 1 K 固定每一级卡诺循环的功 强度与广延 V → λ V V\to\lambda V V → λV 时广延量乘 λ \lambda λ ,强度量不变热力学极限,可分系统才可加
第一定律与熵
这张表是整场考试用得最多的一张。要留意两件事,一是 d \dbar d 与 d \dd d 的区别,功和热不是态函数,只有 U U U 、S S S 、H H H 、A A A 、G G G 是;二是每个式子背后的独立变量,熵方程有两种写法,分别对应 ( T , V ) (T,V) ( T , V ) 和 ( T , P ) (T,P) ( T , P ) 两套独立变量,混用必错。
关系 式子 适用条件 第一定律 Δ U = Δ Q − Δ W \Delta U=\Delta Q-\Delta W Δ U = Δ Q − Δ W 任何变换,不要求准静态 无穷小形式 d U = d Q − d W \dd U=\dbar Q-\dbar W d U = d Q − d W 任何变换,d \dbar d 表示非恰当 准静态功 d W = P d V \dbar W=P\dd V d W = P d V 准静态,且只有 P d V P\dd V P d V 这一种功 热容 C V = ( ∂ U / ∂ T ) V C_V=(\partial U/\partial T)_V C V = ( ∂ U / ∂ T ) V ,C P = ( ∂ H / ∂ T ) P C_P=(\partial H/\partial T)_P C P = ( ∂ H / ∂ T ) P 任何系统,下标是路径 理想气体热容差 C P − C V = N k C_P-C_V=Nk C P − C V = N k 理想气体 热力学基本方程 T d S = d U + P d V T\dd S=\dd U+P\dd V T d S = d U + P d V 可逆变换,或连接两个平衡态 熵方程 T d S = C V d T + T ( ∂ P / ∂ T ) V d V T\dd S=C_V\dd T+T(\partial P/\partial T)_V\dd V T d S = C V d T + T ( ∂ P / ∂ T ) V d V 可逆,独立变量取 ( T , V ) (T,V) ( T , V ) 熵方程 T d S = C P d T − T ( ∂ V / ∂ T ) P d P T\dd S=C_P\dd T-T(\partial V/\partial T)_P\dd P T d S = C P d T − T ( ∂ V / ∂ T ) P d P 可逆,独立变量取 ( T , 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 ( ∂ U / ∂ V ) T = T ( ∂ P / ∂ T ) V − P 由 d S \dd S d S 恰当性得出,任何物质 可逆熵变 S ( B ) − S ( A ) = ∫ A B d Q / T S(B)-S(A)=\int_A^B\dbar Q/T S ( B ) − S ( A ) = ∫ A B d Q / T 沿可逆路径积分 不可逆熵变 S ( B ) − S ( A ) ≥ ∫ A B d Q / T S(B)-S(A)\ge\int_A^B\dbar Q/T S ( B ) − S ( A ) ≥ ∫ A B d Q / T 任意变换,等号只在可逆时 Clausius 定理 ∮ d Q / T ≤ 0 \oint\dbar Q/T\le0 ∮ d Q / T ≤ 0 循环,可逆取等号 绝热系统 Δ S ≥ 0 \Delta S\ge0 Δ S ≥ 0 绝热,可逆取等号 等温膨胀 Δ S = N k ln ( V 2 / V 1 ) \Delta S=Nk\ln(V_2/V_1) Δ S = N k ln ( V 2 / V 1 ) 理想气体,T T T 不变 任意两态 Δ S = C V ln ( T 2 / T 1 ) + N k ln ( V 2 / V 1 ) \Delta S=C_V\ln(T_2/T_1)+Nk\ln(V_2/V_1) Δ S = C V ln ( T 2 / T 1 ) + N k ln ( V 2 / V 1 ) 理想气体,C V C_V C 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 ) z ( ∂ y / ∂ z ) x ( ∂ z / ∂ 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 ( ∂ x / ∂ y ) w ( ∂ y / ∂ z ) w = ( ∂ x / ∂ z ) w 固定同一个 w w w
热力学势与 Maxwell 关系
四个势的记忆法则是自变量成对换。U U U 配 ( S , V ) (S,V) ( S , V ) ,A A A 配 ( T , V ) (T,V) ( T , V ) ,H H H 配 ( S , P ) (S,P) ( S , P ) ,G G G 配 ( T , P ) (T,P) ( T , P ) ,负号出现在换掉的那个变量上。把势写成它自己那对变量的函数(特性函数)之后,其余所有量都能靠微分得到,这就是 Masius 定理。
关系 式子 适用条件 自由能 A = U − T S A=U-TS A = U − T S ,d A = − S d T − P d V \dd A=-S\dd T-P\dd V d A = − S d T − P d V A = A ( T , V ) A=A(T,V) A = A ( T , V ) ,可逆焓 H = U + P V H=U+PV H = U + P V ,d H = T d S + V d P \dd H=T\dd S+V\dd P d H = T d S + V d P H = H ( S , P ) H=H(S,P) H = H ( S , P ) ,可逆Gibbs 势 G = H − T S G=H-TS G = H − T S ,d G = − S d T + V d P \dd G=-S\dd T+V\dd P d G = − S d T + V d P G = G ( T , P ) G=G(T,P) G = G ( T , P ) ,可逆Maxwell 关系 ( ∂ T / ∂ V ) S = − ( ∂ P / ∂ S ) V (\partial T/\partial V)_S=-(\partial P/\partial S)_V ( ∂ T / ∂ V ) S = − ( ∂ P / ∂ S ) V 由 d U \dd U d U 的交叉偏导相等 Maxwell 关系 ( ∂ S / ∂ V ) T = ( ∂ P / ∂ T ) V (\partial S/\partial V)_T=(\partial P/\partial T)_V ( ∂ S / ∂ V ) T = ( ∂ P / ∂ T ) V 由 d A \dd A d A ,用过最多的一条 Maxwell 关系 ( ∂ T / ∂ P ) S = ( ∂ V / ∂ S ) P (\partial T/\partial P)_S=(\partial V/\partial S)_P ( ∂ T / ∂ P ) S = ( ∂ V / ∂ S ) P 由 d H \dd H d H Maxwell 关系 − ( ∂ S / ∂ P ) T = ( ∂ V / ∂ T ) P -(\partial S/\partial P)_T=(\partial V/\partial T)_P − ( ∂ S / ∂ P ) T = ( ∂ V / ∂ T ) P 由 d G \dd G d G ,负号在左边 特性函数 S = − ( ∂ A / ∂ T ) V S=-(\partial A/\partial T)_V S = − ( ∂ A / ∂ T ) V ,P = − ( ∂ A / ∂ V ) T P=-(\partial A/\partial V)_T P = − ( ∂ A / ∂ V ) T A = A ( T , V ) A=A(T,V) A = A ( T , V ) 已知由 A A A 求 C V C_V C V U = A − T ( ∂ A / ∂ T ) V U=A-T(\partial A/\partial T)_V U = A − T ( ∂ A / ∂ T ) V ,C V = − T ( ∂ 2 A / ∂ T 2 ) V C_V=-T(\partial^2A/\partial T^2)_V C V = − T ( ∂ 2 A / ∂ T 2 ) V A = A ( T , V ) A=A(T,V) A = A ( T , V ) 已知等温等容平衡 A A A 不增,平衡时最小d W = 0 \dbar W=0 d W = 0 且 T T T 固定等温等压平衡 G G G 不增,平衡时最小T , P T,P T , P 固定两子系统平衡 ( ∂ A 1 / ∂ V 1 ) T = ( ∂ A 2 / ∂ V 2 ) T (\partial A_1/\partial V_1)_T=(\partial A_2/\partial V_2)_T ( ∂ A 1 / ∂ V 1 ) T = ( ∂ A 2 / ∂ V 2 ) T 即 P 1 = P 2 P_1=P_2 P 1 = P 2 总体积固定,活塞可动
卡诺循环
关系 式子 适用条件 效率 η = W / Q 2 = 1 − Q 1 / Q 2 \eta=W/Q_2=1-Q_1/Q_2 η = W / Q 2 = 1 − Q 1 / Q 2 任何热机,Q 2 Q_2 Q 2 吸热 Q 1 Q_1 Q 1 放热 卡诺效率 η = 1 − T 1 / T 2 \eta=1-T_1/T_2 η = 1 − T 1 / T 2 可逆循环,两个热库,T T T 为绝对温标 卡诺定理 η ≤ 1 − T 1 / T 2 \eta\le1-T_1/T_2 η ≤ 1 − T 1 / T 2 任何热机 可逆循环的热比 Q 2 / T 2 = Q 1 / T 1 Q_2/T_2=Q_1/T_1 Q 2 / T 2 = Q 1 / T 1 可逆,两个热库 绝热方程 V γ P = c o n s t V^\gamma P=\mathrm{const} V γ P = const ,V γ − 1 T = c o n s t V^{\gamma-1}T=\mathrm{const} V γ − 1 T = const 理想气体,准静态绝热,C V C_V C V 常数 绝热指数 γ = 1 + N k / C V \gamma=1+Nk/C_V γ = 1 + N k / C V 理想气体,C V C_V C V 常数 等温吸热 Q 2 = N k T 2 ln ( V b / V a ) Q_2=NkT_2\ln(V_b/V_a) Q 2 = N k T 2 ln ( V b / V a ) 理想气体,等温可逆 卡诺推论 同两热库间所有卡诺热机效率相同 效率与 W W W 、与工作物质都无关
换到 T T T -S S S 平面上看同一个循环,图形会变成一个矩形,这件事本身就值一张图。两个绝热段是等熵的,所以在图上是两条竖线;两个等温段温度不变,是两条横线。
热量在 T T T -S S S 图上是线下的面积,因为 Q = ∫ T d S Q=\int T\dd S Q = ∫ T d S 。顶边下方面积是吸热 Q 2 Q_2 Q 2 ,底边下方面积是放热 Q 1 Q_1 Q 1 ,两者相减就是矩形自己那块灰面积,也就是净功。于是效率一行就出来,η = 1 − Q 1 / Q 2 = 1 − T 1 / T 2 \eta=1-Q_1/Q_2=1-T_1/T_2 η = 1 − Q 1 / Q 2 = 1 − T 1 / T 2 ,连绝热方程都不用碰。这也是 T T T -S S S 图在复习时比 p p p -V V V 图顺手的原因,熵是态函数,绝热线在图上必须是竖线,循环的面积关系一目了然。
相变与 van der Waals
相变这一块,两相共存的条件是 Gibbs 势最小,而 G = m 1 g 1 + m 2 g 2 G=m_1g_1+m_2g_2 G = m 1 g 1 + m 2 g 2 ,对质量分配求变分得 g 1 = g 2 g_1=g_2 g 1 = g 2 。相变在等温等压下进行,对应的判据是 G G G 最小,等温等容才轮到 A A A 。g g g 是单位质量的 Gibbs 势,也叫化学势。
关系 式子 适用条件 相平衡 g 1 ( P , T ) = g 2 ( P , T ) g_1(P,T)=g_2(P,T) g 1 ( P , T ) = g 2 ( P , T ) 一阶相变,等温等压 单位质量微分 d g = − s d T + v d P \dd g=-s\dd T+v\dd P d g = − s d T + v d P 单相内部 Clapeyron 方程 d P d T = Δ s Δ v = l T Δ v \dfrac{\dd P}{\dd T}=\dfrac{\Delta s}{\Delta v}=\dfrac{l}{T\Delta v} d T d P = Δ v Δ s = T Δ v l 一阶相变,沿相界线 潜热 l = T Δ s l=T\Delta s l = T Δ s 一阶相变 一级判据 Δ v ≠ 0 \Delta v\ne0 Δ v = 0 或 Δ s ≠ 0 \Delta s\ne0 Δ s = 0 ,g g g 的一阶导不连续一级相变 二级相变 Δ s = Δ v = 0 \Delta s=\Delta v=0 Δ s = Δ v = 0 ,Clapeyron 失效临界点及以上 三相点 三条 P ( T ) P(T) P ( T ) 相界线的交点 三相共存,P , T P,T P , T 都定死 临界点 相界线的终点 ( T c , P c ) (T_c,P_c) ( T c , P c ) 两相差别消失,一级转二级 van der Waals 方程 ( P + a / V 2 ) ( V − b ) = N k T (P+a/V^2)(V-b)=NkT ( P + a / V 2 ) ( V − b ) = N k T 实际气体的简单模型 vdW 临界量 N k T c = 8 a / 27 b NkT_c=8a/27b N k T c = 8 a /27 b ,P c = a / 27 b 2 P_c=a/27b^2 P c = a /27 b 2 ,V c = 3 b V_c=3b V c = 3 b 由 ( V − V c ) 3 = 0 (V-V_c)^3=0 ( V − V c ) 3 = 0 三重根定出 对应态方程 ( V ˉ − 1 / 3 ) ( P ˉ + 3 / V ˉ 2 ) = 8 T ˉ / 3 (\bar V-1/3)(\bar P+3/\bar V^2)=8\bar T/3 ( V ˉ − 1/3 ) ( P ˉ + 3/ V ˉ 2 ) = 8 T ˉ /3 T ˉ = T / T c \bar T=T/T_c T ˉ = T / T c 等,不含物质常数内压 ( ∂ U / ∂ V ) T = a / V 2 (\partial U/\partial V)_T=a/V^2 ( ∂ U / ∂ V ) T = a / V 2 vdW 气体,推导见后面「Maxwell 关系套用」
van der Waals 的两处修正各有物理来源。硬球排斥让每个分子可用的体积从 V V V 减到 V − b V-b V − b ,b b b 是物质常数;吸引让器壁附近的分子被往回拉,器壁感受到的压强比动理学压强小,减少的量正比于器壁附近分子层与内部的对数,约 N 2 / V 2 N^2/V^2 N 2 / V 2 ,写成 a / V 2 a/V^2 a / V 2 。把修正后的假设 P k i n e t i c V e f f = N k T P_{\mathrm{kinetic}}V_{\mathrm{eff}}=NkT P kinetic V eff = N k T 写出来就是上面那条方程。
图中蓝线是 T < T c T<T_c T < T c 的一条 van der Waals 等温线,在点 1 与点 2 之间它有一段压强随体积上升,压缩率为负,物理上不允许。红线的画法是让两块阴影面积相等,即 Maxwell 等面积构造,用它代替蓝色那一段。红线左端是液相、右端是气相,中间是两相共存,P P P 与 T T T 都保持不变。蓝色曲线上,从点 1 走到极小值这一段的压强低于共存压强,是亚稳的过热液体;从极大值走到点 2 这一段的压强高于共存压强,是亚稳的过冷气体;极小到极大之间那一段就是图里的橙色,压缩率为负,属于不稳定区,既不算液体也不算气体。两个亚稳段的自由能都比共存态高,一有凝结核或者气泡就会跳到红线上。灰虚线是 T > T c T>T_c T > T c 的等温线,单调下降,没有平台,两相也就分不出来了。
相图里的三条线各有来历。固液线、液气线、固气线在一点交汇,那点是三相点,三个相同时共存,P P P 和 T T T 都被定死,所以三相点可以拿来做温标的固定点。液气线的终点是临界点,越过它之后液体和气体的差别消失,转变从一级变成二级。
三条线的斜率都由 Clapeyron 方程决定,符号看 Δ v \Delta v Δ v 。从三相点向左下延伸的那条是固气线,也叫升华线。固液线的 Δ v \Delta v Δ v 很小,所以它几乎竖直,也就意味着改变压强对熔点的影响远小于对沸点的影响。
微正则系综
Γ \Gamma Γ 与 Σ \Sigma Σ 的区别在考卷上很爱出。Σ ( E ) \Sigma(E) Σ ( E ) 是 H < E H<E H < E 那整块体积,Γ ( E ) \Gamma(E) Γ ( E ) 是 E E E 到 E + Δ E+\Delta E + Δ 那层壳,ω = ∂ Σ / ∂ E \omega=\partial\Sigma/\partial E ω = ∂ Σ/ ∂ E 是单位能量的体积。三者在 S = k log S=k\log S = k log 里可以互换,代价只是差一个 O ( log N ) O(\log N) O ( log N ) 的加性项,所以算熵时挑好算的那个。
以 H = p 2 + q 2 H=p^2+q^2 H = p 2 + q 2 为例,H < E H<E H < E 是一张圆盘,半径 E \sqrt{E} E ;壳层就是圆环,宽度对应 Δ \Delta Δ 。Γ \Gamma Γ 是圆环的面积,Σ \Sigma Σ 是整张圆盘的面积,ω \omega ω 是周长。这个例子里圆环面积正比于 Δ \Delta Δ 、与 E E E 无关,圆盘面积正比于 E E E ,取对数只差一个常数;推广到 N N N 个自由度,差别是 O ( log N ) O(\log N) O ( log N ) 量级,不影响热力学量,这就是三者可以互换的来源。
关系 式子 适用条件 壳层体积 Γ ( E ) = ∫ E < H < E + Δ d p d q \Gamma(E)=\int_{E<H<E+\Delta}\dd p\,\dd q Γ ( E ) = ∫ E < H < E + Δ d p d q 微正则,Δ ≪ E \Delta\ll E Δ ≪ E 但 Δ \Delta Δ 是宏观量 累计体积 Σ ( E ) = ∫ H < E d p d q \Sigma(E)=\int_{H<E}\dd p\,\dd q Σ ( E ) = ∫ H < E d p d q 微正则 态密度 ω ( E ) = ∂ Σ ( E ) / ∂ E \omega(E)=\partial\Sigma(E)/\partial E ω ( E ) = ∂ Σ ( E ) / ∂ E ,Γ = ω Δ \Gamma=\omega\Delta Γ = ω Δ 微正则 熵 S ( E , V ) = k log Γ ( E ) S(E,V)=k\log\Gamma(E) S ( E , V ) = k log Γ ( E ) 微正则,Γ \Gamma Γ 是体积需先定单位体积 三式等价 S = k log Γ ≃ k log ω ≃ k log Σ S=k\log\Gamma\simeq k\log\omega\simeq k\log\Sigma S = k log Γ ≃ k log ω ≃ k log Σ 差 O ( log N ) O(\log N) O ( log N ) 的加性常数 温度 ( ∂ S / ∂ E ) V = 1 / T (\partial S/\partial E)_V=1/T ( ∂ S / ∂ E ) V = 1/ T 微正则 压强 P = T ( ∂ S / ∂ V ) E P=T(\partial S/\partial V)_E P = T ( ∂ S / ∂ V ) E 微正则 等价形式 P = − ( ∂ E / ∂ V ) S P=-(\partial E/\partial V)_S P = − ( ∂ E / ∂ V ) S 由 d E = T d S − P d V \dd E=T\dd S-P\dd V d E = T d S − P d V 直接读 基本方程 d E = T d S − P d V \dd E=T\dd S-P\dd V d E = T d S − P d V 只有两个独立参数 广延性 S = S 1 ( E ˉ 1 ) + S 2 ( E ˉ 2 ) S=S_1(\bar E_1)+S_2(\bar E_2) S = S 1 ( E ˉ 1 ) + S 2 ( E ˉ 2 ) ,E ˉ 1 + E ˉ 2 = E \bar E_1+\bar E_2=E E ˉ 1 + E ˉ 2 = E 求和只有一项占主导 平衡条件 ( ∂ S 1 / ∂ E 1 ) = ( ∂ S 2 / ∂ E 2 ) (\partial S_1/\partial E_1)=(\partial S_2/\partial E_2) ( ∂ S 1 / ∂ E 1 ) = ( ∂ S 2 / ∂ E 2 ) 即 T 1 = T 2 T_1=T_2 T 1 = T 2 ,第零定律
能量均分与理想气体
关系 式子 适用条件 均分定理 ⟨ x i ∂ H ∂ x j ⟩ = δ i j k T \left\langle x_i\dfrac{\partial H}{\partial x_j}\right\rangle=\delta_{ij}kT ⟨ x i ∂ x j ∂ H ⟩ = δ ij k T x i x_i x i 是 p p p 或 q q q ,边界项可忽略二次型结果 ⟨ H ⟩ = f k T / 2 \langle H\rangle=fkT/2 ⟨ H ⟩ = f k T /2 H = ∑ A i P i 2 + ∑ B i Q i 2 H=\sum A_iP_i^2+\sum B_iQ_i^2 H = ∑ A i P i 2 + ∑ B i Q i 2 单原子气体 C V = 3 N k / 2 C_V=3Nk/2 C V = 3 N k /2 ,γ = 5 / 3 \gamma=5/3 γ = 5/3 只有平动,f = 3 f=3 f = 3 谐振子 ⟨ H ⟩ = k T \langle H\rangle=kT ⟨ H ⟩ = k T ,C V = N k C_V=Nk C V = N k H = p 2 / 2 m + ω 2 q 2 / 2 H=p^2/2m+\omega^2q^2/2 H = p 2 /2 m + ω 2 q 2 /2 ,f = 2 f=2 f = 2 态数目 Σ ( E ) = C 3 N [ V ( 2 m E ) 3 / 2 ] N \Sigma(E)=C_{3N}[V(2mE)^{3/2}]^N Σ ( E ) = C 3 N [ V ( 2 m E ) 3/2 ] N 理想气体,Ω n ( R ) = C n R n \Omega_n(R)=C_nR^n Ω n ( R ) = C n R n 熵(未加修正) S = N k log [ V ( 4 π m E / 3 N ) 3 / 2 ] + 3 N k / 2 S=Nk\log[V(4\pi mE/3N)^{3/2}]+3Nk/2 S = N k log [ V ( 4 π m E /3 N ) 3/2 ] + 3 N k /2 无 Gibbs 因子 熵(加修正) S = N k log [ ( V / N ) ( 4 π m E / 3 N ) 3 / 2 ] + 5 N k / 2 S=Nk\log[(V/N)(4\pi mE/3N)^{3/2}]+5Nk/2 S = N k log [( V / N ) ( 4 π m E /3 N ) 3/2 ] + 5 N k /2 同种粒子需要 1 / N ! 1/N! 1/ N ! 内能 T = ( ∂ U / ∂ S ) V = 2 U / ( 3 N k ) T=(\partial U/\partial S)_V=2U/(3Nk) T = ( ∂ U / ∂ S ) V = 2 U / ( 3 N k ) ,U = 3 N k T / 2 U=3NkT/2 U = 3 N k T /2 理想气体 压强 P = − ( ∂ U / ∂ V ) S = N k T / V P=-(\partial U/\partial V)_S=NkT/V P = − ( ∂ U / ∂ V ) S = N k T / V 理想气体
配分函数与 MB 分布
关系 式子 适用条件 最可几分布 n ˉ l = g l e − α − β ϵ l \bar n_l=g_le^{-\alpha-\beta\epsilon_l} n ˉ l = g l e − α − β ϵ l 准独立粒子,N N N 很大 归一化 N = e − α Z N=e^{-\alpha}Z N = e − α Z ,即 α = log ( Z / N ) \alpha=\log(Z/N) α = log ( Z / N ) 由 ∑ n ˉ l = N \sum\bar n_l=N ∑ n ˉ l = N 单粒子配分函数 Z ( β , y ) = ∑ l e − β ϵ l = ∫ d p d q e − β ϵ Z(\beta,y)=\sum_le^{-\beta\epsilon_l}=\int\dd p\,\dd q\,e^{-\beta\epsilon} Z ( β , y ) = ∑ l e − β ϵ l = ∫ d p d q e − β ϵ 准独立,ϵ ≤ E \epsilon\le E ϵ ≤ E 的边界可忽略 内能 E = − N ∂ log Z ∂ β E=-N\dfrac{\partial\log Z}{\partial\beta} E = − N ∂ β ∂ log Z 准独立 广义力 Y k = N β ∂ log Z ∂ y k Y_k=\dfrac{N}{\beta}\dfrac{\partial\log Z}{\partial y_k} Y k = β N ∂ y k ∂ log Z y k y_k y k 为外参量压强 P = N β ∂ log Z ∂ V P=\dfrac{N}{\beta}\dfrac{\partial\log Z}{\partial V} P = β N ∂ V ∂ log Z y = V y=V y = V ,要把积分限随 V V V 的依赖算进去温度 β = 1 / k T \beta=1/kT β = 1/ k T 由 d Q = T d S \dbar Q=T\dd S d Q = T d S 定出 熵 S = N k ( log Z − β ∂ log Z ∂ β ) + c o n s t S=Nk\left(\log Z-\beta\dfrac{\partial\log Z}{\partial\beta}\right)+\mathrm{const} S = N k ( log Z − β ∂ β ∂ log Z ) + const 与微正则差 Gibbs 因子 热量 d Q = d E + ∑ k Y k d y k \dbar Q=\dd E+\sum_kY_k\dd y_k d Q = d E + ∑ k Y k d y k 第一定律,y y y 固定时 d Q = ∑ ϵ l d n ˉ l \dbar Q=\sum\epsilon_l\dd\bar n_l d Q = ∑ ϵ l d n ˉ l 自由粒子 Z = V ( 2 π m k T ) 3 / 2 Z=V(2\pi mkT)^{3/2} Z = V ( 2 π mk T ) 3/2 ϵ = p 2 / 2 m \epsilon=p^2/2m ϵ = p 2 /2 m ,无外势理想气体 E = 3 N k T / 2 E=3NkT/2 E = 3 N k T /2 ,P = N k T / V P=NkT/V P = N k T / V ,C V = 3 N k / 2 C_V=3Nk/2 C V = 3 N k /2 由上面那条 Z Z Z 直接求导
涨落
关系 式子 适用条件 相对涨落 ⟨ f 2 ⟩ − ⟨ f ⟩ 2 ⟨ f ⟩ 2 ≪ 1 \dfrac{\langle f^2\rangle-\langle f\rangle^2}{\langle f\rangle^2}\ll1 ⟨ f ⟩ 2 ⟨ f 2 ⟩ − ⟨ f ⟩ 2 ≪ 1 大系统;热力学极限下趋于零 最可几值 ⟨ f ⟩ ≃ f 0 \langle f\rangle\simeq f_0 ⟨ f ⟩ ≃ f 0 上式成立;长程相互作用或连续相变附近失效 能量分配 ⟨ ( δ E 1 ) 2 ⟩ = k T 2 C V \langle(\delta E_1)^2\rangle=kT^2C_V ⟨( δ E 1 ) 2 ⟩ = k T 2 C V 两子系统交换能量,峰宽二阶展开 相对量级 δ E / E ∼ 1 / N \delta E/E\sim1/\sqrt{N} δ E / E ∼ 1/ N 广延系统
必会推导
卡诺定理与卡诺效率
推导分两段,先算理想气体卡诺热机的效率,再证任何热机都超不过它。两段的结论合起来才是卡诺定理。
理想气体卡诺循环由两条等温线和两条绝热线围成,按 a → b → c → d → a a\to b\to c\to d\to a a → b → c → d → a 走,T 2 T_2 T 2 是高温,T 1 T_1 T 1 是低温。
图中 a → b a\to b a → b 与 c → d c\to d c → d 是两条等温线,b → c b\to c b → c 与 d → a d\to a d → a 是两条绝热线,四条线围出的就是卡诺循环。同一体积处两条线的陡峭程度不一样,等温线的斜率是 − P / V -P/V − P / V ,绝热线的是 − γ P / V -\gamma P/V − γ P / V ,比值正好是 γ \gamma γ 。这个差别在图上表现为绝热线离原点更近,也是判断一段曲线属于等温还是绝热的依据。
先推绝热方程。设 P V = N k T PV=NkT P V = N k T 、U = C V T + c o n s t U=C_VT+\mathrm{const} U = C V T + const 且 C V C_V C V 为常数,绝热时 d Q = 0 \dbar Q=0 d Q = 0 ,所以
C V d T = − P d V \begin{align}
C_V\dd T=-P\dd V
\end{align} C V d T = − P d V
对状态方程微分得 P d V + V d P = N k d T P\dd V+V\dd P=Nk\dd T P d V + V d P = N k d T ,消去 d T \dd T d T ,
( 1 + N k C V ) P d V + V d P = 0 ⟹ V γ P = c o n s t , γ = 1 + N k C V \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} ( 1 + C V N k ) P d V + V d P = 0 ⟹ V γ P = const , γ = 1 + C V N k
把 P = N k T / V P=NkT/V P = N k T / V 代进去还可以写成 V γ − 1 T = c o n s t V^{\gamma-1}T=\mathrm{const} V γ − 1 T = const 。这一步只能对理想气体做,因为消 d T \dd T d T 用了 U = U ( T ) U=U(T) U = U ( T ) 。
两条绝热线给出 V b γ − 1 T 2 = V c γ − 1 T 1 V_b^{\gamma-1}T_2=V_c^{\gamma-1}T_1 V b γ − 1 T 2 = V c γ − 1 T 1 和 V d γ − 1 T 1 = V a γ − 1 T 2 V_d^{\gamma-1}T_1=V_a^{\gamma-1}T_2 V d γ − 1 T 1 = V a γ − 1 T 2 ,两式相除得 V b / V a = V c / V d V_b/V_a=V_c/V_d V b / V a = V c / V d ,这个比值后面要用。
两个等温段的热量,等温时 Δ U = 0 \Delta U=0 Δ U = 0 ,所以 Q = ∫ P d V Q=\int P\dd V Q = ∫ P d V ,
Q 2 = N k T 2 ln V b V a , Q 1 = N k T 1 ln V c V d \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} Q 2 = N k T 2 ln V a V b , Q 1 = N k T 1 ln V d V c
代入效率定义,对数部分的比值相同,一起约掉,
η = 1 − Q 1 Q 2 = 1 − T 1 ln ( V c / V d ) T 2 ln ( V b / V a ) = 1 − T 1 T 2 \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} η = 1 − Q 2 Q 1 = 1 − T 2 ln ( V b / V a ) T 1 ln ( V c / V d ) = 1 − T 2 T 1
第二段证一般热机超不过卡诺热机。取一台卡诺热机 C C C 和一台任意热机 X X X ,都工作在 T 2 > T 1 T_2>T_1 T 2 > T 1 之间,记 W = Q 2 − Q 1 W=Q_2-Q_1 W = Q 2 − Q 1 、W ′ = Q 2 ′ − Q 1 ′ W'=Q_2'-Q_1' W ′ = Q 2 ′ − Q 1 ′ 。选正整数 N , N ′ N,N' N , N ′ 使 Q 2 / Q 2 ′ = N ′ / N Q_2/Q_2'=N'/N Q 2 / Q 2 ′ = N ′ / N ,整数取得足够大就能逼近任何比值。
让 C C C 反向运行 N N N 个循环,X X X 正向运行 N ′ N' N ′ 个循环。反向运行的意思是 W , Q 1 , Q 2 W,Q_1,Q_2 W , Q 1 , Q 2 的符号同时改变,这一步是卡诺热机可逆才做得到的。净效果是
W t o t a l = N ′ W ′ − N W ( Q 2 ) t o t a l = N ′ Q 2 ′ − N Q 2 = 0 ( Q 1 ) t o t a l = N ′ Q 1 ′ − N Q 1 \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} W total ( Q 2 ) total ( Q 1 ) total = N ′ W ′ − N W = N ′ Q 2 ′ − N Q 2 = 0 = N ′ Q 1 ′ − N Q 1
高温热库的净吸热为零,于是整个装置只跟低温热库交换热。一圈下来 Δ U = 0 \Delta U=0 Δ U = 0 ,由第一定律 W t o t a l = − ( Q 1 ) t o t a l W_{\mathrm{total}}=-(Q_1)_{\mathrm{total}} W total = − ( Q 1 ) total 。如果 W t o t a l > 0 W_{\mathrm{total}}>0 W total > 0 ,就是从单一热库取热全部变功,违反 Kelvin 表述,所以 W t o t a l ≤ 0 W_{\mathrm{total}}\le0 W total ≤ 0 ,即 ( Q 1 ) t o t a l ≥ 0 (Q_1)_{\mathrm{total}}\ge0 ( Q 1 ) total ≥ 0 ,也就是 N ′ Q 1 ′ ≥ N Q 1 N'Q_1'\ge NQ_1 N ′ Q 1 ′ ≥ N Q 1 。两边除以 N ′ Q 2 ′ = N Q 2 N'Q_2'=NQ_2 N ′ Q 2 ′ = N Q 2 ,
Q 1 Q 2 ≤ Q 1 ′ Q 2 ′ ⟹ 1 − Q 1 Q 2 ≥ 1 − Q 1 ′ Q 2 ′ ⟹ η ≥ η ′ \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} Q 2 Q 1 ≤ Q 2 ′ Q 1 ′ ⟹ 1 − Q 2 Q 1 ≥ 1 − Q 2 ′ Q 1 ′ ⟹ η ≥ η ′
让 X X X 也取卡诺热机,把两者的角色对调再证一遍,会得到反向的不等式,两条合起来取等号,于是所有卡诺热机效率相同,与工作物质和功的大小都无关。这个”与工作物质无关”正是绝对温标的立足点。定义 T 1 / T 2 ≡ 1 − η T_1/T_2\equiv1-\eta T 1 / T 2 ≡ 1 − η ,效率只由两个温度的比值决定;再固定每一级卡诺循环的功为 1 K 1\ \mathrm{K} 1 K 对应的量,单位也定死了。这样定义的 T T T 与理想气体温标 P V = N k T PV=NkT P V = N k T 给出的 T T T 一致。
熵是态函数
要证的是 ∫ A B d Q / T \int_A^B\dbar Q/T ∫ A B d Q / T 在可逆路径上与路径无关。出发点是一条任意循环的 Clausius 不等式
∮ d Q T ≤ 0 \begin{align}
\oint\frac{\dbar Q}{T}\le0
\end{align} ∮ T d Q ≤ 0
可逆取等号。做法是用等温线和绝热线把循环切成很多小卡诺循环。每个小卡诺循环满足 Q j / T j + Q j + 1 / T j + 1 = 0 Q_j/T_j+Q_{j+1}/T_{j+1}=0 Q j / T j + Q j + 1 / T j + 1 = 0 ,把所有小循环加起来,落在内部的分界线被相邻两个小循环各走一次、方向相反,互相抵消,只剩外边界,于是 ∮ r e v d Q / T = 0 \oint_{\mathrm{rev}}\dbar Q/T=0 ∮ rev d Q / T = 0 。如果循环不可逆,同一步骤给出严格小于零,因为不可逆过程是向平衡态弛豫,热失去了做功的能力。
有了 ∮ r e v d Q / T = 0 \oint_{\mathrm{rev}}\dbar Q/T=0 ∮ rev d Q / T = 0 ,取两条从 A A A 到 B B B 的可逆路径 C 1 C_1 C 1 、C 2 C_2 C 2 ,让 C 1 C_1 C 1 接上 C 2 C_2 C 2 的反向,构成一个可逆循环,于是
∫ C 1 d Q T − ∫ C 2 d Q T = 0 ⟹ ∫ C 1 d Q T = ∫ C 2 d Q T \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} ∫ C 1 T d Q − ∫ C 2 T d Q = 0 ⟹ ∫ C 1 T d Q = ∫ C 2 T d Q
积分与路径无关,就可以拿它定义态函数 S S S ,S ( B ) − S ( A ) = ∫ A B d Q / T S(B)-S(A)=\int_A^B\dbar Q/T S ( B ) − S ( A ) = ∫ A B d Q / T ,等价说法是 d S = d Q / T \dd S=\dbar Q/T d S = d Q / T 恰当,1 / T 1/T 1/ T 是积分因子。
两个直接推论。对任意变换,构造一条从 B B B 回到 A A A 的可逆路径,与原来的任意路径拼成循环,Clausius 不等式给出 S ( B ) − S ( A ) ≥ ∫ A B d Q / T S(B)-S(A)\ge\int_A^B\dbar Q/T S ( B ) − S ( A ) ≥ ∫ A B d Q / T ;绝热时右边为零,于是 Δ S ≥ 0 \Delta S\ge0 Δ S ≥ 0 ,等号只在可逆时成立。
四个 Maxwell 关系
来源只有一句话,四个势的全微分都是恰当的,所以交叉偏导相等。以 d A = − S d T − P d V \dd A=-S\dd T-P\dd V d A = − S d T − P d V 为例,把 − S -S − S 看成 T T T 方向的系数、− P -P − P 看成 V V V 方向的系数,
( ∂ 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} ( ∂ V ∂ S ) T = ( ∂ T ∂ P ) V
另外三个照同样的动作做,只是换势、换变量对,
( ∂ T ∂ V ) S = − ( ∂ P ∂ S ) V 由 d U = T d S − P d V ( ∂ S ∂ V ) T = ( ∂ P ∂ T ) V 由 d A = − S d T − P d V ( ∂ T ∂ P ) S = ( ∂ V ∂ S ) P 由 d H = T d S + V d P − ( ∂ S ∂ P ) T = ( ∂ V ∂ T ) P 由 d G = − S d T + V d P \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} ( ∂ V ∂ T ) S ( ∂ V ∂ S ) T ( ∂ P ∂ T ) S − ( ∂ P ∂ S ) T = − ( ∂ S ∂ P ) V = ( ∂ T ∂ P ) V = ( ∂ S ∂ V ) P = ( ∂ T ∂ V ) P 由 d U = T d S − P d V 由 d A = − S d T − P d V 由 d H = T d S + V d P 由 d G = − S d T + V d P
负号的来源是换变量时符号会翻。d U \dd U d U 里 − P d V -P\dd V − P d V 带着负号,交叉偏导就多一个负号;d G \dd G d G 里 − S d T -S\dd T − S d T 带着负号,负号跑到左边。所以背这四条不如每次现场推,两行就能写完,还不会记错符号。
从 S = k log Γ S=k\log\Gamma S = k log Γ 到温度与压强
第一步证 S S S 广延。把系统切成两个固定 N 1 , V 1 , N 2 , V 2 N_1,V_1,N_2,V_2 N 1 , V 1 , N 2 , V 2 的子系统,H = H 1 + H 2 H=H_1+H_2 H = H 1 + H 2 ,忽略相互作用但保留交换能量的能力。整个系统的壳层体积是
Γ ( E ) = ∫ E 1 < H 1 < E 1 + Δ , E 2 < H 2 < E 2 + Δ d p 1 d q 1 d p 2 d q 2 = Γ 1 ( E 1 ) Γ 2 ( E 2 ) \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 ) = ∫ E 1 < H 1 < E 1 + Δ , E 2 < H 2 < E 2 + Δ d p 1 d q 1 d p 2 d q 2 = Γ 1 ( E 1 ) Γ 2 ( E 2 )
麻烦在于 E = E 1 + E 2 E=E_1+E_2 E = E 1 + E 2 的分解不唯一,E 1 E_1 E 1 可以从 0 0 0 变到 E E E 。把 E E E 按 Δ \Delta Δ 切成 E / Δ E/\Delta E /Δ 段,
Γ ( E ) = ∑ i = 0 E / Δ Γ 1 ( E i ) Γ 2 ( E − E i ) \begin{align}
\Gamma(E)=\sum_{i=0}^{E/\Delta}\Gamma_1(E_i)\Gamma_2(E-E_i)
\end{align} Γ ( E ) = i = 0 ∑ E /Δ Γ 1 ( E i ) Γ 2 ( E − E i )
关键是求和里只有一项占主导。设最大的那一项对应 E ˉ 1 + E ˉ 2 = E \bar E_1+\bar E_2=E E ˉ 1 + 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 Γ 1 ( E ˉ 1 ) Γ 2 ( E ˉ 2 ) ≤ Γ ( E ) ≤ Γ 1 ( E ˉ 1 ) Γ 2 ( E ˉ 2 ) ⋅ E /Δ ,取对数后上下界只差 k log ( E / Δ ) k\log(E/\Delta) k log ( E /Δ ) ,而 S ∝ N S\propto N S ∝ N ,这一项是 O ( log N ) O(\log N) O ( log N ) ,可以丢掉。于是
S ( E , V ) = S 1 ( E ˉ 1 , V 1 ) + S 2 ( E ˉ 2 , V 2 ) + 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} S ( E , V ) = S 1 ( E ˉ 1 , V 1 ) + S 2 ( E ˉ 2 , V 2 ) + O ( log N )
主导项的位置由 δ log [ Γ 1 Γ 2 ] = 0 \delta\log[\Gamma_1\Gamma_2]=0 δ log [ Γ 1 Γ 2 ] = 0 配上 δ E 1 + δ E 2 = 0 \delta E_1+\delta E_2=0 δ E 1 + δ E 2 = 0 定出,即
∂ S 1 ∂ E 1 ∣ E ˉ 1 = ∂ S 2 ∂ E 2 ∣ 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} ∂ E 1 ∂ S 1 E ˉ 1 = ∂ E 2 ∂ S 2 E ˉ 2
既然两边相等,就把它定义成 1 / T 1/T 1/ T ,
( ∂ S ∂ E ) V ≡ 1 T \begin{align}
\left(\frac{\partial S}{\partial E}\right)_V\equiv\frac{1}{T}
\end{align} ( ∂ E ∂ S ) V ≡ T 1
这一步同时给出了第零定律 T 1 = T 2 T_1=T_2 T 1 = T 2 。第二步定压强。有两个独立参数时,
d S = ( ∂ S ∂ E ) V d E + ( ∂ S ∂ V ) E d V = 1 T d E + ( ∂ S ∂ V ) E d V \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} d S = ( ∂ E ∂ S ) V d E + ( ∂ V ∂ S ) E d V = T 1 d E + ( ∂ V ∂ S ) E d V
定义 P ≡ T ( ∂ S / ∂ V ) E P\equiv T(\partial S/\partial V)_E P ≡ T ( ∂ S / ∂ V ) E ,上式就变成 d E = T d S − P d V \dd E=T\dd S-P\dd V d E = T d S − P d V 。要检查这个 P P P 是不是真的压强,看两条。它得是强度量,S S S 和 V V V 都广延,偏导是强度量,乘 T T T 还是强度量。它还得和第一章的压强对得上,在 S S S 固定时 d E = − P d V \dd E=-P\dd V d E = − P d V ,所以 P = − ( ∂ E / ∂ V ) S P=-(\partial E/\partial V)_S P = − ( ∂ E / ∂ V ) S ,正是第一章那个压强。
还要说明 S S S 满足第二定律。用 Σ ( E ) \Sigma(E) Σ ( E ) 的写法,孤立系统 N N N 与 E E E 不变、V V V 只增不减,而 Σ ( E ) \Sigma(E) Σ ( E ) 是 V V V 的非减函数(体积变大,H < E H<E H < E 的区域只会变大),所以 S S S 不减。
能量均分定理
要算的是 ⟨ x i ∂ H / ∂ x j ⟩ \langle x_i\partial H/\partial x_j\rangle ⟨ x i ∂ H / ∂ x j ⟩ 在微正则系综下的值,x i x_i x i 取 p i p_i p i 或 q i q_i q i 。按定义
⟨ x i ∂ H ∂ x j ⟩ = 1 Γ ( E ) ∫ E < H < E + Δ d p d q x i ∂ H ∂ x j \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} ⟨ x i ∂ x j ∂ H ⟩ = Γ ( E ) 1 ∫ E < H < E + Δ d p d q x i ∂ x j ∂ H
把壳层积分写成两个累计积分之差,再写成对 E E E 的导数,
⟨ x i ∂ H ∂ x j ⟩ = Δ Γ ( E ) ∂ ∂ E ∫ H < E d p d q x i ∂ H ∂ x j \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} ⟨ x i ∂ x j ∂ H ⟩ = Γ ( E ) Δ ∂ E ∂ ∫ H < E d p d q x i ∂ x j ∂ H
里面的积分做分部积分。注意 ∂ H / ∂ x j = ∂ ( H − E ) / ∂ x j \partial H/\partial x_j=\partial(H-E)/\partial x_j ∂ H / ∂ x j = ∂ ( H − E ) / ∂ x j ,
∫ H < E d p d q x i ∂ ( H − E ) ∂ x j = ∫ H < E d p d q ∂ ∂ x j [ x i ( H − E ) ] − δ i j ∫ H < E d p d q ( 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 < E d p d q x i ∂ x j ∂ ( H − E ) = ∫ H < E d p d q ∂ x j ∂ [ x i ( H − E ) ] − δ ij ∫ H < E d p d q ( H − E )
第一项是全导数,化成 H = E H=E H = E 曲面上的面积分,而在那个曲面上 H − E = 0 H-E=0 H − E = 0 ,所以整项为零。剩下
⟨ x i ∂ H ∂ x j ⟩ = δ i j ω ( E ) ∂ ∂ E ∫ H < E d p d q ( E − H ) = δ i j ω ( 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} ⟨ x i ∂ x j ∂ H ⟩ = ω ( E ) δ ij ∂ E ∂ ∫ H < E d p d q ( E − H ) = ω ( E ) δ ij Σ ( E )
这里用了 ∂ ( E − H ) / ∂ E = 1 \partial(E-H)/\partial E=1 ∂ ( E − H ) / ∂ E = 1 。最后把 Σ / ω \Sigma/\omega Σ/ ω 换成温度,因为 ω = ∂ Σ / ∂ E \omega=\partial\Sigma/\partial E ω = ∂ Σ/ ∂ E ,
Σ ω = ( ∂ log Σ ∂ E ) − 1 = k ( ∂ S ∂ E ) − 1 = k T \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} ω Σ = ( ∂ E ∂ log Σ ) − 1 = k ( ∂ E ∂ S ) − 1 = k T
于是 ⟨ x i ∂ H / ∂ x j ⟩ = δ i j k T \langle x_i\partial H/\partial x_j\rangle=\delta_{ij}kT ⟨ x i ∂ H / ∂ x j ⟩ = δ ij k T 。这一条对任何 Hamiltonian 都成立,只要边界项能丢掉。
要得到 ⟨ H ⟩ = f k T / 2 \langle H\rangle=fkT/2 ⟨ H ⟩ = f k T /2 还得加一个条件,H H H 必须是二次型。设 H = ∑ i A i P i 2 + ∑ i B i Q i 2 H=\sum_iA_iP_i^2+\sum_iB_iQ_i^2 H = ∑ i A i P i 2 + ∑ i B i Q i 2 ,由齐次函数的 Euler 定理
∑ i ( P i ∂ H ∂ P i + Q i ∂ H ∂ Q i ) = 2 H \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} i ∑ ( P i ∂ P i ∂ H + Q i ∂ Q i ∂ H ) = 2 H
两边取系综平均,左边每个自由度贡献 k T kT k T ,共 f f f 个,于是 f k T = 2 ⟨ H ⟩ fkT=2\langle H\rangle f k T = 2 ⟨ H ⟩ ,即 ⟨ H ⟩ = f k T / 2 \langle H\rangle=fkT/2 ⟨ H ⟩ = f k T /2 。单原子理想气体 f = 3 f=3 f = 3 ,U = 3 N k T / 2 U=3NkT/2 U = 3 N k T /2 ;谐振子 H = p 2 / 2 m + ω 2 q 2 / 2 H=p^2/2m+\omega^2q^2/2 H = p 2 /2 m + ω 2 q 2 /2 的两个变量都是二次的,f = 2 f=2 f = 2 ,⟨ H ⟩ = k T \langle H\rangle=kT ⟨ H ⟩ = k T 。
从微正则到 MB 分布
把 μ \mu μ 空间切成 K K K 个小格,第 l l l 格的体积记为 g l g_l g l (最后可以取 1 1 1 ),粒子数 n l n_l n l ,约束是 ∑ n l = N \sum n_l=N ∑ n l = N 与 ∑ ϵ l n l = E \sum\epsilon_ln_l=E ∑ ϵ l n l = E 。给定 { n l } \{n_l\} { n l } 之后,Γ \Gamma Γ 空间里对应的体积正比于把 N N N 个分子分到这 K K K 组里的分法数,
Ω ( { n l } ) ∝ N ! ∏ l n l ! ∏ l g l n l \begin{align}
\Omega(\{n_l\})\propto\frac{N!}{\prod_ln_l!}\prod_lg_l^{n_l}
\end{align} Ω ({ n l }) ∝ ∏ l n l ! N ! l ∏ g l n l
N ! N! N ! 来自分子的排法,n l ! n_l! n l ! 除掉同一格内互换(同一格里分子对 f ( p , q ) f(p,q) f ( p , q ) 的贡献相同,互换不算新状态)。分子数很大,用 Stirling 公式 log n ! ≃ n ( log n − 1 ) \log n!\simeq n(\log n-1) log n ! ≃ n ( log n − 1 ) ,
log Ω ( { n l } ) = N ( log N − 1 ) − ∑ l n l ( log n l − 1 ) + ∑ l n l log g l \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} log Ω ({ n l }) = N ( log N − 1 ) − l ∑ n l ( log n l − 1 ) + l ∑ n l log g l
在两条约束下求极值,引入乘子 α \alpha α 、β \beta β ,
δ [ log Ω − α ∑ l n l − β ∑ l ϵ l n l ] = 0 \begin{align}
\delta\left[\log\Omega-\alpha\sum_ln_l-\beta\sum_l\epsilon_ln_l\right]=0
\end{align} δ [ log Ω − α l ∑ n l − β l ∑ ϵ l n l ] = 0
把 δ n l \delta n_l δ n l 看成独立变分,每个系数都得为零,
− log n ˉ l + log g l − α − β ϵ l = 0 ⟹ n ˉ l = g l e − α − β ϵ 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} − log n ˉ l + log g l − α − β ϵ l = 0 ⟹ n ˉ l = g l e − α − β ϵ l
二阶变分是 − ∑ l ( δ n l ) 2 / n ˉ l < 0 -\sum_l(\delta n_l)^2/\bar n_l<0 − ∑ l ( δ n l ) 2 / n ˉ l < 0 ,所以这个驻点确实是极大。再用 ∑ n ˉ l = N \sum\bar n_l=N ∑ n ˉ l = N 和 ∑ ϵ l n ˉ l = E \sum\epsilon_l\bar n_l=E ∑ ϵ l n ˉ l = E 定 α \alpha α 与 β \beta β 。写成配分函数就是 N = e − α Z N=e^{-\alpha}Z N = e − α Z ,α = log ( Z / N ) \alpha=\log(Z/N) α = log ( Z / N ) 。
β \beta β 的身份靠熵来定。把 n ˉ l \bar n_l n ˉ l 代回 log Ω \log\Omega log Ω ,用 N = e − α Z N=e^{-\alpha}Z N = e − α Z 化掉 α \alpha α ,得到
S = N k ( log Z − β ∂ log Z ∂ β ) + c o n s t \begin{align}
S=Nk\left(\log Z-\beta\frac{\partial\log Z}{\partial\beta}\right)+\mathrm{const}
\end{align} S = N k ( log Z − β ∂ β ∂ log Z ) + const
固定外参量 y y y 对它微分,log Z \log Z log Z 的两项正好抵消,只剩
d S = N k d ( log Z − β ∂ log Z ∂ β ) = k β d E \begin{align}
\dd S=Nk\,\dd\!\left(\log Z-\beta\frac{\partial\log Z}{\partial\beta}\right)
=k\beta\,\dd E
\end{align} d S = N k d ( log Z − β ∂ β ∂ log Z ) = k β d E
要求 d Q = T d S \dbar Q=T\dd S d Q = T d S ,也就是 d S = d E / T \dd S=\dd E/T d S = d E / T (y y y 固定时不做功),比较即得 β = 1 / k T \beta=1/kT β = 1/ k T 。到这一步热力学第二定律在统计力学里被推出来了,前提是第一定律仍然是假设。
典型题型
算熵变
题面给初末态和一条路径,问 Δ S \Delta S Δ S ,有时顺带问环境的熵变和”损失了多少可用能”。判断第一步是这个过程可不可逆。可逆就直接沿路径积 ∫ d Q / T \int\dbar Q/T ∫ d Q / T ;不可逆就不能用路径上的 d Q \dbar Q d Q ,要另设计一条连接同样初末态的可逆路径,因为 S S S 是态函数,走哪条路都行。常用的是先等温后等容,或者直接用 d S = ( C V / T ) d T + ( ∂ P / ∂ T ) V d V \dd S=(C_V/T)\dd T+(\partial P/\partial T)_V\dd V d S = ( C V / T ) d T + ( ∂ P / ∂ T ) V d V 沿一条方便的折线积分。
例,一摩尔理想气体在 300 K 300\ \mathrm{K} 300 K 下自由膨胀到两倍体积,求熵变。自由膨胀 Δ W = 0 \Delta W=0 Δ W = 0 ,绝热 Δ Q = 0 \Delta Q=0 Δ Q = 0 ,所以由第一定律 Δ U = 0 \Delta U=0 Δ U = 0 ,理想气体 U = U ( T ) U=U(T) U = U ( T ) 给出 Δ T = 0 \Delta T=0 Δ T = 0 。实际的路径上 d Q ≡ 0 \dbar Q\equiv0 d Q ≡ 0 ,直接积会得零,那是错的。换一条可逆等温膨胀路径,温度不变,Δ U = 0 \Delta U=0 Δ U = 0 ,气体对外做的功全部来自吸热,
Δ S = ∫ d Q T = 1 T ∫ P d V = N k ln V 2 V 1 = N k ln 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} Δ S = ∫ T d Q = T 1 ∫ P d V = N k ln V 1 V 2 = N k ln 2
取 N k = R ≃ 8.3 J / K Nk=R\simeq8.3\ \mathrm{J/K} N k = R ≃ 8.3 J/K ,得 Δ S ≃ 5.8 J / K \Delta S\simeq5.8\ \mathrm{J/K} Δ S ≃ 5.8 J/K 。环境没变,Δ S e n v = 0 \Delta S_{\mathrm{env}}=0 Δ S env = 0 ,总熵增加了 5.8 J / K 5.8\ \mathrm{J/K} 5.8 J/K ,符合不可逆过程的 Δ S > ∫ d Q / T = 0 \Delta S>\int\dbar Q/T=0 Δ S > ∫ d Q / T = 0 。这 T Δ S T\Delta S T Δ S 就是那句”损失了可用的能量”,同一份熵增如果在可逆过程里,还能换成 5.8 × 300 ≈ 1.7 k J 5.8\times300\approx1.7\ \mathrm{kJ} 5.8 × 300 ≈ 1.7 kJ 的功。
循环效率
题面给 p p p -V V V 图上的一圈,问效率、功或某一支的吸放热。一圈下来 Δ U = 0 \Delta U=0 Δ U = 0 ,所以 W = Q n e t W=Q_{\mathrm{net}} W = Q net ,而 W W W 就是回路围出的面积(顺时针为正)。效率的分母只是吸热的那几段之和,不要把放热段混进去,这是最常见的失分点。若回路由等温线和绝热线组成,先认卡诺,直接用 1 − T 1 / T 2 1-T_1/T_2 1 − T 1 / T 2 ;若不是卡诺,就逐段算 Q Q Q 。
例,工作在 400 K 400\ \mathrm{K} 400 K 与 300 K 300\ \mathrm{K} 300 K 之间的卡诺热机,每循环从高温热库吸热 800 J 800\ \mathrm{J} 800 J 。效率 η = 1 − 300 / 400 = 0.25 \eta=1-300/400=0.25 η = 1 − 300/400 = 0.25 ,做功 W = η Q 2 = 200 J W=\eta Q_2=200\ \mathrm{J} W = η Q 2 = 200 J ,向低温热库放热 Q 1 = Q 2 − W = 600 J Q_1=Q_2-W=600\ \mathrm{J} Q 1 = Q 2 − W = 600 J 。核对 Q 2 / T 2 = 2 J / K Q_2/T_2=2\ \mathrm{J/K} Q 2 / T 2 = 2 J/K 、Q 1 / T 1 = 2 J / K Q_1/T_1=2\ \mathrm{J/K} Q 1 / T 1 = 2 J/K ,比值为一,可逆循环的判据成立。同样两个热库之间的任何其他热机,效率都不超过 0.25 0.25 0.25 ;想提高效率只能抬高 T 2 T_2 T 2 或压低 T 1 T_1 T 1 ,这是卡诺定理的实用含义。
Maxwell 关系套用
题面要一个不好测的偏导,比如 ( ∂ S / ∂ P ) T (\partial S/\partial P)_T ( ∂ S / ∂ P ) T 、( ∂ U / ∂ V ) T (\partial U/\partial V)_T ( ∂ U / ∂ V ) T 、( ∂ H / ∂ P ) T (\partial H/\partial P)_T ( ∂ H / ∂ P ) T 。做法是先认它出现在哪个势的全微分里,用 Maxwell 关系换成可测量(( ∂ V / ∂ T ) P (\partial V/\partial T)_P ( ∂ V / ∂ T ) P 、( ∂ P / ∂ T ) V (\partial P/\partial T)_V ( ∂ P / ∂ T ) V 、C P C_P C P 这类),再把状态方程代进去。不要凭记忆写关系式,现场从 d A \dd A d A 或 d G \dd G d G 推,两行的事。
例,求 van der Waals 气体的 ( ∂ U / ∂ V ) T (\partial U/\partial V)_T ( ∂ U / ∂ 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} ( ∂ V ∂ U ) T = T ( ∂ T ∂ P ) V − P
对 ( P + a / V 2 ) ( V − b ) = N k T (P+a/V^2)(V-b)=NkT ( P + a / V 2 ) ( V − b ) = N k T 求偏导,a a a 、b b b 与 T T T 无关,
( ∂ P ∂ T ) V = N k V − b \begin{align}
\left(\frac{\partial P}{\partial T}\right)_V=\frac{Nk}{V-b}
\end{align} ( ∂ T ∂ P ) V = V − b N k
代回去,
( ∂ U ∂ V ) T = N k T V − b − ( N k T V − b − a V 2 ) = a V 2 \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} ( ∂ V ∂ U ) T = V − b N k T − ( V − b N k T − V 2 a ) = V 2 a
结果是正的,说明等温膨胀时内能增加,因为要克服分子间的吸引做功。a = 0 a=0 a = 0 时回到理想气体的 ( ∂ U / ∂ V ) T = 0 (\partial U/\partial V)_T=0 ( ∂ U / ∂ V ) T = 0 ,也就回到 U = U ( T ) U=U(T) U = U ( T ) 。这道题还顺带给出一个结论,vdW 气体的内能可以写成 U = C V T − a / V + c o n s t U=C_VT-a/V+\mathrm{const} U = C V T − a / V + const ,把它对 T T T 求偏导得 C V C_V C V ,与体积无关。
相变与潜热
题面给潜热、比容差和温度,问相界线的斜率,或者反过来问沸点随压强怎么变。核心只有一条,g 1 = g 2 g_1=g_2 g 1 = g 2 。沿相界线两边都满足这个等式,把 d g = − s d T + v d P \dd g=-s\dd T+v\dd P d g = − s d T + v d P 对两相各写一遍再相减,沿同一条相界线走时两边的 d T \dd T d T 与 d P \dd P d P 相同,于是 Δ s d T = Δ v d P \Delta s\dd T=\Delta v\dd P Δ s d T = Δ v d P ,整理即得 Clapeyron 方程。用 l = T Δ s l=T\Delta s l = T Δ s 还可以写成 l / ( T Δ v ) l/(T\Delta v) l / ( T Δ v ) 。
例,取潜热的量级 l ∼ 2 × 10 6 J / k g l\sim2\times10^{6}\ \mathrm{J/kg} l ∼ 2 × 1 0 6 J/kg 、比容差的量级 Δ v ∼ 1 m 3 / k g \Delta v\sim1\ \mathrm{m^3/kg} Δ v ∼ 1 m 3 /kg 、T = 373 K T=373\ \mathrm{K} T = 373 K ,则
d P d T = l T Δ v ∼ 2 × 10 6 373 × 1 ≈ 5 × 10 3 P a / 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} d T d P = T Δ v l ∼ 373 × 1 2 × 1 0 6 ≈ 5 × 1 0 3 Pa/K
数量级是每开尔文几 kPa,即每升高 1 K 1\ \mathrm{K} 1 K 饱和蒸气压上升几千帕。这个量级说明为什么高压锅里沸点能显著抬高。符号也值得记,气相的熵与比容都比液相大,Δ s > 0 \Delta s>0 Δ s > 0 、Δ v > 0 \Delta v>0 Δ v > 0 ,所以 d P / d T > 0 \dd P/\dd T>0 d P / d T > 0 ,沸点随压强上升。若某处 Δ v \Delta v Δ v 为负(凝固时体积变大),相界线的斜率就反过来,固液线几乎竖直就是 Δ v \Delta v Δ v 很小的结果。
理想气体配分函数
题面给 ϵ ( p ) \epsilon(p) ϵ ( p ) 或 H H H ,要求 E E E 、P P P 、C V C_V C V 、S S S 。路线固定,先算 Z Z Z ,再按表里的公式求导。动量积分是高斯积分,∫ d 3 p e − β p 2 / 2 m = ( 2 π m / β ) 3 / 2 \int\dd^3p\,e^{-\beta p^2/2m}=(2\pi m/\beta)^{3/2} ∫ d 3 p e − β p 2 /2 m = ( 2 π m / β ) 3/2 。
例,ϵ = p 2 / 2 m \epsilon=p^2/2m ϵ = p 2 /2 m ,无外势。位置积分直接给出体积 V V V ,
Z = ∫ d 3 p d 3 q e − β p 2 / 2 m = V ( 2 π m β ) 3 / 2 = V ( 2 π m k T ) 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} Z = ∫ d 3 p d 3 q e − β p 2 /2 m = V ( β 2 π m ) 3/2 = V ( 2 π mk T ) 3/2
取对数 log Z = log V + 3 2 log ( 2 π m k T ) \log Z=\log V+\frac32\log(2\pi mkT) log Z = log V + 2 3 log ( 2 π mk T ) ,再求导,
E = − N ∂ log Z ∂ β = 3 N 2 β = 3 2 N k T P = N β ∂ log Z ∂ V = N k T V \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} E P = − N ∂ β ∂ log Z = 2 β 3 N = 2 3 N k T = β N ∂ V ∂ log Z = V N k T
C V = ( ∂ E / ∂ T ) V = 3 N k / 2 C_V=(\partial E/\partial T)_V=3Nk/2 C V = ( ∂ E / ∂ T ) V = 3 N k /2 。这三个结果和微正则路线算出来的一致,可以互相验算。注意求 P P P 时 Z Z Z 里的 V V V 只出现在一个因子上,如果 ϵ \epsilon ϵ 本身依赖 V V V (比如把 q q q 的积分限放大),就得把那个依赖也算进去,讲义在这一处专门提了一句。
能量均分应用
题面给一个 H H H ,问 ⟨ H ⟩ \langle H\rangle ⟨ H ⟩ 、C V C_V C V 或 γ \gamma γ 。做法是数二次型自由度的个数 f f f ,然后 ⟨ H ⟩ = f k T / 2 \langle H\rangle=fkT/2 ⟨ H ⟩ = f k T /2 、C V = f k / 2 C_V=fk/2 C V = f k /2 、γ = 1 + 2 / f \gamma=1+2/f γ = 1 + 2/ f 。数的时候只数二次项,一次项、四次项、1 / r 1/r 1/ r 项都不算。
例,谐振子 H = p 2 / 2 m + ω 2 q 2 / 2 H=p^2/2m+\omega^2q^2/2 H = p 2 /2 m + ω 2 q 2 /2 。两个变量 p p p 、q q q 都以二次方出现,f = 2 f=2 f = 2 ,所以 ⟨ H ⟩ = 2 ⋅ k T / 2 = k T \langle H\rangle=2\cdot kT/2=kT ⟨ H ⟩ = 2 ⋅ k T /2 = k T ,C V = N k C_V=Nk C V = N k ,γ = 1 + 2 / 2 = 2 \gamma=1+2/2=2 γ = 1 + 2/2 = 2 。这个结果也可以直接验证,p p p 的积分和 q q q 的积分各是高斯型,各自贡献 k T / 2 kT/2 k T /2 。单原子理想气体只有三个平动动量,f = 3 f=3 f = 3 ,C V = 3 N k / 2 C_V=3Nk/2 C V = 3 N k /2 ,γ = 5 / 3 \gamma=5/3 γ = 5/3 。反过来,如果某个坐标只以 q 4 q^4 q 4 的形式出现,那一项不给 k T / 2 kT/2 k T /2 ,这个自由度不能算。
涨落估算
题面问”能不能把两个子系统的能量当成确定值相加”,或者问某个量的相对误差有多大。要算的是峰宽。把 Γ 1 ( E 1 ) Γ 2 ( E − E 1 ) \Gamma_1(E_1)\Gamma_2(E-E_1) Γ 1 ( E 1 ) Γ 2 ( E − E 1 ) 的对数在极值附近展开到二阶,一阶项给 T 1 = T 2 T_1=T_2 T 1 = T 2 ,二阶项给宽度。
log [ Γ 1 ( E 1 ) Γ 2 ( E − E 1 ) ] = log Γ 1 Γ 2 ∣ E ˉ 1 + 1 2 k ( S 1 ′ ′ + S 2 ′ ′ ) ( δ E 1 ) 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} log [ Γ 1 ( E 1 ) Γ 2 ( E − E 1 ) ] = log Γ 1 Γ 2 E ˉ 1 + 2 k 1 ( S 1 ′′ + S 2 ′′ ) ( δ E 1 ) 2 + ⋯
其中 S ′ ′ = ∂ 2 S / ∂ E 2 S''=\partial^2S/\partial E^2 S ′′ = ∂ 2 S / ∂ E 2 。由 ( ∂ S / ∂ E ) V = 1 / T (\partial S/\partial E)_V=1/T ( ∂ S / ∂ E ) V = 1/ T 再求一次导,
∂ 2 S ∂ E 2 = ∂ ∂ E ( 1 T ) = − 1 T 2 ( ∂ T ∂ E ) V = − 1 T 2 C V \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} ∂ E 2 ∂ 2 S = ∂ E ∂ ( T 1 ) = − T 2 1 ( ∂ E ∂ T ) V = − T 2 C V 1
所以 S 1 ′ ′ + S 2 ′ ′ = − ( 1 / T 2 ) ( 1 / C V 1 + 1 / C V 2 ) = − 1 / ( T 2 C V ) S_1''+S_2''=-(1/T^2)(1/C_{V1}+1/C_{V2})=-1/(T^2C_V) S 1 ′′ + S 2 ′′ = − ( 1/ T 2 ) ( 1/ C V 1 + 1/ C V 2 ) = − 1/ ( T 2 C V ) ,用 C V = C V 1 + C V 2 C_V=C_{V1}+C_{V2} C V = C V 1 + C V 2 。分布近似高斯,方差为 − k / ( S 1 ′ ′ + S 2 ′ ′ ) -k/(S_1''+S_2'') − k / ( S 1 ′′ + S 2 ′′ ) ,
⟨ ( δ E 1 ) 2 ⟩ = k T 2 C V \begin{align}
\langle(\delta E_1)^2\rangle=kT^2C_V
\end{align} ⟨( δ E 1 ) 2 ⟩ = k T 2 C V
若子系统 2 远大于 1,C V ≃ C V 1 C_V\simeq C_{V1} C V ≃ C V 1 ,退化成 ⟨ ( δ E 1 ) 2 ⟩ ≃ k T 2 C V 1 \langle(\delta E_1)^2\rangle\simeq kT^2C_{V1} ⟨( δ E 1 ) 2 ⟩ ≃ k T 2 C V 1 ,也就是小系统的能量涨落。相对涨落随 N N N 变小,取理想气体 C V = 3 N k / 2 C_V=3Nk/2 C V = 3 N k /2 、E = 3 N k T / 2 E=3NkT/2 E = 3 N k T /2 ,
δ E E = k T 2 ⋅ 3 N k / 2 3 N k T / 2 = 2 3 N \begin{align}
\frac{\delta E}{E}=\frac{\sqrt{kT^2\cdot3Nk/2}}{3NkT/2}
=\sqrt{\frac{2}{3N}}
\end{align} E δ E = 3 N k T /2 k T 2 ⋅ 3 N k /2 = 3 N 2
N ∼ 10 23 N\sim10^{23} N ∼ 1 0 23 时这个比值是 10 − 12 10^{-12} 1 0 − 12 量级。结论是宏观系统的能量分配几乎确定,所以系综平均可以用最可几值代替,这是整章统计力学能成立的前提。
题型 题面特征 第一动作 算熵变 给初末态,问 Δ S \Delta S Δ S 或总熵 判可逆,不可逆就另设一条可逆路径 循环效率 一圈闭合曲线,问 η \eta η 或 W W W 圈上 Δ U = 0 \Delta U=0 Δ U = 0 ,W W W 取面积,Q 2 Q_2 Q 2 只加吸热段 Maxwell 套用 问一个难测的偏导 找它出现在哪个势的全微分里 相变潜热 给 l l l 、Δ v \Delta v Δ v 、T T T 写 g 1 = g 2 g_1=g_2 g 1 = g 2 ,沿相界线微分 配分函数 给 ϵ ( p ) \epsilon(p) ϵ ( p ) ,问热力学量 先算 Z Z Z ,再对 β \beta β 和 V V V 求导 能量均分 给 H H H ,问 C V C_V C V 或 γ \gamma γ 数二次型自由度 f f f 涨落 问平均值能不能当确定值 把 log [ Γ 1 Γ 2 ] \log[\Gamma_1\Gamma_2] log [ Γ 1 Γ 2 ] 展开到二阶
易错点
先把最容易失分的三块写成散文,剩下的列在表里。
不可逆过程的熵变是第一个大坑。d S = d Q / T \dd S=\dbar Q/T d S = d Q / T 只在可逆时成立,写成 Δ S = ∫ d Q / T \Delta S=\int\dbar Q/T Δ S = ∫ d Q / T 对不可逆过程是错的,因为不可逆时这个积分小于 Δ S \Delta S Δ S 。自由膨胀就是判据,d Q ≡ 0 \dbar Q\equiv0 d Q ≡ 0 然而 Δ S = N k ln 2 \Delta S=Nk\ln2 Δ S = N k ln 2 。办法是换路径,只认初末态。另一个常见漏洞是只算气体,忘了环境,其实要判断的是孤立系(气体加热库)的总熵,自由膨胀时热库没变,总熵才等于气体的熵增。
偏导的固定量是第二个大坑。( ∂ U / ∂ T ) V (\partial U/\partial T)_V ( ∂ U / ∂ T ) V 与 ( ∂ U / ∂ T ) P (\partial U/\partial T)_P ( ∂ U / ∂ T ) P 是两个不同的量,前者是 C V C_V C V ,后者对理想气体等于 C V C_V C V 但一般不等于。同理 ( ∂ S / ∂ V ) T = ( ∂ P / ∂ T ) V (\partial S/\partial V)_T=(\partial P/\partial T)_V ( ∂ S / ∂ V ) T = ( ∂ P / ∂ T ) V 里的下标不能省,省了以后别人没法判断你指的是哪条路径。写偏导时养成每次都写上下标的习惯,考试时省下标扣的分比想不起来公式还多。
C V C_V C V 与 C P C_P C P 何时能互换是第三个。C P − C V = N k C_P-C_V=Nk C P − C V = N k 只对理想气体成立,把它当普遍结论套到 van der Waals 气体上会错。用熵方程时也一样,T d S = C V d T + T ( ∂ P / ∂ T ) V d V T\dd S=C_V\dd T+T(\partial P/\partial T)_V\dd V T d S = C V d T + T ( ∂ P / ∂ T ) V d V 的独立变量是 ( T , V ) (T,V) ( T , V ) ,另一条 T d S = C P d T − T ( ∂ V / ∂ T ) P d P T\dd S=C_P\dd T-T(\partial V/\partial T)_P\dd P T d S = C P d T − T ( ∂ V / ∂ T ) P d P 的独立变量是 ( T , P ) (T,P) ( T , P ) ,两条不能混着用,也不能把一条里的 C V C_V C V 直接换成 C P C_P C P 。
坑 为什么错 正确做法 不可逆过程用 Δ S = ∫ d Q / T \Delta S=\int\dbar Q/T Δ S = ∫ d Q / T 等号只对可逆成立,不可逆时积分偏小 换一条同初末态的可逆路径 只算系统的熵 熵不减是对孤立系说的 系统与热库一起算 偏导不写固定量 不同固定量给出不同的值 每次写全下标 把 C P − C V = N k C_P-C_V=Nk C P − C V = N k 当普遍式 它只在理想气体成立 先看状态方程是不是 P V = N k T PV=NkT P V = N k T Maxwell 关系的负号丢一半 负号来自换变量时的符号翻转 从势的全微分现场推 Γ \Gamma Γ 、Σ \Sigma Σ 、ω \omega ω 混用Γ \Gamma Γ 是壳层体积,Σ \Sigma Σ 是整块体积,ω \omega ω 是态密度算熵时三者可换,算体积时不能换 S = k log Γ S=k\log\Gamma S = k log Γ 直接取对数Γ \Gamma Γ 带 ( d p d q ) N (\dd p\,\dd q)^N ( d p d q ) N 的量纲先规定单位体积,或明确 S S S 只定到差一常数 用 Σ ∝ V N E 3 N / 2 \Sigma\propto V^NE^{3N/2} Σ ∝ V N E 3 N /2 直接得熵 得到 S ∼ N k log V S\sim Nk\log V S ∼ N k log V ,混同种气体时破坏广延 求和中只有一项占主导这一步不能省;同种粒子加 1 / N ! 1/N! 1/ N ! 能量均分乱套 ⟨ x i ∂ H / ∂ x j ⟩ = δ i j k T \langle x_i\partial H/\partial x_j\rangle=\delta_{ij}kT ⟨ x i ∂ H / ∂ x j ⟩ = δ ij k T 是一般的,⟨ H ⟩ = f k T / 2 \langle H\rangle=fkT/2 ⟨ H ⟩ = f k T /2 只对二次型只数二次型自由度 涨落公式里放 C P C_P C P 微正则下 V V V 固定,S ′ ′ = − 1 / ( T 2 C V ) S''=-1/(T^2C_V) S ′′ = − 1/ ( T 2 C V ) 只能用 C V C_V C V 判断违反 Kelvin 时漏看其他变化 等温膨胀从单一热库吸热全变功,但体积变了 先检查”唯一效果”是否成立 相变用 A A A 最小 相变在等温等压下进行 用 G G G 最小,得 g 1 = g 2 g_1=g_2 g 1 = g 2 把 vdW 等温线的上升段当相变平台 那段 ∂ P / ∂ V > 0 \partial P/\partial V>0 ∂ P / ∂ V > 0 ,压缩率为负 用 Maxwell 等面积构造画水平线 E = − N ∂ log Z / ∂ β E=-N\partial\log Z/\partial\beta E = − N ∂ log Z / ∂ β 少负号Z Z Z 随 β \beta β 减小,log Z \log Z log Z 对 β \beta β 的导数为负记住 E E E 必须为正 把 Δ \Delta Δ 当微观小量 Δ ≪ E \Delta\ll E Δ ≪ E 但 Δ \Delta Δ 是宏观量不要把它和 N → ∞ N\to\infty N → ∞ 混为一谈 温标当成温度 P V = F ( T ) PV=F(T) P V = F ( T ) 对任何单值单调 F F F 都成立温度由第零定律定,温标还要定单位与原点
自测清单
能独立答出下面这些,这两章就算过了一遍。答不出的先记下来,回去翻对应的小节。
能说出第零定律的内容,并解释没有它温度为什么不是唯一确定的。
能写出状态方程的一般形式,并说清完备集里参数个数由什么决定。
能判断一个参数是广延还是强度,并说出这样分类在讨论子系统时的用处。
能说清水银温度计这类方案为什么不普适,理想气体温标为什么好一些。
能证明 T V = T P T_V=T_P T V = T P ,并指出哪一步用了 Boyle 定律。
能独立推出 V γ P = c o n s t V^\gamma P=\mathrm{const} V γ P = const 与 V γ − 1 T = c o n s t V^{\gamma-1}T=\mathrm{const} V γ − 1 T = const ,并说明哪一步用了理想气体。
能证卡诺效率是 1 − T 1 / T 2 1-T_1/T_2 1 − T 1 / T 2 ,并把它和绝对温标的定义对上。
能用反证法证卡诺定理,说清取整数 N , N ′ N,N' N , N ′ 的作用和”唯一效果”在哪里用到。
能默写 Kelvin 与 Clausius 两种表述,并能互相反推。
能用无穷多小卡诺循环证 ∮ r e v d Q / T = 0 \oint_{\mathrm{rev}}\dbar Q/T=0 ∮ rev d Q / T = 0 ,并由此说明 S S S 是态函数。
能说清可逆与不可逆过程算 Δ S \Delta S Δ S 的差别,并举出 d Q = 0 \dbar Q=0 d Q = 0 而 Δ S > 0 \Delta S>0 Δ S > 0 的例子。
能把 T d S = d U + P d V T\dd S=\dd U+P\dd V T d S = d U + P d V 改写成两条熵方程,并说明各自的独立变量。
能从四个势的全微分现场推出四条 Maxwell 关系,包括负号的位置。
能用 A A A 最小与 G G G 最小分别处理等温等容和等温等压的平衡。
能写出 Clapeyron 方程,并说明它为什么只对一级相变成立。
能写出 van der Waals 方程,并说清 a a a 与 b b b 各自的物理来源。
能由 ( V − V c ) 3 = 0 (V-V_c)^3=0 ( V − V c ) 3 = 0 三重根定出 T c T_c T c 、P c P_c P c 、V c V_c V c ,并写出对应态方程。
能用 Maxwell 等面积构造解释 vdW 等温线上的平台,分清亚稳段与不物理段。
能写出第三定律,并说明 T → 0 T\to0 T → 0 时 C → 0 C\to0 C → 0 是怎么从 S = ∫ 0 T C d T / T S=\int_0^TC\,\dd T/T S = ∫ 0 T C d T / T 得到的。
能写出等概率先验假设,并说清 Δ \Delta Δ 为什么必须引入。
能说清 Γ \Gamma Γ 、Σ \Sigma Σ 、ω \omega ω 的定义与差别,以及 S S S 用哪一个都行的理由。
能推出 S S S 是广延量,并指出”求和只有一项占主导”这一步为什么不可省。
能用 Γ ( E ) \Gamma(E) Γ ( E ) 定义温度与压强,并证明 P = − ( ∂ E / ∂ V ) S P=-(\partial E/\partial V)_S P = − ( ∂ E / ∂ V ) S 。
能独立走一遍能量均分的积分骨架,包括分部积分和边界项为零。
能说明 ⟨ H ⟩ = f k T / 2 \langle H\rangle=fkT/2 ⟨ H ⟩ = f k T /2 需要二次型,而 ⟨ x i ∂ H / ∂ x j ⟩ = δ i j k T \langle x_i\partial H/\partial x_j\rangle=\delta_{ij}kT ⟨ x i ∂ H / ∂ x j ⟩ = δ ij k T 不需要。
能算出理想气体的 Σ ( E ) \Sigma(E) Σ ( E ) ,并得到 U = 3 N k T / 2 U=3NkT/2 U = 3 N k T /2 与 P V = N k T PV=NkT P V = N k T 。
能说清 Gibbs 佯谬,说明 1 / N ! 1/N! 1/ N ! 修正的作用,以及 N ! N! N ! 与 N N N^N N N 两种取法的区别。
能用 Lagrange 乘子推出 n ˉ l = g l e − α − β ϵ l \bar n_l=g_le^{-\alpha-\beta\epsilon_l} n ˉ l = g l e − α − β ϵ l ,并用二阶变分说明它是极大。
能由 β = 1 / k T \beta=1/kT β = 1/ k T 把 Z Z Z 与 E E E 、P P P 、S S S 连起来,并写出自由粒子的 Z = V ( 2 π m k T ) 3 / 2 Z=V(2\pi mkT)^{3/2} Z = V ( 2 π mk T ) 3/2 。
能估算两个子系统能量分配的涨落,并说明它为什么正比于 k T 2 C V \sqrt{kT^2C_V} k T 2 C V 。
Comment评论