Ti0.33WO3电子结构和太阳辐射屏蔽性能第一性原理研究*

秦京运 舒群威 袁艺 仇伟 肖立华† 彭平 卢国松

1) (贵州理工学院, 贵阳 550003)

2) (湖南大学材料科学与工程学院, 长沙 410082)

节能减排已成为当今社会发展的主题, 对节约能源、提高太阳能的高效综合利用的新型窗用透明隔热材料的理论设计和研究尤其重要.本文采用基于密度泛函理论的计算方法, 研究了六方相三氧化钨Ti掺杂前、后的晶格参数、电子能带结构、形成能和光学性质.研究结果表明, Ti掺杂后晶格体积增大, 系统能量降为负值, 体系具有更好的稳定性; 掺杂后电子能带结构发生很大的变化, 但材料仍保持n型电导率; 随之, 其光学性质也发生改变, 掺杂前h−WO3无近红外吸收性能, 掺杂后的Ti0.33WO3具有很强的近红外吸收性能.在此基础上研究了Ti掺杂h−WO3前、后的太阳辐射屏蔽性能, 掺杂前无太阳辐射屏蔽性能; 掺杂后的Ti0.33WO3薄膜具有可见光高透明、近红外屏蔽的性能.计算结果为Ti掺杂h−WO3在窗用透明隔热材料方面的研究提供了理论依据.

1 引 言

在全球气候变暖的背景下, 节能减排迫在眉睫, 节能材料备受人们青睐.为了实现可见光高透明、近红外辐射强屏蔽, 从而减少室内和车内空调的致冷负荷和汽车尾气的排放, 达到节能减排的目的, 新型环保高透明太阳辐射屏蔽材料的设计和研究迫在眉睫.

迄今为止, 研究者相继开发了掺锑氧化锡(ATO)[1]、掺铟氧化锡 (ITO)[2]、稀土六硼化物(ReB6)[3,4]、氮化钛 (TiN)[5]、二氧化钒 (VO2)[6]、钨青铜类化合物 MxWO3[7−15](M=Ti, Na, K, Rb,Cs; 0 ＜ x ≤ 0.33)等透明隔热窗用近红外吸收材料, 但上述材料各有其优缺点.纳米ATO[1]只能屏蔽光波长1500—2500 nm的近红外辐射, 纳米ITO[2]只能屏蔽光波长1200—2500 nm的近红外辐射, 但对光波长780—1200 nm的高能近红外辐射无能为力.氮化钛[5]是非计量化合物, 它的组成范围为 TiN0.6—TiN1.16, 但其光学性能受N/Ti比率的影响很大.纳米VO2[6]虽然存在从高温金属相到低温绝缘体相的转换特性可用于智能温控窗, 但其透明隔热性能较差, 无法达到高性能窗用透明隔热材料的要求.纳米六硼化物ReB6, 特别是LaB6,虽然可屏蔽掉光波长780—1300 nm的高能近红外辐射[3,4], 但对光波长1300 nm以上的近红外辐射屏蔽性能较差.在诸如此类透明隔热窗用近红外吸收材料中, 由于MxWO3具有其它材料无法比拟的可见光高透明、近红外辐射强吸收屏蔽性能, 研究者对金属离子掺杂六方相三氧化钨(MxWO3)的透明隔热性能进行了大量的实验研究[7−15], 而理论研究较少.Yang等[16,17]采用HSE杂化密度泛函计算了碱金属单掺杂和共掺杂六方相三氧化钨的电子结构和光学性质, 但掺杂含量太低(0.083 %),远少于实验掺杂含量(0.15%— 0.33%); Lee等[18]也采用HSE杂化密度泛函计算了Na0.33WO3和K0.33WO3的光学性质, 但与实验值相差甚远;Yoshio和Adachi[19]采用密度泛函理论(DFT)+U方法研究了铯掺杂六方相三氧化钨的光学性质.Xu等[20]采用广义密度梯度近似(GGA)研究了Cs0.33WO3的光学性质, 理论研究结果与实验结果基本一致.

