用更高精度级别计算，单点能反而更高

wudiazhu · 发表于 Post on 2021-10-27 09:39:30

各位老师好，我使用B3LYP/6-31G**对一个体系做几何优化，优化完成最后一帧能量HF=-2441.709Ha，已收敛且无虚频。取最后一帧用M062X/6-311G**级别做单点，最后读取的能量HF=-2441.100Ha，不太理解为什么更精确的计算级别算得的能量反而更高，这样是合理的吗。下面是我做几何优化的输入和输出文件，做单点的输入和输出文件。谢谢各位老师！

PDA_Cr_m062x_lanl2dz_14_IEFPCM_z_sp.gjf (11.57 KB, 下载次数 Times of downloads: 32)

PDA_Cr_m062x_lanl2dz_14_IEFPCM_z_sp.out (99.74 KB, 下载次数 Times of downloads: 7)

PDA_Cr_b3lyp_lanl2dz_13_IEFPCM_z_loose.gjf (11.62 KB, 下载次数 Times of downloads: 12)

PDA_Cr_b3lyp_lanl2dz_13_IEFPCM_z_loose.zip (1.55 MB, 下载次数 Times of downloads: 4)

snljty · 发表于 Post on 2021-10-27 09:56:04

没有任何直接关系，又不是同一方法下大基组和小基组对比。
顺便，LanL2DZ实在太小，给优化都很勉强，所谓“高精度单点”完全不能接受。看卢老师博文《谈谈赝势基组的选用》http://bbs.keinsci.com/thread-5625-1-1.html如果你要算“高精度单点”，基组统一用def2TZVP即可，算不动至少也用def-TZVP。

Accelerator · 发表于 Post on 2021-10-27 10:08:25

不要给过渡金属体系用M06-2X，喜欢明尼苏达至少也要用M06（虽然也不是好选择）；MN15相比之下好得多
优化B3LYP 单点换成别的杂化泛函也不合适，如果认为另一个泛函精度更高，为什么优化时不用
直接全程PBE0-D3BJ就好了

sobereva · 发表于 Post on 2021-10-27 12:28:15

影响能量的所有因素严格同一时，能量才有可比性，这是最基本常识

wudiazhu · 发表于 Post on 2021-10-27 14:31:36

snljty 发表于 2021-10-27 09:56
没有任何直接关系，又不是同一方法下大基组和小基组对比。
顺便，LanL2DZ实在太小，给优化都很勉强，所谓 ...

谢谢！多谢建议。关于赝势一开始是用的SDD，发现无论如何优化到后面总是震荡，后改成了LanL2TZ，结果又总是在提交任务的时候说读取赝势失败，估计是高斯没有内置这个赝势？遂用了LanL2DZ，知道许多老师说这个赝势已经过时了，但当时想的是我整个体系约70个原子只有1个过渡金属原子使用赝势，可能赝势的质量对计算结果影响不会很大？不知道我这样想对不对。

wudiazhu · 发表于 Post on 2021-10-27 14:38:36

Accelerator 发表于 2021-10-27 10:08
不要给过渡金属体系用M06-2X，喜欢明尼苏达至少也要用M06（虽然也不是好选择）；MN15相比之下好得多
优化B ...

谢谢老师的建议，主要是我做计算的服务器只有12个物理核心，用b3lyp/6-31G**去优化我这个70个原子左右的体系都要接近2天了，所以用更高精度的泛函耗时更高了。感谢您的建议，我去用下PBE0。

wudiazhu · 发表于 Post on 2021-10-27 14:41:08

sobereva 发表于 2021-10-27 12:28
影响能量的所有因素严格同一时，能量才有可比性，这是最基本常识

多谢sob老师提醒，我再补补基本常识。

snljty · 发表于 Post on 2021-10-27 15:06:31

wudiazhu 发表于 2021-10-27 14:31
谢谢！多谢建议。关于赝势一开始是用的SDD，发现无论如何优化到后面总是震荡，后改成了LanL2TZ，结果又总 ...

显然没有内置，有没有内置自己看手册https://gaussian.com/basissets/就知道，根本不需要试。没有内置就去EMSL拷贝定义，看卢老师博文《详解Gaussian中混合基组、自定义基组和赝势基组的输入》和李老师的这个。http://bbs.keinsci.com/thread-4385-1-1.html。http://sobereva.com/60。没有内置不是偷懒干脆不用了的理由。优化震荡和基组没有任何直接关系。我没空仔细看你结构，但是过渡金属的电子结构更难描述，本来就适合稍微大一些的基组。只有一个金属就算你换成def2-TZVP计算量都不会有明显变化，更不应该用这么小连极化函数都没有的基组。

wudiazhu · 发表于 Post on 2021-10-27 15:40:48

snljty 发表于 2021-10-27 15:06
显然没有内置，有没有内置自己看手册https://gaussian.com/basissets/就知道，根本不需要试。没有内置就 ...

好的，多谢老师！受教了！

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] 用更高精度级别计算，单点能反而更高

浏览过的版块