白云母的原子尺度建模及其结构中NH4+和K+扩散迁移的分子机制

《Applied Clay Science》：Atomistic modeling of phengite and molecular mechanisms of NH4 + and K+ diffusional mobility in its structure

　　

白云母（phengite）是一种重要的含水层状硅酸盐矿物，能够将水和钾（K）输送到俯冲板块中直至地球的上地幔。由于其结构能够容纳大量的铵（NH₄⁺），白云母也是俯冲带氮循环的主要贡献者。研究人员采用经典的分子动力学（MD）模拟，结合经验性ClayFF-MOH力场，研究了K-白云母和NH₄-白云母在常温常压以及高温高压条件下的行为。计算得到的K-和NH₄-白云母的结构和弹性性质与现有实验数据吻合良好。然而，模拟的主要焦点是研究驱动K⁺和NH₄⁺阳离子在白云母粘土层间空间中扩散的分子尺度机制。尽管在室温下未观察到阳离子的迁移，但模拟揭示这两种阳离子在高温下均可能发生扩散。此外，高温确保在高达10 GPa的压力下，两种阳离子在白云母晶体内部仍具有扩散迁移能力。另外，粘土层间即使存在少量H₂O，也会降低扩散的活化能垒，从而促进阳离子的迁移。尽管铵离子（NH₄⁺）尺寸较大，但在整个研究的温度和压力范围内，其扩散迁移率始终高于钾离子（K⁺）。

**论文解读：白云母中NH₄⁺与K⁺扩散迁移的原子尺度机制**

**一、研究背景与问题**

白云母（phengite）是一种高压含水层状硅酸盐矿物，在俯冲带中扮演关键角色。它能够在水和钾离子（K⁺）向上地幔输送过程中起载体作用，同时由于结构可容纳大量铵离子（NH₄⁺），白云母也是俯冲带氮循环的重要参与者。然而，俯冲带高温高压条件下白云母中阳离子的迁移行为，尤其是NH₄⁺和K⁺在层间空间中的扩散机制，仍缺乏分子尺度的定量理解。实验手段（如光谱学、衍射）虽能提供宏观信息，但难以直接揭示原子级别的动力学过程。因此，研究人员采用经典分子动力学（MD）模拟，结合改进的ClayFF-MOH力场，系统研究了K-白云母和NH₄-白云母在宽温压范围（室温至1073 K，常压至10 GPa）内的结构、弹性及阳离子扩散行为，旨在阐明俯冲带条件下氨释放与钾迁移的分子机制。该论文发表于《Applied Clay Science》。

**二、主要关键技术方法**

研究人员采用经典分子动力学（MD）模拟，利用改进的ClayFF-MOH力场描述原子间相互作用。模型基于天然2M₁多型白云母的中子衍射结构数据（Mookherjee et al., 2001），通过随机替换Mg²⁺、Fe³⁺、Al³⁺进入八面体和四面体片层构建超胞（64个晶胞），并引入K⁺或NH₄⁺平衡净负电荷。模拟采用NPT和NVT系综，时间步长1 fs，使用LAMMPS和PLUMED软件包。为计算阳离子扩散的活化能垒，应用了偏置元动力学（well-tempered metaDynamics）增强采样方法，以层间阳离子与空位六方环质心的距离为集体变量。此外，还拟合了Birch-Murnaghan状态方程以获取弹性模量和热膨胀系数。

**三、研究结果**

**3.1 白云母的晶格参数与压缩性（Crystallographic lattice parameters and the compressibility of phengite）**
通过MD模拟在无对称约束下计算的晶胞参数与实验数据吻合良好，K-白云母体积偏差小于2.5%，NH₄-白云母偏差在2.4%-4.6%之间。热膨胀模拟显示，高温下计算值与实验值的差异减小（从室温的2%降至873 K的1.2%）。等温压缩性表明，NH₄-白云母比K-白云母更易压缩（结构更柔性），且压缩性随温度升高而增加。在压力高达17 GPa时，K-白云母的模拟结果与实验偏差不超过2%。

**3.2 白云母层间阳离子的局部结构排列（Local structural arrangement of the interlayer cations in phengite）**
通过径向分布函数（RDF）分析，K⁺与表面氧（O_b）的平均距离为3.10 ?，与羟基氧（O_h）的距离为3.90 ?，均与实验一致。高压（10 GPa）下K-O距离缩短至2.81 ?。对于NH₄-白云母，铵氢（H_A）与O_b和O_bts（四面体取代桥氧）的RDF显示强氢键（H键）存在（距离<2.4 ?）。高温下H键减弱，高压下增强。此外，高温下铵氢可与羟基氧形成弱H键，可能促进质子转移和铵分解。

