功能化壳聚糖（–SH/–COO?）选择性吸附Hg(II)、Cd(II)和Pd(II)的密度泛函理论研究：基于ESP/MEP与NCI/RDG分析的计算机理与选择性调控

《Sustainable Chemistry》：Selective Biosorption of Hg(II), Cd(II), and Pd(II) on Functionalized Chitosan (–SH/–COO?): A DFT Study with ESP/MEP and NCI/RGD Analyses Joaquín Hernández-Fernández, Rafael González-Cuello and Rodrigo Ortega-Toro

为解决重金属污染环境中Hg(II)、Cd(II)、Pd(II)选择性去除的难题，研究人员采用密度泛函理论（DFT）探究了功能化壳聚糖（–SH/–COO?等）对三种金属的吸附机理。结果表明，羧酸微环境优先结合Cd(II)，巯基修饰最利于Hg(II)吸附，而天然壳聚糖对Pd(II)的稳定性最强，且吸附主要通过外球配位与非共价协同作用实现。该工作为设计高选择性生物吸附材料提供了电子结构层面的理论指导。

　　

随着工业活动的不断扩张，重金属污染已成为全球性环境挑战，尤其是汞（Hg）、镉（Cd）和钯（Pd）等离子的排放，因其在环境中的持久性、生物累积性和高毒性，对生态系统和人类健康构成严重威胁。这些重金属即使浓度很低，也可能引发神经毒性、肾脏损伤等慢性疾病。在水处理领域，开发高效、低成本且具备高选择性的吸附材料是当前的研究瓶颈。传统吸附材料往往在多种金属离子共存的复杂水体中表现不佳，其去除效率不仅取决于材料的吸附容量，更关键的是活性位点的精细化学特性以及控制金属与聚合物形态的酸碱微环境。因此，从分子层面理解吸附剂与特定金属离子之间的相互作用机制，对于设计下一代高性能吸附材料至关重要。

在此背景下，壳聚糖作为一种来源丰富、可生物降解的天然高分子，因其骨架上的氨基和羟基可作为金属配位点，而被视为极具潜力的生物吸附平台。然而，其性能受诸多因素影响，包括脱乙酰度、构象灵活性以及溶液pH值等，导致吸附选择性难以预测和控制。通过化学修饰，如在壳聚糖上引入硫基（–SH）或羧酸/羧酸根（–COOH/–COO^?）等功能基团，有望“编程”其配位偏好，实现对特定金属的靶向捕获。例如，软酸Hg(II)倾向于与软碱硫基结合，而Cd(II)和Pd(II)则表现出介于共价与静电作用之间的边界行为。尽管实验上已观察到这些趋势，但其背后的微观机理，特别是决定选择性的热力学和电子因素，仍缺乏系统深入的理论阐释。

为了回答上述问题，Joaquín Hernández-Fernández、Rafael González-Cuello 和 Rodrigo Ortega-Toro 在《Sustainable Chemistry》上发表了一项研究，他们采用密度泛函理论（Density Functional Theory, DFT），对天然壳聚糖（CTS）及其多种功能化衍生物（包括CTS–COOH、CTS–COO^?、CTS–NH₃⁺和CTS–SH）选择性吸附Hg(II)、Cd(II)和Pd(II)的过程进行了系统的比较研究。该工作旨在从计算化学角度，揭示不同功能微环境（羧酸、羧酸根、质子化氨基、硫基）对三种金属离子的亲和力差异及其根源，为理性设计功能化壳聚糖吸附剂提供理论框架。

研究团队主要运用了以下几项关键技术方法：1. 计算模型构建：采用壳聚糖四聚体模型，并引入功能化基团，以模拟局部化学环境。2. 密度泛函理论（DFT）计算：使用Gaussian 16软件，主要采用ωB97X-D泛函结合LANL2DZ基组，进行几何结构优化和振动频率分析，确保所得结构为势能面上的真实极小点。3. 显式微水化方案：将金属离子处理为统一的五水合阳离子[M(H₂O)₅]²⁺，以此作为比较不同金属和材料亲和力的参考态，并在此框架下计算结合自由能（ΔG_assoc）。4. 电子密度后处理分析：利用Multiwfn软件对波函数进行分析，生成静电势/分子静电势（Electrostatic Potential/Molecular Electrostatic Potential, ESP/MEP） 图以可视化静电吸引区域，并进行非共价相互作用（Non-Covalent Interaction, NCI） 及约化密度梯度（Reduced Density Gradient, RDG） 分析，以识别和量化稳定复合物的弱相互作用（如氢键、色散作用）。

3.1. 几何优化与理论方法选择（M06-2X-D3 vs. ωB97X-D；LANL2DZ）

研究首先比较了两种泛函（M06-2X-D3和ωB97X-D）对吸附剂模型和五水合金属复合物的计算结果。通过分析总能量、前沿轨道（HOMO/LUMO）及能隙，发现两种方法在定性地描述化学功能化引起的电子结构扰动趋势上保持一致。其中，ωB97X-D泛函能提供更平衡的长程和色散作用描述，因此被选为主力理论方法用于后续金属-生物聚合物复合物的构建与热力学讨论。

3.2. 吸附剂与微水化复合物的电离势与电子亲和能

通过基于Koopmans近似的电离势（I）和电子亲和能（A）分析，研究进一步细化了各体系的电子特征。结果显示，羧酸根模型（CTS–COO^?）在所有吸附剂中具有最高的电子接受倾向，而巯基和质子化氨基模型的电子亲和能则为负值。对于金属水合离子，两种泛函均显示电离势顺序为Cd < Pd < Hg，表明汞复合物最不易电离。

