24686 words
123 minutes
Loading... readings
热力学与统计物理

该笔记主要参考热力学部分——《热学 热力学与统计物理(上册)》(曹烈兆 周子舫),统计物理部分——《热学 热力学与统计物理(下册)》(周子舫 曹烈兆),以及《热力学与统计物理》(汪志诚)。

  

  

  


第1章 热力学的基本规律#

1.1 热力学状态与第零定律#

1.1.1 热力学平衡状态#

  • 要建立热力学这座大厦,第一步是确定我们要研究的对象——平衡态。宏观体系包含极其巨大的粒子数,试图追踪每一个粒子的力学运动是绝望的。热力学的聪明之处在于,它只关注系统在不受外界干扰、经过足够长时间演化后,最终“稳定”下来的状态。这里的“稳定”绝不是微观上的死寂,而是宏观量不再随时间变化,而微观粒子依然在剧烈且无序地运动,即所谓 “宏观静止,微观涨落”

    从非平衡走向平衡所需的时间尺度,我们称为弛豫时间 (τ\tau)。这个概念在物理图像中极其关键。在实际的物理过程里,所有的测量和操作所经历的时间,都必须与弛豫时间进行对比。只有当外界条件变化的时间尺度(tt)远大于 τ\tau 时(tτt \gg \tau),系统才有时间对外界变化做出响应并随时重新达到平衡,我们才能把整个过程视为一连串平衡态的连缀(即准静态过程),此时定义出来的热力学量才是确切且有意义的。

1.1.2 热力学第零定律与温度#

  • 既然定义了单个系统的平衡态,我们自然会进一步考察多个系统之间的相互作用。如果把两个独立的系统进行热接触,它们最终会达到一个共同的新平衡态。那么,如何在接触前,就预判它们之间会不会发生热量交换?这就需要为“热平衡”寻找一个普适的宏观判据。

    为了从逻辑上严密地建立这个判据,物理学引入了热力学第零定律

    “如果两个热力学系统中的每一个都与第三个热力学系统处于热平衡,那么它们也必定互相处于热平衡。”

    从日常经验上看,这似乎是一句理所当然的废话,但在物理逻辑上,它建立了一个严格的数学“等价关系”。这种等价关系的物理本质在于,它强行要求必须存在一个共有的状态函数——温度 (TT)

    第零定律告诉我们,温度不仅仅是感知上的“冷热程度”,更是决定两个系统接触时是否会发生热交换的唯一宏观判据。有了这个定律兜底,我们就可以理直气壮地选取某个特定系统(温度计)的特定物理属性(如水银的体积)作为标准,建立起经验温标。

1.1.3 物态方程与热力学系数#

  • 通过第零定律确立了温度 TT 后,加上我们在力学和几何上熟知的压强 pp 和体积 VV 等宏观参量,我们就有了一套描述系统的语言。对于一个确定的均匀物质系统,其实际的经验告诉我们,只要给定其中几个独立的参量,它的平衡态也就随之唯一确定了。由此引入物态方程(通常写为 f(p,V,T)=0f(p, V, T) = 0)对物质特定平衡状态进行宏观描述。它告诉我们,在这个特定的平衡态下,这几个宏观量是如何具体联系在一起的。

    然而在真实的物理世界里,除了结构简单的气体,对于大多数固体或液体,我们几乎不可能从理论上写出一个完美的、全局适用的物态方程解析式。那在实验室里,引入三个极其重要的宏观实验测量量(热力学系数)

    1. 体胀系数 (Volume Expansivity) α=1V(VT)p\alpha = \frac{1}{V} \left( \frac{\partial V}{\partial T} \right)_p

      压强保持不变(等压过程)的条件下,温度每升高 1K 时,系统体积的相对增加量。它直接反映了物质“热胀冷缩”的剧烈程度。

    2. 等温压缩系数 (Isothermal Compressibility) κT=1V(Vp)T\kappa_T = -\frac{1}{V} \left( \frac{\partial V}{\partial p} \right)_T

      温度保持不变(等温过程)的条件下,压强每增加 1Pa 时,系统体积的相对缩小量。

    3. 压强系数 (Pressure Coefficient) β=1p(pT)V\beta = \frac{1}{p} \left( \frac{\partial p}{\partial T} \right)_V

      体积保持不变(等容过程)的条件下,温度每升高 1K 时,系统内部压强的相对增加量。它反映了如果把物质锁死在一个刚性容器里加热,其内部压强随温度攀升的敏锐程度。

    这三个系数构成了联系物质宏观状态变化的“铁三角”,通过数学上的偏导数循环关系,它们满足严格的制约方程:α=pβκT\alpha = p \beta \kappa_T(推到利用了(Vp)T(pT)V(TV)p=1\left( \frac{\partial V}{\partial p} \right)_{T} \left( \frac{\partial p}{\partial T} \right)_{V} \left( \frac{\partial T}{\partial V} \right)_{p} = -1)。在实际研究中,只要测出其中两个,就能直接算出第三个。在不知道复杂物质确切物态方程的情况下,只要我们在实验室里测出了该物质的这三个表观系数,就能在局部范围内完全掌握它的状态演化规律。

  • 虽然复杂物质极度依赖这些实验系数,但对于最简单的一类物质——气体,我们确实可以建立直观的宏观物态方程模型。最简单的玩具模型是理想气体物态方程pV=nRTpV = nRT),它在极高温度和极低压强下成立,代价是直接抹除了微观细节:把粒子视为没有体积的几何点,并切断了所有相互作用。

    理想气体物态方程的实验推导

    理想气体的宏观物态方程是直接建立在三个控制变量的经典物理实验之上:

    玻意耳定律(等温实验):保持气体温度 TT 和物质的量 nn 不变,压缩或膨胀气体。实验观测发现,气体的体积与压强严格成反比:V1pV \propto \frac{1}{p}

    查理-盖·吕萨克定律(等压实验):保持气体压强 pp 和物质的量 nn 不变,给气体加热。实验观测发现,气体的体积与绝对温度成正比:VTV \propto T

    阿伏伽德罗定律(粒子数实验):在相同的温度 TT 和压强 pp 下,任何气体的体积都直接正比于其包含的分子总数(即物质的量 nn):VnV \propto n

    既然气体的体积 VV 既正比于温度 TT 和物质的量 nn,又反比于压强 pp,我们在数学上可以引入一个普适的比例常数——理想气体常数 RR(由实验精确标定得出),直接将这三个独立的比例关系合并为一个统一的表达式:

    pV=nRTpV = nRT

    但只要稍微靠近真实的物理世界,比如引入范德瓦尔斯方程(p+an2V2)(Vnb)=nRT(p + \frac{an^2}{V^2})(V - nb) = nRT,物理图像立马就丰满了。等式左边多出来的这两个修正项,本质上是对微观粒子间相互作用的“粗视化”补偿:

    内压项修正 (+an2V2+\frac{an^2}{V^2}):源于分子间的长程引力。靠近容器壁的分子被内部的分子往回拉,导致它们撞击器壁的实际压强变小了。因此,我们在宏观测量的压强 pp 上,必须把它损失的这部分给“加回来”,才能描述系统内部真实的受力状态。

    共体积修正 (nb-nb):源于分子间的短程斥力。粒子本身占据了一定的体积,存在不可穿透的硬核,导致气体分子真正能自由活动的空间其实比容器的宏观总体积 VV 要小,所以要减去一个“禁带”体积。

    范德瓦尔斯方程修正项推导

    设理想气体状态方程为 pidVid=nRTp_{id} V_{id} = nRT。实际气体在微观上具有分子本身体积(表现为短程斥力)和分子间相互作用力(表现为长程引力)。

    1. 体积修正(共体积项)

    假设气体分子为直径为 dd 的刚性球,单个分子的真实体积为:

    v0=43π(d2)3=16πd3v_0 = \frac{4}{3}\pi \left(\frac{d}{2}\right)^3 = \frac{1}{6}\pi d^3

    当两个分子发生碰撞时,它们的中心距离最小为 dd。即以任意一个分子中心为球心,半径为 dd 的球状空间内,不能有其他分子中心进入。两个分子共同的“排除体积”(Excluded Volume)为:

    Vexcl=43πd3=8(16πd3)=8v0V_{excl} = \frac{4}{3}\pi d^3 = 8 \left(\frac{1}{6}\pi d^3\right) = 8v_0

    分摊到平均每个分子上的排除体积为:

    vb=12Vexcl=4v0v_b = \frac{1}{2} V_{excl} = 4v_0

    对于 1 mol1 \text{ mol} 气体(NAN_A 个分子),其共体积常量 bb 为:

    b=NAvb=4NAv0b = N_A v_b = 4 N_A v_0

    对于 n moln \text{ mol} 气体,分子实际可自由运动的有效体积 VidV_{id} 为系统总体积 VV 减去所有分子的排除体积:

    Vid=VnbV_{id} = V - nb

    2. 压强修正(内压项)

    处于容器内部的分子受周围分子引力对称,合力为零。但处于器壁表面厚度为分子引力有效作用半径的表面层内的分子,受到内部气体的净引力(指向容器内部)。导致实际气体对器壁的碰撞压强 pp 相比理想压强 pidp_{id} 减小了一个“内压” Δp\Delta p

    p=pidΔp    pid=p+Δpp = p_{id} - \Delta p \implies p_{id} = p + \Delta p

    内压 Δp\Delta p 的大小取决于两个因素:

    1. 表面层内受到内部引力拉扯的分子数密度 ρsurf\rho_{surf}

    2. 内部对表面层分子施加引力的分子数密度 ρint\rho_{int}

    由于在平衡态下,各处分子数密度是均匀的,即 ρsurfnV\rho_{surf} \propto \frac{n}{V}ρintnV\rho_{int} \propto \frac{n}{V}。因此,内压 Δp\Delta p 正比于密度的平方:

    Δpρsurfρint(nV)(nV)=(nV)2\Delta p \propto \rho_{surf} \cdot \rho_{int} \propto \left(\frac{n}{V}\right) \cdot \left(\frac{n}{V}\right) = \left(\frac{n}{V}\right)^2

    引入引力比例常数 aa,得到内压项:

    Δp=an2V2\Delta p = a \frac{n^2}{V^2} pid=p+an2V2p_{id} = p + a \frac{n^2}{V^2}

    pidp_{id}VidV_{id} 的修正表达式代入理想气体状态方程 pidVid=nRTp_{id} V_{id} = nRT,得到范德瓦尔斯方程:

    (p+an2V2)(Vnb)=nRT\left( p + a \frac{n^2}{V^2} \right) (V - nb) = nRT

1.2 热力学第一定律与能量转化#

1.2.1 功与内能#

  • 要改变一个热力学系统的状态,宏观上只有两种手段:做功(如压缩气体、搅拌液体)和传热(如放在火上加热)。功 (WW) 和热量 (QQ) 都是伴随“过程”而生的,它们不仅取决于系统的初末状态,还严重依赖于系统经历的具体路径。在数学上,它们不是全微分(通常用 δW\delta WδQ\delta Q 表示)。

    在热力学中,功不仅局限于机械体积功(pdVp dV)。任何宏观做功过程都可以写成广义力 (YY)广义位移 (dxdx) 的乘积:δW=Ydx\delta W = Y dx。我们需要先计算外界对系统做的功 đWexđW_{ex},则系统对外的元功 δW=đWex\delta W = - đW_{ex}。接下来以电磁介质为例,展示元功的推导过程。在处理电磁介质时,最核心的物理直觉是:必须把“建立真空电磁场消耗的能量”和“真正改变介质内部状态(极化/磁化)所消耗的能量”剥离开来。 只有后者才属于电/磁介质这个热力学系统的“内能”变化。

    1. 电介质的极化功

    物理模型:考虑一个充满各向同性电介质的平行板电容器。极板面积为 AA,间距为 dd,介质体积 V0=AdV_0 = Ad

    外界总电功:当外界电源给极板微小充电 dqdq 时,极板间电势差为 UU。外界做功为:

    đWex=UdqđW_{ex} = U dq

    引入场量:根据电磁学,电势差 U=EdU = Ed,极板自由电荷 q=DAq = DADD 为电位移)。对其微分有 dq=AdDdq = A dD。代入上式:

    đWex=(Ed)(AdD)=V0EdDđW_{ex} = (Ed)(A dD) = V_0 E dD

    剥离真空场能量:介质中的电位移本构关系为 D=ε0E+PD = \varepsilon_0 E + PPP 为极化强度)。代入微分式:

    đWex=V0Ed(ε0E)+V0EdP=d(12ε0E2V0)+V0EdPđW_{ex} = V_0 E d(\varepsilon_0 E) + V_0 E dP = d\left(\frac{1}{2}\varepsilon_0 E^2 V_0\right) + V_0 E dP

    等式右边第一项 d(12ε0E2V0)d(\frac{1}{2}\varepsilon_0 E^2 V_0) 是即使电容器里是真空也必须提供的能量,它储存在真空电场中,与电介质材料的自身热力学状态无关。等式右边第二项 V0EdPV_0 E dP 才是外界电场真正为了抵抗分子内部束缚力、拉扯电荷中心产生极化,从而对电介质材料本身做的功。令系统总电偶极矩 p=V0Pp = V_0 P,则外界对介质做的极化功为 EdpE dp。因此,电介质系统对外所做的元功表达式为:

    δW=Edp\delta W = - E dp
    1. 磁介质的磁化功

    物理模型:考虑一个内部充满磁介质的细长螺线管。截面积为 AA,长度为 ll,体积 V0=AlV_0 = Al,总匝数为 NN

    外界总电磁功:当线圈电流 II 发生变化时,会引起磁通量 Φ=NBA\Phi = NBA 的变化,从而产生感应电动势 E=dΦ/dt\mathcal{E} = -d\Phi/dt。外界电源为了维持电流,必须克服感应电动势做功:

    đWex=EIdt=IdΦ=INAdBđW_{ex} = -\mathcal{E} I dt = I d\Phi = I N A dB

    引入场量:根据安培环路定理,磁场强度 H=NI/lH = NI / l,即 IN=HlIN = Hl。代入上式:

    đWex=(Hl)AdB=V0HdBđW_{ex} = (Hl) A dB = V_0 H dB

    剥离真空场能量:介质中的磁感应强度本构关系为 B=μ0(H+M)B = \mu_0 (H + M)MM 为磁化强度)。代入微分式:

    đWex=V0Hd(μ0H)+V0Hd(μ0M)=d(12μ0H2V0)+μ0V0HdMđW_{ex} = V_0 H d(\mu_0 H) + V_0 H d(\mu_0 M) = d\left(\frac{1}{2}\mu_0 H^2 V_0\right) + \mu_0 V_0 H dM

    同样地,第一项 d(12μ0H2V0)d(\frac{1}{2}\mu_0 H^2 V_0) 是建立真空磁场所需的能量。第二项 μ0V0HdM\mu_0 V_0 H dM 才是外界磁场为了使介质内部微观磁矩发生偏转定向,从而对磁介质材料本身做的功。令系统总磁矩 m=V0Mm = V_0 M,则外界对磁介质做的磁化功为 μ0Hdm\mu_0 H dm。因此,磁介质系统对外所做的元功表达式为:

    δW=μ0Hdm\delta W = - \mu_0 H dm

    但物理学家在实验(焦耳的绝热做功实验)中发现:在一个绝热系统中(没有热量交换,δQ=0\delta Q=0),无论你用什么方式对系统做机械功,只要系统的初末状态确定了,外界对系统所做的总功总是完全相等的。

    这个实验事实意味着,虽然功 WW 本身是依赖路径的,但在“绝热”这个特定限制下,它变成了一个只由初末状态决定的量。在物理逻辑上,这就强制要求系统内部必须存在一个状态函数,它的变化量刚好等于绝热过程中的功。将这个状态函数定义为内能 (UU),其是系统内在的、由微观粒子无规则运动动能和分子间相互作用势能总和所构成的宏观状态函数。

