ASEとQuantum Espressoを使って、六方晶Hfの格子定数aとcを最適化する - Qiita

はじめに Jupyter上で、PythonライブラリであるASE(Atomic Simulation Environment)を介して、第一原理計算ソフトQuantum Espressoで、第一原理計算を実施してみます。 本記事は、書籍『密度半関数理論入門(D.S.ショール，J.A.ステッケル 著、佐々木泰造，末原茂 共訳』(以下、教科書)の例題や課題を参考に、プログラムを作って計算します。 なお、固体DFTの専門家ではありませんので、計算条件、結果等はご自身でご確認/判断ください。 今回は、教科書3章,練習6で紹介されている、 『六方晶Hfの格子定数aとcを最適化する』 問題に取り組んでみます。 前準備 計算に必要になるライブラリや、擬ポテンシャルを用意します。 import os import subprocess import numpy as np import matplotli