Atomistic study of finite temperature properties in ferroelectric BiAlO$_3$

技术摘要：BiAlO3 在有限温度下的原子级性质研究

问题陈述
无铅钙钛矿铁电材料对于可持续技术应用至关重要，然而寻找可行的替代材料仍具挑战性。铋铝酸盐（BiAlO3）已成为一种极具前景的候选材料，通过高压高温技术合成。虽然实验研究已经确立了其菱方（$R3c$）基态和铁电特性，但实验数据与理论预测之间存在显著差异。具体而言，密度泛函理论（DFT）计算（受限于零开尔文）报告的自发极化值（$72.3–82.5 \, \mu\text{C/cm}^2$）与实验测量值（$9–29 \, \mu\text{C/cm}^2$）存在大幅度偏差。此外，现有的计算方法缺乏模拟 BiAlO3 有限温度行为的能力，导致在理解 BiAlO3 自发极化、应变、矫顽场以及铁电（FE）模式与反铁畸变（AFD）模式之间耦合随温度演化过程方面存在空白。

方法论
为了解决这些局限性，作者开发了一种基于第一性原理的原子级模型，用于模拟有限温度下的 BiAl 3。该方法包括以下步骤：

1. 声子分析： 利用 VASP 软件包中的密度泛函微扰理论（DFPT），作者识别了立方相（$Pm\bar{3}m$）中的结构不稳定性。这揭示了对应于 $\Gamma$ 模式（铁电畸变）以及 $R$ 点和 $M$ 点（反铁畸变氧八面体倾斜）的虚频。
2. 有效哈密顿量构建： 基于识别出的不稳定性，构建并参数化了一个有效哈密顿量。该哈密顿量包含铁电局部模式（$u$）、反铁畸变模式（$\omega$，对应于 $R_{25}$ 畸变）和应变（$\eta$）项。至关重要的是，它考虑了这些模式之间的耦合（FE-应变、AFD-应变以及 FE-AFD）。参数由 DFT 计算得出，并以原子单位提供。
3. 蒙特卡洛模拟： 将参数化的哈密顿量应用于单翻转 Metropolis 蒙特卡洛（MC）模拟。采用了具有周期性边界条件的 $20 \times 20 \times 20$ 超胞。模拟采用模拟退火技术，从 1500 K 冷却至 10 K。为了准确捕捉相变附近的竞争不稳定性，在居里温度（$T_C$）附近增加了 MC 步数。
4. 外部扰动： 该模型被扩展用于研究外部电场（沿 [111]、[110] 和 [100] 方向）、静水压力以及单轴/双轴应力的影响。

关键结果

• 基态与相变： 模拟证实了从顺电立方相到铁电菱方（$R3c$）相的结构相变。预测的居里温度为 $1160 \, \text{K}$，且在低温下的自发极化沿 [111] 方向为 $81 \, \mu\text{C/cm}^2$。作者指出，虽然 $T_C$ 高于实验值（$>793 \, \text{K}$），这可能是由于参数化过程中高估了结合能，但极化值与 DFT 基态值一致。
• 模式耦合： 在较宽的温度和电场范围内观察到铁电模式与反铁畸变模式之间的强耦合。AFD 序参数模仿了 FE 模式的趋势，两者均沿 [111] 伪立方方向取向。
• 电场响应： 在沿 [111] 方向的电场作用下，材料表现出单一的回滞回线，证实了 $R3c$ 相。饱和极化（$P_{Sat}$）随电场方向略有变化（$57.8–60.2 \, \mu\text{C/cm}^2$），而矫顽场（$E_C$）表现出显著的各向异性（$2.4–3.8 \, \text{MV/cm}$）。AFD 和应变序参数表现出“蝴蝶型”回线行为，其中应变表现出高饱和度（$0.021$）和显著的剩余应变（$0.018–0.019$）。
• 温度依赖性： 随着温度升高，回滞回线面积、饱和极化、剩余极化和矫顽场均下降。这一趋势与其他铁电材料一致。
• 压力与应力效应：
  • 静水压力： 正静水压力会抑制自发极化和居里温度（变化率为 $-45 \, \text{K/GPa}$），而负压力则产生相反的效果。
  • 单轴应力： 拉伸单轴应力诱导一系列两次相变（立方 $\to$ 四方 $\to$ 菱方），而压缩应力导致单次相变（立方 $\to$ 菱方）。
  • 双轴应力： 趋势相反；拉伸双轴应力导致单次相变，而压缩双轴应力诱导两次相变。

意义与主张
本文声称填补了文献中的关键空白，提供了将 DFT 模拟扩展到 BiAlO3 有限温度的第一种计算方法。通过开发参数化的有效哈密顿量，本研究提供了在实验相关条件下（包括 FE-AFD 耦合的作用以及对机械载荷的响应）理解该材料行为的原子级见解。

作者强调，模拟预测的高饱和应变和剩余应变表明 BiAlO3 是形状记忆应用的有力候选材料。此外，研究表明可以通过静水压力和机械应力调节 BiAlO3 的居里温度和相稳定性，为针对特定器件应用优化材料提供了潜在途径。这项工作强调了超越零开尔文 DFT 计算以充分理解无铅铁电体功能性质的必要性。