phonopy计算模式格林艾森常数

​经验分享

格林艾森常数（Grüneisen parameters）为表示声子非简谐效应的参数。phonopy计算模式格林艾森参数方式如下。

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

以NaCl为例，首先需要进行三个声子谱计算，分别是平衡体积（orig），稍小体积（minus）和稍大体积（plus）的结构，计算时改变晶格参数使得体积获得改变。

计算完成后获得不同体积下的原子间力常数文件FORCE_SETS 或 FORCE_CONSTANTS   （后续命令中需要添加--readfc参数）。

随后通过phonopy-gruneisen（旧版gruneisen）命令获得对应于声子带结构中倒易空间中的路径的模式格林艾森常数（具体设定可根据自身计算情况更改）

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

在倒易网格中的模式格林艾森常数（具体设定可根据自身计算情况更改）：

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 命令导出计算得到的数据，使用方式类似于phonopy-bandplot。

愿有所成

引喻失义   妄自菲薄