计算化学公社

标题: scf总是不收敛应该怎么办？ [打印本页]

作者
Author: 心羽1971411 时间: 2021-4-7 13:59
标题: scf总是不收敛应该怎么办？
本帖最后由心羽1971411 于 2021-4-7 14:00 编辑

老师，您好。我想要使用cp2k计算NaCl-KCl-K2TiF6熔盐体系的Raman光谱。在利用packmol生成pdb文件的时候，向盒子里填充了20Na 40k 40Cl 10Ti 60F；随后使用Multiwfn产生cp2k做结构优化的输入文件，但是在计算过程中总是出现scf不收敛的情况，请问老师是我哪里的设置出现了问题？计算过程中涉及到的文件列于下方。

作者
Author: sobereva 时间: 2021-4-8 03:26
这个初始结构的原子分布太诡异了，和实际相去甚远，SCF必然不好收敛
SCF不收敛先别管它，跑一阵子MD，等原子的分布逐渐自发调节相对合理一些了，到时候可能就容易收敛了
这种体系只用优化出的一个几何结构做计算是没意义的

作者
Author: 心羽1971411 时间: 2021-4-8 08:53

sobereva 发表于 2021-4-8 03:26
这个初始结构的原子分布太诡异了，和实际相去甚远，SCF必然不好收敛
SCF不收敛先别管它，跑一阵子MD，等原 ...

谢谢老师。
老师我还想问一下如果我利用优化后得到的构型然后按照cp2k官网上计算Raman的例子进行计算(How to calculate vibrational spectra from AIMD with CP2K and TRAVIS (Updated 2018).)，也就是按照Massive&Global Equilibration、position trajectory、electron density、Compute the Electromagnetic Properties、Compute the Spectra的顺序计算，这样计算Raman光谱应该问题不大了吧

作者
Author: 喵星大佬 时间: 2021-4-8 11:09

心羽1971411 发表于 2021-4-8 08:53
谢谢老师。
老师我还想问一下如果我利用优化后得到的构型然后按照cp2k官网上计算Raman的例子进 ...

这种东西也只能这么算，而且MD时间还不能太短

作者
Author: 心羽1971411 时间: 2021-4-8 11:29

喵星大佬发表于 2021-4-8 11:09
这种东西也只能这么算，而且MD时间还不能太短

好的，谢谢你。有了你的回复心里总算有点底了。

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