摘要:MgO是具有强极性的离子化合物, 电场诱导MgO吸附H₂是一种有效的储氢方法, 但外加的电场很强, 如何降低所需电场的强度是需要解决的关键问题. 本文在密度泛函理论水平上研究了电场中H₂在(MgO)₄团簇上的吸附性质. 结果表明(MgO)₄能承受强电场并保持立方结构, 可用于电场储氢. 电荷分析表明(MgO)₄在电场中被极化, 其偶极矩增大为场强0.005 a.u. 和0.010 a.u.时的1.67和3.33 Debye. H₂能稳定吸附在单个Mg/O原子上. H₂在Mg上为侧位吸附, 而在O上为端位吸附. 外加电场可提高其吸附强度. 仅需0.010 a. u.的外电场, 就可使H₂在Mg/O上的吸附能由无电场时–0.118/–0.060 eV提高到–0.225/–0.150 eV. 所需电场强度小于较大尺寸的 (MgO)₉团簇, 表明降低团簇尺寸是减少所需电场强度的一种可能方式. 利用QTAIM方法研究了H₂与(MgO)₄间的弱相互作用, 表明电场使团簇及氢分子极化, 从而增强了其间的静电作用. 当团簇尺寸降低时, 更多的原子位于表面, 且具有较低的配位数, 更容易被极化, 因此储氢所需的电场强度更低. 电场中(MgO)₄中最多能吸附16个H₂, 相应的质量密度为16.7 wt%, 表明(MgO)₄是一种可能的电场储氢材料.
关键词: 电场/
(MgO)₄/
储氢/
电子结构

English Abstract

全文HTML

--> --> -->

氢储量丰富、燃烧值高、可再生、清洁无污染, 是化石燃料的理想替代品, 而氢的储运是制约氢能利用的主要因素^[1]. 低温液态(21 K)和高压(35—70 MPa)储氢的成本高, 还存在严重的安全问题^[2]. 研发各种固体储氢材料是解决储氢问题的基本途径^[2,3]. 根据H₂与主体材料的结合强度不同, 固体储氢可粗略分为物理吸附和化学吸附. 物理吸附材料与氢作用较弱(吸附能一般小于–0.1 eV), 需在较低的温度下储存, 且储氢密度较小. 化学吸附材料虽达到很高的储氢密度, 但由于吸附材料与氢较强的化学作用(吸附能一般大于–1 eV), 脱氢较困难, 常需在几百度的高温下进行^[4]. 对理想的储氢材料, 追求的目标是对氢分子产生介于物理吸附和化学吸附之间的结合方式. 根据热力学分析, 在室温和1 MPa 条件下, 理想固体储氢材料吸附单个氢分子的吸附能应为–0.1—–0.2 eV^[5]. 探寻具有适宜吸附强度的储氢方法对寻找室温下储氢材料具有重要意义.
近年来, 电场诱导极性材料储氢逐渐受到人们的关注^[6—11]. 电场极化增强了主体材料与H₂的静电作用, 从而提高对H₂的吸附强度; 而将外电场移除, 又可实现快速脱氢. 通过对电场诱导极性材料(单层BCN, BC₂N, BN, AlN和 C纳米材料等)储氢的研究^{[6—8,10,11]}, 表明电场能显著改善其储氢性能. 然而, 所需外加的电场强度约为0.04—0.05 a.u. (1 a.u. = 5.1 × 10¹⁰ V/m), 如此强的电场, 目前还无法实现. 因此, 如何降低所需的电场强度是电场储氢需要解决的关键问题.
采用极性更强的材料是降低场强的一种可能方式. MgO是具有较强极性的离子化合物. 虽然H₂在MgO晶体表面的吸附很早就受到人们的关注, 而理论与实验研究表明H₂通常在MgO晶体表面形成弱的物理吸附^[12—16]. 纳米材料如团簇因其极大的比表面积和独特的电子结构而具有比相应体相材料更优越的吸附性能. 虽然团簇储氢面临诸多挑战, 如实验研究比较困难; 如何能在温和条件下操作的储氢以及寻找具有高存储密度和优良循环动力学性能的新储氢体系等^[2]. 通过理论研究能在原子、分子水平上探究储氢的微观机制, 从而设计合理的储氢方法, 并通过调控团簇的尺寸及组成, 为寻找适宜的储氢体系及相关的实验研究提供理论指导^[3,17]. Kwapien等^[18]利用从头计算方法研究了H₂在(MgO)_n (n = 6, 8)上的吸附及解离. Chen 等^[19]也在B3LYP/6-31G水平上研究了H₂在 (MgO)₉和(MgO)₁₂幻数团簇上的吸附性质. 虽然 (MgO)_n团簇的储氢能力较MgO晶体有所改善, 但其吸附能仅为–0.03—–0.08 eV, 而将MgO材料置于外电场中, 可进一步提高其储氢性质. 实验上Sun等制备了多孔MgO材料, 发现外加电场后其储氢能力得到一定程度的提高^[9]. 我们之前的理论工作也证实电场能显著改善(MgO)₉储氢性质, 且所需电场强度仅为0.025 a.u.^[20]. 鉴于此, 我们关心若减小团簇尺寸, 能否进一步降低所需外场强度？具有立方结构的(MgO)₄是幻数团簇, 特别稳定, 同时也是MgO岩盐结构的基本组成单元^[21—23], 因此本文研究了电场中(MgO)₄的储氢性质, 结果发现电场能显著改善其储氢性质, 重要的是所需电场小于大尺寸的(MgO)₉团簇, 仅需外加0.010 a.u.的电场就可使其对H₂的吸附能提高到–0.225 eV, 表明减小团簇尺寸是降低电场强度的可能方式. 电场中(MgO)₄中最多能吸附16个H₂, 相应的质量密度为16.7 wt%, 表明(MgO)₄是一种可能的电场储氢材料.

