计算化学公社
标题: G16计算二甲基氧化呋咱生成焓遇见的问题 [打印本页]
作者Author: mapleprinciple 时间: 2022-5-27 08:48
标题: G16计算二甲基氧化呋咱生成焓遇见的问题
本帖最后由 mapleprinciple 于 2022-5-27 12:29 编辑
我按照sob老师的CBS-Q3计算苯甲醛过程来计算二甲基氧化呋咱：分子如下
(, 下载次数 Times of downloads: 19)
输入文件如下
******************************************
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
文件读取：
(, 下载次数 Times of downloads: 39)
计算数据如下：
| 计算结果 | | | |
| | | |
| | | |
| | | |
| | | |
| | | |
| |
| 3768.047KJ/mol 实验值（102KJ/mol） |
我意识到可能是氧化呋咱多重态出现问题，再一次修改输入文件：文件如下
(, 下载次数 Times of downloads: 22)
运算出现错误：
(, 下载次数 Times of downloads: 21)
我不知道整体思路里面细节哪里出现问题，这种化合物应该怎么正确计算
作者Author: sobereva 时间: 2022-5-27 09:05
往前面找更确切的报错信息
有较大概率是优化不收敛，用gview看优化轨迹判断。可以同时写opt并在里面加上帮助收敛的关键词
作者Author: mapleprinciple 时间: 2022-5-27 09:21
想补充以下 这个化合物应该是单重态的~ 到底哪里出了问题
作者Author: mapleprinciple 时间: 2022-5-27 10:43
感谢卢老师的答疑，确实是优化不收敛！主要我现在已经迷糊了，如果二甲化氧化呋咱这个化合物是单重态，结算结果和实验值对不上，肯定说明我计算的方法出现问题，如果是三重态的化合物（还未计算），想呋咱上的N-O键保留，是不是出现自旋极化的现象？又应该怎么处理呢（附上优化最后一步的图片）
作者Author: mapleprinciple 时间: 2022-5-27 10:45
这是优化不收敛，中间的N-O断了~
作者Author: ShuangfeiZhu 时间: 2022-5-27 11:14
计算数据表格中O2分子能量应该是负的
作者Author: mapleprinciple 时间: 2022-5-27 11:42
是负的 我今天早弄表的时候 少打了个负号
作者Author: ionexchangeC 时间: 2022-5-27 11:58
气态的碳原子，自旋多重度应该是3
作者Author: ShuangfeiZhu 时间: 2022-5-27 12:23
我用你的数据算出来是105.5 kJ/mol，你的计算错了
作者Author: mapleprinciple 时间: 2022-5-27 12:28
本帖最后由 mapleprinciple 于 2022-5-27 12:29 编辑
老师，能不能展示一下您的计算过程
作者Author: mapleprinciple 时间: 2022-5-27 12:53
本帖最后由 mapleprinciple 于 2022-5-28 10:20 编辑
生成焓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
是我算错了
作者Author: sobereva 时间: 2022-5-28 01:56
到达了默认优化步数上限没收敛所致。
用最后的结构当初猜重新优化，并且加上calcall
作者Author: mapleprinciple 时间: 2022-5-28 10:33
谢谢卢老师，问题解决了 之前的问题也解决了，是我粗心大意
作者Author: 虚心求教 时间: 2023-3-9 10:54
本帖最后由 虚心求教 于 2023-3-9 11:16 编辑
看了半天终于看出来是数值带错了
| 欢迎光临 计算化学公社 (http://bbs.keinsci.com/) |
Powered by Discuz! X3.3 |