固体表面自由能计算方法综述：从经典理论到实验质量控制

基于接触角测量的理论、方法与系统比较

摘要

固体表面自由能（Surface Free Energy, SFE）是表征固体表面物理化学性质的核心参数，在复合材料、涂层、生物医学、矿物加工等领域具有重要应用价值。本文系统综述了基于接触角测量的表面自由能计算方法，从 Young 方程出发，详细阐述了各类方法的物理基础、数学表述、适用范围及局限性。方法分为两大类：第一类假设固-液界面张力仅取决于固、液表面张力的绝对值，包括 Zisman 临界表面张力法、Neumann 状态方程法和 Chibowski 滞后角法；第二类假设固-液界面张力取决于表面张力中各组分（色散、极性、酸碱等）的贡献，包括 Fowkes 法、Owens-Wendt 法、Wu 调和平均法以及 van Oss-Chaudhury-Good 的 LW-AB 法。此外，本文还专门讨论了高能固体表面自由能的测定方法（两液法）和多孔、粉末、纤维材料的 Washburn 渗透法。通过对微生物细胞表面（147 株菌株）的比较分析，揭示了各方法的内部一致性及其适用范围。进一步，本文针对高能固体（铜、硅、金、铂、不锈钢、铝合金、铝）讨论了各方法的适用性，指出状态方程法在金属体系中的根本局限。在此基础上，本文提出三项实验质量控制观点：(1) 探针液体滴到固体表面后，需采用力学铂金板表面张力及表面清洁度测量仪验证其表面张力是否变化，以判断固体溶出物或表面污染对液体纯度的干扰；(2) 当液滴体积大于 0.5 μL 时，重力、表面粗糙度和结构易导致钉扎，需采用非接触式进液系统形成润湿-去润湿后的稳态接触角；(3) 表面结构影响需通过 3D 形貌测试表征，表面化学多样性可通过可视 UV 分析识别。本文旨在为材料科学、界面化学和生物工程领域的研究者提供从理论选择到实验质量控制的完整指南。

关键词：表面自由能；接触角；Young 方程；LW-AB 法；状态方程；润湿性；探针液体纯度；稳态接触角；表面结构

1. 引言

1805 年，Thomas Young 提出了描述液滴在固体表面平衡状态的力学方程，奠定了表面热力学的理论基础。然而，Young 方程包含两个未知量——固体表面张力 γ_SV 和固-液界面张力 γ_SL，仅通过单次接触角测量无法求解。因此，在过去的两个多世纪中，研究者们发展了大量理论方法来解决这一“一个方程、两个未知数"的经典难题，形成了表面自由能计算的一个丰富而复杂的理论体系。

根据对固-液界面张力性质的不同物理假设，现有方法可大致分为两大流派。第一类方法（常称为单参数法或状态方程法）认为，固-液界面张力 γ_SL 仅是固体表面张力 γ_SV 和液体表面张力 γ_LV 的单一函数，即存在一个普适的状态方程 γ_SL = f(γ_SV, γ_LV)。这类方法的代表性工作包括 Zisman 的临界表面张力法、Neumann 及其合作者的状态方程法以及 Chibowski 基于接触角滞后的计算方法。这类方法的核心优势在于操作简便——理论上仅需一种探测液体即可计算固体的表面自由能。然而，它们的普适性和物理基础一直是学术界激烈争论的焦点。

第二类方法（常称为多参数组分法或表面张力分量法）则认为，表面张力是不同类型分子间相互作用力的加和，固-液界面张力取决于两相中各种相互作用分量（色散、极性、氢键、酸碱等）的匹配情况。Fowkes 在 1960 年代首先提出表面张力可分为色散分量和非色散分量，并引入了几何平均结合规则。随后 Owens 和 Wendt、Wu 以及 van Oss、Chaudhury 和 Good 等人进一步发展了这一思想，形成了 Owens-Wendt 几何平均法、Wu 调和平均法以及 Lifshitz-van der Waals/酸碱（LW-AB）法。这类方法通常需要两种或三种探测液体，但能提供比单参数法丰富得多的表面化学信息（如极性分量、电子给体/受体参数），因此在材料科学、生物界面、胶体与界面化学等领域获得了广泛应用。

此外，针对不同类型的固体材料——低能聚合物、高能金属/氧化物/云母、微生物细胞以及多孔/粉末/纤维材料——研究者还发展了专门的测定技术，如两液法（用于高能固体）和 Washburn 渗透法（用于多孔材料）。

在表面自由能计算领域，有五项经典工作奠定了现代方法的基础。Hansen 对聚合物表面能测量常用方法作了简要概述，重点介绍了 Young 方程、Zisman 图、Fowkes 理论、几何平均与调和平均、LW-AB 法以及状态方程法。Zenkiewicz 系统分析了基于接触角测量的固体 SFE 计算方法，包括 Young 方程、状态方程、表面张力分量划分以及多孔/粉末材料的 Washburn 方法。Schultz、Tsutsumi 和 Donnet 发展了测定高能固体（云母）表面自由能的两液法，分别给出了色散分量和极性分量的测定结果，并与剥离实验对比验证。Sharma 和 Rao 利用 147 株微生物的接触角数据对 Fowkes 法、状态方程法、几何平均法和 LW-AB 法进行了系统比较，评估了各方法的内部一致性，并发现电子给体参数 γ⁻ 可有效区分革兰氏阴阳性菌。Janczuk、Zdziennicka 和 Szymczyk 对界面张力、Gibbs 自由能和接触角的关系进行了全面评述，详细分析了 Neumann 状态方程、Chibowski 滞后角法、Owens-Wendt 法、Wu 法、LW-AB 法以及两液法的理论基础和实验适用性，并给出了大量固体表面张力的计算实例。

本综述按照以下结构展开：第 2 节详细介绍 Young 方程及其热力学基础，并系统阐述三项实验质量控制观点；第 3 节阐述第一类方法（Zisman 法、Neumann 状态方程法、Chibowski 滞后角法），并专门讨论金属体系的适用性；第 4 节阐述第二类方法（Fowkes 法、Owens-Wendt 法、Wu 法、LW-AB 法）；第 5 节介绍两液法；第 6 节介绍 Washburn 渗透法；第 7 节以微生物细胞为例进行比较；第 8 节为总结与方法选择指南；第 9 节为展望。

2. 理论基础：Young 方程及其热力学拓展

2.1 Young 方程的标准形式与力学平衡

Young 方程描述了理想光滑、均匀、刚性、化学惰性且无吸附膜的固体表面上，平衡液滴的力学平衡条件。考虑一个液滴置于固体表面上，三相接触线（固-液-气）处存在三个界面张力：

• γ_SV：固-气界面张力，倾向于收缩固-气界面；

• γ_SL：固-液界面张力，倾向于收缩固-液界面；

• γ_LV：液-气界面张力，沿切线方向作用于液滴轮廓。

在接触线处沿水平方向建立力学平衡：

γ_SV = γ_SL + γ_LV cos θ  (1)

式中 θ 为接触角，定义为液滴在固-气界面一侧的切线角。该方程也可改写为：

γ_LV cos θ = γ_SV – γ_SL  (2)

方程 (1) 的推导假设固体表面是理想化的平坦、刚性和化学均匀的，且忽略了液体蒸气在固体表面的吸附以及重力对液滴形状的影响（对于小液滴可忽略重力）。

2.2 表面张力与表面自由能的热力学等价性

从热力学角度，表面张力 γ 等于在恒定温度、恒定组成和恒定表面积条件下，增加单位表面积所需的 Helmholtz 自由能变化（对于恒容过程）或 Gibbs 自由能变化（对于恒压过程）：

γ = (∂F^s/∂A){T, n_i} = (∂G^s/∂A){T, n_i}  (3)

因此，表面张力在数值上等同于表面自由能。在文献中这两个术语经常互换使用，但严格说来，“表面张力"强调力学概念（单位长度上的力），而“表面自由能"强调热力学概念（单位面积上的能量）。对于液体，两者几乎等价；对于固体，由于弹性约束，通常使用“表面自由能"更为恰当。

2.3 铺展压（Spreading Pressure）的修正

Bangham 和 Razouk（1937）指出，液体蒸气在固体表面的物理吸附会形成吸附膜，降低固-气界面张力。他们引入铺展压 π_e 对 Young 方程进行修正：

γ_SV^0 – π_e – γ_SL = γ_LV cos θ  (4)

其中 γ_SV^0 为无吸附时的“清洁"固体表面张力（通常简记为 γ_SV），π_e = γ_SV^0 – γ_SV，γ_SV 为有吸附膜存在时的固-气界面张力。铺展压可通过吸附等温线积分求得：

π_e = (R T / (V_m A^s)) ∫_0^P (v / p) dp  (5)

式中 V_m 为吸附质的摩尔体积，A^s 为固体比表面积，v 为吸附量（体积），p 为蒸气压力。然而，π_e 的精确测定非常困难，需要对固体进行细致的吸附实验。因此，在实际表面自由能计算中，常根据特定准则判断 π_e 是否可以忽略。Fowkes、Good、Neumann 以及 van Oss 等人普遍认为：当液体表面张力大于固体表面张力（γ_LV > γ_SV）时，液体蒸气在固体表面的吸附不会显著降低固体表面张力（即 π_e ≈ 0）。这一准则对低能聚合物（如 PTFE，γ_SV ≈ 18–20 mN/m）尤其适用，因为水（72.8 mN/m）、甘油（64 mN/m）、甲酰胺（58 mN/m）和大多数有机液体的表面张力均大于聚合物表面张力。因此，在绝大多数聚合物润湿研究中，可直接使用简化形式的 Young 方程（即忽略 π_e）。但对于高能固体（如石英、云母、金属），其清洁表面张力通常远高于液体表面张力，此时 π_e 不可忽略，必须采用两液法（参见第 5 节）或直接测量铺展压。

更一般的修正还包括固-液界面处的吸附膜压力 π_SL。Bialopiotrowicz 和 Janczuk（1990）提出：

