计算化学公社

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

[新手求助] 计算第四周期过渡金属团簇,用哪些基组比较合适?

[复制链接 Copy URL]

46

帖子

0

威望

326

eV
积分
372

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
计算条件一般看了社长的博文,带赝势基组可能比较好,但是对赝势基组不是很了解

想在耗时和精度直接取个恰当的平衡,哪些基组值得考虑

46

帖子

0

威望

326

eV
积分
372

Level 3 能力者

2#
 楼主 Author| 发表于 Post on 2022-11-25 11:40:16 | 只看该作者 Only view this author
算OPT和能量

977

帖子

0

威望

1929

eV
积分
2906

Level 5 (御坂)

偽語言學家+註冊摸魚師

3#
发表于 Post on 2022-11-25 13:45:23 | 只看该作者 Only view this author

如果你考虑耗时,还是上赝势基组吧。
http://sobereva.com/373
ORCA大法好!算得快,还免费,我爷爷的爷爷都说不赖!
山河自落蒼梧月,風雨猶驅草木兵。
san huɔᵊ ki lɔk tʰɔŋ ŋ ᵑgut,fʊŋ ʝi ʝiu kʰui tʰɔu ᵐbʊk pɪŋ。

46

帖子

0

威望

326

eV
积分
372

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2022-11-25 14:40:00 | 只看该作者 Only view this author
chands 发表于 2022-11-25 13:45
如果你考虑耗时,还是上赝势基组吧。
http://sobereva.com/373

ok 谢谢

5万

帖子

99

威望

5万

eV
积分
112492

管理员

公社社长

5#
发表于 Post on 2022-11-25 20:42:56 | 只看该作者 Only view this author
没细节没法说,看下文Q11了解怎么问
在网上求助计算化学问题的时候必须把问题描述得详细、具体、准确、清楚
http://sobereva.com/620http://bbs.keinsci.com/thread-25787-1-1.html
北京科音自然科学研究中心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!

46

帖子

0

威望

326

eV
积分
372

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2022-11-26 16:45:42 | 只看该作者 Only view this author
sobereva 发表于 2022-11-25 20:42
没细节没法说,看下文Q11了解怎么问
在网上求助计算化学问题的时候必须把问题描述得详细、具体、准确、清 ...

谢谢社长提醒
想要计算130个左右的Zn原子团簇和小分子聚合物之间的吸附能,打算用ORCA程序计算,最多能调用的物理核是20个。希望在算的动的情况下精度尽可能的高。看过您一些博文,对于第四周期有赝势的基组有很多可选的,但之前基本都是全电子基组算,没使用过赝势基组。

312

帖子

2

威望

2306

eV
积分
2658

Level 5 (御坂)

7#
发表于 Post on 2022-11-26 18:33:03 | 只看该作者 Only view this author
cxwudi 发表于 2022-11-26 16:45
谢谢社长提醒
想要计算130个左右的Zn原子团簇和小分子聚合物之间的吸附能,打算用ORCA程序计算,最多能 ...

内存也说一说把

46

帖子

0

威望

326

eV
积分
372

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2022-11-26 19:13:19 | 只看该作者 Only view this author

内存可以调用120G

312

帖子

2

威望

2306

eV
积分
2658

Level 5 (御坂)

9#
发表于 Post on 2022-11-26 21:17:42 | 只看该作者 Only view this author
cxwudi 发表于 2022-11-26 19:13
内存可以调用120G

个人愚见def2-QZVPP算单点能,算的动

46

帖子

0

威望

326

eV
积分
372

Level 3 能力者

10#
 楼主 Author| 发表于 Post on 2022-11-26 21:33:50 | 只看该作者 Only view this author
七尺贱 发表于 2022-11-26 21:17
个人愚见def2-QZVPP算单点能,算的动

我也想着单点的时候用全电子,就是整体算OPT的时候,想着看能不能省一点

5万

帖子

99

威望

5万

eV
积分
112492

管理员

公社社长

