请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 1346|回复: 17

[新手求助] 高斯计算的HOMO和LUMO值与文献的差异问题

[复制链接]

21

帖子

0

威望

86

eV
积分
107

Level 2 能力者

发表于 2018-1-12 14:00:50 | 显示全部楼层 |阅读模式
各位大神,请教个问题。用高斯计算有机小分子的HOMO和LUMO,泛函和基组是:rb3lyp/6-31G* ,先进行opt优化,输出fch文件,用Gaussview打开,得到对应的HOMO和LUMO是-0.19845和-0.04850,单位换算得到-5.4026eV 和 -1.3198eV。别人文献中用的相同的基组和泛函,计算的结果不一致,-5.64eV 和 -2.44eV,HOMO相差不大,LUMO相差挺大。

以下两个问题不是太明白:
(1)不明白文献中是否有别的进一步的计算和优化?
(2)opt优化后,是否需要进行NBO计算,然后根据NBO和opt的数据,进一步优化???
谢谢!!!

510

帖子

7

威望

1556

eV
积分
2206

Level 5 (御坂)

发表于 2018-1-12 18:24:43 | 显示全部楼层
1.判断你的结构和文献中是否完全一样
2.不需要

21

帖子

0

威望

86

eV
积分
107

Level 2 能力者

 楼主| 发表于 2018-1-13 10:15:49 | 显示全部楼层
ggdh 发表于 2018-1-12 18:24
1.判断你的结构和文献中是否完全一样
2.不需要

谢谢您的回复!!!结构与文献中一样,就是照着文献中画的,泛函基组也是一样的。实在不太清楚,(1)Gaussview打开fch文件中得到的相应的HOMO和LUMO值是否需要进一步处理?  
(2)是不是用Gaussview画的结构,用“刷子”刷,进行优化,有什么窍门???

172

帖子

0

威望

1763

eV
积分
1935

Level 5 (御坂)

发表于 2018-1-13 17:23:10 | 显示全部楼层
Realzlzl 发表于 2018-1-13 10:15
谢谢您的回复!!!结构与文献中一样,就是照着文献中画的,泛函基组也是一样的。实在不太清楚,(1)Gau ...

照着画也不一定优化后的构型相同,比如顺反构型等
1 不需要
2 不需要,可以先用半经验方法预优化,比如PM6

1万

帖子

25

威望

1万

eV
积分
32295

管理员

公社社长

发表于 2018-1-14 23:49:29 | 显示全部楼层
可能:
1 优化没有收敛到与文献相同的构象
2 文献里的数据是错的
3 你在处理、读取数据时犯了低级错误
4 一些计算细节在文献里没交代,你的计算条件、设定也没有恰与之相符

你的HOMO已经和文献里相差巨大了

死活弄不明白就直接联系作者。
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

21

帖子

0

威望

86

eV
积分
107

Level 2 能力者

 楼主| 发表于 2018-1-15 08:38:20 | 显示全部楼层
赵云跳槽 发表于 2018-1-13 17:23
照着画也不一定优化后的构型相同,比如顺反构型等
1 不需要
2 不需要,可以先用半经验方法预优化,比如 ...

谢谢您的回复!!!理论计算基础比较薄弱,请问PM6是什么?泛函数?

21

帖子

0

威望

86

eV
积分
107

Level 2 能力者

 楼主| 发表于 2018-1-15 08:40:44 | 显示全部楼层
sobereva 发表于 2018-1-14 23:49
可能:
1 优化没有收敛到与文献相同的构象
2 文献里的数据是错的

好的,谢谢sob老师!

1万

帖子

25

威望

1万

eV
积分
32295

管理员

公社社长

发表于 2018-1-15 10:13:54 | 显示全部楼层
Realzlzl 发表于 2018-1-15 08:38
谢谢您的回复!!!理论计算基础比较薄弱,请问PM6是什么?泛函数?

最常用的半经验方法之一
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

21

帖子

0

威望

86

eV
积分
107

Level 2 能力者

 楼主| 发表于 2018-1-15 14:32:29 | 显示全部楼层
sobereva 发表于 2018-1-15 10:13
最常用的半经验方法之一

谢谢sob老师,我再试着计算下。

21

帖子

0

威望

86

eV
积分
107

Level 2 能力者

 楼主| 发表于 2018-1-16 16:10:48 | 显示全部楼层
sobereva 发表于 2018-1-15 10:13
最常用的半经验方法之一

sob老师,PM6/6-311G** 计算显示错误,怎么回事?查询了相关的,没有找到解决办法。

1万

帖子

25

威望

1万

eV
积分
32295

管理员

公社社长

发表于 2018-1-16 17:37:20 | 显示全部楼层
Realzlzl 发表于 2018-1-16 16:10
sob老师,PM6/6-311G** 计算显示错误,怎么回事?查询了相关的,没有找到解决办法。

半经验没法指定基组
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

21

帖子

0

威望

86

eV
积分
107

Level 2 能力者

 楼主| 发表于 2018-1-16 17:42:08 | 显示全部楼层
sobereva 发表于 2018-1-16 17:37
半经验没法指定基组

SOB老师,试验了单独的PM6和PM6/6-311G**,显示错误,刚查到有人b3lyp/6-31g(d):PM6可以计算,不确定这种是否合适???

1万

帖子

25

威望

1万

eV
积分
32295

管理员

公社社长

发表于 2018-1-16 17:44:02 | 显示全部楼层
Realzlzl 发表于 2018-1-16 17:42
SOB老师,试验了单独的PM6和PM6/6-311G**,显示错误,刚查到有人b3lyp/6-31g(d):PM6可以计算,不确定这种 ...

单独写PM6不可能有错误,除非你用的G03,或者其它关键词写错了

显然不能直接写b3lyp/6-31g(d):PM6,如果你要用PM6优化,B3LYP/6-31G*单点,要写成b3lyp/6-31g(d)//PM6
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

21

帖子

0

威望

86

eV
积分
107

Level 2 能力者

 楼主| 发表于 2018-1-16 19:51:32 | 显示全部楼层
sobereva 发表于 2018-1-16 17:44
单独写PM6不可能有错误,除非你用的G03,或者其它关键词写错了

显然不能直接写b3lyp/6-31g(d):PM6,如 ...

谢谢sob老师,可能是我写错了,我再试试。

21

帖子

0

威望

86

eV
积分
107

Level 2 能力者

 楼主| 发表于 2018-1-16 20:27:30 | 显示全部楼层
Realzlzl 发表于 2018-1-16 19:51
谢谢sob老师,可能是我写错了,我再试试。

sob老师,单独用PM6预优化,45秒就算完了,有的分子结构更快,这样正常吗?
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2018-7-18 03:15 , Processed in 0.251404 second(s), 23 queries .

快速回复 返回顶部 返回列表