γ_SV – π_SV – γ_SL + π_SL = γ_LV cos θ  (6)

他们认为在有些体系中 π_SV = π_SL，或者两者均为零，或者两者均大于零但不相等。例如，对于表面活性剂溶液在聚合物上的润湿，即使 γ_LV > γ_SV，表面活性剂分子仍可能沿三相线迁移到固体表面形成吸附层，导致 π_SV > 0。Zdziennicka 和 Janczuk（2020）通过大量实验发现，对于非离子表面活性剂（Triton 系列、烷基糖苷）在 PTFE 上的润湿，铺展压近似等于水表面张力与溶液表面张力差值的一半。

2.4 接触角的类型与滞后现象

实际测量中可区分三种接触角，它们具有不同的物理意义和实际应用：

• 前进角 θ_A：将液滴缓慢增加体积，三相线向前移动（即液滴扩展）时测得的接触角。它代表液滴在干燥清洁固体表面上铺展时的动态接触角，与固-液界面形成过程相关。

• 后退角 θ_R：将液滴缓慢抽吸减少体积，三相线向后回缩（即液滴收缩）时测得的接触角。它代表固-液界面破坏、固体表面被液滴剥离时的动态接触角。

• 平衡角 θ_Y：Young 方程定义的理想平衡接触角，介于 θ_A 和 θ_R 之间（θ_R ≤ θ_Y ≤ θ_A）。

接触角滞后（θ_A > θ_R）的常见原因包括：

• 表面粗糙度（Wenzel 效应）：粗糙度放大了本征润湿性；

• 化学不均匀性（Cassie 效应）：不同化学位点的表面能差异；

• 表面分子重排：高分子链或表面活性剂分子在固-液接触时发生取向变化；

• 吸附膜的形成：液体蒸气在固体表面形成前驱膜；

• 钉扎效应（contact line pinning）：表面缺陷或杂质对三相线的物理钉扎。

Chibowski（2003）利用接触角滞后信息发展了固体表面自由能的计算方法（参见第 3.3 节）。他指出，θ_A 和 θ_R 的测量可以为计算 γ_SV 提供独立于探测液体类型之外的信息。

2.5 理想固体与实际固体的差异

大多数表面自由能计算方法（尤其是 Zisman 法和 Neumann 法）都假设固体表面是理想光滑、化学均匀且惰性的。然而，实际聚合物表面存在不同程度的粗糙度（R_a 通常在纳米到微米级）和化学异质性（如共聚、共混、添加剂析出、氧化等）。这会导致测量得到的表观接触角偏离本征 Young 接触角，进而影响表面自由能的计算精度。因此，在可能的情况下，应使用原子力显微镜（AFM）表征表面形貌，并使用 X 射线光电子能谱（XPS）表征表面化学组成，以辅助解释接触角数据。对于粗糙表面，可采用 Wenzel 方程（cos θ_apparent = r cos θ_Y，r 为粗糙度因子）或 Cassie-Baxter 方程（适用于复合接触）进行修正，但这类修正会引入新的未知参数，使得表面自由能的求解更加复杂。

2.6 实验质量控制的三项核心观点

在进入具体计算方法之前，必须强调：所有表面自由能计算结果的可靠性，首先取决于接触角测量本身的准确性。文献中大量方法争议（尤其是状态方程法在高能固体上的失败）往往源于实验环节的疏忽，而非理论本身的缺陷。基于此，本文提出三项贯穿所有方法的实验质量控制观点。

2.6.1 观点一：探针液体表面张力的原位验证——力学铂金板法

所有表面自由能计算模型均要求已知探针液体的准确表面张力（γ_LV）及其分量（γ_L^d、γ_L^p、γ_L^、γ_L^+、γ_L^-）。这些参数通常从文献中查得（如水在 20°C 时为 72.8 mN/m）。然而，当液滴滴加到固体表面后，以下两种情形会改变液体的实际表面张力：

• 固体溶出物：亲水性固体（如某些盐、氧化物、生物材料）可能向液滴中溶解或释放离子/小分子，改变液体组成。例如，某些矿物粉末与水接触后会释放金属离子，改变水的表面张力。生物细胞表面可能释放多糖或蛋白质，同样会污染液滴。

• 表面污染物：固体表面的有机物或油脂可能被液滴萃取，降低表面张力。即使经过清洗的金属表面，在空气中暴露数分钟后也会吸附一层碳氢化合物，这些污染物可被极性液滴（如水）部分溶解。

验证方法：在接触角测量的同时或单独实验中，将滴加到固体表面并保持相同接触时间（例如 30–60 秒）后的液滴，用微量针筒吸起，采用力学铂金板表面张力及表面清洁度测量仪测定其表面张力，并与原始探针液体的表面张力对比。若偏差超过 ±0.5 mN/m（水的典型测量精度），则表明液体已被污染，接触角数据不可用于表面自由能计算，需对固体进行更严格的清洗或更换探针液体。这一验证应作为所有表面自由能测定实验的必经流程。

2.6.2 观点二：微小液滴与稳态接触角——非接触式进液系统

传统接触角测量通常使用体积为 2–5 μL 的液滴，这是基于普通微注射器的手动或半自动滴液方式。然而，当液滴体积大于 0.5 μL 时，重力（Bond 数 Bo = ρgR²/γ 增大）和表面粗糙度/结构会共同作用，导致三相线钉扎（pinning）。钉扎使得前进角和后退角测量值失真，且无法通过简单的 Young 方程描述。具体来说：

• 重力使液滴形状偏离理想球冠，接触角的几何计算产生误差；

• 液滴体积较大时，三相线在粗糙表面上的运动受到能垒阻碍，液滴呈现“粘滑"行为；

• 前进角和后退角的差值被人为放大，而平衡角难以确定。

解决方案：采用上海梭伦MicroDrop®非接触式进液系统（如压电式或热泡式微喷射技术）产生体积 ≤ 0.5 μL 的微小液滴，并利用润湿-去润湿循环获得稳态接触角。具体步骤：

1. 以非接触方式将亚微升液滴沉积于固体表面，避免针尖对表面的机械干扰；

2. 非接触式高速在材料表面形成液滴过程中，会明显出现润湿与去润湿过程，稳态接触角易形成。

理论依据：当液滴体积足够小时，重力可忽略，液滴形状接近球冠，接触角由表面化学和微观结构决定。此时测得的接触角更接近本征 Young 角。Sharma & Rao（2002）在研究微生物细胞时观察到，即使极小的液滴（约 0.2 μL）仍能获得稳定的平衡阶段（持续 1–3 秒），而较大液滴则立即进入吸收阶段。非接触式进液系统将这一经验系统化，适用于所有固体材料，特别是生物材料、高能固体和粗糙表面。

2.6.3 观点三：表面结构与化学多样性的独立表征——3D 形貌与可视 UV 分析

接触角是固体表面物理结构和化学组成的综合体现。若不区分两者，可能导致表面自由能计算的错误归因。例如，粗糙的亲水表面可能表现为疏水（Cassie 态），而化学均匀的疏水表面经粗糙化后也可能呈现超疏水。因此，仅凭接触角数据无法区分本征润湿性与形貌效应。

结构表征：采用表面 3D 结构测试（如共聚焦显微镜、白光干涉仪或原子力显微镜）获取表面粗糙度参数（Sa、Sq、Ssk、Sku）和分形维数。具体参数含义：

• Sa：算术平均高度，反映平均粗糙度；

• Sq：均方根高度，对尖峰和深谷更敏感；

• Ssk：偏度，正值表示表面以尖峰为主，负值表示以孔洞为主；

• Sku：峰度，反映高度分布的高斯性。

若表面为 Wenzel 态（液体可以填充进粗糙凹陷），需用粗糙度因子 r 修正接触角（cos θ_apparent = r cos θ_Y）；若为 Cassie-Baxter 态（液体下方有空气截留），则需同时知道固-液接触面积分数 f，修正公式为 cos θ_apparent = f (cos θ_Y + 1) – 1。r 和 f 可通过 3D 形貌数据和液体渗透压力估算。

化学多样性表征：采用可视 UV 分析（紫外反射光谱或紫外荧光成像）识别表面化学异质性。不同化学基团（如羟基、羧基、烷基、氟碳基、氨基等）在特定紫外波长下有特征吸收或发射。例如：

• 羟基在 260–280 nm 有弱吸收；

• 芳香族化合物在 250–270 nm 有强吸收；

• 荧光标记的官能团可在紫外激发下发射可见光。

通过 UV 图像可判断表面是否化学均匀。若存在斑点或条纹，表明表面化学组成不均一，则需采用 Cassie 模型或进行局部接触角测量（使用微米级液滴阵列）。对于生物材料，UV 荧光还可用于检测蛋白质、多糖或脂质的分布。

Zenkiewicz（2007）明确指出，表面粗糙度和化学异质性是导致接触角滞后和不同方法结果不一致的主要原因。Janczuk 等（2026）在讨论高能固体（石英）时，强调了表面水膜和羟基分布的影响。将 3D 形貌与 UV 分析纳入实验流程，可使表面自由能计算从“黑箱"回归到物理可解释的层面。

3. 第一类方法：基于 γ_SL = f(γ_SV, γ_LV) 假设的方法

3.1 Zisman 临界表面张力法

Zisman 及其合作者在 1950–1960 年代通过系统实验发现，对于一系列同系物液体（如正构烷烃）在同一种低能固体表面上的接触角，cos θ 与液体表面张力 γ_LV 呈现良好的线性关系。这一经验关系可用下式表示：

cos θ = 1 + b (γ_C – γ_LV)  (7)

当 θ = 0（即 cos θ = 1）时，由上式可得 γ_LV = γ_C。Zisman 将 γ_C 定义为该固体的临界表面张力，并认为 γ_C 是表征固体润湿特征的重要参数。对于 PTFE，用正构烷烃系列测得 γ_C ≈ 18.5 mN/m；对于聚乙烯，γ_C ≈ 31 mN/m；对于聚甲基丙烯酸甲酯（PMMA），γ_C ≈ 39 mN/m。