本文计算在密度泛函理论B3LYP 框架下进行, 电子交换能采用HF和Becke三参数混合形式, 电子相关能的定域部分使用函数VWNⅢ, 而非定域部分使用LYP 的泛函形式^[24,25]. 基组采用CC-PVTZ^[26]. 电场作用通过在哈密顿中引入$ {{\mu}} \cdot {{F}}$项考虑, 其中$ {{\mu}}$为分子偶极矩, $ {{F}}$为外场矢量. 我们优化了场强分别为0, 0.005和0.010 a.u. 时(MgO)₄的构型及一个和多个H₂在(MgO)₄上的稳定吸附结构. 结构优化中力和位移的收敛标准分别是 0.00045 Hartree/?和0.0018 ?. 进行了频率计算以保证优化后的结构没有虚频, 即对应于势能面上的一个局域极小点. 用平均吸附能E_a衡量对H₂的吸附强度, 其定义为

${E_{\rm{a}}} = \left\{ {E\left[ {{{\left( {{\rm{MgO}}} \right)}_4} \cdot n{{\rm{H}}_2}} \right] - E\left[ {{{\left( {{\rm{MgO}}} \right)}_4}} \right] - nE\left[ {{{\rm{H}}_2}} \right]} \right\}/n$

式中$E\left[ {{{\left( {{\rm{MgO}}} \right)}_4} \cdot n{{\rm{H}}_2}} \right]$、$E\left[ {{{\left( {{\rm{MgO}}} \right)}_4}} \right]$及E[H₂]分别为吸附n个H₂的复合体系、(MgO)₄团簇及H₂的总能量. 所有结构的总能量均进行了零点能(ZPE)校正, 并对计算的吸附能进行了基组重叠误差(BSSE)校正. 为研究电场中(MgO)₄吸附体系的稳定性, 我们还进行了Borne-Oppenheimer分子动力学模拟, 模拟温度为300 K, 温度通过Nosé-Hoover方法调节, 动力学时间步长为1 fs, 模拟总时间为3 ps. 对得到的稳定吸附结构采用QTAIM (quantum theory of atoms in molecules)方法进行电子结构分析, 研究了电场中(MgO)₄与H₂的相互作用. 结构优化使用Gaussian09程序完成^[27], 分子动力学模拟采用Dmol³进行^[28], 电子结构分析则采用Multiwfn软件^[29].

