搜索

x

留言板

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

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

基于形变势理论的掺杂计算Sb2Se3空穴迁移率

张冷 张鹏展 刘飞 李方政 罗毅 侯纪伟 吴孔平

张冷, 张鹏展, 刘飞, 李方政, 罗毅, 侯纪伟, 吴孔平. 基于形变势理论的掺杂计算Sb2Se3空穴迁移率. 物理学报, 2024, 73(4): 047101. doi: 10.7498/aps.73.20231406
引用本文: 张冷, 张鹏展, 刘飞, 李方政, 罗毅, 侯纪伟, 吴孔平. 基于形变势理论的掺杂计算Sb2Se3空穴迁移率. 物理学报, 2024, 73(4): 047101. doi: 10.7498/aps.73.20231406
Zhang Leng, Zhang Peng-Zhan, Liu Fei, Li Fang-Zheng, Luo Yi, Hou Ji-Wei, Wu Kong-Ping. Carrier mobility in doped Sb2Se3 based on deformation potential theory. Acta Phys. Sin., 2024, 73(4): 047101. doi: 10.7498/aps.73.20231406
Citation: Zhang Leng, Zhang Peng-Zhan, Liu Fei, Li Fang-Zheng, Luo Yi, Hou Ji-Wei, Wu Kong-Ping. Carrier mobility in doped Sb2Se3 based on deformation potential theory. Acta Phys. Sin., 2024, 73(4): 047101. doi: 10.7498/aps.73.20231406

基于形变势理论的掺杂计算Sb2Se3空穴迁移率

张冷, 张鹏展, 刘飞, 李方政, 罗毅, 侯纪伟, 吴孔平

Carrier mobility in doped Sb2Se3 based on deformation potential theory