然而，随着研究的深入，Zisman 法的局限性逐渐显现：

(1) 线性关系并非普适：当使用非极性液体（正构烷烃）在低能非极性固体（PTFE、聚乙烯）上测量时，cos θ vs. γ_LV 确实呈线性。但当使用极性液体（如水、甘油、甲酰胺）或在极性固体（如 PMMA、尼龙、PET）上测量时，cos θ vs. γ_LV 图呈显著曲线（可用二阶多项式拟合），线性外推会得出错误的 γ_C 值。

(2) γ_C 与液体系列有关：Dann（1970）对 9 种聚合物使用不同系列液体（正构烷烃、混合乙二醇、ASTM 系列等）测得的 γ_C 值差异显著。例如，对于聚对苯二甲酸乙二醇酯（PET），用正构烷烃系列得到 γ_C ≈ 43 mN/m，而用乙醇-水系列得到 γ_C ≈ 37 mN/m。

(3) γ_C 不等于 γ_SV：只有当固-液界面张力 γ_SL = 0 且铺展压 π_e = 0 时，γ_C 才等于 γ_SV。在大多数情况下，γ_SL > 0，因此 γ_C < γ_SV。Kitazaki 和 Hata（1972）进一步将 γ_C 分为三类：γ_C^d（用非极性液体测定，对应色散分量）、γ_C^b（用卤代烃和酯类测定）、γ_C^c（用氢键液体测定）。三类值的差异可达 11 mN/m 以上，说明 γ_C 不是固体的本征参数，而是液体系列和分子间相互作用类型的函数。

适用建议：Zisman 法在涂料、油墨、粘合剂等工业领域仍被广泛用作经验性润湿判据：若液体表面张力小于 γ_C，则该液体可铺展在固体表面；若大于 γ_C，则不能铺展。对于工程快速筛选，这一判据简单有效，但研究者应认识到它的经验本质，不应将 γ_C 等同于真实表面自由能。

3.2 Neumann 状态方程法

Neumann 及其合作者（1974–2000）基于热力学推导，提出固-液界面张力 γ_SL 应满足状态方程形式的普适关系式：γ_SL = f(γ_SV, γ_LV)。他们首先从 Girifalco-Good 方程出发：

γ_SL = γ_SV + γ_LV – 2 φ √(γ_SV γ_LV)  (8)

式中 φ 为分子间相互作用参数（0.5 < φ < 1.15），当两相之间仅存在色散相互作用时 φ=1。Neumann 等人对 Bernett 和 Zisman、Fox 和 Zisman、Ellison 和 Zisman 以及 Shafrin 和 Zisman 等人发表的 8 种低能固体上数十种液体的接触角数据进行了系统分析，发现 φ 与 γ_SL 之间存在线性关系：

φ = 1 – α γ_SL  (9)

根据数据拟合，α = 0.0075 (m²/mJ)，即：

φ = 1 – 0.0075 γ_SL  (10)

将式 (10) 代入式 (8) 并整理，得到状态方程的第一种形式：

γ_SL = [√γ_SV – √γ_LV]² / [1 – 0.015 √(γ_SV γ_LV)]  (11)

结合 Young 方程（忽略铺展压），得：

cos θ = –1 + 2 √(γ_SV/γ_LV) / [1 – 0.015 √(γ_SV γ_LV)]  (12)

由于式 (11) 在 √(γ_SV γ_LV) = 1/0.015 ≈ 66.67 时分母为零，失去物理意义，且该方程对 γ_SV 的解多达三个，难以直接应用。因此 Neumann、Li 和 Kwok 在 1990 年代后期对 Berthelot 结合规则进行了修正。修正的核心思想是：当两相的分子间相互作用能（或表面张力）差异很大时，几何平均会高估相互作用强度，需引入指数衰减因子。他们提出：

φ = exp[–β (γ_LV – γ_SV)²]  (13)

其中 β 为经验常数。通过对 15 种光滑均质固体（包括氟聚合物、聚甲基丙烯酸酯、聚苯乙烯等）上大量液体的接触角数据进行最小二乘拟合，他们得到 β = 0.0001247 (m²/mJ)²。代入 Young-Dupré 方程得：

cos θ = –1 + 2 √(γ_SV/γ_LV) exp[–β (γ_LV – γ_SV)²]  (14)

该方程称为 Neumann-Kwok 状态方程。它避免了分母为零的问题，且对于给定的 θ 和 γ_LV，可通过数值迭代收敛确定 γ_SV。

优点：

• 操作极为简便：仅需一种探测液体的接触角数据（通常用水、甲酰胺或二碘I甲烷）即可求出 γ_SV。

• 数学形式简洁，无需预先知道固体表面张力的任何组分信息。

• 对低能非极性固体（如 PTFE、PE）与其他烷烃或非极性液体的体系，计算结果与 Fowkes 法给出的色散分量比较一致。

争议与局限性：

(1) 热力学基础争议：Morrison（1989, 1991）指出，Neumann 等人的相律分析存在缺陷。对于宏观弯曲液滴（半径 > 1 μm），Laplace 压的影响可以忽略，经典 Gibbs 相律（F = C – P + 2）仍然适用。对于固-液-气三相两组分体系（C=2，P=3），自由度 F=1，即仅温度可独立变化，不可能存在 γ_SL = f(γ_SV, γ_LV) 这样的二元函数。Morrison 认为，状态方程的存在性不能从热力学严格证明。Neumann 团队反驳称，当固-液-气三相中存在弯曲界面（接触线附近曲率较大）时，需采用修正相律，自由度可增至 2，因此状态方程可能存在。但这一争论至今未解决。

(2) β 常数非普适：Janczuk 等人对 PTFE-液体和 PET-液体系统的分析表明，即使在同一固体上，不同液体的 φ-γ_SL 线性关系常数 α*（即式 (9) 的斜率）也不相同。对于 PTFE-水体系，α* ≈ 0.0076（接近 Neumann 值）；对于 PTFE-正构烷烃体系，α* 可低至 0.0039；对于 PTFE-二碘I甲烷体系，α* 甚至为负值（意味着 φ 随 γ_SL 增大而增大，与 Neumann 假设相反）。因此，β 常数高度依赖于液体和固体类型，并非普适常数。

(3) 多液体结果不一致：Zdziennicka 等人（2017）的经典实验表明，对于同一 PTFE 样品，用水（θ=111.1°）、甲酰胺（θ=91.8°）和二碘I甲烷（θ=74.7°）代入式 (14) 计算得到的 γ_SV 分别为 15.76、19.23 和 23.95 mN/m，相差超过 8 mN/m。对于 PMMA，用水（θ=74.25°）、甲酰胺（θ=55.56°）和二碘I甲烷（θ=36.37°）计算得到的 γ_SV 分别为 38.54、38.72 和 42.11 mN/m，虽相对接近但仍不一致。这说明状态方程法对探测液体具有敏感性，违背了“γ_SV 应为固体本征性质"的期望。

(4) 极性体系的偏差：对于含强酸碱相互作用的体系（如 PMMA-水，尼龙-水），状态方程法计算的 γ_SV 系统性低于 Owens-Wendt 法和 LW-AB 法的结果，有时甚至低于 γ_LV，导致 γ_AB 为负数（物理上不可能）。这表明状态方程法本质上假设界面相互作用仅有 Lifshitz-van der Waals 类型，未考虑酸碱相互作用的特异性和不对称性。

适用建议：状态方程法仅适用于低能非极性固体（γ_SV < 30 mN/m）与无特异性相互作用的非极性或弱极性液体（如烷烃、二碘I甲烷、α-溴萘）。对于极性固体（PMMA、尼龙、PET）或极性液体（水、甘油），状态方程法给出的 γ_SV 随探测液体强烈变化，不应使用。在快速估算时，建议用二碘I甲烷作为探测液体，因其表面张力较高且色散分量占主导，结果相对稳定。

3.3 Chibowski 滞后角法

Chibowski（2003）利用前进角 θ_A 和后退角 θ_R 的热力学信息，提出了一种独立于探测液体类型的固体表面自由能计算方法。该方法的出发点是：当液滴在固体表面上前进时，固-气界面张力为 γ_SV（裸固体表面）；当液滴后退时，固体表面已被一层厚吸附膜覆盖，此时的固-气界面张力为 γ_SV_f（覆盖膜的表面张力）。两者的差值即为铺展压 π_e = γ_SV – γ_SV_f。假设两种状态下固-液界面张力 γ_SL 相同（即液滴与固体的相互作用不受吸附膜影响，仅固-气界面改变），且两种状态下的相互作用参数 φ 相同，经过热力学推导得到：

γ_SV = γ_LV * (cos θ_R – cos θ_A) * (1 + cos θ_R)² / [(1 + cos θ_R)² – (1 + cos θ_A)²]  (15)

以及固-液界面张力：

γ_SL = γ_SV – γ_LV cos θ_A  (16)

物理意义：该方法天然地将铺展压 π_e 纳入了计算框架，无需事先假设 π_e = 0。它仅需要同一种液体在固体上的前进角和后退角，避开了不同液体之间参数不一致的麻烦。对于具有显著接触角滞后的体系（如粗糙表面、表面活性剂处理的聚合物、生物材料），滞后角法往往给出比其他方法更接近真实值的 γ_SV。

实例验证：Chibowski 等人对 PTFE、PE、PMMA、PET 等聚合物，以及石英、玻璃等无机材料进行了验证。对于 PTFE-水体系（θ_A ≈ 112°，θ_R ≈ 90°），由式 (15) 计算得 γ_SV ≈ 20.4 mN/m，与 Fowkes 法结果（20.2–20.4 mN/m）非常吻合。对于石英-水体系（θ_A ≈ 45°，θ_R ≈ 20°），计算得 γ_SV ≈ 42 mN/m，与两液法结果相近。

