[論文レビュー] Thermal properties of materials from ab-initio quasi-harmonic phonons

Stefano Baroni, Paolo Giannozzi+1 more|arXiv (Cornell University)|Dec 21, 2011

High-pressure geophysics and materials参考文献 135被引用数 158

ひとこと要約

本論文は、密度汎関数理論（DFT）と準調和近似（QHA）を組み合わせた第一原理的手法を提示し、広い圧力および温度範囲における材料の熱的性質（比熱、エントロピー、フォノン寄与など）を正確に予測する。この手法はゼロ点振動効果を的確に捉え、炭酸塩、酸化物、高圧相において実験的傾向を再現し、地球物理学的条件下で経験的ポテンシャルを上回る性能を示す。

ABSTRACT

This paper gives a short overview of the calculation of thermal properties of materials from first principles, using the Quasi-Harmonic Approximation (QHA). We first introduce some of the thermal properties of interest and describe how they can be calculated in the framework of the QHA; then we briefly recall Density-Functional Perturbation Theory as a tool for calculating phonons from first principles, and present some codes that implement it; finally we review recent applications of first-principle QHA.

研究の動機と目的

• 実験が不可能な極限圧力および極限温度条件下における材料の熱的性質を予測するフレームワークを開発すること。
• 半経験的原子間ポテンシャルの限界を克服し、フォノン寄与を準調和近似を介して第一原理的DFTを用いること。
• 特に融解付近および高圧下における固体の量子的ゼロ点振動および非調和効果を正確に扱うこと。
• 地球内部に存在するケイ酸塩、酸化物、ペロブスカイトなどの地球物理学的に関連する材料をモデル化するための信頼性の高い計算ツールを提供すること。
• グリッドコンピューティングなどのスケーラブルな計算インfraストラクチャを活用した大規模な熱的および構造的挙動のシミュレーションを可能にすること。

提案手法

• 体積または圧力に依存するフォノン周波数を許容することで、温度依存の格子力学を準調和近似（QHA）を用いてモデル化する。
• 密度汎関数摂動論法（DFPT）を用いて、第一原理的電子構造計算から調和フォノン周波数および固有ベクトルを計算する。
• フォノンモードの統合を介して体積の関数としてヘルムホルツ自由エネルギーを計算し、ゼロ点エネルギー補正を含む。
• 標準的な統計力学的関係式を用いて、自由エネルギーから比熱、エントロピー、熱膨張率などの熱力学的性質を導出する。
• フォノンスペクトルおよび熱力学的応答関数の温度および圧力依存性を計算するために、この手法を適用する。
• 大規模並列およびグリッドコンputingアーキテクチャ上でスケーラブルなDFT+DFPT計算を実行するために、Quantum ESPRESSOソフトウェアスイートを用いる。

実験結果

リサーチクエスチョン

• RQ1 準調和近似と第一原理的DFTを組み合わせた手法が、地球物理学的条件下の圧力および温度範囲における材料の熱的性質を正確に予測できるか。
• RQ2 酸化アルミニウム（corundum）やジルコニア（stishovite）のような高圧相の安定性および熱力学的性質に、ゼロ点振動効果がどの程度寄与するか。
• RQ3 カルサイト（calcite）やアラゴナイト（aragonite）のような酸化物において、経験的原子間ポテンシャルが第一原理的手法と比較して、熱的挙動を正しく捉えられていない程度はどの程度か。
• RQ4 低温におけるα-Al₂O₃相の安定化にフォノンが果たす役割は何か。また、高圧多形と比較してその寄与はいかなるものか。
• RQ5 QHA-DFT手法が、SiO₂およびMgSiO₃多形における実験的比熱および熱膨張係数をどの程度正確に再現できるか。

主な発見

• QHA-DFT手法は、α-Quartzおよびstishovite SiO₂の実験的比熱を広い温度範囲で正確に再現し、定量的な一致が得られた。
• ゼロ点振動エネルギーはstishoviteの熱力学的性質に顕著な影響を及ぼし、特に低温域では熱膨張効果を上回る。
• QHAは、低温におけるα-Al₂O₃相がゼロ点振動によって安定化されることを予測しているが、フォノンのみでは室温で安定化されない。
• カルサイト（CaCO₃）においては、高温域ではゼロ点振動が平衡体積に及ぼす寄与が熱膨張よりも顕著であり、その重要性が示された。
• この手法により、α-Al₂O₃相が高圧下でCaIrO₃型およびU₂S₃型多形に転移することが判明し、実験的相図と整合的であった。
• QHA-DFT手法は、Exo-惑星内部に相当する超高圧下におけるSiO₂のピリット型からコタンサイト型への相転移を特定したが、対応する地震的不連続性は観測されなかった。