2

-->

3.1.电场中(MgO)₄的结构

优化得到无外场时(MgO)₄的最稳定结构如图1所示. 该结构为具有T_d对称性的立方结构, 因此所有Mg—O键长均为1.96 ?, 与之前文献报道的结果一致^[21—23]. NPA电荷布居分析表明(MgO)₄中Mg将其3s电子向O 2p轨道转移, 从而形成了Mg阳离子和O阴离子(图1左侧结构图中蓝色代表失电子的原子, 即Mg原子, 而红色小球则代表得电子的O原子). 由于(MgO)₄的高对称性, 每个Mg原子转移的电荷均相同, 为1.493 e, 因此Mg与O通过强的静电相互作用结合.
图 1 不同场强下(MgO)₄的稳定结构(原子颜色与电荷得失相关, 由蓝色至红色表示失电子越多到得电子越多)
Figure1. The stable structures of (MgO)₄ under the electric fields with different intensities(the color is correlated with the gaining or losing of electrons, from blue to red, it represents the variation of from the losing to obtaining electrons).

我们进一步优化了强度分别为0.005 a.u. 和0.010 a.u. 的不同方向的外电场中(MgO)₄的稳定结构, 结果表明当外场方向沿O₁-Mg₆体对角线方向时, 团簇的能量最低, 其稳定结构如图1所示(场强为0.005 a.u. 时 (MgO)₄的结构与场强为0.010 a.u. 时相似, 图中未给出). 一方面电场中(MgO)₄仍保持立方结构, 表明团簇能承受外电场, 可用于电场储氢. 另一方面, 由于Mg为阳离子而O为阴离子, 外加电场后Mg顺电场方向移动而O逆电场方向移动, 因此(MgO)₄的几何结构有一定畸变. 表1给出了不同强度的电场中(MgO)₄的Mg/O原子NPA电荷及Mg—O键长. 由于电子逆着电场方向运动, 导致沿场强方向(即O₁-Mg₆方向)电荷重新布居. 电场中Mg₆, O₃, O₅和 O₇(图1右侧蓝色或篮白色原子)均失电子, 而O₁, Mg₂, Mg₄和Mg₈(图1右侧红色或淡粉原子)得电子. 因此与无电场时相比, Mg和O原子均可分为两类: 一类是电荷绝对值增加的Mg₆/O₁; 另一类是电荷绝对值减少的Mg₂, Mg₄, Mg₈和O₃, O₅, O₇. 第一类Mg₆/O₁原子的电荷由无电场时的1.493 e/–1.493 e增加为场强为0.010 a.u. 时的1.554 e/–1.512 e, 第二类的Mg/O原子的电荷则减少为1.460 e/–1.474 e. 相应地, Mg—O键也可以分为两类: 第一类以Mg₆/O₁为端点的Mg—O键长R_I由无电场时的1.96 ?缩短为场强0.010 a.u. 时的1.93 ?; 而其他Mg—O键长(即第二类R_II)则拉长至1.99 ?. 电场诱导电荷重新分布, 使(MgO)₄极化, 随着电场逐渐增强, 其偶极矩由无电场时的0 Debye 分别增加为场强为0.005 a.u. 和0.010 a.u. 时的1.67和3.33 Debye.

F/a.u.	QMgI/QMgII	QOI/QOII	RI	RII
0	1.493/1.493	–1.493/–1.493	1.96	1.96
0.005	1.476/1.527	–1.484/–1.503	1.94	1.97
0.010	1.460/1.554	–1.474/–1.512	1.93	1.99