Zhang Leng, Zhang Peng-Zhan, Liu Fei, Li Fang-Zheng, Luo Yi, Hou Ji-Wei, Wu Kong-Ping
PDF
HTML
导出引用
  • 硒化锑(Sb2Se3)是一种元素丰富、经济且无毒的太阳电池吸收层材料. 太阳电池的性能在很大程度上取决于载流子的传输特性, 然而在Sb2Se3中, 这些特性尚未得到很好的理解. 通过密度泛函理论和形变势理论, 本文对纯Sb2Se3以及掺杂了As, Bi的Sb2Se3的空穴传输特性进行研究, 计算并分析了影响迁移率的3个关键参数: 有效质量、形变势和弹性常数. 结果显示, 有效质量对迁移率具有最大影响, 掺杂Bi的Sb2Se3表现出最高的平均迁移率. 同时发现, Sb2Se3的空穴迁移率呈现出明显的各向异性, 其中x方向的迁移率远高于y, z方向, 这应该与x方向的原子主要以较强的共价键连接, 而y, z方向以较弱的范德瓦耳斯力连接有关. 载流子传输能力强的方向有助于有效传输和收集光生载流子, 本研究从理论上强调了控制Sb2Se3沿特定方向生长的重要性.
    Antimony selenide (Sb2Se3) is an element-rich, cost-effective, and non-toxic material used as an absorber layer in solar cells. The performance of solar cells is significantly influenced by the transport characteristics of charge carriers. However, these characteristics in Sb2Se3 have not been well understood. In this work, through density functional theory and deformation potential theory, we investigate the hole transport properties of pure Sb2Se3 and As-, Bi-doped Sb2Se3. The incorporation of as element and Bi element does not introduce additional impurity levels within the band gap. However, the band gaps are reduced in both As-Sb2Se3 and Bi-Sb2Se3 due to the band shifts of energy levels. This phenomenon is primarily attributed to the interactions between the unoccupied 4p and 6p states of the doping atoms and the unoccupied 4p states of Se atoms, as well as the unoccupied 5p states of Sb atoms. In this study, we calculate and analyze three key parameters affecting mobility: effective mass, deformation potential, and elastic constants. The results indicate that effective mass has the greatest influence on mobility, with Bi-Sb2Se3 exhibiting the highest average mobility. The average effective mass is highest in As-Sb2Se3 and lowest in Bi-Sb2Se3. The elastic constants of the As- and Bi-doped Sb2Se3 structures show minimal differences compared with that of the intrinsic Sb2Se3 structure. By comparing the intrinsic, As-doped, and Bi-doped Sb2Se3, it is evident that doping has a minor influence on deformation potential energy along various directions. The study reveals that the hole mobility in Sb2Se3 displays significant anisotropy, with higher mobilities observed in the x-direction and the y-direction than in the z-direction. This discrepancy is attributed to stronger covalent bonding primarily in the x- and y-direction, while in the z-direction weaker van der Waals forces is dominant. The directions with enhanced charge carrier transport capability contribute to efficient transfer and collection of photo-generated charge carriers. Therefore, our research theoretically underscores the significance of controlling the growth of antimony selenide along specific directions.
      PACS:
      71.15.Mb(Density functional theory, local density approximation, gradient and other corrections)
      51.50.+v(Electrical properties)
      通信作者: 吴孔平, kpwu@jit.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 61904071, 52002170)、江苏省高校青蓝工程和金陵科技学院高层次人才科研启动金(批准号: jit-rcyj-202001)资助的课题.
      Corresponding author: Wu Kong-Ping, kpwu@jit.edu.cn
    • Funds: Project supported by the National Natural Science Foundation of China (Grant Nos. 61904071, 52002170), the Qing Lan Project of Jiangsu Provincial University of China, and the Jinling Institute of Technology High-level Talent Research Startup Fund (Grant No. jit-rcyj-202001).

    硒化锑(Sb2Se3)材料具有禁带宽度合适、吸收系数大、物相简单、原材料便宜无毒的优势, 非常有潜力制备高效率、低成本的薄膜太阳电池[13]. 在过去的十年中, Sb2Se3太阳能电池的转换效率取得了显著提升, 经历了从不到1%到10.57%[4,5]的快速增长过程. 然而, 与CdTe, Cu(InGa)Se2或混合卤化物钙钛矿器件高于22%的效率相比, Sb2Se3太阳能电池的效率仍有进一步提升的空间. 限制Sb2Se3电池性能的关键因素之一是其较弱的导电性. 导电性主要与载流子浓度和迁移率相关. 迁移率反映了载流子在电场作用下移动的速率, 是衡量半导体导电能力的重要物理指标, 也决定了太阳能电池等电子器件的工作效率. Sb2Se3材料呈带状结构, 由许多一维的(Sb4Se6)n纳米带通过范德瓦耳斯力堆叠而成. 正由于这种结构, Sb2Se3的载流子迁移率呈现出显著的各向异性. 目前, 在实验中很难直接测量载流子沿不同方向的迁移率, 因此只能测量平均迁移率.

    载流子传输能力低的方向会影响光生载流子的传输和收集. 通过调节工艺条件, 将pn结的内建电场方向与载流子传输能力强的方向一致, 有望显著提高载流子的输运效率. 因此, 对Sb2Se3的迁移率进行系统的理论计算具有重要意义. Wang等[6]采用第一性原理密度泛函理论以及玻尔兹曼输运方程, 计算出Sb2Se3在低缺陷浓度情况下, 平均空穴迁移率约为9.55 cm2/(V·s). 然而, 该结果与实验值存在一定差距, 需要进一步的研究来进行验证.

    由于本征Sb2Se3的弱导电性, 许多研究人员通过掺杂来提升其导电性[79]. 这些研究主要关注载流子浓度的提高以增加导电性, 对迁移率的讨论较少, 有关掺杂对Sb2Se3载流子迁移率影响的报道也很少. 又因为Sb2Se3应用在光伏领域以p型导电为主, 因此, 本文重点考察本征态、掺杂砷以及掺杂铋的Sb2Se3x, y, z三个方向上的载流子迁移率. 具体而言, 首先进行标准的密度泛函理论(DFT)计算, 以获得必要的参数, 如总能量、有效质量、形变势和弹性常数等. 接下来, 在室温条件下(T = 300 K), 计算了各个方向上的空穴迁移率, 并对影响输运性质的关键参数进行了讨论.

    本研究采用VASP (Vienna ab initio simulation package)中的第一性原理计算方法进行探究. 我们通过平面波基组展开电子波函数, 并使用广义梯度近似(GGA)下的PBE交换泛函来描述电子与电子之间的交换关联作用[10]. 为了确保计算的准确性, 我们进行了收敛性测试. 通过选择平面波截断能量为400 eV, 使用4 × 3 × 3的k空间取样密度, 能够保证总能量的收敛性. 迭代过程中, 收敛精度选择1 × 10–6 eV/atom. 在布里渊区中, 选取如下特殊K[11], 分别是G(0, 0, 0); X(0.5, 0, 0); Z(0, 0, 0.5); U(0.5, 0, 0.5); Y(0, 0.5, 0); S(0.5, 0.5, 0); T(0, 0.5, 0.5); R(0.5, 0.5, 0.5), 能带结构中的布里渊区路径为G-X-Z-U-Y-S-T-R.

    研究模型结构如图1所示, 总共包含40个原子, 其中包括24个硒原子和16个锑原子. 通过将金属原子(As, Bi)取代锑原子, 构建了金属掺杂的Sb2Se3模型. 在这个计算模型中, 锑原子存在于两个不同的位置, 分别标记为Sb1和Sb2; 而硒原子则分布在3个不同的位置, 标记为Se1, Se2和Se3. 在计算过程中, 我们观察到相比于替换Se原子位置和间隙位置, As原子替代Sb2位置时的形成能最低, Bi原子替代Sb1位置时的形成能最低, 因此本文研究分别选择将As, Bi原子置于Sb2, Sb1位置进行掺杂.

    图 1 Sb2Se3的晶胞\r\nFig. 1. Crystal structure of Sb2Se3 computational model.
    图 1  Sb2Se3的晶胞
    Fig. 1.  Crystal structure of Sb2Se3 computational model.

    根据形变势理论, 结合玻尔兹曼输运理论下的弛豫时间近似和DFT能带结构理论的有效质量近似[12,13], 可得三维材料体系的迁移率为

    $$ {\mu }_{\alpha }^{3{\mathrm{D}}}=\frac{2\sqrt{2{\mathrm{\pi }}}e{{\mathrm{\hslash }}}^{4}{C}_{\alpha }^{3{\mathrm{D}}}}{3{\left({k}_{{\mathrm{B}}}T\right)}^{3/2}{E}_{\alpha }^{2}{m}_{\alpha }^{*5/2}}, $$ (1)

    其中, e代表单电荷电量; $\hbar $代表约化普朗克常数; $ {C}_{\alpha }^{3{\mathrm{D}}} $表示沿载流子输运方向α的弹性常数; kB代表玻尔兹曼常数; T代表温度; Eα是形变势能, m*代表载流子有效质量. 为了直观表示, 本文只计算x, y, z三个方向上的载流子迁移率. 这3个方向上分别对应的弹性常数标记为C11, C22, C33, 对应的形变势标记为E11, E22, E33, 对应的有效质量标记为$ {m}_{xx}^{*} $, $ {m}_{yy}^{*} $, $ {m}_{zz}^{*} $. m*的计算方法为

    $$ {\left(\frac{1}{{m}^{*}}\right)}_{ij}=\frac{1}{{\hslash }^{2}}\frac{{\partial }^{2}E\left(\boldsymbol{k}\right)}{\partial {k}_{i}\partial{k}_{j}},\;\;\;i,\;j=x,\;y,\;z, $$ (2)

    其中, ij分别表示在倒空间的分量, $ E\left(\boldsymbol{k}\right) $是导带底或价带顶附近的能量-波矢关系.

    $ {C}_{\alpha }^{3{\mathrm{D}}} $的计算方法为

    $$ {{C}_{\alpha }^{3{\mathrm{D}}}= \frac{1}{{V}_{0}}\frac{{\partial }^{2}E}{{\partial \left({\text{δ}} l/{l}_{0}\right)}^{2}}\bigg|}_{l={l}_{0}}, $$ (3)

    其中, V0是完全优化后的晶胞体积, l0是沿α方向的晶格参数; δl表示沿l方向的形变; E代表晶胞总能量. Eα的计算方法为

    $$ {E}_{\alpha }=\frac{\partial {E}_{{\mathrm{e}}{\mathrm{d}}{\mathrm{g}}{\mathrm{e}}}}{\partial ({\text{δ}} l/{l}_{0})}, $$ (4)

    其中, $\partial $Eedge是由晶格形变引起的导带底或价带顶的能量偏移.

    平均有效质量$ {\overline{m}}^{*} $可以用(5)式来表示:

    $$ \frac{1}{{\overline{m}}^{*}}=\frac{1}{3}\left(\frac{1}{{m}_{xx}^{*}}+\frac{1}{{m}_{yy}^{*}}+\frac{1}{{m}_{zz}^{*}}\right). $$ (5)

    根据Matthiessen规则, 总迁移率可用(6)式来表示, 其中μtotal是总迁移率, μλ是分模迁移率:

    $$ \frac{1}{{\mu }_{{\mathrm{t}}{\mathrm{o}}{\mathrm{t}}{\mathrm{a}}{\mathrm{l}}}}=\sum _{\lambda }\frac{1}{{\mu }_{\lambda }}. $$ (6)

    本研究采用第一性原理计算来获取不同模型的电子结构、载流子有效质量、形变势能和弹性常数. 在进行能量和性质计算之前, 先对晶体结构进行了优化, 优化过程考虑了范德瓦耳斯修正. 表1列出了优化后的晶格参数, 可以发现, 本征Sb2Se3的晶格参数与文献[14, 15]报道实验值相差不大, 确保了晶体结构的可靠性. 此外, As, Bi原子尺寸分别小于、大于Sb原子, 因此相对应的掺杂了的晶格参数与体积分别小于、大于本征Sb2Se3.

    表 1  优化后的晶格参数
    Table 1.  Optimized lattice parameters.
    文献[14] 文献[15] Sb2Se3 AsSb2Se3 BiSb2Se3
    a 3.98 3.99 3.99 3.97 4.00
    b 11.64 11.63 11.35 11.35 11.36
    c 11.79 11.78 11.65 11.63 11.65
    v3 547.1 546.6 527.3 524.7 528.9
    下载: 导出CSV 
    | 显示表格

    为了深入探讨纯Sb2Se3和掺杂Sb2Se3之间的电子结构、有效质量以及形变势之间的关系, 进行了能带结构和分态密度(PDOS)的计算. 图2(a)展示了未掺杂的Sb2Se3的能带结构和PDOS, 其中费米能级以虚线表示. 可以看出, 计算得到的纯Sb2Se3的带隙为0.83 eV, 价带顶(VBM)和导带底(CBM)主要由Se 4p态和Sb 5p态贡献. Se 4p态在费米能级附近出现局域化, 导致价带比导带更为平坦, 这意味着纯Sb2Se3具有较大的空穴有效质量和较小的迁移率. 这一电子性质与已报道的理论结果相近[11], 初步验证了我们所采用的计算方法的合理性.

    图 2 纯Sb2Se3 (a), As-Sb2Se3 (b), Bi-Sb2Se3 (c)的能带结构及分态密度图. 虚线处是费米能级, 蓝圈处是价带顶(VBM)\r\nFig. 2. Band structures and partial density of states of pure Sb2Se3 (a), As-Sb2Se3 (b), and Bi-Sb2Se3 (c). Dashed line represents the Fermi level and the blue circle represents the valence band maximum (VBM).
    图 2  纯Sb2Se3 (a), As-Sb2Se3 (b), Bi-Sb2Se3 (c)的能带结构及分态密度图. 虚线处是费米能级, 蓝圈处是价带顶(VBM)
    Fig. 2.  Band structures and partial density of states of pure Sb2Se3 (a), As-Sb2Se3 (b), and Bi-Sb2Se3 (c). Dashed line represents the Fermi level and the blue circle represents the valence band maximum (VBM).

    图2(b), (c)显示了As, Bi掺杂的Sb2Se3的能带结构和PDOS. 从能带结构看, 这两种与Sb等电子的元素的掺入并没有在禁带中引入额外的杂质能级. 尽管掺杂产生了一定影响, 但对能带结构的影响相对较小. 具体而言, Bi的掺杂对能带结构的影响要大于As的影响, As-Sb2Se3结构的带隙与为0.82 eV, Bi-Sb2Se3结构的带隙为0.81 eV, 3种结构的VBM均在(0.1111, 0, 0.3889). 此外, 从PDOS的角度来看, As和Bi原子的4p和6p态与导带附近的Sb 5p态之间存在更强的耦合效应. 同样地, As和Bi原子的4p和6p态与价带附近的Se 4p态之间也存在一定的耦合. As-Sb2Se3和Bi-Sb2Se3中能带移动导致的能隙减小, 主要归因于掺杂原子未占据的4p和6p态与Se原子未占据的4p态和Sb原子未占据的5p态之间的相互作用.

    3.3.1   空穴有效质量

    较小的载流子有效质量(<0.5m0, m0是自由电子质量)有利于载流子的输运, 对制备半导体器件很有利. 纯Sb2Se3, As-Sb2Se3与Bi-Sb2Se3的VBM在布里渊区的位置均是(0.1111, 0.000, 0.3889), 采用抛物线拟合方法, 计算了VBM沿着3个方向[100], [010]和[001]的空穴有效质量[16], 结果见表2, 分别由$ {m}_{xx}^{*} $, $ {m}_{yy}^{*} $, $ {m}_{zz}^{*} $表示. 可以发现, 空穴有效质量具有较强的各向异性, 本征Sb2Se3沿着x方向的空穴有效质量较小, 为0.43m0, 这对制备p型半导体有利. 另外, 纯Sb2Se3的有效质量计算结果与Zhang等[17]所计算的0.53m0, 0.42m0, 0.83m0与Wang等[18]所计算的0.85m0, 0.55m0与3m0趋势相似, 表明本计算方法基本可靠.

    表 2  不同Sb2Se3结构的空穴有效质量
    Table 2.  Effective mass of holes for different structures of Sb2Se3.
    m*/m0Sb2Se3As-Sb2Se3Bi-Sb2Se3
    $ {m}_{xx}^{*} $0.430.450.42
    $ {m}_{yy}^{*} $0.880.890.67
    $ {m}_{zz}^{*} $1.081.221.61
    $ {\stackrel{-}{m}}^{*} $0.680.720.67
    下载: 导出CSV 
    | 显示表格

    As掺杂的Sb2Se3沿着各方向的有效质量均比纯Sb2Se3的数值略大, 而Bi掺杂的Sb2Se3沿着x, y方向的有效质量比纯Sb2Se3的小, z方向的有效质量比纯Sb2Se3的大. 3种结构中, As-Sb2Se3的平均有效质量最大, Bi-Sb2Se3的平均有效质量最小.

    3.3.2   形变势与弹性常数

    为了计算形变势能和弹性常数, 对3种结构沿不同方向施加了从–1.5%—1.5%的应变, 步长为0.5%. 弹性常数计算中, 如(3)式所示, V0是完全优化结构的体积, 另一部分$ \dfrac{{\partial }^{2}E}{{\partial \left({\text{δ}} l/{l}_{0}\right)}^{2}} $是体系总能量和所施加应变的二阶导数. 利用VASPKIT[19]图3(a)中总能量和形变的抛物线拟合计算出弹性常数, 见图3(b). 本征Sb2Se3沿[100], [010], [001]方向的弹性常数为85.2, 61.1, 37.6 GPa, 与Silva等[20]算出的81.65, 55.20, 30.9 GPa相近. 3个方向的弹性常数相差较大, 反映了力学性质上的各向异性. As和Bi掺杂了的Sb2Se3结构的弹性常数与本征Sb2Se3结构的数值差别不大.

    图 3 (a)本征Sb2Se3总能量与形变的抛物线拟合, 其决定了弹性常数; (b)不同结构的弹性常数\r\nFig. 3. (a) Parabolic fitting of the intrinsic Sb2Se3 total energy and deformations, determining the elastic constants; (b) calculated elastic constants of different structures.
    图 3  (a)本征Sb2Se3总能量与形变的抛物线拟合, 其决定了弹性常数; (b)不同结构的弹性常数
    Fig. 3.  (a) Parabolic fitting of the intrinsic Sb2Se3 total energy and deformations, determining the elastic constants; (b) calculated elastic constants of different structures.

    形变势是通过在不同方向施加小应变时, 在VBM或CBM线性拟合能量偏移来计算的, 如图4(a)所示. 形变势计算如(4)式所示, 图4(b)呈现了沿[100], [010], [001]方向计算出的3种结构的形变势数值. 3种结构的形变势特点相似, E11E22数值接近, E33的值较小. 通过对比本征、As-Sb2Se3和Bi-Sb2Se3可知, 掺杂对于各个方向的形变势能影响不大.

    图 4 (a)本征Sb2Se3在3个不同方向应变下的价带边缘位置, 实线是线性拟合, 其决定了形变势; (b)不同结构的形变势能\r\nFig. 4. (a) Valence band edge positions of intrinsic Sb2Se3 under strain along three different directions, solid lines represent linear fitting, determining the deformation potential; (b) calculated deformation potential energies of different structures.
    图 4  (a)本征Sb2Se3在3个不同方向应变下的价带边缘位置, 实线是线性拟合, 其决定了形变势; (b)不同结构的形变势能
    Fig. 4.  (a) Valence band edge positions of intrinsic Sb2Se3 under strain along three different directions, solid lines represent linear fitting, determining the deformation potential; (b) calculated deformation potential energies of different structures.
    3.3.3   空穴迁移率

    结合以上空穴有效质量、弹性常数、形变势, 根据(1)式中的关系计算出了空穴迁移率. 表3呈现了不同模型结构的空穴迁移率值. 迁移率主要由3个变量确定, 不同结构下, 弹性常数和形变势之间仅存在微小差异, 可见这两项参数不是影响迁移率的决定性因素. 另一方面, 各个结构、各个方向的空穴有效质量有较大的差异, 且体现在公式上, 有效质量的幂指数是2.5, 因此, 有效质量是影响最终结果的主要因素. 在所研究的3个方向中, 迁移率沿[001]方向最小, 这是因为[001]方向以较弱的范德瓦耳斯力结合为主, 载流子传输时候阻碍较大. 文献[9, 2123]中实验方法测得本征Sb2Se3的平均空穴迁移率约在5.1—45 cm2/(V·s), 本文计算的平均空穴迁移率结果略大. 分析原因可能有两种: 测试样品中存在杂质和缺陷, 因此测试值低估了纯本征Sb2Se3的迁移率; 传统的形变势理论只包括由于纵向声学声子引起的晶胞变形引起的散射, 没有考虑载流子和声子以及载流子之间的相互作用等[24]. 根据(6)式可知, 载流子的整体迁移率小于单一机制引起的迁移率, 因此本文计算值偏大. 此外, 掺As的Sb2Se3的平均迁移率最低, 掺Bi的Sb2Se3的平均迁移率最高; Bi掺杂后, [010]方向的迁移率提高至掺杂前的2倍.

    表 3  三种结构的空穴迁移率
    Table 3.  Hole mobility of Sb2Se3, As-Sb2Se3 and Bi-Sb2Se3 along three principle directions.
    迁移率/(cm2·V–1·s–1)Sb2Se3As-Sb2Se3Bi-Sb2Se3
    μx232.62221.59240.66
    μy32.7831.2066.71
    μz20.0215.118.04
    μavg95.1489.30105.13
    下载: 导出CSV 
    | 显示表格

    通过密度泛函理论和形变势理论研究了纯Sb2Se3和As, Bi掺杂Sb2Se3的晶格参数、电子结构和空穴输运特性. 研究表明, 相较于形变势能和弹性常数, 有效质量对于迁移率的影响最强. Sb2Se3的空穴迁移率具有较强的各向异性, x方向的迁移率远远大于y, z方向的迁移率, 这与[001]方向主要以较弱的范德瓦耳斯力结合为主有关. 载流子传输能力低的方向不利于传输和收集光生载流子, 从理论上反映了Sb2Se3可控方向生长的重要性. 此外, 在本文研究的体系下, Bi-Sb2Se3的平均迁移率最高. 未来的研究需要进一步修正形变势理论, 将其他散射机制考虑入内, 以更准确地计算各项参数对迁移率的影响.

    [1]

    Chen C, Li K H, Tang J 2022 Sol. RRL 6 2200094Google Scholar

    [2]

    Zhang X, Li C, Sun K, Zhou J, Zhang Z 2021 Adv. Energy Mater. 11 2002614Google Scholar

    [3]

    薛丁江, 石杭杰, 唐江 2015 物理学报 64 038406Google Scholar

    Xue D J, Shi H J, Tang J 2015 Acta Phys. Sin. 64 038406Google Scholar

    [4]

    Zhao Y, Wang S, Li C, Che B, Chen X, Chen H, Tang R, Wang X, Chen G, Wang T, Gong J, Chen T, Xiao X 2022 Energy Environ. Sci. 15 5118Google Scholar

    [5]

    Li Z, Liang X, Li G, Liu H, Zhang H, Guo J, Chen J, Shen K, San X, Yu W, Schropp R, Mai Y 2019 Nat. Commun. 10 125Google Scholar

    [6]

    Wang X, Ganose A M, Kavanagh S R, Walsh A 2022 ACS Energy Lett. 7 2954Google Scholar

    [7]

    Spaggiari G, Bersani D, Calestani D, Gilioli E, Gombia E, Mezzadri F, Casappa M, Pattini F, Trevisi G, Rampino S 2022 Int. J. Mol. Sci. 23 15529Google Scholar

    [8]

    Huang M, Lu S, Li K, Lu Y, Chen C, Tang J, Chen S 2022 Sol. RRL 6 2100730Google Scholar

    [9]

    Liang G, Chen X, Ren D, Jiang X, Tang R, Zheng Z, Su Z, Fan P, Zhang X, Zhang Y, Chen S 2021 J. Materiomics 7 1324Google Scholar

    [10]

    Perdew J P, Burke K, Ernzerhof M 1996 Phys. Rev. Lett. 77 3865Google Scholar

    [11]

    Vadapoo R, Krishnan S, Yilmaz H, Marin C 2011 Phys. Status Solidi B 248 700Google Scholar

    [12]

    Bardekn J, Shockley W 1950 Phys. Rev. 80 72Google Scholar

    [13]

    Xi J, Long M, Tang L, Wang D, Shuai Z 2012 Nanoscale 4 4348Google Scholar

    [14]

    El-Sayad E A, Moustafa A M, Marzouk S Y 2009 Physica B 404 1119Google Scholar

    [15]

    Zheng X, Xie Y, Zhu L, Jiang X, Jia Y, Song W, Sun Y 2002 Inorg. Chem. 41 455Google Scholar

    [16]

    Effective Mass Calculator for Semiconductors, Fonari A, Sutton C https://github.com/afonari/emc [2013-3-18

    [17]

    Zhang B, Qian X 2022 ACS Appl. Energy Mater. 5 492Google Scholar

    [18]

    Wang X, Li Z, Kavanagh S R, Ganose A M, Walsh A 2022 Phys. Chem. Chem. Phys. 24 7195Google Scholar

    [19]

    Wang V, Xu N, Liu J C, Tang G, Geng W T 2021 Comput. Phys. Commun. 267 108033Google Scholar

    [20]

    Silva E Lora da, Skelton J M, Rodríguez-Hernández P, Muõz A, Santo M C, Martínez-García D, Vilaplana R, Manjón F J 2022 J. Mater. Chem. C 10 15061Google Scholar

    [21]

    Zhou Y, Leng M, Xia Z, Zhong J, Song H, Liu X, Yang B, Zhang J, Chen J, Zhou K 2014 Adv. Energy Mater. 4 1301846Google Scholar

    [22]

    Madelung O 1996 Semiconductor: Data Handbook (2rd Ed.) (New York: Springer-Verlag Berlin Heidelbergy) p204

    [23]

    Black J, Conwell E, Seigle L, Spencer C 1957 J. Phys. Chem. Solids 2 240Google Scholar

    [24]

    Cheng L, Liu Y 2018 J. Am. Chem. Soc. 140 17895Google Scholar

    期刊类型引用(0)

    其他类型引用(2)

  • 图 1  Sb2Se3的晶胞

    Fig. 1.  Crystal structure of Sb2Se3 computational model.

    图 2  纯Sb2Se3 (a), As-Sb2Se3 (b), Bi-Sb2Se3 (c)的能带结构及分态密度图. 虚线处是费米能级, 蓝圈处是价带顶(VBM)

    Fig. 2.  Band structures and partial density of states of pure Sb2Se3 (a), As-Sb2Se3 (b), and Bi-Sb2Se3 (c). Dashed line represents the Fermi level and the blue circle represents the valence band maximum (VBM).

    图 3  (a)本征Sb2Se3总能量与形变的抛物线拟合, 其决定了弹性常数; (b)不同结构的弹性常数

    Fig. 3.  (a) Parabolic fitting of the intrinsic Sb2Se3 total energy and deformations, determining the elastic constants; (b) calculated elastic constants of different structures.

    图 4  (a)本征Sb2Se3在3个不同方向应变下的价带边缘位置, 实线是线性拟合, 其决定了形变势; (b)不同结构的形变势能

    Fig. 4.  (a) Valence band edge positions of intrinsic Sb2Se3 under strain along three different directions, solid lines represent linear fitting, determining the deformation potential; (b) calculated deformation potential energies of different structures.

    表 1  优化后的晶格参数

    Table 1.  Optimized lattice parameters.

    文献[14] 文献[15] Sb2Se3 AsSb2Se3 BiSb2Se3
    a 3.98 3.99 3.99 3.97 4.00
    b 11.64 11.63 11.35 11.35 11.36
    c 11.79 11.78 11.65 11.63 11.65
    v3 547.1 546.6 527.3 524.7 528.9
    下载: 导出CSV

    表 2  不同Sb2Se3结构的空穴有效质量

    Table 2.  Effective mass of holes for different structures of Sb2Se3.

    m*/m0Sb2Se3As-Sb2Se3Bi-Sb2Se3
    $ {m}_{xx}^{*} $0.430.450.42
    $ {m}_{yy}^{*} $0.880.890.67
    $ {m}_{zz}^{*} $1.081.221.61
    $ {\stackrel{-}{m}}^{*} $0.680.720.67
    下载: 导出CSV

    表 3  三种结构的空穴迁移率

    Table 3.  Hole mobility of Sb2Se3, As-Sb2Se3 and Bi-Sb2Se3 along three principle directions.

    迁移率/(cm2·V–1·s–1)Sb2Se3As-Sb2Se3Bi-Sb2Se3
    μx232.62221.59240.66
    μy32.7831.2066.71
    μz20.0215.118.04
    μavg95.1489.30105.13
    下载: 导出CSV
  • [1]

    Chen C, Li K H, Tang J 2022 Sol. RRL 6 2200094Google Scholar

    [2]

    Zhang X, Li C, Sun K, Zhou J, Zhang Z 2021 Adv. Energy Mater. 11 2002614Google Scholar

    [3]

    薛丁江, 石杭杰, 唐江 2015 物理学报 64 038406Google Scholar

    Xue D J, Shi H J, Tang J 2015 Acta Phys. Sin. 64 038406Google Scholar

    [4]

    Zhao Y, Wang S, Li C, Che B, Chen X, Chen H, Tang R, Wang X, Chen G, Wang T, Gong J, Chen T, Xiao X 2022 Energy Environ. Sci. 15 5118Google Scholar

    [5]

    Li Z, Liang X, Li G, Liu H, Zhang H, Guo J, Chen J, Shen K, San X, Yu W, Schropp R, Mai Y 2019 Nat. Commun. 10 125Google Scholar

    [6]

    Wang X, Ganose A M, Kavanagh S R, Walsh A 2022 ACS Energy Lett. 7 2954Google Scholar

    [7]

    Spaggiari G, Bersani D, Calestani D, Gilioli E, Gombia E, Mezzadri F, Casappa M, Pattini F, Trevisi G, Rampino S 2022 Int. J. Mol. Sci. 23 15529Google Scholar

    [8]

    Huang M, Lu S, Li K, Lu Y, Chen C, Tang J, Chen S 2022 Sol. RRL 6 2100730Google Scholar

    [9]

    Liang G, Chen X, Ren D, Jiang X, Tang R, Zheng Z, Su Z, Fan P, Zhang X, Zhang Y, Chen S 2021 J. Materiomics 7 1324Google Scholar

    [10]

    Perdew J P, Burke K, Ernzerhof M 1996 Phys. Rev. Lett. 77 3865Google Scholar

    [11]

    Vadapoo R, Krishnan S, Yilmaz H, Marin C 2011 Phys. Status Solidi B 248 700Google Scholar

    [12]

    Bardekn J, Shockley W 1950 Phys. Rev. 80 72Google Scholar

    [13]

    Xi J, Long M, Tang L, Wang D, Shuai Z 2012 Nanoscale 4 4348Google Scholar

    [14]

    El-Sayad E A, Moustafa A M, Marzouk S Y 2009 Physica B 404 1119Google Scholar

    [15]

    Zheng X, Xie Y, Zhu L, Jiang X, Jia Y, Song W, Sun Y 2002 Inorg. Chem. 41 455Google Scholar

    [16]

    Effective Mass Calculator for Semiconductors, Fonari A, Sutton C https://github.com/afonari/emc [2013-3-18

    [17]

    Zhang B, Qian X 2022 ACS Appl. Energy Mater. 5 492Google Scholar

    [18]

    Wang X, Li Z, Kavanagh S R, Ganose A M, Walsh A 2022 Phys. Chem. Chem. Phys. 24 7195Google Scholar

    [19]

    Wang V, Xu N, Liu J C, Tang G, Geng W T 2021 Comput. Phys. Commun. 267 108033Google Scholar

    [20]

    Silva E Lora da, Skelton J M, Rodríguez-Hernández P, Muõz A, Santo M C, Martínez-García D, Vilaplana R, Manjón F J 2022 J. Mater. Chem. C 10 15061Google Scholar

    [21]

    Zhou Y, Leng M, Xia Z, Zhong J, Song H, Liu X, Yang B, Zhang J, Chen J, Zhou K 2014 Adv. Energy Mater. 4 1301846Google Scholar

    [22]

    Madelung O 1996 Semiconductor: Data Handbook (2rd Ed.) (New York: Springer-Verlag Berlin Heidelbergy) p204

    [23]

    Black J, Conwell E, Seigle L, Spencer C 1957 J. Phys. Chem. Solids 2 240Google Scholar

    [24]

    Cheng L, Liu Y 2018 J. Am. Chem. Soc. 140 17895Google Scholar

  • [1] 张冷, 沈宇皓, 汤朝阳, 吴孔平, 张鹏展, 刘飞, 侯纪伟. 单轴应变对Sb2Se3空穴迁移率的影响. 物理学报, 2024, 73(11): 117101. doi: 10.7498/aps.73.20240175
    [2] 黄昊, 牛奔, 陶婷婷, 罗世平, 王颖, 赵晓辉, 王凯, 李志强, 党伟. Sb2Se3薄膜表面和界面超快载流子动力学的瞬态反射光谱分析. 物理学报, 2022, 71(6): 066402. doi: 10.7498/aps.71.20211714
    [3] 张小娅, 宋佳讯, 王鑫豪, 王金斌, 钟向丽. In掺杂h-LuFeO3光吸收及极化性能的第一性原理计算. 物理学报, 2021, 70(3): 037101. doi: 10.7498/aps.70.20201287
    [4] 颜俊, 王子毅, 曾若生, 邹炳锁. 零维Sb3+掺杂Rb7Bi3Cl16金属卤化物的三重态自陷激子发射. 物理学报, 2021, 70(24): 247801. doi: 10.7498/aps.70.20211024
    [5] 袁国才, 陈曦, 黄雨阳, 毛俊西, 禹劲秋, 雷晓波, 张勤勇. Mg2Si0.3Sn0.7掺杂Ag和Li的热电性能对比. 物理学报, 2019, 68(11): 117201. doi: 10.7498/aps.68.20190247
    [6] 阮璐风, 王磊, 孙得彦. Sr掺杂对La1-xSrxMnO3/LaAlO3/SrTiO3界面电子结构的影响. 物理学报, 2017, 66(18): 187301. doi: 10.7498/aps.66.187301
    [7] 徐晶, 梁家青, 李红萍, 李长生, 刘孝娟, 孟健. Ti掺杂NbSe2电子结构的第一性原理研究. 物理学报, 2015, 64(20): 207101. doi: 10.7498/aps.64.207101
    [8] 程超群, 李刚, 张文栋, 李朋伟, 胡杰, 桑胜波, 邓霄. B, P掺杂β-Si3N4的电子结构和光学性质研究. 物理学报, 2015, 64(6): 067102. doi: 10.7498/aps.64.067102
    [9] 廖建, 谢召起, 袁健美, 黄艳平, 毛宇亮. 3d过渡金属Co掺杂核壳结构硅纳米线的第一性原理研究. 物理学报, 2014, 63(16): 163101. doi: 10.7498/aps.63.163101
    [10] 曹娟, 崔磊, 潘靖. V,Cr,Mn掺杂MoS2磁性的第一性原理研究. 物理学报, 2013, 62(18): 187102. doi: 10.7498/aps.62.187102
    [11] 吴木生, 徐波, 刘刚, 欧阳楚英. Cr和W掺杂的单层MoS2电子结构的第一性原理研究. 物理学报, 2013, 62(3): 037103. doi: 10.7498/aps.62.037103
    [12] 徐金荣, 王影, 朱兴凤, 李平, 张莉. N掺杂和N-V共掺杂锐钛矿相TiO2的第一性原理研究. 物理学报, 2012, 61(20): 207103. doi: 10.7498/aps.61.207103
    [13] 周传仓, 刘发民, 丁芃, 钟文武, 蔡鲁刚, 曾乐贵. 钶铁矿型MnNb2O6的熔盐法合成、钒掺杂与磁性研究. 物理学报, 2011, 60(4): 048101. doi: 10.7498/aps.60.048101
    [14] 张云, 邵晓红, 王治强. 3C-SiC材料p型掺杂的第一性原理研究. 物理学报, 2010, 59(8): 5652-5660. doi: 10.7498/aps.59.5652
    [15] 徐新发, 邵晓红. Y掺杂SrTiO3晶体材料的电子结构计算. 物理学报, 2009, 58(3): 1908-1916. doi: 10.7498/aps.58.1908
    [16] 杜丽萍, 陈抱雪, 孙 蓓, 陈 直, 邹林儿, 浜中广见, 矶 守. 掺杂As2S8非晶态薄膜波导的光阻断效应. 物理学报, 2008, 57(6): 3593-3599. doi: 10.7498/aps.57.3593
    [17] 金胜哲, 黄祖飞, 明 星, 王春忠, 孟 醒, 陈 岗. 二价金属元素掺杂对LiCoO2体系电子输运性质的影响. 物理学报, 2007, 56(10): 6008-6012. doi: 10.7498/aps.56.6008
    [18] 王先杰, 隋 郁, 千正男, 刘志国, 苗继鹏, 黄喜强, 吕 喆, 朱瑞滨, 程金光, 苏文辉. Fe位Al掺杂对Sr2FeMoO6磁结构和磁输运性质的影响. 物理学报, 2006, 55(2): 849-853. doi: 10.7498/aps.55.849
    [19] 张加宏, 马 荣, 刘 甦, 刘 楣. 掺杂MgCNi3超导电性和磁性的第一性原理研究. 物理学报, 2006, 55(9): 4816-4821. doi: 10.7498/aps.55.4816
    [20] 代月花, 陈军宁, 柯导明, 孙家讹, 胡 媛. 纳米MOSFET迁移率解析模型. 物理学报, 2006, 55(11): 6090-6094. doi: 10.7498/aps.55.6090
  • 期刊类型引用(0)

    其他类型引用(2)

计量
  • 文章访问数:  2532
  • PDF下载量:  89
  • 被引次数: 2
出版历程
  • 收稿日期:  2023-08-31
  • 修回日期:  2023-11-20
  • 上网日期:  2023-11-24
  • 刊出日期:  2024-02-20

/

返回文章
返回