计算化学公社

 找回密码 Forget password
 注册 Register
Views: 25762|回复 Reply: 17
打印 Print 上一主题 Last thread 下一主题 Next thread

[CP2K] CP2K计算中DZVP-MOLOPT-SR-GTH基组的问题

[复制链接 Copy URL]

65

帖子

0

威望

392

eV
积分
457

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
我在CP2K的优化文件中看到了这样的基组,BASIS_SET DZVP-MOLOPT-SR-GTH,泛函是B3LYP,请问这是不是对应B3LYP/6-31G(d,p)?

104

帖子

0

威望

3055

eV
积分
3159

Level 5 (御坂)

被催的人

2#
发表于 Post on 2020-7-11 22:32:31 | 只看该作者 Only view this author
MOLOPT是CP2K的特制基组,一般与GTH赝势联用。DZVP基组的精度相当于6-31G*(据https://www.cp2k.org/_media/even ... ng_basis_pseudo.pdf)。

评分 Rate

参与人数
Participants 2
eV +8 收起 理由
Reason
hebrewsnabla + 5 赞!
thanhtam + 3

查看全部评分 View all ratings

65

帖子

0

威望

392

eV
积分
457

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2020-7-12 16:33:31 | 只看该作者 Only view this author
wth1219 发表于 2020-7-11 22:32
MOLOPT是CP2K的特制基组,一般与GTH赝势联用。DZVP基组的精度相当于6-31G*(据https://www.cp2k.org/_media ...

谢谢回复~

22

帖子

0

威望

331

eV
积分
353

Level 3 能力者

4#
发表于 Post on 2020-10-9 23:50:11 | 只看该作者 Only view this author
从EMSL Basis Set中可能会找到一些熟悉的基组,例如6-31G(d,p)等,另外可以从https://www.basissetexchange.org/找到cp2k格式的基组文件。

36

帖子

0

威望

580

eV
积分
616

Level 4 (黑子)

5#
发表于 Post on 2020-10-10 08:25:32 | 只看该作者 Only view this author

80

帖子

0

威望

897

eV
积分
977

Level 4 (黑子)

6#
发表于 Post on 2020-11-29 20:23:18 | 只看该作者 Only view this author
wth1219 发表于 2020-7-11 22:32
MOLOPT是CP2K的特制基组,一般与GTH赝势联用。DZVP基组的精度相当于6-31G*(据https://www.cp2k.org/_media ...

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

104

帖子

0

威望

3055

eV
积分
3159

Level 5 (御坂)

被催的人

7#
发表于 Post on 2020-12-2 14:35:53 | 只看该作者 Only view this author
elpa 发表于 2020-11-29 20:23
MOLOPT是CP2K的特制基组
请问BASIS_MOLOPT适用于哪些体系?

您好。仅就本人对于CP2K的了解以及浏览过的文献来看,MOLOPT基组的使用没有体系限制。但孤立体系是很少使用CP2K计算的。

80

帖子

0

威望

897

eV
积分
977

Level 4 (黑子)

8#
发表于 Post on 2020-12-9 23:33:30 | 只看该作者 Only view this author
wth1219 发表于 2020-12-2 14:35
您好。仅就本人对于CP2K的了解以及浏览过的文献来看,MOLOPT基组的使用没有体系限制。但孤立体系是很少使 ...

谢谢
cp2k QMMM

5万

帖子

99

威望

5万

eV
积分
112486

管理员

公社社长

9#
发表于 Post on 2021-2-15 05:19:16 | 只看该作者 Only view this author
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的。

评分 Rate

参与人数
Participants 3
eV +12 收起 理由
Reason
挨打祥 + 2 谢谢
乐平 + 5 谢谢分享
wth1219 + 5 受教!

查看全部评分 View all ratings

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

545

帖子

0

威望

3123

eV
积分
3668

Level 5 (御坂)

10#
发表于 Post on 2021-2-15 10:58:38 | 只看该作者 Only view this author

()

本帖最后由 k64_cc 于 2021-2-15 11:00 编辑

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

SR的弥散函数没那么夸张,运算起来快多了,平均下来耗时比DZVP贵一小小点,但是精度比DZVP要好。对于凝聚相体系来说加大弥散本来就没什么必要,平常用SR就行。

3622

帖子

3

威望

1万

eV
积分
18440

Level 6 (一方通行)

第一原理惨品小作坊

11#
发表于 Post on 2021-2-15 14:50:05 | 只看该作者 Only view this author
本帖最后由 卡开发发 于 2021-2-15 14:51 编辑

哪怕是凝聚相或者表面,BSSE的问题还是需要小心的,可以参看(DOI: 10.1021/acs.jpcc.9b03554),需要注意的是有些结果看上去比较好看可能得益于误差抵消。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。
本周忙

113

帖子

0

威望

1504

eV
积分
1617

Level 5 (御坂)

12#
发表于 Post on 2021-7-4 18:40:27 | 只看该作者 Only view this author
由于DZVP-MOLOPT-GTH基组仅涵盖H , C , N , O , F , Si , P , S , Cl元素,而DZVP-MOLOPT-SR-GTH基组覆盖了元素周期表绝大部分元素,
请问计算同一个体系中DZVP-MOLOPT-GTH基组和DZVP-MOLOPT-SR-GTH基组可以混用吗?
两个对易的厄米算子可以有共同本征函数集

5万

帖子

99

威望

5万

eV
积分
112486

管理员

公社社长

13#
发表于 Post on 2021-7-18 02:35:19 | 只看该作者 Only view this 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-SR-GTH
对弥散函数要求比较高的时候才适合特意花更高代价用不带SR的
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

43

帖子

0

威望

147

eV
积分
190

Level 3 能力者

14#
发表于 Post on 2021-7-22 19:53:34 | 只看该作者 Only view this author
sobereva 发表于 2021-7-18 02:35
没特殊情况就都用DZVP-MOLOPT-SR-GTH
对弥散函数要求比较高的时候才适合特意花更高代价用不带SR的

一个体系内DZVP-MOLOPT-SR-GTH和TZVP-MOLOPT-SR-GTH可以混用吗?

432

帖子

11

威望

3428

eV
积分
4080

Level 6 (一方通行)

15#
发表于 Post on 2021-7-22 22:28:39 | 只看该作者 Only view this author
xxq星球 发表于 2021-7-22 19:53
一个体系内DZVP-MOLOPT-SR-GTH和TZVP-MOLOPT-SR-GTH可以混用吗?

可以,这就是混合基组
自由发挥,野蛮生长

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-26 19:43 , Processed in 0.952292 second(s), 22 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list