表1不同场强下(MgO)₄中两类Mg/O原子的NPA 电荷(e)及Mg—O键长(?) (Q_MgI/Q_MgII和Q_OI/Q_OII分别是Mg₂, Mg₄, Mg₈/Mg₆和O₃, O₅, O₇/O₁上的电荷; Mg—O键长R_I = R₁₂ = R₁₄ = R₁₈ = R₆₃ = R₆₅ = R₆₇, 而R_II = R₂₃ = R₂₅ = R₄₃ = R₄₇ = R₈₅ = R₈₇)
Table1.The NPA charges for the two types of Mg/O atoms and the Mg—O distances (Q_MgI/Q_MgII and Q_OI/Q_OII are the charges of Mg₂, Mg₄, Mg₈/Mg₆ and O₃, O₅, O₇/O₁, respectively. The Mg—O distances R_I = R₁₂ = R₁₄ = R₁₈ = R₆₃ = R₆₅ = R₆₇, while R_II = R₂₃ = R₂₅ = R₄₃= R₄₇ = R₈₅ = R₈₇).

2 -->

3.2.电场中单个H₂在(MgO)₄上的吸附

为了确定无外场时H₂在 (MgO)₄上的稳定吸附结构, 我们优化了H₂在 (MgO)₄上的各种可能吸附方式, 包括H₂分别吸附在单个Mg/O原子上方, Mg—O键上方及立方体表面上等. 结果表明H₂只能稳定地吸附在(MgO)₄的单个Mg/O原子上, 与(MgO)₉的结果相似^[20]. 由于(MgO)₄的高对称性, 其表面所有Mg和O原子分别等价, 因此只有两种稳定吸附结构如图2所示. H₂在Mg上为侧位吸附(H₂分子的方向与H₂分子到 Mg离子的连线垂直), 而在O上为端位吸附(H₂分子的方向与H₂分子到 O离子的连线一致).
图 2 不同场强条件下H₂在(MgO)₄上的稳定吸附结构 (绿色球Mg原子; 红色球O原子; 灰色球H原子)
Figure2. The stable structures of H₂ adsorbed on (MgO)₄ under the electric fields with different intensities (the green, red and gray balls are the Mg, O and H atoms, respectively)

电场中 (MgO)₄结构的对称性降低, 所有Mg/O原子不再等价. 如前所述, Mg和O原子均分为两类. 一类是电荷绝对值增加的Mg₆/O₁, 另一类是电荷绝对值减少的Mg₂, Mg₄, Mg₈和O₃, O₅, O₇. 由于外加电场相对较弱, 当H₂分子吸附在电荷上述六个绝对值减少的原子上时, 仍能形成吸附, 但其吸附能仅为–0.03—–0.06 eV(0.010 a.u.), 低于无电场时的吸附能; 而在Mg₁及O₆上时的吸附能高于无电场时的吸附能, 因此这里我们只讨论H₂在这两个原子上的吸附, 相应的吸附结构和吸附能、结构参数及NPA电荷数据分别如图2和表2所示. 图2给出了场强为0.010 a.u.时, H₂在即Mg₆/O₁上的稳定吸附结构(由于场强为0.005 a.u. 时的吸附结构与场强为0.010 a.u. 时的稳定结构相似, 图中未给出). 由于外加电场相对较弱, 与无电场时相似, H₂在Mg和O上的吸附结构仍为侧位吸附和端位吸附. 分子动力学模拟结果表明中这些吸附结构在300 K条件下仍能稳定存在. 以0.010 a.u.时H₂在Mg原子上的吸附结构为例, 其总能量随时间的演化如图3所示, 在整个时间范围内, 总能量以–31458.12 eV为中心波动, 且H₂一直稳定吸附在(MgO)₄上, 表明电场中的吸附体系是稳定的.

Site	F/a. u.	Ea/eV	RH—H/?	RH—Mg/O/?	QH	QH2
H2 on Mg	0	–0.118	0.751	2.217	0.023/0.023	0.046
	0.005	–0.172	0.752	2.179	0.030/0.030	0.060
	0.010	–0.225	0.753	2.093	0.044/0.044	0.088
