搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

一种新型二维TiO2的电子结构与光催化性质

熊子谦 张鹏程 康文斌 方文玉

一种新型二维TiO2的电子结构与光催化性质

熊子谦, 张鹏程, 康文斌, 方文玉
PDF
HTML
导出引用
  • 基于第一性原理计算方法, 设计出了一种新型二维半导体材料TiO2, 并进一步研究了其结构稳定性, 电子结构, 载流子迁移率和光学性质等. 二维TiO2的形成能、声子谱、分子动力学、弹性常数表明, 二维TiO2具有较好的动力学, 热力学和机械稳定性, 具备实验制备的条件, 且能够稳定存在于常温条件下. 电子结构分析表明, 二维TiO2是一种间接带隙半导体, 在GGA+PBE和HSE06算法下的能隙分别为1.19 eV和2.76 eV, 其价带顶和导带底能级分别由Ti-3d和Ti-4s态电子构成, O原子的电子态在费米能级附近贡献很小, 主要分布在深处能级. 载流子迁移率显示, 二维TiO2的迁移率比单层MoS2要小, 其电子和空穴迁移率分别为31.09和36.29 cm2·V–1·s–1. 由于空穴迁移率和电子迁移率的各向异性, 电子-空穴复合率较低, 使得单层TiO2的使用寿命更长, 光催化活性更好. 在应变调控下, 二维TiO2的能隙发生明显响应, 以适用于各种半导体器件的需要. 半导体的带边势和光学性质显示, 在–5%—2%单/双轴应变下, 二维TiO2能够光裂水制H2, 在–5%—5%单/双轴应变下, 能够光裂水制O2, H2O2和O3等. 此外, 二维TiO2对可见光和紫外光具有较高的吸收系数, 说明其在未来光电子器件和光催化材料领域有着潜在的应用前景.
      通信作者: 康文斌, wbkang@hbmu.edu.cn ; 方文玉, jzfangwenyu@163.com
    [1]

    Hashimoto A, Suenaga K, Gloter A, Urita K, Iijima S 2004 Nature 430 870

    [2]

    Liu L, Zhang J, Zhao J, Liu F 2012 Nanoscale 4 5910

    [3]

    Kvashnin D G, Bellucci S, Chernozatonskii L A 2015 Phys. Chem. Chem. Phys. 17 4354

    [4]

    Liu Y, Duan X, Huang Y, Duan X 2018 Chem. Soc. Rev. 47 6388

    [5]

    Gupta S, Kutana A, Yakobson B I 2018 J. Phys. Chem. Lett. 9 2757

    [6]

    Yuan J, Yu N, Xue K, Miao X 2017 Appl. Surf. Sci. 409 85

    [7]

    Lee C, Wei X, Kysar J W, Hone J 2008 Science 321 385

    [8]

    Wang H, Chan T L, Xie Z 2018 Chem. Commun. (Camb) 54 385

    [9]

    Hartman T, Sofer Z 2019 ACS Nano 13 8566

    [10]

    Li P 2019 Phys. Chem. Chem. Phys. 21 11150

    [11]

    Li L, Yang J 2017 Nanotechnology 28 475701

    [12]

    Yuan J, Xie Q, Yu N, Wang J 2017 Appl. Surf. Sci. 394 625

    [13]

    Zhang S, Yan Z, Li Y, Chen Z, Zeng H 2015 Angew. Chem. Int. Ed. Engl. 54 3112

    [14]

    Zhang D, Xiong Y, Cheng J, Chai J, Liu T, Ba X, Ullah S, Zheng G, Yan M, Cao M S 2020 Sci. Bull. 65 138

    [15]

    Cao M S, Shu J C, Wang X X, Wang X, Zhang M, Yang H J, Fang X Y, Yuan J 2019 Ann. Phys.-Berlin 531 1800390

    [16]

    Cao M S, Wang X X, Zhang M, Shu J C, Cao W Q, Yang H J, Fang X Y, Yuan J 2019 Adv. Funct. Mater. 29 1807398

    [17]

    Matta S K, Zhang C, Jiao Y, O'mullane A, Du A 2018 Nanoscale 10 6369

    [18]

    Zhu Y L, Yuan J H, Song Y Q, Wang S, Xue K H, Xu M, Cheng X M, Miao X S 2019 J. Mater. Sci. 54 11485

    [19]

    Sun Y, Cheng H, Gao S, Sun Z, Liu Q, Liu Q, Lei F, Yao T, He J, Wei S, Xie Y 2012 Angew. Chem. Int. Ed. Engl. 51 8727

    [20]

    方文玉, 张鹏程, 赵军, 康文斌 2020 物理学报 69 056301

    Fang W Y, Zhang P C, Zhao J, Kang W B 2020 Acta Phys. Sin. 69 056301

    [21]

    Chen Y, Wang L, Wang W, Cao M S 2017 Appl. Catal., B Environ. 209 110

    [22]

    Zhang Y, Kuwahara Y, Mori K, Yamashita H 2019 Chem. Asian J. 14 278

    [23]

    Liang Z, Sun B, Xu X, Cui H, Tian J 2019 Nanoscale 11 12266

    [24]

    Zeng H, Cui X 2015 Chem. Soc. Rev. 44 2629

    [25]

    Guo Z, Zhou J, Zhu L, Sun Z 2016 J. Mater. Chem. A 4 11446

    [26]

    Zhuang H L, Hennig R G 2013 Chem. Mater. 25 3232

    [27]

    Zhang X, Zhao X, Wu D, Jing Y, Zhou Z 2016 Adv. Sci. 3 1600062

    [28]

    Lacerda A M, Larrosa I, Dunn S 2015 Nanoscale 7 12331

    [29]

    Ishida Y, Motokane Y, Tokunaga T, Yonezawa T 2015 Phys. Chem. Chem. Phys. 17 24556

    [30]

    Yoshida T, Niimi S, Yamamoto M, Nomoto T, Yagi S 2015 J. Colloid Interface Sci. 447 278

    [31]

    Liu X, Chen Z, Cao M S 2019 ACS Appl. Energy Mater. 2 5960

    [32]

    Yang J, Jiang Y L, Li L J, Muhire E, Gao M Z 2016 Nanoscale 8 8170

    [33]

    Song S S, Xia B Y, Chen J, Yang J, Shen X, Fan S J, Guo M L, Sun Y M, Zhang X D 2014 RSC Adv. 4 42598

    [34]

    Liu X, Chen Z, Li W, Cao M S 2017 J. Phys. Chem. C 121 20605

    [35]

    Song X F, Hu L F, Li D H, Chen L, Sun Q Q, Zhou P, Zhang D W 2015 Sci Rep 5 15989

    [36]

    Tao J, Luttrell T, Batzill M 2011 Nat. Chem. 3 296

    [37]

    Xie Q, Yuan J, Yu N, Wang L, Wang J 2017 Comput. Mater. Sci. 135 160

    [38]

    Qiu G, Xiao Q, Hu Y, Qin W, Wang D 2004 J. Colloid Interface Sci. 270 127

    [39]

    Yuan J H, Song Y Q, Chen Q, Xue K H, Miao X S 2019 Appl. Surf. Sci. 469 456

    [40]

    Castellanos-Gomez A, Poot M, Steele G A, Van Der Zant H S, Agrait N, Rubio-Bollinger G 2012 Adv. Mater. 24 772

    [41]

    Yuan J, Yu N, Xue K, Miao X 2017 RSC Adv. 7 8654

    [42]

    Song Y Q, Yuan J H, Li L H, Xu M, Wang J F, Xue K H, Miao X S 2019 Nanoscale 11 1131

    [43]

    Xu L C, Du A, Kou L 2016 Phys. Chem. Chem. Phys. 18 27284

    [44]

    Peng R, Ma Y, He Z, Huang B, Kou L, Dai Y 2019 Nano Lett. 19 1227

    [45]

    Li P, You Z, Haugstad G, Cui T 2011 Appl. Phys. Lett. 98 253105

    [46]

    Wang Y, Ma R, Hu K, Kim S, Fang G, Shao Z, Tsukruk V V 2016 ACS Appl. Mater. Interfaces 8 24962

    [47]

    Li Y, Yu C, Gan Y, Kong Y, Jiang P, Zou D F, Li P, Yu X F, Wu R, Zhao H, Gao C F, Li J 2019 Nanotechnology 30 335703

    [48]

    Liu F, Ming P, Li J 2007 Phys. Rev. B 76 064120

    [49]

    Kudin K N, Scuseria G E, Yakobson B I 2001 Phys. Rev. B 64 235406

    [50]

    Kang J, Sahin H, Peeters F M 2015 Phys. Chem. Chem. Phys. 17 27742

    [51]

    Zhao J, Li Y, Ma J 2016 Nanoscale 8 9657

    [52]

    Yuan J, Yu N, Wang J, Xue K-H, Miao X 2018 Appl. Surf. Sci. 436 919

    [53]

    El Mragui A, Logvina Y, Pinto Da Silva L, Zegaoui O, Esteves Da Silva J C G 2019 Materials 1 2

    [54]

    An X, Hu C, Liu H, Qu J 2018 Langmuir 34 1883

    [55]

    Luican-Mayer A, Zhang Y, Dilullo A, Li Y, Fisher B, Ulloa S E, Hla S W 2019 Nanoscale 11 22351

    [56]

    Yu W, Zhang J, Peng T 2016 Appl. Catal., B Environ. 181 220

    [57]

    Zhong L, Chen X, Qi J 2017 Phys. Chem. Chem. Phys. 19 15388

    [58]

    Hua C, Sheng F, Hu Q, Xu Z A, Lu Y, Zheng Y 2018 J. Phys. Chem. Lett. 9 6695

    [59]

    Sarker H P, Rao P M, Huda M N 2019 ChemPhysChem 20 773

    [60]

    Fang W Y, Li P A, Yuan J H, Xue K H, Wang J F 2019 J. Electron. Mater. 49 959

    [61]

    Liu X, Wang Y, Li F, Li Y 2016 Phys. Chem. Chem. Phys. 18 14638

    [62]

    Xiao J, Long M, Li M, Li X, Xu H, Chan K 2015 Phys. Chem. Chem. Phys. 17 6865

    [63]

    Zhang J, Wageh S, Al-Ghamdi A, Yu J 2016 Appl. Catal., B Environ. 192 101

    [64]

    Mogulkoc A, Mogulkoc Y, Kecik D, Durgun E 2018 Phys. Chem. Chem. Phys. 20 21043

  • 图 1  二维TiO2的晶体结构 (a) 俯视图; (b) 侧视图; (c) K点路径

    Fig. 1.  Crystal structure of monolayer TiO2 (a) top view; (b) side view; (c) K point path.

    图 2  二维TiO2的声子谱

    Fig. 2.  Phonon spectrum of monolayer TiO2.

    图 3  二维TiO2的分子动力学模拟

    Fig. 3.  Molecular dynamics simulation of monolayer TiO2.

    图 4  杨氏模量(左图/蓝色)和泊松比(右图/红色)极坐标图, $\theta $为相对于a轴方向的夹角

    Fig. 4.  Polar coordinates for Young's modulus (left/blue line) and Poisson's ratio (right/red line) $\theta $ is the angle with respect to the a-direction.

    图 5  能带图(左)和分波态密度图(右) (a) 单层P-6M2 TiO2; (b) 单层金红石TiO2; (c) 单层锐钛矿TiO2

    Fig. 5.  Energy band structure (left) and density of states (right) of (a) monolayer P-6M2 TiO2, (b) monolayer rutile TiO2, (c) monolayer anatase TiO2.

    图 6  二维TiO2的差分电荷密度图

    Fig. 6.  The differential charge density of monolayer TiO2.

    图 7  (a) 二维TiO2沿a/b方向的总能量与应变量$\Delta l/l$的关系, 采用二次数据拟合二维结构的平面刚度, 黑色和红色曲线表示沿ab方向的面内刚度; (b), (c)单层TiO2的VBM和CBM随应变量相对真空能级的变化, 采取线性拟合计算形变势

    Fig. 7.  (a) The relation between total energy and the applied strain $\Delta l/l$ along the a/b directions of monolayer TiO2. The quadratic data fitting gives the in-plane stiffness of 2D structures. Black and red curves show the in-plane stiffness along the a and b directions of monolayer TiO2. The shift of VBMs and CBMs for (b-c) monolayer TiO2 with respect to the vacuum energy, as a function of the applied strain along either the a and b direction. The linear fit of the data yields the deformation potential constant.

    图 8  单/双轴应变下能隙变化

    Fig. 8.  Band gap of monolayer TiO2 under uniaxial/biaxial strain, calculated using the HSE06 functional.

    图 9  单双/轴应变下单层TiO2光催化示意图

    Fig. 9.  Schematic diagram of monolayer TiO2 photocatalysis under uniaxial/biaxial strain.

    图 10  光学性质 (a) 介电函数虚部; (b) 光吸收系数

    Fig. 10.  Optical properties: (a) virtual part of dielectric function; (b) absorption coefficient.

    表 1  二维TiO2的结构常数和结合能

    Table 1.  Structure constants and binding energy of monolayer TiO2.

    Materiala/b${\theta _1}$/($^ \circ $)${\theta _2}$/($^ \circ $)l$\sigma $/ÅEf/eV
    TiO22.8991.6668.162.012.26–8.11
    MoS23.1882.5880.742.413.14–7.35
    下载: 导出CSV

    表 2  二维TiO2有效质量${m^ * }$, 形变势常数$E_{_{\rm{d}}}^i$, 弹性常数${C^{{\rm{2 D}}}}$和载流子迁移率${\mu _{{\rm{2 D}}}}$

    Table 2.  Calculated effective mass ${m^ * }$, deformation potential constant $E_{_{\rm{d}}}^i$, elastic modulus ${C^{{\rm{2 D}}}}$, and carrier mobility ${\mu _{{\rm{2 D}}}}$ for monolayer TiO2 along the a ($\zeta \to K$) and b ($\zeta \to G$) directions, where $\zeta $ represents the position of the valence band top and the conduction band bottom.

    Carrier type$m_a^ * $/${m_{\rm{e}}}$$m_b^ * $/${m_{\rm{e}}}$$m_l^ * $/${m_{\rm{e}}}$$\left| {{E_{la}}} \right|$/ eV$\left| {{E_{lb}}} \right|$/ eV$C_a^{2{\rm{D}}}$/ N·m–1$C_b^{{\rm{2 D}}}$/ N·m–1$\mu _a^{2{\rm{D}}}$/ cm2·V–1·s–1$\mu _b^{2{\rm{D}}}$/ cm2·V–1·s–1
    Electrons3.211.392.113.433.3821.2721.2812.9230.75
    Holes4.734.124.411.261.2521.2721.2831.0936.29
    下载: 导出CSV
  • [1]

    Hashimoto A, Suenaga K, Gloter A, Urita K, Iijima S 2004 Nature 430 870

    [2]

    Liu L, Zhang J, Zhao J, Liu F 2012 Nanoscale 4 5910

    [3]

    Kvashnin D G, Bellucci S, Chernozatonskii L A 2015 Phys. Chem. Chem. Phys. 17 4354

    [4]

    Liu Y, Duan X, Huang Y, Duan X 2018 Chem. Soc. Rev. 47 6388

    [5]

    Gupta S, Kutana A, Yakobson B I 2018 J. Phys. Chem. Lett. 9 2757

    [6]

    Yuan J, Yu N, Xue K, Miao X 2017 Appl. Surf. Sci. 409 85

    [7]

    Lee C, Wei X, Kysar J W, Hone J 2008 Science 321 385

    [8]

    Wang H, Chan T L, Xie Z 2018 Chem. Commun. (Camb) 54 385

    [9]

    Hartman T, Sofer Z 2019 ACS Nano 13 8566

    [10]

    Li P 2019 Phys. Chem. Chem. Phys. 21 11150

    [11]

    Li L, Yang J 2017 Nanotechnology 28 475701

    [12]

    Yuan J, Xie Q, Yu N, Wang J 2017 Appl. Surf. Sci. 394 625

    [13]

    Zhang S, Yan Z, Li Y, Chen Z, Zeng H 2015 Angew. Chem. Int. Ed. Engl. 54 3112

    [14]

    Zhang D, Xiong Y, Cheng J, Chai J, Liu T, Ba X, Ullah S, Zheng G, Yan M, Cao M S 2020 Sci. Bull. 65 138

    [15]

    Cao M S, Shu J C, Wang X X, Wang X, Zhang M, Yang H J, Fang X Y, Yuan J 2019 Ann. Phys.-Berlin 531 1800390

    [16]

    Cao M S, Wang X X, Zhang M, Shu J C, Cao W Q, Yang H J, Fang X Y, Yuan J 2019 Adv. Funct. Mater. 29 1807398

    [17]

    Matta S K, Zhang C, Jiao Y, O'mullane A, Du A 2018 Nanoscale 10 6369

    [18]

    Zhu Y L, Yuan J H, Song Y Q, Wang S, Xue K H, Xu M, Cheng X M, Miao X S 2019 J. Mater. Sci. 54 11485

    [19]

    Sun Y, Cheng H, Gao S, Sun Z, Liu Q, Liu Q, Lei F, Yao T, He J, Wei S, Xie Y 2012 Angew. Chem. Int. Ed. Engl. 51 8727

    [20]

    方文玉, 张鹏程, 赵军, 康文斌 2020 物理学报 69 056301

    Fang W Y, Zhang P C, Zhao J, Kang W B 2020 Acta Phys. Sin. 69 056301

    [21]

    Chen Y, Wang L, Wang W, Cao M S 2017 Appl. Catal., B Environ. 209 110

    [22]

    Zhang Y, Kuwahara Y, Mori K, Yamashita H 2019 Chem. Asian J. 14 278

    [23]

    Liang Z, Sun B, Xu X, Cui H, Tian J 2019 Nanoscale 11 12266

    [24]

    Zeng H, Cui X 2015 Chem. Soc. Rev. 44 2629

    [25]

    Guo Z, Zhou J, Zhu L, Sun Z 2016 J. Mater. Chem. A 4 11446

    [26]

    Zhuang H L, Hennig R G 2013 Chem. Mater. 25 3232

    [27]

    Zhang X, Zhao X, Wu D, Jing Y, Zhou Z 2016 Adv. Sci. 3 1600062

    [28]

    Lacerda A M, Larrosa I, Dunn S 2015 Nanoscale 7 12331

    [29]

    Ishida Y, Motokane Y, Tokunaga T, Yonezawa T 2015 Phys. Chem. Chem. Phys. 17 24556

    [30]

    Yoshida T, Niimi S, Yamamoto M, Nomoto T, Yagi S 2015 J. Colloid Interface Sci. 447 278

    [31]

    Liu X, Chen Z, Cao M S 2019 ACS Appl. Energy Mater. 2 5960

    [32]

    Yang J, Jiang Y L, Li L J, Muhire E, Gao M Z 2016 Nanoscale 8 8170

    [33]

    Song S S, Xia B Y, Chen J, Yang J, Shen X, Fan S J, Guo M L, Sun Y M, Zhang X D 2014 RSC Adv. 4 42598

    [34]

    Liu X, Chen Z, Li W, Cao M S 2017 J. Phys. Chem. C 121 20605

    [35]

    Song X F, Hu L F, Li D H, Chen L, Sun Q Q, Zhou P, Zhang D W 2015 Sci Rep 5 15989

    [36]

    Tao J, Luttrell T, Batzill M 2011 Nat. Chem. 3 296

    [37]

    Xie Q, Yuan J, Yu N, Wang L, Wang J 2017 Comput. Mater. Sci. 135 160

    [38]

    Qiu G, Xiao Q, Hu Y, Qin W, Wang D 2004 J. Colloid Interface Sci. 270 127

    [39]

    Yuan J H, Song Y Q, Chen Q, Xue K H, Miao X S 2019 Appl. Surf. Sci. 469 456

    [40]

    Castellanos-Gomez A, Poot M, Steele G A, Van Der Zant H S, Agrait N, Rubio-Bollinger G 2012 Adv. Mater. 24 772

    [41]

    Yuan J, Yu N, Xue K, Miao X 2017 RSC Adv. 7 8654

    [42]

    Song Y Q, Yuan J H, Li L H, Xu M, Wang J F, Xue K H, Miao X S 2019 Nanoscale 11 1131

    [43]

    Xu L C, Du A, Kou L 2016 Phys. Chem. Chem. Phys. 18 27284

    [44]

    Peng R, Ma Y, He Z, Huang B, Kou L, Dai Y 2019 Nano Lett. 19 1227

    [45]

    Li P, You Z, Haugstad G, Cui T 2011 Appl. Phys. Lett. 98 253105

    [46]

    Wang Y, Ma R, Hu K, Kim S, Fang G, Shao Z, Tsukruk V V 2016 ACS Appl. Mater. Interfaces 8 24962

    [47]

    Li Y, Yu C, Gan Y, Kong Y, Jiang P, Zou D F, Li P, Yu X F, Wu R, Zhao H, Gao C F, Li J 2019 Nanotechnology 30 335703

    [48]

    Liu F, Ming P, Li J 2007 Phys. Rev. B 76 064120

    [49]

    Kudin K N, Scuseria G E, Yakobson B I 2001 Phys. Rev. B 64 235406

    [50]

    Kang J, Sahin H, Peeters F M 2015 Phys. Chem. Chem. Phys. 17 27742

    [51]

    Zhao J, Li Y, Ma J 2016 Nanoscale 8 9657

    [52]

    Yuan J, Yu N, Wang J, Xue K-H, Miao X 2018 Appl. Surf. Sci. 436 919

    [53]

    El Mragui A, Logvina Y, Pinto Da Silva L, Zegaoui O, Esteves Da Silva J C G 2019 Materials 1 2

    [54]

    An X, Hu C, Liu H, Qu J 2018 Langmuir 34 1883

    [55]

    Luican-Mayer A, Zhang Y, Dilullo A, Li Y, Fisher B, Ulloa S E, Hla S W 2019 Nanoscale 11 22351

    [56]

    Yu W, Zhang J, Peng T 2016 Appl. Catal., B Environ. 181 220

    [57]

    Zhong L, Chen X, Qi J 2017 Phys. Chem. Chem. Phys. 19 15388

    [58]

    Hua C, Sheng F, Hu Q, Xu Z A, Lu Y, Zheng Y 2018 J. Phys. Chem. Lett. 9 6695

    [59]

    Sarker H P, Rao P M, Huda M N 2019 ChemPhysChem 20 773

    [60]

    Fang W Y, Li P A, Yuan J H, Xue K H, Wang J F 2019 J. Electron. Mater. 49 959

    [61]

    Liu X, Wang Y, Li F, Li Y 2016 Phys. Chem. Chem. Phys. 18 14638

    [62]

    Xiao J, Long M, Li M, Li X, Xu H, Chan K 2015 Phys. Chem. Chem. Phys. 17 6865

    [63]

    Zhang J, Wageh S, Al-Ghamdi A, Yu J 2016 Appl. Catal., B Environ. 192 101

    [64]

    Mogulkoc A, Mogulkoc Y, Kecik D, Durgun E 2018 Phys. Chem. Chem. Phys. 20 21043

  • [1] 谢知, 程文旦. TiO2纳米管电子结构和光学性质的第一性原理研究. 物理学报, 2014, 63(24): 243102. doi: 10.7498/aps.63.243102
    [2] 程旭东, 吴海信, 唐小路, 王振友, 肖瑞春, 黄昌保, 倪友保. Na2Ge2Se5电子结构和光学性质的第一性原理研究. 物理学报, 2014, 63(18): 184208. doi: 10.7498/aps.63.184208
    [3] 王寅, 冯庆, 王渭华, 岳远霞. 碳-锌共掺杂锐钛矿相TiO2 电子结构与光学性质的第一性原理研究 . 物理学报, 2012, 61(19): 193102. doi: 10.7498/aps.61.193102
    [4] 李旭珍, 谢泉, 陈茜, 赵凤娟, 崔冬萌. OsSi2电子结构和光学性质的研究. 物理学报, 2010, 59(3): 2016-2021. doi: 10.7498/aps.59.2016
    [5] 潘凤春, 林雪玲, 陈焕铭. C掺杂金红石相TiO2的电子结构和光学性质的第一性原理研究. 物理学报, 2015, 64(22): 224218. doi: 10.7498/aps.64.224218
    [6] 崔冬萌, 谢泉, 陈茜, 赵凤娟, 李旭珍. Si基外延Ru2Si3电子结构及光学性质研究. 物理学报, 2010, 59(3): 2027-2032. doi: 10.7498/aps.59.2027
    [7] 宋庆功, 刘立伟, 赵辉, 严慧羽, 杜全国. YFeO3的电子结构和光学性质的第一性原理研究. 物理学报, 2012, 61(10): 107102. doi: 10.7498/aps.61.107102
    [8] 程和平, 但加坤, 黄智蒙, 彭辉, 陈光华. 黑索金电子结构和光学性质的第一性原理研究. 物理学报, 2013, 62(16): 163102. doi: 10.7498/aps.62.163102
    [9] 骆最芬, 岑伟富, 范梦慧, 汤家俊, 赵宇军. BiTiO3电子结构及光学性质的第一性原理研究. 物理学报, 2015, 64(14): 147102. doi: 10.7498/aps.64.147102
    [10] 段满益, 周海平, 沈益斌, 陈青云, 丁迎春, 祝文军, 徐 明. 过渡金属与氮共掺杂ZnO电子结构和光学性质的第一性原理研究. 物理学报, 2007, 56(9): 5359-5365. doi: 10.7498/aps.56.5359
    [11] 毕艳军, 郭志友, 孙慧卿, 林 竹, 董玉成. Co和Mn共掺杂ZnO电子结构和光学性质的第一性原理研究. 物理学报, 2008, 57(12): 7800-7805. doi: 10.7498/aps.57.7800
    [12] 杨春燕, 张蓉, 张利民, 可祥伟. 0.5NdAlO3-0.5CaTiO3电子结构及光学性质的第一性原理计算. 物理学报, 2012, 61(7): 077702. doi: 10.7498/aps.61.077702
    [13] 赵佰强, 张耘, 邱晓燕, 王学维. Cu,Fe掺杂LiNbO3晶体电子结构和光学性质的第一性原理研究. 物理学报, 2016, 65(1): 014212. doi: 10.7498/aps.65.014212
    [14] 陈秋云, 赖新春, 王小英, 张永彬, 谭世勇. UO2的电子结构及光学性质的第一性原理研究. 物理学报, 2010, 59(7): 4945-4949. doi: 10.7498/aps.59.4945
    [15] 侯芹英, 孔祥兰, 苏希玉, 齐延华, 支晓芬. Ba0.5Sr0.5TiO3电子结构和光学性质的第一性原理研究. 物理学报, 2009, 58(6): 4128-4131. doi: 10.7498/aps.58.4128
    [16] 李沛娟, 周薇薇, 唐元昊, 张华, 施思齐. CeO2的电子结构,光学和晶格动力学性质:第一性原理研究. 物理学报, 2010, 59(5): 3426-3431. doi: 10.7498/aps.59.3426
    [17] 陈海川, 杨利君. LiGaX2(X=S, Se, Te)的电子结构,光学和弹性性质的第一性原理计算. 物理学报, 2011, 60(1): 014207. doi: 10.7498/aps.60.014207
    [18] 高攀, 柳清菊, 张学军. 氮铁共掺锐钛矿相TiO2电子结构和光学性质的第一性原理研究. 物理学报, 2010, 59(7): 4930-4938. doi: 10.7498/aps.59.4930
    [19] 焦照勇, 郭永亮, 牛毅君, 张现周. 缺陷黄铜矿结构Xga2S4 (X=Zn, Cd, Hg)晶体电子结构和光学性质的第一性原理研究. 物理学报, 2013, 62(7): 073101. doi: 10.7498/aps.62.073101
    [20] 王闯, 赵永红, 刘永. Ga1–xCrxSb (x = 0.25, 0.50, 0.75) 磁学和光学性质的第一性原理研究. 物理学报, 2019, 68(17): 176301. doi: 10.7498/aps.68.20182305
  • 引用本文:
    Citation:
