蒙脱石吸水分子动力学模拟 原创

关键词：蒙脱石；分子动力学；水化机理；层间距；膨胀土

一、研究背景与问题

蒙脱石是膨胀土的主要矿物成分，具有吸水膨胀、失水收缩的特性，在路基、边坡和地基中常引发开裂、隆起等工程病害。传统试验难以在纳米尺度揭示其水化机制。本文采用分子动力学方法，从微观角度揭示了Na基蒙脱石（Na-MMT）的吸水机理。

二、模型构建与仿真方法

基于Wyoming型蒙脱石晶胞（空间群C2/m，a=0.523nm、b=0.906nm、c=0.960nm），构建8晶胞超胞（4a×2b×c），按Loewenstein随机取代规则引入Mg²⁺/Al³⁺取代，化学式为Na6[Si62Al2]Mg4Al28O160(OH)32·nH2O，水分子采用SPC/E模型。在Materials Studio的Forcite模块中使用UFF通用力场，库仑作用取Ewald求和、范德华作用截断半径9Å；采用NPT系综，时间步0.5fs，温度298K，模拟总时长1000ps，前160万步用于平衡，最后200ps做统计平均。

三、关键结果与发现

图1给出蒙脱石超胞吸附模型的XOY表面。几何优化后多面体发生明显畸变，层间Na⁺更贴近硅氧环与均质取代位置，水分子由中部集中分布转为分层分布，说明层间电荷与位置主导着粒子的排布。

图1 蒙脱石超胞吸附模型XOY表面（优化前后对比）

图2为层间距随吸附水分子数的变化。层间距呈阶梯式增长：吸附0–16与32–48个分子时增速明显，分别对应第一、二层水开始吸附；16–32与48–80个分子时增速放缓，表明水层趋于饱和。吸附32、64、96个分子分别对应形成一、二、三层水分子层。

图2 Na-MMT层间距随吸附水分子数呈阶梯式增长

图3为层间Na⁺与水的径向分布函数。gNa⁺–Ow(r)在2.1Å与5Å附近出现双峰，氢键作用范围0.26–0.31nm，表明水分子与粘土层表面原子形成了氢键，这是水化作用的结构基础。

图3 层间Na⁺-Ow径向分布函数RDF

图4为Na⁺与水分子的均方位移MSD。三层水时扩散最强、二层最弱；且Na⁺的迁移始终慢于水分子，源于其与带电粘土表面的强静电作用显著拖慢了金属离子扩散。

图4 Na⁺与水分子均方位移MSD（一/二/三层水）

图5进一步分解2w状态水分子的MSD。z方向扩散几乎为零，受粘土骨架氢键与桥氧限制；粒子主要在层平面内扩散，体现出明显的各向异性。

图5 2w状态水分子MSD各向异性（z方向受限）

图6给出不同水化状态下层间Na⁺、O、H的相对浓度分布。随含水量增加，分布由单峰转为多峰且更分散对称，说明水分子在上下层间分布趋于均匀。

图6 不同水化状态层间Na⁺/O/H相对浓度分布

四、我们提供的服务

① 蒙脱石及粘土矿物分子动力学建模（Materials Studio Forcite / LAMMPS），含力场与模型搭建； ② 水化膨胀与层间距演化的仿真预测，定量给出吸水曲线与膨胀趋势； ③ 水–岩及离子相互作用的力场参数标定与验证； ④ 化学改良（无机盐/有机剂）抑制膨胀方案的仿真比选； ⑤ 膨胀土工程病害机理的微观仿真咨询与报告撰写。 如需针对具体土样或工程场景开展仿真，欢迎与我们联系。

如有需要，欢迎通过公众号"320科技工作室"联系我们。