计算化学公社

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

[Gaussian/gview] 求助周期体系Def2基组的问题

[复制链接 Copy URL]

43

帖子

0

威望

1936

eV
积分
1979

Level 5 (御坂)

跳转到指定楼层 Go to specific reply
楼主
各位同仁,在计算周期性体系的极性时,我用了Def2TZVP基组计算。

但是老板提醒我,“Regardingthe basis, Def2TZVP is a good basis, but usually for PBC calculations I recallpeople used to remove basis functions with small exponents. I don’t know ifthis is done automatically now in Gaussian, but years ago, having smallexponents was a recipe for very slow calculations and/or convergence problems.


我用的g16,请问有没有类似经验的同仁帮忙解释一下?

3814

帖子

4

威望

8005

eV
积分
11899

Level 6 (一方通行)

MOKIT开发者

2#
发表于 Post on 2020-10-11 09:31:03 | 只看该作者 Only view this author
(1)“计算周期性体系的极性”具体是计算什么?偶极矩?偶极矩还得有弥散函数,你这不够。
(2)你用的g16算PBC?究竟是自己截取的有限团簇,还是真的用了周期性(输入文件里有Tv)?
自动做多参考态计算的程序MOKIT

1237

帖子

1

威望

3501

eV
积分
4758

Level 6 (一方通行)

3#
发表于 Post on 2020-10-11 09:34:56 | 只看该作者 Only view this author
宜用pob-TZVP

43

帖子

0

威望

1936

eV
积分
1979

Level 5 (御坂)

4#
 楼主 Author| 发表于 Post on 2020-10-11 09:38:16 | 只看该作者 Only view this author
zjxitcc 发表于 2020-10-11 09:31
(1)“计算周期性体系的极性”具体是计算什么?偶极矩?偶极矩还得有弥散函数,你这不够。
(2)你用的g1 ...

谢谢回复,抱歉没有描述清楚。
(1) 我希望计算的时一个晶体的表面slab模型,计算偶极矩,后面会用def2TZVPD加上弥散函数。
(2) 真周期性,是2维周期即2个TV,计算垂直于TV方向的偶极。

3814

帖子

4

威望

8005

eV
积分
11899

Level 6 (一方通行)

MOKIT开发者

5#
发表于 Post on 2020-10-11 09:43:24 | 只看该作者 Only view this author
本帖最后由 zjxitcc 于 2020-10-11 09:44 编辑
diffcwmol 发表于 2020-10-11 09:38
谢谢回复,抱歉没有描述清楚。
(1) 我希望计算的时一个晶体的表面slab模型,计算偶极矩,后面会用def2 ...

#3楼是正解,pob-TZVP专门给周期性优化过的,应该用这个。sob老师在他博文里解释了《谈谈量子化学中基组的选择》http://sobereva.com/336
真周期性的话,也不要考虑加弥散了,高斯和Crystal都未必能SCF收敛。

自动做多参考态计算的程序MOKIT

43

帖子

0

威望

1936

eV
积分
1979

Level 5 (御坂)

6#
 楼主 Author| 发表于 Post on 2020-10-11 09:47:58 | 只看该作者 Only view this author

非常感谢,去读文献。

43

帖子

0

威望

1936

eV
积分
1979

Level 5 (御坂)

7#
 楼主 Author| 发表于 Post on 2020-10-11 09:48:13 | 只看该作者 Only view this author
zjxitcc 发表于 2020-10-11 09:43
#3楼是正解,pob-TZVP专门给周期性优化过的,应该用这个。sob老师在他博文里解释了《谈谈量子化学中基组 ...

感谢回复。

1237

帖子

1

威望

3501

eV
积分
4758

Level 6 (一方通行)

8#
发表于 Post on 2020-10-11 10:04:09 | 只看该作者 Only view this author
diffcwmol 发表于 2020-10-11 09:47
非常感谢,去读文献。

PBC的话就别加弥散了,到时候又不好收敛

43

帖子

0

威望

1936

eV
积分
1979

Level 5 (御坂)

9#
 楼主 Author| 发表于 Post on 2020-10-11 10:06:30 | 只看该作者 Only view this author
biogon 发表于 2020-10-11 10:04
PBC的话就别加弥散了,到时候又不好收敛

谢谢宝贵的建议,刚仔细看了sob老师的博文,里面写的比较清楚。我再读读文献理解一下。

5万

帖子

99

威望

5万

eV
积分
112551

管理员

公社社长

10#
发表于 Post on 2020-10-11 16:20:17 | 只看该作者 Only view this author
顺带一提,2019年提出的pob-DZVP-rev2和pob-TZVP-rev2在SCF稳定性上比初代pob基组有了很大提升,而且改进了晶格常数计算结果。各种pob基组可以在https://www.chemie.uni-bonn.de/pctc/mulliken-center/software/ssc_basis/basis获得
北京科音自然科学研究中心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

威望

1936

eV
积分
1979

Level 5 (御坂)

11#
 楼主 Author| 发表于 Post on 2020-10-12 05:22:36 | 只看该作者 Only view this author
sobereva 发表于 2020-10-11 16:20
顺带一提,2019年提出的pob-DZVP-rev2和pob-TZVP-rev2在SCF稳定性上比初代pob基组有了很大提升,而且改进了 ...

谢谢老师提醒

本版积分规则 Credits rule

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

GMT+8, 2024-11-27 21:14 , Processed in 0.171413 second(s), 21 queries , Gzip On.

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