计算化学公社

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

[新手求助] 求助涉及碱金属能量问题的基组和泛函选择

[复制链接 Copy URL]

6

帖子

0

威望

17

eV
积分
23

Level 1 能力者

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 dangerK 于 2025-1-25 11:44 编辑

各位老师好,我想计算单电子转移过程中氢原子(H)和碱金属(Li、Na、K)从0价到+1的能量变化,孤立的原子到离子,可以直接用M06-2X/def2-TZVP计算吗?还是有其他更准确的基组和泛函,我算出来能量差有110kcal/mol,感觉有点不太对劲。请教各位老师

139

帖子

0

威望

826

eV
积分
965

Level 4 (黑子)

2#
发表于 Post on 2025-1-25 05:49:16 | 只看该作者 Only view this author
总电子数都不一样,用DFT做这样的计算没有意义。如果你的目的是要计算单个原子的电离能,HOMO轨道的能量可以近似,参见Koopmans定理。

6

帖子

0

威望

17

eV
积分
23

Level 1 能力者

3#
 楼主 Author| 发表于 Post on 2025-1-25 07:08:32 | 只看该作者 Only view this author
ZetaFunction 发表于 2025-1-25 05:49
总电子数都不一样,用DFT做这样的计算没有意义。如果你的目的是要计算单个原子的电离能,HOMO轨道的能量可 ...

这只是半个反应,还有一半是乙醇脱氢,我知道M06-2X/def2-TZVP算有机物的反应是可以的

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

4#
发表于 Post on 2025-1-25 11:17:03 | 只看该作者 Only view this author
别乱拼def2-TZVP
算有机类体系wB97M-V比M06-2X精度明显更好
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272http://bbs.keinsci.com/thread-536-1-1.html

直接把要算的反应式贴出来,免得造成误会。隐藏关键信息只会导致无法得到及时、有效、准确的回复,同时浪费自己和别人的时间。仔细看下文
在网上求助计算化学问题的时候必须把问题描述得详细、具体、准确、清楚
http://sobereva.com/620http://bbs.keinsci.com/thread-25787-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

6

帖子

0

威望

17

eV
积分
23

Level 1 能力者

5#
 楼主 Author| 发表于 Post on 2025-1-25 11:47:53 | 只看该作者 Only view this author
本帖最后由 dangerK 于 2025-1-25 11:49 编辑
sobereva 发表于 2025-1-25 11:17
别乱拼def2-TZVP
算有机类体系wB97M-V比M06-2X精度明显更好
简谈量子化学计算中DFT泛函的选择

好的,感谢sober老师,想请问用wB97M-V或者M06-2X配合def2-TZVP计算碱金属(Li、Na、K)以单电子转移攻击乙醇的能量变化有意义吗?这里面就涉及碱金属0价到+1价,乙醇得一个电子前后的能量变化,这里面的碱金属不知道可不可以当作有机类体系一样处理

887

帖子

3

威望

1670

eV
积分
2617

Level 5 (御坂)

傻傻的木瓜

6#
发表于 Post on 2025-1-25 12:34:26 | 只看该作者 Only view this author
dangerK 发表于 2025-1-25 11:47
好的,感谢sober老师,想请问用wB97M-V或者M06-2X配合def2-TZVP计算碱金属(Li、Na、K)以单电子转移攻击 ...

“能量变化”是电子能变还是焓变或自由能变,碱金属和乙醇是什么状态?这些前提都影响最终结果,但当前描述还是太模糊。

按化学常识,碱金属与乙醇的置换反应总式是M(s) + C2H5OH(l) → C2H5OM(in C2H5OH) + 0.5H2(g)(M代表碱金属,注意物态标注),其中涉及的溶剂合碱金属离子[M(C2H5OH)n]+显然不是“孤立的原子到离子”,可能短暂存在的溶剂化电子[e(C2H5OH)n]-也不是“乙醇得一个电子”这么简单描述。当然如果所有反应物全是气相的孤立原子分子,那结果又完全不同。
√546=23.36664289109

6

帖子

0

威望

17

eV
积分
23

Level 1 能力者

7#
 楼主 Author| 发表于 Post on 2025-1-25 12:44:20 | 只看该作者 Only view this author
Uus/pMeC6H4-/キ 发表于 2025-1-25 12:34
“能量变化”是电子能变还是焓变或自由能变,碱金属和乙醇是什么状态?这些前提都影响最终结果,但当前描 ...

感谢老师,您提到的很有道理,因为碱金属在溶液中的反应性很强,所以我只想考虑气态的反应过程,计算自由能变

887

帖子

3

威望

1670

eV
积分
2617

Level 5 (御坂)

傻傻的木瓜

8#
发表于 Post on 2025-1-25 13:02:18 | 只看该作者 Only view this author
dangerK 发表于 2025-1-25 12:44
感谢老师,您提到的很有道理,因为碱金属在溶液中的反应性很强,所以我只想考虑气态的反应过程,计算自由 ...

全是气相的话,这点小体系可以考虑下热力学组合方法,比如http://bbs.keinsci.com/thread-1132-1-1.html推荐的。注意气相物种缺少溶剂化、形成晶格等显著稳定结构的因素,最终算出来的自由能变结果与现实条件的实验结果不具有可比性。
√546=23.36664289109

6

帖子

0

威望

17

eV
积分
23

Level 1 能力者

9#
 楼主 Author| 发表于 Post on 2025-1-25 13:18:21 | 只看该作者 Only view this author
Uus/pMeC6H4-/キ 发表于 2025-1-25 13:02
全是气相的话,这点小体系可以考虑下热力学组合方法,比如http://bbs.keinsci.com/thread-1132-1-1.html ...

好的,感谢老师,再请教一下,如果我想知道Li0单电子转移之后变成Li+的自由能变,就一个孤立的原子和离子,请问推荐用什么基组和泛函可以更准确描述呢?

887

帖子

3

威望

1670

eV
积分
2617

Level 5 (御坂)

傻傻的木瓜

10#
发表于 Post on 2025-1-26 20:38:52 | 只看该作者 Only view this author
dangerK 发表于 2025-1-25 13:18
好的,感谢老师,再请教一下,如果我想知道Li0单电子转移之后变成Li+的自由能变,就一个孤立的原子和离子 ...

自由能变有电子能(单点能)和热校正量两个部分,准确描述二者需要的计算级别(基组和理论方法)不同,特别是单点能需要比密度泛函更好的级别,其中理论方法可不止密度泛函(比如http://sobereva.com/552的Shermo举例就用了CCSD(T)这种后HF方法),所以当前问题说“基组和泛函”有点局限了。

如果计算资源足够支持,对简单小体系干脆直接全部上8楼链接里提到的热力学组合方法,自动走完整套流程算出所有数据,又能达到很好的精度又省去自己选计算级别了。热力学组合方法G4的原文的补充材料甚至都有锂钠钾电离能的测试。
√546=23.36664289109

6

帖子

0

威望

17

eV
积分
23

Level 1 能力者

11#
 楼主 Author| 发表于 Post on 2025-1-27 00:17:13 | 只看该作者 Only view this author
Uus/pMeC6H4-/キ 发表于 2025-1-26 20:38
自由能变有电子能(单点能)和热校正量两个部分,准确描述二者需要的计算级别(基组和理论方法)不同,特 ...

好的,感谢老师!

本版积分规则 Credits rule

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

GMT+8, 2025-8-14 04:25 , Processed in 0.144345 second(s), 20 queries , Gzip On.

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