| 查看: 4847 | 回复: 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 |
回复此楼» 猜你喜欢
华南师范大学广州市光谱分析与功能探针重点实验室诚招青年英才和博士后岗位研究人员
已经有10人回复
-
200年来人们一直无法破解的难题： 气液临界点状态方程
已经有7人回复
-
物理化学论文润色/翻译怎么收费?
已经有183人回复
-
关于引力延迟下双星系统角动量不守恒的疑惑---求高手解答！
已经有10人回复
-
求标准卡片Cu0.4In0.4Zn1.2S2 PDF#32-0340
已经有2人回复
-
量子点电致发光器件（QLED）：蓝光量子点（QDs）和纳米氧化锌电子传输材料
已经有0人回复
-
需要一个圆偏振光诱导化学反应，请问使用圆偏振光作为光源，应该需要哪些设备搭建
已经有0人回复
-
求ENDNOTE破解软件安装包
已经有4人回复
-
河海大学 （211.双一流高校）钙钛矿光电实验室丁勇团队欢迎有志青年！
已经有3人回复
-
黑龙江省原子能研究院化学课题组招收2026级化学专业或化学相关专业硕士生
已经有0人回复
-
26博士申请
已经有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人回复
雪狼乖乖 
木虫 (正式写手)
- 应助: 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