H2 on O	0	–0.060	0.750	2.365	0.044/–0.068	–0.024
	0.005	–0.101	0.755	2.248	0.069/–0.110	–0.041
	0.010	–0.150	0.763	2.136	0.094/–0.157	–0.063

表2H₂在(MgO)₄上的吸附能、到团簇距离R_H-Mg/O, H—H键长R_H—H及H和H₂的NPA 电荷
Table2.The adsorption energies E_a, H—H bond lengths R_H—H, distances between H₂ and cluster R_H2-Mg/O and NPA charges of H atoms and H₂ for (MgO)₄H₂

图 3 场强为0.010 a.u. H₂吸附在Mg上时总能量随时间的演化(附图中为演化始末的吸附结构)
Figure3. The variation in the energy for the structure of H₂ adsorbed on Mg atom under the electric field with the intensity of 0.010 a. u.(the adsorbed structures are also presented in attached map)

结合图1及图2可见, 吸附H₂后, (MgO)₄团簇的结构未见显著改变, 即吸附过程未显著改变团簇的振动模式, 因此ZPE校正对E_a的影响较小(最大值仅为0.02 eV), 而BSSE的校正则更小(最大值仅为0.01 eV). 表2列出了ZPE及BSSE校正后的平均吸附能E_a, 随着场强增加至0.010 a.u., H₂在Mg/O上吸附能均逐渐增大, 由无电场时的–0.118/–0.060 eV 增加到–0.225/–0.150 eV, 达到了理想的储氢吸附强度. 同时H₂到Mg/O的距离变短, 由无电场时的2.217/2.365 ? 减少至2.093/2.136 ?. 另一方面, 电场中H—H键长变长, 由无电场时的0.751/0.750 ? 增加至0.753/0.763 ?. 上述结果表明电场使H₂与团簇的相互作用增强而H—H键被削弱. 同时, 相同场强条件下, H₂在Mg的吸附能总是大于在O上, 表明H₂在Mg上吸附更稳定.
2 -->

3.3.电场中(MgO)₄与H₂的弱相互作用

为深入理解电场中 (MgO)₄与H₂的相互作用, 我们采用QTAIM的拓扑分析法对其电子结构进行了研究. 拓扑分析中, 电子密度梯度为零的位置称为临界点(CP). 根据电子密度Hessian矩阵负本征值的个数, 临界点可分为四种. 若三个本征值均为负值, 称为核临界点(NCP); 若有两个负本征值, 该临界点常位于两个成键原子之间, 称为键临界点(BCP); 若只有一个负本征值, 该临界点常处于环的中间, 称为环临界点(RCP); 若没有负的本征值, 称为笼临界点(CCP). 拓扑分析得到电场中H₂在Mg/O上吸附结构的临界点如图4所示.
图 4 F = 0.010 a.u. 电场中H₂在Mg/O上稳定吸附结构的临界点 (蓝色球为核临界点, 绿色球为键临界点; 红色球为环临界点)
Figure4. The critical points for the stable adsorption structures of H₂ on Mg/O under the electric field with F = 0.010 a.u.(the blue, green and red balls represent the nuclear, bond and ring critical points, respectively)

