钙钛矿建模.docx

1、钙钛矿建模利用第一性原理对甲胺碘铅进行结构模拟及性能研究摘要钙钛矿型甲胺碘铅薄膜太阳能电池因其结构简单、制备成本低廉等优点吸引了科研工作者的关注。其光电转换效率在近5年内从%迅速提高到15%以上，高于非晶硅的太阳电池效率，被Science评选为2013年十大科学突破之一。本文主要是用Material Studio通过对CH3NH3PbI3钙钛矿建立晶体模型后，利用第一性原理对有机金属钙钛矿CH3NH3PbI3及其延伸结构进行结构模拟、能带分析和性能研究。关键词：第一性原理 ; CH3NH3PbI3钙钛矿; 能带结构一、 绪论钙钛矿太阳电池的结构与原理钙钛矿是指CaTiO3,化学式是CaTiO3

2、，是一种立方晶系的氧化物，是俄国化学家佩罗夫斯基发现的。钙钛矿太阳能电池结构由掺杂氟SnO2(fluorine-oxide,FTO)导电玻璃、电子传输层（ETM）、钙钛矿吸收层（如CH3NH3PbI3、CH3NH3PbCl3、CH3NH3PbBr3等）、空穴传输层（HTM）和金属对电极组成。电子传输层（ETM）多为ZnO、TiO2等，空穴传输层（HTM）多为Spiro-OMeTad、FTAA、H3MT、PEDOT、PASS等固态介质材料。图 钙钛矿太阳能电池结构示意图当钙钛矿吸收层吸收阳光受激后，产生电子空穴对，激子在吸收层运动至ETM/钙钛矿吸收层/HTM界面后发生分离，电子注入ETM层（阳

3、极），空穴注入到HTM(阴极)，最后经外部电路循环在金属对电极复合形成回路电流。图 钙钛矿太阳能电池原理示意图钙钛矿的结构和性能理想钙钦矿结构的化学式为ABX3,属pm-3m空间群,整个结构可以描述为BX6八面体顶点相连形成钙钦矿结构的基本骨架,A位离子填充于所有形成的12配位的空隙中,也可看作A离子与X离子形成立方密堆结构,B位离子填充于八面体空隙中这种结构具有很强的A一X和B一键平衡键长间的不匹配容忍度,可以通过BX6八面体扭转及A或B离子的位移转化为其他低对称空间群结构。而其中的有机金属卤化物钙钛矿作为敏化剂近年来被用于研究敏化太阳电池，得到飞速发展。其化学式为CH3NH3PbX3，其中

4、X=Cl，Br或I。其中A=CH3NH3，B=Pb，X=Cl，Br或I。虽然有机无机杂化钙钛矿的晶体结构和原始钙钛矿的结构类似，但取代其晶胞位置的不是一个原子，而是一个有机的原子团。原子团中含有C-H和N-H键，这些键不够稳定，导致其空间结构也容易倾斜发生改变，而晶体结构的改变伴随着能量的释放，这就是有机无机钙钛矿的导电原理。CH3NH3PbX3钙钛矿作为新型染料敏化太阳能电池的关键原料，自2012以来，在国内外太阳能电池材料领域成为新的、重点研究的方向。而其中CH3NH3PbX3钙钛矿的的晶体结构和光电转换效率是其中的重中之重。CH3NH3PbX3钙钛矿太阳电池不仅具有较高的能量转换效率，而

5、且其核心光电转换材料具有廉价、可溶液制备的特点，便于采用不需要真空条件的卷对卷技术制备，这为CH3NH3PbX3钙钛矿太阳电池的大规模、低成本制造提供可能。不仅如此，CH3NH3PbX3钙钛矿太阳电池还可以制备在柔性衬底上，便于应用在各种柔性电子产品中，例如可穿戴的电子设备、折叠式军用帐篷等。与传统的染料敏化太阳电池相比，CH3NH3PbX3钙钛矿太阳电池不需要液体电解质，不用担心太阳电池的漏液问题。与有机光伏器件相比，CH3NH3PbX3钙钛矿太阳电池的核心光电转换材料是有机-无机杂化材料，材料的耐候性可能会优于有机光伏器件中使用的有机半导体材料。这些优点可能会使钙钛矿太阳电池在实际使用中具

