利用Gaussian计算外重组能的问题

让你变成回忆 · 发表于 Post on 2017-12-13 12:40:29

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

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

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

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

jiangyan · 发表于 Post on 2025-4-14 16:23:38

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

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

铁头头铁 · 发表于 Post on 2024-3-22 10:39:04

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

感谢社长！借楼问下哈

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

My_ComputChem · 发表于 Post on 2018-8-11 09:39:14

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

多谢Sob老师耐心指点！

sobereva · 发表于 Post on 2018-8-11 07:15:05

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

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

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

My_ComputChem · 发表于 Post on 2018-8-10 14:50:57

本帖最后由 My_ComputChem 于 2018-8-10 14:54 编辑

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

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

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

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

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

wty0047 · 发表于 Post on 2018-6-28 22:09:15

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

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

让你变成回忆 · 发表于 Post on 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差距并不大。

lp2014 · 发表于 Post on 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成功计算了外重组能吗？

让你变成回忆 · 发表于 Post on 2018-3-13 19:59:55

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

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

lp2014 · 发表于 Post on 2018-3-13 19:09:41

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

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

lp2014 · 发表于 Post on 2018-3-13 18:56:53

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

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

让你变成回忆 · 发表于 Post on 2017-12-13 14:14:41

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

原来如此，谢谢老师~

sobereva · 发表于 Post on 2017-12-13 13:47:59

1 别读chk文件就完了
2 当前算的只是外重组能，opt是对溶质而言的

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

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

浏览过的版块