由图4可见, 除了Mg/O原子的NCP和各表面的RCP外, 还有对应于Mg—O键、H—H键的BCP, 特别是出现了Mg—H键和O—H键的BCP, 表明吸附的H₂与吸附位的Mg/O原子间存在相互作用. 表3列出了相应BCP的拓扑参数, 分析这些参数可以了解相应BCP的成键性质. 负的拉普拉斯电荷密度(${\nabla ^2}\rho $)和总能密度(H(r))表示共价相互作用, 而正值则表示闭壳层电子间相互作用, 如离子键、范德瓦耳斯作用、氢键和互斥作用等. 电子定域化函数(ELF)也是研究共价键的有力工具, 其值接近1表明强的共价键, 而接近0则对应非共价作用. 由表3数据可见, 不同场强条件下, H—H键的${\nabla ^2}\rho $和H(r)均为负值, 而电子密度ρ较大且ELF接近1, 表明H—H形成强的共价键. 而Mg—O键的${\nabla ^2}\rho $和H(r)均为正值, 而ρ和ELF较小, 表明Mg—O形成闭壳层相互作用. 同时较大的${\nabla ^2}\rho $和H(r)表明Mg—O间形成离子键. H₂与Mg/O原子间${\nabla ^2}\rho $和H(r)也均为正值, 同时ρ和ELF接近0, 表明H₂与Mg/O间也形成闭壳层相互作用, 但较小的${\nabla ^2}\rho $和H(r)表明其间的相互作用较弱, 属于弱相互作用. 另外随着电场的逐渐增强, H₂与Mg/O间ρ, ${\nabla ^2}\rho $, H(r) 和ELF绝对值均增大, 表明电场增强了Mg—H和O—H相互作用, 同时H—H间的ρ和${\nabla ^2}\rho $绝对值减小, 表明H—H键被削弱, 与之前能量分析的结果一致.

F/a. u.	H2 on Mg						H2 on O
	BCP	ρ	${\nabla ^2}\rho $	H(r)	ELF		BCP	ρ	${\nabla ^2}\rho $	H(r)	ELF
0.000	Mg—H	0.013	0.059	0.002	0.029		O—H	0.011	0.038	0.002	0.045
	H—H	0.266	–1.167	–0.292	1.000		H—H	0.265	–1.167	–0.292	1.000
	Mg—O	0.056	0.403	0.005	0.057		Mg—O	0.056	0.403	0.005	0.057
0.005	Mg—H	0.014	0.063	0.002	0.031		O—H	0.015	0.049	0.002	0.060
	H—H	0.260	–1.165	–0.292	1.000		H—H	0.261	–1.142	–0.287	1.000
	Mg—O	0.054—0.059	0.382—0.426	0.004—0.005	0.056—0.059		Mg—O	0.055—0.059	0.390—0.427	0.004—0.005	0.057—0.059
0.010	Mg—H	0.016	0.069	0.002	0.034		O—H	0.019	0.063	0.002	0.080
	H—H	0.255	–1.162	–0.291	0.999		H—H	0.255	–1.102	–0.278	0.999
	Mg—O	0.051—0.061	0.341—0.533	0.004—0.005	0.055—0.060		Mg—O	0.025—0.061	0.358—0.444	0.004—0.005	0.054—0.060

表3H₂在(MgO)₄上吸附结构的拓扑参数
Table3.The topological parameters for the adsorption structures of H₂ on (MgO)₄

约化密度梯度(RDG(r))是研究弱相互作用的有力工具, 其定义为

$RDG\left( r \right) = \frac{1}{{2{{\left( {3{{\text{π}}^2}} \right)}^{1/3}}}}\frac{{\left| {\nabla \rho \left( r \right)} \right|}}{{\rho {{\left( r \right)}^{4/3}}}}.$

对一个体系, 当某些位置的ρ(r)和RDG(r)均较小时, 对应于弱相互作用区域. 电场中(MgO)₄吸附H₂体系的RDG(r)如图5所示. H₂无论吸附在Mg还是O上, 图中均出现两个紫色的弱相互作用区域, 一个是位于(MgO)₄各表面的区域, 另一个则是位于H₂与Mg/O原子中间的位置, 表明H₂与Mg/O间确实存在弱相互作用, 这与拓扑分析的结果一致.
图 5 F = 0.010 a. u. 电场中H₂在Mg/O离子上吸附结构的RDG等能面
Figure5. The RDG isosurface for the adsorption structures of H₂ on Mg/O under the electric field with F = 0.010 a. u.

