钙钛矿建模Word格式.docx

钙钛矿建模Word格式.docx

《钙钛矿建模Word格式.docx》由会员分享，可在线阅读，更多相关《钙钛矿建模Word格式.docx（20页珍藏版）》请在冰豆网上搜索。

钙钛矿太阳能电池结构由掺杂氟SnO2（fluorine-oxide,FTO）导电玻璃、电子传输层（ETM）、钙钛矿吸收层（如CH3NH3PbI3、CH3NH3PbCl3、CH3NH3PbBr3等）、空穴传输层（HTM）和金属对电极组成。

电子传输层（ETM）多为ZnO、TiO2等，空穴传输层（HTM）多为Spiro-OMeTad、FTAA、H3MT、PEDOT、PASS等固态介质材料。

图1.1钙钛矿太阳能电池结构示意图

当钙钛矿吸收层吸收阳光受激后，产生电子空穴对，激子在吸收层运动至ETM/钙钛矿吸收层/HTM界面后发生分离，电子注入ETM层（阳极），空穴注入到HTM（阴极），最后经外部电路循环在金属对电极复合形成回路电流。

图1.2钙钛矿太阳能电池原理示意图

1.2CH3NH3PbX3钙钛矿的结构和性能

理想钙钦矿结构的化学式为ABX3,属pm-3m空间群,整个结构可以描述为BX6八面体顶点相连形成钙钦矿结构的基本骨架,A位离子填充于所有形成的12配位的空隙中,也可看作A离子与X离子形成立方密堆结构,B位离子填充于八面体空隙中"

这种结构具有很强的A一X和B一键平衡键长间的不匹配容忍度,可以通过BX6八面体扭转及A或B离子的位移转化为其他低对称空间群结构。

而其中的有机金属卤化物钙钛矿作为敏化剂近年来被用于研究敏化太阳电池，得到飞速发展。

其化学式为CH3NH3PbX3，其中X=Cl，Br或I。

其中A=CH3NH3，B=Pb，X=Cl，Br或I。

虽然有机无机杂化钙钛矿的晶体结构和原始钙钛矿的结构类似，但取代其晶胞位置的不是一个原子，而是一个有机的原子团。

原子团中含有C-H和N-H键，这些键不够稳定，导致其空间结构也容易倾斜发生改变，而晶体结构的改变伴随着能量的释放，这就是有机无机钙钛矿的导电原理。

CH3NH3PbX3钙钛矿作为新型染料敏化太阳能电池的关键原料，自2012以来，在国内外太阳能电池材料领域成为新的、重点研究的方向。

而其中CH3NH3PbX3钙钛矿的的晶体结构和光电转换效率是其中的重中之重。

CH3NH3PbX3钙钛矿太阳电池不仅具有较高的能量转换效率，而且其核心光电转换材料具有廉价、可溶液制备的特点，便于采用不需要真空条件的卷对卷技术制备，这为CH3NH3PbX3钙钛矿太阳电池的大规模、低成本制造提供可能。

不仅如此，CH3NH3PbX3钙钛矿太阳电池还可以制备在柔性衬底上，便于应用在各种柔性电子产品中，例如可穿戴的电子设备、折叠式军用帐篷等。

与传统的染料敏化太阳电池相比，CH3NH3PbX3钙钛矿太阳电池不需要液体电解质，不用担心太阳电池的漏液问题。

与有机光伏器件相比，CH3NH3PbX3钙钛矿太阳电池的核心光电转换材料是有机-无机杂化材料，材料的耐候性可能会优于有机光伏器件中使用的有机半导体材料。

这些优点可能会使钙钛矿太阳电池在实际使用中具有比染料敏化太阳电池和有机光伏器件更好的性能稳定性和更长的使用寿命[]。

基于上述原因，钙钛矿太阳电池具有非常光明的产业化前景，是现有商业太阳电池最有潜力的竞争者。

1.3选题的意义

随着节能减排的进行，倡导发展可持续社会，发展和利用新能源都是重中之重。

太阳能作为一种安全、无污染、可再生等优点，具有其特殊的优势。

但是目前太阳能电池由于相对煤炭发电成本高，以及转换效率低等缺点限制其发展规模。

所以研究和发现新的材料作为太阳能电池具有重大意义。

钙钛矿近年来作为一种热门的材料，利用当下先进的计算机模拟软件对、相应的钙钛矿材料进行模拟分析，得知它的晶格参数、性能及各种深入研究有很大帮助。

二、CH3NH3PBI3钙钛矿结构的构建

采用基于第一性原理密度泛函理论的平面波赝势方法，对CH3NH3PbI3钙钛矿晶体进行了优化，并对其能带结构，态密度等进行了理论计算。

理论计算主要用到了MS软件中的CASTEP（CambridgeSerialTotalEngeryPackage）模块进行能量计算。

CASTEP是基于密度泛函理论的从头算量子力学程序，利用总能量平面波赝势方法，将离子势用赝势代替，电子波函数通过平面波基组展开，电子-电子的交换和关联作用由局域密度近似（localdensityapproximation,LDA）或广义梯度近似（generalgradientapproximation,GGA）校正，是目前材料科学计算中使用非常广泛的电子结构计算方法。

经过查阅CH3NH3PbI3钙钛矿相关资料，利用MaterialStudio对CH3NH3PbI3钙钛矿进行了建模。

2.1CH3NH3PbI3钙钛矿的主要参数

利用晶体的空间群、晶格参数以及各原子的坐标建立CH3NH3PbI3钙钛矿几何模型。

具体如下表：

