| 查看: 5372 | 回复: 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计算吗？我的结果以哪部分的数据为准呢？ 以上就是我的问题，请各位高手多多指点！！  111.png  HOMO.png  LUMO.png |
回复此楼» 猜你喜欢
中科院西安光机所——促进学术交流，推动产学研融合，分享领域成果，加强校企合作。
已经有0人回复
-
分子动力学 Plumed教程求助
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有119人回复
-
研究生学术新锐奖｜宁波微纳材料研讨会参会通知
已经有2人回复
-
申请2027年材料/化学类博士
已经有6人回复
-
求相关合成条件
已经有12人回复
-
2027春季申博自荐
已经有6人回复
-
青基的材料青椒们，可以把硕士联培生送我这里，要研二男生
已经有63人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请问晶体的理论计算，HOMO and LUMO是什么意思？
已经有4人回复
- gaussian view5.0 算不了HOMO,LUMO轨道？？
已经有16人回复
- 关于HOMO-LUMO gap的疑问？
已经有25人回复
- 六角密排非基面广义层错能计算
已经有12人回复
- Gaussian计算得到的前线分子轨道图能说明什么？
已经有7人回复
- 求助gaussian计算的ECD如何作图
已经有8人回复
- 循环伏安曲线基础问题--氧化峰，还原峰，氧化电位，还原电位及计算HOMO，LUMO的方法
已经有6人回复
- PWscf计算homo lumo或者valence conduction band，问题已解决
已经有14人回复
- 基态的HOMO-LUMO能隙与吸收光谱规律不一致！有这种情况吗？
已经有12人回复
- 【请教】如何根据紫外可见吸收光谱计算半导体带宽及Homo/Lumo
已经有6人回复
- 【求助】用gaussian view 做HOMO,LUMO图
已经有4人回复
- 【求助】Gaussian计算HOMO-LUMO问题
已经有17人回复
- 【求助】HOMO和LUMO值是读取绝对值吗？
已经有4人回复
雪狼乖乖 
木虫 (正式写手)
- 应助: 18 (小学生)
- 贵宾: 0.005
- 金币: 2805.8
- 散金: 567
- 红花: 34
- 帖子: 901
- 在线: 402.2小时
- 虫号: 532524
- 注册: 2008-03-25
- 性别: GG
- 专业: 理论和计算化学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
349924334: 金币+10, ★★★很有帮助, 很有帮助！我刚刚接触量化，以前关于这方面的基础了解的也相对较少。最近看了很多关于这方面的基础，感觉看得越多有些地方就越明白。你的回答也对我有很大的提示，十分感谢！ 2013-12-18 22:00:45
349924334: 金币+10, ★★★很有帮助, 很有帮助！我刚刚接触量化，以前关于这方面的基础了解的也相对较少。最近看了很多关于这方面的基础，感觉看得越多有些地方就越明白。你的回答也对我有很大的提示，十分感谢！ 2013-12-18 22:00:45
|
这三点都有问题。 1点 6-311G+虽然基组更大了，但显然它由于增加了GTO轨道，因此轨道细分更精细了。初看上去好像应该更精确，但事实上，当年6-311g的出现，主要是针对d轨道对s,p轨道的极化影响而产生的，也就是说，对于第三周期和第三周期之后的元素，才比较有效。所以单纯说精度比6-31g是完全不准确的，对于第二周期的元素，6-31g反而更好。 2点 --link1-- 的使用是有技巧的。有三点要注意： 1 chk文件必须在第一步写清楚，因为link1是要读取它的轨道信息的,当然，如果你的link1算的是新的分子，就必须重新写chk文件了； 2 所有的link1 最后一个link 必须要写 %nosave 否则会在最后一个link1计算报错 3 在out文件里，会将所有的计算结果写在一个out文件里，如果不方便，需要其他软件辅助进行切割。 --link1--的计算模板在高斯的test目录里，请用相关搜索文件搜索关键词 link1 3点 opt和freq的结果是有差别的。opt虽然每个点都是单点能，但是如果你仔细比较一下，单纯计算energy 还是有一些出入的，量化计算最好不要想当然，自己算一下才清楚。 |
6楼2013-12-18 21:09:01
雪狼乖乖
木虫 (正式写手)
- 应助: 18 (小学生)
- 贵宾: 0.005
- 金币: 2805.8
- 散金: 567
- 红花: 34
- 帖子: 901
- 在线: 402.2小时
- 虫号: 532524
- 注册: 2008-03-25
- 性别: GG
- 专业: 理论和计算化学
5楼2013-12-17 21:13:56
2楼2013-12-17 13:38:34
3楼2013-12-17 16:33:03
枪下游魂
木虫 (著名写手)
- 应助: 718 (博后)
- 金币: 8336.9
- 散金: 175
- 红花: 80
- 帖子: 1498
- 在线: 542.6小时
- 虫号: 2620258
- 注册: 2013-08-28
- 专业: 理论和计算化学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
349924334: 金币+5, ★有帮助, 谢谢！能具体说一说第二个问题用G09W的操作步骤吗？我现在不清楚我上面表述的操作是不是正确的。 2013-12-17 17:39:47
gmy1990: 金币+2 2013-12-17 20:51:40
感谢参与，应助指数 +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
枪下游魂
木虫 (著名写手)
- 应助: 718 (博后)
- 金币: 8336.9
- 散金: 175
- 红花: 80
- 帖子: 1498
- 在线: 542.6小时
- 虫号: 2620258
- 注册: 2013-08-28
- 专业: 理论和计算化学
7楼2013-12-19 06:41:31
8楼2013-12-19 09:24:21
 |
9楼2013-12-23 13:40:33
dreamyeye
木虫 (著名写手)
- 应助: 39 (小学生)
- 金币: 2029.6
- 散金: 209
- 红花: 7
- 帖子: 1049
- 在线: 357.9小时
- 虫号: 577095
- 注册: 2008-06-22
- 性别: GG
- 专业: 理论和计算化学
10楼2013-12-24 17:17:27