**3.3 白云母的弹性性质与体积模量（Elastic properties of phengite and the bulk modulus）**
利用Birch-Murnaghan方程拟合P-V数据，计算得到室温下K-白云母的体积模量为63.0 GPa，NH₄-白云母为58.1 GPa，后者更易压缩。体积模量随温度线性下降，1073 K时两者均降低15%-20%。整个温度范围内NH₄-白云母的模量始终低于K-白云母。

**3.4 热膨胀系数（Thermal expansion coefficient）**
计算得到的体积热膨胀系数K-白云母为4.4×10^-5 K^-1，NH₄-白云母为4.9×10^-5 K^-1。NH₄-白云母的热膨胀略高，表明热性质主要由层间空间性质而非TOT结构单元决定。

**3.5 白云母层间阳离子迁移（Cation mobility in phengite interlayers）**
**3.5.1 脱水白云母（Dehydrated phengite）**
元动力学模拟计算了阳离子在相邻六方环间跳跃的自由能垒（ΔF）。室温下K⁺和NH₄⁺的ΔF分别为193和188 kJ·mol^-1，无扩散。高温（≥673 K）下扩散发生，常压时K⁺的ΔF为264 kJ·mol^-1，NH₄⁺为246 kJ·mol^-1，后者略低。1073 K、10 GPa时仍可扩散。扩散系数（D）的Arrhenius图显示NH₄⁺的D高于K⁺，但模拟值系统性低于实验值，归因于模型未考虑边缘效应和化学反应。
**3.5.2 水化白云母（Hydrated phengite）**
在873 K、5 GPa条件下，加入少量水（0.5 H₂O/阳离子）显著降低了阳离子扩散的自由能垒，表明层间水促进了阳离子迁移。但定量评估需引入更多集体变量。

**四、讨论与结论**

研究人员利用经典MD模拟系统研究了K-和NH₄-白云母在高温高压下的结构、弹性及阳离子扩散行为。模拟结果与实验在结构和弹性性质上吻合良好，验证了模型的可靠性。研究揭示了阳离子扩散仅在高温下发生，且压力升高会增大活化能垒，但足够的热能仍可克服（如1073 K、10 GPa）。水分子可显著降低扩散能垒。尽管NH₄⁺尺寸较大，其始终表现出比K⁺更高的层间扩散迁移率，这归因于K⁺更强的束缚。

**结论（翻译原文Conclusions部分）**：
在这项工作中，研究人员利用经典分子动力学模拟，结合改进的ClayFF力场（Cygan et al., 2021），研究了K-白云母和NH₄-白云母在高温高压下结构和弹性性质的演化。与实验的良好一致性使研究人员能够探究K⁺和NH₄⁺阳离子在白云母层间空间中的扩散迁移。计算揭示，阳离子扩散仅在升温时发生，在环境条件下未观察到。NH₄⁺扩散的活化能垒与实验中铵损失所确定的活化能吻合良好（Liu et al., 2019b）。模拟结果也证实了两种阳离子扩散系数的Arrhenius温度依赖性。然而，计算得到的扩散系数系统性地低于实验值，这种低估可能是模型局限性的结果，因为计算仅考虑了矿物体相层间的阳离子扩散，未考虑粘土纳米颗粒边缘处层间阳离子的损失（最高可达50%），也未考虑铵阳离子与粘土表面发生化学反应的可能性。尽管压力升高会增加阳离子扩散的活化能垒，从而阻碍迁移，但在足够高的温度下，热能仍可促进扩散，即使在高达10.0 GPa的压力下。白云母层间引入少量H₂O分子会降低阳离子扩散的活化能垒，促进其迁移。然而，尽管铵离子尺寸更大，其层间扩散率始终高于钾离子，这可归因于钾离子在粘土TOT层内更强的束缚或更大的限制。定性而言，模拟结果提供了清晰的分子尺度视角，揭示了含水粘土矿物中氨的释放如何强烈影响其在俯冲带中的脱水作用（例如，Liu et al., 2019b）。然而，对这些机制的详细定量表征需要更广泛的分子模拟，覆盖更宽的温压范围和白云母组成。尽管如此，即使当前的有限模拟已经表明，在俯冲带条件（高达10 GPa和1073 K）下，NH₄⁺在白云母层间的迁移率始终高于K⁺。这意味着一旦铵被结合，在进变质作用过程中可能比钾更容易释放，从而可能解耦N和K循环。此外，即使存在少量层间水，其显著降低活化能垒的作用表明，脱水反应和结构结合水的释放可能触发或加速铵的损失。这些分子尺度发现为解释折返的高压/超高压地体中氮同位素分馏和流体活动元素富集的地球化学观测提供了机理基础。