ljh123 发表于 2025-6-1 15:21
要不改用VASP吧。VASP确实会更稳健些
另外其实你对结构优化不收敛做的尝试是比较少的，还没到放弃的时候 ...

明白了，谢谢老师，我再尝试尝试

sobereva 发表于 2025-6-1 09:12
优先考虑北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/KFP）幻灯片里说的这两条：

谢谢老师

THY 发表于 2025-6-1 15:08
置信半径改成了0.05还是不行，真头疼啊

要不改用VASP吧。VASP确实会更稳健些
另外其实你对结构优化不收敛做的尝试是比较少的，还没到放弃的时候。你之前尝试的OT、对角化、smearing等等，都是针对SCF不收敛做的。结构优化不收敛就去尝试不同的步长、最优化方法。

而且，对于自己没做过的材料，还是表面体系。看你的cif文件也是从数据库下载过来直接切的。这种初始结构，再加上用的是CP2K，是我的话我会先算个100步至少，再看看趋势

ljh123 发表于 2025-5-31 22:33
首先我不是老师哈。我建议把置信半径减小，比如到0.05试试。对角化和OT之类的参数主要是针对SCF计算的。 ...

置信半径改成了0.05还是不行，真头疼啊

kaikai 发表于 2025-6-1 09:04
开OT计算一下试一试

试过了，也还是不行

sobereva 发表于 2025-6-1 09:12
优先考虑北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/KFP）幻灯片里说的这两条：

谢谢老师，我试一下

优先考虑北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/KFP）幻灯片里说的这两条：


另外，别就光看前面20步。起码跑好几十步若没收敛的苗头到时候再说

开OT计算一下试一试

THY 发表于 2025-5-31 20:08
老师，不是没上传全，是前20步是震荡趋势的话，后面可能就没必要算了，对于普通晶体表面结构的优化，一般 ...

首先我不是老师哈。我建议把置信半径减小，比如到0.05试试。对角化和OT之类的参数主要是针对SCF计算的。目前来看应该去调一下结构优化部分的参数。
并且建议增加下结构优化步数。初始结构不太好的话，20步还在震荡还是有可能的。

ljh123 发表于 2025-5-31 20:00
是文件没上传全吗。看起来结构优化只跑了20步啊，这看不出什么趋势。表面结构优化很慢的，跑个几百步看看。 ...

老师，不是没上传全，是前20步是震荡趋势的话，后面可能就没必要算了，对于普通晶体表面结构的优化，一般好像100步以内，肯定是优化完成的，而且能量也肯定是逐渐递减的

是文件没上传全吗。看起来结构优化只跑了20步啊，这看不出什么趋势。表面结构优化很慢的，跑个几百步看看。还有可以尝试降低TRUST_RADIUS

Uus/pMeC6H4-/キ 发表于 2025-5-31 17:12
Mg是典型的在DZVP-MOLOPT-SR-GTH基组/GTH-PBE赝势/GPW组合下对截断能要求极高的元素，现在CUTOFF才400 Ry、 ...

老师您好，根据您的建议，我将cutoff改为1200以后，体系的能量依然在震荡，我将文件更新在1楼中

Mg是典型的在DZVP-MOLOPT-SR-GTH基组/GTH-PBE赝势/GPW组合下对截断能要求极高的元素，现在CUTOFF才400 Ry、SCF收敛限才1.0E-6，这些原子的能量和受力都计算不够准，几何优化自然变形明显、难以收敛。按https://www.cp2k.org/faq:cutoff的标准，Mg需要高达1200 Ry的截断能。如果不想提高截断能那么多，得换用要求更低的GAPW结合全电子基组做计算。