计算化学公社

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

[新手求助] 请教ZIF8掺杂Cu2+的结构优化问题

[复制链接 Copy URL]

5

帖子

0

威望

51

eV
积分
56

Level 2 能力者

本帖最后由 jy1165254033 于 2024-9-2 11:31 编辑

各位老师好,我想研究ZIF8(Zn2+)掺杂Cu2+后的电子结构。先得到初始结构4元环ZIF8,成功收敛,后续用Cu2+ 取代Zn2+结构收敛失败
关键词# opt freq b3lyp/genecp em=gd3bj  基组6-31g(d)/LANL2DZ 尝试过添加SCF=vshift=500也不行


另外我想问一下文献中图E的电荷转移图以及电子数量是怎么得到的。
原文链接https://onlinelibrary.wiley.com/doi/10.1002/adfm.202400523

麻烦各位老师帮忙解答,非常感谢!

202409011132346036..png (189.61 KB, 下载次数 Times of downloads: 1)

文献电荷转移图

文献电荷转移图

202409011117221562..png (96.55 KB, 下载次数 Times of downloads: 1)

Cu2+掺杂

Cu2+掺杂

202409011116506532..png (125.26 KB, 下载次数 Times of downloads: 1)

原始ZIF8

原始ZIF8

202409021054257754..png (63.31 KB, 下载次数 Times of downloads: 2)

简化结构

简化结构

Cu-ZIF8-ring4-2.LOG

110.08 KB, 阅读权限: 60, 下载次数 Times of downloads: 2

Cu-ZIF8-ring4-2.gjf

10.83 KB, 阅读权限: 60, 下载次数 Times of downloads: 0

Cu-Zn-opt-1.LOG

64.68 KB, 阅读权限: 60, 下载次数 Times of downloads: 2

Cu-Zn-opt-1.gjf

7.46 KB, 阅读权限: 60, 下载次数 Times of downloads: 0

ZIF8-ring4-1.zip

1.61 MB, 阅读权限: 60, 下载次数 Times of downloads: 0

25

帖子

0

威望

171

eV
积分
196

Level 3 能力者