1.2.2 热力学第一定律#

  • 在真实的物理世界中,系统不仅可以与外界交换机械功,还可以通过温差直接交换热量,热量的进入显然也会引起系统状态的变化。宏观物理学给出了一个极其坚定的论断——无论是通过做功传递的宏观机械能,还是通过传热传递的微观分子热运动能量,它们在改变系统内能上是完全等效的,且总能量必须守恒。这就是热力学第一定律的物理内核,是能量守恒定律在涉及热现象的宏观过程中的直接表述

    用数学语言来刻画这个物理图景,就是:系统吸收的热量 (δQ\delta Q),一部分用于增加自身的内能 (dUdU),另一部分用于对外做膨胀功 (δW\delta W),即:

    δQ=dU+δW\delta Q = dU + \delta W

    在这个微分方程中,左边的传热和右边的做功都是依赖路径的“过程量”,但神奇的是,当它们按照这个守恒法则组合在一起时,其差值 (dUdU) 却是一个与路径无关的态函数。

1.2.3 焓 (HH) 与热容#

  • 如果把系统锁死在一个刚性容器里(等容过程,dV=0dV = 0),系统对外不做体积功,吸收的热量就完完全全变成了内能的增加:δQV=dU\delta Q_V = dU。但在实验室里,甚至在大自然中,绝大多数的物理和化学过程并不是等容的,而是等压的(比如在敞口烧杯中发生的反应,系统始终承受着恒定的大气压)。系统在吸热的同时往往伴随着体积的膨胀,这意味着它必须推开周围的空气做功(δW=pdV\delta W = p dV)。

    把这个恒定压强下的做功代入第一定律,我们会得到 δQp=dU+pdV\delta Q_p = dU + p dV。因为压强 pp 是常数,它可以被毫无阻碍地塞进微分号里,变成:

    δQp=d(U+pV)\delta Q_p = d(U + pV)

    物理学家为了处理这类极其常见的等压过程,直接把括号里这个由态函数组合而成的新态函数 (U+pV)(U + pV) 定义为焓 (HH)。它的出现完全是为了物理分析上的便利:在等压过程中,系统吸收的热量直接等于其焓的增量 (δQp=dH\delta Q_p = dH)。

    为了描述系统吸热量与温度变化的关系,我们定义了两个宏观可测的热容

    定容热容 (CVC_V)CV=(UT)VC_V = \left(\frac{\partial U}{\partial T}\right)_V

    定压热容 (CpC_p)Cp=(HT)pC_p = \left(\frac{\partial H}{\partial T}\right)_p

    等压升温时,系统吸收的热量不仅要用于增加内能(升温),还要分出一部分能量去抵抗外压做膨胀功,这导致宏观上 CpC_p 必然大于 CVC_V

    理想气体内能表达式的推导

    在焦耳自由膨胀实验中,理想气体向绝热的真空容器中膨胀。因为绝热(δQ=0\delta Q = 0)且向真空膨胀不克服任何阻力(不对外做功,δW=0\delta W = 0),根据热力学第一定律,系统的内能在膨胀前后保持恒定(dU=0dU = 0)。同时,实验精确观测到气体的温度也未发生改变(dT=0dT = 0)。既然体积 VV 经历了剧烈增加,而内能 UU 和温度 TT 却都不变,这在物理上直接导出了一个坚实的结论(焦耳定律):理想气体的内能与体积完全无关,它仅仅是温度的单变量函数,即 U=U(T)U = U(T)

    既然明确了 UU 只是 TT 的函数,我们在定义定容热容 CV=(UT)VC_V = \left(\frac{\partial U}{\partial T}\right)_V 时所用的偏导数,就可以理直气壮地退化为全导数形式:

    CV=dUdT    dU=CVdTC_V = \frac{dU}{dT} \implies dU = C_V dT

    对该微分方程进行积分,即可直接得到理想气体的宏观内能表达式:

    U(T)=CVdT+U0U(T) = \int C_V dT + U_0

    (其中 U0U_0 为积分常数,代表参考态下的内能。由于理想气体的定容热容 CVC_V 在相当宽的温度范围内可视为常数,该表达式在实际计算中通常直接简化为 U=CVT+U0U = C_V T + U_0)。

1.2.4 理想气体的绝热过程#

  • 绝热过程的定义是系统与外界无热交换,即 δQ=0\delta Q = 0。根据热力学第一定律 δQ=dU+δW\delta Q = dU + \delta W,可直接得到微积分方程:

    dU=δW    CVdT=pdVdU = -\delta W \implies C_V dT = -p dV

    气体对外膨胀做功(pdV>0p dV > 0)的能量,完全由消耗系统自身的内能(CVdT<0C_V dT < 0,温度降低)来支付。将理想气体状态方程 pV=nRTpV = nRT 两边取全微分:

    pdV+Vdp=nRdTp dV + V dp = nR dT

    考虑等压升温比等容升温多吸收的热量,刚好等于气体对外做的膨胀功:CpdT=CVdT+pdVC_p dT = C_V dT + p dV,由理想气体状态方程推得等压膨胀功为 pdV=nRdTp dV = nR dT。代入上式并约掉 dTdT,得 CpCV=nRC_p - C_V = nR,利用此关系代入上式消去 nRnR

    pdV+Vdp=(CpCV)dTp dV + V dp = (C_p - C_V) dT

    将上述 dTdT 的表达式代回第一定律的方程 CVdT=pdVC_V dT = -p dV 中并整理等式两边:

    CVVdp+CppdV=0C_V V dp + C_p p dV = 0

    为了求解这个微分方程,将两边同除以 CVpVC_V p V,并引入绝热指数 γ=CpCV\gamma = \frac{C_p}{C_V}

    dpp+γdVV=0\frac{dp}{p} + \gamma \frac{dV}{V} = 0

    对等式两边直接积分:

    lnp+γlnV=常数\ln p + \gamma \ln V = \text{常数}

    最终得到绝热过程的核心特征方程(泊松方程):

    pVγ=常数p V^\gamma = \text{常数}

    对特征方程求导可得绝热线的斜率为 dpdV=γpV\frac{dp}{dV} = -\gamma \frac{p}{V},而等温线的斜率为 pV-\frac{p}{V}。因为 γ>1\gamma > 1,绝热线比等温线更陡峭。在绝热膨胀中,压强的下降由双重因素叠加导致——不仅是因为体积变大(空间变得稀疏),更是因为内能转化为功导致了温度的剧烈下降(分子热运动动能锐减)。

1.3 第二定律与不可逆性#

  • 热力学第一定律确立了能量在转化过程中的守恒关系,但并未对过程的演化方向做出任何限制。在实际的宏观物理世界中,能量的转化和传递具有明确的方向性,第二定律正是对这种宏观不可逆性的严格物理界定。

1.3.1 过程的方向性与第二定律的表述#

  • 自然界中的自发宏观过程均具有单向性。例如,机械功可以自发且无代价地完全转化为热(如摩擦生热),而热量却不能在不产生其他宏观影响的条件下完全转化为功;同样,热量必然自发地从高温物体传向低温物体,而无法自动逆向传递。为了对这种方向性进行确切的物理限定,热力学第二定律给出了两种经典的表述:

    1. 开尔文表述(针对热功转换):不可能从单一热源吸取热量,使之完全变为有用功而不产生其他影响。

      实质:确立了热机效率的绝对上限必然小于 100%。热机在持续输出机械功的过程中,必须向低温热源排放废热。

    2. 克劳修斯表述(针对热量传递):不可能把热量从低温物体传到高温物体而不产生其他影响。

      实质:指明了热量传递的自发方向。制冷机若要实现热量的逆向传递,必须依赖外界对其做功。

    这两种表述在物理上是完全等价的。若其中之一不成立,便可以推出另一个也不成立。它们从不同视角共同确立了宏观热力学过程的不可逆性。

    1. 若开尔文表述不成立,则必然导致克劳修斯表述不成立

    假设存在一台违反开尔文表述的“完美热机”(记为机 A)。该热机能够从单一的高温热源(温度为 T1T_1)吸收热量 QQ,并将其完全转化为宏观机械功 WW(即 W=QW=Q),且不对外界产生任何其他影响。

    现引入一台遵循物理规律的普通制冷机(记为机 B),使其在相同的高温热源 T1T_1 和低温热源 T2T_2 之间工作。利用机 A 输出的功 WW 来直接驱动机 B。机 B 消耗功 WW,从低温热源吸取热量 Q2Q_2,并向高温热源排放热量 Q1Q_1'。根据第一定律的能量守恒,Q1=W+Q2=Q+Q2Q_1' = W + Q_2 = Q + Q_2

    将机 A 与机 B 视为一个闭合的联合系统。分析该系统完成一个循环后的宏观净效应:

    外界净做功:联合系统对外输出功 WW(机 A),同时内部消耗功 WW(机 B),与外界的净功交换为零。 低温热源 T2T_2:被机 B 吸取了热量 Q2Q_2高温热源 T1T_1:向机 A 输出了热量 QQ,同时接收了机 B 排放的热量 Q+Q2Q + Q_2。其净效应为接收了热量 Q2Q_2

    综合来看,该联合系统的唯一最终效果是:自发地将热量 Q2Q_2 从低温热源传到了高温热源,而没有消耗外界的功,也没有引起其他任何状态变化。这直接违反了克劳修斯表述。

    2. 若克劳修斯表述不成立,则必然导致开尔文表述不成立

    假设存在一台违反克劳修斯表述的“完美制冷机”(记为机 C)。该制冷机能够在不消耗任何外界功的条件下,自动从低温热源 T2T_2 吸取热量 Q2Q_2,并将其全部排入高温热源 T1T_1

    现引入一台遵循物理规律的普通热机(记为机 D),使其在相同的两热源间工作。设定机 D 的规模,使其在每个循环中恰好向低温热源 T2T_2 排放热量 Q2Q_2。同时,机 D 从高温热源吸取热量 Q1Q_1,并对外做功 WW。根据第一定律,W=Q1Q2W = Q_1 - Q_2

    将机 C 与机 D 构成一个联合系统。分析其宏观净效应:

    低温热源 T2T_2:机 D 向其排放热量 Q2Q_2,机 C 恰好从其吸取热量 Q2Q_2。该热源的净热量交换为零,状态完全未变。 高温热源 T1T_1:向机 D 输出了热量 Q1Q_1,同时接收了机 C 排放的热量 Q2Q_2。其净效应为向系统输出了热量 Q1Q2Q_1 - Q_2外界净做功:系统对外输出净功 W=Q1Q2W = Q_1 - Q_2

    综合来看,该联合系统的唯一最终效果是:从单一热源 T1T_1 吸收了热量 (Q1Q2)(Q_1 - Q_2),并将其完全转化为有用的机械功 WW,且未对低温热源或外界环境产生任何其他影响。这直接违反了开尔文表述。

1.3.2 卡诺循环与卡诺定理#

  • 在第二定律确立了热功转化的方向性约束后,确定热机在给定高低温热源之间工作的理论效率上限便成为核心问题。为此,卡诺引入了由两个准静态等温过程和两个准静态绝热过程构成的理想可逆循环——卡诺循环

    基于卡诺循环模型及第二定律的开尔文表述,通过逻辑反证法可以严格推导出卡诺定理

    1. 在相同的高温热源和相同的低温热源之间工作的一切实际(不可逆)热机,其效率都不可能大于可逆热机。

    2. 在相同的高温热源和低温热源之间工作的一切可逆热机,其效率都相等,且与工作物质的具体性质无关。

    卡诺定理为热功转换设定了不可逾越的理论极限。其第二条推论表明,可逆热机的效率 η\eta 仅由高、低温热源的状态决定,彻底摆脱了具体热机结构与工作介质属性的制约。

1.3.3 热力学绝对温标的建立#

  • 卡诺定理中“可逆热机效率与工作物质无关”的结论,为建立完全独立于测温物质的温度标尺提供了严格的理论基础。

    根据卡诺定理,在温度为 t1t_1t2t_2 的两个热源间工作的可逆热机,其吸热量 Q1Q_1 与放热量 Q2Q_2 之比 Q1Q2\frac{Q_1}{Q_2} 必然是一个不依赖于具体物质的普适函数:Q1Q2=f(t1,t2)\frac{Q_1}{Q_2} = f(t_1, t_2)

    开尔文利用这一普适关系,直接规定两个热源的绝对温度(记为 T1T_1T2T_2)之比,等于在此两热源间工作的可逆卡诺热机吸收与放出的热量之比,即:

    T1T2=Q1Q2\frac{T_1}{T_2} = \frac{Q_1}{Q_2}

    由此定义的温标称为热力学绝对温标(或开氏温标)。它建立在热力学第二定律的坚实基础之上,彻底排除了经验温标对特定测温物质(如水银、酒精等)物理属性的依赖。进一步地,若将理想气体作为工作物质代入卡诺循环进行推导,可证明热力学绝对温标与理想气体温标在数值上是严格等价的。

1.4 熵与热力学基本方程#

  • 卡诺定理不仅确立了不依赖具体物质的绝对温标,其更深刻的数学推论在于揭示了一个全新的状态函数。本节将从热量传递的普遍规律中严格导出“熵”,并建立宏观热力学的基本微分方程。

1.4.1 克劳修斯等式与熵函数的导出#

  • 在可逆卡诺循环中,工作物质从高温热源 T1T_1 吸热 Q1Q_1,向低温热源 T2T_2 放热 Q2Q_2。根据绝对温标的定义,有 Q1Q2=T1T2\frac{Q_1}{Q_2} = \frac{T_1}{T_2}。若引入热力学通用的符号法则(系统吸热为正,放热为负),将放热量记为 Q2-Q_2,则上式可改写为:

    Q1T1+Q2T2=0\frac{Q_1}{T_1} + \frac{Q_2}{T_2} = 0

    对于任意复杂的闭合可逆循环,可通过引入无数条绝热线,将其分割为无数个微小的可逆卡诺循环。 对所有微小循环求和取极限,可得到克劳修斯等式

    δQrevT=0\oint \frac{\delta Q_{\text{rev}}}{T} = 0

    这表明沿任意可逆闭合路径的积分为零。在数学上,这严格证明了从平衡态 A 到平衡态 B 的积分 ABδQrevT\int_A^B \frac{\delta Q_{\text{rev}}}{T} 其数值完全独立于具体的积分路径,而仅仅取决于系统的初末状态。既然该积分值只与状态有关,物理学便由此严格定义出一个新的状态函数——熵 (SS)。其微元定义式为:

    dS=δQrevTdS = \frac{\delta Q_{\text{rev}}}{T}

    熵的出现,使得过程量热量 δQ\delta Q 能够通过积分因子 1T\frac{1}{T} 转化为全微分 dSdS

