编辑推荐：
本文通过多尺度计算模拟（DFT/MD/MC），揭示了碳纳米锥（CNCs）凭借其独特的凹凸曲面结构，可高效选择性吸附VX与A-234（Novichok）神经毒剂。研究表明，其吸附主要由范德华力驱动，并伴随显著的电子结构变化，展现出在化学防御（如选择性过滤与传感）领域的巨大应用潜力。
引言：化学防御的迫切需求与新材料探索
化学战剂（CWAs）对人类安全构成严重威胁，其中以VX和Novichok（A-234）为代表的有机磷神经毒剂因其极高的毒性、环境稳定性和水解抗性，其检测与消除是当前化学防御领域面临的重大挑战。传统的吸附材料如活性炭、石墨烯等，虽被广泛研究，但其相对均匀的电子景观限制了吸附强度和选择性。而金属有机框架（MOFs）等材料虽能实现强吸附甚至化学吸附，但可能存在湿度敏感、孔道阻塞等问题。因此，开发兼具高效、选择性吸附能力的新材料迫在眉睫。碳纳米锥（CNCs）作为一种具有固有曲率梯度和独特电子不对称性的新型碳纳米材料，在气体分子吸附与传感方面已显示出优异潜力，但其对高持久性神经毒剂的吸附性能尚未被深入探究。本文旨在通过系统的多尺度计算模拟，阐明CNCs对VX和A-234的吸附机理、选择性与潜在应用价值。
研究方法：多层次的计算模拟策略
本研究采用了一套层次化的计算工作流程，以全面评估碳纳米锥对神经毒剂的吸附行为。首先，通过分子对接和蒙特卡洛（MC）模拟进行全局构象搜索，快速探索毒剂分子在CNC表面的可能结合位点与取向。接着，对筛选出的最稳定结构进行密度泛函理论（DFT）几何优化与电子结构计算，以精确获得吸附能、电荷分布、前沿分子轨道（HOMO-LUMO）等关键参数，并分析非共价相互作用（NCI）。最后，利用分子动力学（MD）模拟引入时间维度，在有限温度下评估吸附结构的动态稳定性、扩散动力学及径向分布函数（RDF）。这种从快速筛选到高精度计算，再到动态验证的多尺度方法，旨在弥补单一计算层次的局限，为吸附机制提供更可靠的理论依据。计算模型中定义了两种吸附位点：凹面内表面（P1）和凸面外表面（P2）。
结果与讨论：曲率驱动的选择性吸附与强结合
3.1 吸附能量学与位点选择性
碳纳米锥独特的结构产生了两个截然不同的吸附区域：凹面内表面（P1）和凸面外表面（P2）。计算结果表明，这种曲率差异直接导致了VX和A-234的选择性吸附。VX分子更倾向于结合在凹面P1位点，而A-234分子则与凸面P2位点有更强的亲和力。分子对接结果显示，A-234在凹面位点的吸附能高达-91.4 kcal/mol，而VX在相同位点为-84.1 kcal/mol，表明CNCs的曲率增强了结合能。蒙特卡洛模拟进一步证实了这种选择性：在1:1吸附模型中，A-234的吸附能（Eads）为-76.85 kcal/mol，显著强于VX的-35.46 kcal/mol；在2:10（毒剂分子：CNC簇）模型中，A-234的吸附能更是达到-139.05 kcal/mol，远超VX的-57.15 kcal/mol。这表明CNCs对神经毒剂的捕获偏好顺序为：A-234 >> VX。竞争吸附模拟（引入水分子）显示，即使存在水分子竞争，两种神经毒剂在CNC表面的吸附能仍显著强于水，证明了其在潮湿环境下仍具备选择性吸附潜力。
3.2 吸附机制的本质
为探究强吸附的物理本质，进行了非共价相互作用（NCI）分析。降低密度梯度（RDG）图显示，在sign(λ2)ρ轴负值区域出现的绿色/蓝色峰，表明范德华（vdW）相互作用是CNCs与神经毒剂之间主要的吸附驱动力。分子动力学模拟得出的径向分布函数（RDF）图显示，毒剂分子中关键原子与CNC表面碳原子之间的接触距离约为3.5 ?，这比沙林在石墨烯上的报告距离（3.3 ?）更短，从结构上解释了本体系中更强的吸附能。这些结果共同表明，吸附是一种由范德华力主导的物理/化学吸附过程。
3.3 电子扰动与传感潜力
吸附过程显著改变了CNCs的电子结构，这为其用作传感材料提供了理论基础。静电势（ESP）分析显示，CNCs的尖端呈现亲核性（负ESP），而边缘呈亲电性（正ESP）。前沿分子轨道分析表明，纯净CNCs的HOMO-LUMO能隙（ΔE）为0.937 eV。吸附毒剂分子后，该能隙减小，尤其在凹面P1位点，VX吸附后体系的能隙降至0.882 eV，表明电荷转移程度较高，材料反应性增强。此外，吸附后体系的偶极矩显著增加，P1位点复合物的偶极矩（VX||CNCs约10.71 D，A-234||CNCs约10.54 D）远高于纯净CNCs（6.2854 D），表明体系极性增强。态密度（DOS）分析进一步证实了电子重排，吸附后在费米能级附近出现了新的电子态峰位移。这些因吸附引发的电子结构“指纹”变化，是CNCs具备高灵敏度传感潜力的关键。
3.4 吸附动力学
分子动力学轨迹分析揭示了吸附分子的扩散行为。计算得出的自扩散系数表明，A-234在CNC表面的迁移率（3.63 × 10-2?2/ps）远高于VX（1.7402 × 10-3?2/ps）。这意味着在潜在的CNC基过滤器或传感器中，A-234分子可能在表面相对自由地移动或平衡，而VX分子则被更牢固地束缚，近乎固定。这种动力学行为的差异对于理解毒剂在材料上的捕获与释放过程具有重要意义。
结论与展望
本项计算研究综合运用密度泛函理论、分子动力学和蒙特卡洛模拟，提供了详实的证据，表明碳纳米锥能够作为VX和A-234神经毒剂的高效分子陷阱。其前所未有的曲率结构是赋予其高效分子捕获能力的关键。凹凸分明的表面产生了独特的位点选择性：凹面内表面（P1）优先捕获VX，而凸面外表面（P2）对A-234展现出更强的吸附偏好。这种选择性吸附主要由范德华力驱动，并伴随显著的电子结构扰动（如HOMO-LUMO能隙减小、偶极矩增大、态密度变化），这为开发基于CNCs的化学传感技术奠定了理论基础。尽管模拟在气相中进行，未考虑湿度、竞争吸附等全部环境因素，且报告的吸附能为趋势性指标而非绝对结合自由能，但本研究明确揭示了CNCs作为一种新型吸附材料在化学防御领域（如选择性过滤、高灵敏度检测）的巨大应用潜力。未来的研究工作应聚焦于在更接近真实环境的条件下进行多组分模拟，并通过实验验证这些计算预测，以推动CNC基过滤器和传感器的实际开发。
