计算化学公社

标题: 用不同基组算单点能时个别数据能量差距较大 [打印本页]

作者
Author: Jiangxue 时间: 2024-7-1 16:11
标题: 用不同基组算单点能时个别数据能量差距较大
用ub3lyp/6-31g(d,p)优化，用ub3lyp/6-311++g(d,p)计算单点能，后又提高精度换成了ub3lyp def2tzvp计算单点能。发现换了基组以后计算的能量大部分没有差别，但个别数据出现了14kcal/mol的能量差，后又用ub3lyp def2qzvpp对其进行计算，发现结果与ub3lyp/6-311++g(d,p)计算的能量一致，请问这种情况该怎么办呢？

作者
Author: sobereva 时间: 2024-7-1 17:58
不说具体算的东西问这个毫无意义
倘若是必须带弥散的情况，更贵的def2-TZVP结果很可能还不如6-311++g(d,p)
仔细看
谈谈弥散函数和“月份”基组
http://sobereva.com/119

作者
Author: wzkchem5 时间: 2024-7-1 19:50
检查是不是波函数收敛到了错误的地方

作者
Author: Jiangxue 时间: 2024-7-1 20:31
本帖最后由 Jiangxue 于 2024-7-1 20:36 编辑

sobereva 发表于 2024-7-1 17:58
不说具体算的东西问这个毫无意义
倘若是必须带弥散的情况，更贵的def2-TZVP结果很可能还不如6-311++g(d,p) ...

计算的是一个铁卟啉环和药物分子的加成，输入和输出文件如下，麻烦老师帮我看看
C:\Users\Jiangxue\Desktop

作者
Author: Jiangxue 时间: 2024-7-1 20:36
本帖最后由 Jiangxue 于 2024-7-1 20:49 编辑

Jiangxue 发表于 2024-7-1 20:31
计算的是一个铁卟啉环和药物分子的加成，输入和输出文件如下，麻烦老师帮我看看

输入文件输出文件[/ur ...

(, 下载次数 Times of downloads: 16)

(, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 10)
(, 下载次数 Times of downloads: 2) (, 下载次数 Times of downloads: 0)
(, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 0)

作者
Author: sobereva 时间: 2024-7-2 19:33

Jiangxue 发表于 2024-7-1 20:36

光从gjf文件上没看出有问题
普通泛函结合def2-TZVP算这类体系的电子能量完全够用，没必要考虑其它的

作者
Author: Jiangxue 时间: 2024-7-3 16:54

wzkchem5 发表于 2024-7-1 19:50
检查是不是波函数收敛到了错误的地方

请问怎么检查是不是波函数收敛到了错误的地方呢？是检查计算的opt输出结果是否收敛吗？我的计算结果opt输出能量曲线图降低到了最低点，收敛结果有3个yes，用6-311++g(d,p)计算的结果和def2qzvpp保持一致，但用def2tzvp计算出的结果差了很多，反复检查了输入文件和初始结构，都没发现问题，不知道该如何解决这一问题

作者
Author: wzkchem5 时间: 2024-7-3 20:48

Jiangxue 发表于 2024-7-3 09:54
请问怎么检查是不是波函数收敛到了错误的地方呢？是检查计算的opt输出结果是否收敛吗？我的计算结果opt输 ...

不是，是检查最后一帧的布居分析（电荷和自旋）以及前线轨道，是否符合化学常识

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