局限性：

• 该方法假设前进和后退过程中 φ 相同，但实际可能不同，尤其是当固体表面在液体接触过程中发生化学变化（如水合、离子交换）时。

• 对于低滞后体系（θ_A – θ_R < 5°），式 (15) 分子分母均趋于零，计算误差急剧放大，不适用。

• 要求 θ_A 和 θ_R 的测量精度高（通常需 ±1° 以内），且液滴体积变化速率需足够慢以接近准静态。

• 目前尚未建立适用于多液体组合的扩展形式（即无法同时利用多种液体的滞后信息提高精度）。

尽管存在上述局限，滞后角法在生物材料、粗糙表面、表面活性剂处理材料等难以满足理想固体条件的体系中具有独特优势，值得进一步发展和应用。

3.4 状态方程法在高能金属表面的根本局限性

在进入具体金属数据分析之前，需要重申：所有表面自由能计算方法的可靠性，首先取决于接触角测量本身的准确性。对于高能金属表面，常规实验室空气环境下测量的接触角反映的是污染层或氧化层的表面能，而非金属本征值。状态方程法假设 γ_LV > γ_SV 且 π_e = 0，这对于清洁金属（γ_SV >> γ_LV）不成立。

以下数据来自文献，但必须强调：这些接触角值是在常规实验室空气环境中测得的，固体表面实际覆盖有氧化物、水膜或有机物污染层。

从表中可以清楚地看到，用状态方程法从常规水接触角计算的 γ_SV 值（17–55 mN/m）与实际金属的清洁表面能（1100–2300 mN/m）相差一到两个数量级。这种巨大偏差并非状态方程法本身的数学问题，而是输入数据（接触角）从根本上不具备代表性。因此，强烈不建议使用状态方程法或任何基于空气环境接触角的常规模型来计算金属的真实表面自由能。对于金属，应采用超高真空下的液态金属接触角法或两液法（参见第 5 节），同时必须结合第 2.6 节的实验质量控制措施。

4. 第二类方法：基于表面张力组分的方法

4.1 Fowkes 法及其色散分量的测定

Fowkes（1962, 1964, 1972）是第一个系统提出将表面张力分解为不同分子间相互作用分量的人。他认为表面张力 γ 是各独立分量之和：

γ = γ^d + γ^p + γ^h + γ^i + γ^ + γ^π + γ^e + …  (17)

式中：

• γ^d：色散分量（London dispersion force），存在于所有物质中，来源于瞬时偶极-诱导偶极相互作用。

• γ^p：极性分量（Keesom orientation force），来源于偶极之间的相互作用，仅存在于极性分子中。

• γ^i：诱导分量（Debye induction force），来源于偶极在相邻分子上诱导的偶极。

• γ^h：氢键分量（hydrogen bonding），是一种特殊的定向强相互作用。

• γ^：酸碱分量（Lewis acid-base interaction），广义的电子对给体-受体相互作用，包括氢键。

• γ^π：π 电子分量（芳香体系中的 π-π 堆积）。

• γ^e：静电分量（离子固体表面电荷引起的静电作用）。

Fowkes 指出，尽管表面张力包含多种贡献，但在实际应用中，除色散分量外，其他所有非色散分量可以暂时合并为一个“极性"或“非色散"项（γ^a）。更重要的是，他通过实验和理论推导证明：当两相之间的相互作用仅由色散力主导时（即两相均为非极性，或至少有一相的极性贡献可忽略），色散部分的几何平均结合规则严格成立：

γ_^d = γ_1^d + γ_2^d – 2 √(γ_1^d γ_2^d)  (18)

该式的理论基础是 London 色散常数 C_6^ = √(C_6^ C_6^)（Berthelot 规则），以及分子间距离 d_ = √(d_ d_) 的假设。Fowkes 利用水-饱和烃界面张力实验数据，代入式 (18) 并令 γ_12 = γ_水-烃，γ_1 = γ_水，γ_2 = γ_烃，且设 γ_烃^d = γ_烃，求解得水的色散分量 γ_水^d = 21.8 ± 0.7 mN/m。这一数值至今仍被广泛使用，尽管有更精确的测定（如 Janczuk 等得到 26.85 mN/m）提出异议。

对于固-液体系，结合 Young 方程（忽略铺展压）可得：

γ_LV (1 + cos θ) = 2 √(γ_SV^d γ_LV^d)  (19)

若采用非极性液体（γ_LV = γ_LV^d，γ_LV^p = 0），则固体色散分量可由下式直接计算：

γ_SV^d = γ_LV (1 + cos θ)² / 4  (20)

标准应用实例：对于 PTFE，采用正构烷烃系列（己烷至十六烷），测得接触角随碳链增加而增大（烷烃表面张力从 18.4 到 27.1 mN/m），将每个烷烃的 γ_LV 和 θ 代入式 (20) 得到 γ_SV^d ≈ 20.2–20.4 mN/m，平均 20.3 mN/m，与 PTFE 的表面张力（18–20 mN/m）基本一致。对于聚乙烯（PE），用同样方法得 γ_SV^d ≈ 33.7 mN/m，接近于 PE 的表面张力（约 33–35 mN/m）。这验证了 Fowkes 法的自洽性。

局限性与误区：

• Fowkes 法假设非色散分量的贡献为零或可忽略，对于极性固体（如 PMMA、尼龙、纤维素）不适用，因为 γ_SV^p 贡献显著（可达 10–30 mN/m）。

• 许多研究者错误地将式 (20) 中的 γ_LV 直接用水的表面张力（72.8 mN/m）代入，这会产生荒谬的结果（因为水的 γ_LV^d 仅为 21.8，而非 72.8），必须使用非极性液体（如二碘I甲烷、α-溴萘、正构烷烃）。

• 对于表面张力极低的固体（如氟聚合物，γ_SV ≈ 10–15 mN/m），即使非极性液体也可能部分铺展（θ 接近于 0），需谨慎选择表面张力的烷烃（如己烷）以避免全润湿。

4.2 Owens-Wendt 几何平均法（两液法）

Owens 和 Wendt（1969）将 Fowkes 方法推广到同时包含色散和极性（非色散）相互作用的体系。他们假设固-液界面张力由色散分量和极性分量共同贡献，且两者均采用几何平均结合规则：

γ_SL = γ_S + γ_L – 2 √(γ_S^d γ_L^d) – 2 √(γ_S^p γ_L^p)  (21)

结合 Young 方程：

γ_LV (1 + cos θ) = 2 √(γ_S^d γ_L^d) + 2 √(γ_S^p γ_L^p)  (22)

这是一个关于两个未知数（√γ_S^d 和 √γ_S^p）的线性方程。因此需要两种探测液体：一种非极性（γ_L^p = 0），一种极性（γ_L^p > 0）。通常采用的标准液体对是水 + 二碘I甲烷。二碘I甲烷被当作非极性液体（γ_D = 50.8 mN/m，γ_D^d = 50.8，γ_D^p = 0）。水的参数为 γ_W = 72.8，γ_W^d = 21.8，γ_W^p = 51.0（这些值来自 Fowkes 的水-烃界面张力实验）。将两种液体的接触角 θ_W 和 θ_D 代入式 (22)，得到两个方程：

γ_W (1+cos θ_W) = 2 √(γ_S^d·21.8) + 2 √(γ_S^p·51.0)  (23)
γ_D (1+cos θ_D) = 2 √(γ_S^d·50.8)  (24)

由式 (24) 直接解出 γ_S^d，然后代入式 (23) 解出 γ_S^p。γ_SV = γ_S^d + γ_S^p。

替代液体对：除水和二碘I甲烷外，常用组合还包括：

• 水 + α-溴萘（γ_Br = 44.4，γ_Br^d = 44.4，γ_Br^p = 0）

• 甲酰胺（γ_F = 58.0，γ_F^d = 39.0，γ_F^p = 19.0）+ 二碘I甲烷

• 甘油（γ_G = 64.0，γ_G^d = 34.0，γ_G^p = 30.0）+ 二碘I甲烷

但需注意，不同来源的 γ_L^d 和 γ_L^p 值存在差异，这会导致计算得到的 γ_S^d 和 γ_S^p 波动。

对聚合物的典型结果（Owens & Wendt, 1969；Zdziennicka et al., 2017；Janczuk et al., 2026）：

• PTFE：γ_S^d ≈ 18.5，γ_S^p ≈ 0，γ_SV ≈ 18.5

• PE：γ_S^d ≈ 33.5，γ_S^p ≈ 0–1，γ_SV ≈ 34.0

• PMMA：γ_S^d ≈ 39–42，γ_S^p ≈ 8–12，γ_SV ≈ 48–54

• PET：γ_S^d ≈ 38–41，γ_S^p ≈ 4–8，γ_SV ≈ 44–48

• 尼龙 6：γ_S^d ≈ 35–40，γ_S^p ≈ 10–15，γ_SV ≈ 48–55

优缺点：

• 优点：实验简单，仅需两种液体，计算为初等代数；对大多数中等极性聚合物（PMMA、PET、尼龙、PVC）给出合理且可重复的结果；与 Zisman 临界表面张力相比，提供了极性和色散分量的分离，有助于理解表面化学。

• 缺点：将氢键、偶极、诱导、酸碱等多种相互作用笼统归为单一“极性"项，信息分辨率低；对强酸碱体系（如纤维素、壳聚糖、蛋白质）预测能力差，因为 γ^p 不能反映电子给体/受体的不对称性；当 γ_S^p 计算值为负时（可能由于接触角测量误差或液体参数不匹配），表明方法失效，应重新测量或换用其他液体对。

4.3 Wu 调和平均法

Wu（1971）基于 Hamaker 常数和 Lennard-Jones 势能理论提出了不同的结合规则。他认为，对于色散力，几何平均是合理的（因为 London 色散常数符合 Berthelot 规则）；但对于极性相互作用（偶极-偶极、偶极-诱导偶极），Heitler-London 理论显示更接近调和平均。因此，他提出：

