计算化学公社

标题: 使用混合基组优化高价碘时结构出错 [打印本页]
作者
Author: bmj    时间: 2025-5-15 16:02
标题: 使用混合基组优化高价碘时结构出错
我在优化1-羟基-1,2-苯并碘氧代-3（1H）-酮时，C H O使用6-311g+(d,p)，I使用def2-SVP(d)，优化没有报错且收敛，但是优化后结构会乱掉。输入文件和输出文件在附件中，谢谢各位老师的解答。
作者
Author: zjxitcc    时间: 2025-5-15 16:17
本帖最后由 zjxitcc 于 2025-5-15 16:21 编辑

（1）当你在#p这一行写B3LYP/def2SVP时，坐标底下的自定义基组和赝势就没有用了，不会被读取，因此你gjf文件本身不合理，表明你不知道应该用什么基组/赝势、应该怎么写。
（2）据我所知，不存在def2-SVP(d)，可能是您“自己创造”的基组。高价碘的I原子是超价的，因此需要质量较好的基组（包含角动量较高的极化函数），不适合用def2SVP，使用def2TZVP是一个较好的选择。
（3）不存在6-311G+(d,p)，加号应当写在G前面，即6-311+G(d,p)。这么一个结构优化完全没必要使用弥散函数，H/C/O使用6-311G(d,p)足矣。
（4）这么一个简单的结构优化，无需calcall，有点浪费时间，calcfc即可。

综上所述，一个合理的gjf格式示例如下
  1. %chk=131.chk
  2. %nprocshared=20
  3. %mem=60GB
  4. #p B3LYP/genecp em=GD3BJ opt=calcfc freq

  6. title

  8. 坐标

  10. C H O 0
  11. 6-311G(d,p)
  12. ****
  13. I 0
  14. def2TZVP
  15. ****

  17. I 0
  18. def2
  19. [空行]
  20. [空行]
复制代码
注意gjf文件末尾需要放2~3行空行（已特地标出）。建议您不要看完我的回复 就急着去提交计算任务，而是继续阅读以下几篇Sob老师的博文，稍微系统学习一下，加深印象

详解Gaussian中混合基组、自定义基组和赝势基组的输入
http://sobereva.com/60

谈谈“计算时是否需要加DFT-D3色散校正？”
http://sobereva.com/413
作者
Author: bmj    时间: 2025-5-15 16:26
本帖最后由 bmj 于 2025-5-15 17:18 编辑

好的，谢谢老师建议，我去仔细阅读
老师，后续我按照您给的改了gjf文件，但是结果还是结构不对呢
作者
Author: 18217265596    时间: 2025-5-15 18:35
bmj 发表于 2025-5-15 16:26
好的，谢谢老师建议，我去仔细阅读
老师，后续我按照您给的改了gjf文件，但是结果还是结构不对呢

人家已经和耐心的回复你了
别让别人再问你“结果不对是哪里不对”
作者
Author: tonganlhy    时间: 2025-5-15 19:38
bmj 发表于 2025-5-15 16:26
好的，谢谢老师建议，我去仔细阅读
老师，后续我按照您给的改了gjf文件，但是结果还是结构不对呢

帮你算了一下。自己看看吧。可以的话给我加点ev。
作者
Author: bmj    时间: 2025-5-21 14:45
18217265596 发表于 2025-5-15 18:35
人家已经和耐心的回复你了
别让别人再问你“结果不对是哪里不对”

好的，我下次注意
作者
Author: bmj    时间: 2025-5-21 14:45
tonganlhy 发表于 2025-5-15 19:38
帮你算了一下。自己看看吧。可以的话给我加点ev。

谢谢



欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3