6、有比染料敏化太阳电池和有机光伏器件更好的性能稳定性和更长的使用寿命。基于上述原因，钙钛矿太阳电池具有非常光明的产业化前景，是现有商业太阳电池最有潜力的竞争者。选题的意义随着节能减排的进行，倡导发展可持续社会，发展和利用新能源都是重中之重。太阳能作为一种安全、无污染、可再生等优点，具有其特殊的优势。但是目前太阳能电池由于相对煤炭发电成本高，以及转换效率低等缺点限制其发展规模。所以研究和发现新的材料作为太阳能电池具有重大意义。钙钛矿近年来作为一种热门的材料，利用当下先进的计算机模拟软件对、相应的钙钛矿材料进行模拟分析，得知它的晶格参数、性能及各种深入研究有很大帮助。二、CH3NH3PBI3钙钛矿结

7、构的构建采用基于第一性原理密度泛函理论的平面波赝势方法，对CH3NH3PbI3钙钛矿晶体进行了优化，并对其能带结构，态密度等进行了理论计算。理论计算主要用到了MS软件中的CASTEP(Cambridge Serial Total Engery Package)模块进行能量计算。CASTEP是基于密度泛函理论的从头算量子力学程序，利用总能量平面波赝势方法，将离子势用赝势代替，电子波函数通过平面波基组展开，电子-电子的交换和关联作用由局域密度近似（local density approximation,LDA）或广义梯度近似（general gradient approximation,GGA）校

8、正，是目前材料科学计算中使用非常广泛的电子结构计算方法。 经过查阅CH3NH3PbI3钙钛矿相关资料，利用Material Studio对CH3NH3PbI3钙钛矿进行了建模。钙钛矿的主要参数利用晶体的空间群、晶格参数以及各原子的坐标建立CH3NH3PbI3钙钛矿几何模型。具体如下表：表钙钛矿相关参数空间群P1， 空间群代号为1Lattice type3D Triclinic晶格参数a=埃 b=埃 c=埃建立几何模型建立晶格 首先进入软件，建立3D Atomistic document并建立晶体，选择晶体族群1 P1， 在建立集体空间群之后，接下来输入晶胞的参数，根据之前的参数输入如下图。图晶

9、胞参数设置图添加原子通过Build-Add Atoms添加相应的I，Pb原子和CH3NH3基团，从而形成CH3NH3PbI3结构，虽热只加入了一个I，一个Pb，一个CH3NH3基团，但群的对称操作在晶体中补充了剩余原子。这里，所谓的原子分数坐标是由晶胞所含的结构基元及其在晶胞中的位置决定的。 依次添加完原子后即可得到CH3NH3PbI3钙钛矿几何结构，如下图图钙钛矿几何结构图利用Reflex模块进行粉末衍射图谱我们可以知道任何晶体都是由原子或者原子团在三维空间按相应的空间集合结构组成，所以不同的晶体就会有不同的X射线衍射图。利用MS自带的Reflex模块可以为我们定性、定量的分析的同时，也可以

10、帮我们验证我们的模型是不是准确的，进而进一步研究它的结构。通过如下图选择Reflex模块下的Power Diffraction 即可得到CH3NH3PbI3钙钛矿几何结构X射线粉末衍射图如图：图粉末衍射设置图图粉末衍射图从图中我们可以看到衍射峰强度主要有四处较突出，可绘出下表表衍射峰强度以及相应的角度衍射峰的衍射角衍射峰强度三、对晶体进行优化结构结构优化在CASTEP模块下Calculation时会对给出的晶体结构进行优化，优化的过程是通过不断的迭代计算调整原子的坐标、晶胞的参数使结构的总能量最小化的。通过迭代计算不断减少应力数量级的同时，也使得晶胞内部的应力张量和所施加的外部应力相等，也更加

