[論文レビュー] The generalized quasiharmonic approximation via space group irreducible derivatives

Mark A. Mathis, Amey Khanolkar+10 more|arXiv (Cornell University)|Feb 28, 2022

Acoustic Wave Phenomena Research参考文献 93被引用数 18

ひとこと要約

本論文は、空間群の既約微分を用いて一般化量子調和近似（QHA）を定式化し、結晶内の温度および応力依存の熱力学的性質を正確かつ効率的に計算可能にする。既約2次位移微分と2つのひずみパラメータ化法（グリッド補間とテイラー展開）を活用することで、トリウリアニド（ThO2）における弾性定数および熱膨張の高精度・再現性のある計算が可能となり、第一原理計算および中性子散乱データと照合された。

ABSTRACT

The quasiharmonic approximation (QHA) is the simplest nontrivial approximation for interacting phonons under constant pressure, bringing the effects of anharmonicity into temperature dependent observables. Nonetheless, the QHA is often implemented with additional approximations due to the complexity of computing phonons under arbitrary strains, and the generalized QHA, which employs constant stress boundary conditions, has not been completely developed. Here we formulate the generalized QHA, providing a practical algorithm for computing the strain state and other observables as a function of temperature and true stress. We circumvent the complexity of computing phonons under arbitrary strains by employing irreducible second order displacement derivatives of the Born-Oppenheimer potential and their strain dependence, which are efficiently and precisely computed using the lone irreducible derivative approach. We formulate two complementary strain parametrizations: a discretized strain grid interpolation and a Taylor series expansion in symmetrized strain. We illustrate our approach by evaluating the temperature and pressure dependence of the elastic constant tensor and the thermal expansion in thoria (ThO$_2$) using density functional theory with three exchange-correlation functionals. The QHA results are compared to our measurements of the elastic constant tensor using time domain Brillouin scattering and inelastic neutron scattering. Our irreducible derivative approach simplifies the implementation of the generalized QHA, which will facilitate reproducible, data driven applications.

研究の動機と目的

• 定圧ではなく定応力下で動作する一般化量子調和近似（QHA）を開発すること。
• 任意のひずみ下でのフォノン評価の計算複雑性を克服するため、空間群の既約微分を用いること。
• ひずみ依存の弾性エネルギーおよびダイナミカル行列を計算する実用的で数値的に安定なアルゴリズムを提供すること。
• 材料における温度および応力依存の弾性定数および熱膨張の正確な予測を可能にすること。
• トリウリアニド（ThO2）において、密度汎関数理論および中性子散乱を用いて実験データと照合すること。

提案手法

• 計算コストの低減と対称性の一貫性向上を図るため、空間群の既約微分を用いて一般化QHAを定式化する。
• 効率的なひずみ依存2次核移動微分の計算のため、唯一の既約微分アプローチを採用する。
• 2つのひずみパラメータ化法を適用する：(1) 離散的ひずみグリッド補間と(2) 対称化ひずみ変数におけるテイラー級数展開。
• 第一原理密度汎関数理論を用い、複数の交換相関関数形を用いてひずみの関数としてのダイナミカル行列および弾性エネルギーを計算する。
• パラメータ化されたひずみ依存性を用いて、ヘルムホルツ自由エネルギーからギブス自由エネルギーへのルジャンドル変換を実行する。
• 時間領域ブリューヌール散乱および非弾性中性子散乱測定によるThO2の弾性定数の実験データと照合して結果を検証する。

実験結果

リサーチクエスチョン

• RQ1 対称性に適合した微分を用いて、定応力下での一般化QHAを体系的に定式化する方法は何か？
• RQ2 ひずみ依存フォノンをパラメータ化する際、グリッド補間とテイラー級数展開の相対的利点および数値的安定性は何か？
• RQ3 既約微分の使用が、QHA計算における計算コストを顕著に低減させつつも、精度を保持できるか？
• RQ4 一般化QHAは、ThO2における弾性定数および熱膨張の温度および圧力依存性をどの程度正確に予測できるか？
• RQ5 時間領域ブリューヌール散乱および非弾性中性子散乱を用いた実験測定と比較して、QHAによる弾性定数の予測はどの程度良好か？

主な発見

• 既約微分アプローチにより、ひずみ依存ダイナミカル行列および弾性エネルギーの効率的かつ高精度な計算が可能となり、一般化QHAの実装が簡素化される。
• テイラー級数およびグリッド補間の両パラメータ化法ともに高精度な結果をもたらし、テイラー級数はひずみの所定の次数まで解析的制御を可能にする。
• 本手法は、ThO2における弾性定数テンソルの温度および圧力依存性を正確に予測でき、時間領域ブリューヌール散乱および非弾性中性子散乱による実験測定と良好に一致する。
• 複数の交換相関関数形を用いた計算において、ThO2の熱膨張係数が一貫しており、本手法の堅牢性を示している。
• 対称化ひずみ変数および既約成分の使用により、パラメータ化が結晶の点群対称性を尊重し、不適切なモード混合を回避する。
• 一般化QHAフレームワークにより、実験的に関連のある定応力条件下での再現性があり、データ駆動型の熱力学的性質予測が可能になる。