求助：Co掺杂GDC超胞进行对角化计算时读取OT计算的波函数不收敛

1372133091

各位老师好！我在使用CP2K计算Co掺杂GDC超胞的态密度时，对角化计算过程中无法收敛，DFT+U，磁性都考虑了，仍不收敛，请问有比较好的解决方案吗？
输入、输出文件如下：

liuxp

直接用OT,我那个就是镧系元素多，对角化不收敛用的OT收敛的

1372133091

liuxp 发表于 2024-9-28 18:11
直接用OT,我那个就是镧系元素多，对角化不收敛用的OT收敛的

请问用OT收敛后，我发现.log文件中有相关信息，需要怎么分析呢？

liuxp

1372133091 发表于 2024-9-28 22:24
请问用OT收敛后，我发现.log文件中有相关信息，需要怎么分析呢？

啥意思

1372133091

liuxp 发表于 2024-9-28 22:28
啥意思

emmmm意思是我得到的输出文件中没有.molden文件，该怎么分析，还是说我的输入文件不正确
我添加的指令如下添加之后.log文件中有DOS信息了，但是不知道怎么分析。
如果加了&MO_MOLDEN部分的话，会报错“OT with ADDED_MOS/=0 not implemented
&PRINT
    &MO
       ENERGIES T
       OCCUPATION _NUMBERS T
       &EACH
          QS_SCF 0
       &END EACH
     &END MO

sobereva

开OT时产生molden文件的注意事项下文都说了
详谈使用CP2K产生给Multiwfn用的molden格式的波函数文件
http://sobereva.com/651（http://bbs.keinsci.com/thread-31827-1-1.html）

1372133091

sobereva 发表于 2024-9-28 23:36
开OT时产生molden文件的注意事项下文都说了
详谈使用CP2K产生给Multiwfn用的molden格式的波函数文件
http ...

好的，谢谢Sob老师

liuxp

1372133091 发表于 2024-9-28 23:54
好的，谢谢Sob老师

你这个收敛了吗，用OT的波函数当初猜

1372133091

liuxp 发表于 2024-9-29 15:49
你这个收敛了吗，用OT的波函数当初猜

如果只是针对GDC的3X3X3超胞，能够收敛（见附件），但是我掺Co之后的GDC超胞在对角化时始终无法收敛。
且读取OT计算的波函数在进行对角化计算得到的DOS更合理，同时也能够读取PDOS。
现在发现仅使用OT计算得到的DOS中，带隙变得太宽了，比实验值多了将近一倍了（见下图，两张图的FWHM不一致，请见谅），目前还不知道怎么解决，
我个人的看法：Co掺杂后无法收敛的根本原因还是结构不合理导致的：结构优化后超胞发生了很大的变形（见我之前的帖子，结构优化后变形，http://bbs.keinsci.com/thread-48405-1-1.html），尽管看起来可能只是VESTA显示的问题，但是个人感觉是这个导致的，因为其影响了正常的周期性。至于解决办法现在还没有一个比较好的办法。

liuxp

1372133091 发表于 2024-9-29 17:10
如果只是针对GDC的3X3X3超胞，能够收敛（见附件），但是我掺Co之后的GDC超胞在对角化时始终无法收敛。
...

我最近也在苦恼的这个问题，确实不好解决

1372133091

liuxp 发表于 2024-9-29 19:56
我最近也在苦恼的这个问题，确实不好解决

诶，我有一计，不妨使用杂化泛函试试，我先去试一试

liuxp

1372133091 发表于 2024-9-29 22:30
诶，我有一计，不妨使用杂化泛函试试，我先去试一试

杂化泛函 咱俩的上百原子数 算的动吗