?經(jīng)驗分享

格林艾森常數(shù)（Grüneisen parameters）為表示聲子非簡諧效應的參數(shù)。phonopy計算模式格林艾森參數(shù)方式如下。

http://phonopy.github.io/phonopy/gruneisen.html#phonopy-gruneisen

以NaCl為例，首先需要進行三個聲子譜計算，分別是平衡體積（orig），稍小體積（minus）和稍大體積（plus）的結(jié)構(gòu)，計算時改變晶格參數(shù)使得體積獲得改變。

計算完成后獲得不同體積下的原子間力常數(shù)文件FORCE_SETS 或 FORCE_CONSTANTS ??（后續(xù)命令中需要添加--readfc參數(shù)）。

隨后通過phonopy-gruneisen（舊版gruneisen）命令獲得對應于聲子帶結(jié)構(gòu)中倒易空間中的路徑的模式格林艾森常數(shù)（具體設定可根據(jù)自身計算情況更改）

phonopy-gruneisen orig plus minus --dim="2 2 2" --pa="0 1/2 1/2 1/2 0 1/2 1/2 1/2 0" --band="1/2 1/4 3/4 0 0 0 1/2 1/2 1/2 1/2 0.0 1/2" -p -c POSCAR-unitcell

編輯

在倒易網(wǎng)格中的模式格林艾森常數(shù)（具體設定可根據(jù)自身計算情況更改）：

phonopy-gruneisen orig plus minus --dim="2 2 2" --pa="0 1/2 1/2 1/2 0 1/2 1/2 1/2 0" --mesh="20 20 20" -p -c POSCAR-unitcell --color="RB"

編輯

phonopy-gruneisen 命令選項包括--dim，--mp, --mesh，--band，--pa, --primitive_axis，--readfc，--band_points，--nac，--factor，--nomeshsym，-p，-c，-s, --save，-o等，功能上均與phonopy命令相同，但除上述列出的選項外，可以在phonopy命令上使用的其它選項是不能在phonopy-gruneisen 上使用的。如有必要，可以自行修改phonopy-gruneisen代碼。

如查看幫助，可輸入

phonopy-gruneisen -h

通過 phonopy-gruneisenplot 命令導出計算得到的數(shù)據(jù)，使用方式類似于phonopy-bandplot。

愿有所成

IEchoQ

引喻失義? ?妄自菲薄

IEchoQ