土壤有机质组成、含水量与竞争性Cl?对高岭石界面电双层及微团聚的分子动力学调控机制

《Minerals》：Electric Double Layers in Soils as Basis for Micro-Aggregation

　　

固-液相结合处的电双层（electric double layer, EDL）界面区域内离子行为在表面科学与材料科学中已较为清楚。然而，在土壤科学中并非完全如此；由于吸附于矿物表面的带负电土壤有机质（soil organic matter, SOM）分子本身复杂多样，EDL在该体系中呈现显著复杂性。此方面知识不足导致若干问题：在分子尺度上，何种相互作用机制促成所谓有机-矿物团聚体（organo-mineral aggregates, OMAs）这类土壤结构单元的稳定。研究人员采用经典分子动力学（classical molecular dynamics, MD）模拟，研究高岭石（Kaolinite）与四种标准化土壤有机质样品之间的相互作用；这些样品具有不同芳香族和羧基含量，并由Ca²⁺阳离子使有机质（organic matter, OM）分子聚集并将其稳定于表面。研究人员评估了（i）有机质含量（organic content）、（ii）含水量（water content）和（iii）竞争性Cl^?离子的影响。在低含水量且脂肪族组分较高的土壤中，团簇聚集体更有利；而在高含水量且存在Cl^?离子的条件下，高芳香族组分土壤中的聚集更有利。此外，低酸组分样品因非键相互作用和氢键而更有利于团簇形成。上述环境均体现为各OM样品独特的电荷密度“指纹”（charge-density “fingerprints”）和表面覆盖度（surface coverage）。该工作朝着表征和研究可针对特定局部条件调控的天然土壤迈出一步。

研究背景与问题：离子在固-液界面电双层（electric double layer, EDL；界面电荷与离子分布层）中的行为在表面科学和材料科学中已有较多认识，但在土壤科学中仍不充分。土壤有机质（soil organic matter, SOM）分子带负电、组成复杂且多样，吸附于矿物表面后形成所谓有机-矿物团聚体（organo-mineral aggregates, OMA），其分子尺度稳定机制尚不清楚。土壤由矿物、有机化合物、生物及液气相组成，矿物为微团聚提供基础并影响分散、絮凝、吸附和物质迁移等界面性质；OM是重要碳氮储存库，OMA对改善土壤性质、生物修复和减缓气候变化具有重要意义。已有工作多解析单一组分吸附或构建简化OMA模型，部分真实模型可描述结构但缺乏EDL热力学信息。为此，研究人员以高岭石（Kaolinite）为代表性黏土矿物，系统研究不同有机组成、含水量及竞争性Cl^?条件下OMA的聚集与稳定机制。

方法概况：研究人员从国际腐殖物质学会（IHSS）选取四种标准化样品：Suwannee River III fulvic acid（sfa；3S101F）、Suwannee River III humic acid（sha；3S101H）、Pony Lake fulvic acid（pnl；1R109F）和Leonardite humic acid（leo；1S104H），以VSOMM在中性pH生成OM模型；构建16×16×1高岭石空腔（体积8.2×14.3×0.7 nm³，z向增高10 nm），设置半充水、全充水及全充水加0.1 N CaCl₂三类环境，采用GROMOS 54a8、ClayFF和SPC模型，以GROMACS 2023进行经典MD，生产400 ns并三次重复；以28×48网格计算占据概率θ_ij，并用Gini系数（Gini coefficient, G；空间异质性指标）量化分布，另设benzoate（bzo）对照。

研究结果：电荷密度分布（Charge-Density Profile）：通过Ca²⁺和OM电荷密度剖面，发现近铝醇表面出现z=0.19 nm的内球配合物（inner-sphere complexes, ISC）峰和z=0.42 nm的外球配合物（outer-sphere complexes, OSC）峰；高PC2样品sfa、sha在ISC区峰最高，leo在全充水时无ISC；部分充水强化Ca²⁺桥接，全充水使Ca²⁺迁移至硅氧烷侧；Cl^?吸附于铝醇侧并降低ISC区OM峰，提示其将OM从表面挤出。

有序参数（Order Parameter）：通过水偶极矩有序参数，发现近矿物表面水强极化，z>1 nm后趋于零并呈体相随机取向；部分充水在z=0.82 nm后较全充水更随机；样品间无差异，加CaCl₂不明显影响溶剂有序。