鉴于实验报道Ti掺杂h−WO3形成的Ti0.33WO3也具有太阳辐射屏蔽性能[7], 但其理论研究鲜有报道.因此, 本文采用DFT−GGA方法对Ti0.33WO3的能带结构、形成能、态密度和光学性质进行理论研究, 在此基础上预测其太阳辐射屏蔽性能.

2 晶体结构和计算方法

2.1 晶体结构

六方相三氧化钨(h−WO3)属于六方晶系, 空间群为P6/mmm (No.191), 具有WO6氧八面体六元环和三元环的ab平面内的沿c轴堆积, 形成平行c轴的六边形和三边形的准一维孔道[21].六方相三氧化钨中W原子位于(0.5, 0, 0), O原子位于(0.5, 0, 0.5)和(0.212, 0.424, 0).Ti离子嵌入在WO6八面体六元环内, 位置为(0, 0, 0)[22](如图1所示).

2.2 计算方法

图1 Ti掺杂 h−WO3的 2 × 2 × 1 超晶胞俯视图Fig.1.Top view of the 2 × 2 × 1 supercell of Ti−doped h−WO3.

本文所进行的第一性原理计算是采用基于密度泛函理论的CASTEP软件[23].计算首先采用BFGS算法对六方相三氧化物(h−WO3)的晶体结构进行几何结构优化, 然后建立2 × 2 × 1的超晶胞, 几何结构优化后, Ti原子掺杂在WO6八面体六元环内, 构成Ti0.33WO3晶体结构.在对其进行几何结构优化得到稳定晶体体系的基础上, 进行能带结构、态密度和光学性质的计算.交换关联势采用广义梯度近似(GGA)[24]中的PBE[25]提出的形式, 描述离子实与价电子之间的相互作用势采用超软赝势进行.计算中Ti, W和O各原子的电子价态分别为6s26p1, 5p65d46s2, 2s2p4, 其它轨道电子作为芯电子进行计算.平面波截止能和布里渊区k点的设置分别为500 eV和9 × 9 × 15.迭代过程中的自洽收敛精度为5.0 × 10—6eV/atom, 晶体内应力收敛标准为0.02 GPa, 最大位移为5.0 ×10—4Å以及原子间相互作用力收敛标准为小于0.01 eV/Å.h−WO3的电子能带结构和光学性质计算均采用了剪刀算符修正带隙值.h−WO3的带隙计算值为0.62 eV, 而实验值为2.6—3.25 eV[26−29],因此, 剪刀算符设为2.0 eV.

3 计算结果与讨论

3.1 几何结构优化结果

表1为采用基于第一性原理的广义密度梯度近似方法, 对h−WO3和Ti0.33WO3的进行结构优化及能带结构计算得到的晶格参数和能带间隙.其中, Ti0.33WO3的晶格参数为Ti空位掺杂六方相三氧化钨(h−WO3) 的超晶胞结构几何优化后, 折合成单晶胞的结构参数.计算得到的纯h−WO3的晶格参数和带隙值的计算误差均小于2.0 %, 与实验和理论值基本一致.Ti掺杂后晶格体积比纯h−WO3晶体稍有增大, 这主要是六元环的半径为0.163 nm[30], 而 Ti的离子半径为 0.184 nm[31], 当嵌入的离子半径大于六元环的半径时, 造成晶格畸变, 从而使晶格体积增大.

表1 Ti掺杂六方相WO3前、后的晶格参数、带隙与形成能Table 1.Lattice parameters, band gap and formation energy of pure h−WO3 before and after Ti−doped.

3.2 形成能分析

通过计算Ti掺杂纯h−WO3的形成能可分析其相对稳定性[34,35], 计算公式为

