2025年大学热力学与统计物理试卷

2025年大学热力学与统计物理试卷考试时间：______分钟总分：______分姓名：______一、试述热力学第零定律，并说明由此如何定义温度。在一个由三个可自由移动的、内部能量为U（U>0）且化学成分均匀的物体组成的孤立系统内，发生热平衡过程。已知物体1和物体2的热容量分别为C₁和C₂，初始温度分别为T₁和T₂（T₁>T₂），试求热平衡时的共同温度T及过程的总熵变。二、一定量的理想气体经历一个由等温膨胀、绝热压缩和等温压缩三个过程组成的循环。已知气体在等温膨胀过程开始时的体积为V₁，压强为P₁，温度为T₀；等温膨胀到体积为V₂；绝热压缩回到体积V₁；再经等温压缩回到初始状态(P₁,V₁,T₀)。设气体为单原子分子，试求该循环的效率。三、根据玻尔兹曼分布，一个系统在温度T时，其粒子处在能量为εᵢ的状态的概率为ρᵢ*exp(-εᵢ/kT)，其中ρᵢ是与能量εᵢ无关的状态密度。试由此推导理想气体的内能U=3/2NkT（设气体为单原子分子，只考虑平动自由度）。并说明此结果与能量按自由度均分定理的关系。四、有一巨正则系综，其中包含N个可分辨粒子，每个粒子的单粒子能量可以是ε₁,ε₂,...,ε<0xE2><0x82><0x99>，相应的单粒子配分函数为λᵢ=exp(μ/kt)。试导出系统的内能U和熵S的表达式。其中μ为化学势，k为玻尔兹曼常量，T为温度，V为体积。五、范德瓦尔斯方程(P+a/V²)(V-b)=NkT修正了理想气体状态方程。其中a和b是常数，分别与气体分子间的吸引力和排斥力有关。试推导由范德瓦尔斯方程出发，理想气体的定容热容量Cv（设气体为单原子分子）与温度T的关系。说明与理想气体Cv=(3/2)Nk相比，有何不同。六、考虑一维无限深势阱中运动的粒子，其能量本征值为Eᵢ=(π²ħ²/2mL²)*ι²,ι=1,2,3,...，其中m为粒子质量，L为阱宽，ħ为约化普朗克常量。设系统处于基态，即波函数为ψ₁(x)=sqrt(2/L)*sin(πx/L)(0<x<L)。试计算在0<x<L区间内找到粒子的概率密度。现若系统处于第一激发态，即波函数为ψ₂(x)=sqrt(2/L)*sin(2πx/L)，再计算在x=L/4处找到粒子的概率密度。七、简述序参量的概念及其在相变研究中的作用。以液氦的λ相变为例，说明序参量是如何描述相变的。试卷答案一、解：由热力学第零定律，当三个物体达到热平衡时，它们具有相同的温度T。设物体1、2、3的初始温度分别为T₁、T₂、T₃（T₁>T₂>T₃，假设T₃不影响最终平衡，或可设T₃=T₁或T₂使问题简化），热容量分别为C₁、C₂、C₃。达到平衡时，物体1放热Q₁，物体2放热Q₂，物体3吸热Q₃。由于系统孤立，Q₁+Q₂=-Q₃。由Q=mcΔT，得C₁(T₁-T)+C₂(T₂-T)=-C₃(T-T₃)。若假设C₃=0或T₃即为所求平衡温度T，则T=(C₁T₁+C₂T₂)/(C₁+C₂)。总熵变ΔS=ΔS₁+ΔS₂+ΔS₃=C₁ln(T/T₁)+C₂ln(T/T₂)。（注：若T₃不为T，需将T₃代入上式联立求解，并计算所有物体熵变之和）。二、解：循环效率η=W/Q_H，其中W为循环净功，Q_H为循环中系统吸收的热量。等温膨胀（1→2）：Q₁=W₁=nRT₀*ln(V₂/V₁)，吸收热量。绝热压缩（2→3）：Q₂=0，外界对系统做功W₂=ΔU=nCv(T₃-T₀)=n(3/2)kT₀-n(3/2)kT₃（设单原子气体Cv=3/2k）。等温压缩（3→1）：Q₃=W₃=nRT₃*ln(V₁/V₂)，放出热量。净功W=W₁+W₂=nRT₀*ln(V₂/V₁)+n(3/2)k(T₀-T₃)。净吸热Q_H=Q₁=nRT₀*ln(V₂/V₁)。效率η=W/Q_H=[nRT₀*ln(V₂/V₁)+n(3/2)k(T₀-T₃)]/[nRT₀*ln(V₂/V₁)]=1+(3/2)(T₀-T₃)/T₀*[1/ln(V₂/V₁)]。（此处T₃由绝热方程P₂V₂ᵞ=P₃V₃ᵞ和P₃V₃=nRT₃求出，但为简化表达，未直接代入最终效率公式，实际计算需确定T₃与T₀、T₁的关系）。三、解：理想气体粒子总数为N。粒子处在能量为εᵢ的状态的概率为Pᵢ=ρᵢ*exp(-εᵢ/kT)/Σᵢρᵢ*exp(-εᵢ/kT)=λᵢ/Z，其中λᵢ=exp(μ/kt)=exp(εᵢ/kT)*exp(μ/kT)，μ为化学势，Z为配分函数。内能U=ΣᵢPᵢ*εᵢ=Σᵢ(λᵢ/kT)*εᵢ/Σᵢλᵢ/Z=(1/Z)*kT*Σᵢ(εᵢλᵢ)/Σᵢλᵢ。对理想气体，Z=Σᵢλᵢ=V^(N)*exp(Nμ/kT)=V^(N)*exp(Nεᵢ/kT)*exp(Nμ/kT)=V^(N)*exp(NU/kT)。（推导中用到N=Σᵢgᵢ*Pᵢ，gᵢ为状态i的简并度，此处未明确写出）。代入上式U=kT*(dlnZ/dT)=kT*(d/dT)[NlogV+NlogZ/kT]=kT*[N/V*(dV/dT)-N/kT*(dZ/dT)+N/kT²*Z]。由热力学dU=TdS-PdV，得(∂U/∂T)_V=T(∂S/∂T)_V-P=Cv。又(∂U/∂T)_V=kT*(d/dT)[NlogZ/kT]=k*(N/kT²*Z-N/kT*(dZ/dT))=N/kT*(Z-T(dZ/dT))。令β=1/kT，则U=kT*[N/V*(dV/dT)-N/kT*(dZ/dT)+N/kT²*Z]=N/kT*[Vβ(dV/dT)-Zβ(dZ/dT)+Z]=N/kT*[Vβ(dV/dT)-Zβ(dZ/dT)+Z]=N/kT*[Vβ(dV/dT)-Zβ(dZ/dT)+Z]=N/kT*[Vβ(dV/dT)-Zβ(dZ/dT)+Z]。对单原子理想气体，N/V=n，β(dZ/dT)=-N/kT*(dZ/dT)=-P，代入得U=N*[β(dV/dT)-P+βZ]=N*[β(dV/dT)-P+βZ]=N*[β(dV/dT)-P+βZ]=N*[β(dV/dT)-P+βZ]。利用理想气体状态方程PV=NkT，得U=(3/2)NkT。四、解：巨正则系综包含所有可能粒子数N、能量E和粒子化学势μ的微观状态数。巨配分函数Ω=Σᵢgᵢ*exp[(μN-Eᵢ)/kT]，其中gᵢ为能量为εᵢ的单粒子能级简并度。令x=exp(μ/kT)，则Ω=Σᵢgᵢ*x^N*exp(-εᵢ/kT)=x^N*Σᵢgᵢ*exp(-εᵢ/kT)=x^N*Z，其中Z为正则配分函数。内能U=-(∂lnΩ/∂β)=-(∂ln(x^N*Z)/∂β)，其中β=1/kT。由于β=1/kT，对x求导等于乘以β，对Z求导不变，得U=-N*x*(∂lnZ/∂β)=-N*x*(-kT^2*P)，其中P=-(∂lnZ/∂β)为正则分布下的压强。所以U=NkT^2*(∂lnZ/∂β)=Nμ。（推导中使用了U=-(∂lnZ/∂β)）。熵S=k*(lnΩ-β*(∂lnΩ/∂β))=k*[ln(x^N*Z)-β*d(x^N*Z)/dx]=k*[Nlnx+lnz-β*(N*x*lnx+x^N*dlnz/dx)]。将β=1/kT，dlnz/dx=(1/Z)*(∂Z/∂x)代入，得S=k*[Nlnx+lnz-(1/kT)*(N*x*lnx+x^N*(1/Z)*(∂Z/∂x))]=k*[Nlnx+lnz-Nlnx-x^N*P/(x*Z)]=k*[lnz-x^N*P/(x*Z)]=k*[lnz-P]。（推导中使用了S=k*(lnΩ-β*(∂lnΩ/∂β)），以及P=-(∂lnZ/∂β)）。五、解：由范德瓦尔斯方程P+a/V²=nRT/V-nb，得P=nRT/V-nb-a/V²。定容热容量Cv=(∂U/∂T)_V。内能U=3/2NkT（理想气体部分）+其他项（与a,b有关）。由热力学dU=TdS-PdV，在定容过程(dV=0)下，dU=TdS。因此Cv=T(∂S/∂T)_V。利用熵的微分式dS=nCv(dT)/T+n(∂P/∂T)_V(dV)，在定容下，dS=nCv(dT)/T。对于范德瓦尔斯气体，P=nRT/V-nb-a/V²，故(∂P/∂T)_V=nR/V。代入得dS=nCv(dT)/T+n(nR/V)dV=nCv(dT)/T。定容过程dV=0，故dS=nCv(dT)/T。所以Cv=T(dS/dT)=nCv。（这里似乎直接回到了定容热容量的定义，需要重新思考如何体现a和b的影响）。更准确的方法是，范氏气体内能U=3/2NkT+U₂(V,T)，其中U₂包含了a/V²和-Pb的影响。由热力学dU=TdS-PdV，得dS=(dU/dT)_V*(dT/T)+(dU/dV)_T*(dV/V)。在定容下，dS=(dU/dT)_V*(dT/T)。定容热容量Cv=T(dS/dT)_V=T*[(dU/dT)_V+(dU/dV)_T*(dV/dT)_V/T]=T*[(∂U/∂T)_V+(∂U₂/∂V)_T*(1/T)]。U₂依赖于V和T。如果假设U₂=U₂'(T)+U₂''(V)，则(dU₂/dT)_V=(dU₂'/dT)_V，(dU₂/dV)_T=(dU₂''/dV)_T。由P=nRT/V-nb-a/V²，得(∂U₂/∂V)_T=-Pb-a/V²*(∂V/∂V)_T=-Pb-a/V³。代入得Cv=T*[(∂U/∂T)_V+(-Pb-a/V³)*(1/T)]=Cv_ideal+T*(-Pb/V³-a/V³)。理想气体Cv=3/2Nk。范氏气体Cv=3/2Nk+T*(-Pb/V³-a/V³)。将P=nRT/V-nb-a/V²代入，得Cv=3/2Nk+T*[(-nRT/V+nb+a/V²)b/V³-a/V³]=3/2Nk+T*[-nRTb/V⁴+nb²/V³+ab/V⁵-a/V³]=3/2Nk-nRbT/V+nb²/V³+abT/V⁵-aT/V³。对于高温低密度极限，V很大，a/V⁵和b/V³可忽略，得Cv≈3/2Nk-nRbT/V。对于低温或常温，此项可能与3/2Nk相抵消甚至改变符号。这表明范氏气体的定容热容量随温度和密度变化，与理想气体不同。六、解：概率密度即波函数的模平方。(1)对于基态ψ₁(x)=sqrt(2/L)*sin(πx/L)(0<x<L)，概率密度为|ψ₁(x)|²=(2/L)*sin²(πx/L)。在0<x<L区间内找到粒子的概率为∫[₀ᵐ]|ψ₁(x)|²dx=∫[₀ᵐ](2/L)*sin²(πx/L)dx。利用sin²θ=(1-cos(2θ))/2，得∫[₀ᵐ]sin²(πx/L)dx=∫[₀ᵐ](1-cos(2πx/L))/2dx=(1/L)*[x-(L/2π)*sin(2πx/L)]|₀ᵐ=m/L-(L/2π)*sin(2πm/L)+0-(L/2π)*sin(0)=m/L-(L/2π)*0=m/L。所以概率为m/L=1。（因为积分范围是完整周期）。更标准的计算为∫[₀ᵐ](2/L)*sin²(πx/L)dx=(2/L)*[x/2-(L/4π)*sin(2πx/L)]|₀ᵐ=(2/L)*[m/2-(L/4π)*sin(2πm/L)]=m/L-(1/2π)*sin(2πm/L)。对于整数m，sin(2πm/L)=0，所以概率为m/L。(2)对于第一激发态ψ₂(x)=sqrt(2/L)*sin(2πx/L)(0<x<L)，概率密度为|ψ₂(x)|²=(2/L)*sin²(