纤维素与壳聚糖金属氧化物杂化界面的分子建模表征

《Scientific Reports》：Molecular modeling characterization of cellulose and chitosan metal oxide hybrid interfaces

　　

生物聚合物-金属氧化物杂化物的界面电子结构决定了其稳定性，然而该结构仍未得到充分理解。本文报道了在B3LYP/SDD水平下，对分别以ZnO和V₂O₅单独功能化及其三元组合功能化的纤维素和壳聚糖所构成的10个模型体系进行的密度泛函理论（density functional theory, DFT）研究。根据前沿分子轨道（frontier molecular orbital, FMO）分析，杂化使三元复合材料中的能隙从纤维素的7.106 eV和壳聚糖的6.908 eV分别降低至3.559 eV和3.496 eV，而纤维素-ZnO-V₂O₅的偶极矩从0.536 D增加至26.697 D。反应性描述符的系统性变化包括化学软度（0.141至0.281 eV^-1）和亲电性（0.890至5.418 eV）的增加。利用分子中原子量子理论（Quantum Theory of Atoms in Molecules, QTAIM）进行的拓扑研究表明，从原始聚合物（73–76个键临界点[BCPs]；149–155个总临界点）到具有多达95个BCPs和193个总临界点的复杂杂化物，呈现递进趋势。三元复合物中新的界面金属-氧（M-O）BCPs达到11个（3个Zn-O和8个V-O），证实了异质结的形成。在纤维素-ZnO-V₂O₅中，环临界点从9个增加到22个，同时出现新的笼临界点，表明有机-无机界面处存在环状模式。非共价相互作用分析解析了围绕这些M-O锚点的次级氢键和范德华接触。这种双重共价-非共价锚定机制将M-O界面确定为生物医学、传感和防护涂层应用中生物共轭和表面功能化的关键位点，并为杂化结构稳定性提供了量子化学解释。

**研究背景与科学问题**
纤维素是地球上最丰富的天然结构多糖，是植物细胞壁的主要成分，具有生物相容性、生物可降解性、强力学特性和可再生性等优点；壳聚糖是由D-葡萄糖胺和N-乙酰D-葡萄糖胺单元组成的阳离子多糖，其自由氨基在pH低于6.5时质子化（NH₃⁺），赋予其pH响应性溶解性和抗菌活性。ZnO是一种宽带隙半导体（约3.37 eV），在紫外区具有强吸收；V₂O₅因优异的电子、磁性和催化性能而具有重要技术应用。将多糖与ZnO和V₂O₅等金属氧化物复合可增强其物理化学性能，使其在生物医学、传感和防护涂层等领域具有广阔应用前景。然而，以往的计算研究仅涉及单氧化物杂化物，且多限于前沿轨道描述符分析，缺乏对有机-无机界面的拓扑表征。为此，研究人员开展本研究，旨在系统比较纤维素和壳聚糖与ZnO和V₂O₅形成的二元及三元杂化物的电子结构和界面相互作用，阐明杂化结构稳定性与界面反应性的量子化学机制。该研究成果发表于《Scientific Reports》。

**主要技术方法**
研究人员使用Gaussian 09软件（G09, Revision C.01；美国Gaussian Inc.，于埃及国家研究中心先进科学卓越中心分子建模与光谱实验室运行），在B3LYP杂化泛函结合SDD有效核势基组水平下，对10个模型体系进行密度泛函理论（DFT）计算。几何优化与频率计算同步进行以确保势能极小值。基于前线分子轨道能量计算电离势（IP）、电子亲和势（A）、化学势（μ）、化学硬度（η）、绝对软度（S）和亲电性指数（ω）等全局反应性描述符。采用Multiwfn程序（3.8版）进行QTAIM和NCI拓扑分析，并使用VMD软件（1.9.4版）进行可视化。模型体系包括裸纤维素、裸壳聚糖、裸ZnO、裸V₂O₅及其二元和三元复合物。

