计算化学公社

 找回密码 Forget password
 注册 Register
Views: 12050|回复 Reply: 11
打印 Print 上一主题 Last thread 下一主题 Next thread

[新手求助] 计算水的 homo lumo 值该用什么基组?

[复制链接 Copy URL]

11

帖子

0

威望

49

eV
积分
60

Level 2 能力者

请问下,我用:
#p opt b3lyp/6-311g(2d,2p) scrf=(smd,solvent=water)

来计算水的homo lumo,结果是:
lumo = 1.315 eV
homo = -8.361 eVgap = 9.676 eV


但是我看论坛里,大家给的建议都是用 6-31g 就可以,故我换了下基组:

#p opt b3lyp/6-31g(2d,2p) scrf=(smd,solvent=water)

结果为:
lumo = 2.402 eV
homo = -8.099 eVgap = 10.501 eV



其中 lumo 值差非常多,想问一下大家,哪种合理?谢谢!



87

帖子

0

威望

4014

eV
积分
4101

Level 6 (一方通行)

2#
发表于 Post on 2019-4-17 13:35:55 | 只看该作者 Only view this author
我也有这些疑问,支持一下

6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

3#
发表于 Post on 2019-4-17 13:56:17 | 只看该作者 Only view this author
直接用B3LYP/def2-TZVP

能用得动像样的基组用那么low的干嘛

“都是用 6-31g 就可以”这种说法完全莫名其妙,缺乏前提下这句话是无意义的

空轨道的形状、能量本来就不确定性大,受计算级别影响远大于占据轨道

另外注意HOMO-LUMO gap本来就没有精确不精确一说,因为根本就不是可观测量,没有什么参考标准。用首页google框搜索此论坛能找到大量讨论

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

11

帖子

0

威望

49

eV
积分
60

Level 2 能力者

4#
 楼主 Author| 发表于 Post on 2019-4-17 14:14:22 | 只看该作者 Only view this author
sobereva 发表于 2019-4-17 13:56
直接用B3LYP/def2-TZVP

能用得动像样的基组用那么low的干嘛

谢谢 sob 老师!

除了水,我还想算一下 CF3SO3 阴离子,以及 SO4 硫酸根离子在水溶剂里的 homo lumo,这些都是结构比较简单的,如您所说,我打算用 B3LYP/def2-TZVP 计算。

但是此外,我还想计算 Zn、CF3SO3 离子在水溶剂中的几个配位构型的 homo lumo,结构如下:


Zn·6H2O


Zn CF3SO3 · 5H2O


Zn 2CF3SO3 · 4H2O


对这几个结构来说,感觉用 B3LYP/def2-TZVP 就会比较吃力了。但如果我这里换回 6-311 G,那跟之前小分子的计算基组不一样,是不是又缺乏比较性了呢?
或者我就干脆全部用统一的
#p opt b3lyp/6-311g(2d,2p) scrf=(smd,solvent=water)

您有什么建议吗?

谢谢。

87

帖子

0

威望

4014

eV
积分
4101

Level 6 (一方通行)

5#
发表于 Post on 2019-4-17 14:45:24 | 只看该作者 Only view this author
是不是可以先用6-311G(2d,2p)优化结构,在优化后的结构上再用def2-TZVP 计算HOMO, LUMO

4289

帖子

4

威望

9550

eV
积分
13919

Level 6 (一方通行)

MOKIT开发者

6#
发表于 Post on 2019-4-17 14:49:01 | 只看该作者 Only view this author
waitingseven 发表于 2019-4-17 14:45
是不是可以先用6-311G(2d,2p)优化结构,在优化后的结构上再用def2-TZVP 计算HOMO, LUMO

没有必要这么做。def2-TZVP与6-311G(2d,2p)计算量相当,当然用精度更好的def2-TZVP或def2-TZVPP。
自动做多参考态计算的程序MOKIT

11

帖子

0

威望

49

eV
积分
60

Level 2 能力者

7#
 楼主 Author| 发表于 Post on 2019-4-17 14:50:11 | 只看该作者 Only view this author
zjxitcc 发表于 2019-4-17 14:49
没有必要这么做。def2-TZVP与6-311G(2d,2p)计算量相当,当然用精度更好的def2-TZVP或def2-TZVPP。

谢谢!
我先试试看。

4289

帖子

4

威望

9550

eV
积分
13919

Level 6 (一方通行)

MOKIT开发者

8#
发表于 Post on 2019-4-17 14:55:12 | 只看该作者 Only view this author
jayster 发表于 2019-4-17 14:50
谢谢!
我先试试看。

你列出的三个分子并不大,def2-TZVP完全可以,计算是否吃力看你的机器配置而已。
实际做的时候,优化结构当然不是直接上def2-TZVP了。可以先用def2-SVP优化好结构,在此基础上用def2-TZVP再优化。
自动做多参考态计算的程序MOKIT

11

帖子

0

威望

49

eV
积分
60

Level 2 能力者

9#
 楼主 Author| 发表于 Post on 2019-4-17 15:06:38 | 只看该作者 Only view this author
zjxitcc 发表于 2019-4-17 14:55
你列出的三个分子并不大,def2-TZVP完全可以,计算是否吃力看你的机器配置而已。
实际做的时候,优化结 ...

好的,谢谢!

6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

10#
发表于 Post on 2019-4-17 17:37:26 | 只看该作者 Only view this author
jayster 发表于 2019-4-17 14:14
谢谢 sob 老师!

除了水,我还想算一下 CF3SO3 阴离子,以及 SO4 硫酸根离子在水溶剂里的 homo lumo, ...

写的时候不要偷懒写 6-311 G,会造成极大的误会
以后不要用6-311g(2d,2p),给氢两层p极化根本没必要

def2-TZVP优化不动的体系,用def-TZVP或者与def-TZVP尺寸相近的6-311G**。算单点时候再改成def2-TZVP
参考
谈谈量子化学中基组的选择
http://sobereva.com/336http://bbs.keinsci.com/thread-3545-1-1.html
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

11

帖子

0

威望

49

eV
积分
60

Level 2 能力者

11#
 楼主 Author| 发表于 Post on 2019-4-17 19:27:19 | 只看该作者 Only view this author
sobereva 发表于 2019-4-17 17:37
写的时候不要偷懒写 6-311 G,会造成极大的误会
以后不要用6-311g(2d,2p),给氢两层p极化根本没必要

谢谢 sob 老师,学习了!

1043

帖子

0

威望

4188

eV
积分
5231

Level 6 (一方通行)

12#
发表于 Post on 2019-4-17 20:31:29 来自手机 | 只看该作者 Only view this author
我觉得泛函比基组可能更重要

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-24 16:20 , Processed in 0.232423 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list