请选择 进入手机版 | 继续访问电脑版
第13届北京科音初级量子化学培训班将于10月5~8日于北京举办,请点击此链接查看详情。这是新人一次性正确、完整学习量子化学计算的最好、最快机会,能少走无数弯路,欢迎参加并相互转告!(已报满)

计算化学公社

 找回密码
 现在注册!
查看: 2722|回复: 16

[ORCA] M062x和B3LYP/g计算HOMO和LUMO能级差别巨大

[复制链接]

133

帖子

0

威望

1845

eV
积分
1978

Level 5 (御坂)

发表于 2019-5-25 16:13:30 | 显示全部楼层 |阅读模式
本帖最后由 yaochuang 于 2019-5-25 16:20 编辑

各位老师好,

      最近我用ORCA计算一个包含124个原子的有机分子能级时,发现使用M062x和B3LYP/g计算HOMO和LUMO相差非常大,请大家帮我分析一下是什么原因呢?
      使用 “opt tightscf BLYP D3 GCP(DFT/SVP) def2-SVP def2/J” 优化的分子结构:
       1.PNG
    使用“B3LYP/g def2-tzvp tightscf rijcosx def2/j def2-tzvp/c” 计算HOMO为-5.71 eV,LUMO为-3.73 eV;
    使用“M062x def2-tzvp tightscf rijcosx def2/j def2-tzvp/c” 计算HOMO为-6.64 eV,LUMO为-3.14 eV。
    实验值HOMO为-5.48 eV,LUMO为-3.83 eV。


   之前看到sob老师推荐使用M06-2X,更加适用用于有机分子的计算,所以用ORCA做了尝试,不知道是不是我的计算出现了问题,请大家帮我看看,谢谢!

opt.xyz

8.37 KB, 下载次数: 3

结构文件

E.log

318.11 KB, 下载次数: 3

B3LYP/g输出文件

E-m06.log

318.24 KB, 下载次数: 11

M062x输出文件

179

帖子

0

威望

1708

eV
积分
1887

Level 5 (御坂)

发表于 2019-5-25 16:27:00 | 显示全部楼层
差距不大才不正常。

看来还是老掉牙的b3lyp更好,和我最近的结果类似

2万

帖子

25

威望

3万

eV
积分
58546

管理员

公社社长

发表于 2019-5-25 16:46:36 | 显示全部楼层
“实验值HOMO为-5.48 eV,LUMO为-3.83 eV。”这种描述大错特错,仔细看http://bbs.keinsci.com/thread-13243-1-1.html

算能量等问题合适不代表算gap合适,仔细看诸如
http://bbs.keinsci.com/thread-3053-1-1.html
http://bbs.keinsci.com/thread-4304-1-1.html

HF成份20%左右的泛函通常才适合算HOMO-LUMO gap目的
HF成份50%多的泛函明显算的HOMO-LUMO gap过大


北京科音自然科学研究中心http://www.keinsci.com  致力于计算化学的发展和传播,长期开办最高水准的各种量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训,是提升计算化学研究水平的最佳选择。欢迎加入“北京科音”公众号获取培训最新消息和计算化学资讯!培训相关信息见《北京科音办的培训班FAQ》(http://bbs.keinsci.com/thread-5098-1-1.html)。
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群:思想家公社QQ群1号:18616395,2号:466017436。合计6000人。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一定会被拒绝加入。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

723

帖子

0

威望

1733

eV
积分
2456

Level 5 (御坂)

发表于 2019-5-25 17:11:21 来自手机 | 显示全部楼层
泛函自作用导致的误差。能否尝试下用HF算个单点呢?

444

帖子

0

威望

1679

eV
积分
2123

Level 5 (御坂)

发表于 2019-5-27 13:02:57 | 显示全部楼层
对比能量的绝对值买啥意义,看Gap不就得了。
提供面向计算化学的平台搭建和使用相关的咨询和支持,
wechat / bilibili  abbottcn
Configure, build, and customize the Linux cluster for your computational project.

2899

帖子

27

威望

1万

eV
积分
13567

Level 6 (一方通行)

发表于 2019-5-28 10:58:17 | 显示全部楼层
我猜你这分子是想讨论诸如吸收、发光、激发态电性能等与激发相关的东西,如此大的共轭体系单纯用一两个轨道来分析、对于比较接近的一系列结构、这种分析法儿定性错误的时候比对的多(我知道有很多做器件的人在文章里这么解释,胡乱解释,他们自己也知道错的时候多了去了,他们自己也承认分子轨道图看不懂能说明啥)。

133

帖子

0

威望

1845

eV
积分
1978

Level 5 (御坂)

 楼主| 发表于 2019-5-29 21:42:54 | 显示全部楼层
