计算化学公社

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

[Gaussian/gview] 请问,高斯B3LYP-D3色散校正如何加快几何优化收敛?

[复制链接 Copy URL]

127

帖子

1

威望

1231

eV
积分
1378

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
刚入手高斯,版本是G09 D.01,在使用了6-311G*+优化苯酚结构时,加入GD3校正几何优化了一晚上都没有收敛,去掉了GD3,很快就能收敛了。有什么办法能加快 加入GD3校正的几何优化的收敛速度的方法吗?
另外,GD3校正似乎有解析2二阶导数, 为什么收敛会很难?
输入文件如下:
%chk=C:\Users\Administrator\Desktop\gauss\phenol.chk
# opt freq=(noraman,anharmonic) ub3lyp/6-311+g(d,p)
scrf=(iefpcm,solvent=water) pop=(nbo,savemixed) geom=connectivity
empiricaldispersion=gd3bj

Title Card Required

0 1
C                 -1.80828040   -1.21576763   -0.05801987
C                 -0.71496382   -1.03748118    0.80034763
C                 -0.01425543    0.17616183    0.79975206
C                 -0.40686511    1.21151926   -0.05920929
C                 -1.50018177    1.03323286   -0.91757669
C                 -2.20088994   -0.18041028   -0.91698137
H                 -2.34328608   -2.14241142   -0.05756597
H                 -0.41519818   -1.82799923    1.45618433
H                  0.82051657    0.31228718    1.45513409
H                 -1.79994798    1.82375124   -1.57341273
H                 -3.03566233   -0.31653540   -1.57236295
O                  0.30814201    2.44993132   -0.05981535
H                 -0.08593261    3.04664546   -0.70031315

1 2 1.5 6 1.5 7 1.0
2 3 1.5 8 1.0
3 4 1.5 9 1.0
4 5 1.5 12 1.0
5 6 1.5 10 1.0
6 11 1.0
7
8
9
10
11
12 13 1.0
13

5万

帖子

99

威望

5万

eV
积分
112492

管理员

公社社长

2#
发表于 Post on 2014-12-9 04:24:26 | 只看该作者 Only view this author
D3(BJ)没有解析二阶导数,虽然你做freq不报错,但实际上是不精确的。D3(零阻尼)才有二阶导数。所以要做freq的话应当用零阻尼。看你的输入文件貌似还要做非谐振计算,这种情况最好就别用色散校正,要考虑色散的话就用M062X这样的。
然而,不管有没有解析二阶导数,opt如果你不写calcfc/calcall,都是一样的,并不会直接计算二阶导数。

加了D3(BJ)优化不收敛巧合因素居多,并没有D3(BJ)会导致普遍性难以收敛的情况。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

127

帖子

1

威望

1231

eV
积分
1378

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2014-12-9 06:18:27 | 只看该作者 Only view this author
非常感谢 sob! 明白了。

127

帖子

1

威望

1231

eV
积分
1378

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2014-12-9 06:46:05 | 只看该作者 Only view this author
重新用GD3(去掉bj) 算了一遍,完全和sob大大说的一样, 收敛飞快,不到10分钟就算完了几何优化部分

105

帖子

0

威望

867

eV
积分
974

Level 4 (黑子)

5#
发表于 Post on 2014-12-9 11:04:00 | 只看该作者 Only view this author
Shannon 发表于 2014-12-9 06:46
重新用GD3(去掉bj) 算了一遍,完全和sob大大说的一样, 收敛飞快,不到10分钟就算完了几何优化部分

你好,我想问一下你是在溶液中(水)优化溶质的构型吗,你是加了色散矫正,是否像这种情况一定得加呢?还有就是这个empiricaldispersion=gd3bj类似的关键词在高斯手册上有吗?好像没看见。
谢谢。

5万

帖子

99

威望

5万

eV
积分
112492

管理员

公社社长

6#
发表于 Post on 2014-12-9 14:54:51 | 只看该作者 Only view this author
tjuchan 发表于 2014-12-9 11:04
你好,我想问一下你是在溶液中(水)优化溶质的构型吗,你是加了色散矫正,是否像这种情况一定得加呢?还 ...


用不用溶剂模型和是否用色散校正没直接关系。

empiricaldispersion=gd3bj关键词见此贴
DFT-D色散校正的使用
http://sobereva.com/210
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

127

帖子

1

威望

1231

eV
积分
1378

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2014-12-10 02:06:51 | 只看该作者 Only view this author
tjuchan 发表于 2014-12-9 11:04
你好,我想问一下你是在溶液中(水)优化溶质的构型吗,你是加了色散矫正,是否像这种情况一定得加呢?还 ...

手册上没有,俺是把sob大大的博客基本看完了,偶然翻到这个方法的

5万

帖子

99

威望

5万

eV
积分
112492

管理员

公社社长

8#
发表于 Post on 2014-12-10 02:41:34 | 只看该作者 Only view this author
Shannon 发表于 2014-12-10 02:06
手册上没有,俺是把sob大大的博客基本看完了,偶然翻到这个方法的

网上的pdf版手册都是老版本的,DFT-D3又是从D.01才开始加入的,所以下载的手册看不到。官网上的最新版手册里有
http://www.gaussian.com/g_tech/g_ur/k_dft.htm
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

12

帖子

0

威望

39

eV
积分
51

Level 2 能力者

9#
发表于 Post on 2023-10-18 11:33:14 | 只看该作者 Only view this author
Shannon 发表于 2014-12-10 02:06
手册上没有,俺是把sob大大的博客基本看完了,偶然翻到这个方法的

请问sob大大的博客在哪看

1万

帖子

0

威望

7402

eV
积分
18171

Level 6 (一方通行)

10#
发表于 Post on 2023-10-18 16:05:24 | 只看该作者 Only view this author
dio. 发表于 2023-10-18 04:33
请问sob大大的博客在哪看

http://sobereva.com/
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

本版积分规则 Credits rule

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

GMT+8, 2024-11-27 00:51 , Processed in 0.180082 second(s), 21 queries , Gzip On.

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