计算化学公社

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

[综合交流] 简单量子化学问题答疑专帖(旧版,已关帖)

   关闭 [复制链接 Copy URL]

146

帖子

0

威望

1329

eV
积分
1475

Level 4 (黑子)

7021#
发表于 Post on 2018-6-22 08:43:35 | 只看该作者 Only view this author
老师,我在做IRC路径的时候,关键词IRC=(maxpoints=50,calcfc,stepsize=5),最后任务没能够正常结束,产生了47帧图像,而且第一步的点就掉下去了。这样应该怎样解决呢?

IRC.jpg (85.9 KB, 下载次数 Times of downloads: 17)

IRC.jpg

IRC1.jpg (54.95 KB, 下载次数 Times of downloads: 13)

IRC1.jpg

6万

帖子

99

威望

6万

eV
积分
125190

管理员

公社社长

7022#
 楼主 Author| 发表于 Post on 2018-6-22 12:16:13 | 只看该作者 Only view this author
刘小英雄steven 发表于 2018-6-22 08:43
老师,我在做IRC路径的时候,关键词IRC=(maxpoints=50,calcfc,stepsize=5),最后任务没能够正常结束,产生 ...

看具体报错

常见问题解决方法这里说了
在Gaussian中计算IRC的方法和常见问题
http://sobereva.com/400
北京科音自然科学研究中心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

541

帖子

3

威望

6622

eV
积分
7223

Level 6 (一方通行)

7023#
发表于 Post on 2018-6-22 23:32:59 | 只看该作者 Only view this author
本帖最后由 ABetaCarw 于 2018-6-23 11:58 编辑

老师好。我最近实验发现 石墨烯对某种农药有较大的吸附能力 但是对另一种农药吸附效果很差 甚至可以判断为没有吸附效果。吸附染料时候有相同的情况出现 对某种染料吸附很好 对别的不好。这个可以通过什么样的计算去模拟?算弱相互作用 对于说明吸附牢固程度有意义,对于吸附容量 意义不大吧。
恍惚月余,深谙人与人之间的差距。以后还应努力学习,才能与强者比肩。

230

帖子

0

威望

790

eV
积分
1020

Level 4 (黑子)

7024#
发表于 Post on 2018-6-23 10:02:16 | 只看该作者 Only view this author
老师您好,最近在计算nbo和单点的时候经常出现scf不收敛问题,想问您一下:

1)我算nbo时出现scf不收敛问题,用了scf=xqc就一直在做L508的计算,十分耗时。感觉有可能与加的弥散函数有关,想问您一下,可不可以用您帖子里说的novaracc,int=acc2e=12, conver=6 等关键词,会不会对最后结果产生影响呢?

2)计算单点能量的时候也是因为加了弥散函数所以不好收敛,已经加了novaracc,int=acc2e=12, conver=6 等关键词,效果仍然不好。能不能使用conver=5 或者conver=4 呢,不知道这样的结果能否使用呢?谢谢老师。

6万

帖子

99

威望

6万

eV
积分
125190

管理员

公社社长

7025#
 楼主 Author| 发表于 Post on 2018-6-23 12:37:05 | 只看该作者 Only view this author
Jameslzy 发表于 2018-6-23 10:02
老师您好,最近在计算nbo和单点的时候经常出现scf不收敛问题,想问您一下:

1)我算nbo时出现scf不收敛 ...

当前问题和你做不做NBO没有丝毫关系

1 会
如博文所示,绝对不要优先考虑xqc

2
解决SCF不收敛问题的方法
http://sobereva.com/61
里面那么多方法,干嘛优先想着降低收敛限?降到4、5没法用
北京科音自然科学研究中心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

6万

帖子

99

威望

6万

eV
积分
125190

管理员

公社社长

7026#
 楼主 Author| 发表于 Post on 2018-6-23 12:38:20 | 只看该作者 Only view this author
ABetaCarw 发表于 2018-6-22 23:32
老师好。我最近实验发现 石墨烯对某种农药有较大的吸附能力 但是对另一种农药吸附效果很差 甚至可以判断为 ...

吸附强度可以算相互作用能
如果牵扯到吸附容量,那就得做分子动力学模拟了
北京科音自然科学研究中心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

威望

814

eV
积分
892

Level 4 (黑子)

7027#
发表于 Post on 2018-6-23 15:14:14 | 只看该作者 Only view this author
Sob老师好。1.  如果计算基态双自由基体系已经考虑了对称性破缺,那么用CCSD(T)后HF方法计算单点能的时候,可以直接读取之前优化时自旋极化单重态结构的.chk文件算单点能吧?CCSD(T)算单点的时候,参考态波函数还要用对称破缺的方式再用一次U结合guess=mix再给出stable测试判断是否是稳定的波函数吗?
2. 如果一个闭壳层分子与双自由基发生加成反应,生成的TS或者中间体是不是也应当考虑自旋极化单重态呢?
谢谢Sob老师啦,周末愉快!
自己决定自己的DNA