对于低能表面（如聚合物），固-液界面张力由色散分量的调和平均和极性分量的调和平均共同组成：

γ_SL = γ_S + γ_L – 4 γ_S^d γ_L^d/(γ_S^d+γ_L^d) – 4 γ_S^p γ_L^p/(γ_S^p+γ_L^p)  (25)

结合 Young 方程：

γ_LV (1+cos θ) = 4 γ_S^d γ_L^d/(γ_S^d+γ_L^d) + 4 γ_S^p γ_L^p/(γ_S^p+γ_L^p)  (26)

对于高能表面（如汞、玻璃、金属、石墨），Wu 建议色散分量用几何平均，极性分量仍用调和平均（混合平均）：

γ_SL = γ_S + γ_L – 2 √(γ_S^d γ_L^d) – 4 γ_S^p γ_L^p/(γ_S^p+γ_L^p)  (27)

对应的 Young 方程为：

γ_LV (1+cos θ) = 2 √(γ_S^d γ_L^d) + 4 γ_S^p γ_L^p/(γ_S^p+γ_L^p)  (28)

Wu 调和平均法通常也需要两种液体（一非极性，一极性）。但求解过程中需要解二次方程组（因为 γ_S^d 和 γ_S^p 以线性与非线性混合形式出现），与 Owens-Wendt 法相比计算稍繁琐。

比较研究：Dalal（1987）对 12 种聚合物使用 6 种液体的接触角数据，比较了几何平均和调和平均两种方法。他发现在大多数情况下，两种方法计算出的 γ_SV 非常接近（差值 < 2 mN/m），但几何平均的数据拟合优度（R²）略高于调和平均。Janczuk 和 Bialopiotrowicz（1989）对 PTFE、PE、PET 和 PMMA 进行了系统研究，发现调和平均法对二碘I甲烷-水组合会给出 γ_S^d > γ_L（物理上不可能，因为色散分量不能超过总表面张力），因此他们认为几何平均更可靠。

当前共识：调和平均法的理论依据较弱（Heitler-London 理论主要针对气体分子，凝聚态中的极性相互作用更复杂），且实验验证并未显示调和平均优于几何平均。因此，Owens-Wendt 几何平均法在聚合物润湿领域占据主导地位。

4.4 van Oss-Chaudhury-Good 的 LW-AB 法

4.4.1 物理基础：将极性相互作用拆分为酸-碱不对称相互作用

van Oss、Chaudhury 和 Good（1986–1988）对 Fowkes 和 Owens-Wendt 的方法进行了最深刻的改进。他们认为，将表面张力简单划分为“色散"和“极性"两个分量过于粗糙，因为极性相互作用中，氢键和 Lewis 酸碱相互作用具有本质的不对称性——酸只能与碱反应，而不能与另一个酸发生强特异性作用。因此，他们提出将表面张力分解为：

γ = γ^ + γ^  (29)

其中：

• γ^：Lifshitz-van der Waals 分量，包括色散（London）、偶极-偶极（Keesom）和偶极-诱导偶极（Debye）三种电动力学相互作用。根据 Lifshitz 的宏观连续介质理论，这三种相互作用的衰减长度和温度依赖性相似，可以合并处理。在凝聚相中，Keesom 和 Debye 贡献通常小于表面张力的 2%，因此 γ^ 在数值上接近 Fowkes 的 γ^d，但物理含义更广。

• γ^：Lewis 酸碱分量，来源于电子对给体（碱）与电子对受体（酸）之间的短程相互作用（作用范围 < 5 Å，远小于 LW 作用范围）。

酸碱分量进一步分解为两个独立的参数：电子受体参数 γ^+（Lewis 酸度，空轨道接受电子对的能力）和电子给体参数 γ^-（Lewis 碱度，孤对电子给出的能力）。两者通过几何平均相关联：

γ^ = 2 √(γ^+ γ^-)  (30)

这一表达式基于 Small 和 Kollman 关于短程力结合规则的理论，以及 Drago 在溶液热力学中提出的 E 和 C 方程（但简化为单参数几何平均）。

不对称相互作用规则：两相之间的酸碱相互作用自由能变化为：

ΔG_^ = –2 ( √(γ_1^+ γ_2^-) + √(γ_1^- γ_2^+) )  (31)

注意到交叉项的出现：相 1 的酸（γ_1^+）只与相 2 的碱（γ_2^-）相互作用，相 1 的碱（γ_1^-）只与相 2 的酸（γ_2^+）相互作用。这是与 Owens-Wendt 极性项（对称的 2 √(γ_1^p γ_2^p)）的本质区别。

LW 部分的结合规则：假设 γ^ 分量符合 Fowkes 的几何平均（因为 LW 相互作用符合 Berthelot 规则）：

γ_^ = γ_1^ + γ_2^ – 2 √(γ_1^ γ_2^)  (32)

结合酸碱部分，总界面张力为：

γ_ = γ_1 + γ_2 – 2 √(γ_1^ γ_2^) – 2 ( √(γ_1^+ γ_2^-) + √(γ_1^- γ_2^+) )  (33)

结合 Young 方程，得：

γ_LV (1+cos θ) = 2 ( √(γ_S^ γ_L^) + √(γ_S^+ γ_L^-) + √(γ_S^- γ_L^+) )  (34)

该方程包含三个未知数：γ_S^、γ_S^+ 和 γ_S^-。因此，需要三种探测液体：一种非极性（γ_L^+ = γ_L^- = 0），两种极性且已知 γ_L^、γ_L^+、γ_L^- 的液体。

4.4.2 标准液体的参数标定

为了求解式 (34)，必须首先知道探测液体的 γ_L^、γ_L^+ 和 γ_L^-。van Oss 等人采用了以下标定策略：

(1) 非极性液体：二碘I甲烷（CH₂I₂）和 α-溴萘（C₁₀H₇Br）被认为仅具有 LW 相互作用（γ^=0，即 γ^+=γ^-=0）。它们的表面张力通过悬滴法精确测量：γ_D = 50.8 mN/m，γ_Br = 44.4 mN/m。因此 γ_D^=50.8，γ_D^+=γ_D^-=0；γ_Br^=44.4，γ_Br^+=γ_Br^-=0。

(2) 参比液体水：水的 γ_W = 72.8，γ_W^ 通过水-饱和烃界面张力实验确定为 21.8 mN/m（Fowkes 法）。因此 γ_W^=72.8–21.8=51.0。为了将 γ_W^ 分解为 γ_W^+ 和 γ_W^-，van Oss 等人最初假设水是“双极性"且 γ_W^+ = γ_W^-，因此 γ_W^+ = γ_W^- = √(γ_W^/2) = √(51.0/2) = 25.5 mN/m。这个假设基于水的自电离特性和对称的氢键网络，但受到争议（一些研究认为水的 γ_W^+/γ_W^- 比值在 0.3 到 3 之间）。

(3) 其他极性液体（甘油、甲酰胺、乙二醇、DMSO 等）：采用三种方法之一确定其参数 (a) 通过测定它们与二碘I甲烷或 α-溴萘的界面张力，利用式 (34)（已知固体为非极性，即 γ_S^+=γ_S^-=0）反推 γ_L^+ 和 γ_L^-；(b) 测定它们在水（已知 γ_W^+、γ_W^-）和已知非极性固体（如 PTFE 或 PMMA）上的接触角，联立方程求解；(c) 测定它们与另一种已标定极性液体的界面张力。

常用液体的标准参数（van Oss 传统值）如下：

Janczuk 等人（1993, 2017）通过 PTFE 接触角重新测定了水的 γ^=26.85 mN/m（高于传统的 21.8 mN/m），并据此重新计算了甘油、甲酰胺等液体的酸碱参数。这些修正值在物理上可能更准确，但改变后所有固体参数需重新计算，目前尚未成为通用标准。为便于与文献比较，建议继续使用 van Oss 传统值。

4.4.3 数学求解与病态条件

式 (34) 可写成矩阵形式 A·x = b，其中 x = [√γ_S^, √γ_S^+, √γ_S^-]ᵀ，A 为 3×3 矩阵，其行由 [√γ_L^, √γ_L^-, √γ_L^+] 组成（注意顺序与 x 匹配），b = [γ_L (1+cos θ)/2] 为列向量。

Hollander（1995）系统研究了该线性方程组的病态问题。他定义液体 i 的酸-碱比值 Q_i = γ_i^- / γ_i^+。对于两种极性液体，若它们的 Q 值非常接近（即 ΔQ = |Q₁ – Q₂| 很小），则矩阵 A 的两行近似线性相关，条件数 Cond(A) 极大，导致求解时 √γ_S^+ 和 √γ_S^- 对接触角的微小误差极其敏感，可能出现负的平方根（无物理意义）。

通过数值模拟，Hollander 得出以下建议：

• ΔQ ≥ 15：方程组条件良好（Cond(A) < 10），结果稳定可靠。

• 5 ≤ ΔQ < 15：条件数中等（Cond(A)=10–100），结果可能仍有较大误差，需谨慎。

• ΔQ < 5：病态严重（Cond(A) > 100），结果不可靠，应避免使用。

常用液体对的 ΔQ 值（基于 van Oss 传统参数）：

• 水-甲酰胺：17.37 – 1.0 = 16.37

• 水-甘油：14.64 – 1.0 = 13.64（可接受，但略低于 15）

• 水-乙二醇：24.48 – 1.0 = 23.48（理论上更优，但乙二醇表面张力较低）

• 甲酰胺-甘油：17.37 – 14.64 = 2.73（禁止使用）

条件数的具体数值（Volpe & Siboni, 1997）：

• 水-甲酰胺：Cond(A) ≈ 3.5，

• 水-甘油：Cond(A) ≈ 4.9–5.4，良好

• 水-乙二醇：Cond(A) ≈ 5.8–9.8，中等

• 甲酰胺-甘油：Cond(A) ≈ 106–151，极差（严禁使用）

