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

技術要約：BiAlO3における有限温度特性の原子論的研究

問題提起
鉛フリーのペロブスカイト強誘電体は、持続可能な技術応用において極めて重要であるが、鉛ベースの材料に代わる実行可能な選択肢の探索は依然として困難である。ビスマスアルミネート（BiAlO3）は、高圧・高温技術を通じて合成される有望な候補として浮上している。実験的研究により、その菱面体晶（$R3c$）基底状態と強誘電性が確立されている一方で、実験データと理論予測の間には大きな乖離が存在する。具体的には、絶対零度（0 K）に限定される密度汎関数理論（DFT）計算では、自発分極の値が（$72.3–82.5 \, \mu\text{C/cm}^2$）と報告されており、これは実験による測定値（$9–29 \, \mu\text{C/cm}^2$）と大幅に異なっている。さらに、既存の計算手法はBiAlO3の有限温度における挙動をモデル化する能力を欠いており、自発分極、歪み、抗電界、および強誘電（FE）モードと反強電位（AFD）モードの結合の温度進化に関する理解に空白を残している。

手法
これらの限界に対処するため、著者らはBiAlO3を有限温度でシミュレートするための第一原理に基づく原子論的モデルを開発した。その手法は以下のステップを含む：

1. フォノン解析： VASPパッケージ内の密度汎関数摂動論（DFPT）を用い、著者らは立方晶（$Pm\bar{3}m$）相における構造的不安定性を特定した。これにより、$\Gamma$モード（強誘電歪み）および$R$点および$M$点（反強電位酸素八面体の傾斜）に対応する虚数振動数が明らかになった。
2. 有効ハミルトニアンの構築： 特定された不安定性に基づき、有効ハミルトニアンのパラメータ化を行った。このハミルトニアンには、強誘電局所モード（$u$）、反強電位モード（$\omega$、$R_{25}$歪みに対応）、および歪み（$\eta$）の項が含まれる。決定的なことに、これはモード間の結合（FE-歪み、AFD-歪み、およびFE-AFD）を考慮している。パラメータはDFT計算から導出され、原子単位で提供されている。
3. モンテカルロ・シミュレーション： パラメータ化されたハミルトニアンは、単一反転メトロポリス・モンテカルロ（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}$）。負の圧力はこれとは逆の効果を持つ。
  • 一軸応力： 引張一軸応力は、一連の2つの相転移（立方晶 $\to$ 正方晶 $\to$ 菱面体晶）を誘起するのに対し、圧縮応力は単一の転移（立方晶 $\to$ 菱面体晶）をもたらす。
  • 二軸応力： トレンドは逆転し、引張二軸応力は単一の転移を導くが、圧縮二軸応力は2つの転移を誘起する。

意義と主張
本論文は、BiAlO3の有限温度における挙動を拡張するための初の計算的アプローチを提供することで、文献における重要な空白を埋めるものであると主張している。パラメータ化された有効ハミルトニアンを開発することにより、本研究は、実験的に関連のある条件下（FE-AFD結合の役割や機械的荷重への応答を含む）での材料の挙動に関する原子論的な洞察を提供する。

著者らは、シミュレーションによって予測された高い飽和歪みおよび残留歪みが、BiAlO3が形状記憶用途への有望な候補となり得ることを示唆していると強調している。さらに、本研究は、静水圧および機械的応力を通じてBiAlO3のキュリー温度および相安定性を調整できることを示しており、特定のデバイス応用に向けて材料を最適化するための潜在的な経路を提示している。この研究は、鉛フリー強誘電体の機能的特性を完全に理解するためには、絶対零度のDFT計算を超えていく必要があることを強調している。