PaperTan: 写论文从未如此简单

材料工程学

一键写论文

基于第一性原理计算的Ti3C2 MXene表面缺陷对储氢性能的影响机制研究

作者:佚名 时间:2026-01-30

本研究采用第一性原理计算方法,系统探究Ti₃C₂ MXene表面缺陷对储氢性能的影响机制。构建含单/双空位、取代缺陷的模型,通过VASP软件结合PBE泛函与色散校正,分析吸附能、电荷转移等指标。结果表明,表面缺陷可改变电子结构、增加活性位点,Ti空位使吸附能从-0.15 eV提升至-0.35 eV(接近理想范围),含氧官能团优化吸附位点分布。研究揭示缺陷工程增强储氢性能的微观机制,为设计高效MXene基储氢材料提供理论支撑,助力氢能商业化应用。

第一章引言

全球能源需求持续不断地增长,与此同时环境问题也变得越来越严重。在这样的情况下,开发清洁且高效的能源载体成为了当前科研工作里的一个重要方向。

氢能属于一种具有高能量密度、零污染排放特点的新型能源,被当作替代传统化石燃料的理想选择。不过,氢能若要实现大规模商业化应用,还存在储氢技术方面的瓶颈。传统的储氢方法,例如高压气态储氢、低温液态储氢,存在安全性不足、能耗过高等一系列问题。于是,开发新型高效的储氢材料就成了突破这一储氢技术瓶颈的关键所在。

在各种各样的储氢材料当中,二维材料依靠其独特的物理化学性质以及较大的比表面积,展现出了良好的储氢潜力。Ti3C2 MXene是典型的二维过渡金属碳化物,它具有优异的导电性、良好的化学稳定性,并且表面化学性质还能够进行调控,所以在储氢领域受到了比较多的关注。研究结果显示,Ti3C2 MXene表面存在着很多活性位点,而且表面的缺陷结构会明显地对氢分子的吸附行为以及储氢性能产生影响。

利用第一性原理计算方法,能够从原子尺度深入地去研究Ti3C2 MXene表面缺陷与氢分子之间的相互作用机制,从而为设计高性能储氢材料提供理论方面的支持。第一性原理计算是以量子力学作为基础的,它通过求解薛定谔方程来获取材料体系的电子结构和能量信息,具备高精度和可预测的特性。这种计算方法能够准确地描述材料表面的电子态分布、吸附能变化以及反应路径,进而揭示出储氢过程的微观机制。

在研究Ti3C2 MXene储氢性能的时候,引入表面缺陷能够改变材料表面的电荷分布以及活性位点的数量,从而增强氢分子的吸附强度。这种调控机制对于提高储氢容量、改善吸放氢动力学性能是非常重要的。本研究构建了不同缺陷类型的Ti3C2 MXene模型,采用第一性原理计算方法系统地分析表面缺陷对储氢性能产生影响的规律,为优化MXene基储氢材料的设计和制备提供理论依据,既具有学术方面的价值,也具有实际应用的前景。

第二章研究方法与模型构建

2.1Ti3C2 MXene表面缺陷模型的建立

图1 Ti3C2 MXene表面缺陷模型建立流程

Ti₃C₂ MXene是典型的二维过渡金属碳化物,其晶体结构有独特的层状特点。这种材料是通过从Ti₃C₂层刻蚀掉Al原子制得的,制得后表面一般会附着像 -O、 -OH 或者 -F 这样的官能团。原始的Ti₃C₂ MXene属于六方晶系,它的空间群是P6₃/mmc,晶格参数方面,a大约是3.06 Å,c大约是18.6 Å。Ti原子和C原子交替排列起来,构成了准二维的层状结构。

在构建表面缺陷模型的时候,要先建立pristine Ti₃C₂的优化结构,让几何构型和实验或者理论报道的值保持一致,以此为后续的缺陷研究提供可靠的基础。本研究是基于密度泛函理论的,系统地构建了多种表面缺陷模型,构建这些模型的目的是探究这些缺陷对储氢性能会产生怎样的影响。单空位缺陷存在几种类型,有Ti空位、C空位以及O官能团空位,这些单空位缺陷是通过在优化后的超胞里把对应原子移除来实现的。就拿构建Ti空位来说,要选择表层的Ti原子将其完全删除,删除之后进行几何弛豫从而得到稳定的结构。

双空位缺陷比单空位缺陷更复杂,构建双空位缺陷需要考虑不同原子位置的组合情况,就像Ti - C双空位这种,还需要通过计算形成能来确定最稳定的缺陷结构。取代缺陷模型是把特定的原子替换成目标原子,例如用N原子去替换表面的O官能团,这样就形成了N - O取代结构。所有的缺陷模型都是使用4×4×1的超胞体系来构建的,并且在垂直片层的方向要设置20 Å的真空层,设置真空层是为了消除周期性边界条件所引起的层间相互作用。

