计算化学公社

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

[Gaussian/gview] 请教关于Linux系统下使用optDFTw进行w值的计算过程

[复制链接 Copy URL]

6万

帖子

99

威望

6万

eV
积分
125136

管理员

公社社长

16#
发表于 Post on 2017-1-6 22:44:05 | 只看该作者 Only view this author
yezhonghua 发表于 2017-1-6 21:31
sober 老师,我优化文献上一个结构,文献值w=0.183;但是我优化了14圈了,J=0.002,依然没有收敛;但是w ...

可以确定。如果反复不收敛,可以在这个小区间里再用scan方式确定
北京科音自然科学研究中心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

430

帖子

0

威望

784

eV
积分
1214

Level 4 (黑子)

17#
 楼主 Author| 发表于 Post on 2017-1-7 20:52:26 | 只看该作者 Only view this author
sobereva 发表于 2017-1-6 22:44
可以确定。如果反复不收敛,可以在这个小区间里再用scan方式确定

谢谢sober老师,主要是我的体系150左右的原子,严格按照您的收敛规则来的话,耗时比较长。我用scan跑一下试试。
另外,还有出现情况是开始7圈是收敛的,第八圈的N-1.gjf是不收敛的;程序接着就跑第9圈了;这样的情况是否应该停止,重新对原gjf文件进行SCF收敛方面进行处理。谢谢!

430

帖子

0

威望

784

eV
积分
1214

Level 4 (黑子)

18#
 楼主 Author| 发表于 Post on 2017-1-7 21:12:15 | 只看该作者 Only view this author
冰释之川 发表于 2017-1-6 22:03
由于Sob老师的程序采用的是基于Brent算法来最小化J^2的,一次次迭代有高低起伏是正常的,你等结果跑完即 ...

谢谢,我发现开始跑8圈没有问题,第九圈的N-1.gjf突然说无法收敛,然后程序跑第10圈去了,这个不影响结果么?我感觉不收敛的情况下,不值得去跑第10圈吧。结果跑了20个小时都没有结束。

6万

帖子

99

威望

6万

eV
积分
125136

管理员

公社社长

19#
发表于 Post on 2017-1-7 21:46:55 | 只看该作者 Only view this author
yezhonghua 发表于 2017-1-7 20:52
谢谢sober老师,主要是我的体系150左右的原子,严格按照您的收敛规则来的话,耗时比较长。我用scan跑一下 ...


应该停止

如果能定位到0.1828-0.1830这样0.0001左右的精度,就已经足够了,没必要继续弄下去,可以scan一下再做确认确实是这个位置附近。
北京科音自然科学研究中心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

78

帖子

0

威望

260

eV
积分
338

Level 3 能力者

20#
发表于 Post on 2017-4-1 19:51:08 | 只看该作者 Only view this author
冰释之川 发表于 2017-1-6 22:03
由于Sob老师的程序采用的是基于Brent算法来最小化J^2的,一次次迭代有高低起伏是正常的,你等结果跑完即 ...

请问一下 使用OPTDFTw计算时    w收敛线为0.0001,输出文件中看是否收敛是看J值,还是J^2那?    最后跑了16圈J=0.002,是不是代表没收敛,就在输入文件中加入scf=(tight,maxcycle=500),这样做对吗?

1102

帖子

18

威望

6703

eV
积分
8165

Level 6 (一方通行)

計算化学の社畜

21#
发表于 Post on 2017-4-1 20:03:05 | 只看该作者 Only view this author
阳光总在风雨后 发表于 2017-4-1 19:51
请问一下 使用OPTDFTw计算时    w收敛线为0.0001,输出文件中看是否收敛是看J值,还是J^2那?    最后跑 ...

看J^2是否达到最小值了。推荐用J^2指标
Stand on the shoulders of giants

6万

帖子

99

威望

6万

eV
积分
125136

管理员

公社社长

22#
发表于 Post on 2017-4-2 01:32:57 | 只看该作者 Only view this author
阳光总在风雨后 发表于 2017-4-1 19:51
请问一下 使用OPTDFTw计算时    w收敛线为0.0001,输出文件中看是否收敛是看J值,还是J^2那?    最后跑 ...


scf=(tight,maxcycle=500)毫无意义。看此文
常见的多余的和被滥用的Gaussian关键词
http://sobereva.com/331
没收敛就继续等,几十圈都收敛不了,如果收敛精度已经足够满足需要了可以直接停掉,或者尝试换不同的初猜w值
北京科音自然科学研究中心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

78

帖子

0

威望

260

eV
积分
338

Level 3 能力者