表2.1CH3NH3PbI3钙钛矿相关参数

空间群

P1，空间群代号为1

Latticetype

3DTriclinic

晶格参数

a=8.8埃b=8.8埃c=13.0475埃

2.2建立几何模型

2.2.1建立晶格

首先进入软件，建立3DAtomisticdocument并建立晶体，选择晶体族群1P1，在建立集体空间群之后，接下来输入晶胞的参数，根据之前的参数输入如下图。

图2.1晶胞参数设置图

2.2.2添加原子

通过Build-AddAtoms添加相应的I，Pb原子和CH3NH3基团，从而形成CH3NH3PbI3结构，虽热只加入了一个I，一个Pb，一个CH3NH3基团，但群的对称操作在晶体中补充了剩余原子。

这里，所谓的原子分数坐标是由晶胞所含的结构基元及其在晶胞中的位置决定的。

依次添加完原子后即可得到CH3NH3PbI3钙钛矿几何结构，如下图

图2.2CH3NH3PbI3钙钛矿几何结构图

2.3利用Reflex模块进行粉末衍射图谱

我们可以知道任何晶体都是由原子或者原子团在三维空间按相应的空间集合结构组成，所以不同的晶体就会有不同的X射线衍射图。

利用MS自带的Reflex模块可以为我们定性、定量的分析的同时，也可以帮我们验证我们的模型是不是准确的，进而进一步研究它的结构。

通过如下图3.4选择Reflex模块下的PowerDiffraction即可得到CH3NH3PbI3钙钛矿几何结构X射线粉末衍射图如图3.5：

图2.2粉末衍射设置图

图2.3粉末衍射图

从图中我们可以看到衍射峰强度主要有四处较突出，可绘出下表

表2.3衍射峰强度以及相应的角度

衍射峰的衍射角

衍射峰强度

13.56

81.79

14.32

99.87

23.67

16.89

24.63

20.82

三、对晶体进行优化结构

3.1结构优化

在CASTEP模块下Calculation时会对给出的晶体结构进行优化，优化的过程是通过不断的迭代计算调整原子的坐标、晶胞的参数使结构的总能量最小化的。

通过迭代计算不断减少应力数量级的同时，也使得晶胞内部的应力张量和所施加的外部应力相等，也更加接近真实的材料的结构。

3.1.1参数设置

通过在菜单栏下Modules下的CASTEP中选择calculation，在下拉菜单栏中选择GeometryOptional,选择精度为coarse，在Electronic栏中点击More,在里面设置charge为0.15。

点击Run运行，就会出现相应的工作进度，在运行时以及完成后，我们可以看到会有文件生成，有代表优化结构步骤的图、优化过程记录运行状态的文档、优化后的结构图以及包含能量、压力、应力的图等。

3.1.2运行结果显示

在运行晶体结构优化后，我们可以看到里面文件包含的关于能量收敛的图如下，我们可以看到里面包含着一些能量的变化、外界压力、迭代次数位移的变化等信息。

（a）

（b）

图3.1优化过程能量图

从图中我们可以发现，在优化的过程中表现了收敛，在文档中我们发现图下图的一部分数据也表明我们的优化过程是收敛。

图3.2文件数据

3.2能量计算

经过几何优化后，我们得到了较为真实可靠的ch3nh3pbi3型钙钛矿晶体结构，接下来我们继续利用CASTEP模块继续对优化后的结构进行能量计算，同时也可以在分析的过程中出现态密度图。

3.2.1参数设置

如下图所示，我们选择Module中的CASTEP中的calculation，但是此时我们应该选择下拉单中的Energy进行能量计算。

图3.3能量计算参数设置

3.2.2运行结构显示

我们也可以对利用CASTEP中的Analysis进行分析获得态密度图和能带结构。

通过勾选Analysis对话框中的”ShowDos”即可在一个途中显示出来，如下图：

图3.4态密度和能带结构图

从图中我们可以的到：

1、能带禁带宽度为1.441ev,符合钙钛矿太阳电池吸收层光学性能的要求。

2、根据我们所学的半导体知识，价带和导带相交时，为金属。

如果不想交久为半导体或者绝缘体。

从图中，导带是高对称点附近类似开口向上的抛物线形状的能带，我们发现导带和价带并没有相交，而且Eg并没有那么大，所以ph3nh3pbi3是半导体，这也符合我们把半导体作为太阳电池材料。

在经过能量计算运行结束后，我们可以看到如下两个图，

图3.5能量收敛图

从图中可以看到随着计算次数的不断增加，能量在某一值附近，即对晶体进行的计算是收敛的。

3.3计算结果分析

根据运行完成后的文件，我们可以打开后缀名为CASTEP的文件，里面包含了很多能量的相关结论。

如下表举例所示，每经过一次迭代后，都会显示相应次数后的总能量：

表3.1迭代后总能量

迭代次数

CASTEP文件给出的能量

14

Convergedin14iterationstoatotalenergyof-12.4182eV

16

Convergedin16iterationstoatotalenergyof-311.6138eV

20

Convergedin20iterationstoatotalenergyof-1655.7917eV

21

Convergedin21iterationstoatotalenergyof-262.5918eV

最终我们也可以从CASTEP文件中得出能量结果，如下表

表3.2CH3NH3PBI3能量结果表

ElectronicMinimizationParameters

Method

Treatingsystemasmetallicwithdensitymixingtreatmentofelectrons

totalenergy/atomconvergencetol

0.2000E-05eV

eigen-energyconvergencetolerance

0.7934E-06eV

Fermienergyconv