柠檬酸功能化碳纳米管对稀土元素吸附驱动传输的位点特异性选择性：一项基于第一性原理的计算研究

《The Journal of Physical Chemistry C》：Site-Specific Selectivity of Citrate-Functionalized Carbon Nanotubes for Adsorption-Driven Transport of Rare Earth Elements: An Ab Initio Computational Study

　　

稀土元素（REEs）的分离因其配位化学与水合行为高度相似而仍是开放难题。钕（Nd）与钇（Y）构成典型模型体系：Nd是废弃磁体中优先回收目标，Y在离子半径与配位行为上模拟重稀土，常作为Nd共污染物。基于膜的分离被视为替代传统溶剂萃取的高效低耗途径，而开发具特定选择性的新型膜涂层依赖于带定制结合位点的纳米结构。研究人员采用密度泛函理论（DFT）阐明Y3+与Nd3+在柠檬酸功能化碳纳米管（CNT-Cit）上的吸附行为，以评估该功能化纳米结构作选择性膜涂层的潜力。研究人员系统计算吸附能（ΔEads）、结合能（BECOO–）、形变能（ΔEdef）及吸附自由能，量化金属配合物随结合拓扑变化的吸附与稳定性。研究确立CNT-Cit上REEs结合的结构—能量关系：吸附强度与配位几何强烈依赖金属离子与功能化拓扑。具体而言，管壁位点经柠檬酸中心羧基的吸附对Nd为放能（exergonic）、对Y为吸能（endergonic）；管尖位点吸附过强不利于离子传输。管壁柠檬酸基团对Nd与Y的一放一吸表明CNT-Cit对Nd优于Y具选择性。虽未直接模拟组装CNT网络中位点间跳跃势垒或通量，研究为设计选性CNT基网络作REE回收膜涂层提供分子级准则。

研究背景与动机：稀土元素（Rare Earth Elements, REEs）支撑异相催化、永磁体、荧光粉与光学器件等技术，但从二次资源回收受限于REEs化学行为相近、离子尺寸与水合结构差异细微，传统溶剂萃取、离子交换与沉淀路线试剂密集且二次废物多。吸附策略因可模块化与低影响受关注，而高选择性而非仅高容量是可持续分离的核心瓶颈。钕（Nd）与钇（Y）被选为模型体系：Nd属轻稀土（LREE）源自废磁体，Y离子半径与配位行为模拟钇（Ho–Er）等重稀土（HREE），二者分离可检验吸附剂在LREE/HREE范围内的选择性控制能力。碳纳米材料具高比表面与可修饰界面，柠檬酸为多羧基配体，提供多个氧供体与构象柔性，共价锚定于CNT后引入稳定螯合位点。在膜电容去离子（MCDI）或电渗析（ED）中，功能化CNT涂层可作阳离子交换膜（CEM），其分离机制更可能走表面介导传输（surface-diffusion）：离子在功能位点吸附，经互联CNT网络hopping迁移，再由对侧释放，而非单纯尺寸排阻。研究人员据此用第一性原理计算评估CNT-Cit在水相下位点特异性吸附热力学是否足以产生REE选择性，而不直接模拟膜内通量。论文发表于《The Journal of Physical Chemistry C》。

主要技术方法：研究人员使用Northwest Computational Chemistry Package（NWChem）套件，采用B3LYP交换关联泛函与6-311G*基组描述C/H/O，Nd与Y采用Stuttgart RSC 1997有效核势（ECPs），引入Grimme DFT-D4色散校正；水相以COnductor-like Screening Model（COSMO）连续介质模拟（介电常数80.2），金属首层水壳显式建模。碳纳米管取zigzag (9,0)单壁CNT（直径0.71 nm，长约2.5 nm，216碳），端氢饱和；侧壁功能化于空位缺陷处酯化柠檬酸侧羧基得CNT_wall-Cit，管尖功能化得CNT_end-Cit。采用静电嵌入簇方法：量子区含柠檬酸及吸附点附近碳，外区CNT原子以羟基化模型静电势（ESP）点电荷固定。吸附能ΔE_ads依优化复合物、游离CNT-Cit与六水合离子能差计算；结合能BE_COO–取碎片能差；形变能ΔE_def分CNT-Cit与金属水合物贡献；吸附自由能近似计平移与构型熵，假设吸附不释放水壳水分子故忽略脱水利熵。

研究结果：

3.1 CNT–柠檬酸模型几何：原始CNT端氢饱和，C–H约1.10 ?。CNT_wall–OH借移除碳造空位使柠檬酸酯化，锚定C–O_Cit键1.39 ?，C_anchor–O–C_Cit角118°；CNT_end-Cit相应1.40 ?与117°，酯C–O取向平行管轴，侧壁取向垂直指向空位。侧壁柠檬酸侧羧基与羟基氧距管壁碳2.8–3.1 ?，管尖因向外伸展无此接触。

3.2 CNT-Cit金属配合物几何：柠檬酸均作双齿配体仅经羧基氧配位。Y配合物配位数8，保留满水壳；Nd配合物多置换1–3个内层水至次层氢键网络（如CNT_wall-Cit_c-Nd换3水）。Y–O_water 2.3–2.6 ?短于Nd–O_water 2.1–2.9 ?分布，符Y³⁺半径小电荷密度高。但羧基键长反转：Y–O_COO对称双齿2.2–2.4 ?，bite角55–58°；Nd呈不对称双齿，仅CNT_end-Cit_l-Nd近对称（1.89、1.94 ?，bite角66.5°），余者一短一长。

3.3 吸附、结合与形变能：CNT_wall-Cit_c位点Nd的ΔE_ads = –11.19 kcal/mol（放能），Y为+7.84 kcal/mol（吸能）；CNT_wall-Cit_l二者皆正（Nd +10.01，Y +33.96）；CNT_end-Cit_c Y –2.08、Nd +8.41；CNT_end-Cit_l Nd –61.45、Y –26.91 kcal/mol，吸附过强不利传输。BE_COO–上Nd（–175至–235 kcal/mol）远强于Y（–63至–101 kcal/mol）；ΔE_def中Nd水合物形变（154–196 kcal/mol）显著高于Y（65–95 kcal/mol），说明Nd为成键付出更大几何重排代价却被更强静电补偿。管壁中心羧基位点因而构成Nd放能/Y吸能的判别界面。

讨论与结论翻译：研究人员确立CNT-Cit上REEs结合的结构—能量关系，显示吸附强度与配位几何强烈依赖金属离子与功能化拓扑。管壁位点经柠檬酸中心羧基的吸附对Nd为放能、对Y为吸能，而管尖位点吸附过强不宜离子传输。管壁柠檬酸对Nd与Y的一放一吸预示CNT-Cit对钕优于钇具选择性。尽管未直接模拟组装CNT网络中位点间跳跃势垒或通量，该研究为设计选择性CNT基网络用作REE回收膜涂层提供了分子级准则。

意义：论文以Nd/Y为探针，从DFT层面证明羧基拓扑（管壁中心vs侧位vs管尖）可翻转三价稀土吸附热力学符号，指出表面扩散机制下“可逆放能”位点（CNT_wall-Cit_c）是膜涂层选性来源，为重稀土模拟离子Y与轻稀土Nd的识别提供定量能量判据，避开直接动力学模拟而给出涂层设计的热力学边界。