因此，最佳液体组合是：非极性液体（二碘I甲烷）+ 水 + 甲酰胺。若甲酰胺不可得或与固体反应，水+甘油组合勉强可用，但需接受较大误差。

4.4.4 全局优化与最小二乘扩展

为了克服病态问题和单个液体组合带来的偶然性，Volpe 和 Siboni（1997）建议采用全局非线性优化方法。他们将所有固体-液体对的接触角数据（多个固体，多个液体）和所有液体的表面张力参数（γ_L, γ_L^, γ_L^+, γ_L^-）以及固体的未知参数（γ_S^, γ_S^+, γ_S^-）同时拟合，通过最小化如下目标函数：

U = Σ_^ (γ_ – (γ_^ + 2√(γ_^+ γ_^-)))²

• Σ_^ (γ_ 的约束项)²

• Σ_^ Σ_^ [ γ_ (1+cos θ_)/2 – ( √(γ_^ γ_^) + √(γ_^+ γ_^-) + √(γ_^- γ_^+) ) ]²

通过这种方式，可以利用所有数据的信息，减少对单一液体组合的依赖。该方法的难点在于需要大量高质量接触角数据，且初始值选择很重要（可能收敛到局部极小）。目前尚未成为常规实践，但代表了 LW-AB 法未来发展的方向。

4.5 探针液体分量准确性的重要性及验证

所有组分法都依赖探针液体已知的分量值。文献中同一液体（如水）的分量值存在显著差异。这直接导致不同实验室对同一种固体计算出的 γ_S^p 或 γ_S^- 相差可达 10–20 mN/m。

建议：

• 尽量采用 van Oss 传统值以保持与大多数文献可比；

• 在计算报告中必须明确给出所采用的液体参数来源；

• 对于关键体系，使用至少两组不同极性液体（如水+甲酰胺 vs 水+甘油）进行计算，若结果不一致（差值 > 10%），则需采用第 2.6.1 节的方法检查液体是否被污染。

5. 高能固体表面自由能的测定：两液法（Schultz 法）

高能固体（如云母、新鲜玻璃、金属、石墨、氧化物）的清洁表面张力通常很高（从几百到几千 mJ/m²），远大于水及大多数有机液体的表面张力。按照 Young 方程，液滴会在其上全铺展（θ=0），无法通过常规空气-液滴-固体系测量接触角。Schultz、Tsutsumi 和 Donnet（1976‑1977）针对这一难题，开发了两液法（Solid‑Liquid‑Liquid, SLL），也称为“液下接触角法"。

5.1 基本原理与方程式

将高能固体浸没在一种液体 L₂ 中（通常为正构烷烃，如己烷、辛烷、癸烷、十六烷），然后在固体表面（在 L₂ 环境下）放置另一种与 L₂ 不混溶的液体 L₁ 的液滴（通常为水）。此时三相系统为：固体 S – 液体 L₁（水） – 液体 L₂（烃）。Young 方程在该体系中的形式为：

γ_ = γ_ cos θ_ + γ_  (35)

其中 γ_ 为固体-烃界面张力，γ_ 为固体-水界面张力，γ_ 为水-烃界面张力，θ_ 为水在固体上（在烃环境中）的接触角。

根据 Fowkes 的几何平均假设，并将烃视为仅具有色散分量的非极性液体（γ_H = γ_H^d，且 I_^P = 0，即无非色散相互作用），有：

γ_ = γ_S + γ_H – 2 √(γ_S^d γ_H)  (36)
γ_ = γ_S + γ_W – 2 √(γ_S^d γ_W^d) – I_^P  (37)

其中下标 H 代表烃，W 代表水，I_^P 为水-固体之间的非色散相互作用能（包括氢键、酸碱、偶极等）。水-烃界面张力 γ_ 由液体性质决定，可以通过测量获得（或查表）。

将式 (36) 和 (37) 代入式 (35)，消去 γ_S（固体总表面张力），整理得到：

γ_W – γ_H + γ_ cos θ_ = 2 √(γ_S^d) ( √(γ_W^d) – √(γ_H) ) + I_^P  (38)

关键点：式 (38) 中左边各项均可由实验测定（γ_W、γ_H 已知，γ_ 查表或测量，θ_ 测量），右边有两个未知数：√γ_S^d 和 I_^P。如果测量水在一系列不同正构烷烃（不同 γ_H）中的接触角，以 [γ_W – γ_H + γ_ cos θ] 为纵坐标，以 [√γ_W^d – √γ_H] 为横坐标作图，应得到一条直线，斜率为 2 √γ_S^d，截距为 I_^P。

这就是 Schultz 两液法测定高能固体色散分量的核心公式。

5.2 实验细节与结果（以云母为例）

Schultz 等人使用新鲜解离的云母（muscovite mica）作为模型高能固体。解离在正构烷烃中进行以避免空气污染，然后直接在该烃中测量水的接触角。

测得的数据（部分）：

γ_W = 72.6 erg/cm²，γ_W^d = 21.6 erg/cm²。

将数据代入式 (38) 作图，得直线斜率 = 2 √γ_S^d = 10.95，故 γ_S^d = (10.95/2)² = (5.475)² ≈ 30.0 erg/cm²。截距 I_^P = 100 erg/cm²。

因此，云母的色散分量 γ_S^d = 30.0 erg/cm²，非色散（极性+酸碱）相互作用能为 100 erg/cm²。总表面自由能 γ_S = γ_S^d + γ_S^P，其中 γ_S^P 需通过极性液体进一步测定。

5.3 极性分量的测定（两液法扩展）

将多种极性液体（L₂ = 二氯甲烷、氯仿、硝基甲烷、硝基Y乙烷、醇类、苯等）代替烃，测量水在其中的接触角。由于 L₂ 具有非零的极性分量，式 (38) 需修改为包含 I_^P 项。已知 γ_S^d 和 I_^P 后，可计算出每种 L₂ 对应的 I_^P。

发现 I_^P 与 √γ_^P 成良好线性关系：

I_^P = 2 √(γ_S^P γ_^P)  (39)

斜率给出 γ_S^P = 90 erg/cm²。因此云母的总表面自由能 γ_S = 30 + 90 = 120 erg/cm²。

验证：Bailey 和 Price 通过云母剥离实验（剥离不匹配的云母片，避免库仑力贡献）测得 γ_S = 120 erg/cm²，与两液法结果一致。云母-水界面张力计算值 γ_SW = 42 erg/cm²，而剥离法测得包含库仑贡献的 γ_SW = 107 erg/cm²，扣除库仑部分（约 50 erg/cm²）后与 42 erg/cm² 接近。这验证了两液法的正确性。

5.4 两液法的适用条件与推广

成功应用的条件：

• L₁ 和 L₂ 必须互不相溶或溶解度极低（通常 < 0.1%），否则液滴体积会因溶解而减小，接触角随时间衰减。

• L₁ 必须能够置换固体表面上的 L₂ 吸附层，即 L₁‑L₂ 界面张力应足够高，且 L₁ 对固体的亲和力大于 L₂ 对固体的亲和力。Schultz 导出了一个热力学判据：I_^P > 2 √(γ_S^d) (√γ_S – √γ_H)(√γ_H – √γ_W^d)。对于云母-水-烃体系，该条件满足。

• 固体必须是化学稳定的，不与 L₁ 或 L₂ 发生反应、溶胀或溶解。

应用范围：

• 无机晶体：云母、石英、方解石、萤石、硫化物矿物等。

• 金属：金、银、铜、铝、钢（需注意氧化层）。

• 碳材料：石墨、金刚石、碳纤维。

• 玻璃、陶瓷。

• 高能聚合物（如聚酰亚胺）虽然 γ_SV 在 40‑50 mN/m，仍可用两液法提高精度。

局限性：

• 需要多种烃（至少 4 种）进行线性拟合，实验工作量较大。

• 对某些固体（如金属），烃的吸附膜可能很稳定，水不能置换，导致接触角测量值偏大，γ_S^d 被低估。

• 不适用于多孔材料（液体会渗入孔道）或可溶胀聚合物。

5.5 两液中探针液体纯度的控制

在两液法中，水与烃直接接触，烃中微量杂质（如表面活性剂）会显著降低水-烃界面张力 γ_WH，从而影响式 (38) 的截距和斜率。因此，必须用第 2.6.1 节的方法分别测量水和烃的表面张力，以及水-烃界面张力，确保与文献值偏差 < 1 mN/m。

6. 多孔、粉末、纤维材料的 Washburn 渗透法

对于常规测角法无法应用的样品——多孔固体（如滤纸、陶瓷、岩石）、粉末（如药物粉末、颜料、催化剂）、颗粒（如离子交换树脂、砂粒）以及纤维束（如棉、尼龙、碳纤维）——Washburn 渗透法是实际可行的选择。该法基于液体在毛细管或多孔介质中的自发渗透动力学。

6.1 Washburn 方程的推导

考虑一个半径为 r 的圆柱形毛细管，垂直插入液体中。液体由于表面张力的驱动而上升，上升速度受粘滞力平衡。忽略重力（对于短时间渗透）和惯性项，由 Poiseuille 方程和 Laplace 压力差可得：

dx/dt = (r γ_L cos θ) / (4 η x)  (40)

其中 x 为渗透深度，t 为时间，η 为液体粘度。积分得：

x² = (r γ_L cos θ) t / (2 η)  (41)

实际测量中，通常记录渗透质量 m 而非深度 x。假设固体为均匀多孔介质，具有 n 个平行毛细管，每个半径为 r，则：

m = n π r² ρ x  (42)

代入式 (41) 得：

m² = (n² π² r^5 ρ² γ_L cos θ) t / (2 η)  (43)

整理为：

cos θ = (m² / t) * (2 η) / (n² π² r^5 ρ² γ_L)  (44)

定义仪器常数 C = 1/(n² π² r^5)，该常数仅取决于固体样品的堆积结构和孔道几何。但由于 n 和 r 难以独立测定，通常通过已知 θ=0 的参考液体（全润湿）来标定 C。对于参比液体（cos θ=1），有：

