求助：有限位移法计算三阶力常数文件产生的POSCAR文件过多

查看数: 2149 | 评论数: 9 | 收藏 Add to favorites 0

关灯 | 提示：支持键盘翻页<-左右->

发布时间: 2024-5-4 16:11

结合Phonopy和Phono3py算拉曼，需要使用有限位移法计算三阶力常数，我用的原胞有26个原子：H S Cd 1 13 12 ，但是我用Phono3py产生POSCAR--*的时候有两万多个POSCAR，这怕不是要算到服务器炸掉，我看其他人都是只有 ...

回复 Reply

ymeng 发表于 Post on 2024-5-13 11:35:39

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

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

936039846 发表于 Post on 2024-5-7 20:16:30

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

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

卡开发发 发表于 Post on 2024-5-7 13:23:55

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

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

ymeng 发表于 Post on 2024-5-7 12:26:47

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

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

ymeng 发表于 Post on 2024-5-7 12:25:42

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

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

936039846 发表于 Post on 2024-5-6 09:47:26

PHONO3PY临近原子考虑了全部，用shengbte选取合理的截断半径会少一点

卡开发发 发表于 Post on 2024-5-5 06:28:54

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

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

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

ymeng 发表于 Post on 2024-5-4 21:12:35

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

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

乐平 发表于 Post on 2024-5-4 16:13:09

体系的对称性太低了吧

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register