计算化学公社

标题: CP2K计算中DZVP-MOLOPT-SR-GTH基组的问题 [打印本页]

作者
Author:
CCCC    时间: 2020-7-11 20:21
标题: CP2K计算中DZVP-MOLOPT-SR-GTH基组的问题
我在CP2K的优化文件中看到了这样的基组,BASIS_SET DZVP-MOLOPT-SR-GTH,泛函是B3LYP,请问这是不是对应B3LYP/6-31G(d,p)?

作者
Author:
wth1219    时间: 2020-7-11 22:32
MOLOPT是CP2K的特制基组,一般与GTH赝势联用。DZVP基组的精度相当于6-31G*(据https://www.cp2k.org/_media/even ... ng_basis_pseudo.pdf)。
作者
Author:
CCCC    时间: 2020-7-12 16:33
wth1219 发表于 2020-7-11 22:32
MOLOPT是CP2K的特制基组,一般与GTH赝势联用。DZVP基组的精度相当于6-31G*(据https://www.cp2k.org/_media ...

谢谢回复~
作者
Author:
tuzhidingdong    时间: 2020-10-9 23:50
从EMSL Basis Set中可能会找到一些熟悉的基组,例如6-31G(d,p)等,另外可以从https://www.basissetexchange.org/找到cp2k格式的基组文件。

作者
Author:
haibeih    时间: 2020-10-10 08:25
用推荐的基组就好呀~
https://www.cp2k.org/exercises:2018_uzh_cmest:faq
作者
Author:
elpa    时间: 2020-11-29 20:23
wth1219 发表于 2020-7-11 22:32
MOLOPT是CP2K的特制基组,一般与GTH赝势联用。DZVP基组的精度相当于6-31G*(据https://www.cp2k.org/_media ...

MOLOPT是CP2K的特制基组
请问BASIS_MOLOPT适用于哪些体系?

作者
Author:
wth1219    时间: 2020-12-2 14:35
elpa 发表于 2020-11-29 20:23
MOLOPT是CP2K的特制基组
请问BASIS_MOLOPT适用于哪些体系?

您好。仅就本人对于CP2K的了解以及浏览过的文献来看,MOLOPT基组的使用没有体系限制。但孤立体系是很少使用CP2K计算的。
作者
Author:
elpa    时间: 2020-12-9 23:33
wth1219 发表于 2020-12-2 14:35
您好。仅就本人对于CP2K的了解以及浏览过的文献来看,MOLOPT基组的使用没有体系限制。但孤立体系是很少使 ...

谢谢
作者
Author:
sobereva    时间: 2021-2-15 05:19
MOLOPT与一般的量化领域的基组有极大的差异。DZVP-MOLOPT-SR-GTH与6-31(d,p)虽然表面上看都是带极化函数的2-zeta基组,但有极大不同,体现在
(1)MOLOPT系列都是和GTH赝势结合使用的赝势基组
(2)MOLOPT系列都是完全广义收缩基组,而6-31G(d,p)是片段收缩。这导致基于前者在没有特意考虑广义收缩的程序里做计算,比如在Multiwfn中做波函数分析,耗时会高非常非常多
(3)DZVP-MOLOPT-GTH是带有指数很小的高度弥散的高斯函数的,带SR(short-range)的DZVP-MOLOPT-SR-GTH虽然没有了指数那么小的高斯函数,但仍然有半弥散程度的函数,这是6-31G(d,p)不具备的。DZVP-MOLOPT-SR-GTH的最外层的GTF的指数比6-311G**的还要小。
(4)由于上述(3)的特征,DZVP-MOLOPT-SR-GTH,特别是DZVP-MOLOPT-GTH,在算需要弥散函数的偶极矩计算问题上明显好于6-31G(d,p),BSSE问题也小得多得多。非要更贴切地对比,DZVP-MOLOPT-GTH应该和6-31+G(d,p)比。

CP2K里做GPW计算基本上能用MOLOPT都建议用MOLOPT,支持的元素比较全,气相和凝聚相都适用。一般用SR版就可以。如果对偶极矩、弱相互作用问题比较重视,用更贵的不带SR的。

作者
Author:
k64_cc    时间: 2021-2-15 10:58
标题: ()
本帖最后由 k64_cc 于 2021-2-15 11:00 编辑

DZVP-MOLOPT系列是他们用小分子特别tune过的,在分子体系中达到了很高的精度,但是算起来很慢,特别慢,非常慢(原文献里也稍微提了一下这事)

SR的弥散函数没那么夸张,运算起来快多了,平均下来耗时比DZVP贵一小小点,但是精度比DZVP要好。对于凝聚相体系来说加大弥散本来就没什么必要,平常用SR就行。
作者
Author:
卡开发发    时间: 2021-2-15 14:50
本帖最后由 卡开发发 于 2021-2-15 14:51 编辑

哪怕是凝聚相或者表面,BSSE的问题还是需要小心的,可以参看(DOI: 10.1021/acs.jpcc.9b03554),需要注意的是有些结果看上去比较好看可能得益于误差抵消。
作者
Author:
nusiew    时间: 2021-7-4 18:40
由于DZVP-MOLOPT-GTH基组仅涵盖H , C , N , O , F , Si , P , S , Cl元素,而DZVP-MOLOPT-SR-GTH基组覆盖了元素周期表绝大部分元素,
请问计算同一个体系中DZVP-MOLOPT-GTH基组和DZVP-MOLOPT-SR-GTH基组可以混用吗?
作者
Author:
sobereva    时间: 2021-7-18 02:35
nusiew 发表于 2021-7-4 18:40
由于DZVP-MOLOPT-GTH基组仅涵盖H , C , N , O , F , Si , P , S , Cl元素,而DZVP-MOLOPT-SR-GTH基组覆盖了 ...

没特殊情况就都用DZVP-MOLOPT-SR-GTH
对弥散函数要求比较高的时候才适合特意花更高代价用不带SR的
作者
Author:
xxq星球    时间: 2021-7-22 19:53
sobereva 发表于 2021-7-18 02:35
没特殊情况就都用DZVP-MOLOPT-SR-GTH
对弥散函数要求比较高的时候才适合特意花更高代价用不带SR的

一个体系内DZVP-MOLOPT-SR-GTH和TZVP-MOLOPT-SR-GTH可以混用吗?
作者
Author:
丁越    时间: 2021-7-22 22:28
xxq星球 发表于 2021-7-22 19:53
一个体系内DZVP-MOLOPT-SR-GTH和TZVP-MOLOPT-SR-GTH可以混用吗?

可以,这就是混合基组

作者
Author:
xxq星球    时间: 2021-7-23 10:21
丁越 发表于 2021-7-22 22:28
可以,这就是混合基组

谢谢!
作者
Author:
zhaoxiaoni    时间: 2023-1-4 17:17
丁越 发表于 2021-7-22 22:28
可以,这就是混合基组

这两个基组在精度上差别大吗
作者
Author:
丁越    时间: 2023-1-4 18:29
zhaoxiaoni 发表于 2023-1-4 17:17
这两个基组在精度上差别大吗

一个是2-zeta基组,另一个是3-zeta基组,肯定是后者精度高,相应很耗时。




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3