计量
  • 文章访问数:  613
  • PDF下载量:  53
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-04-29
  • 修回日期:  2020-05-20
  • 上网日期:  2020-05-20
  • 刊出日期:  2020-08-20

一种新型二维TiO2的电子结构与光催化性质

  • 1. 湖北医药学院公共卫生与管理学院, 十堰 442000
  • 2. 湖北省南水北调水源区生物医药研发检测共享平台, 十堰 442000
  • 通信作者: 康文斌, wbkang@hbmu.edu.cn ; 方文玉, jzfangwenyu@163.com

摘要: 基于第一性原理计算方法, 设计出了一种新型二维半导体材料TiO2, 并进一步研究了其结构稳定性, 电子结构, 载流子迁移率和光学性质等. 二维TiO2的形成能、声子谱、分子动力学、弹性常数表明, 二维TiO2具有较好的动力学, 热力学和机械稳定性, 具备实验制备的条件, 且能够稳定存在于常温条件下. 电子结构分析表明, 二维TiO2是一种间接带隙半导体, 在GGA+PBE和HSE06算法下的能隙分别为1.19 eV和2.76 eV, 其价带顶和导带底能级分别由Ti-3d和Ti-4s态电子构成, O原子的电子态在费米能级附近贡献很小, 主要分布在深处能级. 载流子迁移率显示, 二维TiO2的迁移率比单层MoS2要小, 其电子和空穴迁移率分别为31.09和36.29 cm2·V–1·s–1. 由于空穴迁移率和电子迁移率的各向异性, 电子-空穴复合率较低, 使得单层TiO2的使用寿命更长, 光催化活性更好. 在应变调控下, 二维TiO2的能隙发生明显响应, 以适用于各种半导体器件的需要. 半导体的带边势和光学性质显示, 在–5%—2%单/双轴应变下, 二维TiO2能够光裂水制H2, 在–5%—5%单/双轴应变下, 能够光裂水制O2, H2O2和O3等. 此外, 二维TiO2对可见光和紫外光具有较高的吸收系数, 说明其在未来光电子器件和光催化材料领域有着潜在的应用前景.

