基于密度泛函理论的氢氧化镁（101）表面改性机理研究

李立清， 周润， 龙慧婷， 孔会民， 邹来禧*， 罗仙平， 邵羽南， 杨国飞

（1. 江西理工大学材料冶金化学学部，江西 赣州 341000；2. 西部矿业集团有限公司青海省高原矿物加工工程与综合利用重点实验室， 西宁 810000；3. 江西同益高分子材料科技有限公司，江西 赣州 341000； 4. 江西宏特绝缘材料有限公司，江西 抚州 344000）

氢氧化镁（Mg（OH）2）因兼具无毒、无味、无腐蚀、阻燃、消烟、热分解温度高（340～450 ℃）和环保等优点，已成为高分子材料阻燃领域中极具前途的阻燃剂之一[1]。但是，目前研发的氢氧化镁阻燃剂还存在如下不足[2]：①产品粒径大，阻燃效率低，实际应用时需大量填充，严重损害材料的力学性能。②产品粒度不均匀，团聚现象严重，不利于在高分子材料中分散。③产品表面极性强，与有机高分子的相容性较差。④以水氯镁石制备的氢氧化镁，氯含量高，作为阻燃剂会产生HCl 等有害气体。为弥补上述不足，研究者从制备方法和表面改性方面，对氢氧化镁的形貌、粒径和表面极性等进行了研究。

在改性方面，目前大部分研究集中于改性剂的选择，在制备过程中，通过表面活性剂（硬脂酸钠等）、偶联剂（有机硅烷等）以及表面聚合包覆（聚合物等）降低表面极性和增加与有机高分子的相容性。PIPEROPOULOS 等[3]研究了表面活性剂十六烷基三甲基溴化铵（CTAB）对氢氧化镁结构、形貌和物理性能的影响，结果表明，当CTAB浓度为2.0 mmol/L时，促进形成了分离良好的氢氧化镁颗粒，降低了氢氧化镁的平均粒径；ZHU 等[4]采用钛酸酯偶联剂改善氢氧化镁与聚对苯二甲酸乙二醇酯（PET）的相容性，结果表明，直径为0.1～0.5 μm 的PET 与改性纳米Mg（OH）2颗粒（改性剂含量：10%、20%和30%）相容性良好；YANG 等[5]以合成的碳微球包覆氢氧化镁提高其阻燃性能，结果表明，氢氧化镁和碳微球（CMSs）对PET 具有协同阻燃作用。虽然这些方法可以在一定程度上改善氢氧化镁的粒径分布和表面极性，但对氢氧化镁生长机制、粒径控制机理和改性剂机理的研究较少，开发存在一定的盲目性。

密度泛函理论（Density Functional Theory，DFT）近几年在材料合成、模拟计算等方面的应用研究取得了较大进展，成为计算凝聚态物理、计算量子化学和计算材料科学的重要根基和中心技术，其应用范围越来越广[6-8]。本研究基于DFT 方法，选取了3 种常用的氢氧化镁改性剂：油酸[9]、硅烷偶联剂G-570（γ-甲基丙烯酰氧基丙基三甲氧基硅烷）以及甲基丙烯酸甲酯（MMA）[10]，计算模拟了这3 种物质在氢氧化镁表面的吸附过程，探究氢氧化镁（101）表面的改性机理，并通过吸附能、态密度、差分电荷密度和Mulliken布居电荷分析界面之间的相互作用。

本研究从理论上分析改变氢氧化镁表面极性、提高氢氧化镁与有机材料相容性的改性机理，再结合实际生产，以期制备高分散、高附加值的氢氧化镁阻燃材料。我国盐湖资源丰富，富含大量镁资源，而作为盐湖资源开发废弃物之一的水氯镁石产量很高，大量堆积的水氯镁石不仅极大地浪费了资源，还会形成“镁害”，威胁当地生态环境。因此，利用盐湖废弃水氯镁石制备高性能氢氧化镁阻燃材料，不仅可以保护盐湖地区生态环境，还能增加经济效益，实现废弃水氯镁石的综合利用，对于盐湖资源的高效开发和高值化利用具有重要意义。

1 计算方法和模型建立

1.1 计算方法

