计算化学公社

 找回密码 Forget password
 注册 Register
Views: 421|回复 Reply: 6
打印 Print 上一主题 Last thread 下一主题 Next thread

[新手求助] 使用混合基组优化高价碘时结构出错

[复制链接 Copy URL]

10

帖子

0

威望

51

eV
积分
61

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
我在优化1-羟基-1,2-苯并碘氧代-3(1H)-酮时,C H O使用6-311g+(d,p),I使用def2-SVP(d),优化没有报错且收敛,但是优化后结构会乱掉。输入文件和输出文件在附件中,谢谢各位老师的解答。

660ad490a9057406245c2d5f2a89c73.png (43.48 KB, 下载次数 Times of downloads: 15)

660ad490a9057406245c2d5f2a89c73.png

输出文件.rar

2.03 MB, 下载次数 Times of downloads: 4

131.gjf

1.11 KB, 下载次数 Times of downloads: 12

输入文件

4104

帖子

4

威望

8869

eV
积分
13053

Level 6 (一方通行)

MOKIT开发者

2#
发表于 Post on 2025-5-15 16:17:05 | 只看该作者 Only view this author
本帖最后由 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

  5. title

  6. 坐标

  7. C H O 0
  8. 6-311G(d,p)
  9. ****
  10. I 0
  11. def2TZVP
  12. ****

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

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

谈谈“计算时是否需要加DFT-D3色散校正?”
http://sobereva.com/413
自动做多参考态计算的程序MOKIT

10

帖子

0

威望

51

eV
积分
61

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2025-5-15 16:26:06 | 只看该作者 Only view this author
本帖最后由 bmj 于 2025-5-15 17:18 编辑

好的,谢谢老师建议,我去仔细阅读
老师,后续我按照您给的改了gjf文件,但是结果还是结构不对呢

214

帖子

0

威望

889

eV
积分
1103

Level 4 (黑子)

4#
发表于 Post on 2025-5-15 18:35:14 | 只看该作者 Only view this author
bmj 发表于 2025-5-15 16:26
好的,谢谢老师建议,我去仔细阅读
老师,后续我按照您给的改了gjf文件,但是结果还是结构不对呢

人家已经和耐心的回复你了
别让别人再问你“结果不对是哪里不对”

310

帖子

0

威望

461

eV
积分
771

Level 4 (黑子)

5#
发表于 Post on 2025-5-15 19:38:06 | 只看该作者 Only view this author
bmj 发表于 2025-5-15 16:26
好的,谢谢老师建议,我去仔细阅读
老师,后续我按照您给的改了gjf文件,但是结果还是结构不对呢

帮你算了一下。自己看看吧。可以的话给我加点ev。

I131.rar

320.43 KB, 下载次数 Times of downloads: 5

out文件

10

帖子

0

威望

51

eV
积分
61

Level 2 能力者

6#
 楼主 Author| 发表于 Post on 2025-5-21 14:45:07 | 只看该作者 Only view this author
18217265596 发表于 2025-5-15 18:35
人家已经和耐心的回复你了
别让别人再问你“结果不对是哪里不对”

好的,我下次注意

10

帖子

0

威望

51

eV
积分
61

Level 2 能力者

7#
 楼主 Author| 发表于 Post on 2025-5-21 14:45:30 | 只看该作者 Only view this author
tonganlhy 发表于 2025-5-15 19:38
帮你算了一下。自己看看吧。可以的话给我加点ev。

谢谢

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-14 06:49 , Processed in 0.159289 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list