English Abstract

    • 石墨烯自从2004年被成功制备以来[1], 便吸引了众多学者的兴趣. 随着研究的深入, 越来越多的新型二维材料被报道出来, 例如氧化石墨烯[2], h-BN[3], 过渡金属硫化物(transition metal dichalcogenides, TMDCs)[4]. 各族化合物也相继被报道, 如III族的硼烯[5]和铝烯[6], IV族的石墨烯[7]、硅烯[8]、锗烯[9]和锡烯[10], V族的黑磷烯[11]、砷烯[12]和锑烯[13]等, 这极大丰富了二维材料家族. 二维材料因其独特的结构[14-16], 如表面积大、载流子迁移路径短, 导致电子-空穴的复合率很低, 因而具有较好的光催化性能, 成为理想的光催化剂[17]. 例如, Zhu等[18]通过计算表明, 二维硅硫族化合物SiX(X = S, Se, Te)具有优异的光催化性质, 且通过双轴应变能够灵活调节材料的能隙和带边势, 是一种很有前途的水分解光催化剂. Sun等[19]通过机械剥离二维SnS2, 它在1.0 V电压下产生的光电流密度为2.75 mA·cm–2, 几乎是体相SnS2的72倍. 方文玉等[20]通过H, F修饰单层GeTe, 使得GeTe的载流子复合率降低, 且其光催化性能明显增强. 其他新型纳米结构光催化材料[21]也被陆续报道, 如二维MoS2[22], C3N4[23], 过渡金属硫化物(TMDCs)[24]、过渡金属碳化物或氮化物MXenes[25]、III族卤化物[26]和MnPSe3[27]也已在理论上预测用于光催化剂应用.

      锐钛矿相TiO2是一种常见的光催化材料, 其能隙约为3.23 eV[28], 在太阳光照射下能够光裂水制氢和析氧[29], 受到很多学者的青睐[30,31]. 关于二维TiO2的研究近期也有诸多报道[32], 比如Song等[33]研究TiO2纳米片对肝脏和脾脏毒性反应, 发现纳米结构TiO2是一种安全的生物材料. Liu等[34]以CdSxSe1-x为中间层, 在TiO2纳米管上制备出了多层阶梯式结构, 能够显著提升其光电电流强度和量子效率, 为光电化学水分解提供了一条全新的思路. Song等[35]通过椭圆偏振光谱的非接触无损检测方法, 证明二维TiO2是可以作为半导体器件中栅氧化物. Tao等[36]通过钛金属氧化方法在金红石TiO2(011)表面制备纳米结构TiO2, 能够将TiO2的能隙降低至2.1 eV, 其光催化活性显著增强.

      2017年, Xie等[37]提出了一种新型二维材料设计方法, 即在保持原晶体结构的前提下, 可以用相邻族元素替换的方式设计出新材料. 受此启发, Ti和Mo同属过渡金属, 且原子半径Mo(1.30 Å)和Ti(1.32 Å)、电负性Mo(3.9)和Ti(3.45)都比较接近. 而O和S是同族元素, 电子排列相似. 因此, 我们设想分别用Ti和O去替换单层MoS2中的Mo和S原子, 设计出一种与单层MoS2相同结构(P-6M2)的二维TiO2. 目前关于P-6M2 结构的单层TiO2尚未见报道, 本文基于第一性原理计算方法, 系统分析二维P-6M2 TiO2的稳定性、电子结构、载流子迁移率和光学性质等, 为实验上制备出新型光电子材料提供理论指导.

    • 计算所选取的二维TiO2与单层MoS2具有相同的结构, 其对称群为P-6M2. 为了避免层与层之间弱范德华力的影响, 沿z轴方向设置了20 Å的真空层, 晶胞结构如图1所示. 文中所有计算都是基于CASTEP[38]软件来完成的, 交换相关能用Perdew-Burke-Ernzerhof(PBE)泛函中的广义梯度近似(GGA)表示. 各原子选取的价电子组态分别为Ti-3p63d24s2, O-2s2sp4. 计算均在倒格矢中进行, 其中平面波截断能设置为Ecut = 450 eV, 采用20 × 20 × 1 Monkhorst-Pack 特殊K点对全布里渊区求和. 自洽场计算中, 能量的收敛精度设置为1.0 × 10–6 eV·atom–1, 内应力收敛精度0.001 GPa, 各原子受力不大于0.001 eV·nm–1, 最大位移收敛标准1.0 × 10–3 nm. 为了避免GGA+PBE算法导致能隙偏小的问题[39], 选取Heyd-Scuseria-Ernzerhof(HSE06)泛函杂化函数计算材料的能带结构. 此外, 选取模守恒赝势和有限差分方法计算声子谱, 采用NVT正则标准对5×5×1的超晶胞结构进行分子动力学模拟(AIMD), 计算各种体系时均进行结构优化.

      图  1  二维TiO2的晶体结构 (a) 俯视图; (b) 侧视图; (c) K点路径

      Figure 1.  Crystal structure of monolayer TiO2 (a) top view; (b) side view; (c) K point path.

    • 表1是二维TiO2自由弛豫后的结构常数和结合能Ef等, 其中a/b为晶格常数, ${\theta _1}$${\theta _2}$分别表示Ti—O—Ti和O—Ti—O的键角, l为Ti与O原子间的键长, $\sigma $为TiO2优化后单层的厚度buckling height, 所有参数均在图1中做了标注. 为了便于比较, 同时计算了单层MoS2各项结构参数, 各项参数与Castellanos-Gomez等[40]计算结果基本相符. 从表1可以看出, 二维TiO2的各项参数与单层MoS2很接近, 其中晶格常数a/b, xz平面内O—Ti—O原子间夹角${\theta _2}$、键长l以及buckling height $\sigma $均小于单层MoS2. 分析认为, 虽然Mo(1.30 Å)原子半径较Ti(1.32 Å)原子稍小, 但是S(1.04 Å)原子半径比O(0.66 Å)大得多, 导致MoS2的上述参数都比TiO2要大. 同时, 由于TiO2$\sigma $值和键长l更小, 从而导致其xy面内的原子夹角${\theta _1}$较MoS2更大.

      Materiala/b${\theta _1}$/($^ \circ $)${\theta _2}$/($^ \circ $)l$\sigma $/ÅEf/eV
      TiO22.8991.6668.162.012.26–8.11
      MoS23.1882.5880.742.413.14–7.35

      表 1  二维TiO2的结构常数和结合能

      Table 1.  Structure constants and binding energy of monolayer TiO2.

      对于一种新型材料, 实验的可行性是一个重要指标, 形成能能够表征材料实验合成的难易程度, 其计算公式为[39]

      ${E_{\rm{f}}} = \dfrac{{E_{\rm{t}}} - {E_{{\rm{Ti}}}} - 2{E_{\rm{O}}}}{3},$

      式中${E_{\rm{f}}}$表示形成能, ${E_{\rm{t}}}$表示优化后TiO2原胞的能量, ${E_{{\rm{Ti}}}}$${E_{\rm{O}}}$则分别表示单个Ti和O原子的能量. 从表1可以看出, 单层TiO2和MoS2的形成能分别为–8.11和–7.35 eV, 即单层TiO2的制备是一种放热反应, 且相对于单层MoS2更容易在实验室合成.

    • 通常情况下, 可以通过声子谱的软化程度来判断二维材料的动力学稳定性[41]. 图2是二维TiO2的声子谱计算结果, 由于二维TiO2原胞共有3个原子, 因此其总共有9支振动模式, 包括3支低频的声学模和6支高频的光学模. 从图2中可以看出, TiO2声子谱都没有虚频, 因此, 二维TiO2具有较好的动力学稳定性.

      图  2  二维TiO2的声子谱

      Figure 2.  Phonon spectrum of monolayer TiO2.

      为了进一步研究二维TiO2的热力学稳定性, 对TiO2的5 × 5 × 1超晶胞结构进行分子动力学模拟计算[42], 计算中选取NVT正则标准, 温度设置为500 K, 模拟时间5 ps, 步长2 fs. 计算结果如图3所示, 可以发现, 二维TiO2在加热过程中总能量在增加, 但是增加幅度非常小(约4.5 eV), 材料的温度一直在500 K上下起伏. 更重要的是, TiO2超晶胞在整个过程中没有出现明显的原子重构及键断裂现象, 这表明二维TiO2具有较好的热稳定性, 能够在常温下稳定存在.

      图  3  二维TiO2的分子动力学模拟

      Figure 3.  Molecular dynamics simulation of monolayer TiO2.

      此外, 对于一个稳定的二维材料也应该具有稳定的机械稳定性, 通常情况下可以用Born–Huang准则[43]来表征, 即需满足${C_{11}}{C_{22}} - C_{12}^2 > 0$${C_{66}} > 0$, 这里的${C_{ij}}$指的是弹性常数. 通过计算二维TiO2的弹性常数显示, ${C_{11}} = {C_{22}}$=77.18 GPa, ${C_{12}}$=8.92 GPa, ${C_{66}}$=34.13 GPa, 很显然满足Born–Huang准则. 综上分析, 二维TiO2同时具有较好的动力学、热力学和机械稳定性.

    • 作为一种新的二维材料, 研究其力学性质可以带来新的认识和应用的可能性. 为此, 根据二维TiO2的弹性常数进一步计算了其在xy平面内的杨氏模量$Y(\theta )$和泊松比$\upsilon (\theta )$($\theta $是相对于a轴的夹角), 计算公式满足[44]:

      $Y(\theta ) = \frac{{{C_{11}}{C_{12}} - C_{12}^2}}{{{C_{11}}{{\sin }^4}\theta + \left(\dfrac{{{C_{11}}{C_{12}} - C_{12}^2}}{{{C_{66}}}} - 2{C_{12}}\right){{\sin }^2}\theta {{\cos }^2}\theta + {C_{22}}{{\cos }^4}\theta }},$

      $\upsilon (\theta ) = \frac{{{C_{12}}{{\sin }^4}\theta - \left({C_{11}} + {C_{22}} - \dfrac{{{C_{11}}{C_{12}} - C_{12}^2}}{{{C_{66}}}}\right){{\sin }^2}\theta {{\cos }^2}\theta + {C_{12}}{{\cos }^4}\theta }}{{{C_{11}}{{\sin }^4}\theta + \left(\dfrac{{{C_{11}}{C_{12}} - C_{12}^2}}{{{C_{66}}}} - 2{C_{12}}\right){{\sin }^2}\theta {{\cos }^2}\theta + {C_{22}}{{\cos }^4}\theta }}.$

      力学性质计算结果如图4所示, 由于空间反演对称性, 二维TiO2的杨氏模量和泊松比表现出明显的各向异性, 但沿a轴和b轴方向上的力学性质相同, 且在a/b轴方向达到最大值, 在与a轴45°和135°夹角处取最小值. 对比常见二维材料的力学性质发现, 二维TiO2的杨氏模量最高可达76.15 GPa, 明显小于石墨烯(1.0 TPa)[45]、单层MoS2(300 GPa)[40]和氧化石墨烯(200 GPa)[46]的杨氏模量, 但大于黑磷烯(41.3 GPa)[11], 稍小于二维TnSe (101.37 GPa)[47]. 对于二维TiO2的泊松比, 最大可达0.116, 小于石墨烯(0.175)[48]、h-BN(0.211)[49]和单层MoS2(0.25)[50], 且大于8-Pmmn硼烯[41] (0.042—0.062)的泊松比. 因此, 二维TiO2是一种常见的正泊松比材料, 具有中等偏小的杨氏模量, 因此在柔性电子器件中具有潜在的应用前景.

      图  4  杨氏模量(左图/蓝色)和泊松比(右图/红色)极坐标图, $\theta $为相对于a轴方向的夹角

      Figure 4.  Polar coordinates for Young's modulus (left/blue line) and Poisson's ratio (right/red line) $\theta $ is the angle with respect to the a-direction.

    • 为考虑电子自旋极化(SOC)[51]对能带结构的影响, 分别计算了GGA+PBE算法加上自旋极化和不加自旋极化能带结构, 结果显示两种能带图几乎没有区别(能隙仅相差0.002 eV), 因此可以认为自旋极化对体系的计算结果几乎没有影响, 可以不用考虑. 为了比较这种P-6M2结构的TiO2与常见TiO2结构 (金红石相和锐钛矿相) 性质的不同, 分别计算出了单层金红石(rutile TiO2)和单层锐钛矿(anatase TiO2)的能带结构和分波态密度, 如图5所示. 图中分别描绘了GGA+PBE (红色)和HSE06 (蓝色)两种算法下的能带结构, 以及HSE06算法下的分波态密度. 观察能带图可以发现, 两种算法下的能级色散曲线的变化趋势基本一致. 通常情况下, GGA+PBE会低估半导体的能隙, 而HSE06计算的结果则能与实验值很好地符合[52]. 通过图5(a)可知, P-6M2 TiO2计算能隙分别为1.19 eV (GGA+PBE)和2.76 eV (HSE06), 因此2.76 eV可以视为二维TiO2的真实能隙, 这比锐钛矿TiO2的能隙(3.2 eV)[53]更小. 由于价带顶和导带底不在同一点, P-6M2 TiO2属于间接带隙半导体. 比对图5(b)图5(c)发现, 单层rutile TiO2和anatase TiO2的计算能隙分别高达5.02和4.70 eV (HSE06), 这明显高于体相rutile TiO2 (3.0 eV)和体相anatase TiO2 (3.2 eV)的能隙[54]. 单从能隙宽度来看, P-6M2 TiO2属于窄带隙半导体, 它能够对太阳光谱产生很好的响应. 而单层rutile TiO2和anatase TiO2能隙均大于4 eV (接近5 eV), 属于宽带隙半导体, 很难对太阳光谱产生响应. 因此, 相比较而言, P-6M2 TiO2更适用于光催化材料, 因此后文中将主要对单层P-6M2 TiO2进行分析. 对照分波态密度可知, 在–3—6 eV的能量范围内, 能级曲线主要由Ti原子外的电子贡献, O原子的贡献很少, 仅在能级深处(–2.5 eV)处有少量的分布. 此外, 半导体的价带顶主要由Ti-3 d组成, 导带底则是由Ti-4 s组成. 态密度在费米面上下出现两个孤立的峰值, 这与价带顶和导带底各自形成一条孤立的能级相符合; 而在价带深处(–3—–1.5 eV)和导带远处(4.5—6.0 eV)则出现相对较宽的峰值, 与之对应的则是能级曲线在以上两个范围内的局域化程度明显加剧, 并且发生能级简并现象. 由于半导体对光电性质起决定性作用的主要是费米面附近电子的分布, 半导体费米面上下出现峰值, 这将允许更多的载流子在价带顶和导带底发生跃迁而产生光吸收现象, 这显然对其光催化性质是有利的.

      图  5  能带图(左)和分波态密度图(右) (a) 单层P-6M2 TiO2; (b) 单层金红石TiO2; (c) 单层锐钛矿TiO2

      Figure 5.  Energy band structure (left) and density of states (right) of (a) monolayer P-6M2 TiO2, (b) monolayer rutile TiO2, (c) monolayer anatase TiO2.

      为研究晶体中原子间的成键特征以及电子转移情况, 选取Ti和O原子在同一平面内的差分电荷密度[55], 如图6所示. 图6中差分电荷密度线保持同一刻度, 刻度条的数值表示差分电荷密度等高线数值的大小, 其中蓝色表示失去电子, 红黄色表示得到电子, 颜色越深表示电子转移得越多. 很明显, Ti主要表现为失去电子, O则得到电子, 即电子由低电负性的Ti(3.45)转移至高电负性的O(7.54)原子, 因此, Ti—O之间主要以离子键形式存在.

      图  6  二维TiO2的差分电荷密度图

      Figure 6.  The differential charge density of monolayer TiO2.

      载流子迁移率是半导体材料的重要指标, 通过计算价带顶光生空穴和导带底光生电子的迁移率, 可以很直观地反映电子和空穴的复合率[56], 合理评价半导体的光催化活性. 对于二维半导体材料, 可以通过Bardeen和Shockley提出的形变势理论[57]来表达, 如果把二维材料看作均匀的介质, 其载流子迁移率的计算公式如下[58]:

      ${\mu _{{\rm{2D}}}} = \frac{{e{\hbar ^3}{C^{_{{\rm{2D}}}}}}}{{{k_{\rm{B}}}T{m^ * }m_l^{\rm{*}}{{(E_l^i)}^2}}},$

      式中, $e$为电子带电量, $\hbar $为约化普朗克常数, ${k_{\rm{B}}}$为玻尔兹曼常数, ${m^ * }$为载流子有效质量(可以根据公式${m^ * } = {\hbar ^2}{\left( {\dfrac{{{\partial ^2}{E_0}(k)}}{{\partial {k^2}}}} \right)^{ - 1}}$计算[59], ${E_0}(k)$为载流子迁移方向上价带顶或导带底能级), ${m_l}^ * $为平均有效质量(计算公式为$m_l^ * = \sqrt {m_a^ * m_b^ * } $, $m_a^ * $$m_b^ * $分别为沿a轴和b轴方向上的有效质量), T为温度(300 K), ${E_l}$为形变势[60](计算公式为$E_l^i = \dfrac{{\Delta {E_i}}}{{\Delta l/{l_0}}}$, 这里$\Delta {E_i}$为第i条能级在压缩或者拉伸应变下能量的变化, $\Delta l$为压缩或者拉伸长度, ${l_0}$为晶胞初始长度, 设置的应变范围为$\varepsilon $(–1%—1%), 应变步长为0.5%), ${C^{{\rm{2 D}}}}$为平面刚度, 计算公式为[61]

      ${C^{{\rm{2D}}}} = {{{{({\partial ^2}E} / {\partial {\varepsilon ^2}}})} / {{S_0}}},$

      式中$E - {E_0}$, ${S_0}$${{\Delta l}}/{{{l_0}}}$分别表示晶胞的总能量的变化, xy平面面积, 以及应变压缩或者拉伸比例. 图7 (a)为二维TiO2总能量随应变量的变化, 采取二次曲线拟合计算平面刚度; 图7(b)图7(c)表示价带顶(VBM)和导带底(CBM)相对真空能级的变化, 采取线性拟合计算材料的形变势, 计算中均是采取3 × 3 × 1超晶胞结构进行. 根据(4)式和(5)式可以计算出二维TiO2的载流子迁移率等有关性质, 如表2所列.

      Carrier type$m_a^ * $/${m_{\rm{e}}}$$m_b^ * $/${m_{\rm{e}}}$$m_l^ * $/${m_{\rm{e}}}$$\left| {{E_{la}}} \right|$/ eV$\left| {{E_{lb}}} \right|$/ eV$C_a^{2{\rm{D}}}$/ N·m–1$C_b^{{\rm{2 D}}}$/ N·m–1$\mu _a^{2{\rm{D}}}$/ cm2·V–1·s–1$\mu _b^{2{\rm{D}}}$/ cm2·V–1·s–1
      Electrons3.211.392.113.433.3821.2721.2812.9230.75
      Holes4.734.124.411.261.2521.2721.2831.0936.29

      表 2  二维TiO2有效质量${m^ * }$, 形变势常数$E_{_{\rm{d}}}^i$, 弹性常数${C^{{\rm{2 D}}}}$和载流子迁移率${\mu _{{\rm{2 D}}}}$

      Table 2.  Calculated effective mass ${m^ * }$, deformation potential constant $E_{_{\rm{d}}}^i$, elastic modulus ${C^{{\rm{2 D}}}}$, and carrier mobility ${\mu _{{\rm{2 D}}}}$ for monolayer TiO2 along the a ($\zeta \to K$) and b ($\zeta \to G$) directions, where $\zeta $ represents the position of the valence band top and the conduction band bottom.

      图  7  (a) 二维TiO2沿a/b方向的总能量与应变量$\Delta l/l$的关系, 采用二次数据拟合二维结构的平面刚度, 黑色和红色曲线表示沿ab方向的面内刚度; (b), (c)单层TiO2的VBM和CBM随应变量相对真空能级的变化, 采取线性拟合计算形变势

      Figure 7.  (a) The relation between total energy and the applied strain $\Delta l/l$ along the a/b directions of monolayer TiO2. The quadratic data fitting gives the in-plane stiffness of 2D structures. Black and red curves show the in-plane stiffness along the a and b directions of monolayer TiO2. The shift of VBMs and CBMs for (b-c) monolayer TiO2 with respect to the vacuum energy, as a function of the applied strain along either the a and b direction. The linear fit of the data yields the deformation potential constant.

      图7可以看出, 二维TiO2沿ab两个方向上的总能量、CBM和VBM变化趋势非常相似, 因而两个方向上的平面刚度(21.27, 21.28 N·m–1)和形变势(3.42, 3.38 eV)差异很小, 如表2所列. 计算载流子有效质量发现, 二维TiO2空穴质量明显大于电子质量, 这说明价带顶附近能级较导带底附近能级更加平坦, 对照图5也证实了这一点. 根据(4)式可知, 正是有效质量的不同导致载流子迁移率表现出明显的各向异性, 沿a, b方向上的电子迁移率为19.92和31.09 cm2·V–1·s–1, 空穴迁移率分别为30.75和36.29 cm2·V–1·s–1, 这比单层MoS2的迁移率(200 cm2·V–1·s–1)[62]要小. 这说明二维TiO2载流子迁移率相对较小, 并且电子和空穴的迁移率表现出明显差异, 即两者的分离效率和异步化程度较高, 导致电子和空穴的复合率较低[63], 这样的半导体更耐用, 光催化活性更好.

      为进一步研究弹性形变对二维TiO2电子结构的影响, 分别计算了在–5%—5%的面内压缩/拉伸应变下TiO2的能隙变化, 计算选取HSE06函数, 结果如图8所示. 很明显, 二维TiO2的能隙对弹性应变有着明显响应, 其能隙随压缩/拉伸应变单调增大/减小, 并且在沿a/b轴双轴方向应变下能隙变化得更快, 这说明二维TiO2的能隙可以通过形变进行调控, 以适用于各种电子器件的需求.

      图  8  单/双轴应变下能隙变化

      Figure 8.  Band gap of monolayer TiO2 under uniaxial/biaxial strain, calculated using the HSE06 functional.

    • 前文分析结果表明二维TiO2属于间接窄带隙半导体, 且通过单/双轴应变可以有效调控其能隙, 为此进一步分析了二维TiO2在不同应变下的光催化解水能力. 一般情况下, 参照标准氢电极的氧化还原势, 半导体的氧化还原电势为[56]:

      $CB = \chi - {E_{\rm e}} - 0.5{E_g},$

      $VB = {E_{\rm c}} + {E_g},$

      其中CBVB分别表示导带底与价带顶的氧化还原电势, $\chi $表示材料的电负性, 它可以通过计算各原子密立根电负性几何平均值得到. Ee为以氢为标准时自由电子的电势(4.50 eV), Eg为半导体能隙. 已知Ti和O原子的电负性分别为3.45和7.54, 可导TiO2电负性$\chi $为5.81. 通过(6)式和(7)式计算得到二维TiO2在不同应变下的氧化还原势, 并将其与光催化水的部分产物的电势作比较, 如图9所示.可以看出, TiO2的导带底的电势随着压缩/拉伸形变的增加而减小/增加, 价带顶的电势则随着压缩/拉伸形变的增加而增加/减小, 且在–5%—2%的单/双轴应变下其电势小于0 eV, 即低于光裂水制氢的还原势(0.0 eV), 因此在这个应变范围内, 二维TiO2能够光裂水制H2. 值得注意的是, 在整个应变范围内(–5%—5%), TiO2价带顶的电势很大(大于2.5 eV), 具有很强的氧化性, 这说明二维TiO2还能够光催化水析出O2, H2O2和O3等, 而这些产物都有较强的氧化能力, 可以直接氧化很多污染物, 这使得TiO2有望成为一种新型二维光催化材料.

      图  9  单双/轴应变下单层TiO2光催化示意图

      Figure 9.  Schematic diagram of monolayer TiO2 photocatalysis under uniaxial/biaxial strain.

      采用HSE06泛函杂化函数计算二维TiO2的介电函数, 可以根据Kramers-Kronig的色散关系和直接跃迁概率的定义推导出其吸收系数等光学常数, 计算关系可以表达为[64]:

      $\varepsilon (\omega ) = {\varepsilon _1}(\omega ) + {\rm{i}}{\varepsilon _2}(\omega ),$

      $\alpha (\omega ) = \sqrt 2 \omega {[{({\varepsilon _1}^2(\omega ) + {\varepsilon _2}^2(\omega ))^{1/2}} - {\varepsilon _1}(\omega )]^{1/2}},$

      其中$\varepsilon (\omega )$表示介电函数, ${\varepsilon _1}(\omega )$${\varepsilon _2}(\omega )$分别表示实部和虚部, $\omega $是频率, $\alpha (\omega )$为光吸收系数. 根据(8)式和(9)式得到介电函数的虚部和吸收系数随波长的变化, 如图10所示, 这里分别选取面内沿x轴方向xx (蓝色)和沿z轴方向zz (红色)两个光线入射方向. 从图中可以看出, 二维TiO2对光谱的响应范围内(100—700 nm), ${\varepsilon _2}(\omega )$$\alpha (\omega )$的变化趋势保持一致, 在200 nm附近出现峰值, 在400 nm以上的分布很少, 这说明其对中长波紫外线及短段可见光具有一定的吸收效果. 值得注意的是, 吸收系数对应的吸收边约为480 nm, 对应的光子能量约为2.58 eV, 与TiO2的能隙(2.76 eV)比较接近. 这主要是由于吸收边与光学带隙直接相关, 当光子能量达到能隙范围, 半导体便能产生光-电子耦合效应, 激发电子在价带与导带间跃迁. 此外, 吸收系数达到了105的数量级, 与硅的吸收系数[54]相当, 说明二维TiO2作为一种高效的光吸收材料有望应用于太阳能电池和其他光电器件中.

      图  10  光学性质 (a) 介电函数虚部; (b) 光吸收系数

      Figure 10.  Optical properties: (a) virtual part of dielectric function; (b) absorption coefficient.

    • 基于第一性原理计算方法, 设计出了一种新型二维TiO2, 通过计算其结合能、声子谱、分子动力学和弹性常数表明, 二维TiO2具备实验合成的条件, 同时拥有较好的动力学、热力学和机械稳定性, 能够在常温条件下稳定存在. 电子结构分析表明, 二维TiO2的计算能隙为2.67 eV, 为一种间接带隙半导体. 其迁移率比单层MoS2要小, 电子和空穴迁移率分别为31.09和36.29 cm2·V–1·s–1 即二维TiO2的电子和空穴的复合率较低, 具备较好的光催化活性. 在单/双轴应变下, 二维TiO2能隙随压缩/拉伸应变量单调增大/减小, 即可以通过外力应变有效调控其能隙, 以适应不同电子器件的需求. 分析其光催化性质表明, 二维TiO2具备较好的光催化效果, 当施加–5%—2%的单/双轴应变时, 其还原势低于0 eV, 能够光裂水制H2, 而在–5%—5%的应变范围内, 二维TiO2的氧化势高于2.5 eV, 能够光裂水析O2, H2O2和O3等. 另一方面, 二维TiO2的吸收系数达到105的数量级, 光吸收范围覆盖中长波紫外线及可见光波段, 因此二维TiO2作为一种高效的光吸收材料有望应用于太阳能电池和其他光电器件中.

      感谢国家超级计算深圳中心计算资源的支持.

参考文献 (64)

目录

    /

    返回文章
    返回