静电势（Electrostatic Potential）：通过静电势剖面，发现OM与Ca²⁺的离子贡献近似Poisson–Boltzmann势；水贡献影响z<0.2 nm处Ca²⁺势斜率，z>1 nm处水势与离子势相互抵消使总势为零，支持远离界面呈体相行为。

跨EDL的复杂结构（Complex Structures Across the EDL）：通过统计Ca²⁺形成ISC+OSC的总数，发现其与主成分PC1（主要反映芳香度）的相关性高于PC2（主要反映酸度/羧基组成）；高芳香性leo有利于形成阳离子桥；部分充水Ca²⁺位点最多，全充水不相关；Cl^?改变趋势并使低芳香性pnl有利于桥接OM。

非键相互作用（Non-Bonding Interactions）：通过非键成对相互作用能，发现OM-Ca²⁺与PC2强相关，OM-OM与PC1和PC2强相关，且以静电贡献为主；低PC2样品有利于阳离子桥，低芳香含量有利于OM聚集；pnl高脂肪族通过疏水基团包埋、羧基暴露于水，Ca²⁺经羧基间阳离子桥连接聚集体。

氢键（Hydrogen Bonds）：通过氢键数量分析，发现OM–矿物氢键与PC2较强相关，OM-OM氢键与PC1和PC2相关；PC2越高，全充水时与表面氢键越多、吸附能力增强；部分充水因阳离子桥主导而氢键形成减少；高PC1的OM-OM氢键增加，说明高脂肪族通过氢键聚集程度较低；Cl^?增加氢键数并有利于OM聚集。

表面占据（Surface Occupancy）：通过占据概率热图和Gini系数，发现sfa在部分充水中近表面OM固定、少量可动；bzo部分充水虽经Ca²⁺桥接但不形成团簇且分布均匀，全充水无特异吸附位点；低水时随芳香性增加，平均占据概率θ增加、G降低，低芳香样品低覆盖且局部位点团簇，高芳香样品覆盖更广且较均匀；高水时随芳香性增加，θ降低、G增加，趋势反转；Cl^?降低表面覆盖并促进团簇。

讨论部分总结：OM的多阴离子性质使其在铝醇表面强吸附，Ca²⁺主要经OSC桥接OM与矿物，但ISC区域峰强反映少量高稳定Ca²⁺的桥接作用。不同SOM类型导致吸引与排斥相互作用不同，形成独特EDL指纹。高芳香样品在低水时有利于阳离子桥和最大表面分布；高PC2样品虽阳离子桥较弱，但可由较多氢键补偿而形成中等团簇。高水促进Ca²⁺向硅氧烷侧迁移，驱动OM阴离子部分向铝醇侧并增加OM–矿物氢键，同时因铝醇侧Ca²⁺减少而促进OM分散；高水体系平均θ较大。Cl^?替换部分吸附OM、参与EDL屏蔽但不明显扰动水偶极，并反转芳香性依赖、促进团簇；移除Cl^?后OM覆盖可逆恢复约98%。

研究结论翻译：OM吸附于矿物表面主要由静电现象驱动。不同有机组成的标准化土壤样品在固-液界面表现出独特静电“指纹”，并因EDL形成导致吸附物种的团簇与分布。矿物表面占据强烈受SOM芳香含量影响；虽未直接评估，结果提示酸性组成也发挥作用。Ca²⁺等阳离子在聚集和吸附机制中必不可少，形成阳离子桥，使OM分子聚集并固定于高岭石表面；氢键通过暴露于环境的带电物种维持OM之间及OM与矿物表面的连接。含水量通过极化离子物种影响空腔内OM吸附机制：低水有利于大OM团簇及经阳离子桥吸附；高水促进空腔内极化，驱动阳离子至硅氧烷侧，促进超分子结构分解并使OM分布于矿物表面，最终形成矿物界面质子与OM分子间氢键。相反，高水空腔中可通过增加溶液离子强度促进OM团簇，重组聚集体并以单价阴离子占据矿物表面；或通过增加SOM脂肪含量，将脂肪族基团包埋于OM结构内部并暴露可溶基团以稳定OM聚集体。据研究人员所知，这是首次以原子分辨率描述精炼土壤模型EDL界面，并通过改变溶剂环境和有机输入组成进行研究。该行为提示特定OM-OM和OM–矿物相互作用在OMA形成与稳定过程中随环境变化。这种调控水平代表土壤建模的重要改进，可为特定土壤的生物修复、碳封存或全球变化缓解靶向方法提供科学基础。