应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » Gaussian新手计算HOMO-LUMO的数值跟文献所给出数值不一样
52/1返回列表
查看: 4494  |  回复: 11
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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
回复此楼

» 猜你喜欢
1楼2013-12-17 10:49:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dreamyeye

木虫 (著名写手)
后面的机组是优化和算频率的，前面的是优化好以后算单点能的

发自小木虫Android客户端
赞一下 回复此楼
12楼2016-02-20 16:52:22
 已阅   回复此楼   关注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, ★有帮助, 谢谢！能具体说一说第二个问题用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 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