计算化学公社

标题: 利用Gaussian计算外重组能的问题 [打印本页]

作者
Author: 让你变成回忆 时间: 2017-12-13 12:40
标题: 利用Gaussian计算外重组能的问题
sob老师您好，按照您的帖子《乱谈外重组能的计算》一文中，我按照你给的四个步骤计算溶液下的外重组能遇到如下问题：

（1）在第（2）步的时候，通过去读第一步中的chk文件，但最后提示错误：Inconsistent number of charges on the checkpoint file. 然后任务终止在scf的第一步；

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

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

作者
Author: sobereva 时间: 2017-12-13 13:47
1 别读chk文件就完了
2 当前算的只是外重组能，opt是对溶质而言的

作者
Author: 让你变成回忆 时间: 2017-12-13 14:14

sobereva 发表于 2017-12-13 13:47
1 别读chk文件就完了
2 当前算的只是外重组能，opt是对溶质而言的

原来如此，谢谢老师~

作者
Author: lp2014 时间: 2018-3-13 18:56

让你变成回忆发表于 2017-12-13 14:14
原来如此，谢谢老师~

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

作者
Author: lp2014 时间: 2018-3-13 19:09

让你变成回忆发表于 2017-12-13 14:14
原来如此，谢谢老师~

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

作者
Author: 让你变成回忆 时间: 2018-3-13 19:59

lp2014 发表于 2018-3-13 18:56
你好，想请教一下你现在可以正常计算吗？不读chk文件意思是所有四步都是基于log文件的结构再计算吗？谢谢 ...

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

作者
Author: lp2014 时间: 2018-3-13 20:37

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

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

作者
Author: 让你变成回忆 时间: 2018-3-13 21:09

lp2014 发表于 2018-3-13 20:37
我记得之前G09好像能正常算，但现在用G16就会出现‘Inconsistent number of charges on the checkpoint f ...

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

作者
Author: wty0047 时间: 2018-6-28 22:09

让你变成回忆发表于 2018-3-13 21:09
过了挺久了，都忘了。
当时我应该是改用G09了。因为不涉及激发态的频率分析，G09和G16差距并不大。

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

作者
Author: My_ComputChem 时间: 2018-8-10 14:50
本帖最后由 My_ComputChem 于 2018-8-10 14:54 编辑

sobereva 发表于 2017-12-13 13:47
1 别读chk文件就完了
2 当前算的只是外重组能，opt是对溶质而言的

Sob老师，见此帖有如下两个问题：

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

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

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

作者
Author: sobereva 时间: 2018-8-11 07:15

My_ComputChem 发表于 2018-8-10 14:50
Sob老师，见此帖有如下两个问题：

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

1 没法考虑
2 本来用的就是隐式溶剂模型，不可能优化溶剂分子。
直接用隐式溶剂模型算基态本身就对应于溶剂快慢部分都已经弛豫，即平衡溶剂
不了解这些知识的话看
Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法
http://sobereva.com/314
里面2.5、3.6节

作者
Author: My_ComputChem 时间: 2018-8-11 09:39

sobereva 发表于 2018-8-11 07:15
1 没法考虑
2 本来用的就是隐式溶剂模型，不可能优化溶剂分子。
直接用隐式溶剂模型算基态本身就对应于 ...

多谢Sob老师耐心指点！

作者
Author: 铁头头铁 时间: 2024-3-22 10:39

sobereva 发表于 2017-12-13 13:47
1 别读chk文件就完了
2 当前算的只是外重组能，opt是对溶质而言的

感谢社长！借楼问下哈

既然四个文件输入结构是一样的，那么我做单电子转移，用氧化态还是还原态结构？
还是说氧化态结构算出来的是氧化态的外重组能，还原态结构算出来的是还原态外重组能？

作者
Author: jiangyan 时间: 2025-4-14 16:23

wty0047 发表于 2018-6-28 22:09
**** 作者被禁止或删除内容自动屏蔽 ****

您好，请问您是怎么解决的，我也遇到了一i样的问题

欢迎光临计算化学公社 (http://bbs.keinsci.com/)