如何知道这两个弱相互作用的类型？根据QTAIM理论, RCP电子密度Hessian矩阵的第二个本征值(λ₂)为正值, 而BCP则为负值. 同时弱相互作用强度与ρ(r)相关, 如范德瓦耳斯作用具有很小的ρ(r), 而位阻效应等互斥作用则具有相对较大的ρ(r). 因此λ₂的正负号(记为[λ₂])与ρ(r)的乘积(记为[λ₂]ρ(r))可作为确定弱相互作用类型的指标函数. 利用[λ₂]ρ(r)结合RDG(r), 不仅可以定位体系中弱相互作用的位置, 同时也能确定弱相互作用的类型. 利用[λ₂]ρ(r)值对RDG(r)等能面进行着色, 形成填色图, 只需考察弱相互作用区域的颜色, 即可确定其类型. 蓝色表示氢键、卤键等具有较强吸引的相互作用; 绿色表示范德瓦耳斯吸引作用; 而红色则表示位阻效应对应的互斥作用. 电场中(MgO)₄吸附H₂体系的RDG(r)等能面填色图如图6所示. 由图可见, 位于(MgO)₄表面的区域为红色, 表明其间是互斥作用, 而H₂与Mg/O原子间则是绿色的弱相互作用区域, 表明其间是范德瓦耳斯吸引作用.
图 6 F = 0.010 a. u. 电场中H₂在Mg/O离子上吸附结构的填色RDG等能面
Figure6. The color-filled map of RDG isosurface for the adsorption structures of H₂ on Mg/O under the electric field with F = 0.010 a. u.

范德瓦耳斯作用本质是静电作用, 因此我们对(MgO)₄吸附H₂体系的NPA电荷进行了计算并列于表2. H₂吸附在Mg上, 两个H原子均带正电荷, H₂上的净电荷为正, 表明电荷由H₂分子向团簇转移. 当H₂吸附在O上时, 近O离子端的H原子带正电荷, 另一端H原子带负电荷, 表明H₂被极化. 同时H₂上的净电荷为负, 表明电荷由团簇向H₂分子转移. 随着电场的增强, H₂吸附在Mg/O上时的电荷绝对值均单调增加, 表明电场促进了H₂与团簇间的电荷转移. 另一方面, 我们还计算了电场中(F = 0.010 a.u.)与无电场时H₂吸附在Mg/O上的电荷密度差分, 如图7所示. 图中的暖色/冷色和实线/虚线分别表示电子增加和减小的区域. 当H₂吸附在Mg上时, Mg与H₂间电子密度增加, 而相反的另一侧电子密度减小, 表明电场中H₂被极化; 当H₂吸附在O上时, 极化效应更加明显, O与H₂间电子密度减少, 而相反的另一侧电子密度增加. 这表明电场中氢分子被有效极化. 综合上述分析, 一方面电场促进了H₂与(MgO)₄间的电荷转移, 另一方面电场使团簇及H₂被有效极化, 增强了其间的静电作用, 从而提高了(MgO)₄对H₂的吸附. 另外, 在相同的场强条件下, 当H₂吸附在Mg上时, 转移的电荷量总是大于吸附在O上, 因此H₂吸附在Mg上时更稳定.
图 7 H₂在Mg (a)和O (b)原子上吸附结构的电荷密度差分
Figure7. The color-filled contours of charge density difference for adsorption structures of H₂ on Mg (a) and O (b), respectively

2 -->

3.4.电场中多个H₂在(MgO)₄上的吸附