10#
发表于 Post on 2024-9-5 08:44:59 | 只看该作者 Only view this author
sobereva 发表于 2024-9-4 12:02
没具体输入文件没法判断
都当周期性结构计算了,哪还有悬键的事(除非你就是要计算ZIF-8的表面而非体相 ...

好的,谢谢社长,我先编译popt版试试

5万

帖子

99

威望

5万

eV
积分
112354

管理员

公社社长

9#
发表于 Post on 2024-9-4 12:02:42 | 只看该作者 Only view this author
艾小野 发表于 2024-9-3 20:04
谢谢sob老师,没想到可以得到sob老师的回复。其实我之前也得到了社长老师的电荷选取指点。实操是,Materi ...

没具体输入文件没法判断
都当周期性结构计算了,哪还有悬键的事(除非你就是要计算ZIF-8的表面而非体相)
1000原子纯泛函计算对于CP2K来说太简单了,也耗不了什么内存

照下文编译,用popt版,不会有什么问题
CP2K第一性原理程序在CentOS中的简易安装方法
http://sobereva.com/586http://bbs.keinsci.com/thread-21608-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!

25

帖子

0

威望

171

eV
积分
196

Level 3 能力者

8#
发表于 Post on 2024-9-3 20:04:40 | 只看该作者 Only view this author
sobereva 发表于 2024-9-3 11:33
如果目的是做分子动力学模拟,推荐算REPEAT电荷,看下文的介绍,下文给了产生相应文件的例子
使用Multiw ...

谢谢sob老师,没想到可以得到sob老师的回复。其实我之前也得到了社长老师的电荷选取指点。实操是,Material Studio 构建了3层ZIF-8结构(848个原子,氢原子饱和悬挂键),再用Multiwfn生成cp2k的输入文件来算ZIF周期的REPEAT电荷,但是cp2k运行后报错a non-zero exit code. Per user-direction, the job has been aborted. 看了公社内其他社员帖子说是内存泄漏,需要换另一个版本的MPI。但我还是怀疑应该是我算的原子太多,我看了楼主帖子是只用了四元环来算,所以想请教对于算周期性结构电荷来说,用的结构层数并不会存在很大影响吗?烦请各位前辈不吝赐教。

5

帖子

0

威望

51

eV
积分
56

Level 2 能力者

7#
 楼主 Author| 发表于 Post on 2024-9-3 11:38:41 | 只看该作者 Only view this author
sobereva 发表于 2024-9-3 11:35
就是常规的SCF不收敛,看下文
解决SCF不收敛问题的方法
http://sobereva.com/61

感谢sob老师!蹲到回复了哈哈哈~我再去学习一下

5万

帖子

99

威望

5万

eV
积分
112354

管理员

公社社长

6#
发表于 Post on 2024-9-3 11:35:27 | 只看该作者 Only view this author
jy1165254033 发表于 2024-9-2 11:20
感谢sob老师耐心解答,第一次发帖就有回复有点激动已经上传了输入文件,电荷数是我根据这个帖子中 h ...

就是常规的SCF不收敛,看下文
解决SCF不收敛问题的方法
http://sobereva.com/61

带单电子情况和带的净电荷完全是两码事。严重怀疑你的整体的净电荷设的不对

评分 Rate

参与人数
Participants 1
eV +4 收起 理由
Reason
jy1165254033 + 4 谢谢

查看全部评分 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!

5万

帖子

99

威望

5万

eV
积分
112354

管理员

公社社长

5#
发表于 Post on 2024-9-3 11:33:36 | 只看该作者 Only view this author
艾小野 发表于 2024-9-2 21:50
您好,我最近想用cp2k来算一下ZIF-8的电荷,看了您回sob老师的帖子有点拿不准,我想问一下您是怎么求ZIF的 ...

如果目的是做分子动力学模拟,推荐算REPEAT电荷,看下文的介绍,下文给了产生相应文件的例子
使用Multiwfn非常便利地创建CP2K程序的输入文件
http://sobereva.com/587http://bbs.keinsci.com/thread-21668-1-1.html
如果是其它目的,看下文
使用Multiwfn对周期性体系计算Hirshfeld(-I)、CM5和MBIS原子电荷
http://sobereva.com/712http://bbs.keinsci.com/thread-46269-1-1.html

以上原子电荷在下文里都有介绍
一篇深入浅出、完整全面介绍原子电荷的综述文章已发表!
http://sobereva.com/714http://bbs.keinsci.com/thread-46067-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!

25

帖子

0

威望

171

eV
积分
196

Level 3 能力者

4#
发表于 Post on 2024-9-2 21:50:52 | 只看该作者 Only view this author
您好,我最近想用cp2k来算一下ZIF-8的电荷,看了您回sob老师的帖子有点拿不准,我想问一下您是怎么求ZIF的电荷分布的呢?望您赐教。

5

帖子

0

威望

51

eV
积分
56

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2024-9-2 11:20:13 | 只看该作者 Only view this author
sobereva 发表于 2024-9-2 07:35
每次自觉说报错或者上传输入输出文件
仔细看
在网上求助计算化学问题的时候必须把问题描述得详细、具体、 ...

感谢sob老师耐心解答,第一次发帖就有回复有点激动已经上传了输入文件,电荷数是我根据这个帖子中 http://bbs.keinsci.com/forum.php?mod=viewthread&tid=26131 配体2-甲基咪唑配体带单个电子,12个配体,4个二价金属离子,因此设置-4。原始ZIF8我先用的xtb预优化了一下,然后再用gaussian算的,至于后面的Cu2+掺杂结构,我是直接替换了一个Zn2+,后面又算了一个简化结构同样收敛失败,不知道问题出在哪里,烦请老师指点迷津。

5万

帖子

99

威望

5万

eV
积分
112354

管理员

公社社长

2#
发表于 Post on 2024-9-2 07:35:36 | 只看该作者 Only view this author
每次自觉说报错或者上传输入输出文件
仔细看
在网上求助计算化学问题的时候必须把问题描述得详细、具体、准确、清楚
http://sobereva.com/620http://bbs.keinsci.com/thread-25787-1-1.html
另外,这种体系用CP2K优化好得多得多
并且搞清楚净电荷设置对不对

仔细看下文
谈谈怎么考察、计算、分析化学体系的电子密度
http://sobereva.com/715http://bbs.keinsci.com/thread-46380-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!

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

GMT+8, 2024-11-24 03:03 , Processed in 0.291158 second(s), 26 queries , Gzip On.

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