其中Edoped和Epure分别表示空位掺杂模型与纯h−WO3超胞模型的总能量; µTi为Ti原子的化学势;n为掺杂进入纯h−WO3的Ti原子的个数.形成能的计算结果列于表1.结果表明, Ti空位掺杂纯h−WO3的形成能为—2.359 eV, 负的形成能表明Ti空位掺杂纯h−WO3体系具有良好的稳定性.

3.3 能带结构和态密度分析

3.3.1 纯h−WO3能带结构和态密度分析

为方便比较和分析空位掺杂对纯h−WO3能带结构和态密度的影响, 先计算了纯h−WO3的能带结构和态密度如图2所示.从图2(a)的能带结构计算结果可知, 纯h−WO3带隙较大, 导带底和价带顶分别位于Brillouin区高对称点Γ和点A处,为间接带隙半导体.其间接带隙计算值为0.62 eV,经剪刀算符2.0 eV修剪后间接带隙值为2.62 eV.由图2纯h−WO3的总态密度可知, 价带和价带顶态密度则主要由O−2p组成, 导带和导带底的态密度主要是W−5d态构成.从分态密度可知, 纯h−WO3的价带和导带都主要由局域化的O−2p和W−5d轨道电子构成并发生杂化, 表明O原子与W原子之间共价键比较明显, 这一特性与文献的计算结果一致.

3.3.2 Ti0.33WO3能带结构和态密度分析

图2 纯h−WO3的能带结构和态密度: (a)能带结构; (b)总态密度; (c) W的分态密度; (d) O1的态密度; (e) O2的分态密度Fig.2.Energy band structure and DOS of pure h−WO3: (a) Energy band structure; (b) TDOS of h−WO3; (c) PDOS of W;(d) PDOS of O1; (e) PDOS of O2.

图3 Ti掺杂h−WO3(Ti0.33WO3)能带结构和态密度: (a)能带结构; (b)总态密度; (c) Ti的分态密度; (d) W的分态密度;(e) O1的态密度; (f) O2的分态密度Fig.3.Energy band structure and DOS of Ti0.33WO3: (a) Energy band structure; (b) TDOS of Ti0.33WO3; (c) PDOS of Ti;(d) PDOS of W; (d) PDOS of O1; (e) PDOS of O2.

Ti空位掺杂h−WO3能带结构和态密度如图3所示.对比图3与图2的分态密度可知, Ti0.33WO3价带主要由 O−2p、W−5d 和 Ti−6s态组成, 导带主要W−5d和Ti−6p态组成, 尤为重要的是W−5d越过费米能级EF进入价带, 价带顶和导带底均由W−5d态构成.这表明Ti掺杂h−WO3主要扮演着捐献电子的角色.从图3(a)可知, 价带顶和导带底的位置没有改变; 但价带整体下移, 导带底越过费米能级EF形成简并态, 并且导带和价带之间仍存在一定的间隙, 这表明Ti空位掺杂构成的Ti0.33WO3晶体仍保持n型电导率.h−WO3掺杂前后费米能级的变化表明晶体从半导体向导体转变,这必将导致材料的光学性质的改变[36,37].后文中的光学性质的计算结果将进一步印证这一结论.

3.4 光学性质

为研究Ti掺杂前后h−WO3的光学性质, 对Ti空位掺杂前后六方相h−WO3的介电函数进行了计算研究.在计算h−WO3的光学性质时, 首先采用剪刀算符对其电子结构进行了修正, 计算光学性质时也采用相同剪刀算符进行修正, 本文h−WO3的电子结构和光学性质计算的剪刀算符取值为 2.0 eV.介电函数表达式为 ε(ω) = ε1(ω) +iε2(ω), 相关计算公式如下[38]:

式中, C, V代表导带、价带; BZ表示第一布里渊区; ℏ 是狄拉克常量; K是倒格矢; |a·MC,V(K)|2是动量矩阵元; EC(K), EV(K)代表导带、价带本征能级.根据克拉默斯−克勒尼希(Kramers−Kronig)色散关系及电子跃迁概率的定义利用介电函数的实部 ε1(ω)与虚部 ε2(ω)可推导出固体的反射率R(ω)和吸收系数 α(ω).

