Relaxation approach to quantum-mechanical modeling of ferroelectric and antiferroelectric phase transitions

技術要約：強誘電および反強誘電相転移の量子力学的モデリングへの緩和アプローチ

問題提起
強誘電および反強誘電材料は、温度および電場によって誘起される特徴的なヒステリシスループを示す。これらの現象の伝統的なモデリングは、古典的な枠組みに依存しており、通常、転移を古典的な領域内で発生するアレニウス型の活性化プロセスとして扱う。これらのモデルは、メタステーブル（準安定）状態にある系が、熱活性化を介してエネルギー障壁を乗り越えなければならないと仮定しており、転移速度は $e^{-E_A/k_BT}$ に比例すると想定している。

しかし、著者らは、第一原理シミュレーションにこれらの伝統的な手法を適用した際の決定的な失敗を指摘している。具体的には、古典的なモデルは非物理的に高い保磁界およびスイッチング電場を予測する（例：PbTiO3に対して実験値の約175 kV/cmに対し、約3018 kV/cmと予測する）。また、実験的な0 Kでのヒステリシスループを再現することにも失敗する。この不一致は、古典的な遷移状態理論において、エネルギー障壁が消失するまで（寿命が無限大となるまで）系がメタステーブル状態に閉じ込められる一方で、実際の系は有限の電場でスイッチングを行うという点に起因する。さらに、標準的な密度汎関数理論（DFT）シミュレーションは、相空間を制限する小さなスーパーセルサイズ（例：ドメイン形成を妨げる）に限定されており、多くの場合、現実の試料における障壁を低下させる欠陥を無視して、純粋な材料を想定している。SrTiO3やBaTiO3のような材料においては量子効果（零点エネルギーなど）が役割を果たしていることが示唆されてきたが、量子力学的な取り扱いがヒステリシスループやスイッチング電場に関する特定の論争を解決できるかどうかは、未だ確立されていなかった。

手法
著者らは、根本的に異なる枠組みを提案している。すなわち、相転移を古典的な活性化イベントとしてではなく、量子力学的領域における緩和プロセスとしてモデリングすることである。彼らは、以下の3つの具体的なモデルを開発し、比較検討している。

1. 基底状態緩和（GSR）モデル:

   • 系は時間発展する量子状態 $|\psi(t)\rangle$ によって記述される。
   • 励起状態（メタステーブル状態）から基底状態への緩和をシミュレートするために、エネルギー固有値 $E_n$ に虚数成分を含める：$E'_n = E_n - i\hbar\gamma_n/2$。ここで $\gamma_n$ は緩和率である。
   • 緩和率は、すべての低エネルギー状態に関連するボーア周波数の和に比例すると定義される：$\gamma_n = k_{scale} \sum_{E_m < E_n} (E_n - E_m)/\hbar$。
   • 分極は、分極演算子の期待値として計算される。このモデルは0 Kに限定される。

2. 密度行列ベース（DMB）モデル:

   • 有限温度を組み込むために、著者らは密度行列形式 $\rho(t)$ を利用する。
   • 時間発展は、単なる基底状態ではなく、平衡密度行列 $\rho_{eq}$（カノニカルアンサンブル）へと系を駆動する修正されたリウヴィル方程式に従う。
   • この方程式は、ボルツマン因子 $e^{-(E_m - E_n)/k_BT}$ によって正規化された、ポピュレーションの減衰と獲得の両方を考慮している。
   • 物理的な妥当性を確保するために、トレース保存項が追加されている。

3. リンドブラッド演算子（LO）モデル:

   • このアプローチは、熱浴との相互作用を記述するために、リンドブラッド演算子 $L_{kj}$ を用いたマスター方程式を用いて緩和を実装する。
   • この手法は、密度行列のトレースを明示的に保存し、詳細釣合いの条件を満たす。これにより、状態間の転移速度の比が平衡時におけるボルツマン分布と一致することを保証する。

主要な入力およびパラメータ

• エネルギーランドスケープ: PbTiO3、PbZrO3、PbHfO3、CsGeBr3、HfO2、およびAlNを含む材料について、DFT（VASP）を用いて計算される。エネルギープロファイル $U(P)$ は、均一な分極反転（ナノドメインの核生成）を表す、各相間の線形補間から導出される。
• 質量と電荷: 極性モードの有効質量は、DFT計算による調和振動数から導出される。ボルン有効電荷（$Z^*$）は、分極をスーパーセルの歪みに変換する。
• 緩和率: 解析的な表現（フェルミの黄金律）が存在するが、緩和係数 $k_{scale}$ を通じて実験データ（IR/Raman線幅）にフィッティングされる。
• ダイナミクス: 主なシミュレーションでは、固有のダイナミクス（ボーア周波数）は緩和タイムスケール上で平均化されるため、その影響はほぼ無視できるとしてオフにされている。

結果

• 反強誘電体（PbZrO3）: GSRモデルは反強誘電ヒステリシスループを正常に再現し、0 Kにおいて実験値に近いスイッチング電場（315 kV/cm）を予測した。これは古典的モデルでは達成不可能な結果である。DMBおよびLOモデルは、これを有限温度（10–1000 K）へと拡張する。
  • DMBモデルは、低温度において段階的な分極進化と矩形ヒステリシスループを予測する。
  • LOモデルは、指数関数的な進化と傾斜のあるヒステリシスループを予測し、PbZrO3のループの具体的な形状に対するより良い一致を示す。
• 強誘電体（PbTiO3）: これらのモデルは強誘電ヒステリシスループを再現し、非ペロブスカイト（AlN、HfO2）がペロブスカイトよりも著しく高い保磁界を示す傾向を正しく予測する。
• 温度依存性: モデルは、温度上昇に伴うヒステリシスループの細り（スリミング）を捉え、反強誘電から常誘電挙動への遷移を記述する。
• 共存電場の補正: 著者らは、DFTが共存電場（極性と反極性相のエネルギーが等しくなる電場）を過小評価していることを指摘している。シミュレーションのループ幅を実験と一致させるために、反強誘電体におけるシミュレーションで使用されるAC電場周波数を調整した。

意義および主張
本論文は、（反）強誘電相転移を量子力学的な緩和プロセスとして扱うことが、従来の古典的なアレニウス型モデルの限界を克服することを主張している。主な意義は、効率的な第一原理ベースのシミュレーションを実行することで、古典的な熱活性化が不可能な0 Kにおいても、実験と良好に一致するスイッチing電場とヒステリシス形状を得られる点にある。

著者らは、この枠組みが、長らく純粋に古典的であると考えられてきた転移における量子力学の本質的な役割を実証していると断言している。彼らは、このアプローチが広範な相転移（磁性、弾性、マルチフェロイック、電子的転移）や、量子トンネル効果および化学反応速度を伴う他の現象にも適用可能であることを示唆している。本研究は、相転移を理解するための新しい視点を提供するものであり、特に「型破りな」仮定である緩和駆動型の量子進化が、強誘電材料のモデリングにおける不一致を解決するために必要であることを強調している。