为研究电场中(MgO)₄的储氢质量密度, 我们优化了场强为0.010 a.u.时多个H₂在(MgO)₄上的吸附结构. 由于H₂只能吸附在单个Mg/O原子上方, 因此首先优化了多个H₂吸附在(MgO)₄中单个Mg/O上的结构, 发现每个Mg/O原子上最多能吸附2个H₂分子. 因此我们又优化了16个H₂吸附在(MgO)₄上的结构(如图8所示), 当继续增加H₂分子数目时, 则不能形成稳定的吸附结构, 表明电场中(MgO)₄最多能吸附16个H₂, 相应的质量密度为16.7 wt%, 达到了美国能源部制定的质量密度需达到6 wt%的要求, 表明(MgO)₄是一种优良的电场储氢材料. 吸附的16个H₂的键长为0.747—0.752 ?, 到吸附的Mg/O原子间的距离为2.316—2.958 ?, 由于先吸附的H₂对后吸附H₂的位阻效应, 相应每个H₂分子的平均吸附能也减小为–0.09 eV.
图 8 F = 0.010 a.u. (MgO)₄吸附16个H₂的稳定结构(绿色球为Mg原子, 红色球为O原子; 灰色球为H原子)
Figure8. The adsorption structures of 16 H₂ on (MgO)₄ under the electrical field with the intensity of 0.010 a. u. (the green, red and grey balls are Mg, O and H atoms, respectively).

电场诱导极性材料储氢是一种具有潜力的储氢方法, 如何降低所需的电场强度是需要解决的关键问题. Zhou等^[8]提出采用富含阴离子的强极性材料是一种可能的方式. 对电场中(MgO)₄储氢行为的研究表明, 只需外加较弱的电场就可显著提升其对H₂的吸附性能, 所加外电场的强度小于较大尺寸的(MgO)₉团簇, 表明降低团簇尺寸也是减少外电场强度的一种可能方式. 如获得同样的吸附能(–0.23 eV), 对(MgO)₄仅需外加0.010 a.u.的电场, 而对(MgO)₉, 外加电场需达到0.250 a.u.^[19]. 为进一步证实该结论, 我们在同样的理论水平上, 对尺寸更小的(MgO)₂在电场中吸附H₂的结构进行了优化, 结果表明达到同样的吸附能, 所需电场确实更低, 仅为0.005 a.u.. 电子结构分析表明电场使(MgO)₄及H₂被有效极化, 即极化效应是提高(MgO)₄对H₂的吸附性能的重要原因. 当团簇尺寸降低时, 更多的原子位于表面, 而这些位于表面的原子配位数较低((MgO)₉, (MgO)₄和(MgO)₂中Mg/O原子的最高配位数分别为5, 3和2), 更容易被极化, 因此仅需较弱的电场就可显著改善其储氢性质.
将小尺寸 (MgO)_n幻数团簇担载于轻质基底材料上形成复合体系, 可能进一步提升其储氢性能, 并具有更复杂、更丰富的物理过程. 另外一些稳定的幻数团簇具有原子特性, 称为超原子. 以这些超原子为单元组装形成的新型复合材料, 仍能保持超原子单元本身优良的物化性能^[30]. 因此以小尺寸 (MgO)_n幻数团簇为单元形成的复合物或组装材料, 仍可能具有较好的电场储氢性能, 而其中涉及的尺寸和电场效应、存储密度及储氢机理等也是值得进一步研究的课题.

本文在密度泛函理论水平上研究了电场中H₂在(MgO)₄上的吸附行为. 结果表明电场能提高H₂在(MgO)₄部分位置上的吸附强度. 仅需外加0.010 a.u. 的电场, 就可使H₂在Mg/O上的吸附能由无电场时–0.118/–0.060 eV分别提高到–0.225/–0.150 eV. 所加外电场的强度小于较大尺寸的(MgO)₉团簇, 因此降低团簇尺寸是降低外电场强度的一种可能方式. 电子结构分析表明, 电场中H₂与(MgO)₄形成的弱相互作用为范德瓦耳斯静电吸引. 一方面电场促进了H₂与(MgO)₄间的电荷转移, 另一方面电场使团簇及H₂被有效极化, 增强了其间的静电作用, 从而提高了(MgO)₄对H₂的吸附. 电场中(MgO)₄最多能吸附16个H₂, 相应的质量密度为16.7 wt%, 表明(MgO)₄可作为电场储氢材料.
感谢国家超级计算机中心在广州和深圳提供的计算资源.