3.3. 最终优化结构及主导配位模式（N/O/S）

几何优化结果表明，在大多数研究的复合物中，金属离子的第一水合壳层（五个水分子）都得以完整保留，未发生水分子被壳聚糖杂原子（N, O, S）系统性的内界配体取代。例如，在CTS–M、CTS–COOH–M、CTS–NH₃⁺–M和CTS–SH–M系列中，金属与生物聚合物的最短接触距离均较长，表明稳定化主要源于外球结合，即通过静电吸引、水分子介导的耦合以及第二球层的重组来实现。

一个关键的例外出现在CTS–COO^?–Cd复合物中，其中一个羧酸根氧原子与Cd(II)的距离缩短至2.658 ?，显示出内界配位的萌芽特征，但金属的第一水合壳层仍然保留，表明这是一种部分内球或溶剂共享的排列方式。相比之下，Hg(II)和Pd(II)在CTS–COO^?上仍保持外球结合模式。

3.4. 生物吸附热力学：亲和力与选择性（ΔE_assoc/ΔG_assoc）

结合自由能（ΔG_assoc）计算揭示了清晰的金属依赖性亲和规律：羧酸（–COOH）微环境对Cd(II)的稳定化作用最强；巯基化（–SH）微环境为Hg(II)提供了最有利的结合；而天然壳聚糖（CTS）对Pd(II)的计算稳定化作用最强。从选择性角度看，天然CTS对Pd(II)表现出强选择性；CTS–COOH是“Cd特异性”的吸附剂；CTS–SH对Hg(II)和Pd(II)共选择性，并能有效排除Cd(II)；而CTS–COO^?和CTS–NH₃⁺则表现出对多种金属的中等亲和力。

3.5. 选择性的静电起源：ESP/MEP分析

ESP/MEP图谱显示，不同功能化显著改变了吸附剂表面的静电景观。CTS–COO^?在羧酸根周围产生了连续、强烈的负电势区域，是强大的静电吸引子。CTS–COOH则表现出更具方向性的偶极特征。CTS–SH显示出高各向异性和可极化的硫微环境特征。CTS–NH₃⁺则因质子化而呈现电正性区域。在形成复合物后，金属阳离子附近出现强电正性区域，而锚定位点的杂原子处负电势增强，这种静电指纹的分离是复合物形成的标志。研究指出，选择性不仅取决于静电势的强弱，更与静电作用的方向性、位点的可极化性以及第二球层的预组织有关。

3.6. 酸碱模拟：质子化状态与电荷一致性

通过对比不同质子化状态（如–COOH与–COO^?，–NH₂与–NH₃⁺）的模型，研究证实质子化状态通过改变给体位点可用性、静电方向性和局部水合组织来调节金属亲和力，而不仅仅是改变形式电荷。这解释了为何质子化的羧酸（–COOH）模型对Cd(II)的结合自由能优于去质子化的羧酸根（–COO^?）模型。

3.7. 次级相互作用与微水化的作用：NCI/RDG

NCI/RDG分析表明，所有复合物的稳定化都强烈依赖于色散作用和微环境堆积（对应于sign(λ₂)ρ ≈ 0的贡献），同时存在来自适度吸引相互作用（sign(λ₂)ρ < 0）的贡献。这说明全局的复合物稳定化源于多种弱相互作用的协同平衡，包括氢键、色散接触和局部重组，而非单一的强配位键。这一发现与几何优化结果（外球结合主导）一致，强调了第二球层组织在选择性吸附中的关键作用。

3.8. 与实验文献的对比

计算结果与部分实验观察定性一致。例如，实验证实硫基化壳聚糖对Hg(II)的吸附增强，羧甲基化壳聚糖能改善对Cd(II)等二价金属离子的吸附性能。然而，研究也指出，由于计算采用了简化的微水化阳离子模型，且真实水溶液中的金属形态更为复杂，因此计算结果应视为可在受控条件下进行实验验证的机制假说。

结论与讨论

本研究在统一的显式微水化框架下，通过DFT计算系统比较了功能化壳聚糖对Hg(II)、Cd(II)和Pd(II)的选择性吸附。主要结论是，功能化以金属特异性和非单调的方式重塑了壳聚糖的吸附特性：羧酸（–COOH）微环境最有利于Cd(II)，巯基化（–SH）微环境对Hg(II)的结合最有利，而天然壳聚糖对Pd(II)的计算稳定化作用最强。结构分析表明，大多数情况下金属保留其第一水合壳层，主导的稳定机制是外球结合，并辅以诱导极化和协同的第二球层相互作用，而非系统的内界配体取代。ESP/MEP分析揭示了决定选择性的异质且具方向性的静电吸引区域，而NCI分析则证实稳定化源于色散和适度吸引性弱相互作用的累积效应。

这项研究的重要意义在于，它从电子结构层面提供了一个可比较的理论框架，将功能化壳聚糖的化学修饰（引入–SH/–COO^?等基团）与对特定重金属离子的选择性吸附能力直接关联起来。研究指出，设计高选择性吸附剂时，不应仅仅关注引入基团的电荷或硬度/软度，还需综合考虑其创造的局部微环境的静电方向性、可极化性以及组织第二球层弱相互作用网络的能力。这些计算获得的机理见解和选择性趋势，为未来在明确pH、离子强度和配体条件下，定向设计和实验验证高性能壳聚糖基吸附材料提供了重要的理论指导和可检验的假设。尽管计算的绝对能量值可能受到水化处理方式和模型简化程度的影响，但研究所揭示的相对亲和力层次和化学上截然不同的微环境效应，是指导后续实验研究和材料开发的坚实理论基础。