理论价值填补分子尺度机制空白突破传统研究 “重宏观、轻微观” 的局限首次通过分子动力学MD模拟将轻木生物炭吸附微塑料的机制从 “宏观吸附容量描述” 推进到 “原子级作用力量化”明确范德华力80~90%为核心驱动力、静电作用次要贡献、氢键可忽略的作用规律解决了 “多作用力耦合下主导机制不明” 的关键科学问题。建立 “结构 - 性能” 理论关联通过构建生物炭全原子模型BC300 用萘单元、BC700 用 HBC 单元量化了 “表面含氧官能团含量→相互作用能强度→吸附容量”“孔隙结构发达程度→扩散系数→吸附动力学” 的构效关系为碳基吸附材料的分子设计提供理论范式。验证 MD 模拟的可靠性实现宏观实验吸附容量顺序 PSPAPE与微观模拟相互作用能顺序 PSPAPE的高度吻合证明 MD 模拟可作为 “预测吸附性能、解析作用机制” 的有效工具为后续微塑料吸附领域的理论研究提供方法参考。应用价值指导生物炭定向优化明确 BC300 因高含氧官能团更适合 “高容量吸附”BC700 因高扩散速率更适合 “快速动态处理”为实际污染治理中 “按需选择生物炭热解温度” 提供依据同时提示 “强化表面 π 电子位点如保留芳香结构” 可提升 PS 类微塑料去除效率为生物炭改性提供靶向方向。降低污染治理成本基于轻木低成本可再生原料的研究基础结合 MD 模拟对 “吸附机制的精准解析”可减少传统 “试错式改性” 的实验成本加速生物炭从实验室走向实际水体微塑料修复的产业化进程。文献介绍近日内蒙古农业大学刘文静教授团队材料科学与艺术设计学院 / 内蒙古沙生灌木纤维化及能源开发利用重点实验室在《Separation and Purification Technology》分离纯化领域权威期刊IF9.1JCR Q1 区发表关于 “轻木生物炭吸附微塑料的分子动力学模拟机制” 研究。研究核心是通过宏观吸附实验与分子动力学MD模拟深度耦合首次从原子尺度阐明轻木生物炭300℃/700℃热解对聚苯乙烯PS、聚酰胺PA、聚乙烯PE三种典型微塑料的吸附驱动力、扩散规律及热力学本质破解了长期以来生物炭吸附微塑料的 “分子黑箱”为高效碳基吸附剂的设计与应用提供关键理论支撑。背景介绍微塑料MPs作为 “环境隐形杀手”已广泛存在于海洋、陆地、大气中通过食物链累积威胁生态稳定与人类健康开发高效去除技术迫在眉睫。现有微塑料去除技术中吸附法因 “操作简便、无二次污染、速率快” 成为主流但核心瓶颈在于吸附机制不明导致吸附剂优化无章可循。生物炭虽因 “低成本、可再生、表面可调” 成为潜力吸附剂但过往研究仅能通过宏观实验如吸附等温线、动力学模型推测 “可能存在范德华力 / 静电作用”无法量化各作用力贡献占比更无法解释 “为何 300℃热解生物炭吸附容量更高而 700℃生物炭扩散更快” 等构效关系问题。传统实验手段如 XPS、BET仅能捕捉 “吸附前后的静态变化”无法观测 “微塑料分子在生物炭孔隙中运动、结合的动态过程”—— 而分子动力学MD模拟作为 “原子级显微镜”可通过构建全原子模型、计算相互作用能、追踪分子轨迹从动态视角解析吸附机制成为突破这一研究瓶颈的核心技术。理论研究方法MD 模拟核心设计本研究的核心亮点在于 **“精准建模 关键参数量化”**从理论计算角度构建了可复现、可推广的研究框架具体设计如下图 1 两种生物炭模型的康诺利孔隙表面及展示单壁碳纳米管引入的中间过程1. 生物炭全原子模型构建匹配实验特征建模方法采用 Segarra-GlandtSG模型基于 Materials StudioMS软件构建 50×50×50 Å 的模拟单元以 “芳香烃单元” 为基本 building blockBC300 用萘2 个芳香环、BC700 用六苯并苯HBC6 个芳香环通过调整单元密度与芳香环数量调控孔隙结构。官能团校正通过 Perl 脚本引入 C-C、C-O、CO 官能团比例严格匹配 XPS 实验结果BC300C-O 20%、CO 4%BC700C-O 8%、CO 4%确保模型与实际生物炭化学特征一致。结构优化先通过 DMol³ 模块进行密度泛函理论DFT原子级优化再用 Forcite 模块的 NPT 系综Andersen 控温、Berendsen 控压退火处理使模型达到热力学平衡。图2 300℃热解轻木生物炭BC300与微塑料溶液的固液双相模型。a聚苯乙烯PSb聚酰胺PAc聚乙烯PE2. 微塑料模型与固液两相体系构建微塑料建模采用全原子模型比粗粒化模型更精准基于 PS、PA、PE 的分子结构绘制单体通过 MS 聚合功能连接 10 个单体形成线性链匹配实验中 18μm 微塑料的分子尺度特征。溶剂相设计按宏观实验浓度调整水分子比例构建 50×50×10 Å 的水相单元通过 “Build Layers” 工具将生物炭模型与水 - 微塑料体系组合引入真空层形成固液两相界面模拟实际水体环境。3. MD 模拟参数与关键指标计算力场与系综选用 COMPASS 力场适用于有机高分子与碳材料体系采用 NPT 系综温度 25~60℃、压力 1atm施加周期性边界条件消除边界效应。核心参数计算相互作用能EinterEvdwEelec量化范德华力与静电作用的贡献占比氢键分析基于 “键长 3.8 Å、键角 90°” 标准统计氢键数量与强度吸附等压热Qst通过公式计算判断吸附类型40kJ/mol 为物理吸附扩散系数D从分子均方位移MSD推导反映微塑料在生物炭中的运动速率。结果和讨论从 MD 模拟视角解读核心发现吸附驱动力范德华力主导静电作用辅助图3 两种生物炭与微塑料MPs的相互作用能a300℃热解轻木生物炭BC300b700℃热解轻木生物炭BC700MD 模拟显示两种生物炭与微塑料的总相互作用能Einter均为负值表明自发吸附且范德华力Evdw占比 80~90%是吸附的核心驱动力。原因解析生物炭的芳香环结构与微塑料分子尤其是 PS 的苯环形成 π-π 堆叠同时碳骨架与微塑料的疏水区域产生疏水相互作用共同贡献范德华力静电作用的特殊性虽生物炭与微塑料均带负电Zeta 电位实验验证但分子尺度下局部电荷不均仍产生微弱静电吸引贡献 10~20% 相互作用能且 BC300 因高含氧官能团C-O/CO的极性静电作用略强于 BC700。图4 模拟过程中氢键键长与键角的分布概率。a、b300℃热解轻木生物炭BC300c、d700℃热解轻木生物炭BC700此外氢键分析显示图 4仅 PA 与生物炭的含氧官能团存在极低概率氢键键长≈3.8 Å、键角 120°对总相互作用能贡献 1%可忽略不计。构效关系官能团决定 “吸附容量”孔隙决定 “扩散速率”图5 C1s 轨道的 X 射线光电子能谱XPS峰拟合结果。a300℃热解轻木生物炭BC300b700℃热解轻木生物炭BC700图 6 两种生物炭与微塑料MPs的相互作用能a300℃热解轻木生物炭BC300b700℃热解轻木生物炭BC700吸附容量差异的分子机制BC300 的 Einter 绝对值如与 PS 作用能≈-120kJ/mol显著高于 BC700≈-80kJ/mol对应宏观实验中 BC300 吸附容量更高 —— 核心原因是 BC300 保留更多 C-O/CO 官能团增强了与微塑料的范德华力和静电作用图 5、图 6。图7 体系内的均方位移与扩散系数。a300℃热解轻木生物炭BC300b700℃热解轻木生物炭BC700扩散速率差异的原子解释BC700 的扩散系数如 PS 在 BC700 上 D≈1.85×10⁻⁵cm²/s是 BC300D≈1.57×10⁻⁵cm²/s的 1.18 倍图 7—— 因 700℃热解使生物炭形成石墨微晶孔隙结构更发达比表面积≈429.65m²/g是 BC300 的 60 倍降低了微塑料分子的迁移壁垒且 PS 因与生物炭表面 π-π 作用强扩散速率最快PA 次之PE 最慢。3. 热力学本质物理吸附主导温度影响有限图8 摩尔吸附热随温度的变化曲线。a300℃热解轻木生物炭BC300b700℃热解轻木生物炭BC700吸附等压热Qst计算结果显示图 8两种生物炭对三种微塑料的 Qst 均 40kJ/molBC300 的 Qst≈28kJ/molBC700≈12kJ/mol证明吸附过程为物理吸附无化学键形成易脱附再生。温度效应随温度升高25~60℃Qst 略有上升BC300 从 25kJ/mol 升至 28kJ/mol因温度升高为微塑料分子提供更多动能帮助克服吸附能垒 —— 但 Einter 仅提升 5~8%说明温度对吸附性能影响有限范德华力仍为绝对主导这与宏观实验中 “温度升高吸附容量略有增加” 的结果完全吻合。4. 微塑料吸附顺序PSPAPE模拟与实验完美匹配图9 环境因素对吸附容量的影响。a、b微塑料浓度c、d吸附时间e、f吸附温度MD 模拟中图6两种生物炭对微塑料的相互作用能顺序均为 PSPAPE如 BC300 与 PS 作用能≈-120kJ/molPA≈-100kJ/molPE≈-80kJ/mol与图9宏观实验中 “吸附容量 PSPAPE” 的顺序一致。分子层面解释PS 含芳香环与生物炭形成强 π-π 作用PA 含极性酰胺基-CONH-但生物炭表面疏水削弱了相互作用PE 为非极性聚合物无功能性位点仅靠弱疏水作用结合 —— 这一解释为 “不同类型微塑料的选择性去除” 提供了理论依据。总结本研究通过MD 模拟与宏观实验的深度耦合实现了三大突破① 首次量化轻木生物炭吸附微塑料的作用力贡献范德华力 80~90% 主导② 建立 “生物炭分子结构官能团 / 孔隙→微观参数相互作用能 / 扩散系数→宏观性能吸附容量 / 速率” 的关联模型③ 验证了全原子 MD 模拟在微塑料吸附领域的可靠性为该领域理论研究提供可复现的方法学框架。展望当前研究仍存在 “理想化模型” 局限① 未考虑纳米塑料比表面积更大作用力更强的吸附机制② 未引入 pH、离子强度、溶解有机质等实际水体因素如高 pH 可能削弱静电作用芳香有机质可能竞争 π-π 位点。未来方向应聚焦① 开展 “量子化学DFT 分子动力学MD 介观模拟” 的多尺度研究解析更复杂的作用机制② 结合多因素实验如不同 pH、离子强度构建 “贴近实际环境” 的 MD 模型推动理论研究向产业化应用转化。文章信息原文标题Molecular dynamics simulation of microplastics adsorption mechanisms on balsa biochar标题中文翻译轻木生物炭吸附微塑料的分子动力学模拟机制期刊Separation and Purification Technology分离纯化领域 Q1 区IF9.1DOIRedirecting通讯作者单位内蒙古农业大学材料科学与艺术设计学院内蒙古沙生灌木纤维化及能源开发利用重点实验室呼和浩特010018