请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 5680|回复: 11

[Gaussian/gview] 求助,高斯B3LYP/PCM计算双原子自由基SCF不收敛

[复制链接]

127

帖子

1

威望

1213

eV
积分
1360

Level 4 (黑子)

发表于 2014-10-27 00:35:16 | 显示全部楼层 |阅读模式
大家好,有个问题想请教。 用B3LYP/g-311++g**/SCRF(PCM)计算双原子的 自由基,比如OH自由基和NO自由基,总是SCF不收敛。而计算大点的分子,比如NO2自由基,就收敛了。不知道这个是什么原因。
另外,用了一组sob先辈的解决方案,把SCF改成SCF=(noincfock,conver=6,maxcyc=200) iop(5/22=20),请问这个计算结果可靠吗?
另外,查了下高斯的官网, iop(5/22=20) 指令的意思是 Energy-based mixing,可否解释下这个指令 是什么意思?
高斯input文件:
#p B3LYP/Integral(Grid=UltraFine)/6-311++G(d,p)
SCF=(noincfock,conver=6,maxcyc=200) iop(5/22=20)
SCRF(solvent=water)
OPT
GEOM(PrintInputOrient)
NoSymm

NO Geometry Optimization

0 2
N -4.447667 1.287769 -0.000000
O -2.973329 1.417097 0.000000

1万

帖子

25

威望

2万

eV
积分
42611

管理员

公社社长

发表于 2014-10-27 00:57:14 | 显示全部楼层
你写了很多关键词没什么用。我用下面这个输入文件很顺利地就算完了
#p B3LYP/6-311+G(d,p) scrf=PCM opt

NO Geometry Optimization

0 2
N -4.447667 1.287769 -0.000000
O -2.973329 1.417097 0.000000

没必要用5/22
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。达5000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

127

帖子

1

威望

1213

eV
积分
1360

Level 4 (黑子)

 楼主| 发表于 2014-10-27 01:08:27 | 显示全部楼层
sobereva 发表于 2014-10-27 00:57
你写了很多关键词没什么用。我用下面这个输入文件很顺利地就算完了
#p B3LYP/6-311+G(d,p) scrf=PCM opt
...

可能是程序bug吧,我用的是高斯03,刚才试了下用先辈发的input, 迭代最大次数超过了还是没收敛。

最后一轮迭代的输出信息是这个样子的。

Cycle 128  Pass 1  IDiag  1:
With in-vacuo wavefunction:
Error on total polarization charges =  0.01054
E= -129.825476662981     Delta-E=        0.000169721250 Rises=F Damp=F
DIIS: error= 5.42D-03 at cycle 128 NSaved=  15.
NSaved=15 IEnMin= 2 EnMin= -129.829138030045     IErMin= 3 ErrMin= 1.64D-03
ErrMax= 5.42D-03 EMaxC= 1.00D-01 BMatC= 1.33D-03 BMatP= 7.14D-05
IDIUse=3 WtCom= 9.91D-01 WtEn= 9.28D-03
Large coefficients: NSaved= 15 BigCof=   28.21 CofMax=   10.00 Det=-4.07D-11
Inversion failed.  Reducing to 14 matrices.
Rare condition: small coef for last iteration:  0.000D+00
Coeff-Com:  0.653D-01 0.562D-02 0.288D+01-0.325D+01-0.345D+01 0.974D+01
Coeff-Com: -0.200D+01 0.890D-01 0.136D+01 0.326D+00-0.193D+01-0.232D+01
Coeff-Com: -0.401D+01 0.349D+01
Coeff-En:   0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00
Coeff-En:   0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00
Coeff-En:   0.000D+00 0.100D+01
Coeff:      0.653D-01 0.562D-02 0.288D+01-0.325D+01-0.345D+01 0.974D+01
Coeff:     -0.200D+01 0.890D-01 0.136D+01 0.326D+00-0.193D+01-0.232D+01
Coeff:     -0.401D+01 0.349D+01
Gap=     0.119 Goal=   None    Shift=    0.000
Gap=     0.190 Goal=   None    Shift=    0.000
RMSDP=8.84D-04 MaxDP=1.38D-02 DE= 1.70D-04 OVMax= 5.81D-02

>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(UB+HF-LYP) =  -129.825476663     A.U. after  129 cycles
             Convg  =    0.8839D-03             -V/T =  2.0073
             S**2   =   0.8419
