针对CsPbI3基钙钛矿太阳能电池的全面数值与模拟研究
本次向您介绍的研究工作由M. Khalid Hossain（通讯作者，隶属于孟加拉国原子能委员会原子能研究所以及日本九州大学跨学科工程科学研究生院先进能源工程科学系）领衔，与Mirza Humaun Kabir Rubel（拉杰沙希大学材料科学与工程系）、G. F. Ishraque Toki（东华大学材料科学与工程学院）、Intekhab Alam（卡尔加里大学机械与制造工程系）、Md. Ferdous Rahman（贝古姆·罗克亚大学电气与电子工程系）和H. Bencherif（阿尔及利亚国家高等可再生能源、环境与可持续发展学院）合作完成。该研究以论文《Effect of Various Electron and Hole Transport Layers on the Performance of CsPbI3‑Based Perovskite Solar Cells: A Numerical Investigation in DFT, SCAPS‑1D, and wxAMPS Frameworks》的形式，于2022年11月14日在线发表在ACS Omega期刊上（2022年，第7卷，第43210-43230页）。
本研究属于材料科学与光伏技术领域。其学术背景源于21世纪日益增长的能源需求与化石燃料带来的环境污染问题。钙钛矿太阳能电池（Perovskite Solar Cells, PSCs）作为第三代光伏技术，因其制造工艺简单、成本低廉而备受关注。其中，无机CsPbI3钙钛矿因其与杂化钙钛矿相似的光电特性以及显著改善的热稳定性，被视为极具潜力的光吸收材料。然而，CsPbI3基PSCs的实际性能仍未达到杂化PSCs的高水平，且其光电性能易受环境湿度和温度影响，从黑色的钙钛矿相转变为黄色的非钙钛矿相。为了获得高性能、低成本且稳定的CsPbI3基PSCs，选择合适的电子传输层（Electron Transport Layer, ETL）和空穴传输层（Hole Transport Layer, HTL）至关重要。这两个层在器件的可靠性和性能中扮演关键角色，它们影响电荷的提取、传输以及复合过程。传统上，研究者多依赖实验试错来寻找最佳材料组合，但这种方法耗时且昂贵。因此，本研究旨在通过先进行第一性原理计算，再结合专业太阳能电池模拟软件进行系统性数值筛选，从理论层面探究并优化CsPbI3与不同ETL/HTL的组合，以期快速识别出能实现高效率的器件结构，为后续实验制备提供明确的指导。研究目标具体包括：1) 利用密度泛函理论（Density Functional Theory, DFT）详细分析立方相CsPbI3的结构、电子和光学性质；2) 利用SCAPS-1D软件，系统地筛选多种ETL和HTL与CsPbI3的组合，找出最佳性能的配对；3) 对最优的几个器件结构，深入分析其厚度、串联/并联电阻、工作温度等参数对性能的影响，并研究其内部物理过程；4) 使用另一款模拟软件wxAMPS对SCAPS-1D的结果进行交叉验证。
本研究的工作流程可分为三个核心步骤，环环相扣，从材料基础物性计算到器件性能的全面优化与评估。
第一步是材料层面的第一性原理计算。 研究对象是立方相（空间群Pm-3m）的CsPbI3单晶结构。研究团队采用剑桥序列总能量包（Cambridge Serial Total Energy Package, CASTEP）代码，在DFT框架下进行了一系列计算。这里使用了广义梯度近似（Generalized Gradient Approximation, GGA）下的Perdew-Burke-Ernzerhof（PBE）泛函来处理交换关联能。研究中首先对晶体结构进行几何优化，以获得稳定的构型和准确的晶格常数。计算中设置了520 eV的截断能，并采用了12×12×12的k点网格进行布里渊区采样。对于电子电荷密度图和费米面拓扑的计算，则使用了更密的17×17×17 k点网格以确保精度。这一步骤是后续器件模拟的基石，其核心任务是获取CsPbI3的准确本征属性。具体而言，他们计算了能带结构以确定其带隙类型和大小，计算了态密度（Density of States, DOS）以分析不同原子轨道对电子结构的贡献，绘制了电子电荷密度图以研究化学键性质，描绘了费米面拓扑以判断材料的电子输运特性，并计算了多种光学性质（如介电函数、折射率、吸收系数等）以评估材料的光学响应。所有DFT计算均基于第一性原理，无需依赖实验参数，旨在从量子力学层面深入理解CsPbI3的内在特性。
第二步是器件结构的系统性筛选与初步优化。 这是本研究的核心环节。研究对象不再是单一材料，而是由多层薄膜构成的完整平面异质结PSC器件结构，基本构型为：前电极（ITO）/ ETL / CsPbI3吸收层 / HTL / 背电极（Au）。研究团队设定了8种ETL材料（包括PCBM、TiO₂、ZnO、C60、IGZO、SnO₂、WS₂、CeO₂）和12种HTL材料（包括Cu₂O、CuSCN、CuSbS₂、Spiro-MeOTAD、V₂O₅、CBTS、CFTS、P3HT、PEDOT:PSS、NiO、CuO、CuI）。通过排列组合，理论上研究了8 × 12 = 96种不同的器件构型。这一庞大的筛选工作是通过专业的太阳能电池模拟软件SCAPS-1D（一维太阳能电池电容模拟器）完成的。该软件通过数值求解半导体器件中的泊松方程和载流子连续性方程，能够精确模拟器件在光照和偏压下的电流-电压（J-V）特性等性能参数。在模拟中，研究者为每一层材料（ITO、ETLs、CsPbI3、HTLs）输入了基于文献报道和DFT计算结果的详细物理参数，例如厚度、带隙、电子亲和能、介电常数、态密度、载流子迁移率、掺杂浓度等（具体数值见原文表1和表2）。为了更接近实际器件，他们还在ETL/CsPbI3和CsPbI3/HTL两个界面处引入了特定密度的缺陷态（见原文表3）。所有模拟均在标准测试条件（AM1.5G光谱，1000 W/m²光照强度，300 K温度，1 MHz频率）下进行。此步骤的目标是在保持CsPbI3吸收层厚度、掺杂密度等参数不变的前提下，单纯通过改变ETL和HTL的材料种类，找出能产生最高光电转换效率（Power Conversion Efficiency, PCE）的配对组合。SCAPS-1D软件的高效性使得在短时间内评估所有96种组合成为可能，这是实验方法难以企及的。
第三步是对最优器件结构的深入分析与验证。 在完成初步筛选后，研究聚焦于性能最好的六种器件结构（均使用CBTS作为HTL，分别搭配PCBM、TiO₂、ZnO、C60、IGZO、WS₂作为ETL）。针对这六个“优胜”结构，研究进行了更深入的参数化研究，以评估其性能的鲁棒性和优化潜力。具体包括：1) 厚度优化：系统地改变CsPbI3吸收层（400-2200 nm）和各ETL（50-500 nm）的厚度，通过绘制等高线图来观察其对开路电压（Voc）、短路电流密度（Jsc）、填充因子（FF）和PCE的影响。2) 电阻分析：研究串联电阻（Rs，从0到6 Ω·cm²）和并联电阻（Rsh，从10到10⁷ Ω·cm²）对器件性能参数的影响。3) 温度稳定性：考察工作温度从275 K升至450 K时，各性能参数的变化趋势。4) 内部物理过程分析：模拟了这六个器件的能带排列、电容-电压（C-V）特性、载流子产生与复合速率分布、量子效率（QE）谱等，以深入理解其工作原理和性能差异的根本原因。5) 交叉验证：为了确保SCAPS-1D模拟结果的可靠性，研究者使用另一款广泛认可的器件模拟软件wxAMPS（Widget-Provided Analysis of Microelectronic and Photonic Structures）对六个最优器件的性能进行了重复模拟，并将结果进行了对比。
本研究的主要结果如下：
在DFT计算部分，研究获得了立方相CsPbI3的详细材料特性。优化后的晶格常数为6.407 Å，与实验值高度吻合，且形成能为负值（-6.905 eV/原子），证实了其结构稳定性。电子结构计算显示，CsPbI3具有直接带隙，采用GGA-PBE方法计算得到的带隙值为1.483 eV（注：此值通常因PBE泛函的局限而低估，文中用于器件模拟的带隙采用了更接近实验值的1.694 eV）。态密度分析表明，费米能级附近的电子态主要来源于Pb-5d轨道，其次是Cs-6s轨道，表明Pb原子在光电过程中起主导作用。电子电荷密度图证实了Pb和I之间较强的共价键特性，而Cs和I之间则主要为离子键。费米面拓扑结构同时显示出电子型和空穴型费米面，表明CsPbI3具有多能带特征。光学性质计算给出了介电函数、折射率、光吸收系数等一系列参数，其中吸收系数曲线表明该材料在可见光到紫外区域有良好的光吸收能力。
在SCAPS-1D筛选部分，研究得出了关键结论。在96种组合中，性能最优的器件结构为 ITO/TiO₂/CsPbI₃/CBTS/Au。该结构模拟得到的PCE为17.90%，对应的Voc为0.997 V，Jsc为21.075 mA/cm²，FF高达85.21%。研究发现，当与性能优异的CBTS HTL搭配时，TiO₂、ZnO、WS₂和IGZO等ETL都能使器件获得接近18%的PCE。而CeO₂和SnO₂作为ETL时性能普遍较差。从能带排列角度看，高性能的ETL（如TiO₂、ZnO）其导带底与CsPbI3的导带底匹配良好，有利于电子提取；同时，高性能的HTL（如CBTS、CuSCN）其价带顶与CsPbI3的价带顶匹配良好，有利于空穴提取，从而减少了界面处的载流子复合损失。
在深入分析部分，研究揭示了各因素对最优器件性能的影响规律。厚度优化结果显示，对于TiO₂、ZnO、WS₂等ETL，存在一个最佳的吸收层厚度范围（通常在1000-2000 nm）和ETL厚度范围（通常较薄，小于250 nm），在此范围内器件效率最高（最高模拟效率可达18.40%左右）。串联电阻Rs的增大会显著降低FF和PCE，但对Voc和Jsc影响较小。并联电阻Rsh的增大则能显著提升所有性能参数，特别是当Rsh从较低值增加到约10³ Ω·cm²时，提升效果明显，之后趋于饱和。温度升高会导致Voc、FF和PCE的显著下降，主要归因于缺陷密度增加和载流子复合加剧。C-V分析表明所有器件在零偏压下处于耗尽状态。载流子产生和复合速率分布图显示，在吸收层靠近ETL一侧的区域（约0.9-1.0 μm深度）是光生载流子产生和复合最强烈的区域。J-V和QE曲线进一步证实了使用TiO₂、ZnO和WS₂作为ETL的器件具有优越的性能。
交叉验证结果显示，wxAMPS软件模拟得到的性能参数（PCE等）普遍略高于SCAPS-1D的结果（例如，ITO/TiO₂/CsPbI₃/CBTS/Au结构在wxAMPS中PCE为19.75%），但两者的趋势完全一致，最佳材料组合的排序也相同，这有力地支持了SCAPS-1D模拟结果的合理性和可靠性。
本研究得出的核心结论是：通过结合DFT计算和SCAPS-1D/wxAMPS数值模拟，成功地系统评估了CsPbI3的材料性质，并从大量可能的ETL/HTL组合中，高效地筛选出有望实现高效率的无机PSC器件结构。其中，以无机、地壳储量丰富的铜钡硫代锡酸盐（CBTS）作为HTL，与传统的TiO₂或ZnO等ETL搭配，构成了性能优异的器件构型（尤其是ITO/TiO₂/CsPbI₃/CBTS/Au）。该结构模拟效率达到17.90%，并且通过参数优化显示其有潜力达到18.4%以上。研究还全面评估了该器件对厚度、电阻和温度的敏感性，为其实际制备提供了重要的设计参考。最终，该模拟研究为实验上制备低成本、高稳定性和高效率的CsPbI3基全无机PSC指明了有前途的材料组合方向，避免了耗时的实验试错。
本研究的亮点和创新性在于：1) 系统性：首次同时将如此多种类的ETL（8种）和HTL（12种）与CsPbI3进行组合筛选，覆盖面广，结论更具普适性和参考价值。2) 多尺度模拟结合：将原子尺度的DFT计算（揭示材料本征特性）与器件尺度的物理模拟（预测宏观性能）有机结合，形成了从材料到器件的完整研究链条。3) 高效的研究范式：展示了利用先进模拟软件快速、低成本地筛选高性能光伏器件材料组合的强大能力，为新材料和新结构的探索提供了高效的研究范式。4) 突出了新型HTL材料：重点推介了CBTS这种相对较新、具有可调带隙和良好光吸收能力的无机HTL材料，并证实了其与CsPbI3搭配的优越性，为替代昂贵或不稳定的有机HTL提供了新选择。5) 深入的分析与验证：不仅给出了最佳组合，还对其进行了全面的参数敏感性分析和内部物理过程剖析，并使用了两种独立模拟软件进行结果比对，增强了研究的严谨性和可信度。这项工作的价值不仅在于为CsPbI3基PSCs找到了有潜力的优化路径，其采用的研究方法也为其他新型光伏材料的器件设计与性能预测提供了可借鉴的模板。作者在文末建议，未来可以将机器学习引入此类研究，以进一步加速最优材料组合的搜索过程。