计算化学公社

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

[Gaussian/gview] 使用def2-SVPD泛函计算单点始终不收敛

[复制链接 Copy URL]

83

帖子

0

威望

647

eV
积分
730

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
老师们好,想请教一下,基于图片的分子我想求其偶极距,基于采用b3lyp 6-31g(d)优化好的结构 采用def2-SVPD基组算单点 但始终不收敛。我看了帖子,尝试过以下方法:(1)SCF(novaracc,noincfock) int=acc2e=12  (2)scf(novaracc,conver=6)  (3)scf(novaracc,conver=6)   (4)SCF=noDIIS  (5)不加溶剂  (6)SCF=Fermi  这是目前尝试过的方法但仍不收敛。
看了sob的帖子 一些方法比如QC不建议使用  想请教一下各位有类似经验吗?我采取什么方法可以让其收敛呢?下一步准备采用的是换基组方法  步步为营6-31G*->6-311G**->def2-TZVP>def2-SVPD  可以吗?还有什么改进方法推荐吗?

202209100949505830..png (86.51 KB, 下载次数 Times of downloads: 13)

202209100949505830..png

202209100949303902..png (21.67 KB, 下载次数 Times of downloads: 11)

202209100949303902..png

PNDIT-F3N-Br.png (39.12 KB, 下载次数 Times of downloads: 12)

PNDIT-F3N-Br.png

83

帖子

0

威望

647

eV
积分
730

Level 4 (黑子)

32#
 楼主 Author| 发表于 Post on 2022-9-19 08:02:28 | 只看该作者 Only view this author
zjxitcc 发表于 2022-9-18 13:50
你也可以尝试6-31+G(d,p),6-311+G(d,p)之类的,不一定非要def2-SVPD吧。。

好的 谢谢老师

83

帖子

0

威望

647

eV
积分
730

Level 4 (黑子)

31#
 楼主 Author| 发表于 Post on 2022-9-19 08:02:18 | 只看该作者 Only view this author
scf 发表于 2022-9-16 12:55
几何优化用过色散修正吗?做些构象搜索?有可能分子处于一个相对比较多参考态的状态SCF收敛困难,但这个位 ...

好的谢谢老师

4104

帖子

4

威望

8869

eV
积分
13053

Level 6 (一方通行)

MOKIT开发者

30#
发表于 Post on 2022-9-18 13:50:08 | 只看该作者 Only view this author
兰兰爱学习 发表于 2022-9-16 11:20
请教一下各位老师 基于1L 采取的不收敛方法 又根据老师的建议采取读取之前优化成功的 check 文件里的几何结 ...

你也可以尝试6-31+G(d,p),6-311+G(d,p)之类的,不一定非要def2-SVPD吧。。
自动做多参考态计算的程序MOKIT

194

帖子

0

威望

3694

eV
积分
3888

Level 5 (御坂)

29#
发表于 Post on 2022-9-16 12:55:04 | 只看该作者 Only view this author
几何优化用过色散修正吗?做些构象搜索?有可能分子处于一个相对比较多参考态的状态SCF收敛困难,但这个位置可能不是全局极小值。虽然qc比较慢,没有办法的时候也可以试试

83

帖子

0

威望

647

eV
积分
730

Level 4 (黑子)

28#
 楼主 Author| 发表于 Post on 2022-9-16 11:20:34 | 只看该作者 Only view this author
请教一下各位老师 基于1L 采取的不收敛方法 又根据老师的建议采取读取之前优化成功的 check 文件里的几何结构和电子结构作为初始猜测 进行计算 仍是不收敛  请问是此体系不适用于这个基组吗?还是学生的操作有问题? 以下是采取初猜方法的报错输出文件。

ljl-SVPD-a-sp-PNDIT-F3N-Br.log

196.79 KB, 下载次数 Times of downloads: 6

83

帖子

0

威望

647

eV
积分
730

Level 4 (黑子)

27#
 楼主 Author| 发表于 Post on 2022-9-15 14:16:29 | 只看该作者 Only view this author
zjxitcc 发表于 2022-9-15 14:08
我在24L回复的意思是仅把关键词那一行中的sp去掉。但你不仅把关键词那一行中的sp去掉了,还把Title行 ...

好的老师明白了 谢谢老师

4104

帖子

4

威望

8869

eV
积分
13053

Level 6 (一方通行)

MOKIT开发者

26#
发表于 Post on 2022-9-15 14:08:45 | 只看该作者 Only view this author
兰兰爱学习 发表于 2022-9-15 13:39
老师 我改了以后报错显示gjf格式有问题 请问是哪块我写错了吗?

我在24L回复的意思是仅把关键词那一行中的sp去掉。但你不仅把关键词那一行中的sp去掉了,还把Title行给删了(因为你原来Title行里也含有sp字符)。要保留Title行,但Title行内容可以随便写。

如果你不懂啥是Title行,可以对比22L与25L的gjf文件
自动做多参考态计算的程序MOKIT

83

帖子

0

威望

647

eV
积分
730

Level 4 (黑子)

25#
 楼主 Author| 发表于 Post on 2022-9-15 13:39:38 | 只看该作者 Only view this author
