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

计算化学公社

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

[Gaussian/gview] 利用Gaussian计算外重组能的问题

[复制链接]

205

帖子

0

威望

1116

eV
积分
1321

Level 4 (黑子)

发表于 2017-12-13 12:40:29 | 显示全部楼层 |阅读模式
本帖最后由 让你变成回忆 于 2017-12-13 12:57 编辑

sob老师您好,按照您的帖子《谈谈外重组能的计算》一文中,我按照你给的四个步骤计算溶液下的外重组能遇到如下问题:

(1)在第(2)步的时候,通过去读第一步中的chk文件,但最后提示错误:Inconsistent number of charges on the checkpoint file. 然后任务终止在scf的第一步;

(2)在第三步中,要使溶剂快、慢部分均驰豫,为什么不是在电荷为1,自旋多重度为2的时候进行opt任务?

内重组能的计算那篇博文我follow过,能全部走通,这个外重组能计算就走不通了。 肯定老师指正错误。

1万

帖子

25

威望

1万

eV
积分
32905

管理员

公社社长

发表于 2017-12-13 13:47:59 | 显示全部楼层
1 别读chk文件就完了
2 当前算的只是外重组能,opt是对溶质而言的
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

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

205

帖子

0

威望

1116

eV
积分
1321

Level 4 (黑子)

 楼主| 发表于 2017-12-13 14:14:41 | 显示全部楼层
sobereva 发表于 2017-12-13 13:47
1 别读chk文件就完了
2 当前算的只是外重组能,opt是对溶质而言的

原来如此,谢谢老师~

65

帖子

0

威望

279

eV
积分
344

Level 3 能力者

发表于 2018-3-13 18:56:53 | 显示全部楼层

你好,想请教一下你现在可以正常计算吗?不读chk文件意思是所有四步都是基于log文件的结构再计算吗?谢谢

65

帖子

0

威望

279

eV
积分
344

Level 3 能力者

发表于 2018-3-13 19:09:41 | 显示全部楼层

sob老师,如果不读chk文件的话,1到2步的话,如何读取溶剂的慢部分呢?

205

帖子

0

威望

1116

eV
积分
1321

Level 4 (黑子)

 楼主| 发表于 2018-3-13 19:59:55 | 显示全部楼层
lp2014 发表于 2018-3-13 18:56
你好,想请教一下你现在可以正常计算吗?不读chk文件意思是所有四步都是基于log文件的结构再计算吗?谢谢 ...

这个问题我仔细研究过,好像是G09和G16在处理这个问题上有区别。按照sob老师文中的步骤,G09能正常运行,但是G16却不行。问了下好像G16有些变化。

65

帖子

0

威望

279

eV
积分
344

Level 3 能力者

发表于 2018-3-13 20:37:30 | 显示全部楼层
让你变成回忆 发表于 2018-3-13 19:59
这个问题我仔细研究过,好像是G09和G16在处理这个问题上有区别。按照sob老师文中的步骤,G09能正常运行, ...

我记得之前G09好像能正常算,但现在用G16就会出现‘Inconsistent number of charges on the checkpoint file.’这个错误。那如果第二步中分子结构不读chk文件,从log中获得,你用G16成功计算了外重组能吗?

205

帖子

0

威望

1116

eV
积分
1321

Level 4 (黑子)

 楼主| 发表于 2018-3-13 21:09:53 | 显示全部楼层
lp2014 发表于 2018-3-13 20:37
我记得之前G09好像能正常算,但现在用G16就会出现‘Inconsistent number of charges on the checkpoint f ...

过了挺久了,都忘了。  
当时我应该是改用G09了。因为不涉及激发态的频率分析,G09和G16差距并不大。

1

帖子

0

威望

33

eV
积分
34

Level 2 能力者

发表于 2018-6-28 22:09:15 | 显示全部楼层
让你变成回忆 发表于 2018-3-13 21:09
过了挺久了,都忘了。  
当时我应该是改用G09了。因为不涉及激发态的频率分析,G09和G16差距并不大。

你好同学,你算出来的外重组能有没有出现最后算出来的外重组能的值为0的呢?我是严格按照步骤来的,最后第一部和第四步的E值一样,第二步和第三步的E值一样,好忧伤。。

4

帖子

0

威望

49

eV
积分
53

Level 2 能力者

发表于 4 天前 | 显示全部楼层
本帖最后由 My_ComputChem 于 2018-8-10 14:54 编辑
sobereva 发表于 2017-12-13 13:47
1 别读chk文件就完了
2 当前算的只是外重组能,opt是对溶质而言的

Sob老师,见此帖有如下两个问题:

1. 按您给出的四步算得的外重组能应该只是电子能,如何将溶剂重排的熵变也考虑进来呢?

2. 在第三步中也只是读取第二步的chk文件进行的单点计算,并未对溶剂分子的分布进行优化,何以见得溶剂的快、慢部分皆已弛豫?

问题可能比较低级,还请您不吝赐教。多谢!



外重组能计算步骤.png

1万

帖子

25

威望

1万

eV
积分
32905

管理员

公社社长

发表于 3 天前 | 显示全部楼层
My_ComputChem 发表于 2018-8-10 14:50
Sob老师,见此帖有如下两个问题:

1. 按您给出的四步算得的外重组能应该只是电子能,如何将溶剂重排的 ...

1 没法考虑
2 本来用的就是隐式溶剂模型,不可能优化溶剂分子。
直接用隐式溶剂模型算基态本身就对应于溶剂快慢部分都已经弛豫,即平衡溶剂
不了解这些知识的话看
Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法
http://sobereva.com/314
里面2.5、3.6节
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

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

4

帖子

0

威望

49

eV
积分
53

Level 2 能力者

发表于 3 天前 | 显示全部楼层
sobereva 发表于 2018-8-11 07:15
1 没法考虑
2 本来用的就是隐式溶剂模型,不可能优化溶剂分子。
直接用隐式溶剂模型算基态本身就对应于 ...

多谢Sob老师耐心指点!
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

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

GMT+8, 2018-8-14 18:39 , Processed in 0.215521 second(s), 26 queries .

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