计算化学公社

标题: 求助：Co掺杂GDC超胞进行对角化计算时读取OT计算的波函数不收敛 [打印本页]

作者
Author: 1372133091 时间: 2024-9-28 17:05
标题: 求助：Co掺杂GDC超胞进行对角化计算时读取OT计算的波函数不收敛
各位老师好！我在使用CP2K计算Co掺杂GDC超胞的态密度时，对角化计算过程中无法收敛，DFT+U，磁性都考虑了，仍不收敛，请问有比较好的解决方案吗？
输入、输出文件如下：

作者
Author: liuxp 时间: 2024-9-28 18:11
本帖最后由 liuxp 于 2024-9-28 18:12 编辑

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

作者
Author: 1372133091 时间: 2024-9-28 22:24

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

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

作者
Author: liuxp 时间: 2024-9-28 22:28

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

啥意思

作者
Author: 1372133091 时间: 2024-9-28 22:49
本帖最后由 1372133091 于 2024-9-28 22:58 编辑

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

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

作者
Author: 1372133091 时间: 2024-9-28 23:54

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

好的，谢谢Sob老师

作者
Author: liuxp 时间: 2024-9-29 15:49

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

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

作者
Author: 1372133091 时间: 2024-9-29 17:10
本帖最后由 1372133091 于 2024-9-29 17:26 编辑

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显示的问题，但是个人感觉是这个导致的，因为其影响了正常的周期性。至于解决办法现在还没有一个比较好的办法。

作者
Author: liuxp 时间: 2024-9-29 19:56

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

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

作者
Author: 1372133091 时间: 2024-9-29 22:30

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

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

作者
Author: liuxp 时间: 2024-9-30 08:11

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

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

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