计算化学公社

标题: cp2k进行振动分析，输入文件中使用不同NPROC_REP，得到的虚频数目不同 [打印本页]

作者
Author: 。。。 时间: 2024-10-21 22:26
标题: cp2k进行振动分析，输入文件中使用不同NPROC_REP，得到的虚频数目不同
各位老师好，我对优化后的Pt6团簇进行振动分析，使用的服务器是128核，只修改了NPROC_REP的数目但是得到的虚频数目却完全不同，我也尝试使用更严格的EPS_SCF和EPS_DEFAULT,但是虚频始终无法消除。我感觉如果只是修改NPROC_REP数目就会对频率的结果造成影响的话，应该是我的输入文件存在问题？但是目前我暂时没有头绪，希望各位老师指导。

作者
Author: sobereva 时间: 2024-10-22 00:43
把SCF收敛精度设得更严以减小随机性的影响。并且注意检查副本的out文件看SCF都收敛没有

作者
Author: 。。。 时间: 2024-10-22 14:03

sobereva 发表于 2024-10-22 00:43
把SCF收敛精度设得更严以减小随机性的影响。并且注意检查副本的out文件看SCF都收敛没有

老师，我把SCF精度设置地更严格，但是算出来频率还是存在多个虚频，我也检查了副本的out文件，发现有一两次的SCF没有收敛。我后续把输入文件中LSD删掉再进行了一次振动分析，得到的频率中就不存在虚频，并且out文件中的SCF均收敛，所以老师，我之前的计算中出现虚频是因为LSD的原因吗

作者
Author: sobereva 时间: 2024-10-22 15:24

。。。发表于 2024-10-22 14:03
老师，我把SCF精度设置地更严格，但是算出来频率还是存在多个虚频，我也检查了副本的out文件，发现有一两 ...

如果体系本来就是闭壳层，就不该开LSD（等价于UKS），否则额外浪费很多时间，往往SCF也更困难
但如果实际是开壳层，就必须开LSD，并且遇到问题解决什么问题，且确保体系收敛到了真正的基态（比如做波函数稳定性测试）
对于仅仅算孤立的Pt6，强烈不建议用CP2K，应当用Gaussian等主流量子化学程序，能提供波函数稳定性测试功能，而且DFT还有解析Hessian

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