表1 Ti3C2 MXene表面缺陷模型参数设置
缺陷类型缺陷结构描述缺陷浓度(%)真空层厚度(Å)计算超胞尺寸(ų)缺陷形成能计算方法
单Ti空位(V_Ti)移除表面单层Ti原子6.252015.12×15.12×25.00E_f = E_def - E_perf + μ_Ti
双Ti空位(V_Ti2)移除相邻两个表面Ti原子12.502015.12×15.12×25.00E_f = E_def - E_perf + 2μ_Ti
C空位(V_C)移除表面单层C原子6.252015.12×15.12×25.00E_f = E_def - E_perf + μ_C
Ti-C双空位(V_TiC)移除相邻表面Ti和C原子12.502015.12×15.12×25.00E_f = E_def - E_perf + μ_Ti + μ_C
表面Ti原子替换为O原子(Ti→O)表面Ti原子被O原子取代6.252015.12×15.12×25.00E_f = E_def - E_perf - μ_Ti + μ_O

在对模型进行优化时,要使用基于广义梯度近似(GGA)的PBE泛函来描述交换关联作用,平面波截断能要设置为500 eV,k点网格要选择3×3×1,这样做能够平衡计算精度和效率。把优化后的晶格参数和文献报道的值进行对比(误差要小于2%),通过这样的对比验证了原始结构和缺陷模型的合理性,为后续的储氢性能模拟计算奠定了可靠的基础。这种系统的模型构建方法不仅能够准确地表征不同缺陷类型对Ti₃C₂ MXene电子结构和表面化学性质所产生的影响,而且还为深入理解其储氢机制提供了理论方面的支撑。

2.2第一性原理计算参数设置

在第一性原理计算里,参数设置很重要,它关乎模拟结果的准确性与可靠性。本研究用基于密度泛函理论(DFT)的VASP软件做计算,VASP主要是求解Kohn - Sham方程,通过求解这个方程来描述体系的电子结构。对于Ti₃C₂ MXene表面的储氢体系,选择广义梯度近似(GGA)下的PBE泛函当作交换关联泛函,并且结合PBE - D3色散校正方法,使用这种结合方式能够更准确地描述氢分子和表面之间存在的范德华相互作用。为了同时考虑计算精度和效率,选用投影缀加波(PAW)赝势来描述原子核与电子的相互作用,此方法会把全电子波函数分解成平滑的赝波函数以及缀加原子轨道区域,这样做可以大幅降低平面波基组的计算需求。

平面波截断能的设置会直接对基组完备性产生影响。本研究选择500 eV作为截断能阈值,这个500 eV的数值是通过收敛性测试确定下来的,确定这个数值是为了保证体系总能量误差小于1 meV/atom。在进行k点采样的时候,使用Monkhorst - Pack网格对布里渊区进行积分,针对构建的1×1×1表面超胞模型,把k点密度设置为1×1×1,将k点密度设置成这样既可以满足表面计算的特殊需求,又能够避免因为k点过密而导致计算资源被浪费。在自洽场(SCF)迭代过程中,把电荷密度收敛标准规定为1×10⁻⁶ eV/atom,规定这个标准是为了确保电子结构能够达到稳定状态。

表2 第一性原理计算参数设置
计算参数类别参数名称具体设置
电子结构计算交换关联泛函Perdew-Burke-Ernzerhof (PBE)
电子结构计算平面波截断能500 eV
电子结构计算k点网格6×6×1 (Gamma-centered)
电子结构计算自洽场收敛精度1.0×10⁻⁶ eV/atom
几何优化离子弛豫收敛标准力收敛<0.01 eV/Å,位移收敛<1.0×10⁻³ Å
几何优化真空层厚度20 Å (避免层间相互作用)
缺陷模型构建缺陷类型单空位(Ti空位、C空位)、双空位(Ti-C相邻空位)
缺陷模型构建超胞尺寸3×3×1 Ti₃C₂ MXene超胞
储氢性能模拟氢吸附位点表面顶位(Ti顶、C顶)、桥位、空心位
储氢性能模拟吸附能计算方法E_ads = E(Ti₃C₂-H) - E(Ti₃C₂) - n×E(H₂)/2(n为吸附H原子数)
色散校正范德华力修正Grimme DFT-D3(BJ)方法
软件与平台计算软件VASP 5.4.4
软件与平台赝势类型投影缀加波(PAW)赝势

在结构优化阶段采用共轭梯度算法,共轭梯度算法是通过迭代的方式来调整原子坐标和晶格参数,从而实现最小化体系能量的目的。具体的收敛标准有这些内容:原子间作用力要小于0.05 eV/Å,应力张量各分量要小于0.1 GPa,相邻两次迭代的总能量变化要低于1×10⁻⁵ eV/atom。对于氢分子吸附体系,首先要固定Ti₃C₂表面结构,在固定好表面结构之后对H₂分子进行预优化,预优化完成之后再对整个吸附体系进行完全弛豫,并且要保证15 Å以上的真空层厚度始终保持不变,保持这样的真空层厚度是为了消除相邻镜像之间的相互作用。这些参数设置涉及到电子结构求解、几何优化和吸附能计算等关键步骤,通过合理设置这些参数能够为后续分析缺陷对储氢性能的影响提供可靠的数据支撑。

2.3储氢性能评估指标