**研究结果**
*前沿分子轨道（FMO）分析*
裸纤维素和裸壳聚糖的HOMO-LUMO能隙最宽（分别为7.106和6.908 eV）。ZnO的能隙最窄（2.402 eV），而V₂O₅的能隙（6.405 eV）接近宽隙生物聚合物，其HOMO（-12.478 eV）和LUMO（-6.073 eV）均远低于其他所有结构。杂化后能隙一致收窄，三元复合物纤维素-ZnO-V₂O₅和壳聚糖-ZnO-V₂O₅的能隙分别降至3.559和3.496 eV，降幅约50%。偶极矩从裸纤维素的0.536 D大幅增至纤维素-ZnO-V₂O₅的26.697 D，但壳聚糖-ZnO-V₂O₅的偶极矩（5.965 D）仅略高于裸壳聚糖（5.878 D）。分子静电势（MESP）图显示金属中心（Zn或V）为最正电性位点，配位氧和羟基氧周围为负电性区域，表明金属-氧界面是亲核攻击和氢键供体的优先位点。

*非共价相互作用（NCI）分析*
约化密度梯度（RDG）分析显示，裸生物聚合物的等值面主要由糖苷键间和赤道羟基网络的绿色弱相互作用区域主导。复合物形成后，在所有杂化物的多糖-金属氧化物界面均出现新的明确等值面，表现为低RDG的蓝色强吸引特征，与QTAIM识别的新界面键临界点一致，表明多糖并非简单物理吸附，而是通过共价和非共价双重机制锚定于金属氧化物相，这对杂化材料在水性或生理介质中的胶体和界面稳定性具有直接意义。

*QTAIM拓扑分析*
复合物形成后总临界点（CPs）数量从裸纤维素的149个和裸壳聚糖的155个增加至二元和三元复合物的181–193个。裸纤维素和裸壳聚糖分别显示73和76个BCPs，平均拉普拉斯量为负值（-0.673和-0.669 a.u.），表明以共价键为主的网络。纤维素-ZnO-V₂O₅的BCPs从73增至94，环临界点（RCPs）从9增至22。三元复合物中界面M-O BCPs达到11个（3个Zn-O和8个V-O），证实了异质结的形成。V-O键的拉普拉斯量（+1.075 a.u.）和电子密度（0.212 a.u.）显著大于Zn-O键（+0.690 a.u.，0.137 a.u.），与钒的较高正式氧化态和键级一致。三元复合物均呈现最高的总CPs数（各193个），反映最大的结构和电子复杂性。

*反应性描述符与结构-性能关联*
杂化使化学软度从0.141 eV^-1（纤维素）和0.145 eV^-1（壳聚糖）增加至0.281和0.286 eV^-1，同时系统性地降低电离势和化学硬度。亲电性指数在裸氧化物中最高（ZnO为13.831 eV，V₂O₅为13.433 eV），在杂化物中为3.473–5.949 eV。能隙收缩（约50%）使光学吸收阈值移向可见光区并将LUMO能量（-2.822至-2.604 eV）置于有利于光激发电子注入的区间；高偶极矩（最高26.697 D）有利于极性分析物的吸附和电荷转移检测；低硬度（1.780–1.748 eV）和高软度则有利于电极/电解质界面的电荷转移动力学。

**讨论与结论**
综合FMO、MESP、NCI、QTAIM和全局反应性描述符分析，研究人员得出了一致的量子化学图像：金属氧化物片段主导了每个复合物的前沿轨道和静电行为；O-连接在拓扑上产生了一组新的键和环临界点，通过强极化的M-O配位键及周围较弱的氢键/范德华接触将两相结合。该双重共价-非共价锚定机制在三元复合物中最为显著，为功能化纤维素/壳聚糖-金属氧化物平台的结构完整性和界面反应性提供了量子化学解释，并将金属-氧（M-O）界面确定为表面功能化和生物共轭的原子级靶点。研究人员同时指出研究局限性：有限单体模型的计算能隙未收敛到固态（周期性）数值，未包含显式溶剂或生物介质，B3LYP/SDD水平未显式考虑色散作用（但ωB97X-D和B3LYP-D3(BJ)基准计算维持了所述趋势），描述符与性能之间的关联为相关性而非因果性。未来研究将纳入显式水和生理环境，扩展至周期性DFT和分子动力学模拟，并在生物传感和抗菌应用中对预期锚定行为进行实验验证。