关键词：分子模拟；元动力学；CO2矿化；ReaxFF；水膜
一、文章简要介绍
这篇Nature Communications文章研究二氧化碳在润湿矿物表面的矿化机制，来自法国、英国和香港理工大学的联合团队，2025年发表，DOI：10.1038/s41467-025-65794-w。矿物碳化是潜在的十亿吨级碳汇，水泥混凝土每年能吸收大量CO2，但自然风化碳化太慢。文章用分子模拟，把氢氧化钙介孔里CO2矿化的完整路径一步步算清楚。
矿化发生的场所很特别：潮湿环境下，矿物表面会吸附一层纳米厚的水膜，它就像一个纳米反应器。实验发现相对湿度超过50%碳化才明显加快。研究者选了氢氧化钙（水泥水化的主要钙源）做对象，重点盯住它的基面（001面）。
图1 界面CO2矿化的完整模拟路径：水吸附、离子溶解、CO2溶剂化、碳物种转化、离子配对
图1把整条路径串起来看：先吸水成膜，再溶出钙离子，然后CO2进入水膜、转化成碳酸根，最后和钙离子配对成碳酸钙。五个步骤环环相扣。
二、分子模拟步骤
步骤一：建模型。氢氧化钙原胞扩成超胞，沿基面劈开造出5纳米狭缝孔。为了模拟真实矿物表面的粗糙度，做了两套模型：一套平整表面，一套挖掉六个氢氧化钙单元形成台阶缺陷。水膜厚度按相对湿度90%用巨正则蒙特卡洛（GCMC）定。
图2 氢氧化钙介孔吸附水膜：相对湿度约90%时形成约1纳米厚水膜（GCMC模拟）
图2是水吸附的结果：相对湿度大于30%才形成单层水膜，90%时约三层、1纳米厚。这个厚度下的矿化活性最好。
步骤二：选力场。吸附水用Clay-FF力场配SPC/E水模型；反应过程用ReaxFF反应力场，它能描述化学键的断裂与生成。所有反应模拟在NVT系综300开尔文下跑，时间步0.25飞秒，先弛豫1纳秒，再生产3到10纳秒。
步骤三：加元动力学加速采样。矿化反应是稀有事件，普通模拟等不到。文章用well-tempered元动力学（WT-MTD），在PLUMED软件里给体系加偏置势，沿反应坐标采样，把自由能面算出来。
图3 Ca2+溶解的自由能面：台阶位与完整位溶解路径对比（元动力学模拟）
图3是钙离子溶解的自由能面。台阶位的钙只连四个羟基，比完整位少两个约束，能探索的空间大得多。台阶位要先爬0.49电子伏的能垒到亚稳态，之后一步0.09电子伏就能完全脱附，形成(CaOH)+水合离子。
图4 CO2从气相进入水膜的溶剂化路径与自由能变化
图4看CO2怎么进入水膜：从气相穿过气-水界面到膜内，再到矿物表面，每穿透一层水都要付出能量代价。
三、关键结果与发现
第一个结果：CO2的转化路径分两种。水合路径在气-水界面更容易发生，CO2直接和水反应先生成碳酸；羟基化路径则是CO2与溶解出来的OH-反应，能垒高很多，多发生在矿物表面附近。碳酸再逐步去质子化变成碳酸氢根、碳酸根。
图5 CO2水合与碳酸去质子化的自由能面：协同与分步两种水合方式
图5是碳物种转化的自由能面：中性和酸性膜里CO2主要走水合路线，碱性膜里CO2更愿意潜进水里找OH-反应。去质子化的能垒在薄膜里比体相水高，界面环境确实会改变反应。
第二个结果：离子配对三条路都能走。碳酸钙成核从钙离子和碳酸根配对开始。强碱膜里CO3直接结合Ca2+；中性膜里Ca2+和HCO3一边结合一边脱质子；酸性膜里则从H2CO3开始。不管哪种环境，碳酸钙总能生成。
图6 CaCO3离子配对的三种路径：Ca2+-CO3、Ca2+-HCO3、Ca2+-H2CO3
图6是三条配对路径的自由能面。Ca2+和CO3结合只须爬0.19电子伏，放出0.24电子伏，配对时挤掉六个水分子，熵增让反应更顺。
第三个结果：pH是把双刃剑。酸性利于钙溶解（实验里pH=1比pH=7快一个数量级），碱性利于CO2转化，两者互相竞争。近中性到微碱的水膜两边都兼顾，矿化效率最高。实验也验证：中性雾环境下碳化程度最大，强碱反而生成大量无定形碳酸钙。
图7 pH调控矿化机理总览与碳化实验验证：近中性条件矿化程度最高
图7把pH的作用串成完整图景：溶解抬升碱度、带动碳转化，碳转化又消耗碱度、回头带动溶解，离子配对再把产物沉淀掉，形成一个能量上划算的循环，让矿化持续进行。CO2 speciation是限速步骤，不是钙溶解。
文章还提示：矿化路径和实验活化能趋势吻合，说明这套从水吸附到离子配对的框架可以推广到硅酸盐等其他含钙矿物，为加速碳化固碳提供设计依据。
四、我们提供的服务
我们可提供以下模拟服务：①LAMMPS分子动力学与反应力场模拟；②GCMC吸附等温线计算；③元动力学增强采样与自由能面分析；④矿物-水界面反应机理研究。如有类似模拟需求，欢迎私信聊聊。