24小时热门版块排行榜    

查看: 5292  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 澳大利亚 Murdoch University 全奖博士招生(3个名额)地质化工冶金领域 +80 AI8RaGaPaSCR 2026-08-07 84/4200 2026-08-14 10:45 by rkvMDKApspAV
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 by 龙-樱
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +4 majunge000 2026-08-11 6/300 2026-08-13 16:50 by 四季常青藤
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +22 Lanmanbaby 2026-08-09 36/1800 2026-08-12 22:50 by sdfapple719
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +4 工厂打螺丝 2026-08-08 4/200 2026-08-12 09:18 by hansi2025
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
信息提示
请填处理意见