图4 (a)为Ti掺杂h−WO3前后的介电函数虚部ε2(ω), 虚部与带间跃迁有关, 取决于导带与价带的跃迁, 可通过能带结构与态密度解释.从图4 (a)可以看出, h−WO3的介电函数阈能约为2.6 eV, 这是h−WO3的基本吸收边, 当光子能量高于此值时,虚部曲线迅速上升.纯h−WO3的介电函数有两个特征峰, 峰1和峰2分别对应光子能量为5.8 eV和 11.4 eV.峰 1对应 O−2p态电子向 W−5d态直接跃迁所致; 峰2是下价带能级电子向上导带能级跃迁所致.Ti原子空位掺杂后, Ti0.33WO3在低能0.1 eV处增加新的最大特征峰, 这是价带顶W−5d态向导带底W−5d态跃迁所致, 其它两个特征峰向低能方向移动.

采用第一性原理计算出来的材料的反射谱R(ω)和吸收谱α(ω)代入公式[39]:

(其中d为薄膜厚度)可计算材料致密薄膜和涂层薄膜的太阳辐射屏蔽性能.需要指出的是材料致密薄膜的理论透光率是根据反射谱和吸收谱共同作用计算出来的结果; 材料涂层薄膜的理论透光率是只根据吸收谱计算出来的结果, 因为涂层薄膜的反射率极低[40].图5为200 nm厚h−WO3和Ti0.33WO3薄膜的理论透过率.从图5(a)中可以看出h−WO3致密薄膜对太阳辐射从紫外−可见光到近红外的辐射屏蔽作用很弱; 但其掺杂后(Ti0.33WO3)对可见光有很高的透光率, 在光波长711 nm处最大透光率为61.9%; 而对近红外辐射有很强的屏蔽作用,在光波长1607 nm处最小透光率为2.0%, 两者差值约为60%.从图5(b)中掺杂h−WO3前后涂层薄膜的理论透过率可以看出, h−WO3涂层薄膜对太阳辐射无屏蔽作用; 而掺杂后Ti0.33WO3涂层的理论透过率曲线呈“吊铃”型分布, 在可见光678 nm处透光率最大值可达82.8%, 在近红外1524 nm处透光率最小值为5.0%.两者相差77.8%.从以上分析可以看出, h−WO3掺杂前对太阳辐射几乎无屏蔽作用, 而Ti掺杂后的Ti0.33WO3无论是致密薄膜还是涂层薄膜对太阳辐射均有很好的可见光高透明、近红外强屏蔽作用.涂层薄膜的理论透光率与文献[7]实验结果基本一致.

图4 h−WO3和Ti0.33WO3的光学性质: (a)介电函数虚部; (b)反射谱; (c)吸收谱Fig.4.Optical performance of h−WO3 andTi0.33WO3: (a) Imaginary part of the dielectric function; (b) reflectivity; (c) absorption spectrum.

图5 h−WO3和Ti0.33WO3薄膜的理论透过率: (a)致密薄膜; (b)涂层薄膜Fig.5.Theoretical transmittance of h−WO3 andTi0.33WO3 films: (a) The compacted film; (b) the coated film.

4 结 论

本文采用第一性原理对h−WO3及掺杂后形成的Ti0.33WO3材料的晶格参数、形成能、电子结构和光学性质进行了计算和分析, 并预测了其薄膜的理论透光率.计算结果表明: 1) Ti掺杂h−WO3后晶格体积增大, 系统能量为负值, 晶体具有更好的结构稳定性; 2) Ti掺杂半导体h−WO3后晶体能带结构发生变化, 构成的Ti0.33WO3晶体保持n型电导率; 3) h−WO3掺杂前无近红外吸收/反射性能,而掺杂后形成的Ti0.33WO3具有很强的近红外吸收/反射性能, 可用作可见光高透明、近红外强吸收屏蔽的窗用透明隔热材料.