11、接近真实的材料的结构。参数设置通过在菜单栏下Modules下的CASTEP 中选择calculation，在下拉菜单栏中选择 Geometry Optional ,选择精度为coarse，在Electronic栏中点击More ,在里面设置charge为。点击Run 运行，就会出现相应的工作进度，在运行时以及完成后，我们可以看到会有文件生成，有代表优化结构步骤的图、优化过程记录运行状态的文档、优化后的结构图以及包含能量、压力、应力的图等。运行结果显示在运行晶体结构优化后，我们可以看到里面文件包含的关于能量收敛的图如下，我们可以看到里面包含着一些能量的变化、外界压力、迭代次数位移的变化等信息。(

12、a)(b)图优化过程能量图从图中我们可以发现，在优化的过程中表现了收敛，在文档中我们发现图下图的一部分数据也表明我们的优化过程是收敛。图 文件数据能量计算经过几何优化后，我们得到了较为真实可靠的ch3nh3pbi3型钙钛矿晶体结构，接下来我们继续利用CASTEP模块继续对优化后的结构进行能量计算，同时也可以在分析的过程中出现态密度图。参数设置如下图所示，我们选择Module 中的CASTEP中的calculation ，但是此时我们应该选择下拉单中的Energy 进行能量计算。图能量计算参数设置运行结构显示我们也可以对利用CASTEP中的Analysis 进行分析获得态密度图和能带结构。通过勾

13、选Analysis对话框中的”Show Dos”即可在一个途中显示出来，如下图：图态密度和能带结构图从图中我们可以的到：1、能带禁带宽度为,符合钙钛矿太阳电池吸收层光学性能的要求。2、根据我们所学的半导体知识，价带和导带相交时，为金属。如果不想交久为半导体或者绝缘体。从图中，导带是高对称点附近类似开口向上的抛物线形状的能带，我们发现导带和价带并没有相交，而且Eg并没有那么大，所以ph3nh3pbi3是半导体，这也符合我们把半导体作为太阳电池材料。在经过能量计算运行结束后，我们可以看到如下两个图，（a）（b）图能量收敛图从图中可以看到随着计算次数的不断增加，能量在某一值附近，即对晶体进行的计算是

14、收敛的。计算结果分析根据运行完成后的文件，我们可以打开后缀名为CASTEP的文件，里面包含了很多能量的相关结论。如下表举例所示，每经过一次迭代后，都会显示相应次数后的总能量：表迭代后总能量迭代次数CASTEP文件给出的能量14Converged in 14 iterations to a total energy of eV16Converged in 16 iterations to a total energy of eV20Converged in 20 iterations to a total energy of eV21Converged in 21 iterations to a

15、 total energy of eV最终我们也可以从CASTEP文件中得出能量结果，如下表表 CH3NH3PBI3能量结果表Electronic Minimization ParametersMethodTreating system as metallic with density mixing treatment of electronstotal energy / atom convergence toleigen-energy convergence toleranceFermi energy convergence tolerance从上表中，我们可以得到总能量是，本征吸收层收敛容

16、差是，费米能量收敛容差是。得到的能量数据可以更好的帮助我们了解CH3NH3PBI3钙钛矿材料，同时也为下面和参杂、替换相应的钙钛矿做对比，具有很大的意义。四、CH3NH3PbI3钙钛矿几何结构延伸正如前面所说的理想钙钦矿结构的化学式为ABX3，其中 A、B 是指金属元素，X 是指非金属元素。大约有超过 20 多种元素可以占据 A位置，例如 Ba、Sr、Pb、K、Sc、Y 以及镧系元素等；而几乎有50 余种元素可以占据 B 位置；此外，能占据 X 位置的元素不止是O，卤族元素 F、Cl、Br、I 等以及其他一些非金属元素 N 等都可以占据 X 位置形成钙钛矿结构。可以利用B离子代替B离子会形成A

17、BxB1-XX3型固溶体化合物,最为常见的B位离子无序化合物为铅基钙钦矿固溶体化合物,如节，把CH3NH3PbI3的Pb替换成为Sn，形成CH3NH3SnI3钙钛矿结构。同理也可以利用X离子代替X离子会形成ABXxX3-X型固溶体化合物,按照钙钛矿比例去替换相应原子，如节，把I离子与Cl离子按1:2参杂，形成CH3NH3PbI2Cl钙钛矿。用软件模拟分析相应的晶体结构。的掺杂对钙钛矿电池的能带结构的影响通过MS软件中Edit中Select Atom 选择Pb原子，再利用Modify中的Modify Element中选择用Sn替换Pb原子，得到如下图所示的CH3NH3SnI3钙钛矿结构，进而来研

