关于计算硫化铜体系的参数选取问题

kirSTY

本人在计算硫化铜晶体001晶面结构优化的时候遇到了SCF不收敛的问题，在阅读社长的“SCF不收敛的处理方法”那篇帖子后进行了一系列调整和测试。
目前发现自旋极化的参数选取对收敛性影响较大，此外想请教各位大佬：

1.Cu表现为+2价时有单电子，针对CuS体系的计算是否需要开UKS？
2.了解到CuS体系好像是强相关体系，是否需要DFT+U？
3.自旋极化目前经过计算测试发现取5的时候相对容易收敛，可否认为这个体系就应该选取5？
4.计算方面取原胞+3，3，1的K点不收敛，取超胞+OT可以收敛，但是优化后的结构与原始结构相差较大(出现错层)，个人觉得还是计算参数没有取对，关于这个希望各位大佬提供一点思路？
5.其他是否有需要注意的地方，比如基组和泛函的选取，甚至是更换计算软件？目前计算软件是CP2K，选用的是multiwfn默认的PBE泛函和DZVP基组

谢谢各位大佬！

sobereva

我没写过叫“SCF不收敛的处理方法”的帖子
帖子标题必须准确，并且给出具体链接，就像写论文时要恰当引用一样
CP2K中遇到SCF难收敛时的解决方法
http://sobereva.com/665（http://bbs.keinsci.com/thread-37196-1-1.html）

sobereva

1 如果你确认是开壳层就需要。但是光写个UKS没用，必须打破初猜波函数对称性，参考论坛里CP2K的磁性体系计算的讨论，我回复过很多次

2 这么典型的体系文献里早就计算过，参考文献

3 SCF是否容易收敛和什么状态是基态毫无必然关系。显然要对比不同磁性状态的能量

4 SCF收敛问题博文里该说的都说了

5 基组名必须写得非常准确，DZVP-GTH和DZVP-MOLOPT-SR-GTH差异大了去了，后者精度和耗时远高于前者。不需要换程序

kirSTY

好的谢谢社长！