11#
发表于 Post on 2022-11-27 01:03:37 | 只看该作者 Only view this author
cxwudi 发表于 2022-11-26 16:45
谢谢社长提醒
想要计算130个左右的Zn原子团簇和小分子聚合物之间的吸附能,打算用ORCA程序计算,最多能 ...

说具体化学组成
Zn的数目很多和很少的情况答案截然不同
北京科音自然科学研究中心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!

46

帖子

0

威望

326

eV
积分
372

Level 3 能力者

12#
 楼主 Author| 发表于 Post on 2022-11-27 20:48:58 | 只看该作者 Only view this author
sobereva 发表于 2022-11-27 01:03
说具体化学组成
Zn的数目很多和很少的情况答案截然不同

264个Zn原子。小分子聚合物12个C,1个N,1个S,1个Na,1个O,16个H。高斯VIEW作图模型在附件里。

Zn_complex.rar

116.46 KB, 下载次数 Times of downloads: 4

1万

帖子

0

威望

7402

eV
积分
18171

Level 6 (一方通行)

13#
发表于 Post on 2022-11-27 21:25:21 | 只看该作者 Only view this author
cxwudi 发表于 2022-11-27 13:48
264个Zn原子。小分子聚合物12个C,1个N,1个S,1个Na,1个O,16个H。高斯VIEW作图模型在附件里。

首先,你这个Zn团簇是实验上就那么大,还是说实验上是宏观尺寸的Zn,抠出一个团簇来近似整个Zn表面?如果是后者的话,建议用embedded cluster方法,参见https://pubs.acs.org/doi/10.1021/acs.inorgchem.9b00994,这样QM区很可能可以取得比现在这样小一些,精度还更高。此外还可以考虑用ONIOM,高级别用DFT,低级别用GFN2-xTB。
其次,建议反复确认是不是想研究表面上完全没有被氧化的Zn。现实中的Zn表面有(相比原子尺寸而言)很厚的氧化层,所以实验测的Zn表面的吸附很可能相当于是吸附在了ZnO或Zn(OH)2上。只有用plasma、惰性气氛下球磨等手段,得到的才是干净的Zn表面,计算用干净的Zn表面才有意义。
再其次,这个分子会不会和Zn表面发生化学反应,也是值得考虑的。
这几个问题确认没问题以后,才可以考虑基组的问题。可以参照http://sobereva.com/373,结合自己是否能算得动来选择。可以先用LANL2DZ预优化看看,如果已经快到可用计算资源的上限了,那就接着用LANL2DZ优化,但是单点得用至少LANL2TZ(f)、def2-TZVP才放心;如果LANL2DZ结构优化尚有余力,那就尝试LANL2TZ。
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

46

帖子

0

威望

326

eV
积分
372

Level 3 能力者

14#
 楼主 Author| 发表于 Post on 2022-11-27 21:50:59 | 只看该作者 Only view this author
wzkchem5 发表于 2022-11-27 21:25
首先,你这个Zn团簇是实验上就那么大,还是说实验上是宏观尺寸的Zn,抠出一个团簇来近似整个Zn表面?如果 ...

谢谢老师指导。实验时Zn表面是打磨过的,保证尽量无氧化的情况下。模型是参见团簇模型抠的。还有个小问题就是用混合基组,对非Zn原子用全电子基组是不是会更好一点?

1万

帖子

0

威望

7402

eV
积分
18171

Level 6 (一方通行)

15#
发表于 Post on 2022-11-27 22:47:54 | 只看该作者 Only view this author
cxwudi 发表于 2022-11-27 14:50
谢谢老师指导。实验时Zn表面是打磨过的,保证尽量无氧化的情况下。模型是参见团簇模型抠的。还有个小问题 ...

对,但是有一定概率算不动。所以不妨赝势基组和全电子基组都试一下,全电子基组能算得动就用全电子基组,反之用赝势基组
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

本版积分规则 Credits rule

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

GMT+8, 2024-11-27 00:06 , Processed in 2.410579 second(s), 24 queries , Gzip On.

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