[論文レビュー] Ab Initio bulk free energy surface of proper ferroelectrics

Pinchen Xie, Yixiao Chen+4 more|ArXiv.org|May 24, 2022

Machine Learning in Materials Science被引用数 6

ひとこと要約

本研究では、密度汎関数理論（DFT）データで訓練されたディープポテンシャル（DP）機械学習ポテンシャルを用いて、PbTiO₃のab initioマルチスケールモデルを提案し、グローバルな分極の関数として自由エネルギー表面を正確に捉える。この手法は非調和的効果および非線形分極応答を含み、実験的エンタルピー、自発分極、比熱、誘電感受率と優れた一致を示す。一方、強化されたサンプリングにより、Curie温度を超えても持続的な局所双極子を有する一次遷移、順序-無秩序転移が明らかになった。

ABSTRACT

We report a systematic and accurate approach for deriving the bulk free energy surface (FES), a function of temperature, polarization, and strain, from the first-principles density functional theory (DFT) of proper ferroelectrics. The core of our approach is the metadynamics algorithm that extracts the polarization dependence of the FES from all-atom molecular dynamics simulations without an a priori ansatz. The rest of the FES is derived from the metadynamics trajectories that span the relevant phase space. We demonstrate our approach in the case of lead titanate. The errors across the phase transition, due to DFT numerics, all-atom molecular dynamics, and free energy evaluation by enhanced sampling, can be systematically controlled and are of the order of 1meV/atom. The accuracy of the resulting ab initio FES is only limited by the adopted functional approximation of DFT.

研究の動機と目的

• 線形近似や調和限界を越えて、フェロエレクトリック材料の全自由エネルギー表面を捕捉する第一原理に基づくマルチスケールモデルの構築。
• DFTデータで訓練された機械学習ポテンシャルを用いて、PbTiO₃における非線形的および非調和的寄与を含む分極応答を正確に記述すること。
• 強化されたサンプリングを用いて、PbTiO₃におけるフェロエレクトリック転移の性質を調査し、特にCurie温度を超えて局所双極子が持続するかを検討すること。
• 自発分極、比熱、誘電感受率を含む実験的熱力学的および誘電的性質と、モデルを比較して妥当性を検証すること。
• シミュレートされた自由エネルギー表面に基づいて、単一ドメインPbTiO₃結晶のランダウ＝デヴォンシアー型理論を確立すること。

提案手法

• PbTiO₃のポテンシャルエネルギー表面および分極表面を、第一原理DFTデータで訓練したディープポテンシャル（DP）モデルを用いて高精度に表現すること。
• グローバル分極の集団変数空間における自由エネルギー表面を探索するため、バイアス要因を用いたウェルテンペラードメタダイナミクスを用いる。
• 等温等圧条件の下で位相空間をサンプリングするため、時間刻み0.5 fsおよび周期的境界条件を用いたNPT-MDシミュレーションを実施する。
• 短い第一原理MD軌道から訓練およびテストデータセットを生成し、高精度なモデル訓練に十分な1835配置の低次元データセットを構築する。
• グローバル分極を、ベリー位相およびワニエ中心位置に基づく現代的分極理論により定義し、局所双極子を補助変数として用いる。
• 標準誤差1 eÅ（3×3×3スーパーセルで）の誤差指標を用いて、DPモデルのDFTデータに対する妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 DFTデータで訓練されたディープポテンシャルモデルは、グローバル分極の関数としてPbTiO₃の自由エネルギー表面をどの程度正確に再現できるか？
• RQ2 非調和的および非線形分極効果を完全に含めた場合、PbTiO₃におけるフェロエレクトリック転移の性質はいかなるものか？
• RQ3 Curie温度を超えて局所双極子が持続するか。これは、転移が順序-無秩序的性質を示す証拠となるか？
• RQ4 シミュレートされた自由エネルギー表面は、単一ドメイン結晶のランダウ＝デヴォンシアー型記述をどの程度支持するか？
• RQ5 モデルから計算された熱力学的性質—自発分極、比熱、誘電感受率—は、実験測定値とどの程度一致するか？

主な発見

• 3×3×3スーパーセルにおけるグローバル双極子予測で、標準誤差約1 eÅを達成し、極めて歪んだ構造に対しても高い精度を示した。
• 静的ボーンチャージを用いた線形モデルは、DPモデルの2倍の誤差を示し、相境界付近で顕著な外れ値を示した。これは、非線形的・非摂動的分極モデル化の必要性を確認した。
• ウェルテンペラードメタダイナミクスによる強化されたサンプリングにより、2ナノ秒後に反強誘電相と強誘電相の自由エネルギー差が収束し、確率的誤差は1ユニットセルあたり約0.01 meVであった。
• 自由エネルギー表面には、反強誘電相と強誘電相に対応する2つの明確な極小値が存在し、その障壁の高さと形状から一次転移が支持された。
• Curie温度を超えても局所双極子が持続することが観察され、これは転移が純粋に変位的モデルとは対照的に、順序-無秩序的性質を示す証拠となった。
• シミュレートされたエンタルピー、自発分極、比熱、誘電感受率は、単結晶の実験的測定値と優れた定量的一致を示した。