计算化学公社

标题: 请教计算溶剂化结构和各溶剂分子HOMO LUMO的计算级别？ [打印本页]

作者
Author: Hanna110 时间: 2026-5-2 21:22
标题: 请教计算溶剂化结构和各溶剂分子HOMO LUMO的计算级别？
打算计算各种溶剂化结构 (例如Na+-DME-PF6-)、溶剂分子、电解质盐、PF6-(阴离子) 的homo lumo level, 网上参考文献用的计算方法不一，自己尝试的不同计算方法算出来的差别也很大，目前打算采用如下方法用高斯计算homo lumo, 大家帮忙看看行不行
step 1: 对于溶剂化结构、溶剂分子、电解质盐均考虑pcm，用b3lyp-d3/6-311+g** opt freq；对于PF6-, 真空下 b3lyp-d3/6-311+g** opt freq
step 2: 用step1结构，考虑smd, 用b3lyp-d3/6-311+g** 算能量，从out文件读取homo lumo level

作者
Author: sobereva 时间: 2026-5-2 21:24
“自己尝试的不同计算方法算出来的差别也很大”这种描述毫无意义，也不说尝试了什么方法、结果各是什么，注意看下文
在网上求助计算化学问题的时候必须把问题描述得详细、具体、准确、清楚、完整
http://sobereva.com/620（http://bbs.keinsci.com/thread-25787-1-1.html）

没理由唯独PF6-做opt freq不带溶剂模型。若无特殊原因，计算级别都统一

既然单点和opt freq的级别一样（至于IEFPCM换成SMD莫名其妙，毫无意义，当前关心的完全不是体系总能量，SMD没任何额外好处），根本毫无必要单独做单点，直接从opt freq任务输出文件里读取轨道能量计算HOMO-LUMO gap就完事了，或者用Multiwfn载入opt freq任务的fch文件后进入主功能0从文本窗口直接读取HOMO-LUMO gap，下文说了
使用Multiwfn观看分子轨道
http://sobereva.com/269（http://bbs.keinsci.com/thread-462-1-1.html）

当前计算级别没显著问题
算这些小分子带色散校正完全多余（虽然无害）

作者
Author: Hanna110 时间: 2026-5-3 11:15

sobereva 发表于 2026-5-2 21:24
“自己尝试的不同计算方法算出来的差别也很大”这种描述毫无意义，也不说尝试了什么方法、结果各是什么，注 ...

1. Na+solvent体系和Na+-solvent-无机阴离子体系用b3lyp/6-311+g**算homo lumo的时候是否用smd或者pcm会对结果产生影响吗？用smd，read是模拟了更加贴近真实的溶剂环境，我用这个的初衷是想算溶剂化结构在真实溶剂环境下的homo lumo. 本来opt也用smd但结构不收敛，所以才换成pcm
2. 之前对比过M062X/def2tzvp和b3lyp/6-311+g**,均考虑smd, 算出来的homo数值差别高达2eV, 综合272博文和文献最终还是决定采用b3lyp算homo lumo的具体数值

作者
Author: sobereva 时间: 2026-5-3 13:17

Hanna110 发表于 2026-5-3 11:15
1. Na+solvent体系和Na+-solvent-无机阴离子体系用b3lyp/6-311+g**算homo lumo的时候是否用smd或者pcm会 ...

1 看下文，弄清楚SMD和IEFPCM是什么关系、什么时候才适合用SMD。SMD的非极性部分对轨道能量、电子结构根本没有直接影响，比scrf默认的IEFPCM毫无额外优点
谈谈隐式溶剂模型下溶解自由能和体系自由能的计算
http://sobereva.com/327（http://bbs.keinsci.com/thread-3345-1-1.html）

2 仔细看下文
正确地认识分子的能隙(gap)、HOMO和LUMO
http://sobereva.com/543（http://bbs.keinsci.com/thread-16758-1-1.html）

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