这里会显示出您选择的修订版和当前版本之间的差别。
两侧同时换到之前的修订记录前一修订版后一修订版 | 前一修订版上一修订版两侧同时换到之后的修订记录 | ||
adf:phononandthermal [2017/05/08 10:43] – liu.jun | adf:phononandthermal [2024/05/17 21:57] – [如何计算声子谱、热力学性质] liu.jun | ||
---|---|---|---|
行 1: | 行 1: | ||
======如何计算声子谱、热力学性质====== | ======如何计算声子谱、热力学性质====== | ||
- | BAND使用高精度数值基、Slater基组,因此精度非常高,但效率较低。声子谱的计算量较大,因此首先推荐使用ATK软件进行计算,参考:[[atk: | + | BAND使用高精度数值基、Slater基组,因此精度非常高,但效率较低。声子谱的计算量较大,因此首先推荐使用QATK、Quantum ESPRESSO软件进行计算,参考:[[atk: |
+ | 2024版以后,有BANDlicense,就可以使用AMS中完全集成化的Quantum Espresso。建议用户直接使用Quantum Espresso计算即可:[[adf: | ||
=====参数设置===== | =====参数设置===== | ||
- | BAND的泛函很丰富,不过这里我们使用PBE泛函。BAND的基组比ADF同级别基组精度更高(因为不仅使用了同级Slater基组,还插入了数值基),因此本例使用DZP基组就足够了。冻心不影响振动性质,因此使用Frozen core: Large。BAND会自动根据用户精度需求设置对应的k点密度,此处用了Numerical quality: Normal。并自动根据晶体的对称性,设置计算声子谱的k空间路径。参数设置如下图所示: | + | BAND的泛函很丰富,不过这里我们使用PBE泛函。BAND的基组比ADF同级别基组精度更高(因为不仅使用了同级Slater基组,还插入了数值基),因此本例使用DZP基组就足够了。冻心不影响振动性质,因此使用Frozen core: Large。参数设置如下图所示: |
{{ : | {{ : | ||
+ | {{ : | ||
+ | |||
+ | 点击Task: | ||
+ | |||
+ | {{ : | ||
=====结果查看===== | =====结果查看===== | ||
- | ADF LOGO > Spectra打开谱图窗口,点击菜单栏Spectra > Phonon Spectrum,显示声子谱: | + | ADF LOGO > BandStructure,显示声子谱: |
{{ : | {{ : | ||
- | 图中有三个窗口:左边的窗口显示的是k空间,及其高对称点在k空间的位置;中间窗口显示的是沿着这些高对称点的声子谱函数,下方path显示的是声子谱的k点路径(注意字母G代表高对称点Γ);右边显示点是声子态密度。 | + | 图中有三个窗口:左边的窗口显示的是k空间,及其高对称点在k空间的位置;中间窗口显示的是沿着这些高对称点的声子谱函数,下方path显示的是声子谱的k点路径(注意字母G代表高对称点Γ),也就是声子谱的横坐标,与左边窗口可以对照着看;右边显示点是声子态密度。 |
ADF LOGO > Output打开输出文件的文本内容,拖动窗口滚轮到如下位置: | ADF LOGO > Output打开输出文件的文本内容,拖动窗口滚轮到如下位置: | ||
行 47: | 行 53: | ||
</ | </ | ||
显示从0K到1000K的热力学性质。 | 显示从0K到1000K的热力学性质。 | ||
- | |||
- | 用户也可以通过修改*.run文件中的关键词,修改温度范围: | ||
- | <code bash> | ||
- | PhononConfig | ||
- | SuperCell | ||
- | 2 0 0 | ||
- | 0 2 0 | ||
- | 0 0 2 | ||
- | SubEnd | ||
- | MinTempKelvin 298 | ||
- | MaxTempKelvin 398 | ||
- | NumTemperatures 101 | ||
- | End | ||
- | </ | ||
- | 其中, | ||
- | <code bash> | ||
- | MinTempKelvin 298 | ||
- | MaxTempKelvin 398 | ||
- | NumTemperatures 101 | ||
- | < | ||
- | 设定温度范围为298K-398K,每1K的热力学都要计算,因此要计算101种温度。 | ||
- | |||
- | 其他更详细的设置,包括超胞、步长、对称性等,参考:https:// | ||
- |