计算化学公社

标题: 用BP86算NiH-三重态时波函数跳动，导致三重态的能量高，输入文件加啥限制波h... [打印本页]

作者
Author: 掌上明珠 时间: 2021-5-4 10:29
标题: 用BP86算NiH-三重态时波函数跳动，导致三重态的能量高，输入文件加啥限制波h...
本帖最后由掌上明珠于 2021-5-4 10:33 编辑

用BP86算NiH-三重态时波函数跳动，导致三重态的能量高，输入文件加啥限制波函数不跳动
Initial guess orbital symmetries:
Alpha Orbitals:
   Occupied  (SG) (SG) (PI) (PI) (SG) (SG) (PI) (PI) (SG) (SG)
               (DLTA) (DLTA) (PI) (PI) (SG) (SG)
   Virtual (PI) (PI) (SG) (SG) (PI) (PI) (DLTA) (DLTA) (SG)
               (SG) (PI) (PI) (SG) (PI) (PI) (SG) (DLTA) (DLTA)
               (PI) (PI) (SG) (SG) (SG) (?A) (?A) (?A) (?A) (SG)
               (DLTA) (DLTA) (?A) (?A) (SG) (DLTA) (DLTA) (PI)
               (PI) (PI) (PI) (SG) (SG) (DLTA) (DLTA) (PI) (PI)
               (SG) (?B) (?B) (?C) (?C) (DLTA) (DLTA) (?C) (?C)
               (?B) (?C) (?C) (DLTA) (DLTA) (?C) (?C) (SG)
Beta  Orbitals:
   Occupied  (SG) (SG) (PI) (PI) (SG) (SG) (PI) (PI) (SG) (SG)
               (DLTA) (PI) (PI) (SG)
   Virtual (DLTA) (SG) (PI) (PI) (SG) (SG) (DLTA) (PI) (PI)
               (DLTA) (SG) (SG) (PI) (PI) (SG) (PI) (PI) (SG)
               (DLTA) (DLTA) (PI) (PI) (SG) (SG) (SG) (?A) (?A)
               (?A) (?A) (SG) (DLTA) (DLTA) (?A) (?A) (SG) (DLTA)
               (DLTA) (PI) (PI) (PI) (PI) (SG) (SG) (DLTA) (DLTA)
               (PI) (PI) (SG) (?B) (?C) (?C) (DLTA) (?B) (DLTA)
               (?C) (?C) (?B) (DLTA) (?C) (?C) (DLTA) (?C) (?C)
               (SG)

Orbital symmetries:
Alpha Orbitals:
   Occupied  (SG) (SG) (PI) (PI) (SG) (SG) (PI) (PI) (SG) (SG)
               (DLTA) (DLTA) (PI) (PI) (SG) (SG)
   Virtual (PI) (PI) (SG) (SG) (PI) (PI) (DLTA) (DLTA) (SG)
               (SG) (PI) (PI) (SG) (PI) (PI) (SG) (DLTA) (DLTA)
               (PI) (PI) (SG) (SG) (SG) (?A) (?A) (?A) (?A) (SG)
               (DLTA) (DLTA) (?A) (?A) (SG) (DLTA) (DLTA) (PI)
               (PI) (PI) (PI) (SG) (SG) (DLTA) (DLTA) (PI) (PI)
               (SG) (?B) (?B) (?C) (?C) (DLTA) (DLTA) (?C) (?C)
               (?B) (?C) (?C) (DLTA) (DLTA) (?C) (?C) (SG)
Beta  Orbitals:
   Occupied  (SG) (SG) (PI) (PI) (SG) (SG) (PI) (PI) (SG) (SG)
               (DLTA) (PI) (PI) (SG)
   Virtual (DLTA) (SG) (PI) (PI) (SG) (SG) (DLTA) (PI) (PI)
               (DLTA) (SG) (SG) (PI) (PI) (SG) (PI) (PI) (SG)
               (DLTA) (DLTA) (PI) (PI) (SG) (SG) (SG) (?A) (?A)
               (?A) (?A) (SG) (DLTA) (DLTA) (?A) (?A) (SG) (DLTA)
               (DLTA) (PI) (PI) (PI) (PI) (SG) (SG) (DLTA) (DLTA)
               (PI) (PI) (SG) (?B) (?C) (?C) (DLTA) (?B) (DLTA)
               (?C) (?C) (?B) (DLTA) (?C) (?C) (DLTA) (?C) (?C)

