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

ymeng

结合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

乐平

体系的对称性太低了吧


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

ymeng

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

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

卡开发发

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

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

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

936039846

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

ymeng

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

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

ymeng

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

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

卡开发发

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

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

936039846

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

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

ymeng

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

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