pymatgenで状態図を描く

pymatgenで状態図を描く方法の備忘録。

ついでにMaterials Project Rest を使う方法も調べた。

Materials Project API key の登録

まず、Materials Project Rest を使うにはユーザー登録してAPIキーを作成しなければならない。

登録は以下のリンクから行える: https://www.materialsproject.org/dashboard

コマンドラインで以下を実行すると.pmgrc.yamlに自分のAPIキーを設定できて、後でいちいちAPIキーをベタ書きしなくてよくなる

pmg config -a PMG_MAPI_KEY your_api_key

状態図を描く

pymatgen.analysis.phase_diagramを使うと、基準となる物質(単体でも化合物でもよい)の組成に関してギブス自由エネルギー密度のconvex hullを描くことができる(すごい)。 凸包の計算自体はscipy.spatial.ConvexHullでやっているみたい。

gist.github.com

参照リンク