230

帖子

0

威望

790

eV
积分
1020

Level 4 (黑子)

7028#
发表于 Post on 2018-6-23 15:15:06 | 只看该作者 Only view this author
sobereva 发表于 2018-6-23 12:37
当前问题和你做不做NBO没有丝毫关系

1 会

谢谢老师解答,那我再试试博文里的其他方法。我还想问一下,就是计算nbo时如果一个关键词用scf(novaracc,noincfock)int=acc2e=12,另一个不用,这两个结果会有比较大的差别吗,如果有的话我就把结构再重新用关键词刷一下。

239

帖子

0

威望

2333

eV
积分
2572

Level 5 (御坂)

7029#
发表于 Post on 2018-6-23 15:59:39 | 只看该作者 Only view this author
老师好,我现在在做一个BDE的计算
1.在分别优化后的片段上做振动分析freq,计算得到Thermal correction to Enthalpy,高斯B3LYP-D3BJ/6-31g**
2.在ORCA中计算单点,得到Total Energy,PWPB95-D3BJ/def2-TZVP
这样两步操作得到焓值Enthalpy(Thermal correction to Enthalpy+Total Energy),BDE=Enthalpy(A+B-AB),以上计算思路没有问题吧
还有就是在ORCA中,在当前计算条件下,自动进行了MP2的计算,这没有问题吧
麻烦老师了

38

帖子

0

威望

1553

eV
积分
1591

Level 5 (御坂)

7030#
发表于 Post on 2018-6-24 12:15:43 | 只看该作者 Only view this author
本帖最后由 xiaojingling016 于 2018-6-24 12:20 编辑
dreamyeye 发表于 2018-6-14 19:37
碘离子就不需要优化,算个单点能就行,可以把出错的输入和输出文件贴上来看看。

其实我想用的是它的SCF能量,输入文件和输出如下:

QQ图片20180624121745.png (5.83 KB, 下载次数 Times of downloads: 8)

输入

输入

s0.log

7.26 KB, 下载次数 Times of downloads: 0

输出

6万

帖子

99

威望

6万

eV
积分
125190

管理员

公社社长

7031#
 楼主 Author| 发表于 Post on 2018-6-24 12:25:47 | 只看该作者 Only view this author
柒月小鱼 发表于 2018-6-23 15:59
老师好,我现在在做一个BDE的计算
1.在分别优化后的片段上做振动分析freq,计算得到Thermal correction to  ...

没问题
双杂化泛函本身就有MP2这一步,肯定得算,只不过此时的MP2部分是基于普通泛函的轨道做的
北京科音自然科学研究中心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

6万

帖子

99

威望

6万

eV
积分
125190

管理员

公社社长

7032#
 楼主 Author| 发表于 Post on 2018-6-24 12:29:01 | 只看该作者 Only view this author
Jameslzy 发表于 2018-6-23 15:15
谢谢老师解答,那我再试试博文里的其他方法。我还想问一下,就是计算nbo时如果一个关键词用scf(novaracc ...

前两个不影响结果,最后一个仅对大体系有略微影响,一般可以忽略
北京科音自然科学研究中心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

6万

帖子

99

威望

6万

eV
积分
125190

管理员

公社社长

7033#
 楼主 Author| 发表于 Post on 2018-6-24 12:55:56 | 只看该作者 Only view this author
QuantumHan 发表于 2018-6-23 15:14
Sob老师好。1.  如果计算基态双自由基体系已经考虑了对称性破缺,那么用CCSD(T)后HF方法计算单点能的时候, ...

1 如果chk里已经是对称破缺波函数了,算CCSD(T)的时候guess=mix就完了,stable和mix都不写

2 应当考虑。即便不是对称破缺态,这么算完结果自然而然也是闭壳层的结果,不会有错。
北京科音自然科学研究中心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

239

帖子

0

威望

2333

eV
积分
2572

Level 5 (御坂)

7034#
发表于 Post on 2018-6-24 14:57:44 | 只看该作者 Only view this author
sobereva 发表于 2018-6-24 12:25
没问题
双杂化泛函本身就有MP2这一步,肯定得算,只不过此时的MP2部分是基于普通泛函的轨道做的

好的 谢谢老师
您说了就有谱了

230

帖子

0

威望

790

eV
积分
1020

Level 4 (黑子)

7035#
发表于 Post on 2018-6-24 15:21:16 | 只看该作者 Only view this author
sobereva 发表于 2018-6-24 12:29
前两个不影响结果,最后一个仅对大体系有略微影响,一般可以忽略

明白了,多谢老师解答

本版积分规则 Credits rule

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

GMT+8, 2026-2-26 12:55 , Processed in 0.186326 second(s), 22 queries , Gzip On.

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