全文HTML
--> --> -->CdTe | CdS | *CdTe[25] | *CdS[25] | #CdTe[26] | #CdS[27] | |
a/? | 6.642 | 4.214 | 6.646 | 4.212 | 6.481 | 4.140 |
b /? | 6.642 | 4.214 | 6.646 | 4.212 | 6.481 | 4.140 |
c /? | 6.642 | 6.850 | 6.646 | 6.858 | 6.481 | 6.720 |
注: *为理论值, #为实验值. |
表1CdTe和CdS晶格参数优化结果
Table1.Optimal lattice parameters of CdTe and CdS.
本研究采用基于密度泛函理论(DFT)的Material Studio软件的CASTEP模块进行计算. 计算过程中采用投影缀加平面波(PAW)方法描述赝势[28-30], 广义梯度近似(GGA-PBE)泛函确定交换关联能[31]. 平面波基的截止能量设置为400 eV, 布里渊区的Monkhorst-Pack网格设置为4 × 4 × 1. 由于Mn原子的存在, 计算过程中考虑了自旋轨道耦合(SOC)因素[32]. 所有原子均采用Broyden-Fletcher-Goldfarb-Shanno (BFGS)方法进行几何优化, 直到每个原子上的离子力小于0.02 eV/A, 总能量变化小于1.0 × 10–6 eV.
3.1.CdS (002)和CdMnTe (111)表面模型和电子性质
33.1.1.CdS(002)和CdMnTe(111)表面模型
薄膜生长实验发现[33,34], CdMnTe薄膜太阳能电池的CdS和CdMnTe的择优生长面分别为(002)和(111)晶面, 因此计算中选取(2 × 2)-CdS (002)表面(如图1(a)所示)和(图 1 表面模型的建立 (a) CdS (002)表面模型; (b) CdMnTe (111)-A1表面模型; (c) CdMnTe (111)-B1表面模型
Figure1. Building of the surface models: (a) Model of the CdS (002) surface; (b) A1 surface model of CdMnTe (111); (c) B1 surface model of CdMnTe (111).
结构弛豫后的CdS (002)和CdMnTe (111)表面原子位置和键长的变化如图2所示. 图2(a)中CdS (002)表面的Cd原子沿z轴向上移动了0.035 ?, S原子沿z轴向下移动了0.052 ?. 图2(b)中CdMnTe (111)-A1表面的Cd原子沿z轴向下移动了0.021 ?, Te原子沿z轴向上移动了0.013 ?. 图2(c)中B1表面模型的Cd原子沿z轴向下移动了0.073 ?, Te原子沿z轴向上移动了0.040 ?. 层间距离大小和层间距离变化的计算结果见表2. 从表2中的数据可以看出, 弛豫后第1层与第2层, 第2层与第3层之间的原子键长与层间距离的变化比第3层与第4层, 第4层与第5层之间的变化大, 说明层间收缩或扩张主要发生在表面, 内部原子键长与层间距离基本保持不变.
图 2 弛豫后的表面原子位置和层间距离的变化 (a) CdS (002)表面模型; (b) CdMnTe (111)-A1表面模型; (c) CdMnTe (111)-B1表面模型
Figure2. Variations of atomic positions and interlayer spacing after relaxation: (a) Model of the CdS (002) surface; (b) A1 surface model of CdMnTe (111); (c) B1 surface model of CdMnTe (111).
Surface | items | d1-2/? | d2-3/? | d3-4/? | d4-5/? |
CdMnTe-A1 | dCd-Te | 2.829 | 2.837 | 2.831 | 2.831 |
dMn-Te | 2.605 | ||||
?d | 0.138 | 0.089 | 0.000 | 0.000 | |
CdMnTe-B1 | dCd-Te | 2.791 | 2.841 | 2.831 | 2.831 |
dMn-Te | 2.697 | ||||
?d | 0.182 | 0.093 | 0.000 | 0.000 | |
CdS | dCd-S | 2.534 | 2.538 | 2.536 | 2.536 |
?d | 0.120 | 0.109 | 0.000 | 0.000 |
表2弛豫后CdMnTe (111)和CdS (002)表面层间距离的变化
Table2.Variations of interlayer spacing of CdMnTe (111) and CdS (002) surfaces after relaxation.
3
3.1.2.CdS (002)和CdMnTe (111)表面的电子性质
CdS (002)和CdMnTe (111)表面模型的总态密度及其第1层和第5层的局域态密度如图3所示. 图3(a)所示的CdS (002)表面模型总态密度中费米能级附近出现了额外的电子态, 即表面态, 从第1层和第5层的局域态密度可知, 该表面态主要是由CdS表面第1层中S-3p, Cd-4p, Cd-5s轨道引起的. 第5层由于距离表面较远, 受到表面的影响较小, 表现出较好的体相特性. 表面态出现的主要原因是表面原子结构的周期性遭到破坏, 未能形成成对电子而导致电子密度的增加. 这些新的表面电子态对CdS表面的光学性质具有一定的影响.图 3 表面模型总态密度和第1层及第5层局域态密度 (a) CdS (002); (b) CdMnTe (111)-A1; (c) CdMnTe (111)-B1
Figure3. Total density of states and the local density of states of the first layer and the fifth layer: (a) CdS (002); (b) A1 of CdMnTe (111); (c) B1 of CdMnTe (111).
图3(b)代表了CdMnTe (111)-A1表面模型的总态密度和局域态密度, 其总态密度中费米能级附近也出现了新的表面态. 该表面态是由CdMnTe表面第1层中Mn-3d和Te-5p引起的, 远离表面的内部在费米能级附近没有出现新的电子态. 图3(c)所示的B1表面模型中总态密度的费米能级附近也出现了由第1层中Mn-3d和Te-5p轨道引起的表面态, 但其峰值略低.
2
3.2.CdS/CdMnTe异质结模型和电子性质
33.2.1.CdS/CdMnTe异质结模型
由CdS (002)和CdMnTe (111)表面模型构建的CdS/CdMnTe异质结模型如图4所示, 模型A2为Mn原子位于界面层(Layer 1)的异质结模型(Mn原子位置同模型A1), B2为Mn原子位于中间层(Layer 2)的异质结模型(Mn原子位置同模型B1). 异质结包含77个原子, 化学式为H7S18Cd34MnTe15, 晶格参数为a = b = 8.15 ?, 晶格失配度约为3.5%. 同样在第5层外添加20 ?的真空层以减小周期性边界条件对异质结模型的影响, 采用氢原子钝化最外层原子的悬挂键以减小表面态对异质结界面的影响. 另外, 令界面两端第1, 2原子层自由弛豫, 第3—5原子层固定. 采用单点能量计算方法来获得合适的界面距离, 如图5所示. 可以看出, 当界面距离为2.26 ?时, 模型A2系统总能量最低; 当界面距离为2.55 ?时, 模型B2系统总能量最低.图 4 CdS/CdMnTe异质结模型 (a) Mn原子位于界面层的异质结模型(模型A2); (b) Mn原子位于中间层的异质结模型(模型B2)
Figure4. CdS/CdMnTe heterojunction models: (a) Mn atom being in the interface layer (model A2); (b) Mn atom in the internal layer (model B2).
图 5 模型A2和B2的总能量随界面间距的变化
Figure5. Variation of total energy dependent on the interface distance of the models A2 and B2, respectively.
采用BFGS方法分别对两种模型进行结构优化, 优化后的模型如图6所示. 图6(a)和图6(b)分别表示模型A2弛豫前后原子键长的变化, 模型A2弛豫后界面上原子键长为dCd—Te = 2.902—2.994 ?(弛豫前2.871 ?), dMn—Te = 2.740—2.720 ? (弛豫前2.651 ?), dCd—S = 2.531—2.666 ?(弛豫前2.497 ?); 界面附近原子键长为dCd—Te = 2.841—2.852 ?(弛豫前2.969 ?), dCd—S = 2.502—2.547 ? (弛豫前2.542 ?), 远离界面处的原子键长和位置基本保持不变. 图6(c)和图6(d)分别表示模型B2弛豫前后原子键长的变化, 弛豫后界面上的原子键长为dCd—Te = 2.871—2.942 ?(弛豫前2.832 ?), dCd—S = 2.541—2.579 ?(弛豫前2.503 ?); 界面附近原子的键长为dCd—Te = 2.813—2.814 ?(弛豫前2.906 ?), dMn—Te = 2.696 ?(弛豫前2.661 ?), dCd—S = 2.543—2.584 ?(弛豫前2.759 ?), 同样远离界面处的原子键长和位置也基本保持不变.
图 6 弛豫前后界面及其附近原子键长的变化 (a)弛豫前的模型A2; (b)弛豫后的模型A2; (c)弛豫前的模型B2; (d)弛豫后的模型B2
Figure6. Bond lengths in/near of the interface before and after relaxation: (a) Before relaxation of A2; (b) after relaxation of A2; (c) before relaxation of B2; (d) after relaxation of B2.
弛豫前后界面上原子键长的变化是由界面处原子的种类和结构的差异造成的. 在两种异质结模型中, 模型B2弛豫前后界面附近原子键长变化的幅度偏大, 说明模型B2界面处原子间作用力较大, 界面原子的结合更加稳定, 更益于形成稳定的界面.
根据CdS/CdMnTe异质结界面结合能的大小可以判断其结构的稳定性, 若结合能为负值, 则异质结模型是稳定的, 且负值越大模型越稳定[35,36]. 界面结合能的计算公式为
3
3.2.2.CdS/CdMnTe异质结的电子性质
下面通过异质结的总态密度、局域态密度、原子的分态密度和差分电荷密度研究异质结的电子性质. 图7(a)和图7(b)分别表示模型A2和B2的总态密度和局域态密度, 可以看出, 异质结A2和B2在费米能级附近约–0.15—0.20 eV范围内均出现了新的电子态, 即界面态, 并且模型A2中费米能级附近出现了更高的界面态. 界面态在载流子传输过程中可以充当俘获中心, 阻止光生载流子通过费米能级, 导致带隙减小, 开路电压和短路电流降低. 局域态密度分析发现, 异质结的界面态主要是由CdMnTe (111)第5层引起的, 而CdS (002)和CdMnTe (111)第1层的费米能级附近均没有界面态. 也就是说, 这两个异质结界面不会捕获载流子, 避免了费米能级在界面处被钉扎. 这可能是因为界面处原子的紧密结合促进了共价键的形成, 消除了界面悬挂键而使电子达到了饱和状态, 故CdS/CdMnTe异质结界面处不显示界面态.图 7 CdS/CdMnTe异质结总态密度与局域态密度 (a)模型A2; (b)模型B2
Figure7. Total density of states and local density of states of the CdS/CdMnTe heterojunction: (a) Model A2; (b) model B2.
图8表示CdS/CdMnTe异质结界面处原子的分态密度, 图8(a)说明模型A2界面处所有的原子在费米能级附近均没有其他电子态. 在–5 eV附近, Mn-3 d和Te-5p轨道发生微弱杂化, 在0.5—1.5 eV处, Mn-3p和Te-5p轨道也发生了较弱的杂化, 在–3—–1 eV处, Mn-3p与Cd-3p轨道发生部分杂化, 在1.5 —2.5 eV处, Mn-4s, Cd-3p与S-3s电子发生较小的轨道杂化. 这些原子轨道的杂化作用促进了界面处Mn与Te, Mn与S, S与Cd原子之间的成键. 对于模型B2, 如图8(b)所示, 在–5—–4.5 eV与–3.5—–1.5 eV处, Te-5p和S-3p发生强烈的轨道杂化, 在1—2 eV处, Cd-4p和S-3s轨道发生杂化, 在–5 eV附近, Cd-5s和S-3p轨道发生微弱的杂化. 此外, H原子1s轨道在费米能级处出现了界面态, 说明H原子的加入对CdS/CdMnTe异质结有不利的影响.
图 8 CdS/CdMnTe异质结界面处原子的分态密度 (a)模型A2; (b)模型B2
Figure8. Part density of states of different atoms in the CdS/CdMnTe interface: (a) Model A2; (b) model B2.
界面上原子轨道杂化能力的强弱是影响界面结合程度的关键因素, 异质结B2具有比异质结A2更强的界面原子轨道杂化能力, 说明异质结B2具有更加稳定的结构.
异质结沿z轴方向的电荷密度差可以由下式表示:
垂直于z轴方向的电荷重新分布如图9所示, 图中蓝色区域表示电荷累积, 黄色区域表示电荷耗尽. 从图9(a)和图9(c)所示的差分电荷示意图中可以看出, 电荷的重新分配主要发生在CdS/CdMnTe界面处, 这是由于界面处原子轨道的杂化作用促进了界面的结合, 远离界面处的原子由于CdMnTe与CdS两相之间的相互作用较弱几乎没有电荷转移. 沿z轴方向的平面平均电荷密度差代表电荷的重新分布, 正值表示电子累积, 负值表示电子耗尽. 通过界面处偶极矩变化的定量计算, 表明电子主要通过界面从CdMnTe侧转移到CdS侧. 此外, 通过Mulliken电荷计算出异质结A2电荷转移量为0.89, 异质结B2电荷转移量为0.91, 此结果表明异质结B2具有较强的电荷转移能力.
图 9 电荷重新分布图 (a), (c)异质结A2与B2的差分电荷示意图和相应平面差分电荷密度曲线, 黄色区域代表电荷消耗, 蓝色区域代表电荷积累, 等值面值为0.00103 e/?3; (b), (d) A2与B2沿z轴方向的平均静电势
Figure9. Distribution diagram of charges: (a), (c) Charge density difference and the corresponding planar differential charge density curve of the A2 model and B2 model, respectively. The yellow region represents charge depletion, the bule region indicates charge accumulation, the isosurface value is 0.00103 e/?3. (b), (d) the average electrotactic potential difference of the A2 model and B2 model, respectively, along the direction of z axis.
图9(b)和图9(d)分别表示CdS/CdMnTe异质结A2和B2模型界面间的静电势. 计算结果显示, 异质结B2在界面处的电位更低, 表明电子更容易在此界面区域转移和积累. Guo Y等[37]研究表明, 在CdS/CdMnTe异质结界面处将发生载流子的有效分离, 较大的导带偏移有助于增强电子注入CdS层的能力. 因为异质结B2具有更高的界面结合能和更强的电荷转移能力, 这将有助于减小CdS/CdMnTe价带偏移而增加导带偏移, 对促进载流子的分离具有积极的意义.
2
3.3.CdS (002), CdMnTe (111)表面及其异质结的光学性质
研究CdS (002), CdMnTe (111)表面及其异质结的光学性质对于扩展多结串联电池的光谱响应范围具有重要意义[38]. 介电函数图 10 表面模型和异质结模型的吸收光谱 (a) CdS (002), CdMnTe (111)-B1表面和异质结B2模型; (b)异质结A2和B2模型
Figure10. Absorption spectra of surface models and heterojunction models: (a) Surface CdS (002), surfaces of CdMnTe (111)-B1, and B2 model of CdMnTe/CdS heterojunction model; (b) heterojunction A2 model and B2 model.