liyuanhe211 发表于 2019-5-28 10:58
我猜你这分子是想讨论诸如吸收、发光、激发态电性能等与激发相关的东西,如此大的共轭体系单纯用一两个轨道 ...

是的,我主要就是讨论这些内容,请问对于这种体系您有没有好的方法和建议呢?

133

帖子

0

威望

1845

eV
积分
1978

Level 5 (御坂)

 楼主| 发表于 2019-5-29 21:43:37 | 显示全部楼层
abin 发表于 2019-5-27 13:02
对比能量的绝对值买啥意义,看Gap不就得了。

由于需要讨论载流子传输,和能级的值还是有一定的关系。

133

帖子

0

威望

1845

eV
积分
1978

Level 5 (御坂)

 楼主| 发表于 2019-5-29 21:44:06 | 显示全部楼层
sobereva 发表于 2019-5-25 16:46
“实验值HOMO为-5.48 eV,LUMO为-3.83 eV。”这种描述大错特错,仔细看http://bbs.keinsci.com/thread-1324 ...

谢谢sob老师,学习了。

133

帖子

0

威望

1845

eV
积分
1978

Level 5 (御坂)

 楼主| 发表于 2019-5-29 21:44:33 | 显示全部楼层
granvia 发表于 2019-5-25 17:11
泛函自作用导致的误差。能否尝试下用HF算个单点呢?

谢谢建议,我下来试一下。

9

帖子

0

威望

135

eV
积分
144

Level 2 能力者

发表于 2019-6-23 11:12:38 | 显示全部楼层
实验测的HOMO值,一般是通过循环伏安法测氧化电位,然后与参比做比较得出的数值,这种方法科学嘛?

2万

帖子

25

威望

3万

eV
积分
58546

管理员

公社社长

发表于 2019-6-24 02:45:51 | 显示全部楼层
bryan-chen 发表于 2019-6-23 11:12
实验测的HOMO值,一般是通过循环伏安法测氧化电位,然后与参比做比较得出的数值,这种方法科学嘛?

实验根本测不到HOMO,仔细看此贴http://bbs.keinsci.com/thread-13243-1-1.html
你提到的那种做法都是搞实验的对理论稀里糊涂的人长期以讹传讹的说法,令无数量化外行信以为真。
北京科音自然科学研究中心http://www.keinsci.com  致力于计算化学的发展和传播,长期开办最高水准的各种量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训,是提升计算化学研究水平的最佳选择。欢迎加入“北京科音”公众号获取培训最新消息和计算化学资讯!培训相关信息见《北京科音办的培训班FAQ》(http://bbs.keinsci.com/thread-5098-1-1.html)。
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群:思想家公社QQ群1号:18616395,2号:466017436。合计6000人。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一定会被拒绝加入。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

9

帖子

0

威望

135

eV
积分
144

Level 2 能力者

发表于 2019-7-13 09:14:13 | 显示全部楼层
sobereva 发表于 2019-6-24 02:45
实验根本测不到HOMO,仔细看此贴http://bbs.keinsci.com/thread-13243-1-1.html
你提到的那种做法都是搞 ...

老师,在OLED器件制备中,又需要能级匹配,如果HOMO-LUMO不可测,可否将各个有机分子材料在同一泛函/基组下进行计算,得到HOMO-LUMO(或者HOMO-1,LUMO+1)然后比较其大小,看是否能级匹配?

2万

帖子

25

威望

3万

eV
积分
58546

管理员

公社社长

发表于 2019-7-14 01:17:43 | 显示全部楼层
bryan-chen 发表于 2019-7-13 09:14
老师,在OLED器件制备中,又需要能级匹配,如果HOMO-LUMO不可测,可否将各个有机分子材料在同一泛函/基组 ...

从实用角度上来说,或许勉强可以。
北京科音自然科学研究中心http://www.keinsci.com  致力于计算化学的发展和传播,长期开办最高水准的各种量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训,是提升计算化学研究水平的最佳选择。欢迎加入“北京科音”公众号获取培训最新消息和计算化学资讯!培训相关信息见《北京科音办的培训班FAQ》(http://bbs.keinsci.com/thread-5098-1-1.html)。
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群:思想家公社QQ群1号:18616395,2号:466017436。合计6000人。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一定会被拒绝加入。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

16

帖子

0

威望

503

eV
积分
519

Level 4 (黑子)

发表于 5 天前 | 显示全部楼层
sobereva 发表于 2019-5-25 16:46
“实验值HOMO为-5.48 eV,LUMO为-3.83 eV。”这种描述大错特错,仔细看http://bbs.keinsci.com/thread-1324 ...

sob老师,用B3LYP 优化完结构,是直接从opt的结果看HOMO-LUMO就行,还是需要再用HF方法算一下单点,再从单点读HOMO-LUMO?
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

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

GMT+8, 2020-9-20 21:21 , Processed in 0.166694 second(s), 32 queries .

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