计算化学公社

标题: 请教各位老师关于有机小分子与锂离子的结合能计算问题 [打印本页]

作者
Author: Fencros7y 时间: 2024-7-28 23:48
标题: 请教各位老师关于有机小分子与锂离子的结合能计算问题
本帖最后由 Fencros7y 于 2024-7-31 14:49 编辑

各位老师好。

最近我想尝试一下计算一系列大环分子和锂离子的结合能，以前有个师兄给我分享过一篇卟啉分子和锂离子结合能计算的文献，如附图所示。

我希望请教各位老师如下问题：

1、如果我的计算只涉及到的原子只有C、H、N、O和Li+，那么应该选取什么泛函和基组比较合适呢？文献使用的是B3LYP/6-311++G(d,p)，这种涉及阳离子的体系需要加这样的弥散函数吗？

2、关于溶剂化环境，我使用的是和文献一致的DOL/DME（二氧戊环/乙二醇二甲醚）=1:1的溶剂，文献中写的是“dielectric constant=7.1/7.2，solvent radius of 3.71/4.19 Å”。

那么根据sob老师的 http://sobereva.com/327 这篇文章里的介绍，我是否可以写scrf(solvent=generic,read)，然后定义eps=7.15，epsinf=1.9，而溶剂半径不需要理会。请问老师我的上述理解对吗？

作者
Author: sobereva 时间: 2024-7-29 03:01
1 6-311+G(2d,p)是性价比极高的选择，这点从J. Phys. Chem. A 2023, 127, 7023−7035的表1上可见。算弱相互作用能除非用>=4-zeta基组否则总应当加弥散，和是否是阳离子没关系。下文说了
谈谈弥散函数和“月份”基组
http://sobereva.com/119

2 对

作者
Author: Fencros7y 时间: 2024-7-29 16:55

sobereva 发表于 2024-7-29 03:01
1 6-311+G(2d,p)是性价比极高的选择，这点从J. Phys. Chem. A 2023, 127, 7023−7035的表1上可见。算 ...

感谢sob老师的回复！

作者
Author: 吹梦到西洲 时间: 2024-10-12 10:19
我想请问一下优化和振动分析使用B3LYP(D3/BJ)/6-311+G(2d,p), 后来您能量计算使用了什么级别呢?用M062X(D3)/ma-tzvpp是否合适?

欢迎光临计算化学公社 (http://bbs.keinsci.com/)