23#
发表于 Post on 2017-4-5 08:31:52 | 只看该作者 Only view this author
sobereva 发表于 2017-4-2 01:32
scf=(tight,maxcycle=500)毫无意义。看此文
常见的多余的和被滥用的Gaussian关键词
http://sobereva. ...

作业跑了14圈就结束了   您说的继续跑是什么意思那?  是这个圈数可以自己设定吗?

1102

帖子

18

威望

6703

eV
积分
8165

Level 6 (一方通行)

計算化学の社畜

24#
发表于 Post on 2017-4-5 09:05:04 | 只看该作者 Only view this author
阳光总在风雨后 发表于 2017-4-5 08:31
作业跑了14圈就结束了   您说的继续跑是什么意思那?  是这个圈数可以自己设定吗?

首先你要说清楚是 SCF不收敛,还是w优化迭代在振荡。这是两码事
Stand on the shoulders of giants

430

帖子

0

威望

784

eV
积分
1214

Level 4 (黑子)

25#
 楼主 Author| 发表于 Post on 2017-5-23 09:03:21 | 只看该作者 Only view this author
冰释之川 发表于 2017-1-5 17:02
1. w值的优化肯定是基于优化好的基态构型坐标来找的。换句话说,你要在该结构下做一系列不同w值的单点能 ...

sober老师,请教一下我在基态优化好的基础上进行w值搜索,开始的N.gjf和N-1.gjf SCF都能收敛,但是N+1.gjf无法收敛,导致计算中断。说明我的基态优化结果应该是没有问题的。
因为我是在优化好基态的基础上进行w优化的,所以我需要重新优化基态结构么?还是对已经优化好的基态进行稍微的结构调整?谢谢sober老师

1102

帖子

18

威望

6703

eV
积分
8165

Level 6 (一方通行)

計算化学の社畜

26#
发表于 Post on 2017-5-23 11:39:45 | 只看该作者 Only view this author
yezhonghua 发表于 2017-5-23 09:03
sober老师,请教一下我在基态优化好的基础上进行w值搜索,开始的N.gjf和N-1.gjf SCF都能收敛,但是N+1.gj ...

结构不变,SCF不收敛参考《解决SCF不收敛问题的方法》http://sobereva.com/61
Stand on the shoulders of giants

430

帖子

0

威望

784

eV
积分
1214

Level 4 (黑子)

27#
 楼主 Author| 发表于 Post on 2017-5-23 13:36:23 | 只看该作者 Only view this author
冰释之川 发表于 2017-5-23 11:39
结构不变,SCF不收敛参考《解决SCF不收敛问题的方法》http://sobereva.com/61

谢谢!我用的表头是#p LC-wPBE/6-31g(d)  SCF(novaracc,noincfock) int=acc2e=12
基本上算是考虑到SCF不收敛的情况了。
我现在改成#p LC-wPBE/6-31g(d) nosymm SCF(novaracc,noincfock) int=acc2e=12
如果不行的话,我可否成尝试改变一下优化后的结构?

6万

帖子

99

威望

6万

eV
积分
125136

管理员

公社社长

28#
发表于 Post on 2017-5-23 13:42:25 | 只看该作者 Only view this author
yezhonghua 发表于 2017-5-23 13:36
谢谢!我用的表头是#p LC-wPBE/6-31g(d)  SCF(novaracc,noincfock) int=acc2e=12
基本上算是考虑到SCF不 ...

博文里明确说了,这三个通常只对因为弥散函数导致的不收敛管用
绝对不要以为加了这三个就万事大吉了
北京科音自然科学研究中心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

1102

帖子

18

威望

6703

eV
积分
8165

Level 6 (一方通行)

計算化学の社畜

29#
发表于 Post on 2017-5-23 13:46:42 | 只看该作者 Only view this author
本帖最后由 冰释之川 于 2017-5-23 13:48 编辑
yezhonghua 发表于 2017-5-23 13:36
谢谢!我用的表头是#p LC-wPBE/6-31g(d)  SCF(novaracc,noincfock) int=acc2e=12
基本上算是考虑到SCF不 ...

你w优化本身就是基于优化好的几何构型来执行的,什么叫“改变一下优化后的结构”,几何结构不能再动来动去了?如果你使用了某种方法SCF不收敛,那么换其他方法试试。
Stand on the shoulders of giants

430

帖子

0

威望

784

eV
积分
1214

Level 4 (黑子)

30#
 楼主 Author| 发表于 Post on 2017-5-23 21:14:38 | 只看该作者 Only view this author
冰释之川 发表于 2017-5-23 13:46
你w优化本身就是基于优化好的几何构型来执行的,什么叫“改变一下优化后的结构”,几何结构不能再动来动 ...

谢谢大神,我按照sober老师的帖子再试试

本版积分规则 Credits rule

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

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

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