作者
Author: zjxitcc 时间: 2021-5-4 10:55
本帖最后由 zjxitcc 于 2021-5-4 10:58 编辑

不用管啥“波函数跳动”，当前关键词（严重）不合理。举一种输入文件的合理示例

%chk=NiO.chk
%mem=64GB
%nprocshared=32
#p UBP86/def2TZVPP nosymm stable=opt
title
-1 3
Ni 0.0 0.0 0.0
H 0.0 0.0 1.59077
--Link1--
%chk=NiO.chk
%mem=64GB
%nprocshared=32
#p opt=calcfc freq UBP86/def2TZVPP nosymm guess=read geom=allcheck

复制代码

咋一步就结构优化收敛了。。你这是别人优化好的结构是吧。。。

PS：这是阴离子，最好用ma-def2-TZVP基组。

作者
Author: 掌上明珠 时间: 2021-5-4 18:13

zjxitcc 发表于 2021-5-4 10:55
不用管啥“波函数跳动”，当前关键词（严重）不合理。举一种输入文件的合理示例

咋一步就结构优化收敛了 ...

您好老师，NIH负离子应该是单重态能量高，我现在用BP86泛函def2TZVPP基组总是三重态能量高，我现在是想着用啥关键词是使三重态能量降低？

作者
Author: sobereva 时间: 2021-5-5 00:17

掌上明珠发表于 2021-5-4 18:13
您好老师，NIH负离子应该是单重态能量高，我现在用BP86泛函def2TZVPP基组总是三重态能量高，我现在是想着 ...

不要自说自话。先把其它可能存在不合理性的地方解决了再说，看2L的回复

作者
Author: 掌上明珠 时间: 2021-5-5 09:07
我现在用着二楼的输入文件进行计算，能量并没有变化

作者
Author: zjxitcc 时间: 2021-5-5 15:36
本帖最后由 zjxitcc 于 2021-5-5 15:39 编辑

掌上明珠发表于 2021-5-5 09:07
我现在用着二楼的输入文件进行计算，能量并没有变化

“NIH负离子应该是单重态能量高，我现在用BP86泛函def2TZVPP基组总是三重态能量高”这问题与该帖问题不同，属于另一个问题。我们先认定你在一开始提出的问题已被解决。

接着解决你现在说的“NIH负离子应该是单重态能量高”这个问题。该问题影响因素繁多，以你目前提供的（极为）有限的信息，不足以让回答者提供一针见血的答案。首先咱们排除“想着用啥关键词是使三重态能量降低”这个馊主意。大概有几种可能的原因（我们先假设你说的“NIH负离子应该是单重态能量高”这话是对的）：

（1）你算单重态时输入文件有误，例如你同样用了一些不合理的关键词（如果你想解决这个问题，毫无疑问得晒出你的计算文件，没有文件至少也得有截图）；
（2）可能你比较的对象有误，例如应该分别取单重态局域极小点与三重态局域极小点的Gibbs自由能进行比较。你可能没做到这点。
（3）可能是当前泛函不足以精准描述NIH负离子的单-三态 gap，得用更高级的泛函，或弃用DFT改用多参考方法。
在排除（1）和（2）的问题前，请不要考虑（3）。

作者
Author: 掌上明珠 时间: 2021-5-6 09:27
这是按照您的进行修改，进行的计算

作者
Author: zjxitcc 时间: 2021-5-6 10:22
本帖最后由 zjxitcc 于 2021-5-6 10:23 编辑

掌上明珠发表于 2021-5-6 09:27
这是按照您的进行修改，进行的计算

（1）要对他人回复时，请使用“回复”功能，否则对方收不到提醒；
（2）我只指出你在1L中三重态（严重）不合理的计算关键词，没说过单重态像你这么写关键词。UDFT单重态单点计算时需加上guess=mix，看Sob老师这篇博文了解基本知识《谈谈片段组合波函数与自旋极化单重态》http://sobereva.com/82
改正之后还有问题请按照我在6L的回复继续排查。

作者
Author: 掌上明珠 时间: 2021-5-6 11:56

zjxitcc 发表于 2021-5-6 10:22
（1）要对他人回复时，请使用“回复”功能，否则对方收不到提醒；
（2）我只指出你在1L中三重态（严重） ...

好的，谢谢老师。

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