C = (m_ref² / t) / (ρ_ref² γ_) * (2 η_ref)  (45)

注意：不同文献中常将 C 合并进一个复合常数 B = C/2，则式 (44) 变为：

cos θ = (m² / t) * (η / (ρ² γ_L)) * B  (46)

其中 B 需通过参考液体标定。

6.2 实验装置与步骤

现代商用仪器采用力传感器（或高精度天平）记录随时间增加的渗透质量。典型步骤：

1. 将粉末或纤维装入底部带筛网的样品管中，压实至恒定堆积密度（通常通过固定填充质量和敲击次数控制）。

2. 用全润湿的参比液体（如正己烷或正庚烷，对大多数材料 θ=0）测量 m²/t，标定出仪器-样品组合常数 B。

3. 用目标测量液体（如水、甲酰胺、二碘I甲烷等）重复实验，得到 m²/t，代入式 (46) 计算 cos θ，进而得到 θ。

4. 使用第 3 节或第 4 节中的表面自由能方法（将 θ 视为前进角或平均接触角）计算 γ_SV 及其分量。

注意事项：

• 样品的填充状态必须高度可重复，因为 B 对孔隙率极其敏感（r^5 项）。建议每次测量后重新填充，并做 3‑5 次重复。

• 仅当 θ < 90° 时液体能自发渗入，所以该法不适用于疏水材料（θ > 90°）。对于疏水粉末，需用低表面张力液体（如乙醇-水混合液）并外推到水。

• 液体渗透可能引起粉末团聚或溶胀（如高吸水性聚合物），导致孔道堵塞，m²/t 曲线偏离直线（后期斜率下降）。应只取初始线性段（渗透质量平方与时间成比例）的数据。

• 由于孔道弯曲度和非圆柱形结构，Washburn 方程给出的 θ 是表观平均值，不等于在理想平坦表面上的 Young 接触角。对于粗糙和多孔材料，该方法仍需谨慎解释。

6.3 Washburn 法中的稳态接触角问题

Washburn 法记录的是液体渗入孔道的动态过程，其 m²/t 曲线往往非线性。为获得稳态接触角，应使用微小液滴（体积 < 0.5 μL）沉积在粉末柱顶部，待液滴渗入后测量渗透时间。对于可溶胀材料，需结合 3D 形貌（测量渗入前后样品体积变化）和 UV 分析（检测渗出液中是否含固体成分）来验证。

6.4 应用实例与局限

成功应用：

• 药物粉末（如乳糖、微晶纤维素）的润湿性评估，用于预测片剂崩解和溶出行为。

• 矿物粉末（如石英、方解石、煤粉）在浮选剂中的润湿性，用于浮选工艺优化。

• 碳纤维和玻璃纤维的表面处理效果评估，用于复合材料界面设计。

• 滤纸和纺织品（棉、羊毛、聚酯）的吸水性和抗污性测试。

主要局限：

• 无法区分前进角与后退角（渗透过程是动态的，但 Washburn 方程给出的是单一 cos θ 值）。

• 对粉末的粒径分布和形状敏感，不同批次可能 B 值不同，需频繁标定。

• 当液体与固体之间发生化学吸附或离子交换时，渗透过程中表面性质可能改变，违反恒定 θ 的假设。

• 对于亲水但可溶胀的材料（如木材、纸张、某些水凝胶），渗透初期 θ 可能真实反映润湿性，但后期结构变化导致解释困难。

尽管如此，Washburn 渗透法仍是目前工业界和研究实验室处理粉末和多孔材料表面自由能的标准化方法。

7. 微生物细胞表面自由能的系统比较研究

Sharma 和 Rao（2002）的工作是表面自由能计算方法比较研究的一个标准化的方法。他们利用 van der Mei 等人累积发表的 147 株微生物分离株的接触角数据，对 Fowkes 法、状态方程法、几何平均法、LW-AB 法进行了全面比较，并以革兰氏染色分类和 XPS 表面化学分析作为外部验证。

7.1 接触角数据的统计特征

147 株微生物包括：140 株细菌（其中 111 株革兰氏阳性，29 株革兰氏阴性）和 7 株酵母（Candida species）。接触角测量采用薄膜法（细菌沉积在醋酸纤维素滤膜上，干燥至生理相关含水量），使用四种探测液体：水（W）、甲酰胺（F）、二碘I甲烷（D）、α-溴萘（Br）。

统计结果：

• θ_W：范围从约 0°（某些亲水乳酸菌）到 120°（高度疏水链球菌），平均值约 55°，分布广泛。

• θ_F：范围从 7° 到约 90°，平均值约 45°，分布较水窄。

• θ_D：范围从 36° 到约 72°，平均值 53°，标准偏差 9.3%，95% 集中在 36‑72°。

• θ_Br：范围从 20° 到约 60°，平均值 37°，标准偏差 11.2%，95% 集中在 20‑60°。

关键观察：极性液体（水和甲酰胺）的接触角分布宽，说明微生物细胞表面的极性/酸碱性质差异显著；而非极性液体（二碘I甲烷和 α-溴萘）的接触角分布非常集中，说明色散分量 γ^ 相对恒定，符合“生物表面 γ^ ≈ 40 mN/m ± 10%"的规律。

7.2 色散分量（γ^）的测定：Fowkes 法与状态方程法的比较

使用二碘I甲烷接触角，通过 Fowkes 法（式 20）计算 γ^（即 γ_S^d）。使用 α-溴萘接触角同样计算。两种液体计算出的 γ^ 值存在系统偏差：γ^(D) 平均 32.7 mN/m，γ^(Br) 平均 35.6 mN/m，但两者在微生物间呈线性相关（R=0.89）。这表明色散分量对非极性液体的选择有一定依赖性，但相对顺序保持一致。

使用状态方程法（式 14）以二碘I甲烷计算 γ_SV，得平均值 34.5 mN/m；以 α-溴萘计算得平均值 41.5 mN/m。两者相差较大，且与 Fowkes 法的 γ^ 不一致（状态方程值普遍高 5‑10 mN/m）。状态方程法以水计算 γ_SV 时，低 θ_W 的菌株（亲水）给出高 γ_SV（~70 mN/m），高 θ_W 的菌株（疏水）给出低 γ_SV（~20 mN/m），且当 θ_W > 66° 时出现 γ_SV < γ^，导致 γ^AB 为负。这是状态方程法在极性体系中的致命缺陷。

7.3 几何平均法（Owens-Wendt）的内部一致性检验

采用 α-溴萘测定 γ_S^d = 35.6 mN/m（固定），再用水或甲酰胺分别计算 γ_S^p。

• 用水：γ_S^p vs. θ_W 呈良好二次关系（R²=0.937），γ_S^p 始终 ≥ 0，γ_SV = γ_S^d+γ_S^p = 35–75 mN/m，随 θ_W 增大而递减。

• 用甲酰胺：γ_S^p vs. θ_F 的二次拟合 R² 仅 0.514，说明甲酰胺作为极性探测液体效果很差。这是因为甲酰胺的 γ_F^p（19 mN/m）比水的 γ_W^p（51 mN/m）小得多，对 γ_S^p 的灵敏度低，加上接触角测量误差，导致结果分散。

结论：几何平均法应使用水作为极性液体，二碘I甲烷或 α-溴萘作为非极性液体。甲酰胺不适合作为几何平均法的极性液体（但可作为 LW-AB 法的第二极性液体）。

7.4 几何平均法的最小二乘拟合（四种液体同时使用）

将水、甲酰胺、二碘I甲烷、α-溴萘四种液体的接触角同时代入式 (22) 进行最小二乘回归，求解 γ_S^d 和 γ_S^p。结果：

• 得到的 γ_S^d 与直接用二碘I甲烷的 Fowkes 法结果高度一致（R²=0.96，斜率接近 1），说明最小二乘法“偏好"二碘I甲烷的信息。

• 得到的 γ_S^p 与用水+α-溴萘配对结果一致，与用甲酰胺+α-溴萘配对结果无关。这进一步表明甲酰胺不是几何平均法的合适液体。

因此，若决定使用几何平均法，建议只用水+二碘I甲烷，不必强行使用甲酰胺。

7.5 LW-AB 法的内部一致性

使用二碘I甲烷或 α-溴萘作为非极性液体，水+甲酰胺作为两种极性液体，分别计算 γ_S^、γ_S^+、γ_S^-。

(1) γ_S^+（电子受体）：90% 以上的菌株 γ_S^+ < 5 mN/m，许多小于 0.1 mN/m。改变非极性液体（D vs. Br）后，γ_S^+ 值的相关系数 R=0.74，尚可但不WM。极小值区域（~0）两种液体结果一致。

(2) γ_S^-（电子给体）：这是 LW-AB 法中最稳健的参数。两种非极性液体给出的 γ_S^- 值几乎一致（R² ≈ 0.99，斜率 ≈ 1），且截距接近 0。γ_S^- 的分布呈双峰：第一峰在 0–25 mN/m（约占 38% 菌株），第二峰在 35–65 mN/m（约占 57% 菌株）。仅少数菌株落在 25–40 mN/m 区间。

(3) γ_S^：两种非极性液体给出的 γ_S^ 值相关性良好（R=0.89），但系统性偏差仍然存在（Br 值比 D 值平均高 3 mN/m）。这可能是由于二碘I甲烷和 α-溴萘的 γ_L^ 标定值不同所致。

(4) γ_S^ 和 γ_S^：γ_S^ 对非极性液体选择的敏感性中等（R=0.82），γ_S^ 的敏感性稍高（R=0.85）。

7.6 革兰氏阴阳性菌的区分能力

将 γ_S^- 与革兰氏染色结果相关联，发现：

• 革兰氏阳性菌（111 株）：γ_S^- 分布较分散，38.7% 在 0–25 区间，54.9% 在 35–65 区间，6.4% 在 25–35 区间。