zjxitcc 发表于 2022-9-15 12:39
去掉sp,这没用。“0 1”与“H     0”间要空一行。

老师 我改了以后报错显示gjf格式有问题 请问是哪块我写错了吗?

baocuo.zip

4.2 KB, 下载次数 Times of downloads: 8

4104

帖子

4

威望

8869

eV
积分
13053

Level 6 (一方通行)

MOKIT开发者

24#
发表于 Post on 2022-9-15 12:39:50 | 只看该作者 Only view this author
兰兰爱学习 发表于 2022-9-15 09:10
好的老师已经上传在下方  请问老师输入文件是哪部分出错了导致的报错呢?
/home/user5/ljl/jiarongji/  ...

去掉sp,这没用。“0 1”与“H     0”间要空一行。
自动做多参考态计算的程序MOKIT

83

帖子

0

威望

647

eV
积分
730

Level 4 (黑子)

23#
 楼主 Author| 发表于 Post on 2022-9-15 09:12:19 | 只看该作者 Only view this author
乐平 发表于 2022-9-14 15:58
你提交任务就知道对不对了。

能正常运行起来就是对的。

好的老师 已经提交任务发现报错 麻烦请教一下老师是哪里出错了呢?

baocuo.zip

4.2 KB, 下载次数 Times of downloads: 4

83

帖子

0

威望

647

eV
积分
730

Level 4 (黑子)

22#
 楼主 Author| 发表于 Post on 2022-9-15 09:10:07 | 只看该作者 Only view this author
zjxitcc 发表于 2022-9-14 19:14
不要复制粘贴(很容易出现内容不全、内容繁杂、空格、空行不对等情况),直接上传输入文件简单了事。要是 ...

好的老师已经上传在下方  请问老师输入文件是哪部分出错了导致的报错呢?
/home/user5/ljl/jiarongji/ 是文件所在服务器的位置   读取的文件是优化好的ljl-b-a-opt-PNDIT-F3N-Br.chk 文件

baocuo.zip

4.2 KB, 下载次数 Times of downloads: 4

4104

帖子

4

威望

8869

eV
积分
13053

Level 6 (一方通行)

MOKIT开发者

21#
发表于 Post on 2022-9-14 19:14:46 | 只看该作者 Only view this author
兰兰爱学习 发表于 2022-9-14 08:28
老师请问这样改对吗?
%chk=/home/user5/ljl/jiarongji/ljl-b-a-opt-PNDIT-F3N-Br.chk
%nprocshared=24 ...

不要复制粘贴(很容易出现内容不全、内容繁杂、空格、空行不对等情况),直接上传输入文件简单了事。要是体系需要保密,可以仅把坐标删去,其他部分保留。
自动做多参考态计算的程序MOKIT

1104

帖子

0

威望

3954

eV
积分
5058

Level 6 (一方通行)

20#
发表于 Post on 2022-9-14 15:58:48 | 只看该作者 Only view this author
兰兰爱学习 发表于 2022-9-14 08:28
老师请问这样改对吗?
%chk=/home/user5/ljl/jiarongji/ljl-b-a-opt-PNDIT-F3N-Br.chk
%nprocshared=24 ...

你提交任务就知道对不对了。

能正常运行起来就是对的。

不能正常运行就会报错。

83

帖子

0

威望

647

eV
积分
730

Level 4 (黑子)

19#
 楼主 Author| 发表于 Post on 2022-9-14 08:28:17 | 只看该作者 Only view this author
本帖最后由 兰兰爱学习 于 2022-9-14 08:30 编辑
乐平 发表于 2022-9-13 10:55
1)没看到图……

2)不需要下载 .chk 文件啊,你只有告诉 g16 你的 .chk 文件的路径和文件名就行了啊
...

老师请问这样改对吗?
%chk=/home/user5/ljl/jiarongji/ljl-b-a-opt-PNDIT-F3N-Br.chk
%nprocshared=24
%mem=24GB
#p sp  B3LYP/gen geom=check guess=read     EMPIRICALDISPERSION=GD3BJ  scrf=(pcm,solvent=Acetonitrile)  

sp

0 1
H     0
S    3   1.00
     13.0107010              0.19682158D-01
      1.9622572              0.13796524
      0.44453796             0.47831935
S    1   1.00
      0.12194962             1.0000000
P    1   1.00
      0.8000000              1.0000000
P    1   1.00
      0.11704099050          1.0000000
****
C     0
S    5   1.00
   1238.4016938              0.54568832082D-02
    186.29004992             0.40638409211D-01
     42.251176346            0.18025593888
     11.676557932            0.46315121755
      3.5930506482           0.44087173314
S    1   1.00
      0.40245147363          1.0000000
S    1   1.00
      0.13090182668          1.0000000
S    1   1.00
      0.67053540256D-01      1.0000000
P    3   1.00
      9.4680970621           0.38387871728D-01
      2.0103545142           0.21117025112
      0.54771004707          0.51328172114
P    1   1.00
      0.15268613795          1.0000000
D    1   1.00
      0.8000000              1.0000000
D    1   1.00
      0.11713185140          1.0000000
****
~~~~·
谢谢老师指教

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

GMT+8, 2025-8-15 09:38 , Processed in 0.214314 second(s), 24 queries , Gzip On.

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