KE= 1.288843157076D+02 PE=-3.456293013007D+02 EE= 6.689657962702D+01

1万

帖子

25

威望

2万

eV
积分
42611

管理员

公社社长

发表于 2014-10-27 01:15:52 | 显示全部楼层
我用的G09 B.01
G09的SCF收敛方法比G03处理得好
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。达5000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

225

帖子

0

威望

1090

eV
积分
1315

Level 4 (黑子)

发表于 2015-1-10 22:10:53 | 显示全部楼层
sobereva 发表于 2014-10-27 00:57
你写了很多关键词没什么用。我用下面这个输入文件很顺利地就算完了
#p B3LYP/6-311+G(d,p) scrf=PCM opt
...

弱弱的问一下,这是进行溶剂化模型下的优化吗?但是没有指定溶剂。
抱歉,我没做过溶剂中的结构优化,只是问一下。

127

帖子

1

威望

1213

eV
积分
1360

Level 4 (黑子)

 楼主| 发表于 2015-1-10 22:28:25 | 显示全部楼层
虎王 发表于 2015-1-10 22:10
弱弱的问一下,这是进行溶剂化模型下的优化吗?但是没有指定溶剂。
抱歉,我没做过溶剂中的结构优化,只 ...

scrf默认是溶剂是水。如果模拟的是水环境下的反应,有时候忘了写溶剂 也没啥关系,

225

帖子

0

威望

1090

eV
积分
1315

Level 4 (黑子)

发表于 2015-1-10 22:44:46 | 显示全部楼层
Shannon 发表于 2015-1-10 22:28
scrf默认是溶剂是水。如果模拟的是水环境下的反应,有时候忘了写溶剂 也没啥关系,

原来如此,非常感谢。

1528

帖子

3

威望

5604

eV
积分
7192

Level 6 (一方通行)

Ab Initio Amateur

发表于 2015-1-10 23:00:56 | 显示全部楼层
如果只是指定PCM,G09和G03的处理方式好像是有点差异的。
Performs a reaction field calculation using the integral equation formalism model (IEFPCM). This is the default. Some details of the formalism and the implementation have changed with respect to Gaussian 03, as described in [Scalmani10]. IEFPCM is a synonym for PCM.
我不确定是不是这个原因,换成G09试试。
满招损,谦受益。热衷于理论和方法研究水平不高但欢迎讨论。

127

帖子

1

威望

1213

eV
积分
1360

Level 4 (黑子)

 楼主| 发表于 2015-1-10 23:21:35 | 显示全部楼层
卡开发发 发表于 2015-1-10 23:00
如果只是指定PCM,G09和G03的处理方式好像是有点差异的。
Performs a reaction field calculation using t ...

恩,谢谢。 g09是可以收敛的。 这个问题出现之后查过文献,g09默认的溶剂化计算方式改变了,有文献报道说精度提升了不少。

16

帖子

0

威望

81

eV
积分
97

Level 2 能力者

发表于 2019-4-15 15:39:26 | 显示全部楼层
我是一个新手,我有个问题想请教一下您,就是关于自由基优化的问题,我看您NO自由基的电荷和多重度与NO中性分子的是一样的,那这不就是想当于一样的吗?怎么体现出来你在优化的是自由基,我看您的关键词就是日常优化的关键词差不多,多以我一直不理解,优化自由基,是怎么来确定它的自旋多重度的,怎么体现是在优化自由基而不是中性分子呢?

1万

帖子

25

威望

2万

eV
积分
42611

管理员

公社社长

发表于 2019-4-16 00:11:36 | 显示全部楼层
TJNina 发表于 2019-4-15 15:39
我是一个新手,我有个问题想请教一下您,就是关于自由基优化的问题,我看您NO自由基的电荷和多重度与NO中性 ...

NO中性分子 = NO自由基
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。达5000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

16

帖子

0

威望

81

eV
积分
97

Level 2 能力者

发表于 2019-4-16 15:49:41 | 显示全部楼层
sobereva 发表于 2019-4-16 00:11
NO中性分子 = NO自由基

谢谢老师,人真的很好
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2019-6-20 16:16 , Processed in 0.184708 second(s), 23 queries .

快速回复 返回顶部 返回列表