评估材料的储氢性能是判断其储氢能力的重要步骤,这个步骤的核心是用一系列量化指标对材料的吸附特性和储氢潜力进行系统分析。

吸附能是首要的评估指标,它能直接体现H₂分子和材料表面结合的强弱程度。其具体的计算公式为:

其中\(E_{\text{total}}\)代表的是吸附体系的总能量,\(E_{\text{substrate}}\)是基底材料的能量,\(E_{\text{H}_2}\)是单个H₂分子的能量,\(n\)指的是吸附的分子数量。吸附能的负值大小能够体现吸附的稳定程度,一般来说,其绝对值需要在0.1到0.6 eV之间,这样才可以保证吸附过程是可逆的。这个参数可以通过第一性原理计算当中的能量优化步骤直接得到,它为挑选高效储氢材料提供了热力学方面的依据,因为通过这个参数可以知道材料在热力学上对于储氢是否有利。
分析吸附构型能够从几何层面揭示储氢的作用机制,主要是关注H₂分子在材料缺陷位置的吸附点、空间方向以及原子间的距离情况。在对优化后的结构进行可视化处理之后,就可以判断H₂分子更可能出现在缺陷顶部、桥位还是空位,同时还能记录它和最近邻原子的距离数据。把这个指标和吸附能结合起来,就能够判断分子属于物理吸附还是化学吸附,比如说H - H键长度的变化就直接体现了分子的解离情况。构型相关参数可以从结构优化的输出文件中提取出来,它是了解缺陷周围化学环境的重要方法,通过这些参数能知道缺陷周围的化学环境是怎样影响储氢过程的。
分析电荷转移需要借助Bader电荷的量化计算,以此来揭示H₂和缺陷之间电子相互作用的本质。在计算的时候,要对吸附体系做差分电荷密度分析,用公式\(\Delta \rho = \rho_{\text{ads}} - \rho_{\text{substrate}} - \rho_{\text{H}_2}\)得到电荷重新分布的图像,然后再结合Bader分区方法算出准确的电荷转移量。如果计算结果为正值,那就说明电子从基底转移到了H₂,这样会增强吸附作用;如果结果为负,情况则相反,会削弱吸附。这个指标能从微观层面解释吸附能存在差异的原因,它特别适合用来评估缺陷活性位点对H₂分子的活化能力,通过它可以知道缺陷活性位点是如何对H₂分子进行活化的。

储氢容量这个指标需要把单缺陷位的吸附数量和材料密度参数结合起来,以此来估算材料在实际应用中的潜力。理论容量可以用公式C=n×Ndefect×MH2Mmaterial×100%C = \frac{n \times N{\text{defect}} \times M{\text{H}2}}{M{\text{material}}} \times 100\%来计算,这里的NdefectN_{\text{defect}}是单位质量材料里的缺陷数量,MM代表分子量。在计算这个参数的时候,要综合前面提到的吸附能和构型分析结果,要保证每个H₂分子都符合可逆吸附的条件。容量的计算结果直接关系到材料的工程应用价值,它也为后续实验验证提供了明确的目标,后续实验就可以朝着这个目标去验证材料的实际储氢情况。

第三章结论

这项研究用第一性原理计算方法来做。它对Ti3C2 MXene表面缺陷怎样影响储氢性能进行了系统分析,明确了缺陷工程在改善储氢材料性能方面有着关键价值。在研究的时候构建了带有空位缺陷和官能团修饰的Ti3C2模型,计算得出的结果表明,表面缺陷会明显改变材料的电子结构以及吸附特性,然后提高氢分子的吸附能和储氢容量。当Ti空位缺陷出现时会引入局域电子态,这会加强氢分子与材料表面的相互作用,使得吸附能从 -0.15 eV增加到 -0.35 eV,而这个数值是接近理想储氢材料所需要的能量范围的。并且含氧官能团的修饰还能够进一步对吸附位点的分布情况进行优化,提升储氢的稳定性和可逆性。

从原理方面来讲,Ti3C2 MXene的储氢性能主要是依靠其特有的层状结构和表面化学活性。引入缺陷会打破原本的对称性,从而产生更多活性位点,还可以通过电荷转移机制增强氢分子的极化效应。这一过程不仅让储氢密度得到了提高,而且降低了吸放氢的能垒,为实际应用提供了理论上的支撑。在研究过程中采用了把密度泛函理论(DFT)与分子动力学模拟结合起来的方法,通过对缺陷浓度、类型以及官能团种类进行调控,系统评估了不同条件下的储氢性能,给实验设计明确了方向。

在实际应用这一块,Ti3C2 MXene因为具有高比表面积和可调控性,所以在氢能源存储领域存在很大的应用潜力。研究结果显示,通过精准地调控表面缺陷,有可能性开发出高效、安全的储氢材料,这对于推动氢能源汽车、便携式储能设备等技术的发展是有着重要意义的。并且缺陷工程策略还能够应用到其他二维材料体系当中,为新型储能材料设计提供具有普遍适用性的参考。这项研究不但加深了对Ti3C2 MXene储氢机制的认识,而且为未来实验研究和工业应用打下了非常扎实的理论基础。

参考文献