1.4.2 热力学基本方程#

  • 将热力学第一定律 δQ=dU+δW\delta Q = dU + \delta W 与熵的微元定义相结合。对于一个只做可逆体积功的简单可压缩系统,有 δW=pdV\delta W = p dVδQ=TdS\delta Q = T dS。代入第一定律,即得到热力学基本微分方程

    dU=TdSpdVdU = T dS - p dV

    该方程虽然是借助“可逆过程”的条件(δQ=TdS\delta Q = T dSδW=pdV\delta W = p dV)推导出来的,但由于方程中包含的物理量(U,T,S,p,VU, T, S, p, V全部都是状态函数,状态函数的微小变化仅取决于相邻的初末平衡态。因此,该基本方程适用于简单系统相邻两个平衡态之间的任何过程(包括不可逆过程)。它是连接热力学第一定律与第二定律的核心枢纽。

    为了将抽象的熵函数具体化,我们以理想气体为例,推导其宏观熵的解析表达式。这不仅是热力学计算的基础,也能直观地展示态函数如何通过状态参量来完全确定。

    1. 以 (T,V)(T, V) 为独立变量的推导

    从热力学基本微分方程出发:

    TdS=dU+pdV    dS=dUT+pTdVT dS = dU + p dV \implies dS = \frac{dU}{T} + \frac{p}{T} dV

    将理想气体的两个核心属性代入该式:

    • 焦耳定律:内能仅为温度的函数,dU=nCVdTdU = nC_V dT

    • 物态方程pV=nRT    pT=nRVpV = nRT \implies \frac{p}{T} = \frac{nR}{V}

    代入后得到熵的微元微分式:

    dS=nCVdTT+nRdVVdS = nC_V \frac{dT}{T} + nR \frac{dV}{V}

    假定在所讨论的温度范围内,定容热容 CVC_V 为常数,对等式两边直接积分,得到以 T,VT, V 为自变量的熵函数表达式:

    S(T,V)=nCVlnT+nRlnV+S0S(T, V) = nC_V \ln T + nR \ln V + S_0

    2. 以 (T,p)(T, p) 为独立变量的推导

    在许多实际过程中(如等压过程),选取温度和压强作为独立变量更为方便。我们借助焓的微元定义 dH=dU+pdV+VdpdH = dU + p dV + V dp,将基本方程改写为:

    TdS=dHVdp    dS=dHTVTdpT dS = dH - V dp \implies dS = \frac{dH}{T} - \frac{V}{T} dp

    同理,代入理想气体的焓变性质(dH=nCpdTdH = nC_p dT)与物态方程(VT=nRp\frac{V}{T} = \frac{nR}{p}):

    dS=nCpdTTnRdppdS = nC_p \frac{dT}{T} - nR \frac{dp}{p}

    在定压热容 CpC_p 为常数的假设下,积分可得:

    S(T,p)=nCplnTnRlnp+S0S(T, p) = nC_p \ln T - nR \ln p + S'_0

    式中的 S0S_0S0S'_0 是积分常数。在经典宏观热力学框架内,熵只能确定到相差一个任意常数,我们只能计算过程前后的熵变 (ΔS\Delta S)。熵的绝对数值零点,必须等到引入热力学第三定律和量子统计物理后才能严格标定。上述简洁的对数表达式严格依赖于“热容是常数”这一近似。如果温度变化跨度极大(激发了分子的振动自由度),热容将成为温度的函数 CV(T)C_V(T),此时必须退回保留积分号的原始形式:S=nCV(T)TdT+nRlnV+S0S = \int \frac{nC_V(T)}{T} dT + nR \ln V + S_0

    S(T,V)S(T, V) 表达式中,体积项前是 正号 (++)。物理图像:体积膨胀,分子活动空间增大,系统微观状态数增多(混乱度增加),导致熵增。在 S(T,p)S(T, p) 表达式中,压强项前是 负号 (-)。物理图像:压强增大意味着系统被压缩,分子活动空间受限,混乱度降低,导致熵减。

    由于熵是态函数,不论理想气体经历的是可逆过程还是极度复杂的不可逆过程(如绝热自由膨胀),计算其初末态的熵变 ΔS\Delta S 时,只需将初态和末态的参量代入上述公式直接相减即可,完全不需要考虑中间的实际路径。

1.4.3 克劳修斯不等式与熵增加原理#

  • 对于包含不可逆过程的实际循环,根据卡诺定理,其效率必定小于相同热源间的可逆卡诺循环。经过严格的数学推演,可将克劳修斯等式推广为包含不可逆过程的克劳修斯不等式

    δQT0\oint \frac{\delta Q}{T} \le 0

    (等号适用于全可逆循环,小于号适用于存在不可逆性的循环)。

    若系统经任意过程(可逆或不可逆)从初态 A 演化至终态 B,结合熵的定义,可由克劳修斯不等式导出:

    SBSAABδQTS_B - S_A \ge \int_A^B \frac{\delta Q}{T}

    其微元形式为:

    dSδQTdS \ge \frac{\delta Q}{T}

    熵增加原理

    将上述微元关系应用于孤立系统(系统与外界无物质和能量交换,即必定有 δQ=0\delta Q = 0),不等式直接化简为:

    dS0dS \ge 0

    孤立系统中发生的一切自发过程,必然朝着熵增加的方向进行;当系统的熵达到特定宏观约束下的极大值时,过程停止,系统进入稳定的热力学平衡态。熵增加原理为自然界宏观演化的方向性和限度,提供了绝对且严密的数学判据。

1.5 辅助特性函数的初步引入#

  • 熵增加原理(dS0dS \ge 0)为宏观过程的方向性提供了绝对判据,但其严格的适用前提是孤立系统(即系统的内能 UU 和体积 VV 保持不变)。在实际的物理和化学过程中,系统往往并非孤立,而是处于等温、等容或等压等特定约束条件下。为了在这些常见约束下寻找系统平衡的判据,需要引入新的热力学状态函数。

1.5.1 自由能 (FF) 与等温等容系统的判据#

  • 考察一个与恒温热源(温度为 TT)接触的系统。根据克劳修斯不等式,系统在微小过程中吸收的热量 δQ\delta Q 满足:

    δQTdS\delta Q \le T dS

    结合热力学第一定律 δQ=dU+δW\delta Q = dU + \delta W(其中 δW\delta W 为系统对外做功),可得:

    dU+δWTdSdU + \delta W \le T dS

    由于热源温度 TT 恒定(等温过程),TdST dS 可改写为 d(TS)d(TS)。移项整理得到:

    d(UTS)δWd(U - TS) \le -\delta W

    为此,定义一个新的状态函数——自由能(Helmholtz Free Energy, FF

    FUTSF \equiv U - TS

    代入上式,得到等温过程的核心不等式:

    dFδWdF \le -\delta W

    在等温过程中,系统对外所做的功,永远不可能超过其自由能的减少量。自由能的减少量是在该条件下系统做功的理论上限(可逆过程取等号)。自由能 FF 即为内能中能够“自由”转化为对外做功的那部分能量。

    等温等容判据:若系统进一步被限制在等容条件下(dV=0dV = 0),且无其他形式的功,则 δW=0\delta W = 0。上述不等式化简为:

    dF0dF \le 0

    在等温等容条件下,系统内部发生的一切不可逆过程,必然朝着自由能减少的方向进行。当自由能达到极小值时,系统进入稳定的平衡态。

1.5.2 吉布斯函数 (GG) 与等温等压系统的判据#

  • 实验室中最常见的条件是系统同时与恒温热源(温度 TT)和恒压环境(大气压 pp)接触。在此条件下,系统对外的体积功为 δW=pdV\delta W = p dV。代入等温过程的不等式 d(UTS)δWd(U - TS) \le -\delta W,有:

    d(UTS)pdVd(U - TS) \le -p dV

    由于环境压强 pp 恒定,pdVp dV 可改写为 d(pV)d(pV)。移项合并所有全微分项,得到:

    d(UTS+pV)0d(U - TS + pV) \le 0

    为此,定义另一个新的状态函数——吉布斯函数(Gibbs Free Energy, GG

    GUTS+pV=HTSG \equiv U - TS + pV = H - TS

    代入上式,得到等温等压过程的核心不等式:

    dG0dG \le 0

    等温等压判据:在等温等压条件下,系统自发进行的过程必然导致吉布斯函数减小;当吉布斯函数达到极小值时,系统处于热力学平衡态。

1.5.3 判据的绝对对应关系#

  • 孤立系统(恒 U,VU, V\rightarrow SS 极大

    等温等容系统(恒 T,VT, V\rightarrow 自由能 FF 极小

    等温等压系统(恒 T,pT, p\rightarrow 吉布斯函数 GG 极小

    切忌在条件不匹配时乱用判据。例如,在绝热膨胀过程中强行使用吉布斯函数极小值作为判据是严重的物理错误。判据的有效性严格绑定于其推导时所引入的外部约束。

  

  

  


第2章:热力学函数与应用#

  • 在第一章中,我们确立了热力学基本方程 dU=TdSpdVdU = TdS - pdV。但该方程在实际应用中存在一个缺陷:它包含了内能 UU 和熵 SS 这样无法在实验中直接测量的“绝对量”。为了将理论推导与实验观测(如压强 pp、体积 VV、温度 TT 及热容)对接,我们必须引入一套严格的宏观数学转换工具。

2.1 热力学基本函数与特性的数学推演#

2.1.1 勒让德变换与四个基本热力学函数#

  • 热力学基本方程 dU=TdSpdVdU = TdS - pdV 表明,内能 UU 的“自然变量”是 (S,V)(S, V)。但在实验室中,熵 SS 是极其难以控制的变量,我们更习惯控制温度 TT 和压强 pp

    为了实现独立变量的转换而不丢失系统的热力学信息,数学上采用勒让德变换(Legendre Transformation)。其核心思想是:若要把函数 f(x)f(x) 的自变量转换为其导数 y=fxy = \frac{\partial f}{\partial x},只需构造新函数 g=fxyg = f - xy 即可。利用该变换,我们将包含热量交换的项 TSTS 与包含做功的项 pVpV 进行组合,严格定义出宏观热力学的四个基本函数(及其全微分式):

    1. 内能 (Internal Energy) U(S,V)U(S, V)

      dU=TdSpdVdU = TdS - pdV
    2. 焓 (Enthalpy) H(S,p)U+pVH(S, p) \equiv U + pV

      dH=TdS+VdpdH = TdS + Vdp
    3. 自由能 (Helmholtz Free Energy) F(T,V)UTSF(T, V) \equiv U - TS

      dF=SdTpdVdF = -SdT - pdV
    4. 吉布斯函数 (Gibbs Free Energy) G(T,p)UTS+pVG(T, p) \equiv U - TS + pV

      dG=SdT+VdpdG = -SdT + Vdp

2.1.2 麦克斯韦关系#

  • 四个基本函数的全微分不仅确立了自然变量,还直接导出了热力学中最强大的数学等式。对于任意状态函数(全微分)dz=Mdx+Ndydz = M dx + N dy,数学上必然满足二阶混合偏导数相等:(My)x=(Nx)y\left(\frac{\partial M}{\partial y}\right)_x = \left(\frac{\partial N}{\partial x}\right)_y。将此纯数学定理直接应用于上述四个基本函数的微分表达式,立即得到四大麦克斯韦关系式

    1. dUdU(TV)S=(pS)V\left(\frac{\partial T}{\partial V}\right)_S = -\left(\frac{\partial p}{\partial S}\right)_V

    2. dHdH(Tp)S=(VS)p\left(\frac{\partial T}{\partial p}\right)_S = \left(\frac{\partial V}{\partial S}\right)_p

    3. dFdF(SV)T=(pT)V\left(\frac{\partial S}{\partial V}\right)_T = \left(\frac{\partial p}{\partial T}\right)_V

    4. dGdG(Sp)T=(VT)p\left(\frac{\partial S}{\partial p}\right)_T = -\left(\frac{\partial V}{\partial T}\right)_p

    麦克斯韦关系是热力学理论走向实验的绝对桥梁。特别是第3和第4个关系式,它们成功地将包含熵 SS 的不可测偏导数,完美地转化为了由 p,V,Tp, V, T 构成的完全可测偏导数。

  • 雅可比行列式恒等式在热力学中,雅可比行列式有一个极为方便的性质,它可以将任何偏导数转化为行列式的商:

    (ux)y=(u,y)(x,y)\left( \frac{\partial u}{\partial x} \right)_y = \frac{\partial(u, y)}{\partial(x, y)}

    利用这一工具,我们可以得到一个极其重要的恒等式:

    (T,S)(p,V)=1\frac{\partial(T, S)}{\partial(p, V)} = 1

    证明:

    直接取由自由能 dFdF 导出的麦克斯韦关系:

    (SV)T=(pT)V\left(\frac{\partial S}{\partial V}\right)_T = \left(\frac{\partial p}{\partial T}\right)_V

    利用雅可比行列式的链式乘法法则,将等式两边同乘 (T,V)(p,V)\frac{\partial(T, V)}{\partial(p, V)},即得证:

    (T,S)(p,V)=1\frac{\partial(T, S)}{\partial(p, V)} = 1

    该恒等式表明,从温熵坐标 (T,S)(T, S) 变换到压容坐标 (p,V)(p, V) 时,其雅可比行列式(即变换的模)严格等于 1。这不仅意味着在平衡态热力学的可逆过程中相空间的“面积”是守恒的,更证明了 (T,S)(T, S)(p,V)(p, V) 作为两对共轭变量,在热力学结构中具有完全对等的物理地位。

2.1.3 特性函数及其应用#

  • 如果把某一个热力学函数表示为特定的一组独立变量的函数,使得不用借助于任何其他的函数(如由实验测定的物态方程或热容量),只要对该函数进行偏导数运算,就可以求得均匀系的全部热力学函数,从而把系统的全部热力学性质完全确定下来,那么这个函数就被称为在该组特定独立变量下的特性函数。常见的热力学特性函数有:

    内能的特性函数:U(S,V)U(S, V)

    焓的特性函数:H(S,p)H(S, p)

    自由能的特性函数:F(T,V)F(T, V)

    吉布斯函数的特性函数:G(T,p)G(T, p)

    1.以 (T,V)(T, V) 为独立变量的表述(自由能表示)

    在等温过程中,最方便的特性函数是自由能 F(T,V)F(T, V)。一旦获知 F(T,V)F(T, V) 的解析式,系统的所有性质可按如下逻辑导出:

    基本参数:

    根据 dF=SdTpdVdF = -SdT - pdV,直接求一阶偏导:

    熵:S=(FT)VS = -\left( \frac{\partial F}{\partial T} \right)_V

    压强:p=(FV)Tp = -\left( \frac{\partial F}{\partial V} \right)_T (此即系统的状态方程)

    能量参数:

    内能:U=F+TS=FT(FT)VU = F + TS = F - T\left( \frac{\partial F}{\partial T} \right)_V

    焓:H=U+pV=FT(FT)VV(FV)TH = U + pV = F - T\left( \frac{\partial F}{\partial T} \right)_V - V\left( \frac{\partial F}{\partial V} \right)_T

    特性参量:

    吉布斯函数:G=F+pV=FV(FV)TG = F + pV = F - V\left( \frac{\partial F}{\partial V} \right)_T

    定容热容:CV=(UT)V=T(2FT2)VC_V = \left( \frac{\partial U}{\partial T} \right)_V = -T\left( \frac{\partial^2 F}{\partial T^2} \right)_V

    2.以 (U,V)(U, V) 为独立变量的表述(熵表示)

    根据热力学基本方程的另一种形式 dS=1TdU+pTdVdS = \frac{1}{T}dU + \frac{p}{T}dV,可知在以 (U,V)(U, V) 为独立变量时,S(U,V)S(U, V) 是该系统的特性函数(属于熵表示,Entropy Representation)。

    已知 S(U,V)S(U, V) 后,推导逻辑如下:

    基本参数:

    温度:根据 1T=(SU)V\frac{1}{T} = \left( \frac{\partial S}{\partial U} \right)_V,可得 T=[(SU)V]1T = \left[ \left( \frac{\partial S}{\partial U} \right)_V \right]^{-1}

    压强:根据 pT=(SV)U\frac{p}{T} = \left( \frac{\partial S}{\partial V} \right)_U,可得 p=T(SV)U=(S/V)U(S/U)Vp = T \left( \frac{\partial S}{\partial V} \right)_U = \frac{(\partial S / \partial V)_U}{(\partial S / \partial U)_V}

    其他热力学函数:

    利用已求得的 TTpp,可以直接代入定义式求出其他态函数:

    F=UTSF = U - TS

    H=U+pVH = U + pV

    G=UTS+pVG = U - TS + pV

  • 特性函数的“全能性”严格依赖于变量的选择。例如,虽然 UU 是状态函数,但如果将其表达为 U(T,V)U(T, V),它就失去了特性函数的地位。因为从 U(T,V)U(T, V) 出发,我们只能得到热容信息,却无法唯一地推导出系统的状态方程(压强 pp)。

    假设已知系统的内能函数为 U(T,V)U(T, V),试着去求该系统的状态方程 p(T,V)p(T, V)

    借助能态方程(由麦克斯韦关系导出):

    (UV)T=T(pT)Vp\left(\frac{\partial U}{\partial V}\right)_T = T\left(\frac{\partial p}{\partial T}\right)_V - p

    将上式整理为关于未知量 pp 的偏微分方程:

    (pT)V1Tp=1T(UV)T\left(\frac{\partial p}{\partial T}\right)_V - \frac{1}{T}p = \frac{1}{T}\left(\frac{\partial U}{\partial V}\right)_T

    等式两边同乘积分因子 1T2\frac{1}{T^2},左边凑成全微分形式:

    [T(pT)]V=1T2(UV)T\left[ \frac{\partial}{\partial T} \left( \frac{p}{T} \right) \right]_V = \frac{1}{T^2}\left(\frac{\partial U}{\partial V}\right)_T

    在保持体积 VV 不变的条件下,对温度 TT 进行偏积分:

    pT=1T2(UV)TdT+f(V)\frac{p}{T} = \int \frac{1}{T^2}\left(\frac{\partial U}{\partial V}\right)_T dT + f(V)

    得到压强的最终表达式:

    p(T,V)=T1T2(UV)TdT+Tf(V)p(T, V) = T \int \frac{1}{T^2}\left(\frac{\partial U}{\partial V}\right)_T dT + T \cdot f(V)

    由于我们进行的是偏积分,产生的积分常数不再是一个具体的数值,而是一个完全未知的、关于体积 VV 的任意函数 f(V)f(V)。这表明,即使我们完全掌握了 U(T,V)U(T, V) 的解析式,也无法唯一确定压强 pp 的表达式(即无法给出明确的物态方程)。因此,U(T,V)U(T, V) 无法仅通过偏导数就给出系统的全部信息,它失去了特性函数的地位

2.1.4 广义热容关系#

  • 广义功与变量的物理对应

    为了展现热力学框架的极致普适性,我们推导任意共轭参量 (x,y)(x, y) 下的热容关系。在普遍的热力学系统中,系统在微小过程中的“做功”项(δW\delta W)总可以抽象为两个共轭变量的乘积:

    δW=ydx\delta W = y dx

    xx(广义位移):通常代表系统的某种广延量(如体积 VV、总电偶极矩 pep_e、总磁矩 mm)。

    yy(广义力):通常代表促使系统状态改变的强度量(如压强 pp、外加电场 EE、外加磁场 HH)。

    (注:以最常见的流体膨胀功 pdVp dV 为例,我们通常直接映射 xVx \leftrightarrow Vypy \leftrightarrow p)

    基于此,设系统保持广义位移 xx 恒定时的热容为 CxC_x,保持广义力 yy 恒定时的热容为 CyC_y。其热力学严格定义为:

    Cx=T(ST)x,Cy=T(ST)yC_x = T\left(\frac{\partial S}{\partial T}\right)_x, \quad C_y = T\left(\frac{\partial S}{\partial T}\right)_y

    1. 广义热容之差 (CyCxC_y - C_x)

    选取温度 TT 和广义位移 xx 作为独立变量,将熵写为 S(T,x)S(T, x),取全微分:

    dS=(ST)xdT+(Sx)TdxdS = \left(\frac{\partial S}{\partial T}\right)_x dT + \left(\frac{\partial S}{\partial x}\right)_T dx

    等式两边同除以 dTdT,并在保持广义力 yy 恒定的条件下取偏导数:

    (ST)y=(ST)x+(Sx)T(xT)y\left(\frac{\partial S}{\partial T}\right)_y = \left(\frac{\partial S}{\partial T}\right)_x + \left(\frac{\partial S}{\partial x}\right)_T \left(\frac{\partial x}{\partial T}\right)_y

    等式两边同乘温度 TT,代入热容定义,直接得到广义热容差的普遍公式

    CyCx=T(Sx)T(xT)yC_y - C_x = T \left(\frac{\partial S}{\partial x}\right)_T \left(\frac{\partial x}{\partial T}\right)_y

    V,pV, p 代入上式,并调用麦克斯韦关系 (SV)T=(pT)V\left(\frac{\partial S}{\partial V}\right)_T = \left(\frac{\partial p}{\partial T}\right)_V,公式立刻退化为我们熟知的形式:

    CpCV=T(pT)V(VT)p=TVα2κTC_p - C_V = T \left(\frac{\partial p}{\partial T}\right)_V \left(\frac{\partial V}{\partial T}\right)_p = \frac{T V \alpha^2}{\kappa_T}

    CyC_y 之所以不等于 CxC_x,是因为在等广义力 yy(如等压)升温时,系统必须改变其广义位移 xx(如体积膨胀)。这一变化不仅涉及克服内部相互作用势能的改变,还需要对外输出广义功(δW=ydx\delta W = y dx),从而必须吸收额外的热量。

    2. 广义热容之比 (Cy/CxC_y / C_x)

    利用微积分中的循环关系定理 (AB)C(BC)A(CA)B=1\left(\frac{\partial A}{\partial B}\right)_C \left(\frac{\partial B}{\partial C}\right)_A \left(\frac{\partial C}{\partial A}\right)_B = -1,将偏导数展开:

    对于 (S,T,y)(S, T, y) 系统: (ST)y=(Sy)T(yT)S\left(\frac{\partial S}{\partial T}\right)_y = - \left(\frac{\partial S}{\partial y}\right)_T \left(\frac{\partial y}{\partial T}\right)_S

    对于 (S,T,x)(S, T, x) 系统: (ST)x=(Sx)T(xT)S\left(\frac{\partial S}{\partial T}\right)_x = - \left(\frac{\partial S}{\partial x}\right)_T \left(\frac{\partial x}{\partial T}\right)_S

    将两者相比,并利用链式法则化简 (S/y)T(S/x)T=(xy)T\frac{(\partial S / \partial y)_T}{(\partial S / \partial x)_T} = \left(\frac{\partial x}{\partial y}\right)_T 以及 (y/T)S(x/T)S=(yx)S\frac{(\partial y / \partial T)_S}{(\partial x / \partial T)_S} = \left(\frac{\partial y}{\partial x}\right)_S

    CyCx=(S/T)y(S/T)x=(xy)T(yx)S\frac{C_y}{C_x} = \frac{(\partial S / \partial T)_y}{(\partial S / \partial T)_x} = \left(\frac{\partial x}{\partial y}\right)_T \left(\frac{\partial y}{\partial x}\right)_S

    整理得到广义热容比的普遍公式

    CyCx=(yx)S(yx)T\frac{C_y}{C_x} = \frac{\left(\frac{\partial y}{\partial x}\right)_S}{\left(\frac{\partial y}{\partial x}\right)_T}

    V,pV, p 代入上式:

    CpCV=(pV)S(pV)T\frac{C_p}{C_V} = \frac{\left(\frac{\partial p}{\partial V}\right)_S}{\left(\frac{\partial p}{\partial V}\right)_T}

    由于等温压缩系数 κT=1V(Vp)T\kappa_T = -\frac{1}{V}\left(\frac{\partial V}{\partial p}\right)_T,绝热压缩系数 κS=1V(Vp)S\kappa_S = -\frac{1}{V}\left(\frac{\partial V}{\partial p}\right)_S,该比值直接映射为力学参量之比:

    CpCV=κTκS\frac{C_p}{C_V} = \frac{\kappa_T}{\kappa_S}

    热学响应之比(热容比)与力学响应之比(压缩系数比)存在着完美的对应关系。由于必然有 Cp>CVC_p > C_V,因此在相同的体积压缩量下,绝热压缩引起的压强攀升一定比等温压缩更剧烈(κT>κS\kappa_T > \kappa_S)。

2.2 获取低温的物理过程与原理#

2.2.1 节流膨胀与焦耳-汤姆孙效应(等焓制冷)#

  • 气体在持续的压力差下,通过多孔塞(或节流阀)从高压侧 (p1p_1) 缓慢渗漏到低压侧 (p2p_2) 的过程,称为节流膨胀。假设多孔塞及管道完全绝热(Q=0Q = 0)。

    分析特定质量气体在通过节流阀前后的能量转化:

    外界对气体做推挤功:W1=p1V1W_1 = p_1 V_1

    气体膨胀对外界做功:W2=p2V2W_2 = -p_2 V_2

    根据第一定律 ΔU=W1+W2\Delta U = W_1 + W_2,有 U2U1=p1V1p2V2U_2 - U_1 = p_1 V_1 - p_2 V_2

    移项整理得:U1+p1V1=U2+p2V2U_1 + p_1 V_1 = U_2 + p_2 V_2

    这证明了节流膨胀前后的初末态满足 H1=H2H_1 = H_2。在热力学处理中,该过程被抽象为一个等焓过程

  • 为了定量评估节流后的温度变化,定义焦耳-汤姆孙系数 μJT=(Tp)H\mu_{JT} = \left(\frac{\partial T}{\partial p}\right)_H。我们利用上一节的雅可比行列式性质进行极简推导:

    μJT=(Tp)H=(T,H)(p,H)=(T,H)(T,p)(p,H)(T,p)=(Hp)T(HT)p\mu_{JT} = \left(\frac{\partial T}{\partial p}\right)_H = \frac{\partial(T, H)}{\partial(p, H)} = \frac{\frac{\partial(T, H)}{\partial(T, p)}}{\frac{\partial(p, H)}{\partial(T, p)}} = -\frac{\left(\frac{\partial H}{\partial p}\right)_T}{\left(\frac{\partial H}{\partial T}\right)_p}

    分母即为定压热容 CpC_p。对于分子,将基本方程 dH=TdS+VdpdH = TdS + Vdp 在等温条件下对压强求导:

    (Hp)T=T(Sp)T+V\left(\frac{\partial H}{\partial p}\right)_T = T\left(\frac{\partial S}{\partial p}\right)_T + V

    调用由吉布斯函数 dGdG 导出的麦克斯韦关系 (Sp)T=(VT)p\left(\frac{\partial S}{\partial p}\right)_T = -\left(\frac{\partial V}{\partial T}\right)_p

    (Hp)T=T(VT)p+V\left(\frac{\partial H}{\partial p}\right)_T = -T\left(\frac{\partial V}{\partial T}\right)_p + V

    代回原式,得到焦耳-汤姆孙系数的普适解析式:

    μJT=1Cp[T(VT)pV]=VCp(Tα1)\mu_{JT} = \frac{1}{C_p} \left[ T \left(\frac{\partial V}{\partial T}\right)_p - V \right] = \frac{V}{C_p} (T\alpha - 1)

    (其中 α=1V(VT)p\alpha = \frac{1}{V}\left(\frac{\partial V}{\partial T}\right)_p 为体胀系数)。

    1. 理想气体无效应:对于理想气体 V=nRT/pV = nRT/p,必然有 Tα=1T\alpha = 1,导致 μJT=0\mu_{JT} = 0。理想气体内能仅依赖于温度,等焓膨胀既不改变内能,也不改变温度。

    2. 真实气体的制冷与转化:真实气体的 μJT\mu_{JT} 受分子间相互作用主导。

      μJT>0\mu_{JT} > 0 时,压强下降 (dp<0dp<0) 必然导致温度下降 (dT<0dT<0),实现制冷。(微观机制:分子间距变大,克服分子间吸引力做功消耗了分子动能)。

      μJT<0\mu_{JT} < 0 时,压强下降导致温度上升,即致热。(微观机制:气压极高时斥力占优,膨胀时斥力做正功转化为动能)。

      μJT=0\mu_{JT} = 0 得到的 TpT-p 关系即为转化曲线。在工程液化中,必须先将气体预冷至转化曲线内部,方可利用节流阀实现持续降温。

2.2.2 绝热可逆膨胀(等熵制冷)#

  • 若气体在绝热气缸中缓慢推动无摩擦活塞对外做功,过程严格可逆且绝热(δQrev=0\delta Q_{\text{rev}} = 0),这是一个等熵过程

    定义等熵膨胀降温系数 μS=(Tp)S\mu_S = \left(\frac{\partial T}{\partial p}\right)_S。同样利用雅可比行列式:

    μS=(T,S)(p,S)=(T,S)(T,p)(p,S)(T,p)=(Sp)T(ST)p\mu_S = \frac{\partial(T, S)}{\partial(p, S)} = \frac{\frac{\partial(T, S)}{\partial(T, p)}}{\frac{\partial(p, S)}{\partial(T, p)}} = -\frac{\left(\frac{\partial S}{\partial p}\right)_T}{\left(\frac{\partial S}{\partial T}\right)_p}

    分母为 CpT\frac{C_p}{T}。分子再次调用麦克斯韦关系 (Sp)T=(VT)p=Vα\left(\frac{\partial S}{\partial p}\right)_T = -\left(\frac{\partial V}{\partial T}\right)_p = -V\alpha

    μS=TVαCp\mu_S = \frac{T V \alpha}{C_p}
  • 比较两种制冷机制的系数差值:

    μSμJT=TVαCpVCp(Tα1)=VCp>0\mu_S - \mu_{JT} = \frac{TV\alpha}{C_p} - \frac{V}{C_p}(T\alpha - 1) = \frac{V}{C_p} > 0

    由于 μS>μJT\mu_S > \mu_{JT},在经历相同的压强降幅(Δp<0\Delta p < 0)时,绝热可逆膨胀的降温幅度永远大于节流膨胀。从能量守恒的角度看这是必然的:节流膨胀气体仅克服内聚力做功;而绝热膨胀气体不仅要克服内聚力,还要向外输出宏观机械功,导致其内能消耗更为剧烈。

2.2.3 磁制冷:顺磁体的绝热去磁#

  • 气体液化技术在 1K 左右遭遇物理极限(氦气的蒸气压变得极低)。要进入 mK(毫开尔文)级别的极低温领域,体积膨胀已失效,物理学转而利用固体内部的电子磁矩系统

    设系统总磁矩为 mm,外加磁场强度为 HH。磁场对顺磁体系统做功为 δW=Hdm\delta W = H dm(注:当系统磁矩增加时,外界磁场对其做正功)。结合第一定律 δQ=dUδW\delta Q = dU - \delta W,热力学基本方程改写为:

    dU=TdS+HdmdU = TdS + Hdm

    为了利用麦克斯韦关系,定义磁介质的“吉布斯函数” G=UTSHmG = U - TS - Hm。取全微分并代入基本方程:

    dG=SdTmdHdG = -SdT - mdH

    由此立刻得到对应的麦克斯韦关系:(SH)T=(mT)H\left(\frac{\partial S}{\partial H}\right)_T = \left(\frac{\partial m}{\partial T}\right)_H

    • 绝热去磁的核心是考察等熵条件下减弱磁场引起的温度变化,即求 (TH)S\left(\frac{\partial T}{\partial H}\right)_S

    利用雅可比行列式:

    (TH)S=(SH)T(ST)H=TCH(mT)H\left(\frac{\partial T}{\partial H}\right)_S = -\frac{\left(\frac{\partial S}{\partial H}\right)_T}{\left(\frac{\partial S}{\partial T}\right)_H} = -\frac{T}{C_H} \left(\frac{\partial m}{\partial T}\right)_H

    根据顺磁体的居里定律,在恒定磁场下,温度升高导致热运动加剧,磁化强度必然减弱,即 (mT)H<0\left(\frac{\partial m}{\partial T}\right)_H < 0

    代入上式,得出绝对结论:(TH)S>0\left(\frac{\partial T}{\partial H}\right)_S > 0。绝热状态下,减小磁场(dH<0dH < 0),必然导致温度下降(dT<0dT < 0)。

2.3 热力学第三定律与负温度#

2.3.1 能斯特定理与绝对零度不可达#

  • 能斯特在总结大量低温化学反应和相变实验后发现,随着温度趋近于绝对零度,系统不仅熵的绝对值在减小,等温过程中的熵变(ΔS\Delta S)也严格趋于零。其热力学数学表述为:对于任何处于内部平衡态的凝聚系,当温度趋于绝对零度时,系统经历任何可逆等温过程的熵变趋于零。即:

    limT0ΔS=0\lim_{T \to 0} \Delta S = 0

    这就是热力学第三定律的能斯特表述。普朗克随后将其进一步强化为绝对熵的标定:当 T0T \to 0 时,所有凝聚系的绝对熵趋于零(limT0S=0\lim_{T \to 0} S = 0)。在微观统计物理中,这对应着系统基态的非简并性(即微观状态数 Ω=1\Omega = 1S=kBlnΩ=0S = k_B \ln \Omega = 0)。

  • 能斯特定理在逻辑上直接宣判了绝对零度是宏观热力学过程的禁区。

    回顾上一节的磁制冷循环:制冷依赖于“高场等温磁化(熵减小)”和“绝热去磁(等熵降温)”的交替。如果能斯特定理成立,当系统逼近 0 K 时,任意两个不同外场参数(如 H1H_1H2H_2)下的熵-温度曲线(STS-T 曲线)在 T=0T=0 处必然收敛于同一点(因为等温熵变 ΔS0\Delta S \to 0)。无论在何种高场下开始绝热退磁(沿水平等熵线向左移动),由于 STS-T 曲线在底部必定汇聚并变平,等熵线永远无法在有限次循环内与 T=0T=0 的纵轴相交

    这就引出了第三定律的等价表述:不可能通过有限的手续将任何系统冷却到绝对零度。

  • 第三定律不仅限制了制冷机的下限,还给所有物质的低温热力学响应系数下了数学铁律。结合第二章的麦克斯韦关系,当 T0T \to 0 时:

    1. 热容趋于零:由于热容定义为 Cx=T(S/T)xC_x = T(\partial S / \partial T)_x。在 T0T \to 0 时,熵的绝对值趋于常数(0),其对温度的导数不发散,因此必然有 limT0CV=0\lim_{T \to 0} C_V = 0limT0Cp=0\lim_{T \to 0} C_p = 0。这意味着在极低温下,极其微小的热量扰动就能引起巨大的温度剧变。

    2. 体胀系数趋于零:由麦克斯韦关系 (S/p)T=Vα(\partial S/\partial p)_T = -V\alpha。因为 T0T \to 0 时等温熵变 ΔS0\Delta S \to 0,即熵对压强完全不敏感((S/p)T0(\partial S/\partial p)_T \to 0),必然有 limT0α=0\lim_{T \to 0} \alpha = 0。物质在绝对零度附近将完全丧失热胀冷缩的力学性质。

2.3.2 负温度的统计物理本质#

  • 在热力学中,温度的最严密定义是 1T=(SU)V,N\frac{1}{T} = \left(\frac{\partial S}{\partial U}\right)_{V, N}。要使温度 TT 出现负值,也就是说,随着系统总能量 EE 的增加,系统可占据的微观状态数必须减少。负温度系统必须满足一个严苛的先决条件:系统的总能量必须有严格的上限,否则能量越高系统越混乱,无法达到负温度状态。

    对于处于热力学平衡的系统,处于能量为 εi\varepsilon_i 的微观状态上的粒子数 nin_i,服从玻尔兹曼分布律

    niexp(εikBT)n_i \propto \exp\left(-\frac{\varepsilon_i}{k_B T}\right)

    或者考察高能级(ε2\varepsilon_2)与低能级(ε1\varepsilon_1)的粒子数之比:

    n2n1=exp(ε2ε1kBT)\frac{n_2}{n_1} = \exp\left(-\frac{\varepsilon_2 - \varepsilon_1}{k_B T}\right)

    这个简单的指数公式,决定了不同温度下达到平衡时系统的总内能(U=niεiU = \sum n_i \varepsilon_i)大小。

    设想一个总能量受限的双能级自旋系统(低能级 ε1\varepsilon_1,高能级 ε2\varepsilon_2):

    1. 常规状态(正温度 T>0T > 0 且能量较低):

      由于 T>0T > 0,指数项必定小于 1。因此 n2<n1n_2 < n_1高能级上的粒子数永远少于低能级。系统极其“守规矩”,绝大多数粒子沉寂在低能态。此时总内能 UU 处于较低水平。如果系统继续吸收热量,粒子逐渐向高能级跃迁,混乱度(熵 SS)增加,表现为 TT 逐渐升高。

    2. 极限混乱状态(T±T \to \pm\infty 且能量居中):

      当系统吸收了足够多的能量,使得 n2=n1n_2 = n_1(高低能级粒子数完全相等)时,系统达到了绝对的最混乱状态,熵 SS 达到抛物线的顶点。

      代入玻尔兹曼公式:n2n1=1    exp(ΔεkBT)=1\frac{n_2}{n_1} = 1 \implies \exp\left(-\frac{\Delta \varepsilon}{k_B T}\right) = 1。这要求分母 T±T \to \pm\infty

    3. 粒子数反转与负温度(T<0T < 0 且能量极高):

      如果系统继续被强行注入能量,高能级上的粒子数将超过低能级,即 n2>n1n_2 > n_1。这在物理上被称为粒子数反转(Population Inversion)(这也是激光产生的原理)。

      回到玻尔兹曼公式,由于 n2n1>1\frac{n_2}{n_1} > 1,指数必然大于 0。由于 Δε\Delta \varepsilonkBk_B 都是正数,唯一的数学可能是:温度 TT 必须变为负值

      此时,随着能量 UU 继续增加(更多粒子跑到高能级),系统反而变得越来越有序,混乱度 SS 开始下降。这就完美呼应了宏观的热力学定义 1T=SU<0\frac{1}{T} = \frac{\partial S}{\partial U} < 0

    通过玻尔兹曼分布,我们可以极其直观地看到:

    1. 任何正温度系统,无论温度多高,其高能级的粒子数总是小于低能级粒子数。

    2. 任何负温度系统,其微观前提就是粒子数反转,即高能级粒子数大于低能级粒子数。

    由于内能是所有粒子能量的总和(U=niεiU = \sum n_i \varepsilon_i),一个大量粒子占据高能级的系统,其总内能绝对、严格地大于任何由低能级主导的正常系统。

    既然负温度系统蕴含极高的内能,当它与普通正温度系统接触时,负温度系统必然充当“高温热源”向外释放热量(这也满足第二定律使联合系统总熵增加的要求)。因此,负温度并非比绝对零度更冷,而是跨越了正无穷后,比任何正温度都更热的极端状态。真实的宏观温度高低演化排序(从极冷到极热、内能从低到高)应严格表述为:

    +0 K+300 K+ K K300 K0 K+0 \text{ K} \rightarrow +300 \text{ K} \rightarrow +\infty \text{ K} \equiv -\infty \text{ K} \rightarrow -300 \text{ K} \rightarrow -0 \text{ K}

  

  

  


第3章 相变#

  • 在实际的热力学过程中,物质经常以多个相(如固、液、气相)共存的形式出现,并在特定条件下发生相与相之间的转变。本章将基于热力学基本方程和平衡判据,讨论复相平衡的普遍规律及相变的宏观物理特征。由于相变通常在等温等压条件下进行,吉布斯函数极小值判据是分析相变问题的基础。

3.1 热力学平衡与稳定性判据#

  • 在讨论相变与复相平衡之前,需确立系统处于平衡状态的热力学判据。这涉及热力学势的一阶变分(平衡条件)和二阶变分(稳定性条件)。特别需要指出的是,当选取不同类型的热力学参量(广延量与强度量)作为独立变量时,系统稳定性所对应的数学特性存在显著差异。

    考虑置于等温等压环境中的封闭系统,其状态由温度 TT、压强 pp 以及系统内部参量 ξ\xi(例如某相的物质的量、化学反应进度或序参量)共同描述。根据热力学第二定律,系统在此约束下演化时,吉布斯函数 G(T,p;ξ)G(T, p; \xi) 趋于减小。

    1. 平衡条件(一阶变分):系统达到平衡时,吉布斯函数对内部参量的一阶偏导数为零:
    (Gξ)T,p=0\left(\frac{\partial G}{\partial \xi}\right)_{T,p} = 0
    1. 稳定条件(二阶变分):仅一阶变分为零不足以保证平衡态的稳定性。为了使系统能够抵抗微小涨落,平衡态对应的必须是局部极小值。因此,GG 作为内部参量 ξ\xi 的函数,在平衡点处应满足凸性条件:
    (2Gξ2)T,p>0\left(\frac{\partial^2 G}{\partial \xi^2}\right)_{T,p} > 0
  • 当系统处于平衡态时,内部参量 ξ\xi 不再是独立变量,而是环境参量(如 T,pT, p)的函数。此时,热力学势仅依赖于其各自的自然变量。系统稳定性的微观涨落要求,在不同自然变量空间中表现为不同符号的二阶偏导数。

    1. 广延量空间的凸性

      对于以内能 U(S,V)U(S, V) 为代表的势函数,其自然变量均为广延量。热力学稳定性要求总内能对任何局部的熵或体积涨落具有恢复力,这在数学上等价于要求 UU 是其自变量 SSVV 的凸函数。其二阶变分要求 d2U=dTdSdpdV>0d^2U = dTdS - dpdV > 0。由此导出系统的基本响应函数必须为正:

      热容要求:(2US2)V=TCV>0    CV>0\left(\frac{\partial^2 U}{\partial S^2}\right)_V = \frac{T}{C_V} > 0 \implies C_V > 0

      压缩系数要求:(2UV2)S=1VκS>0    κS>0\left(\frac{\partial^2 U}{\partial V^2}\right)_S = \frac{1}{V\kappa_S} > 0 \implies \kappa_S > 0

    1. 强度量空间的凹性

      当通过勒让德变换引入强度量作为独立变量时,热力学势的凹凸性会发生相应改变。以吉布斯函数 G(T,p)=UTS+pVG(T, p) = U - TS + pV 为例,对其自然变量 TTpp 求二阶偏导数:

      对温度的二阶偏导:

    (2GT2)p=(ST)p=CpT\left(\frac{\partial^2 G}{\partial T^2}\right)_p = -\left(\frac{\partial S}{\partial T}\right)_p = -\frac{C_p}{T}

      由于热力学稳定性要求定压热容 Cp>0C_p > 0,因此必然得到 (2GT2)p<0\left(\frac{\partial^2 G}{\partial T^2}\right)_p < 0

      对压强的二阶偏导:

    (2Gp2)T=(Vp)T=VκT\left(\frac{\partial^2 G}{\partial p^2}\right)_T = \left(\frac{\partial V}{\partial p}\right)_T = -V\kappa_T

      同理,由于等温压缩系数 κT>0\kappa_T > 0,必然得到 (2Gp2)T<0\left(\frac{\partial^2 G}{\partial p^2}\right)_T < 0

    系统热力学稳定性所施加的物理约束表明:平衡态的热力学势对其广延参量呈现凸性(二阶导数为正),而对其强度参量呈现凹性(二阶导数为负)。因此,尽管系统在趋向平衡时吉布斯函数随内部参量 ξ\xi 取极小值,但平衡态的吉布斯函数 Geq(T,p)G_{eq}(T, p) 作为一个多变量曲面,在 (T,p)(T, p) 坐标系中是一个凹函数。这是勒让德变换的必然代数结果。

3.2 复相平衡条件与相图#

  • 设一封闭系统内存在同一物质的两个相(记为相 1 和相 2),系统处于恒温(TT)和恒压(pp)环境中。根据热力学第二定律,在等温等压条件下,系统内部发生的任何自发过程必然导致总吉布斯函数 GG 减小;当系统达到热力学平衡时,总吉布斯函数取极小值,即:

    dG=0dG = 0

    系统的总吉布斯函数为各相吉布斯函数之和(G=G1+G2G = G_1 + G_2)。由于相变过程中伴随着物质在两相间的转移,需要引入广义热力学基本方程:

    dG=SdT+Vdp+μ1dn1+μ2dn2dG = -S dT + V dp + \mu_1 dn_1 + \mu_2 dn_2

    其中,nin_i 为相 ii 的物质的量,μi=(Gini)T,p\mu_i = \left(\frac{\partial G_i}{\partial n_i}\right)_{T,p} 为相 ii 的化学势(对于纯物质,即为摩尔吉布斯函数)。

  • 在等温等压条件下(dT=0,dp=0dT = 0, dp = 0),系统的状态演化需满足:

    dG=μ1dn1+μ2dn20dG = \mu_1 dn_1 + \mu_2 dn_2 \le 0

    由于封闭系统总物质的量守恒,必有 dn1=dn2dn_1 = -dn_2。代入上式得到:

    (μ2μ1)dn20(\mu_2 - \mu_1) dn_2 \le 0

    根据该不等式,可得出以下结论:

    1. 自发过程的方向:若 μ1>μ2\mu_1 > \mu_2,则必定有 dn2>0dn_2 > 0。即物质总是自发地从化学势较高的相转移到化学势较低的相。

    2. 平衡条件:当相变停止、两相达到共存的平衡态时,dG=0dG = 0。由于 dn2dn_2 是相互独立的微小变分,要使 (μ2μ1)dn2=0(\mu_2 - \mu_1) dn_2 = 0 恒成立,必然要求括号内为零。

    综合热学和力学条件,两相平衡必须同时满足以下三个等式:

    1. 热平衡条件T1=T2T_1 = T_2

    2. 力学平衡条件p1=p2p_1 = p_2 (忽略气液等界面的表面张力影响)

    3. 相变平衡条件μ1(T,p)=μ2(T,p)\mu_1(T, p) = \mu_2(T, p)

  • 相变平衡条件 μ1(T,p)=μ2(T,p)\mu_1(T, p) = \mu_2(T, p) 是决定系统状态的核心方程。

    两相共存线:对于确定的物质,μ1\mu_1μ2\mu_2TTpp 的已知函数。方程 μ1(T,p)=μ2(T,p)\mu_1(T, p) = \mu_2(T, p)pTp-T 平面上确定了一条曲线。这条曲线即为两相共存线(如汽化线、熔解线)。在曲线上的任意状态点,两相可以平衡共存;在曲线两侧的区域,化学势较低的单一相成为唯一的稳定相。

    三相点:若系统中有三个相平衡共存,必须同时满足:

    μ1(T,p)=μ2(T,p)=μ3(T,p)\mu_1(T, p) = \mu_2(T, p) = \mu_3(T, p)

    这构成了一个包含两个独立方程的方程组,对应两个未知数 (T,p)(T, p)。一般情况下,该方程组只有唯一解。因此,在 pTp-T 相图上,三条两相共存线必然交汇于一个确定的点,即三相点。在三相点处,系统的温度和压强被完全确定,系统自由度为零。

3.3 一阶相变与克拉珀龙方程(Clapeyron Equation)#

  • 在相平衡点(两相共存线上),相 1 和相 2 的化学势必须严格相等:μ1(T,p)=μ2(T,p)\mu_1(T, p) = \mu_2(T, p)。即化学势在相变过程中是连续的

    但是,化学势的一阶偏导数通常并不连续。根据热力学基本方程 dμ=sdT+vdpd\mu = -s dT + v dp(其中 ssvv 分别为摩尔熵和摩尔体积),我们有:

    (μT)p=s,(μp)T=v\left(\frac{\partial \mu}{\partial T}\right)_p = -s, \quad \left(\frac{\partial \mu}{\partial p}\right)_T = v

    如果在相变发生时,系统状态跨越共存线,虽然 μ\mu 连续,但其一阶偏导数发生跃变,即:

    s1s2,v1v2s_1 \neq s_2, \quad v_1 \neq v_2

    这种相变被称为一阶相变。此时有:

    摩尔熵的突变(Δs=s2s10\Delta s = s_2 - s_1 \neq 0)意味着相变过程中系统需要吸收或放出热量。在等温等压下,这部分热量被称为摩尔相变潜热 L=T(s2s1)L = T(s_2 - s_1)

    摩尔体积的突变(Δv=v2v10\Delta v = v_2 - v_1 \neq 0)意味着相变伴随着宏观的体积膨胀或收缩。

    固-液熔化、液-气汽化、固-气升华以及不同固体结晶形态之间的转变,均属于典型的一阶相变。

  • pTp-T 相图上,两相共存线不仅是一条几何曲线,它的斜率 dpdT\frac{dp}{dT} 蕴含着深层的热力学必然性。假设系统沿两相共存线发生状态改变,两相必须分别保持化学势平衡:

    μ1(T,p)=μ2(T,p)\mu_1(T, p) = \mu_2(T, p)

    在两相共存线上取微分有:

    dμ1=dμ2d\mu_1 = d\mu_2

    代入化学势的全微分表达式 dμ=sdT+vdpd\mu = -s dT + v dp

    s1dT+v1dp=s2dT+v2dp-s_1 dT + v_1 dp = -s_2 dT + v_2 dp

    重新排列各项,分离出 ppTT 的微分:

    (v2v1)dp=(s2s1)dT(v_2 - v_1) dp = (s_2 - s_1) dT dpdT=s2s1v2v1\frac{dp}{dT} = \frac{s_2 - s_1}{v_2 - v_1}

    利用相变潜热定义 L=T(s2s1)L = T(s_2 - s_1) 替换掉极难直接测量的摩尔熵差,我们得到克拉珀龙方程

    dpdT=LT(v2v1)=LTΔv\frac{dp}{dT} = \frac{L}{T(v_2 - v_1)} = \frac{L}{T \Delta v}
  • 克拉珀龙方程是一阶相变中最核心的数学方程,它将相图上的几何特征(共存线斜率)与物质的宏观热力学响应(潜热、体积差)进行了完美的绝对绑定。

    1. 正常物质的熔化(如二氧化碳、绝大多数金属)

      熔化过程吸收热量,因此潜热 L>0L > 0;温度恒定为正 T>0T > 0;固态熔化为液态时体积膨胀,因此 Δv=vlvs>0\Delta v = v_l - v_s > 0

      代入克拉珀龙方程可得 dpdT>0\frac{dp}{dT} > 0

      推论:正常物质的固-液共存线在相图上具有正斜率。外界压强增大,会导致熔点升高。

    2. 水/冰体系的反常膨胀(复冰现象)

      冰熔化成水同样吸热(L>0L > 0),但由于冰的特殊氢键晶体结构,固态的摩尔体积竟然大于液态的摩尔体积(即冰浮在水面上),因此 Δv=vlvs<0\Delta v = v_l - v_s < 0

      代入方程可得绝对结论:dpdT<0\frac{dp}{dT} < 0

      推论:冰的固-液共存线斜率为负。这意味着外界施加的压强越大,冰的熔点反而越低。这正是溜冰鞋能够滑行的热力学基础(冰刀下方的极高局部压强使得冰在冰点以下局部熔化为水,形成润滑层)。

3.4 范德瓦尔斯等温线与麦克斯韦等面积法则#

  • 实际的气液一阶相变在等温条件下是一个等压过程(pvp-v 图上为水平线段)。范德瓦尔斯方程引入了分子体积和分子间吸引力,试图给出气液两相的统一描述,但在临界温度以下其解析形式与物理稳定性条件存在冲突。

    范德瓦尔斯方程为:

    (p+av2)(vb)=RT\left(p + \frac{a}{v^2}\right)(v - b) = RT

    将压强 pp 表示为摩尔体积 vv 的函数,在 pvp-v 图上绘制等温线。

    1. T>TcT > T_c 时,(pv)T<0\left(\frac{\partial p}{\partial v}\right)_T < 0,等温线单调下降,系统仅存在单一气相。

    2. T<TcT < T_c 时,等温线呈现非单调的极值结构。在极小值和极大值之间的区段,满足:

      (pv)T>0\left(\frac{\partial p}{\partial v}\right)_T > 0

    根据热力学稳定性判据,系统的等温压缩系数 κT=1v(vp)T\kappa_T = -\frac{1}{v}\left(\frac{\partial v}{\partial p}\right)_T。在上述区段内 κT<0\kappa_T < 0,系统处于力学不稳定状态,微小的密度涨落将导致系统自发坍缩或膨胀。因此,该波浪线区段在自然界中不能稳定存在。

  • 为了修正非物理的振荡,确定两相共存的实际压强 p0p_0,必须引入复相平衡的化学势判据。

    设饱和液相状态为 A(vl,p0)A(v_l, p_0),饱和气相状态为 C(vg,p0)C(v_g, p_0)。两相平衡共存的必要条件是化学势相等:

    μ(T,p0)A=μ(T,p0)C\mu(T, p_0)_A = \mu(T, p_0)_C

    利用化学势的全微分 dμ=sdT+vdpd\mu = -s dT + v dp。在等温条件下(dT=0dT = 0),沿理论上的范德瓦尔斯等温线(包括不稳定的振荡段)从态 AA 积分至态 CC

    μCμA=ACdμ=ACvdp\mu_C - \mu_A = \int_A^C d\mu = \int_A^C v dp

    代入平衡条件 μC=μA\mu_C = \mu_A,得出:

    ACvdp=0\int_A^C v dp = 0

    几何意义:在 pvp-v 图上,实际相变的水平直线 p=p0p = p_0 必须切割范德瓦尔斯等温线,使得水平线与理论曲线之间围成的上方区域面积(液相侧)与下方区域面积(气相侧)严格相等。此即为确定相变压强的麦克斯韦等面积法则

  • 随着温度升高,等温线上的极大值与极小值相互逼近。当 T=TcT = T_c 时,两极值点重合,波浪结构消失,液相和气相的摩尔体积差 Δv0\Delta v \to 0

    临界点是临界等温线上的拐点,且该点切线水平。其数学条件为压强对体积的一阶和二阶偏导数同时为零:

    (pv)T=Tc=0,(2pv2)T=Tc=0\left(\frac{\partial p}{\partial v}\right)_{T=T_c} = 0, \quad \left(\frac{\partial^2 p}{\partial v^2}\right)_{T=T_c} = 0

    将范德瓦尔斯方程代入上述条件,解得临界常数:

    vc=3b,pc=a27b2,Tc=8a27Rbv_c = 3b, \quad p_c = \frac{a}{27b^2}, \quad T_c = \frac{8a}{27Rb}

    在临界点处,由于 (pv)T=0\left(\frac{\partial p}{\partial v}\right)_T = 0,等温压缩系数发散:κT\kappa_T \to \infty。这意味着流体在临界点的微小压强扰动会导致巨大的密度涨落,一阶相变的宏观边界在此终结。

3.5 连续相变与临界现象#

3.5.1 埃伦费斯特相变分类与二阶相变方程#

  • 为了将不伴随潜热和体积突变的相变纳入严格的热力学框架,埃伦费斯特提出按化学势(或摩尔吉布斯函数 gg)偏导数出现不连续性的最低阶次来对相变进行分类。

    设系统在温度 TT 和压强 pp 下发生相变(由相 1 转变为相 2),两相平衡的核心条件始终是摩尔吉布斯函数连续:

    g1(T,p)=g2(T,p)g_1(T, p) = g_2(T, p)
    1. 一阶相变gg 连续,但 gg 的一阶偏导数不连续。

      由于 (gT)p=s\left(\frac{\partial g}{\partial T}\right)_p = -s(gp)T=v\left(\frac{\partial g}{\partial p}\right)_T = v

      一阶偏导数跃变意味着存在摩尔熵差 Δs0\Delta s \neq 0(即潜热 L=TΔs0L = T\Delta s \neq 0)和摩尔体积差 Δv0\Delta v \neq 0

    2. 二阶相变gg 连续,且 gg 的一阶偏导数也连续,但 gg 的二阶偏导数发生有限的跃变

      一阶导数连续意味着:

      s1=s2    Δs=0(无相变潜热)s_1 = s_2 \implies \Delta s = 0 \quad (\text{无相变潜热}) v1=v2    Δv=0(无体积跃变)v_1 = v_2 \implies \Delta v = 0 \quad (\text{无体积跃变})

      二阶导数跃变对应着系统的三个宏观响应函数发生阶跃:

      定压摩尔热容:cp=T(2gT2)p    Δcp0c_p = -T\left(\frac{\partial^2 g}{\partial T^2}\right)_p \implies \Delta c_p \neq 0

      等温压缩系数:κT=1v(2gp2)T    ΔκT0\kappa_T = -\frac{1}{v}\left(\frac{\partial^2 g}{\partial p^2}\right)_T \implies \Delta \kappa_T \neq 0

      体胀系数:α=1v2gTp    Δα0\alpha = \frac{1}{v}\frac{\partial^2 g}{\partial T \partial p} \implies \Delta \alpha \neq 0

  • 对于一阶相变,我们利用 g1=g2g_1 = g_2 推导出了决定相界斜率的克拉珀龙方程。

    对于二阶相变,由于一阶偏导数在相界线上处处相等,我们可以利用 s1(T,p)=s2(T,p)s_1(T, p) = s_2(T, p)v1(T,p)=v2(T,p)v_1(T, p) = v_2(T, p) 建立新的相界斜率方程。

    设系统沿二阶相变共存线发生微小位移 (dT,dp)(dT, dp),两相的熵和体积必须保持同步变化:

    ds1=ds2    d(Δs)=0ds_1 = ds_2 \implies d(\Delta s) = 0 dv1=dv2    d(Δv)=0dv_1 = dv_2 \implies d(\Delta v) = 0
    1. 由熵的连续性推导第一方程:

      Δs\Delta s 视为 (T,p)(T, p) 的函数,取全微分:

      d(Δs)=(ΔsT)pdT+(Δsp)Tdp=0d(\Delta s) = \left(\frac{\partial \Delta s}{\partial T}\right)_p dT + \left(\frac{\partial \Delta s}{\partial p}\right)_T dp = 0

      利用热力学关系 (sT)p=cpT\left(\frac{\partial s}{\partial T}\right)_p = \frac{c_p}{T},以及麦克斯韦关系 (sp)T=(vT)p=vα\left(\frac{\partial s}{\partial p}\right)_T = -\left(\frac{\partial v}{\partial T}\right)_p = -v\alpha。代入上式:

      ΔcpTdTvΔαdp=0\frac{\Delta c_p}{T} dT - v \Delta \alpha dp = 0

      移项整理,得到埃伦费斯特第一方程

      dpdT=ΔcpTvΔα\frac{dp}{dT} = \frac{\Delta c_p}{T v \Delta \alpha}
    2. 由体积的连续性推导第二方程:

      同理,对 Δv\Delta v 取全微分:

      d(Δv)=(ΔvT)pdT+(Δvp)Tdp=0d(\Delta v) = \left(\frac{\partial \Delta v}{\partial T}\right)_p dT + \left(\frac{\partial \Delta v}{\partial p}\right)_T dp = 0

      利用定义 (vT)p=vα\left(\frac{\partial v}{\partial T}\right)_p = v\alpha(vp)T=vκT\left(\frac{\partial v}{\partial p}\right)_T = -v\kappa_T。代入上式:

      vΔαdTvΔκTdp=0v \Delta \alpha dT - v \Delta \kappa_T dp = 0

      移项整理,得到埃伦费斯特第二方程

      dpdT=ΔαΔκT\frac{dp}{dT} = \frac{\Delta \alpha}{\Delta \kappa_T}
  • 埃伦费斯特方程(dpdT=ΔcpTvΔα=ΔαΔκT\frac{dp}{dT} = \frac{\Delta c_p}{T v \Delta \alpha} = \frac{\Delta \alpha}{\Delta \kappa_T})在数学上结构严密,它预言了在二阶相变点,热容和压缩系数只是发生一个有限高度的阶跃

    然而,当物理学家对超导转变、液氦的正常流体-超流体转变(λ\lambda 相变)以及铁磁体的居里点进行高精度热力学测量时,发现了致命的矛盾:在逼近真实相变温度 TcT_c 时,cpc_pα\alphaκT\kappa_T 并非发生有限的阶跃,而是趋于无穷大(发散)。在数学上,如果 Δcp\Delta c_p \to \inftyΔα\Delta \alpha \to \infty,埃伦费斯特方程将变成 \frac{\infty}{\infty} 的不定型,完全失去了描述相界行为的能力。这一物理危机表明:埃伦费斯特的纯宏观泰勒展开式分类法,在面临临界点附近巨大的微观涨落时彻底失效。物理学必须抛弃“有限阶跃”的图像,转而寻找能够处理“无穷大发散”的新理论。

3.5.2 朗道二级相变理论#

  • 为了从微观对称性机制上解释二阶相变,朗道提出了一套基于自由能极值与稳定性分析的唯象理论。相变过程伴随着物质宏观对称性的改变。高温相通常对称性较高(无序态),低温相通常对称性较低(有序态)。引入序参量 η\eta 定量描述系统偏离高对称态的程度。规定在高温无序相中 η=0\eta = 0,在低温有序相中 η0\eta \neq 0

    在相变点 TcT_c 附近,η\eta 是一个极小量。假设系统的吉布斯函数 G(T,p,η)G(T, p, \eta)TcT_c 附近解析,可按 η\eta 的幂级数展开。由于系统在高对称相具有反演对称性(η\etaη-\eta 物理状态等价),展开式中奇数次幂项系数必为零:

    G(T,p,η)=G0(T,p)+A(T,p)η2+B(T,p)η4+G(T, p, \eta) = G_0(T, p) + A(T, p)\eta^2 + B(T, p)\eta^4 + \cdots

    为了保证系统在任意温度下处于热力学稳定状态(即 η\eta \to \infty 时自由能 GG 不至于发散至负无穷而导致系统崩溃),必须有最高次项系数 B(T,p)>0B(T, p) > 0

  • 系统在热力学平衡时,自由能必须取极小值。求一阶和二阶偏导数:

      平衡条件:Gη=2η(A+2Bη2)=0\frac{\partial G}{\partial \eta} = 2\eta (A + 2B\eta^2) = 0

      稳定条件:2Gη2=2A+12Bη2>0\frac{\partial^2 G}{\partial \eta^2} = 2A + 12B\eta^2 > 0

    考察高对称相(η=0\eta = 0)的稳定性。将 η=0\eta = 0 代入二阶导数得:2Gη2η=0=2A(T,p)\left.\frac{\partial^2 G}{\partial \eta^2}\right|_{\eta=0} = 2A(T, p)

    1. 在高温区 (T>TcT > T_c):实验事实表明,系统平衡于高对称相。即方程必须仅有唯一的实数解 η=0\eta = 0

      η=0\eta = 0 代入稳定条件,要求该极小值点满足 2Gη2η=0=2A>0\left.\frac{\partial^2 G}{\partial \eta^2}\right|_{\eta=0} = 2A > 0。因此严格要求:A(T,p)>0A(T, p) > 0

    2. 在低温区 (T<TcT < T_c):系统发生对称性破缺,平衡于低对称相。即平衡方程除了 η=0\eta=0 之外,必须存在非零的实数解 η0\eta \neq 0

      由平衡条件得出 A+2Bη2=0A + 2B\eta^2 = 0,解得:

      η2=A2B\eta^2 = -\frac{A}{2B}

      由于物理上要求序参量 η\eta 必须为实数(即 η2>0\eta^2 > 0),且前提条件已定 B>0B > 0,因此要使该非零解存在,代数上严格要求:

      (注:此时将非零解代入稳定条件,得到二阶导数为 2A+12B(A2B)=4A2A + 12B(-\frac{A}{2B}) = -4A。既然 A<0A < 0,则 4A>0-4A > 0,证明该非零解确实是一个稳定的极小值点。)

    3. 在临界点 (T=TcT = T_c):作为高温相与低温相的边界,系数 A(T,p)A(T, p) 必须连续地从正值过渡到负值。

      因此,发生二级相变的绝对必要条件为:A(Tc,p)=0A(T_c, p) = 0

  • 既然确定了 A(Tc,p)=0A(T_c, p) = 0,在非常逼近相变点的温度区间内,将系数 A(T,p)A(T, p) 围绕 TcT_c 进行泰勒展开,保留最低阶非零项(线性项):

    A(T,p)(AT)Tc(TTc)A(T, p) \approx \left( \frac{\partial A}{\partial T} \right)_{T_c} (T - T_c)

    a=(AT)Tca = \left( \frac{\partial A}{\partial T} \right)_{T_c},即得到:

    A(T,p)=a(TTc)A(T, p) = a(T - T_c)

    同时,当 T>TcT > T_cA>0A > 0,必然要求系数 a>0a > 0

  • A=a(TTc)A = a(T - T_c) 代回平衡条件 2η(A+2Bη2)=02\eta (A + 2B\eta^2) = 0 进行求解:

    1. 高温相 (TTcT \ge T_c)A0A \ge 0,方程仅有唯一的实数解 η=0\eta = 0

    2. 低温相 (T<TcT < T_c)A<0A < 0,除了不稳定的 η=0\eta=0 外,方程出现两个对称的非零实数解(代表极小值点):

    η=±A2B=±a(TcT)2B\eta = \pm \sqrt{\frac{-A}{2B}} = \pm \sqrt{\frac{a(T_c - T)}{2B}}

    这在数学上完美预言了序参量随温度下降呈现指数为 β=1/2\beta = 1/2 的连续增长。

  • 将低温相的非零解 η2=a(TcT)2B\eta^2 = \frac{a(T_c - T)}{2B} 代回自由能表达式:

    G(T<Tc)=G0a2(TcT)24BG(T < T_c) = G_0 - \frac{a^2(T_c - T)^2}{4B}

    求熵 s=(GT)ps = -\left(\frac{\partial G}{\partial T}\right)_p

    s(T<Tc)=s0a2(TcT)2Bs(T < T_c) = s_0 - \frac{a^2(T_c - T)}{2B}

    T=TcT = T_c 时,Δs=0\Delta s = 0,符合二级相变无潜热特征。

    求定压热容 cp=T(sT)pc_p = T\left(\frac{\partial s}{\partial T}\right)_p

    cp(T<Tc)=cp0+a2T2Bc_p(T < T_c) = c_{p0} + \frac{a^2 T}{2B}

    而高温相 cp(T>Tc)=cp0c_p(T > T_c) = c_{p0}。因此在相变点,热容发生大小为 Δcp=a2Tc2B\Delta c_p = \frac{a^2 T_c}{2B}有限跃变。朗道理论在平均场近似下,从微观对称性出发,极其自洽地推导出了埃伦费斯特二阶相变的全部宏观特征。

3.5.3 涨落的影响与临界指数#

  • 朗道理论本质上是一种“平均场近似”,它假设系统内的序参量是均匀分布的,忽略了热力学涨落(即局部状态偏离宏观平均值的现象)。

    通常情况下涨落的宏观影响可以忽略。但当系统极其逼近相变点 TcT_c 时,决定系统恢复力的系数 A(T)0A(T) \to 0。此时,微小的局部热涨落会被无限放大并达到宏观尺度。

    这种剧烈的涨落彻底破坏了平均场假设。因此,朗道理论预测的热容等参量的“有限跃变”,在真实的实验中变为了趋于无穷大的“发散”。

  • 朗道理论的泰勒展开在临界区失效,物理学转而采用幂律关系来描述发散行为。

    定义无量纲的约化温度 t=TTcTct = \frac{T - T_c}{T_c}。在临界区(t1|t| \ll 1),系统的宏观热力学量服从如下渐近关系:

      序参量 η\etaη(t)β(T<Tc)\eta \propto (-t)^\beta \quad (T < T_c)

      定压热容 cpc_pcptαc_p \propto |t|^{-\alpha}

      等温压缩系数 κT\kappa_TκTtγ\kappa_T \propto |t|^{-\gamma}

    这些指数(α,β,γ\alpha, \beta, \gamma)即为临界指数。真实系统的实验测定值偏离了朗道理论的预言(例如实验中 α>0\alpha > 0 表现为热容发散,而朗道预言 α=0\alpha = 0 即有限跃变)。

    实验表明,临界指数具有深刻的普适性:具有完全不同微观分子间作用力的系统(如气液系统与铁磁自旋系统),只要它们的空间维数和序参量对称性相同,其测得的临界指数就完全一致。

    此外,各临界指数之间并非相互独立,而是满足严格的热力学代数约束(称为标度律),例如:

    α+2β+γ=2\alpha + 2\beta + \gamma = 2

    这表明在连续相变的临界点附近,主导系统宏观奇异行为的不再是微观相互作用的细节,而是系统的全局对称结构。

  

  

  


第4章 多元系复相平衡与化学平衡#

  • 在前三章中,我们处理的主要是处于封闭状态的单组分系统。本章将热力学框架扩展到包含多个组分,且组分的物质的量可以发生变化(由于相变或化学反应)的系统。

4.1 多元系的热力学基本方程与吉布斯-杜亥姆方程#

  • 设系统由 kk 种不同的化学组分构成,各组分的物质的量分别为 n1,n2,,nkn_1, n_2, \dots, n_k

    当系统的组分发生微小变化时,由于引入或移出粒子会改变系统的总能量,热力学基本方程必须引入与粒子数相关的化学功项 i=1kμidni\sum_{i=1}^k \mu_i dn_i

    四个核心热力学函数的全微分方程推广为:

    dU=TdSpdV+i=1kμidnidU = TdS - pdV + \sum_{i=1}^k \mu_i dn_i dH=TdS+Vdp+i=1kμidnidH = TdS + Vdp + \sum_{i=1}^k \mu_i dn_i dF=SdTpdV+i=1kμidnidF = -SdT - pdV + \sum_{i=1}^k \mu_i dn_i dG=SdT+Vdp+i=1kμidnidG = -SdT + Vdp + \sum_{i=1}^k \mu_i dn_i

    由上述方程可得出组分 ii 的化学势 μi\mu_i 的严密数学定义,最常用的是基于吉布斯函数的定义:

    μi=(Gni)T,p,nji\mu_i = \left(\frac{\partial G}{\partial n_i}\right)_{T, p, n_{j \neq i}}

    物理意义:在等温、等压且其他组分物质的量不变的条件下,增加一摩尔组分 ii 所引起的系统总吉布斯函数的增量。

  • 温度 TT、压强 pp 是强度量,而吉布斯函数 GG 和物质的量 nin_i 是广延量。

    在等温等压下,如果将系统中所有组分的物质的量同时放大 λ\lambda 倍,系统的总吉布斯函数也必将放大 λ\lambda 倍。这在数学上表明 GG 是关于各组分 nin_i一阶齐次函数

    G(T,p,λn1,,λnk)=λG(T,p,n1,,nk)G(T, p, \lambda n_1, \dots, \lambda n_k) = \lambda G(T, p, n_1, \dots, n_k)

    根据数学上的欧拉齐次函数定理,一阶齐次函数满足:

    G=i=1kni(Gni)T,p,njiG = \sum_{i=1}^k n_i \left(\frac{\partial G}{\partial n_i}\right)_{T,p,n_{j \neq i}}

    代入化学势的定义,得到多元系吉布斯函数的绝对积分表达式:

    G=i=1kniμiG = \sum_{i=1}^k n_i \mu_i

    这表明,多元系的吉布斯函数严格等于各组分化学势与其物质的量的加权和。

  • 对上述积分形式 G=niμiG = \sum n_i \mu_i 取全微分:

    dG=i=1kμidni+i=1knidμidG = \sum_{i=1}^k \mu_i dn_i + \sum_{i=1}^k n_i d\mu_i

    将其与广义热力学基本方程 dG=SdT+Vdp+μidnidG = -SdT + Vdp + \sum \mu_i dn_i 联立,两式相减,消去 dGdGμidni\sum \mu_i dn_i 项,得到:

    i=1knidμi=SdT+Vdp\sum_{i=1}^k n_i d\mu_i = -SdT + Vdp

    此即为吉布斯-杜亥姆方程

    推论与物理意义

    1. 系统的强度状态参量(T,p,μ1,,μkT, p, \mu_1, \dots, \mu_k)并不是完全独立的。它们的变化受到吉布斯-杜亥姆方程的严格约束。

    2. 在等温等压过程中(dT=0,dp=0dT=0, dp=0),方程退化为:

      i=1knidμi=0\sum_{i=1}^k n_i d\mu_i = 0

      或者除以总物质的量 n=nin = \sum n_i,使用摩尔分数 xi=ni/nx_i = n_i/n 表达:

      i=1kxidμi=0\sum_{i=1}^k x_i d\mu_i = 0

      这表明在恒温恒压下的混合物中,如果某一成分的化学势升高,必然导致其他某些成分的化学势降低,以维持总体的热力学制约关系。

    对于单组分理想气体,在等温条件下 dμ=vdp=RTpdpd\mu = v dp = \frac{RT}{p} dp。积分可得单组分气体的化学势:

    μ(T,p)=μ0(T)+RTlnp\mu(T, p) = \mu^0(T) + RT \ln p

    其中 μ0(T)\mu^0(T)p=1p=1 标准大气压下的标准化学势。

    对于多元理想气体混合物,根据道尔顿分压定律,组分 ii 在混合物中的行为如同它单独占据整个容器时的状态。因此,组分 ii 的化学势仅取决于温度及其分压 pip_i

    μi(T,p)=μi0(T)+RTlnpi\mu_i(T, p) = \mu_i^0(T) + RT \ln p_i

    又因为分压 pi=xipp_i = x_i pxix_i 为该组分的摩尔分数,pp 为系统总压),可将其改写为:

    μi(T,p)=μi0(T)+RTlnp+RTlnxi\mu_i(T, p) = \mu_i^0(T) + RT \ln p + RT \ln x_i

    该方程将化学势清晰地分离为了温度依赖项、总压依赖项以及组分浓度(摩尔分数)依赖项。

4.2 独立组分数与吉布斯相律 (Gibbs Phase Rule)#

  • 在明确了多元系的热力学基本方程后,我们需要解决一个宏观热力学中最普遍的拓扑问题:对于一个包含多个组分和多个相的系统,它究竟拥有多少个可以独立改变的强度状态参量(如温度、压强、浓度)?这个问题的最终答案即为吉布斯相律。

    设系统由 cc 个独立组分构成,共存在 ϕ\phi 个相(记为 α,β,,ϕ\alpha, \beta, \dots, \phi)。当系统达到热力学平衡时,必须同时满足以下三个条件:

    1. 热平衡:所有相的温度相等。
    T(α)=T(β)==T(ϕ)=TT^{(\alpha)} = T^{(\beta)} = \dots = T^{(\phi)} = T
    1. 力学平衡:忽略表面张力和重力场,所有相的压强相等。
    p(α)=p(β)==p(ϕ)=pp^{(\alpha)} = p^{(\beta)} = \dots = p^{(\phi)} = p
    1. 相变与化学平衡每一种组分在各个相中的化学势必须严格相等。对于任意组分 ii (i=1,2,,ci = 1, 2, \dots, c):
    μi(α)(T,p,x1(α),)=μi(β)(T,p,x1(β),)==μi(ϕ)(T,p,x1(ϕ),)\mu_i^{(\alpha)}(T, p, x_1^{(\alpha)}, \dots) = \mu_i^{(\beta)}(T, p, x_1^{(\beta)}, \dots) = \dots = \mu_i^{(\phi)}(T, p, x_1^{(\phi)}, \dots)

    (其中 xi(α)x_i^{(\alpha)} 为组分 iiα\alpha 相中的摩尔分数)。

  • 在实际化学系统中,总物种数(记为 KK)不一定等于独立决定系统组成所需的变量数。如果系统中各物种之间发生化学反应,或存在其他特定的比例约束,组分数会随之减少。

    独立组分数 cc 的计算公式为:`

    c=KRRc = K - R - R'

    KK:系统中的总化学物种数。

    RR:物种之间能够达到化学平衡的独立化学反应数目

    RR':同一相中各物种浓度之间的其他独立限制条件(例如电中性条件,或由初始投料比决定的特定浓度比例)。

  • 自由度 (ff) 定义为在保持系统相数 ϕ\phi 不变的前提下,系统可以独立改变的强度参量(T,pT, p 及各相浓度)的数目。

    根据代数方程理论,自由度等于系统中总强度参量的数目减去这些参量必须满足的独立约束方程数目

    1. 总强度参量数目

      整个系统有一个共同的温度 TT 和压强 pp(共2个)。

      对于每一个相,其组成由 cc 个组分的摩尔分数决定。由于满足归一化条件 i=1cxi=1\sum_{i=1}^c x_i = 1,每个相有 c1c - 1 个独立的摩尔分数。

      因此,ϕ\phi 个相的总强度参量数目为:2+ϕ(c1)2 + \phi(c - 1)

    2. 独立约束方程数目

      复相平衡条件要求每一组分在各相中的化学势相等。对于单一组分 ii,连等式 μi(α)=μi(β)==μi(ϕ)\mu_i^{(\alpha)} = \mu_i^{(\beta)} = \dots = \mu_i^{(\phi)} 包含了 ϕ1\phi - 1 个独立的方程式。

      系统共有 cc 个独立组分,因此化学势等式构成的总约束方程数目为:c(ϕ1)c(\phi - 1)

    3. 计算自由度 ff

      将上述两项相减:

      f=[2+ϕ(c1)]c(ϕ1)f = [2 + \phi(c - 1)] - c(\phi - 1) f=2+ϕcϕcϕ+cf = 2 + \phi c - \phi - c\phi + c f=cϕ+2f = c - \phi + 2

    这就是极其著名的吉布斯相律 (Gibbs Phase Rule)。式中的“2”代表温度和压强这两个外界物理变量。

  • 由于自由度 ff 必须是非负整数(f0f \ge 0),相律对系统可能存在的相数给出了绝对的上限:

    ϕc+2\phi \le c + 2

    单组分系统 (c=1c=1)

    f=3ϕf = 3 - \phi

    由于 f0f \ge 0,最大相数 ϕmax=3\phi_{\max} = 3。当存在三个相(固、液、气)时,f=0f = 0。这证明了单组分系统的三相点在 pTp-T 图上是一个唯一确定的无自由度点,不可能出现四相共存。

    二元系统 (c=2c=2)

    f=4ϕf = 4 - \phi

    最大相数 ϕmax=4\phi_{\max} = 4。如果保持压强恒定(去除了1个自由度),则条件自由度 f=cϕ+1=3ϕf^* = c - \phi + 1 = 3 - \phi。在恒压下的二元合金相图中,最多只能观察到三相共存(如共晶反应点)。

4.3 化学平衡条件与质量作用定律#

  • 在多元相平衡中,我们假设物质在各相间转移但不发生种类的改变。本节将吉布斯极小值判据应用于发生化学反应的系统,推导出决定反应进行程度的普适法则。

    设一封闭系统中发生如下一般形式的化学反应:

    ν1A1+ν2A2++νkAk=0\nu_1 A_1 + \nu_2 A_2 + \dots + \nu_k A_k = 0

    简写为 i=1kνiAi=0\sum_{i=1}^k \nu_i A_i = 0。其中 AiA_i 代表化学物种,νi\nu_i 为化学计量数(规定:生成物 νi>0\nu_i > 0,反应物 νi<0\nu_i < 0)。

    引入反应进度(Extent of Reaction, λ\lambda。当反应发生极小的推进 dλd\lambda 时,各组分物质的量的变化必定与其化学计量数成严格的比例关系:

    dni=νidλdn_i = \nu_i d\lambda

    在等温等压条件下,系统自发演化的方向必须满足吉布斯函数减小(dG0dG \le 0)。代入多元系全微分方程:

    dG=i=1kμidni=(i=1kνiμi)dλ0dG = \sum_{i=1}^k \mu_i dn_i = \left( \sum_{i=1}^k \nu_i \mu_i \right) d\lambda \le 0
    1. 反应自发方向

      νiμi<0\sum \nu_i \mu_i < 0,必然要求 dλ>0d\lambda > 0,反应自发向正方向进行。

      νiμi>0\sum \nu_i \mu_i > 0,必然要求 dλ<0d\lambda < 0,反应自发向逆方向进行。

    2. 化学平衡的绝对条件:当系统达到彻底的热力学平衡时,总吉布斯函数处于极小值,微小变分 dG=0dG = 0。由于 dλd\lambda 可正可负,括号内的项必须严格为零:

      i=1kνiμi=0\sum_{i=1}^k \nu_i \mu_i = 0

    这就是推导一切化学平衡规律的根本方程。它表明在平衡态下,生成物化学势的加权和等于反应物化学势的加权和。

  • 将理想气体混合物组分的化学势公式 μi(T,p)=μi0(T)+RTlnpi\mu_i(T, p) = \mu_i^0(T) + RT \ln p_i 代入平衡条件 νiμi=0\sum \nu_i \mu_i = 0 中:

    i=1kνi[μi0(T)+RTlnpi]=0\sum_{i=1}^k \nu_i \left[ \mu_i^0(T) + RT \ln p_i \right] = 0

    展开并重新组合各项:

    RTi=1kνilnpi=i=1kνiμi0(T)RT \sum_{i=1}^k \nu_i \ln p_i = -\sum_{i=1}^k \nu_i \mu_i^0(T)

    利用对数运算性质 νilnpi=ln(piνi)\sum \nu_i \ln p_i = \ln \left( \prod p_i^{\nu_i} \right),并定义标准摩尔反应吉布斯函数 ΔG0(T)=νiμi0(T)\Delta G^0(T) = \sum \nu_i \mu_i^0(T),上式变为:

    RTln(i=1kpiνi)=ΔG0(T)RT \ln \left( \prod_{i=1}^k p_i^{\nu_i} \right) = -\Delta G^0(T)

    对两边取指数,我们得到了极其重要的质量作用定律(Law of Mass Action)

    i=1kpiνi=Kp(T)\prod_{i=1}^k p_i^{\nu_i} = K_p(T)

    其中,Kp(T)=exp(ΔG0(T)RT)K_p(T) = \exp\left(-\frac{\Delta G^0(T)}{RT}\right) 定义为压强平衡常数

    物理推论

    1. 质量作用定律表明,在达到化学平衡时,各参与气体分压的幂乘积是一个绝对的常数。

    2. 因为 ΔG0(T)\Delta G^0(T) 仅仅是温度 TT 的函数,所以 KpK_p 仅仅是温度的函数,它与系统的总压强、反应容器体积以及反应物的初始投料比例毫无关系

  • 在实际工程计算中,我们更关心物质的浓度(摩尔分数 xix_i)。利用分压与总压的关系 pi=xipp_i = x_i p,代入质量作用定律:

    i=1k(xip)νi=(i=1kxiνi)pνi=Kp(T)\prod_{i=1}^k (x_i p)^{\nu_i} = \left( \prod_{i=1}^k x_i^{\nu_i} \right) p^{\sum \nu_i} = K_p(T)

    定义摩尔分数平衡常数 Kx=i=1kxiνiK_x = \prod_{i=1}^k x_i^{\nu_i},反应前后的气体摩尔数变化为 Δν=νi\Delta \nu = \sum \nu_i。移项可得:

    Kx=Kp(T)pΔνK_x = K_p(T) \cdot p^{-\Delta \nu}

    物理推论

    1. 虽然 KpK_p 仅依赖于温度,但摩尔分数平衡常数 KxK_x 却依赖于系统的总压强 pp

    2. 当反应是一个气体分子数增加的过程(Δν>0\Delta \nu > 0)时,增大总压 pp 会导致 pΔνp^{-\Delta \nu} 减小,从而使 KxK_x 减小,意味着平衡向逆方向(体积缩小方向)移动。这就是我们在下一节探讨勒夏特列原理的数学基础。

4.4 反应热与化学平衡的移动#

  • 在前一节中,我们导出了决定化学平衡态的核心常数 Kp(T)K_p(T)。本节将利用经典热力学基本方程,定量考察外界条件(温度和压强)的变化如何导致平衡态的移动,并严格推导出勒夏特列原理的数学表达。

    在等温等压且不作非体积功的条件下,系统发生化学反应所吸收的热量 QpQ_p 等于系统的焓变 ΔH\Delta H。定义标准摩尔反应焓(反应热):

    ΔH0(T)=i=1kνiHi0(T)\Delta H^0(T) = \sum_{i=1}^k \nu_i H_i^0(T)

    为了寻找平衡常数 KpK_p 与温度的关系,我们需要借助吉布斯函数 GG 与焓 HH 之间的基本联系。由 G=HTSG = H - TS,两边同除以 TT 并对 TT 求恒压偏导:

    [(G/T)T]p=1T(GT)pGT2\left[ \frac{\partial (G/T)}{\partial T} \right]_p = \frac{1}{T}\left(\frac{\partial G}{\partial T}\right)_p - \frac{G}{T^2}

    代入 (GT)p=S\left(\frac{\partial G}{\partial T}\right)_p = -S,得到:

    [(G/T)T]p=STHTST2=HT2\left[ \frac{\partial (G/T)}{\partial T} \right]_p = -\frac{S}{T} - \frac{H - TS}{T^2} = -\frac{H}{T^2}

    这就是著名的吉布斯-亥姆霍兹方程。将其应用于标准态下的化学反应,即有:

    ddT(ΔG0T)=ΔH0T2\frac{d}{dT} \left( \frac{\Delta G^0}{T} \right) = -\frac{\Delta H^0}{T^2}

    (注:因为标准态仅为温度的函数,偏微商变为全微商)。

  • 根据质量作用定律,压强平衡常数 KpK_p 满足:

    ΔG0=RTlnKp    ΔG0T=RlnKp\Delta G^0 = -RT \ln K_p \implies \frac{\Delta G^0}{T} = -R \ln K_p

    将其代入上述吉布斯-亥姆霍兹方程中:

    RdlnKpdT=ΔH0T2-R \frac{d \ln K_p}{dT} = -\frac{\Delta H^0}{T^2}

    整理得到决定化学平衡随温度移动的根本法则——范霍夫等压方程

    dlnKpdT=ΔH0RT2\frac{d \ln K_p}{dT} = \frac{\Delta H^0}{RT^2}

    勒夏特列原理的热力学证明

    (1) 温度对平衡移动的影响

    由范霍夫等压方程可知,KpK_p 随温度的变化率完全由反应热 ΔH0\Delta H^0 的符号决定(由于 RT2>0RT^2 > 0 恒成立):

    • 吸热反应 (ΔH0>0\Delta H^0 > 0)

    dlnKpdT>0\frac{d \ln K_p}{dT} > 0。升高温度 (dT>0dT > 0),必导致 lnKp\ln K_p 增大,平衡常数变大,平衡向**正方向(产物方向)**移动。

    • 放热反应 (ΔH0<0\Delta H^0 < 0)

    dlnKpdT<0\frac{d \ln K_p}{dT} < 0。升高温度 (dT>0dT > 0),必导致 lnKp\ln K_p 减小,平衡常数变小,平衡向**逆方向(反应物方向)**移动。

    (2) 压强对平衡移动的影响

    在 4.3 节中我们导出了摩尔分数平衡常数 KxK_x 与总压 pp 的关系:Kx=Kp(T)pΔνK_x = K_p(T) p^{-\Delta \nu},其中 Δν=νi\Delta \nu = \sum \nu_i 是反应前后气体分子总数的相对增量。

    对上式两边取对数,并在等温条件下对压强 pp 求导:

    lnKx=lnKp(T)Δνlnp\ln K_x = \ln K_p(T) - \Delta \nu \ln p (lnKxp)T=Δνp\left( \frac{\partial \ln K_x}{\partial p} \right)_T = -\frac{\Delta \nu}{p}

    由于 p>0p > 0 恒成立,摩尔分数组成的移动方向完全由 Δν\Delta \nu 的符号决定:

    • 分子数增加的反应 (Δν>0\Delta \nu > 0)

    (lnKxp)T<0\left( \frac{\partial \ln K_x}{\partial p} \right)_T < 0。增大系统总压 (dp>0dp > 0),导致 KxK_x 减小,平衡向**逆方向(分子数减少的方向)**移动。

    • 分子数减少的反应 (Δν<0\Delta \nu < 0)

    (lnKxp)T>0\left( \frac{\partial \ln K_x}{\partial p} \right)_T > 0。增大系统总压 (dp>0dp > 0),导致 KxK_x 增大,平衡向**正方向(也是分子数减少的方向)**移动。

    • 分子数不变的反应 (Δν=0\Delta \nu = 0)

    压强变化不影响 KxK_x,平衡不发生移动。