谢谢老师，我重新画一下初始结构试试

龚湘姣 发表于 2022-8-7 16:14
老师，那我这个可以从哪些方面考虑解决几何优化不收敛呢？现在2.0E-5的精度SCF都收敛不了的话，SCF收敛限 ...

那个季鏻盐的结构是饱和的，应该不是长链分子的一段，而且季鏻盐XY方向都伸出盒子外面。你开了XY方向的PBC，很可能本来的季鏻盐跟PBC重复过来的季鏻盐挨得很近，这样结构就不合理了。建议把季鏻盐包括在盒子里，并检查一下PBC条件下是否有原子挨得很近的情况。

sobereva 发表于 2022-8-6 03:33
不了解CP2K不要乱改。Multiwfn产生的CP2K几何优化输入文件里设的EPS_SCF 2.0E-06对大多数体系都是很合适的 ...

老师，那我这个可以从哪些方面考虑解决几何优化不收敛呢？现在2.0E-5的精度SCF都收敛不了的话，SCF收敛限设得更严（同时适当减小EPS_DEFAULT）还需要尝试吗

Eudaimonia 发表于 2022-8-5 18:54
SCF收敛限放的太松了，会导致几何优化收敛反复震荡，先提高一个数量级再看看

老师，计算结果出来了，改成2.0E-05后SCF就收敛不了了：
Leaving inner SCF loop after reaching   200 steps.


  Electronic density on regular grids:       -539.9999998622        0.0000001378
  Core density on regular grids:              539.9999999787       -0.0000000213
  Total charge density on r-space grids:        0.0000001165
  Total charge density g-space grids:           0.0000001165

  Overlap energy of the core charge distribution:               0.00415062713215
  Self energy of the core charge distribution:              -3730.21242659181553
  Core Hamiltonian energy:                                   1020.37104852121581
  Hartree energy:                                            1295.47022113939875
  Exchange-correlation energy:                               -304.36239004881151
  Dispersion energy:                                           -0.40372187706263

  Total energy:                                             -1719.13311822994274

*** WARNING in qs_scf.F:539 :: SCF run NOT converged ***

谢谢老师的提醒，当时是因为用2.0E-06没有收敛，有同学建议降低精度算了一下。结果就出现了上面的情况。

不了解CP2K不要乱改。Multiwfn产生的CP2K几何优化输入文件里设的EPS_SCF 2.0E-06对大多数体系都是很合适的值。几何优化不收敛当然要用几何优化不收敛的解决方法，改SCF方面设置毫无意义，而且你还把SCF计算精度改糙了，只会妨碍几何优化收敛。从SCF收敛层面增加几何优化收敛的可能性那也是把SCF收敛限设得更严（需同时适当减小EPS_DEFAULT）。

Eudaimonia 发表于 2022-8-5 18:54
SCF收敛限放的太松了，会导致几何优化收敛反复震荡，先提高一个数量级再看看

好的，感谢！我现在就试试

SCF收敛限放的太松了，会导致几何优化收敛反复震荡，先提高一个数量级再看看