计算化学公社

 找回密码 Forget password
 注册 Register
楼主 Author: 3044305615
打印 Print 上一主题 Last thread 下一主题 Next thread

[Gaussian/gview] 计算scf不收敛,错误如下,求助。

[复制链接 Copy URL]

739

帖子

0

威望

1715

eV
积分
2454

Level 5 (御坂)

16#
发表于 Post on 2015-1-2 21:29:41 来自手机 | 只看该作者 Only view this author
用这句话试试,对我来说,如果结构没有画错,就是神器

739

帖子

0

威望

1715

eV
积分
2454

Level 5 (御坂)

17#
发表于 Post on 2015-1-2 21:31:38 来自手机 | 只看该作者 Only view this author
maxcycle=200,xqc
试试这个,好用有时候

739

帖子

0

威望

1715

eV
积分
2454

Level 5 (御坂)

18#
发表于 Post on 2015-1-3 08:55:35 | 只看该作者 Only view this author
本帖最后由 lao7 于 2015-1-3 08:59 编辑

上面,手机不太好回复。发了两次

67

帖子

0

威望

1057

eV
积分
1124

Level 4 (黑子)

19#
发表于 Post on 2015-1-5 09:00:33 | 只看该作者 Only view this author
sobereva 发表于 2015-1-2 12:50
这样得到的几何结构可能不准确(虽然一般差得不会太多),为了准确,最好还是去掉之后继续优化一下


      谢谢sob,不过我试了一个150个原子的双核Cr体系,用的方法和基组是M06/6-31G(d)/Lanl08(f) 做优化,发现用 scf=conver=6 完成了优化,在得到结构基础上,去掉scf=conver=6 的关键词,结果又不收敛了。 (第一个link 502 过程SCF收敛了,但两个Displacement判据不收敛,优化几圈后一个link502 SCF收敛失败)。
      于是我把结构做一个频率分析,准备看看存在的虚频是什么样的,好根据振动方向来微调结构,结果发现频率分析中没有虚频,这似乎说明得到的结构已经是一个稳定结构了,那这个几何结构可以使用么?(freq 作业中没有使用放宽scf收敛标准)

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

20#
发表于 Post on 2015-1-5 09:30:38 | 只看该作者 Only view this author
没有虚频不代表是极小点结构。收敛精度不高的情况依然可以没有虚频。只要结构处在极小点的二次元区域内就都是没有虚频的。

是否能用取决于你对结构要求的精度。不要求很准的话,勉强能用。