• 革兰氏阴性菌（29 株）：γ_S^- 分布明显偏向高值，仅 24.1% 在 0–25 区间，69.0% 在 35–65 区间，6.9% 在 25–35 区间。

统计检验（t 检验）显示革兰氏阴性菌的平均 γ_S^- 显著高于革兰氏阳性菌（p < 0.01）。这一差异与 XPS 表面化学分析一致：革兰氏阴性菌的 N/C、O/C、P/C 比值较低，说明细胞壁中肽聚糖和蛋白质含量较少，而富含脂多糖（LPS）的外膜使表面呈强 Lewis 碱性。

实际意义：γ_S^- 可作为快速区分不同类型细菌的表面化学指纹，无需繁琐的基因测序或生化鉴定。这一发现对生物膜研究、口腔微生物学、病原菌检测和矿物生物浮选具有潜在应用价值。

7.7 微生物细胞测量中的特殊注意事项

微生物细胞表面柔软、含水且易变形，常规接触角测量更易受钉扎和液体污染影响。Sharma & Rao（2002）采用的薄膜法和干燥至“生理相关含水量"是有效措施。在此基础上，建议：

• 使用微小液滴（≤ 0.2 μL），避免液滴被细胞层吸收；

• 用力学铂金板法检查水/甲酰胺在细胞膜上停留 30 秒后的表面张力，判断溶出物影响；

• 用 3D 共聚焦显微镜表征细胞层厚度和均匀性，用 UV 荧光检测表面蛋白质或多糖的分布。

7.8 各方法的最终评分

基于 147 株菌的分析，Sharma 和 Rao（2002）给出以下评分（0–5 分）：

综合建议：

• 对于生物材料（细胞、蛋白质、多糖、生物组织），优先选用 LW-AB 法（水+甲酰胺+二碘I甲烷），并报告 γ_S^+ 和 γ_S^-。

• 若只能进行两种液体的测量，则用水+二碘I甲烷的几何平均法，并接受 γ_S^p 仅为粗略的极性指标。

• 避免使用状态方程法处理极性固体或极性液体数据，因其易产生负的 γ_S^。

• Fowkes 法仅适用于非极性固体（如 PTFE、PE、PP）。

8. 总结与方法选择指南

8.1 方法选择决策树

根据固体类型和实验条件，建议按以下决策树选择方法：

1. 样品为多孔、粉末、纤维 → 必须使用 Washburn 渗透法（第 6 节）。若样品可压制成光滑平整片（如某些药物粉末），则可尝试常规测角法。

2. 样品为高能固体（γ_SV > 45 mN/m，且为无孔致密体） → 优先选择两液法（第 5 节）。若两液法设备不具备，可尝试状态方程法（用二碘I甲烷），但结果可靠性较低。

3. 样品为低能聚合物（γ_SV = 15–45 mN/m） → 进入子决策：

• 若只需快速估算润湿性（工业质控）→ Zisman 法（用正构烷烃系列测定 γ_C）。

• 若已知固体为非极性（PTFE、PE、PP）→ Fowkes 法（用二碘I甲烷，式 20）即得 γ_SV。

• 若固体为中等极性（PMMA、PET、尼龙、PVC）且希望得到极性分量 → Owens-Wendt 几何平均法（水+二碘I甲烷）。

• 若希望深入表征酸碱性质（如生物材料、含官能团聚合物）→ LW-AB 法（水+甲酰胺+二碘I甲烷），并确保 ΔQ ≥ 15。

• 若接触角滞后显著（θ_A – θ_R > 15°）→ Chibowski 滞后角法可作为交叉验证。

4. 样品为生物材料（细胞、蛋白质、多糖） → 强烈推荐 LW-AB 法（水+甲酰胺+二碘I甲烷），额外报告 γ_S^- 并与其他生物学表征（如表面电荷、XPS）关联。

8.2 常见错误与避免方法

1. 错误：用纯水在 PTFE 上的接触角通过 Owens-Wendt 法计算 γ_SV，却不测量非极性液体。
   纠正：Owens-Wendt 法需要两种液体，单独用水只能给出一个方程，无法解出 γ_S^d 和 γ_S^p。需同时测二碘I甲烷接触角。

2. 错误：将 Fowkes 法的式 (20) 中 γ_LV 直接代入水的表面张力 72.8。
   纠正：水的 γ_LV^d 仅为 21.8，应使用非极性液体（如二碘I甲烷、α-溴萘、正构烷烃）的 γ_LV 代入。若只有水，则需先查得 γ_W^d=21.8，代入式 (19) 计算。

3. 错误：使用状态方程法处理 PMMA-水体系，并接受 γ_SV=38.5 mN/m，却不检验其他液体的结果。
   纠正：状态方程法对极性体系结果不可靠。应用几何平均法或 LW-AB 法交叉验证。

4. 错误：使用甲酰胺和甘油作为 LW-AB 法的两种极性液体（ΔQ=2.73 < 15）。
   纠正：应换为水+甲酰胺（ΔQ=16.37）或水+甘油（ΔQ=13.6，勉强），或水+乙二醇（ΔQ=23.5，但乙二醇表面张力较低，需确保 θ>0）。

5. 错误：在多孔粉末上测量接触角时，直接使用测角法滴加液滴。
   纠正：液体会瞬间渗入，无法形成稳定液滴。必须采用 Washburn 渗透法。

8.3 关键液体参数推荐值汇总

注：对于 LW-AB 法，水的参数（γ_W^+=γ_W^-=25.5）虽有争议，但仍是事实上的标准，便于不同研究之间比较。若采用 Janczuk 修正值（γ_W^=26.85），则需同时修正所有其他液体的参数，并重新计算固体参数，导致与文献数据不可比。因此，除非整个研究领域共同改变标准，建议继续使用 van Oss 传统值。

8.4 实验质量控制清单

8.5 各材料推荐方法总表

9. 展望

尽管表面自由能的计算方法已经发展了近两个世纪，但仍有许多开放性问题待解决。

9.1 标准液体参数的绝对标定
目前所有方法都依赖于参比液体（特别是水）的 γ^、γ^+、γ^- 值的绝对标定。van Oss 的传统标定（γ_W^+=γ_W^-=25.5）缺乏理论依据，仅是方便假设。未来应结合第一性原理计算（DFT、分子动力学模拟）和独立实验（如反气相色谱、表面力仪、石英晶体微天平）来建立更可靠的标度。

9.2 动态表面自由能
目前的静态方法忽略了接触角的时间依赖性（液滴形成后数秒至数分钟内的变化），而这对于表面可重构的体系（如聚合物、蛋白质层、细胞膜）至关重要。发展时间分辨接触角技术，结合表面自由能的时间演化分析，将提供更真实的界面信息。

9.3 粗糙度和异质性的 Wenzel-Cassie 修正
大多数方法假设固体表面光滑均匀。对于实际粗糙表面，需结合 Wenzel 或 Cassie-Baxter 模型修正接触角。但这会引入新的参数（粗糙度因子 r、面积分数 f），使得表面自由能求解更复杂。机器学习方法或许可以同时从大量接触角数据中反演出本征 Young 角、粗糙度和化学组成。

9.4 高温和特殊环境下的测量
表面自由能对温度敏感（dγ/dT 通常为 –0.05 到 –0.2 mN/(m·K)）。大多数文献数据在室温（20–25°C）下测量，对于高温熔体（聚合物熔体、金属熔体）、高温涂层、低温生物保存等应用，需要扩展测量方法至变温条件。关于高温熔体表面张力技术，上海梭伦的真实液滴®/RealDrop®/TrueDrop®高温熔体表面张力和接触角测量仪SL200HT已经实现，采用超高温注射泵测试悬滴或约束停滴表面张力，测试熔体表面张力及其三明治效应。

9.5 无接触光学方法
近年来发展了通过液滴振荡（悬滴或坐滴）的形变分析，无需接触角直接反推表面张力；或将接触角测量与椭偏光谱、和频光谱（SFG）结合，获取表面化学信息。这些方法有望避免接触角测量中的人为误差（基线选择、液滴轮廓拟合），但尚处于实验室阶段。

9.6 开放数据库与标准操作流程
建议建立标准化的表面自由能数据库，收集不同实验室使用标准方法测得的参考聚合物（PTFE、PE、PMMA、PET）的 γ_SV 及各分量，作为内部质量控制标准。同时制定详细的样品制备（清洗、干燥、表面老化）和接触角测量（液滴体积、环境湿度、温度、前进/后退程序）的标准操作流程，以降低实验室间的偏差。

10. 结论

本文基于五篇经典文献，系统整合了固体表面自由能计算的主要方法，包括 Zisman 临界表面张力法、Neumann 状态方程法、Chibowski 滞后角法、Fowkes 法、Owens-Wendt 几何平均法、Wu 调和平均法、van Oss 等 LW-AB 法、两液法（Schultz 法）以及 Washburn 渗透法。通过对微生物细胞表面（147 株菌株）的比较分析，揭示了各方法的内部一致性及其适用范围：LW-AB 法信息丰富、内部一致性最好；几何平均法（水+二碘I甲烷）可作为简便替代；状态方程法不推荐用于极性体系。

针对高能金属表面，本文指出状态方程法及其他常规空气环境接触角模型均不适用，因为常规条件下测得的接触角反映的是污染层或氧化层的表面能，而非金属本征值。对于金属，应采用超高真空下的液态金属接触角法或两液法。

在此基础上，本文提出了三项实验质量控制观点：(1) 探针液体滴加后表面张力的原位验证（力学铂金板法）；(2) 微小液滴稳态接触角的获取（非接触式进液系统，液滴体积 ≤ 0.5 μL）；(3) 表面结构与化学多样性的独立表征（3D 形貌 + 可视 UV 分析）。将这些措施与理论方法选择相结合，可以显著提高表面自由能计算的准确性和可重复性，尤其在金属、生物材料和多孔体系等挑战性材料中。