小木虫

当前位置：首页 > 第一原理 >phonopy 为什么有限位移法得到的二阶力常数不完整？

phonopy 为什么有限位移法得到的二阶力常数不完整？

作者 Sahira: 来源: 小木虫 250 5 举报帖子

第一个问题：
同时用有限位移法和DFPT法计算了声子谱，得到的力常数文件不一致，有限位移法得到的不完整。
共有96组，有限位移法只有12组，DFPT法96组都有

第二个问题：
为什么不管什么体系，二阶力常数都是96组？我算了三个体系，二阶力常数文件都是96组。

截屏2020-09-04 上午10.36.18.png

截屏2020-09-04 上午10.34.53.png 返回小木虫查看更多

今日热帖

精华评论

漫天飘雪

根据重复性，删掉了一些，如果你想得到完整的，加入一个参数。
FULL_FORCE_CONSTANTS
FULL_FORCE_CONSTANTS = .TRUE. is used to compute full supercell constants matrix. The default setting is .FALSE.. By .TRUE. or .FALSE., the array shape becomes (n_patom, n_satom, 3, 3) or (n_satom, n_satom, 3, 3), respectively. The detail is found at FORCE_CONSTANTS and force_constants.hdf5.
Sahira

引用回帖:
2楼: Originally posted by 漫天飘雪 at 2020-09-04 11:16:19
根据重复性，删掉了一些，如果你想得到完整的，加入一个参数。
FULL_FORCE_CONSTANTS
FULL_FORCE_CONSTANTS = .TRUE. is used to compute full supercell constants matrix. The default setting is .FALSE.. By ...

谢谢，但是从官网的手册来看，FORCE_CONSTANTS是从FORCE_SETS中产生的，用有限位移法得到的FOECE_SETS就不完整，所以提取的力常数也不完整。官网上写的是VASP还是得用DFPT方法得到完整的力常数。

截屏2020-09-04 下午9.17.26.png

，
漫天飘雪

它是完整的，只是由于重复性，少写了一部分。你可以通过phonopy的命令，把你认为的“完整”FORCE_CONSTANTS写出来。需要加入以下两个参数：
FULL_FORCE_CONSTANTS = .TRUE.
WRITE_FORCE_CONSTANTS = .TRUE.
Sahira

更新：漫天飘雪说的对，是可以的
Sahira

引用回帖:
4楼: Originally posted by 漫天飘雪 at 2020-09-04 21:32:23
它是完整的，只是由于重复性，少写了一部分。你可以通过phonopy的命令，把你认为的“完整”FORCE_CONSTANTS写出来。需要加入以下两个参数：
FULL_FORCE_CONSTANTS = .TRUE.
WRITE_FORCE_CONSTANTS = .TRUE.

加了这两个参数之后是产生了完整的力常数，非常感谢！