张思(1985-), 女, 硕士研究生. E-mail:zhangsi0825@163.com
采用第一性原理方法, 模拟了应变对四方相BaTiO3(BTO)空位形成能的影响, 并计算了应变和双空位共同作用下BaTiO3体系的极化. 计算结果表明, 应变改变了空位形成能, 且在-3%到3%的应变范围内, 氧空位的形成能逐渐减小, 亦即压应变抑制了氧空位的形成, 这与Ti-O键的变化有关. 此外压应变在改变双空位形成能的同时, 增大了体系的极化和极化偏移, 对于1A和2A组合的双空位即使在-3%的应变下其极化偏移也只有1 μC/cm2和 0.3 μC/cm2, 而极化偏移是印记效应的表现, 因此应变对印记效应的影响很小, 可以忽略.
The first-principle methods are used to calculate the effects of in-plane strain on the formation energies of neutral vacancies. It is shown that the formation energies of neutral O vacancy increase with compressive strain, while decrease for tensile strain. In other words, the O vacancy can be suppressed by compressive strain, which is contributed to the variation of the Ti-O bond. The formation energies of Ba-O divacancies and polarization of the tetragonal BaTiO3 are studied under different strain conditions. The strain changes the formation energies of Ba-O divacancies and the polarization of the systems containing Ba-O divacancies. More over, the results show that the polarization offset increases gradually when the strain changes from tensile to compressive strain. In other words the polarization imprint of BaTiO3 increase under compressive strain, but the changes are slight and negligible.
作为一种典型的钙钛矿结构铁电材料, 钛酸钡在室温下具有四方相结构, 且自发极化可在外电场作用下翻转. 利用这一特性可以制造性能优越的铁电存储器(FeRAM). 由于铁电材料在微电子领域的广泛应用, 使得铁电薄膜成为目前研究的主要方向之一[ 1].
在薄膜制备过程中, 由于铁电薄膜和衬底之间的晶格差异以及热膨胀系数的不同, 会在薄膜中引入应变, 使得铁电氧化物薄膜的性质与块体材料有很大差异. 同时, 应变与极化之间存在强烈的耦合作用, 使得应变对体系结构有很大影响, 甚至引起相变[ 2]. 此外, 应变对铁电材料的相变温度、压电和热电效应等性能参数有着明显影响, 并在理论和实验研究中受到广泛的关注. 例如室温下应变能够在顺电相SrTiO3中引发铁电相变[ 3]. 当沉积在DyScO3上的BaTiO3薄膜的双轴向压应变达到1.3%时, 其相变温度提高500℃, 剩余极化也至少提高2.5倍以上[ 4]. Ederer等[ 5]通过计算研究了应变对铁电材料极化的影响, 结果显示, 压应变增大体系的自发极化, 拉应变减小体系的自发极化.
与此同时, 在薄膜和器件制备过程中不可避免会引入缺陷, 而晶格缺陷的产生与铁电器件的失效行为密切相关[ 6], 例如氧空位的有序化是引起BaTiO3晶体中铁电疲劳的重要原因[ 7]. 此外, 印记效应也是铁电失效的一种重要表现形式, 印记效应是指电滞回线沿电压轴的偏移, 进而导致无法正确区分上下极化状态[ 8]. 目前解释印记效应有两种模型: 偶极矩模型[ 9]和屏蔽模型[ 10]. 同时, 对于立方相BaTiO3缺陷的计算也有不少报道. Moriwake计算了中性空位的形成能[ 11], Lee等[ 12]在此基础上计算了不同离子价态的形成能. 由于缺陷和应变对铁电器件性能都有较大的影响, 因而研究应变对缺陷的影响以及应变和缺陷对体系极化的影响具有重要意义. 基于以上原因, 本工作利用第一性原理方法, 通过改变四方相BaTiO3晶格常数来模拟应变, 分析了应变对中性空位和钡氧双空位形成能的影响以及应变和钡氧双空位共同作用下体系极化的变化.
以四方相BaTiO3作为研究对象, 采用基于密度泛函理论的第一性原理方法, 利用广义梯度近似(GGA)处理电子之间的交换关联能[ 13, 14]. 考虑5s25p66s2、3s23p64s23d2、2s22p4电子作为Ba、Ti、O的价电子, 其他电子则作为芯电子. 由于PBEsol泛函计算钙钛矿结构铁电材料的理论结果更接近实验值[ 15], 采用PBEsol泛函进行计算, 并用共轭梯度(CG)算法对原胞体积和离子位置进行弛豫. 首先对体系进行收敛测试, 结果表明采用500 eV和8×8×8的 k点网格已经达到精度要求. 计算得到四方相BaTiO3晶格常数 a=0.3974 nm、 c=0.4055 nm, 结果列于表1, 并与实验和其他理论计算结果进行了比较. Roman等[ 15]的结果表明采用PBEsol泛函计算结果与实验值符合得更好, 而本课题组的计算结果与其基本一致. 这为下一步采用PBEsol泛函得到的晶格结构进行形成能计算提供了依据.
![]() | 表1 四方相BaTiO3晶格常数、原子位置沿 z方向的偏移量、体积和自发极化 Table 1 Structural properties and macroscopic polarization of the tetragonal BaTiO3 |
利用改变 a、 b轴晶格常数并弛豫 c轴晶格常数的方法来模拟平面应变, 空位则通过移除格点上的原子来实现, 其形成能定义式采用Lee等[ 12]的形式.对于中性氧空位的形成能由于其接近导带底, 采用Padilla等的方法对其进行修正[ 19, 20]. 各原子的化学势可以通过Ba-Ti-O体系的“相图”[ 12]及求解相图所对应的各材料的低温稳定相的总能得到[ 20, 21].
VBa- VO双空位通过移除格点上的Ba离子和O离子得到, 并采用Zhang和Northrup[ 22]的方法来计算 VBa- VO的形成能. 由于第二近邻 VBa- VO双空位的形成能在几种组合中最低[ 23], 本文采用PBEsol利用计算得到的晶格常数计算了应变对第二近邻 VBa- VO形成能的影响. 在此基础上利用现代极化理论计算含有 VBa- VO双空位体系的极化, 并把移除了Ba、O离子的四方中心对称结构作为初始构型, 通过Berry相方法计算体系的极化值[ 23].
图1为BTO自发极化随应变的变化, 图中直线利用Ederer等的极化增量与应变关系得到[ 5], 其公式为:
![]() | (1) |
其中 c33、 c31和 ceff为压电常数, n为泊松比. 通过比较可以发现, 本文结果与Ederer等的计算结果基本一致. 如图2所示, c轴应变与 ab轴平面应变呈线性变化关系, 即在平面应变从压应变到拉应变的范围内, c轴应变由拉应变逐渐变为压应变; 而体积则在-3%到3%的应变范围内呈线性增大.
![]() | 图1 BTO自发极化随应变的变化, 直线由式(1)得到Fig. 1 Variation of spontaneous polarizations with in-plane strain, line is calculated using Eq.(1) |
利用计算得到的晶格常数对空位形成能进行计算. 空位的形成能与原子化学势有关, 在不同的相平衡点空位形成能也不同. 图3给出了不同相平衡点的空位形成能, 在氧化环境下Ti空位的形成能最低是主要缺陷; 而从A点到G点即从氧化环境到还原环境的过程中, O空位逐渐成为主要缺陷, 形成能不断降低. 其变化趋势与Lee等[ 12]立方相的结果一致, 且空位形成能的计算结果与采用PBE泛函计算得到的结果基本一致[ 21].
为了模拟应变下空位形成能的变化, 对体系施加-3%到3%的平面应变, 其定义为 εs=(a1-a0)/a0, 其中 a0、 a1表示应变前后 a、 b晶格常数, εs为体系应变. 图4给出了应变下 VBa、 VO、 VTi空位相平衡点E点的形成能, 结果表明在压应变到拉应变的整个变化范围内 VO形成能逐渐增小, 与体积随应变变化趋势相反. 换言之 VO形成能随体积的增大而减小, 随体积的减小而增大; 而 VBa形成能在应变作用下的变化, 具有和 VO相反的变化趋势; VTi的形成能从拉应变到1%的压应变变化范围内逐渐增大, 在-1%处达到最大值, 当压应变继续增大其形成能开始下降.
![]() | 图4 Ba、Ti和O中性空位形成能Fig. 4 Variation of neutral vacancy formation energy of Ba, Ti and O with in-plane strain |
三种空位形成能随应变的不同变化趋势与BaTiO3中三种原子的结合方式有关, 在BaTiO3中Ti与O原子存在较强的共价杂化作用, 在压应变下 c轴方向的Ti-O作用不断加强, 使得O空位更难脱离体系. 而Ti与下方O原子的作用不断减弱, 当压应变达到1%时已基本无共价杂化作用, 所以其形成能在-1%的应变下开始下降. 对于Ba在Ti-O之间共价杂化不断加强的同时, Ba-O之间作用减弱, 使得在压应变下Ba的形成能随压应变的增大而减小. 本文的计算结果表明, 可以利用压应变来抑制氧空位的形成.
在3×3×3的超胞格点中移除Ba离子和O离子来模拟 VBa- VO双空位, 从而计算其形成能. 已有研究结果表明, 对于不同组合 VBa- VO, 第二近邻的形成能最低[ 23]. 基于此, 本文中只计算了应变对第二近邻 VBa- VO双空位缺陷形成能的影响. 图5给出了应变下 VBa- VO双空位的形成能, 其中1A、1B、2A和2B分别对应第二近邻 VBa- VO双空位的四种不同组合. 在拉应变到压应变的变化范围内, 1B与1A形成能的变化趋势基本相同, 都是随拉应变的增大而减小, 随压应变的增大先增大后减小, 当应变达到-1%时形成能开始减小. 对于2B与2A其形成能在由拉应变到压应变的范围内逐渐增大, 其变化幅值约为0.6 eV. 由图5可知, 1B与1A以及2B与2A相互间的形成能差值随压应变的增大而增大, 随 着拉应变的增大而减小, 能量差值反映了体系结构的不对称性, 压应变在增大其差值的同时, 增大了体系印记, 但是其影响很小.
由于应变和双空位对体系极化都有影响, 因此在形成能的基础上计算了应变下含有双空位体系的极化, 如图6所示. 由图6可知, 应变下缺陷结构的极化随应变的变化趋势与完整晶格随应变的变化趋势相同. 且顺极化方向的 VBa- VO双空位组合增大了体系极化, 逆极化方向 VBa- VO双空位组合减小了体系极化. 但是应变下1A、2A组合极化值较完整晶格极化值的偏移量减小了, 而1B、2B则具有增大的趋势. 如果把顺极化方向组合的极化看做+ Pr, 而把逆极化方向组合的极化看做
![]() | (2) |
利用式(2)得到应变与极化偏移的关系如图7所示, 在压应变作用下1A和2A所对应的极化偏移都随压应变的增大而增加, 随拉应变的增大而减小. 同时, 有序的 VBa- VO双空位引起印记效应的产生[ 23], 则由图7可知压应变增大了有序的 VBa- VO双空位引起的印记效应, 但是其影响很小. 1A结构在-3%的应变条件下增大了1μC/cm2, 而2A只有0.3μC/cm2, 相对于应变对体系极化的影响可以忽略.
![]() | 图6 含有 VBa- VO双空位体系极化随应变的变化Fig. 6 Variation of polarization with in-plane strain for systems cantaining VBa- VO divacancy |
采用第一性原理计算了应变对四方相钛酸钡缺陷和极化的影响, 应变改变了空位形成能, 这与钛酸钡中存在较强的共价杂化有关. 氧空位作为还原环境下的主要缺陷, 由于其形成能随压应变的增大而增大, 从而可以利用压应变来抑制氧空位的形成. 通过计算 VBa- VO双空位形成能, 我们发现虽然不同组合形成能随应变的变化趋势不同, 但是1A-1B以及2A-2B相互之间的能量差都随压应变的增大而增大, 使得体系结构的不对称性增大, 从而增大了体系的印记. 此外, 应变下含有双空位的体系极化在拉应变到压应变的范围内线性增大, 同时1A和2A结构对应的极化偏移也具有相同的变化趋势. 由于极化偏移是印记效应的主要表现, 所以压应变增大了有序钡氧双空位体系的印记效应, 但是其影响较小可以忽略.
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|
[13] |
|
[14] |
|
[15] |
|
[16] |
|
[17] |
|
[18] |
|
[19] |
|
[20] |
|
[21] |
|
[22] |
|
[23] |
|