本研究采用DFT 理论计算方法对氢氧化镁改性机理进行分析。DFT 被广泛用于描述分子之间的相互作用[11]。利用Materials Studio 软件Adsorption Locator 模块中的模拟退火法确定吸附的最佳位置，选择由原子集定义的表面区域，最大吸附距离设置为5 Å（1 Å = 0.1 nm）；几何优化则在DMol3模块中进行，采用PBE（Perdew-Burke-Ernzerhof）开发的广义梯度近似（Generalized Gradient Approximatio，GGA）泛函计算交换关联能，采用DFT Semi-core Pseudopots 处理原子的内层电子，基组为DNP 双重数值基组[12]。布里渊区k点网格设置为1 × 1 × 1，自洽收敛为2 × 10-5Ha，最大应力为0.005 Ha/Å，最大位移为0.005 Ha[13]。为了评估界面与分子间的相互作用和电荷转移情况，进行了吸附能、态密度、差分电荷密度和布居数的分析[12,14]。

1.2 模型建立

油酸、硅烷偶联剂G-570和MMA模型如图1所示：图1（a）为优化后的油酸分子，其分子式为C18H34O2；图1（b）为优化后的硅烷偶联剂G-570，其分子式为C10H20O5Si；图1（c）为优化后的MMA，其分子式为C5H8O2。

图1 结构优化后的吸附剂模型：（a） 油酸； （b） 硅烷偶联剂G-570； （c） MMAFig.1 Adsorbent model after structural optimization:（a） oleic acid； （b） silane coupling agent G-570； （c） MMA

氢氧化镁（101）面模型（图2（b））：氢氧化镁模型由Materials Pjoject 网站下载后优化得到（图2（a）），氢氧化镁（101）面模型如图2（b）所示，优化后晶格常数与文献[15]相比较吻合（表1），其误差低于2%，这表明优化后的结构模型与实际相符合,本计算方法以及所设置实验参数合理，结果可靠。由于氢氧化镁晶体在（101）面有微观内应变，具有表面极性大、表面能高、晶粒易团聚的特点[16]，是氢氧化镁改性研究的重点内容，因此构建了一层6 × 6 × 1 氢氧化镁（101）表面模型，真空层厚度设置为30 Å，以避免不必要的潜在相互作用力[17]；结构优化后，氢氧化镁（101）面的H-O-Mg 和O-Mg-O 键角分别为125.3°和110.1°、层间距离为2.335 Å，与文献[18]一致，计算结果合理。

表1 氢氧化镁晶格常数对比Table 1 Comparison of lattice constants of Mg（OH）2

图2 结构优化后（a） 氢氧化镁模型和（b）氢氧化镁 （101）面模型Fig.2 After structural optimization: （a） Mg（OH）2 model；（b） Mg（OH）2 （101） surface model

2 结果与讨论

2.1 稳定性分析

氢氧化镁表面吸附改性剂的稳定性可由吸附能进行判断，若吸附能为负值，则表明吸附为放热反应，吸附能绝对值越大，吸附越稳定[19]，计算公式如式（1）所示[20]：

式（1）中：Ea为吸附能；EMH/abs为氢氧化镁（101）表面吸附改性剂后的能量；EMH为氢氧化镁（101）表面的总能量；Eabs为表面改性剂的总能量。

通过吸附定位确定3种改性剂在氢氧化镁（101）表面的最佳吸附位点，即吸附能最低的模型（图3）；然后将所得模型进行几何优化，获得硅烷偶联剂G-570、MMA 和油酸吸附在氢氧化镁（101）面的最优结构（图4）；计算得到吸附3 种改性剂时的吸附能，结果见表2。

表2 氢氧化镁（101）表面吸附改性剂吸附能计算结果Table 2 Calculation results of adsorption energy of magnesium hydroxide （101） surface adsorption modifier单位：eV

图3 计算得到改性剂最佳吸附位点模型：（a） 吸附油酸； （b）吸附MMA； （c）吸附硅烷偶联剂G-570Fig.3 Optimal adsorption site model of the modifier: （a） adsorption oleic acid； （b）adsorption MMA；（c） adsorption silane coupling agent G-570

图4 优化后的氢氧化镁（101）表面吸附模型： （a） 油酸； （b） MMA； （c）硅烷偶联剂G-570Fig.4 Optimized surface adsorption model of magnesium hydroxide （101） : （a） oleic acid；（b） MMA； （c） silane coupling agent G-570

如图4（a）所示，吸附后油酸与氢氧化镁（101）表面相邻原子之间的距离为2.164、2.056 Å，硅烷偶联剂G-570、MMA与氢氧化镁（101）表面相邻原子之间的距离约为2.5 Å。由于相互作用的强度与距离成反比[21]，所以氢氧化镁（101）表面吸附油酸体系的结构稳定性更强。结合表2 数据分析，相比硅烷偶联剂G-570和MMA，氢氧化镁（101）表面吸附油酸体系的吸附能更低，稳定性更好。