18、究CH3NH3SnI3钙钛矿几何结构的性能、能量及其性质。图钙钛矿几何结构粉末衍射图我们依然用MS软件中的Reflex模块进行粉末衍射图如下：图粉末衍射图从图中我们可以看到衍射峰强度主要有两处较突出，可绘出下表表衍射峰强度以及相应的角度衍射峰的衍射角衍射峰强度几何优化几何优化的运行结果如下图所示：（a)(b)图优化过程能量图能量计算对利用CASTEP中的Analysis 进行分析获得态密度图和能带结构。通过勾选Analysis对话框中的”Show Dos”即可在一个途中显示出来，如下图：图态密度和能带结构图从图中可知其的禁带宽度为,能量运行结束后，得到如下两个图：（a)(b)图能量收敛图从图中

19、可以看到随着计算次数的不断增加，能量在某一值附近，即对晶体进行的计算是收敛的。计算结果分析同样根据运行完成后的文件，我们可以打开后缀名为CASTEP的文件，里面包含了很多能量的相关结论。如下表举例所示，每经过一次迭代后，都会显示相应次数后的总能量：表迭代后总能量迭代次数CASTEP文件给出的能量14Converged in 14 iterations to a total energy of eV17Converged in 17 iterations to a total energy of eV21Converged in 21 iterations to a total energy o

20、f eV最终我们也可以从CASTEP文件中得出能量结果，如下表表 CH3NH3SNI3能量结果表Electronic Minimization ParametersMethodTreating system as metallic with density mixing treatment of electronstotal energy / atom convergence toleigen-energy convergence toleranceFermi energy convergence tolerance从上表中，我们可以得到总能量是，本征吸收层收敛容差是，费米能量收敛容差是。掺杂

21、CH3NH3PbI3钙钛矿的几何机构及其分析通过MS软件中 选择四个I原子，再利用Modify中的Modify Element中选择用Cl替换I原子，得到如下图所示的CH3NH3PbI2Cl钙钛矿结构。图钙钛矿几何结构粉末衍射图我们依然用MS软件中的Reflex模块进行粉末衍射图如下：图粉末衍射图从图中我们可以看到衍射峰强度主要有三处较突出，可绘出下表表衍射峰强度以及相应的角度衍射峰的衍射角衍射峰强度几何优化几何优化的运行结果如下图所示：（a)(b)图优化过程能量图能量计算对利用CASTEP中的Analysis 进行分析获得态密度图和能带结构。通过勾选Analysis对话框中的”Show Do

22、s”即可在一个途中显示出来，如下图：图态密度和能带结构图从图中可知其的禁带宽度为,能量运行结束后，得到如下两个图：（a)(b)图能量收敛图从图中可以看到随着计算次数的不断增加，能量在某一值附近，即对晶体进行的计算是收敛的。计算结果分析同样根据运行完成后的文件，我们可以打开后缀名为CASTEP的文件，里面包含了很多能量的相关结论。如下表举例所示，每经过一次迭代后，都会显示相应次数后的总能量：表迭代后总能量迭代次数CASTEP文件给出的能量16Converged in 16 iterations to a total energy of eV19Converged in 19 iterations

23、 to a total energy of eV20Converged in 20 iterations to a total energy of eV21Converged in 21 iterations to a total energy of eV最终我们也可以从CASTEP文件中得出能量结果，如下表表 CH3NH3PbI2Cl1能量结果表Electronic Minimization ParametersMethodTreating system as metallic with density mixing treatment of electronstotal energy /

24、atom convergence toleigen-energy convergence toleranceFermi energy convergence tolerance从上表中，我们可以得到总能量是，本征吸收层收敛容差是，费米能量收敛容差是。五、结论（1）CH3NH3BX3型钙钛矿材料很多，可以通过替换、参杂等变换得到，但是他们的基本结构性质差不多。（2）我们可以发现以上模拟分析的三种除去本征吸收层收敛容差的能量是不相同的之外，总能量是相同的，费米能量收敛容差也是相同的。（3）在得到结构的X射线粉末衍射图时，我们也可以通过做实验的到相应的钙钛矿材料得出相应的图谱，进行验证，获得更精准的结论。在得到的结构上得到的结论也说明我们的计算是收敛的，得到的相应结果也说明其相应能量、光学属性。六、致谢时间飞逝，转眼大学生活即将结束，在大学不仅学到了专业知识，也学会了成长。大学期间也得到了很多老师朋友同学的帮助，在此向他们表示衷心的感谢！首先感谢我的父母，在经济和精神上给我帮助和支持。其次，感谢我的论文指导老师钟炜对我的教导和帮助。也感谢我的班主任章娟四年来的关心和帮助。最后感谢老师和同学们的指导和关心。谢谢！参考文献