计算化学公社

标题: 求助：有限位移法计算三阶力常数文件产生的POSCAR文件过多 [打印本页]

作者
Author: ymeng 时间: 2024-5-4 16:11
标题: 求助：有限位移法计算三阶力常数文件产生的POSCAR文件过多
结合Phonopy和Phono3py算拉曼，需要使用有限位移法计算三阶力常数，我用的原胞有26个原子：H S Cd 1 13 12 ，但是我用Phono3py产生POSCAR--*的时候有两万多个POSCAR，这怕不是要算到服务器炸掉，我看其他人都是只有几百个，为什么我会产生几万个POSCAR?调整哪个参数可以减少POSCAR的数量？command: phono3py -d --dim-fc2="2 2 1" --dim="1 1 1" -c POSCAR-unit

作者
Author: 乐平 时间: 2024-5-4 16:13
体系的对称性太低了吧

对称性越高，POSCAR-* 的数目就越少

作者
Author: ymeng 时间: 2024-5-4 21:12

乐平发表于 2024-5-4 16:13
体系的对称性太低了吧

谢谢前辈，那提高对称性要在结构优化的阶段就得考虑了吧。。。要强行提高对称性吗？

作者
Author: 卡开发发 时间: 2024-5-5 06:28

ymeng 发表于 2024-5-4 21:12
谢谢前辈，那提高对称性要在结构优化的阶段就得考虑了吧。。。要强行提高对称性吗？

体系确有对称性那就提高一下，没有那也没办法。原理上说有对称性的情况程序能识别的等价原子多，这样可以减少一些扰动产生的结构。

另外有vasp_raman.py的脚本可以试试，这么推荐的原因是，这样允许使用DFPT来进行二阶导数计算，但具体谁能算的更少点我不好说，具体有什么坑我不太确定。

作者
Author: 936039846 时间: 2024-5-6 09:47
PHONO3PY临近原子考虑了全部，用shengbte选取合理的截断半径会少一点

作者
Author: ymeng 时间: 2024-5-7 12:25

卡开发发发表于 2024-5-5 06:28
体系确有对称性那就提高一下，没有那也没办法。原理上说有对称性的情况程序能识别的等价原子多，这样可以 ...

我已经用DFPT算完二阶力常数了，但是我看帖子说三阶力常数必须用有限位移法来计算，他们之间的关系我不是很懂

作者
Author: ymeng 时间: 2024-5-7 12:26

936039846 发表于 2024-5-6 09:47
PHONO3PY临近原子考虑了全部，用shengbte选取合理的截断半径会少一点

好的，看来必须要用shengbte了，phono3py每次都是两万多POSCAR

作者
Author: 卡开发发 时间: 2024-5-7 13:23

ymeng 发表于 2024-5-7 12:25
我已经用DFPT算完二阶力常数了，但是我看帖子说三阶力常数必须用有限位移法来计算，他们之间的关系我不是 ...

因为程序没解析三阶导数，三阶导需要用二阶导差分出来，当然理论上说也可以通过能量和一阶导数（即力）去构造，但不管怎样体系一大这样要计算的构型点就非常多。

作者
Author: 936039846 时间: 2024-5-7 20:16

ymeng 发表于 2024-5-7 12:26
好的，看来必须要用shengbte了，phono3py每次都是两万多POSCAR

也可以试一试alamode的lasso法（类似机器学习）计算三阶力常数

作者
Author: ymeng 时间: 2024-5-13 11:35

936039846 发表于 2024-5-6 09:47
PHONO3PY临近原子考虑了全部，用shengbte选取合理的截断半径会少一点

您好，我在安装shengbte的时候，修改arch.make文件的时候MPIFC=mpiifort这个部分一直报错，这个部分是做什么的？

欢迎光临计算化学公社 (http://bbs.keinsci.com/)