| 查看: 4932 | 回复: 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 |
回复此楼» 猜你喜欢
湖北文理学院材料与化工专硕还有调剂名额，欢迎大家推荐考生！
已经有1人回复
-
中南林科大化工院专硕第三批调剂
已经有2人回复
-
物理化学论文润色/翻译怎么收费?
已经有191人回复
-
26年储能、电池方向博士申请
已经有2人回复
-
大连工业大学杰青/长江团队-生物质材料-储能电池方向招收2026级博士生
已经有1人回复
-
欢迎硕博咨询
已经有26人回复
-
求审搞机会
已经有4人回复
-
昆明理工大学冶能院离子液体冶金课题组招收博士生
已经有2人回复
-
锂电池正极材料前驱体法制备正极材料
已经有1人回复
-
深圳理工大学-湖南大学招收电化学方向联培博士生1名
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请问晶体的理论计算，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人回复
枪下游魂
木虫 (著名写手)
- 应助: 718 (博后)
- 金币: 8336.9
- 散金: 175
- 红花: 80
- 帖子: 1498
- 在线: 542.6小时
- 虫号: 2620258
- 注册: 2013-08-28
- 专业: 理论和计算化学
7楼2013-12-19 06:41:31
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
5