G16计算二甲基氧化呋咱生成焓遇见的问题

mapleprinciple · 发表于 Post on 2022-5-27 08:48:34

本帖最后由 mapleprinciple 于 2022-5-27 12:29 编辑

我按照sob老师的CBS-Q3计算苯甲醛过程来计算二甲基氧化呋咱：分子如下

输入文件如下

******************************************

Gaussian16: EM64W-G16RevA.03 17-Dec-2016

26-May-2022

******************************************

%nprocshared=1

Will use upto 1 processors via shared memory.

%mem=4GB

%chk=F:\A\iso-reaction\CBS.chk

------------------------------

#geom=connectivity #p cbs-qb3

------------------------------

-------------------

Title CardRequired

-------------------

SymbolicZ-matrix:

Charge = 0 Multiplicity = 1

文件读取：

计算数据如下：

计算结果	电荷数	自旋多重度	能量（Hartree）
二甲基氧化呋咱	0	1	-415.226006
C（g）	0	3	-37.783017
N2（g）	0	1	-109.395156
O2（g）	0	3	-150.161299
H2（g）	0	1	-1.162778
最终生成焓	0.345857（hartree）
换算+石墨的升发焓	3768.047KJ/mol 实验值（102KJ/mol）

我意识到可能是氧化呋咱多重态出现问题，再一次修改输入文件：文件如下

运算出现错误：

我不知道整体思路里面细节哪里出现问题，这种化合物应该怎么正确计算

sobereva · 发表于 Post on 2022-5-27 09:05:12

往前面找更确切的报错信息
有较大概率是优化不收敛，用gview看优化轨迹判断。可以同时写opt并在里面加上帮助收敛的关键词

mapleprinciple · 发表于 Post on 2022-5-27 09:21:51

想补充以下这个化合物应该是单重态的~ 到底哪里出了问题

mapleprinciple · 发表于 Post on 2022-5-27 10:43:37

sobereva 发表于 2022-5-27 09:05
往前面找更确切的报错信息
有较大概率是优化不收敛，用gview看优化轨迹判断。可以同时写opt并在里面加上帮 ...

感谢卢老师的答疑，确实是优化不收敛！主要我现在已经迷糊了，如果二甲化氧化呋咱这个化合物是单重态，结算结果和实验值对不上，肯定说明我计算的方法出现问题，如果是三重态的化合物（还未计算），想呋咱上的N-O键保留，是不是出现自旋极化的现象？又应该怎么处理呢（附上优化最后一步的图片）

mapleprinciple · 发表于 Post on 2022-5-27 10:45:45

这是优化不收敛，中间的N-O断了~

ShuangfeiZhu · 发表于 Post on 2022-5-27 11:14:58

计算数据表格中O2分子能量应该是负的

mapleprinciple · 发表于 Post on 2022-5-27 11:42:24

ShuangfeiZhu 发表于 2022-5-27 11:14
计算数据表格中O2分子能量应该是负的

是负的我今天早弄表的时候少打了个负号

ionexchangeC · 发表于 Post on 2022-5-27 11:58:20

气态的碳原子，自旋多重度应该是3

ShuangfeiZhu · 发表于 Post on 2022-5-27 12:23:46

mapleprinciple 发表于 2022-5-27 11:42
是负的我今天早弄表的时候少打了个负号

我用你的数据算出来是105.5 kJ/mol，你的计算错了

mapleprinciple · 发表于 Post on 2022-5-27 12:28:29

本帖最后由 mapleprinciple 于 2022-5-27 12:29 编辑

ShuangfeiZhu 发表于 2022-5-27 12:23
我用你的数据算出来是105.5 kJ/mol，你的计算错了

老师，能不能展示一下您的计算过程

mapleprinciple · 发表于 Post on 2022-5-27 12:53:42

本帖最后由 mapleprinciple 于 2022-5-28 10:20 编辑

mapleprinciple 发表于 2022-5-27 12:28
老师，能不能展示一下您的计算过程

生成焓H=E（呋咱）-4*E（C(g)）-E（N2（g））-3*E（H2（g））-E（O2（g））
=-415.226006-4*（-37.783017）-（-109.395156）-3*（-1.62778）-（-150.161299）
=0.345857（hartree）
转换：0.345857*627.5095 *4.184+4*H（石墨蒸发焓）=3768.047466KJ/MOL

是我算错了

sobereva · 发表于 Post on 2022-5-28 01:56:31

mapleprinciple 发表于 2022-5-27 10:45
这是优化不收敛，中间的N-O断了~

到达了默认优化步数上限没收敛所致。
用最后的结构当初猜重新优化，并且加上calcall

mapleprinciple · 发表于 Post on 2022-5-28 10:33:36

sobereva 发表于 2022-5-28 01:56
到达了默认优化步数上限没收敛所致。
用最后的结构当初猜重新优化，并且加上calcall

谢谢卢老师，问题解决了之前的问题也解决了，是我粗心大意

虚心求教 · 发表于 Post on 2023-3-9 10:54:26

本帖最后由虚心求教于 2023-3-9 11:16 编辑

看了半天终于看出来是数值带错了

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

[新手求助] G16计算二甲基氧化呋咱生成焓遇见的问题

浏览过的版块