第零步 · 先看地图
1. 一个装置,两个差,四种效应
前两份都在搭工具。这一份第一次动手用。
先看一个真实现象,它比任何公式都更能说明这一节在讲什么:
拿一块多孔的塞子(烧结玻璃、陶瓷片、甚至一叠滤纸都行),两边灌上稀盐水,然后用力把水从一边压向另一边。
这时如果在两边插上电极、接一个电压表——表针会动。你只是推了水,却量到了电压。
反过来也成立:接上电源加一个电压,水会自己流起来。你只是通了电,液体却被推着走了。
这两件事分别叫流致电势(也叫流动电势,英文 streaming potential,本份两个名字通用)和电渗,它们和另外两个亲戚合称电动效应。这一节要回答的问题是:
- 这些效应为什么会存在?(答案:液体和电荷由同一批载体搬运——贴着孔壁那层带电的液体,这就是交叉效应)
- 它们之间有没有关系?(答案:有,而且是被昂萨格关系锁死的定量关系)
装置:在第 2 份那个装置上再加一个差
场景和第 2 份「情形三」几乎一样:两个容器 I 和 II,中间有一带小孔的隔板(那块多孔塞子)。区别只有一点——这次多了一个压强差。
图 1 电动效应的装置。$T^{\mathrm{I}} = T^{\mathrm{II}}$,但 $p^{\mathrm{I}} \ne p^{\mathrm{II}}$、$\varphi^{\mathrm{I}} \ne \varphi^{\mathrm{II}}$,并设带电粒子为单组元。
条件写清楚就是三行:
$$T^{\mathrm{I}} = T^{\mathrm{II}},\qquad p^{\mathrm{I}} \ne p^{\mathrm{II}},\qquad \varphi^{\mathrm{I}} \ne \varphi^{\mathrm{II}}$$
注意这个装置的关键设计:温度是相同的。
这意味着上一份辛苦算出来的那个力 $\nabla\!\left(\frac{1}{T}\right)$ 在这一节完全不出场。剩下的两个力来自压强差和电势差(第 4 段会看到,它们严格的形式是 $\dfrac{\Delta p}{T}$ 和 $\dfrac{\Delta\varphi}{T}$——那个 $T$ 不能丢)。
这样安排是有道理的:它是能凑出「两个力」的最简单系统。既然要演示昂萨格关系,就至少需要两个力;而温度一旦也不相等,就还要多管一路热流,还没入门就先被复杂度压垮。先在最小的系统里把套路走通,第 4 份再回头处理带温差的情形。
8.1–8.2 · 已讲四条经验定律 → 力与流的框架
8.3 · 已讲熵产生率:力和流已定准
8.4 · 本份电动效应:四个效应压成两个
8.5热电效应:三个系数只需两个
8.6远离平衡:非平衡态相变
这一份的路线只有四步,而且每一步都是把前两份的现成工具搬过来用:
| 步 | 做什么 | 用到的工具 | 段号 |
| 1 | 把压强差翻译成化学势差 | $\left(\dfrac{\partial\mu}{\partial p}\right)_{T} = v$ | 第 2 段 |
| 2 | 算熵产生率 | 第 2 份情形三的现成结果 | 第 3 段 |
| 3 | 读出力与流,写唯象方程 | 判据 $\dfrac{\mathrm{d}_iS}{\mathrm{d}t}=\sum J X$ | 第 4 段 |
| 4 | 解出四个效应,再用昂萨格砍掉一半 | $L_{12}=L_{21}$ | 第 5–7 段 |
①「力」与「流」是什么。系统偏离平衡就会出现梯度(温度差、浓度差、电势差、压强差……),有梯度就有流(热流、物质流、电流……)。非平衡热力学把前者统称力 $X$,后者统称流 $J$。近平衡时两者成线性关系:
$$J_j = \sum_i L_{ji}X_i$$
非对角项 $L_{12}$、$L_{21}$ 就是交叉效应——一种力引起另一种流。这一节四个效应,全部是交叉效应。
② 力和流不能乱选,必须从熵产生率里读出来。判据是:
$$\frac{\mathrm{d}_i S}{\mathrm{d}t} = \sum_{j} J_j X_j$$
把熵产生率整理成「若干项相加、每项两个因子相乘」,每一项里的两个因子就是一对配对的流和力。$\mathrm{d}_i S$ 的下标 $i$ 是 internal,指系统内部自己产生的熵(区别于从外界流进来的熵流 $\mathrm{d}_e S$),它永远 $\ge 0$——这就是热力学第二定律。
③ 昂萨格倒易关系。照上面那个判据定出来的系数,满足
$$L_{ij} = L_{ji}$$
它把独立系数减半,还能从一个效应预言另一个。这一节就是这句话的第一个实战案例。
另外用到第 2 份的一个量:电化学势 $\tilde\mu = \mu + q\varphi$。带电粒子同时是「物质」又是「电荷」,化学势管前者、电势管后者,两笔账必须合成一笔才有意义。$\tilde\mu$ 处处相等,才叫平衡。
落点 · 主干这一份要用的就是这三条
$$\frac{\mathrm{d}_i S}{\mathrm{d}t} = \sum_j J_j X_j \quad\Longrightarrow\quad J_j = \sum_i L_{ji}X_i \quad\text{且}\quad L_{ij}=L_{ji}$$
起点 · 主干里那个现象
$$\text{推水} \;\longrightarrow\; \text{量到电压}$$
热力学本身不回答这个问题——它只管「有没有关系、关系多大」,不管机制。但知道机制会让整节课踏实很多,所以这里说一下。
关键在固液界面。玻璃、陶瓷这类材料泡在水里,表面基团会电离,留下一层固定的负电荷。溶液里的正离子被吸过来,在贴近壁面处形成一层可以流动的正电荷。这一对合称双电层。
于是孔道里的液体不再是电中性均匀的:靠近壁面的那一薄层,净带正电。接下来两件事就都顺理成章了:
| 你做什么 | 发生什么 | 叫什么 |
| 用压强推水 | 水把那层正电荷一起拖走 → 一边积正电、一边积负电 → 出现电压 | 流动电势 |
| 加一个电压 | 电场推动那层正电荷 → 正电荷靠黏滞拽着整团水走 | 电渗 |
两件事共用同一个载体:贴壁那层带电的液体。这正是第 2 份里说过的判断法则——两个过程有没有交叉效应,取决于它们是不是共用同一批载体。共用,所以耦合;耦合,所以昂萨格关系把它们锁在一起。
但请注意热力学的立场:上面这套双电层图像,热力学一个字也用不上。它照样能推出四个效应之间的定量关系,而且不管孔壁是玻璃、陶瓷还是生物膜,结论完全一样。
这就是这套方法真正厉害的地方:它把「机制」和「关系」彻底分开了。机制交给界面化学去研究,关系由热力学一次性锁定——换个材料,$L_{ij}$ 的数值会变,但 $L_{12}=L_{21}$ 岿然不动。
第一步 · 把压强翻译成化学势
2. 熵产生率里没有「压强」这一项
动手第一步就撞上一个麻烦,值得慢慢说清楚。
第 2 份算出来的熵产生率,力是化学势差(或电化学势差)。可这一节的装置里,两边的差别是压强和电势。压强差不在那个式子里,怎么办?
办法是:不去改公式,而是把压强差翻译成化学势差。
因为化学势本来就是压强的函数——压强一变,化学势跟着变。所以「两边压强不同」这句话,本身就等价于「两边化学势不同」。
课本把这件事写成一行:化学势是三个变量的函数
$$\mu^{\alpha} = \mu^{\alpha}\big(T,\ p^{\alpha},\ n^{\alpha}\big),\qquad \alpha = \mathrm{I},\ \mathrm{II}$$
其中 $n^{\alpha} = \dfrac{1}{v^{\alpha}}$,而 $v^{\alpha} = \dfrac{V^{\alpha}}{N^{\alpha}}$,$N^{\alpha}$ 为物质的量,$v^{\alpha}$ 为比体积,$n^{\alpha}$ 为粒子数密度。(「粒子数密度」这个叫法课本写得不严谨:既然 $N^{\alpha}$ 是物质的量,$1/v^{\alpha}$ 的单位就是 $\mathrm{mol/m^3}$,严格说是摩尔浓度,与「每立方米几个粒子」差一个阿伏伽德罗常数。不影响本节任何结论——$n$ 在后面根本没被用到。)
于是结论是:由于两个容器的压强不相等,而化学势又是压强的函数,所以实际上就是两个容器中粒子的化学势不等,即 $\mu^{\mathrm{I}} \ne \mu^{\mathrm{II}}$。再加上本来就有的电势差——这个情况就和第 2 份的「情形三」完全相同了。
翻译的汇率是多少?
知道「压强差会造成化学势差」还不够,我们需要定量的换算。这里用到热力学的一条基本关系:
$$\Big(\frac{\partial\mu}{\partial p}\Big)_{T} = v$$
翻译成大白话:压强每升高一点,化学势就升高「一个比体积」那么多。
由于偏离平衡态不远(压强差很小),可以直接用这个斜率做线性近似,把微分换成差分:
$$\Delta\mu = v\,\Delta p,\qquad \Delta p = p^{\mathrm{I}} - p^{\mathrm{II}}$$
「偏离平衡态不远」这五个字不是客套话,它是这一步成立的前提。
严格的关系是微分式 $\mathrm{d}\mu = v\,\mathrm{d}p$,其中 $v$ 本身还随压强变化。只有当 $\Delta p$ 足够小、$v$ 在这段区间里几乎是常数时,才能把它提出来写成 $\Delta\mu = v\Delta p$。
这和第 1 份里「流正比于力」是同一个近似的两次使用——整套线性非平衡热力学,都活在这个「近」字上。
化学势 $\mu$:在温度、压强不变的条件下,往系统里再加一点点物质,吉布斯函数 $G$ 增加多少。
$$\mu = \Big(\frac{\partial G}{\partial N}\Big)_{T,\,p}$$
单位 $\mathrm{J/mol}$。最好用的直觉是把它类比成温度和压强:
- 温度不等 $\Rightarrow$ 热会流动,从高温流向低温;
- 压强不等 $\Rightarrow$ 体积会变化,从高压推向低压;
- 化学势不等 $\Rightarrow$ 物质会流动,从高 $\mu$ 流向低 $\mu$。
三者是同一张桌子上的三个位置。第 2 份已经把最后一条严格证明出来了(限于单组元)。
比体积 $v$:每单位物质占多大体积。
$$v = \frac{V}{N}$$
本节里 $N$ 是物质的量(单位 mol),所以 $v$ 的单位是 $\mathrm{m^3/mol}$,也就是通常说的摩尔体积。它是密度的倒数一类的量:东西越松散,$v$ 越大。
| 物质(常温常压) | 大致的 $v$ |
| 水 | $1.8\times10^{-5}\ \mathrm{m^3/mol}$(18 mL/mol) |
| 理想气体 | $2.4\times10^{-2}\ \mathrm{m^3/mol}$(24 L/mol) |
气体的 $v$ 比液体大三个数量级——这就是为什么同样的压强差,对气体造成的化学势差要大得多。
落点 · 主干正在用它们的地方
$$\Delta\mu = v\,\Delta p$$
起点 · 就是上面主干里的这个关系
$$\Big(\frac{\partial\mu}{\partial p}\Big)_{T,\,n} = v$$
第 1 步 · 单组元时,化学势就是「每摩尔的吉布斯函数」
对只有一种物质的系统(本节正是这样,课本明说了「设带电粒子为单组元」),吉布斯函数与化学势的关系特别简单:
$$G = \mu N \qquad\Longleftrightarrow\qquad \mu = \frac{G}{N}$$
为什么?分两步。第一步,$G$ 是广延量(物质翻倍,$G$ 翻倍),所以它一定能写成 $G = N\,g(T,p)$,其中 $g$ 是某个强度量——把 $N$ 摩尔看成 $N$ 份一模一样的一摩尔,总量就是每份的 $N$ 倍。
第二步,这个 $g$ 到底是不是 $\mu$?用化学势的定义验一下:
$$\mu = \Big(\frac{\partial G}{\partial N}\Big)_{T,p} = \frac{\partial (N g)}{\partial N} = g$$
是的。这一步不能省——光说「$G$ 广延、$\mu$ 强度」推不出 $G=\mu N$,广延性只告诉你 $G$ 能写成 $N$ 乘上某个强度量,至于那个强度量是不是 $\mu$,只有用定义验过才知道。
第 2 步 · 写出吉布斯函数的微分
这是热力学里的标准结果(由 $G = U - TS + pV$ 与吉布斯方程直接得到):
$$\mathrm{d}G = -S\,\mathrm{d}T + V\,\mathrm{d}p \qquad (\mathrm{d}N = 0)$$
它说的是:$G$ 的两个「自然变量」正是 $T$ 和 $p$——恰好就是我们这个装置控制的两个量。注意括号里的条件:物质的量不变时才只有这两项,否则还要加上 $+\mu\,\mathrm{d}N$。
第 3 步 · 两边同除以 $N$,变成「每摩尔」的版本
令 $s = S/N$(比熵)、$v = V/N$(比体积)。严格的做法是:由 $G=\mu N$ 得 $\mathrm{d}G = N\,\mathrm{d}\mu + \mu\,\mathrm{d}N$,与开系的 $\mathrm{d}G = -S\mathrm{d}T + V\mathrm{d}p + \mu\,\mathrm{d}N$ 相减,$\mu\,\mathrm{d}N$ 消掉,剩下 $N\,\mathrm{d}\mu = -S\,\mathrm{d}T + V\,\mathrm{d}p$(这就是吉布斯–杜安关系),再两边除以 $N$:
$$\mathrm{d}\mu = -s\,\mathrm{d}T + v\,\mathrm{d}p$$
最后固定温度(本节两容器温度相同,$\mathrm{d}T = 0$),只剩一项:
终点 · 回到主干的这个式子
$$\Big(\frac{\partial\mu}{\partial p}\Big)_{T} = v$$
顺手记住这个式子的另一半,第 4 份要用:$\left(\dfrac{\partial\mu}{\partial T}\right)_{p} = -s$。
也就是说,压强通过 $v$ 影响化学势,温度通过 $-s$ 影响化学势。这一节只用到前一半(因为温度相同);到了热电效应,用的就是后一半了。
起点 · 主干里的这组定义
$$v^{\alpha} = \frac{V^{\alpha}}{N^{\alpha}},\qquad n^{\alpha} = \frac{1}{v^{\alpha}}$$
| 符号 | 叫什么 | 单位 | 是不是强度量 |
| $V$ | 体积(整个容器的) | $\mathrm{m^3}$ | 否,广延量 |
| $v$ | 比体积(每摩尔的) | $\mathrm{m^3/mol}$ | 是 |
| $n = 1/v$ | 浓度(每立方米有多少摩尔) | $\mathrm{mol/m^3}$ | 是 |
错误的理解:把 $\left(\frac{\partial\mu}{\partial p}\right)_T = v$ 里的 $v$ 当成容器的体积 $V$。
为什么错:量纲立刻就对不上。$\mu$ 的单位是 $\mathrm{J/mol}$,$p$ 的单位是 $\mathrm{Pa} = \mathrm{J/m^3}$,所以
$$\Big[\frac{\partial\mu}{\partial p}\Big] = \frac{\mathrm{J/mol}}{\mathrm{J/m^3}} = \frac{\mathrm{m^3}}{\mathrm{mol}}$$
正是比体积的单位,不是体积的单位。更根本的理由:等式左边 $\mu$ 是强度量,右边就必须也是强度量;而 $V$ 是广延量,放进去这个式子会随容器大小而变,那就荒唐了。
课本这里还有一处更要紧的记号松懈:$n$ 根本不是独立变量。对单组元流体,强度状态由 $(T,p)$ 就完全定死了,浓度 $n$ 由物态方程 $n = n(T,p)$ 跟着决定。所以写 $\mu(T,p,n)$ 是过定的,而 $\left(\frac{\partial\mu}{\partial p}\right)_{T,n}$ 这种「固定 $T$ 和 $n$ 再对 $p$ 求导」的操作严格说没有定义——正确的写法就是主干里的 $\left(\frac{\partial\mu}{\partial p}\right)_{T} = v$。
课本把 $n$ 列进去,暗示的其实是「溶质溶在溶剂里」的多组元图像(那里浓度确实是独立的,$v$ 也要改读成偏摩尔体积)。这和下一段要说破的那处理想化是同一个病根。
另一处叫法上的松懈:课本把 $n = 1/v$ 称作「粒子数密度」,但既然 $N$ 是物质的量(mol),$1/v$ 的单位就是 $\mathrm{mol/m^3}$,严格说是摩尔浓度而不是「每立方米几个粒子」。
两者只差一个阿伏伽德罗常数,不影响本节任何结论($n$ 在后面根本没被用到,它只是化学势那三个自变量之一)。知道有这回事就行。
第二步 · 算熵产生率
3. 把现成的结果搬过来,代一次
翻译做完了,现在可以直接套用第 2 份的现成结果——因为「化学势差 + 电势差」这个局面,和那里的「情形三」一模一样。
第 2 份情形三算出的熵产生率是(单组元,去掉求和号):
$$\frac{\mathrm{d}_i S}{\mathrm{d}t} = -\Big(\frac{\mu^{\mathrm{I}} + q\varphi^{\mathrm{I}}}{T} - \frac{\mu^{\mathrm{II}} + q\varphi^{\mathrm{II}}}{T}\Big)\cdot\frac{\mathrm{d}N^{\mathrm{I}}}{\mathrm{d}t}$$
括号里正是电化学势差:$\tilde\mu = \mu + q\varphi$,这个式子就是 $-\Delta\!\left(\dfrac{\tilde\mu}{T}\right)\cdot\dfrac{\mathrm{d}N^{\mathrm{I}}}{\mathrm{d}t}$。把它按「化学势」和「电势」拆成两项:
$$\frac{\mathrm{d}_i S}{\mathrm{d}t} = -\frac{\Delta\mu}{T}\cdot\frac{\mathrm{d}N^{\mathrm{I}}}{\mathrm{d}t} - \frac{q\,\Delta\varphi}{T}\cdot\frac{\mathrm{d}N^{\mathrm{I}}}{\mathrm{d}t}$$
最后一步:把上一段的翻译 $\Delta\mu = v\,\Delta p$ 代进第一项,压强差就正式进场了:
$$\frac{\mathrm{d}_i S}{\mathrm{d}t} = -\frac{\Delta p}{T}\cdot v\,\frac{\mathrm{d}N^{\mathrm{I}}}{\mathrm{d}t} - \frac{\Delta\varphi}{T}\cdot q\,\frac{\mathrm{d}N^{\mathrm{I}}}{\mathrm{d}t}$$
(8.4.1)
盯住 (8.4.1) 的形状,它已经把答案摆在脸上了。右边是两项相加,每一项都是两个因子相乘——这正是判据 $\dfrac{\mathrm{d}_i S}{\mathrm{d}t} = \sum_j J_j X_j$ 要求的样子。
所以这个系统有两个力、两个流。下一段只需要把它们从式子里读出来、起个名字。整节最难的部分到这里就结束了。
起点 · 主干直接搬来用的那个式子
$$\frac{\mathrm{d}_i S}{\mathrm{d}t} = -\Delta\Big(\frac{\tilde\mu}{T}\Big)\cdot\frac{\mathrm{d}N^{\mathrm{I}}}{\mathrm{d}t}$$
这一框是给「不想欠账」的读者准备的。不看它完全不影响往下读——上面那个式子是第 2 份的成品,可以当已知条件用。
第 1 步 · 对每一室写吉布斯方程,再加起来
每一室自己处于平衡态(局部平衡假设),所以对每一室都成立 $T\,\mathrm{d}S^{\alpha} = \mathrm{d}U^{\alpha} + p^{\alpha}\mathrm{d}V^{\alpha} - \mu^{\alpha}\mathrm{d}N^{\alpha}$。两室相加:
$$T\,\mathrm{d}S = \sum_{\alpha}\Big(\mathrm{d}U^{\alpha} + p^{\alpha}\mathrm{d}V^{\alpha}\Big) - \sum_{\alpha}\mu^{\alpha}\mathrm{d}N^{\alpha}$$
请注意 $p\,\mathrm{d}V$ 项还在。它不能像第 2 份「情形二」那样以「体积不变」为由划掉——两边压强不等、又要有净流动,对不可压缩液体来说体积就不可能都不变,实际装置正是靠活塞或液柱维持这个压差的。
第 2 步 · 关键:那两项 $p\,\mathrm{d}V$ 并没有消失,而是并进了热量
整个体系的第一定律是「内能的增加 = 吸的热 $-$ 对外做的功」,而这里对外做的功正是推动两个活塞的功 $W = p^{\mathrm{I}}\mathrm{d}V^{\mathrm{I}} + p^{\mathrm{II}}\mathrm{d}V^{\mathrm{II}}$:
$$\mathrm{d}U^{\mathrm{I}} + \mathrm{d}U^{\mathrm{II}} = \delta Q - W$$
移项就看出来了:第 1 步那个括号整体就是热量——
$$\sum_{\alpha}\Big(\mathrm{d}U^{\alpha} + p^{\alpha}\mathrm{d}V^{\alpha}\Big) = \delta Q$$
所以 $p\,\mathrm{d}V$ 不是被「抵消」掉的,而是被吸收进了「热量」这个名字里。再用 $\mathrm{d}N^{\mathrm{II}} = -\mathrm{d}N^{\mathrm{I}}$ 把 II 那一项换掉:
$$\mathrm{d}S = \underbrace{\frac{\delta Q}{T}}_{\text{熵流}} - \underbrace{\frac{1}{T}\big(\mu^{\mathrm{I}} - \mu^{\mathrm{II}}\big)\mathrm{d}N^{\mathrm{I}}}_{\text{熵产生}}$$
第一项是与外界换热带来的(跨过体系边界)= 熵流;第二项是体系内部两容器之间倒腾粒子造成的 = 熵产生。形式和第 2 份完全一样,只是那里的「$\mathrm{d}Q$」在这里要理解成含了活塞功在内的真正热量。
第 3 步 · 粒子带电,再加一笔电势能的账
有了电场之后,内能上要加一项 $q\varphi^{\alpha}\mathrm{d}N^{\alpha}$;交换粒子时这部分能量也跟着迁移。于是熵产生那一项里,$\mu$ 处处换成 $\tilde\mu = \mu + q\varphi$:
终点 · 回到主干搬来的那个式子
$$\frac{\mathrm{d}_i S}{\mathrm{d}t} = -\Big(\frac{\tilde\mu^{\mathrm{I}}}{T} - \frac{\tilde\mu^{\mathrm{II}}}{T}\Big)\frac{\mathrm{d}N^{\mathrm{I}}}{\mathrm{d}t}$$
为什么可以「原样搬来」而不用重算:因为第 2 份情形三对装置的要求只有两条——两容器温度相同、粒子带电。本节的装置这两条都满足,至于压强等不等,那个推导从头到尾没有用到。
压强是通过化学势间接进来的(上一段那步翻译),而不是作为一个新的独立项——$p\,\mathrm{d}V$ 那些项已经被并进「热量」里去了,留在熵产生那一项里的只有 $\mu$。这正是上一段那步翻译的全部价值。
起点 · 主干里定下的约定
$$\Delta p = p^{\mathrm{I}} - p^{\mathrm{II}}$$
本节全程是「I 减 II」,即前者减后者,$\Delta\varphi = \varphi^{\mathrm{I}} - \varphi^{\mathrm{II}}$ 同理。这是个纯粹的约定,本身没有对错,但一旦定下就必须从头贯彻到尾。
为什么它特别要紧:因为第 5–7 段要算的四个效应,全部是比值,比如 $\left(\dfrac{\Delta\varphi}{\Delta p}\right)_{J_e=0}$。中途换一次方向,分子分母就各带一个负号——
- 两个都换,比值不变(还好);
- 只换一个,比值整体变号(完蛋)。
而 (8.4.9) 那两条结论恰恰是关于符号的(一正一负)。方向搞混,结论就直接读反。
一个能自查的办法:算完之后,用物理直觉验一下符号。比如「把水从 I 压向 II」应当对应 $\Delta p = p^{\mathrm{I}} - p^{\mathrm{II}} > 0$,而物质确实往 II 走。如果你算出的流方向和这个对不上,多半就是某一步把 $\Delta$ 写反了。
这个习惯值得养成:任何带方向的结果,算完都用一个极端情形试一下。比第二遍检查代数有效得多。
第三步 · 读出力与流
4. 两个力,两个流,四个系数
(8.4.1) 已经是「两项相加、每项两因子相乘」的形状。按判据,每一项里的两个因子就是一对配对的流和力。给它们起名字:
$$J_e = -q\,\frac{\mathrm{d}N^{\mathrm{I}}}{\mathrm{d}t}\quad(\text{电流}),\qquad J_m = -v\,\frac{\mathrm{d}N^{\mathrm{I}}}{\mathrm{d}t}\quad(\text{物质流})$$
这里必须说破课本的一处理想化,否则后面会自相矛盾。
照上面的定义,两个流都正比于
同一个 $\dfrac{\mathrm{d}N^{\mathrm{I}}}{\mathrm{d}t}$,于是它们的比值恒为常数:
$$\frac{J_e}{J_m} = \frac{q}{v} = \text{常数}$$
严格照这个模型走,两个流
并不独立——$J_e = 0$ 就意味着 $\dfrac{\mathrm{d}N^{\mathrm{I}}}{\mathrm{d}t}=0$,于是 $J_m$ 也为零。
那么第 1 段讲的「开路推水、却量到电压」就根本不可能发生了。
再往下追一层,会看得更清楚。把 $J_e \equiv \frac{q}{v}J_m$ 这个恒等式代进 (8.4.2)(8.4.3),会解出 $L_{11} = \frac{q}{v}L_{21}$、$L_{12} = \frac{q}{v}L_{22}$,于是系数矩阵
退化成秩 1:
$$\det L = L_{11}L_{22} - L_{12}L_{21} = 0$$
再用上 $L_{12}=L_{21}$,代进 (8.4.5) 就得到 $\left(\frac{\Delta\varphi}{\Delta p}\right)_{J_e=0} = -\frac{v}{q}$——一个
只含 $v$ 和 $q$、与材料完全无关的普适数。而实测的流动电势系数明明随材料和浓度变化好几个数量级。更要命的是,此时 $J_e=0$ 其实等价于 $v\Delta p + q\Delta\varphi = 0$,也就是 $\Delta\tilde\mu = 0$——
那是平衡态、力本身为零,根本没有什么「效应」可言。
病根在「设带电粒子为单组元」这个简化上。真实的电动效应至少需要
两个组元:体积流主要由
中性的溶剂承担,电流由双电层里那
一点点净离子电荷承担——两者由不同的东西搬运,才可能各走各的。
课本这样设,是为了让 $\Delta\mu = v\Delta p$ 和 $G=\mu N$ 这两个只有单组元才有的简单关系能用。
好消息是:四个效应之间的关系式 (8.4.5)–(8.4.9),在多组元的严格处理中形式完全不变(前提是
两侧是同一种溶液、只差压强和电势;若两侧浓度也不同,就还要多出一对共轭量,变成 $3\times3$ 的问题)。把体积流 $J_V=\sum_k v_k J_k$ 与电流 $I=F\sum_k z_k J_k$ 当作两个流,熵产生率照样整理成 $\frac{\Delta p}{T}J_V + \frac{\Delta\varphi}{T}I$ 两项,而且新的系数矩阵仍然对称——
这一层才是它们能和真实实验对上的原因。顺带一提,(8.4.9) 左边那条关系在昂萨格之前就被实验发现了,叫
Saxén 关系(1892 年)。
所以下面照着算就行,但心里要清楚:两个流的独立性来自「溶剂和离子是两拨人」,而不是来自这个单组元式子。
代回 (8.4.1),它就变成了标准形状:
$$\frac{\mathrm{d}_i S}{\mathrm{d}t} = J_m\cdot\frac{\Delta p}{T} + J_e\cdot\frac{\Delta\varphi}{T}$$
| 流 $J$ | 配对的力 $X$ | 物理意思 |
| 第一对 | $J_e = -q\,\dfrac{\mathrm{d}N^{\mathrm{I}}}{\mathrm{d}t}$ | $\dfrac{\Delta\varphi}{T}$ | 每秒有多少电荷穿过隔板 |
| 第二对 | $J_m = -v\,\dfrac{\mathrm{d}N^{\mathrm{I}}}{\mathrm{d}t}$ | $\dfrac{\Delta p}{T}$ | 每秒有多少体积的液体穿过隔板 |
有了力和流,线性唯象方程就可以写出来了。注意是每个流对两个力都有响应:
$$J_e = L_{11}\,\frac{\Delta\varphi}{T} + L_{12}\,\frac{\Delta p}{T}$$
(8.4.2)
$$J_m = L_{21}\,\frac{\Delta\varphi}{T} + L_{22}\,\frac{\Delta p}{T}$$
(8.4.3)
课本在这里加了一句很重的话:只有像上面那样写出的流与力的两个方程,系数之间才有昂萨格关系:
$$L_{12} = L_{21}$$
(8.4.4)
「只有像上面那样写出的」——这半句话是整节的地基。
但要说准:
并不是随便动一下就会失效。比如嫌 $T$ 碍事,把
两个力同时乘上 $T$(改用 $\Delta\varphi$ 和 $\Delta p$ 当力、流不动),系数矩阵只是
整体除以 $T$,$L_{12}=L_{21}$
照样成立。
真正会出事的是「拆散配对」——只对其中一对动手。比如只把电那一对的力改成 $\Delta\varphi$,物质那一对仍用 $\dfrac{\Delta p}{T}$:
$$L'_{12} = L_{12},\qquad L'_{21} = \frac{L_{21}}{T}$$
改的是第
一个力,动的就是系数矩阵的第
一列($L_{11}$ 与 $L_{21}$);$L_{12}$ 是「力 2 对流 1 的响应」,跟力 1 怎么写无关,
纹丝不动。一个除了 $T$、一个没除,原本相等的两个立刻就不等了。
这正是第 2 份那条:力和流必须
成对地从熵产生率里读出来。照这个办法读出来的一套,昂萨格关系才
有依据;拆散配对之后,就没有任何东西保证它了。而这一节最后的成果(四个效应压成两个)
全靠 (8.4.4)。
四个系数里,$L_{11}$、$L_{22}$ 是「老实人」——电势差引起电流、压强差引起物质流,都是意料之中。真正有意思的是 $L_{12}$ 和 $L_{21}$:
| 系数 | 它描述的因果 | 身份 |
| $L_{11}$ | 电势差 $\to$ 电流 | 老熟人:欧姆定律 |
| $L_{22}$ | 压强差 $\to$ 物质流 | 老熟人:达西定律(液体渗流) |
| $L_{12}$ | 压强差 $\to$ 电流 | 交叉项:推水居然出电 |
| $L_{21}$ | 电势差 $\to$ 物质流 | 交叉项:通电居然推水 |
这一节剩下的全部内容,就是把这四个系数翻译成四个可以实际测量的量,然后看看 (8.4.4) 会带来什么。
起点 · 主干里的定义
$$J_e = -q\,\frac{\mathrm{d}N^{\mathrm{I}}}{\mathrm{d}t},\qquad J_m = -v\,\frac{\mathrm{d}N^{\mathrm{I}}}{\mathrm{d}t}$$
负号不是算出来的,是定出来的。它只是在选「流的正方向」。
$\dfrac{\mathrm{d}N^{\mathrm{I}}}{\mathrm{d}t} > 0$ 的意思是「容器 I 里的物质在增加」,也就是东西流进 I、方向是 II → I。而我们希望「流为正」表示从 I 流向 II(和 $\Delta = \mathrm{I} - \mathrm{II}$ 的方向一致),所以加个负号把方向掰过来。
怎么确认这个选择是对的?看它有没有让式子自洽。代回熵产生率:
$$\begin{aligned}\frac{\mathrm{d}_i S}{\mathrm{d}t} &= -\frac{\Delta p}{T}\Big(v\frac{\mathrm{d}N^{\mathrm{I}}}{\mathrm{d}t}\Big) - \frac{\Delta\varphi}{T}\Big(q\frac{\mathrm{d}N^{\mathrm{I}}}{\mathrm{d}t}\Big)\\[3pt] &= \frac{\Delta p}{T}\,J_m + \frac{\Delta\varphi}{T}\,J_e\end{aligned}$$
两个负号都被吸收掉了,剩下干干净净的「力 × 流」求和,正是判据要的形状。这就是选这个方向的理由——它让式子最漂亮。
再做一次物理自查:设 $q > 0$(正离子),把液体从 I 压向 II,即 $\Delta p > 0$。此时物质离开 I,$\dfrac{\mathrm{d}N^{\mathrm{I}}}{\mathrm{d}t} < 0$,于是 $J_m = -v\dfrac{\mathrm{d}N^{\mathrm{I}}}{\mathrm{d}t} > 0$。
「压强差为正,物质流也为正」——方向对上了。而且此时 $J_m\dfrac{\Delta p}{T} > 0$:在没有电势差($\Delta\varphi = 0$)的情形下,这就是全部的熵产生,确实为正,第二定律满足。两个检查都过,说明约定自洽。(两个力都在时要看两项之和,单看一项不足以下结论。)
顺带交代一下 $q$ 的身份:既然 $N$ 是物质的量(mol),为了让 $q\Delta\varphi$ 的单位是 $\mathrm{J/mol}$、$J_e = -q\dot N$ 的单位是安培,这里的 $q$ 必须是每摩尔的电荷量($\mathrm{C/mol}$,即 $zF$,$F$ 是法拉第常数),而不是单个粒子的电荷。上面说「$q>0$ 即正离子」时,指的也是这个量的符号。
起点 · 主干里的定义
$$J_m = -v\,\frac{\mathrm{d}N^{\mathrm{I}}}{\mathrm{d}t}$$
数一下单位:$v$ 是 $\mathrm{m^3/mol}$,$\dfrac{\mathrm{d}N}{\mathrm{d}t}$ 是 $\mathrm{mol/s}$,乘起来是
$$[J_m] = \frac{\mathrm{m^3}}{\mathrm{mol}}\cdot\frac{\mathrm{mol}}{\mathrm{s}} = \frac{\mathrm{m^3}}{\mathrm{s}}$$
立方米每秒——这是体积流量。所以 $J_m$ 名叫「物质流」,量的其实是每秒有多少体积的液体穿过隔板。
这不是随手为之,而是刻意的选择,有三个好处:
- 它是实验室里直接能读的量。接一根量管看液面走多快,或者干脆称重,就得到 $\mathrm{m^3/s}$。而「每秒多少摩尔」还得先知道摩尔质量、再换算。
- 它让力变成压强。力与流相乘必须给出熵产生率;配 $\Delta p$(单位 $\mathrm{J/m^3}$)的流,量纲上就必须是 $\mathrm{m^3/s}$。是「力要用压强」这个愿望,反过来决定了流必须用体积。
- 它让昂萨格关系变得可以直接验证。下一段你会看到,正因为 $J_m$ 是体积流,$\left(\dfrac{J_m}{J_e}\right)$ 和 $\left(\dfrac{\Delta\varphi}{\Delta p}\right)$ 才会有完全相同的单位——两个八竿子打不着的实验,数值可以直接对比。
这里体现了一条通用的经验:$X$ 和 $J$ 怎么分家,是有讲究的。只要乘积仍等于熵产生率,你可以把一个因子挪到另一个上去(比如把 $v$ 从流挪进力)。挑哪一种,标准是「哪一对更接近仪器读数」。
这一节挑的这一对,让两个因子恰好是压强差和体积流量——压力表和流量计各读一个,谁都不用换算。
第四步 · 四个效应(上)· 掐断一个流
5. 不让它流,看它憋出多大的差
(8.4.2)(8.4.3) 是两个方程、四个量(两个力、两个流)。要从里面读出可测的结论,办法是每次固定一个量为零,看剩下的三个之间是什么关系。
四种固定方式,正好给出四个效应。课本按顺序编号为 (1)–(4),但按你固定的是「流」还是「力」来分两组,会好记得多:
| 这一段(掐断一个流) | 下一段(归零一个力) |
| 令 $J_e = 0$:流致电势 | 令 $\Delta p = 0$:电渗效应 |
| 令 $J_m = 0$:电渗压 | 令 $\Delta\varphi = 0$:流致电流 |
| 测的是「憋出多大的差」(力/力) | 测的是「两个流的比」(流/流) |
先看这一组:把一条通路堵死,另一个力就会自己长出来。
效应(1):流致电势——推水,憋出电压
令 $J_e = 0$(电路开路,不许电流通过),代入 (8.4.2):
$$0 = L_{11}\frac{\Delta\varphi}{T} + L_{12}\frac{\Delta p}{T}$$
两边的 $T$ 约掉,解出电势差与压强差的比:
$$\Big(\frac{\Delta\varphi}{\Delta p}\Big)_{J_e = 0} = -\frac{L_{12}}{L_{11}}$$
(8.4.5)
此时 $J_e = 0$,但 $\Delta\varphi \ne 0$,是由于压力差产生的物质流而引起的电位差。
这就是第 1 段开头那个现象:你只是推水,却量到了电压。之所以要「开路」,正是因为电压表本身几乎不让电流通过——所以实验上测到的,恰好就是 $J_e = 0$ 这个条件下的值。
效应(4):电渗压——加电压,憋出压强
反过来,令 $J_m = 0$(把液体的通路堵死,不许液体净流动),代入 (8.4.3):
$$0 = L_{21}\frac{\Delta\varphi}{T} + L_{22}\frac{\Delta p}{T}$$
$$\Big(\frac{\Delta p}{\Delta\varphi}\Big)_{J_m = 0} = -\frac{L_{21}}{L_{22}}$$
(8.4.8)
此时物质流为零,即 $J_m = 0$,压力差是由单位电势引起的。加上电压,液体想走却走不了,于是在堵头处憋出一个压强差。
图 2 这一组的共同套路:堵死一条路,另一个力就会自己长出来,长到刚好把这条路上的流顶成零为止。
两个效应的结构完全对称,值得并排记:
(8.4.5) 堵电路 $\Rightarrow$ 分母是 $L_{11}$(电那一行的对角元);
(8.4.8) 堵液路 $\Rightarrow$ 分母是 $L_{22}$(物质那一行的对角元)。
规律:你堵谁,被堵那个流的对角元就是分母,同一行的交叉元就是分子,前面带负号。(这条口诀只管这一组「堵流」的两个效应;下一段那两个是流与流之比,没有负号,套上去会错。)
负号是怎么来的?不是「力反抗力」——那两个量连量纲都不同,谈不上谁反抗谁。真正在互相抵消的是两股流:以 (8.4.5) 为例,压强差通过交叉项 $L_{12}$ 驱动出一股电流,憋起来的电势差再通过对角项 $L_{11}$ 驱动出一股反向电流,两股恰好抵消,净电流才为零。负号记录的是「这两股流必须反向」。
顺带一句:$L_{12}$ 本身的正负是材料性质,不是普适的。孔壁带负电还是带正电,会让整个比值反号。换句话说,这个符号本身就携带着材料的信息——界面化学正是反过来利用这一点,用测到的流动电势去反推孔壁的带电状况。
起点 · 主干的式 (8.4.5),原样照抄
$$\Big(\frac{\Delta\varphi}{\Delta p}\Big)_{J_e = 0} = -\frac{L_{12}}{L_{11}}$$
先看这个比值的单位,它本身就很有意思:
$$\Big[\frac{\Delta\varphi}{\Delta p}\Big] = \frac{\mathrm{V}}{\mathrm{Pa}} = \frac{\mathrm{J/C}}{\mathrm{N/m^2}} = \frac{\mathrm{N\cdot m/C}}{\mathrm{N/m^2}} = \frac{\mathrm{m^3}}{\mathrm{C}}$$
立方米每库仑。请记住这个单位,第 7 段会用它做一件漂亮的事。
代一组实际数量级。电解质溶液流过玻璃多孔塞,实测的流动电势系数随溶液浓度变化很大,大致在 $10^{-8}\sim10^{-6}\ \mathrm{m^3/C}$ 之间——溶液越稀、电导越低,这个系数越大(所以要把信号做大,就用稀溶液)。这里取偏低的一端(约 $0.1\ \mathrm{mol/L}$ 的量级):
$$\Big|\frac{\Delta\varphi}{\Delta p}\Big| \approx 2\times10^{-8}\ \mathrm{m^3/C}$$
加 $1\ \mathrm{bar}$(约一个大气压)的压差,即 $\Delta p = 10^{5}\ \mathrm{Pa}$:
$$|\Delta\varphi| \approx 2\times10^{-8}\times10^{5} = 2\times10^{-3}\ \mathrm{V} = 2\ \mathrm{mV}$$
几毫伏。不算大,但用普通电压表完全测得到——这也是为什么这个效应在十九世纪就被发现了。
它在工程上是个真麻烦,不只是实验室现象:输油管道里的油高速流过管壁,同样会积累流动电势。油是绝缘的,电荷跑不掉,越积越高,最后放电打火——这是油罐车和输油管道着火的一个真实原因。
所以行业规范里会限制流速,并要求良好接地与跨接。知道了 (8.4.5),就知道接地在这个框架里改的是什么:不是材料,而是边界条件——把「开路、$J_e=0$」换成「允许 $J_e \ne 0$」。
不过要说准:接地和跨接主要防的是孤立导体(罐车、法兰、浮顶)积起电位后打火;至于绝缘油自己体相里的电荷,接地帮不上忙,只能靠液体自身的电导慢慢弛豫掉(时间常数约为 $\varepsilon/\sigma$)。
顺带一提,行业里还有一条独立的措施是加抗静电剂,那才是真正去动材料本身(提高液体电导率,也就是提高 $L_{11}$,让积累起来的电荷更快弛豫掉)。两条措施一条改边界、一条改材料,正好对应这个式子里的两个位置。
起点 · 主干里的这个效应
$$\Big(\frac{\Delta\varphi}{\Delta p}\Big)_{J_e = 0} = -\frac{L_{12}}{L_{11}}$$
在界面化学里,这个式子被反过来用:不是拿它算电势,而是拿测到的电势去反推材料表面的性质。
要反推的那个量叫 zeta 电位($\zeta$)——大致是双电层里「能跟着液体一起走的那部分」与「固定在壁上的那部分」之间的电势差。它是判断胶体稳不稳定的核心指标:
- $|\zeta|$ 大 $\Rightarrow$ 颗粒之间互相排斥强 $\Rightarrow$ 悬浮液稳定,不结块;
- $|\zeta|$ 小 $\Rightarrow$ 排斥不过范德瓦耳斯吸引 $\Rightarrow$ 颗粒絮凝、沉降。
这件事在很多行业里是硬指标:涂料会不会分层、牛奶会不会结块、水厂加多少絮凝剂、药物微粒能不能久放——判据都是 $\zeta$。
而 $\zeta$ 没法直接测(你没法把一根探针插进几纳米厚的双电层)。于是就用宏观的电动效应去反推:测流动电势,或者测电渗流量,代进理论公式解出 $\zeta$。
所以在那边它解决了什么问题:它把一个纳米尺度、测不到的界面性质,变成了一个用压力表和电压表就能测的宏观量。
这也再次说明本节这套关系的价值:正因为四个效应之间有确定的关系,才可以挑其中最好测的那一个去做实验,再换算到其余的。下一段的电渗,正是实践中更常用的那一个。
第四步 · 四个效应(下)· 归零一个力
6. 只留一个力,看两个流的比
另一组的套路正好相反:不堵任何一条路,而是把其中一个力调成零,让系统只受另一个力驱动,然后看两个流的比值。
效应(2):电渗效应——通电,水就流
令 $\Delta p = 0$(两边压强相同),(8.4.2)(8.4.3) 各只剩一项:
$$J_e = L_{11}\frac{\Delta\varphi}{T},\qquad J_m = L_{21}\frac{\Delta\varphi}{T}$$
两式相除,$\dfrac{\Delta\varphi}{T}$ 整个约掉:
$$\Big(\frac{J_m}{J_e}\Big)_{\Delta p = 0} = \frac{L_{21}}{L_{11}}$$
(8.4.6)
此时压差为零,即 $\Delta p = 0$,这个比值是由于单位电流引起的物质流。
大白话:每通过 1 库仑电荷,能带动多少体积的液体。这就是「通电推水」,也就是电渗泵的工作原理。
效应(3):流致电流——推水,水就带电荷走
令 $\Delta\varphi = 0$(两边电势相同,比如把两个电极短路),同样各剩一项:
$$J_e = L_{12}\frac{\Delta p}{T},\qquad J_m = L_{22}\frac{\Delta p}{T}$$
$$\Big(\frac{J_e}{J_m}\Big)_{\Delta\varphi = 0} = \frac{L_{12}}{L_{22}}$$
(8.4.7)
此时电势差为零,即 $\Delta\varphi = 0$,但有电流,是由单位体积流引起的电流。
请把效应(1)和效应(3)并排看一眼,它们非常容易混。两个都是「推水」,但
电路的状态完全不同:
| (1)流致电势 | (3)流致电流 |
| 电路 | 开路(接电压表) | 短路(把两电极接通) |
| 条件 | $J_e = 0$ | $\Delta\varphi = 0$ |
| 你测到 | 电压 $\Delta\varphi$ | 电流 $J_e$ |
同一个装置、同一个动作,只因为电路接法不同,测到的就是两个不同的效应。「条件」那一行是关键:开路是「不许有电流」,短路是「不许有电压」——恰好互补。
起点 · 主干的式 (8.4.6),原样照抄
$$\Big(\frac{J_m}{J_e}\Big)_{\Delta p = 0} = \frac{L_{21}}{L_{11}}$$
先看单位。$J_m$ 是体积流($\mathrm{m^3/s}$),$J_e$ 是电流($\mathrm{C/s}$),所以
$$\Big[\frac{J_m}{J_e}\Big] = \frac{\mathrm{m^3/s}}{\mathrm{C/s}} = \frac{\mathrm{m^3}}{\mathrm{C}}$$
又是立方米每库仑——和上一段流动电势的单位一模一样。这不是巧合,第 7 段揭晓。
代数字。取和上一段同一个量级 $\dfrac{J_m}{J_e} = 2\times10^{-8}\ \mathrm{m^3/C}$,通 1 毫安电流:
$$\begin{aligned}J_m &= 2\times10^{-8}\ \mathrm{m^3/C}\times 10^{-3}\ \mathrm{C/s}\\[2pt] &= 2\times10^{-11}\ \mathrm{m^3/s}\end{aligned}$$
换成好懂的单位($1\ \mathrm{m^3} = 10^{12}\ \mathrm{nL}$):
$$J_m = 20\ \mathrm{nL/s} = 1.2\ \mathrm{\mu L/min}$$
每分钟一微升多一点。听起来微不足道,但这正是微流控芯片需要的流量档位——在一根头发丝粗细的通道里,这已经是很可观的流速了。
电渗泵有一个机械泵永远比不了的优点:它没有任何运动部件。没有活塞、没有叶轮、没有阀门,只有两个电极和一段多孔材料。
而且流速正比于电流,调电压就能精确调流量,还不产生脉动(机械泵的活塞往复会让流量一跳一跳)。所以芯片实验室、微型燃料电池的冷却、乃至一些植入式给药装置,用的都是它。
代价是:只能泵导电的液体(要有离子才有双电层),而且耗电不低——毕竟它是靠拖拽那薄薄一层电荷,把整团水「捎带」着走的。
起点 · 主干里的电渗效应
$$\Big(\frac{J_m}{J_e}\Big)_{\Delta p = 0} = \frac{L_{21}}{L_{11}}$$
不是一回事,但它们是同一枚硬币的两面,都源于双电层加电场。
| 电渗 electroosmosis | 电泳 electrophoresis |
| 谁不动 | 固体(孔壁)不动 | 液体整体不动 |
| 谁在动 | 液体整团流过去 | 颗粒 / 分子在液体里走 |
| 典型用途 | 电渗泵、微流控输液 | DNA 电泳、蛋白质分离 |
「谁不动」这一行是分辨它们的唯一标准。物理机制上,两者其实是同一个过程在不同参考系下的样子:都是「带电的那一层被电场推动,靠黏滞把邻居带着走」。区别只在于你把哪一方钉死。
错误的理解:「电泳就是带电颗粒被电场推着走,所以跟本节的交叉效应没关系。」
为什么错:这样理解会漏掉一件在实验室里天天要处理的事——做电泳时,电渗同时也在发生。毛细管电泳里,管壁的双电层会驱动整根管子里的液体一起流动,所有样品都被这股流「托着」跑。
结果就是:你测到的迁移速度 = 颗粒自己的电泳速度 + 背景电渗速度。不把电渗扣掉,测出的电泳迁移率就是错的。所以毛细管电泳要么专门涂层抑制电渗,要么加一个中性标记物去标定它。
一个记忆钩子:「渗」字说的是液体穿过多孔介质(渗透、渗流);「泳」字说的是东西在液体里游。看字就知道谁在动。
第五步 · 收网 · 这一节真正的成果
7. 四个效应,其实只有两个
四个效应都算出来了。把它们摆在一起:
| 效应 | 固定什么 | 结果 | 式号 |
| (1)流致电势 | $J_e = 0$ | $\left(\dfrac{\Delta\varphi}{\Delta p}\right)_{J_e=0} = -\dfrac{L_{12}}{L_{11}}$ | (8.4.5) |
| (2)电渗效应 | $\Delta p = 0$ | $\left(\dfrac{J_m}{J_e}\right)_{\Delta p=0} = \dfrac{L_{21}}{L_{11}}$ | (8.4.6) |
| (3)流致电流 | $\Delta\varphi = 0$ | $\left(\dfrac{J_e}{J_m}\right)_{\Delta\varphi=0} = \dfrac{L_{12}}{L_{22}}$ | (8.4.7) |
| (4)电渗压 | $J_m = 0$ | $\left(\dfrac{\Delta p}{\Delta\varphi}\right)_{J_m=0} = -\dfrac{L_{21}}{L_{22}}$ | (8.4.8) |
如果没有昂萨格关系,这四个效应就是独立的——四个不同的实验,测四个互不相干的数,谁也管不着谁。
但 (8.4.4) 说 $L_{12} = L_{21}$。把它代进去,(1) 和 (2) 的右边就只差一个负号,(4) 和 (3) 也一样:
$$\Big(\frac{\Delta\varphi}{\Delta p}\Big)_{J_e=0} = -\Big(\frac{J_m}{J_e}\Big)_{\Delta p=0},\qquad \Big(\frac{\Delta p}{\Delta\varphi}\Big)_{J_m=0} = -\Big(\frac{J_e}{J_m}\Big)_{\Delta\varphi=0}$$
(8.4.9)
于是4 个不可逆过程只有两个是独立的。课本管这叫「不可逆过程热力学得到的典型结果」。
图 3 $L_{12}=L_{21}$ 把四个效应两两配对。横向相连的两个,测一个就知道另一个。
请体会 (8.4.9) 左边那一条到底有多离谱。它把两个看起来毫无关系的实验捆在了一起:
左边:不通电,拿压力表和电压表——推水,量它憋出多少电压。
右边:不加压,拿电源和流量计——通电,量它推出多少液体。
一个实验里根本没有电源,另一个实验里根本没有压差。两台仪器、两个操作、两种物理直觉。可它们测出来的数,除了一个负号必须严格相等。
这不是近似、不是经验规律,是热力学保证的。而且它对任何多孔材料都成立——玻璃、陶瓷、生物膜、页岩,全都一样。
起点 · 主干的 (8.4.9) 左边一条,原样照抄
$$\Big(\frac{\Delta\varphi}{\Delta p}\Big)_{J_e=0} = -\Big(\frac{J_m}{J_e}\Big)_{\Delta p=0}$$
先把第 5、6 段留下的那个伏笔收掉:为什么两边单位相同?
| 量 | 单位 |
| 左边 | $\dfrac{\Delta\varphi}{\Delta p}$ | $\dfrac{\mathrm{V}}{\mathrm{Pa}} = \dfrac{\mathrm{J/C}}{\mathrm{N/m^2}} = \dfrac{\mathrm{m^3}}{\mathrm{C}}$ |
| 右边 | $\dfrac{J_m}{J_e}$ | $\dfrac{\mathrm{m^3/s}}{\mathrm{C/s}} = \dfrac{\mathrm{m^3}}{\mathrm{C}}$ |
两边都是 $\mathrm{m^3/C}$。这不是巧合——一个等式如果两边量纲不同,它根本不可能成立。所以量纲相同是 (8.4.9) 能存在的必要条件,也是你事后检查有没有抄错的第一道关。
现在做一次真正的预言。假设你手上只有一台电源和一个流量计,做了电渗实验,测得:
$$\Big(\frac{J_m}{J_e}\Big)_{\Delta p=0} = 2\times10^{-8}\ \mathrm{m^3/C}$$
(对应第 6 段那个数:1 mA 推出 20 nL/s。)
那么不必再做任何实验,(8.4.9) 直接给出流动电势系数:
$$\Big(\frac{\Delta\varphi}{\Delta p}\Big)_{J_e=0} = -2\times10^{-8}\ \mathrm{m^3/C}$$
换算成实验室里读得到的话:每加一个大气压,两端会出现约 2 mV 的电压(负号只表示两者反号,具体哪一端为正取决于孔壁带什么电)。
$$|\Delta\varphi| = 2\times10^{-8}\times10^{5} = 2\ \mathrm{mV}$$
这正是第 5 段实例里那个数——两段的数字本来就是同一个数,现在你知道为什么了。
反过来用,它还是一把尺子。如果你把两个实验都做了,理论说结果必须相等;那么两次测量对不对得上,就成了一次对整套装置的独立校验。
对不上,说明装置有系统误差(电极极化了?塞子堵了?温度没控住?)。一条理论关系,反过来成了实验的标尺——这是昂萨格关系在实验室里最常被使用的方式,甚至比拿它省实验更常见。
起点 · 主干的结论
$$4 \text{ 个效应} \;\xrightarrow{\;L_{12}=L_{21}\;}\; 2 \text{ 个独立}$$
先数系数。矩阵 $L$ 有 4 个元素;昂萨格关系给出 1 个约束,独立的只剩 3 个:$L_{11}$、$L_{22}$,以及公共的 $L_{12}=L_{21}$。
再数效应。四个效应都是比值,各自只含两个系数的商:
| 效应 | 比值 |
| (1) 与 (2) | 都是 $\dfrac{L_{12}}{L_{11}}$(差一个负号) |
| (3) 与 (4) | 都是 $\dfrac{L_{12}}{L_{22}}$(差一个负号) |
所以四个效应实际上只携带两个独立的信息:$\dfrac{L_{12}}{L_{11}}$ 和 $\dfrac{L_{12}}{L_{22}}$。「哪两个」没有唯一答案——你可以挑 (1) 和 (3)、也可以挑 (2) 和 (4),只要不是同一对里的两个就行。
实践中怎么挑?挑最好测的那两个。比如挑 (2) 电渗(只需电源和流量计,信号大)和 (4) 电渗压(只需电源和压力表),把 (1)(3) 全部换算出来。
顺带说清一件事:为什么四个效应凑不齐三个独立系数?因为它们全是
比值,比值天然丢掉了「整体大小」这一维。
要把 $L_{11}$、$L_{22}$、$L_{12}$ 三个都定下来,还得再测
一个绝对量——比如量一下这块塞子的电导,或者量一下它的渗透率。
但这里有个陷阱:测这两个量时,必须说清另一条路是通的还是堵的。
| 测什么 | 另一条路 | 实际测到的是 |
| 电导 | 液路通($\Delta p=0$) | $L_{11}$ ✓ |
| 电导 | 液路堵($J_m=0$) | $L_{11}-\dfrac{L_{12}^2}{L_{22}}$ |
| 渗透率 | 电路短路($\Delta\varphi=0$) | $L_{22}$ ✓ |
| 渗透率 | 电路开路($J_e=0$) | $L_{22}-\dfrac{L_{12}^2}{L_{11}}$ |
后两行都比前两行
小,小掉的正是交叉效应那一块。其中最后一行有个名字叫
电黏滞效应:开路时量到的渗透率偏低,
看起来就像液体的黏度变大了——其实液体没变,只是它一边流一边在憋电压,而憋出来的电压又反过来把它往回推。
这和第 2 份末尾说过的是同一件事:说「热导率」之前必须先交代另一个流被固定成什么。这里是同一条教训的第二次出现。
「两个比值 + 一个绝对量 = 三个独立系数」,账就平了。
还有一点容易被误读:(8.4.9) 那两条式子,彼此并不是两条新信息。每一条单独拿出来都
等价于 $L_{12}=L_{21}$ 这一句话——它们只是同一个事实的两种测法(一种从「堵电路 vs 不加压」测,另一种从「堵液路 vs 短路」测)。所以做完其中一条的验证,另一条不会再提供新的检验力。
本份小结
这一份的逻辑链条是:装置是两容器夹一块多孔隔板,温度相同、压强不等、电势不等(这是能凑出「两个力」的最简系统)→ 先把压强差翻译成化学势差,靠的是 $\left(\frac{\partial\mu}{\partial p}\right)_T = v$,近平衡下写成 $\Delta\mu = v\Delta p$ → 于是局面与第 2 份的「情形三」相同,熵产生率直接搬来,代入后得 (8.4.1),右边恰好是两项「力×流」→ 读出两个流 $J_e=-q\dot N^{\mathrm I}$、$J_m=-v\dot N^{\mathrm I}$ 与两个力 $\frac{\Delta\varphi}{T}$、$\frac{\Delta p}{T}$,写出唯象方程 (8.4.2)(8.4.3),此时且仅此时 $L_{12}=L_{21}$ → 依次令 $J_e=0$、$\Delta p=0$、$\Delta\varphi=0$、$J_m=0$,得到四个效应 (8.4.5)–(8.4.8) → 代入昂萨格关系,得 (8.4.9):四个不可逆过程只有两个是独立的。
下一份(8.5 热电效应)会做同一件事,但难度上一个台阶:那里温度不再相同,力换成了温度梯度和电势梯度,而且要处理的是三个热电系数(泽贝克、佩尔捷、汤姆孙)之间的关系——最后同样由昂萨格关系收网,给出开尔文关系式。