不收敛有很多办法解决,几何优化不收敛见《量子化学计算中帮助几何优化收敛的常用方法》(http://hi.baidu.com/sobereva/item/596c242087e820cedcf69a9d
SCF不收敛见《解决SCF不收敛问题的方法》(http://hi.baidu.com/sobereva/item/68b8d2866c0223ebe596e0dc
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

48

帖子

0

威望

107

eV
积分
155

Level 3 能力者

21#
 楼主 Author| 发表于 Post on 2015-1-5 20:42:29 | 只看该作者 Only view this author
lao7 发表于 2015-1-2 21:31
maxcycle=200,xqc
试试这个,好用有时候

好的 谢谢 我其实没太明白xqc的作用

67

帖子

0

威望

1057

eV
积分
1124

Level 4 (黑子)

22#
发表于 Post on 2015-1-6 10:34:39 | 只看该作者 Only view this author
sobereva 发表于 2015-1-5 09:30
没有虚频不代表是极小点结构。收敛精度不高的情况依然可以没有虚频。只要结构处在极小点的二次元区域内就都 ...

多谢sob回复,我有一个结构做几何优化的时候, 位移判据始终不收敛,力判据则是早早就收敛了,看起来是进入了一个势能平坦区,结构整体刚性还挺强的,没有大的柔链结构。那这种情况是缩小步长在周围找找还是怎么想法跳出这个势能平坦区呢?

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

23#
发表于 Post on 2015-1-6 10:41:04 | 只看该作者 Only view this author
nkallwar 发表于 2015-1-6 10:34
多谢sob回复,我有一个结构做几何优化的时候, 位移判据始终不收敛,力判据则是早早就收敛了,看起来是进 ...

平坦区里面也有极小点,不过这个极小点可能是不显著的。
先试试gdiis或缩小步长上限
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

48

帖子

0

威望

107

eV
积分
155

Level 3 能力者

24#
 楼主 Author| 发表于 Post on 2015-1-15 10:41:07 | 只看该作者 Only view this author
本帖最后由 3044305615 于 2015-1-15 10:53 编辑
sobereva 发表于 2014-12-25 22:58
就是由于你用了xqc,到了64步时普通方式迭代没收敛,自动切换到了qc方式收敛。qc方式每一步很慢,得耐心等 ...


sob 老师  如果两个原子距离过近(限域) 而导致排斥能高 计算出的能量已经超出高斯给的范围了 而导致的不收敛该怎么办?
输出文件太长了 不知道哪部分 该粘贴

命令行用了iop(5/13=1) empiricaldispersion=gd3
Variable       Old X    -DE/DX   Delta X   Delta X   Delta X     New X
                                 (Linear)    (Quad)   (Total)
    R1        2.75816       NaN   0.00000       NaN   0.00000   2.75816
    R2        4.96577       NaN   0.00000       NaN   0.00000   4.96577
    R3        4.10904       NaN   0.00000       NaN   0.00000   4.10904
    R4        4.96577       NaN   0.00000       NaN   0.00000   4.96577
    R5        4.96577       NaN   0.00000       NaN   0.00000   4.96577
    R6        4.31456       NaN   0.00000       NaN   0.00000   4.31456
    R7        4.86419       NaN   0.00000       NaN   0.00000   4.86419
    R8        4.31456       NaN   0.00000       NaN   0.00000   4.31456
    R9        4.86418       NaN   0.00000       NaN   0.00000   4.86418
   R10        4.31455       NaN   0.00000       NaN   0.00000   4.31455
   R11        4.31454       NaN   0.00000       NaN   0.00000   4.31454
   R12        4.31454       NaN   0.00000       NaN   0.00000   4.31454
   R13        4.86418       NaN   0.00000       NaN   0.00000   4.86418
   R14        4.31455       NaN   0.00000       NaN   0.00000   4.31455
   R15        4.77969       NaN   0.00000       NaN   0.00000   4.77969
   R16        4.22340       NaN   0.00000       NaN   0.00000   4.22340
   R17        4.77969       NaN   0.00000       NaN   0.00000   4.77969
   R18        4.22340       NaN   0.00000       NaN   0.00000   4.22340
还有这个表示
RMSDP=7.07D-01 MaxDP=8.51D+01 DE=      NaN OVMax= 1.00D+00

以及文件末尾
RFO could not converge Lambda in  999 iterations.
Linear search skipped for unknown reason.




6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

25#
发表于 Post on 2015-1-15 15:19:52 | 只看该作者 Only view this author
3044305615 发表于 2015-1-15 10:41
sob 老师  如果两个原子距离过近(限域) 而导致排斥能高 计算出的能量已经超出高斯给的范围了 而导致 ...

绝对不要用5/13=1,SCF都没正常收敛,之后计算的受力可能会很邪乎
先用gview观察初始结构,把严重不合理的诸如距离过近的地方调整一下(比如用扫把型图标)
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

48

帖子

0

威望

107

eV
积分
155

Level 3 能力者

26#
 楼主 Author| 发表于 Post on 2015-1-16 16:20:46 | 只看该作者 Only view this author
sobereva 发表于 2015-1-15 15:19
绝对不要用5/13=1,SCF都没正常收敛,之后计算的受力可能会很邪乎
先用gview观察初始结构,把严重不合理 ...

sob 老师 我发现我的错误是log刚开始的grimme 色散的值 就是超出高斯给的范围了,为什么呢?是不是这个值没有,才导致后来的不收敛啊,您能帮忙解答一下么?但是d3不是针对1-94号元素么?

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

27#
发表于 Post on 2015-1-16 16:23:27 | 只看该作者 Only view this author
3044305615 发表于 2015-1-16 16:20
sob 老师 我发现我的错误是log刚开始的grimme 色散的值 就是超出高斯给的范围了,为什么呢?是不是这个值 ...

去掉D3再试自然就明白了。如果能正常优化完毕(此时肯定已经消除了原子间过近的接触),再重新用D3继续优化。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

52

帖子

0

威望

1208

eV
积分
1260

Level 4 (黑子)

28#
发表于 Post on 2016-10-11 15:54:07 | 只看该作者 Only view this author
sob老师,您好!使用scf=xqc进行几何构型优化,非得把关键词stable写到输入文件(gjf)中吗?

谢谢老师!

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

29#
发表于 Post on 2016-10-11 16:55:05 | 只看该作者 Only view this author
huangtf 发表于 2016-10-11 15:54
sob老师,您好!使用scf=xqc进行几何构型优化,非得把关键词stable写到输入文件(gjf)中吗?

谢谢老师!

没必要
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

8

帖子

0

威望

51

eV
积分
59

Level 2 能力者

30#
发表于 Post on 2017-7-21 09:51:50 | 只看该作者 Only view this author
SCF=IntRep请问sob老师,这个关键词的作用?

本版积分规则 Credits rule

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

GMT+8, 2026-2-22 03:13 , Processed in 0.234533 second(s), 20 queries , Gzip On.

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