24小时热门版块排行榜    

查看: 5457  |  回复: 11
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

349924334

捐助贵宾 (著名写手)

[求助] Gaussian新手计算HOMO-LUMO的数值跟文献所给出数值不一样 已有4人参与

本人刚接触Gaussian一周左右,使用的是Window版本G09W。高斯使用指南和FAQ也都读过2遍了。
我学习Gaussian主要是看HOMO-LUMO,现在觉得大概入门。所以在文献中找一个结构进行了计算,结构如图。
但是我运算出的结果数值上跟文献有些许偏差,HOMO-LUMO图形基本一致。
我的方法是先用Chem3D计算最小能量预优化,RMS=0.0001,计算完后另存为PDB;在G09W中Utilities--NewZMat转化为Z-矩阵后File--Modify,分别输入:
%chk=D:\G09W\test\test111.chk
%mem=1000MB
# opt freq b3lyp/6-31g(d) geom=connectivity
然后运行,大概跑了3天后Processing Completed.
我计算出来的HOMO=-5.14eV,LUMO=-1.90eV;文献中的数值是HOMO=-5.73eV,LUMO=-2.09eV。
问题是
1.文献中的方法是6-311+g(d,p)//6-31g(d),它的意思是先用6-31g(d)优化,再用6-311+g(d,p)计算频率吗?我计算的偏差是不是就是因为计算方法不同导致的?
2.如果要在G09W中建立连续任务怎么操作?可不可以进行如下的操作:转换为Z-矩阵后File--Modify,然后输入
%chk=D:\G09W\test\test111.chk
%mem=1000MB
# opt b3lyp/6-31g(d) geom=connectivity
再之后Addition Step: # freq b3lyp/6-311+g(d,p)。 运行。
不知道这种操作方式是否可行?是不是就是连续任务?在Addition Step里还需要补全% Section那一项吗?
3.我在查看OUT文件时发现有两段 Population analysis using the SCF density的计算,这里的数据不就应该是HOMO-LUMO的能量数值吗?为什么会有两段数据呢?是第一部分进行的opt优化,第二部分进行的freq计算吗?我的结果以哪部分的数据为准呢?
以上就是我的问题,请各位高手多多指点!!
Gaussian新手计算HOMO-LUMO的数值跟文献所给出数值不一样
111.png


Gaussian新手计算HOMO-LUMO的数值跟文献所给出数值不一样-1
HOMO.png


Gaussian新手计算HOMO-LUMO的数值跟文献所给出数值不一样-2
LUMO.png
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

枪下游魂

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
349924334: 金币+5, 有帮助, 谢谢!能具体说一说第二个问题用G09W的操作步骤吗?我现在不清楚我上面表述的操作是不是正确的。 2013-12-17 17:39:47
gmy1990: 金币+2 2013-12-17 20:51:40
1、我觉得可以这么理解。6-311+g(d,p)的精度会比6-31g高不少。
2、可用--Link1--来分割,%的内容是要重新输入的,并用geom=allcheck来读入优化好的结构
3、opt会算,freq最后会再检查一遍结构是否合理。我的理解,如果你opt和freq的方法完全一样,这两个结果应该没有差别。
4楼2013-12-17 16:51:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

marson

银虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
349924334: 金币+5, 有帮助, 谢谢!那如果我也按照6-311+g(d,p)//6-31g(d)这种方法计算的话结果就可能跟文献一致了吧?图形差不多一致是不是就意味着趋势是与文献相同的呢?数值相差就是精度的差别呗? 2013-12-17 14:55:24
gkf高: 金币+2 2013-12-17 15:44:31
想要跟文献中的数值一样的话一定要严格按照文献中的方法跟步骤,另外,高斯版本也要一致~不过差别应该是精度,对趋势是没有大的影响的~
思想有多远,我们就能走多远~
2楼2013-12-17 13:38:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

marson

银虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
349924334: 金币+5, ★★★很有帮助, 谢谢! 2013-12-17 17:34:20
我看过有对比不同分子的HOMO及LUMO轨道值,不同的方法可能会造成数值的不同,但分子间的大小顺序趋势是相通的~另外,想要重复文献中的数值必须严格按照文献中的方法基组及版本,如果仍对应不上就可能要考虑文献中他们数据造假的问题了~
思想有多远,我们就能走多远~
3楼2013-12-17 16:33:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

雪狼乖乖

木虫 (正式写手)


【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
349924334: 金币+5, 有帮助 2013-12-18 07:56:55
有两个问题你需要注意:
1  请注意文章中是否注明了这个分子的对称性,一般HOMO-LUMO的误差都是分子对称性不一致造成的;
2  请换一下你的头像,不然我老想点你的头像拿金币。
5楼2013-12-17 21:13:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 s3fFTmArrBt6 2026-09-13 3/150 2026-09-16 18:39 by rRvxEj3zorCy
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 5/250 2026-09-16 16:12 by i3KBRD19ERkG
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +6 LwdutQ8HoqWP 2026-09-13 6/300 2026-09-16 15:24 by i3KBRD19ERkG
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 QUjhNVAcOSff 2026-09-13 6/300 2026-09-16 14:27 by i3KBRD19ERkG
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-16 08:51 by SjQobnC04j06
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-16 07:12 by LH5NkK5Ud4bh
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-16 05:51 by LH5NkK5Ud4bh
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-16 05:27 by LH5NkK5Ud4bh
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 6F5UbRU2I5hL 2026-09-14 5/250 2026-09-16 04:51 by LH5NkK5Ud4bh
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-16 03:23 by DgNGHc3h5tPl
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +7 3n8v2C8RimXI 2026-09-13 7/350 2026-09-16 02:35 by DgNGHc3h5tPl
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-16 02:14 by DgNGHc3h5tPl
[考博] 申博发邮件 +4 Lmengk 2026-09-15 4/200 2026-09-16 00:05 by 木心想继续深造
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 20:59 by zNcEhTcH7ihq
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 20:14 by zNcEhTcH7ihq
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 0pYnUiPDfdnk 2026-09-14 3/150 2026-09-15 19:35 by CgyNCDVNhGVg
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 5/250 2026-09-15 14:58 by CgyNCDVNhGVg
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 07:50 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 07:02 by 5BDX0d0WFp7t
[考博] 西北工业大学材料学院电化学传感与催化课题组招收2027年推免研究生、博士生、博士后 +3 马方园 2026-09-10 4/200 2026-09-14 09:29 by Chem张zz
信息提示
请填处理意见