2.2 态密度分析

为了进一步揭示吸附质和氢氧化镁之间的相互作用的本质，对吸附后的氢氧化镁进行了基于能带理论的分波态密度分析（Partial Density of States，PDOS）。若相邻原子之间的分态密度在相同能量处出现了尖峰，可用于直观说明吸附质和氢氧化镁（101）表面相邻原子之间的作用强弱[22]。

吸附质和氢氧化镁（101）表面相邻原子之间的PDOS分析如图5所示，即3种吸附体系中Mg原子和O 原子的杂化峰重叠情况。由图5 可见，在-25～5 eV之间，这3 种吸附体系的Mg 原子和O 原子的s、p 轨道均存在杂化峰，其中，吸附油酸体系的Mg317 和O323 之间的杂化峰数量最多，表明氢氧化镁吸附油酸的相互作用力最强，这也与前述实验结果一致。

图5 吸附质与氢氧化镁（101）表面相邻原子之间的PDOS图示:（a）和（b）吸附油酸；（c）吸附MMA;（d）吸附硅烷偶联剂G-570Fig.5 PDOS diagram between the adsorbent and adjacent atoms on the surface of magnesium hydroxide （101）：（a） and （b）adsorption oleic acid； （c）adsorption MMA； （d） adsorption silane coupling agent G-570

2.3 差分电荷密度分析

为了更加直观地表征改性剂和氢氧化镁（101）表面之间的相互作用力，对吸附后的电子密度差进行了分析，以揭示吸附过程中电子密度的变化，其计算公式如式（2）所示:

式（2）中：Δρ为吸附前后的电子密度差；ρabs+MH为吸附后的总电子密度；ρabs为油酸的电子密度；ρMH为氢氧化镁（101）表面的电子密度。

图6所示为氢氧化镁（101）表面3种吸附体系的差分电荷密度图，电荷减少和电荷密集的区域分别表示为蓝色和红色，等值面为0.05 e/Å。由图6（a）可见，油酸末端羧基O 原子周围电荷比较密集，蓝色区域主要在氢氧化镁（101）表面Mg 原子周围，并且油酸与氢氧化镁表面之间电子转移明显，说明在吸附过程中，电子从氢氧化镁（101）表面转移到油酸末端羧基，形成较强的Mg-O 相互作用；而MMA（图6（b））和硅烷偶联剂G-570（图6（c））与氢氧化镁表面之间的电子转移并不明显，说明并未形成较强的相互作用力。

图6 氢氧化镁（101）表面3种吸附体系的差分电荷密度分析图：（a）油酸；（b）MMA；（c）硅烷偶联剂G-570Fig.6 Differential charge density analysis diagram of three adsorption systems on magnesium hydroxide （101） surface:（a） oleic acid；（b） MMA； （c） silane coupling agent G-570

2.4 Mulliken布居电荷分析

Mulliken 布居电荷分析有助于更清晰地研究吸附质在氢氧化镁（101）表面的电子转移情况[23-24]，具体分析如表3所列。

表3 氢氧化镁（101）表面吸附体系布居电荷分析Table 3 Population charge analysis of magnesium hydroxide （101） surface adsorption system

由表3可知，3种吸附体系中氢氧化镁（101）表面的Mg原子均失去电子、吸附质的O原子均得到电子，并且与其他2种体系相比，吸附油酸体系的得失电子数最多、电荷转移数最多、相互作用力最大。因此，油酸可稳定吸附在氢氧化镁（101）表面，得益于油酸和氢氧化镁之间能够形成较强的Mg-O 相互作用，这也与上述态密度、差分电荷密度和吸附能的分析结果一致。

3 结 论

基于密度泛函理论（DFT）的理论计算，对常用改性剂油酸、硅烷偶联剂G-570 以及MMA 在氢氧化镁（101）表面的吸附情况进行了模拟，研究了改性剂对氢氧化镁表面的改性机理。吸附能、态密度、差分电荷密度和Mulliken 布居电荷等参数的分析结果表明，3种改性剂分子均可吸附在氢氧化镁（101）表面，以降低表面能，其中，油酸的吸附能最低，可在氢氧化镁（101）表面稳定吸附，并且吸附后油酸与氢氧化镁（101）表面之间发生了明显的电荷转移，形成Mg-O相互作用。改性剂通过这种相互作用稳定吸附在氢氧化镁表面，改变了氢氧化镁表面极性，提高氢氧化镁与高分子材料的相容性。本研究为提高氢氧化镁阻燃材料的性能研究提供了理论参考。