计算化学公社

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

[CP2K] MOF表面吸附Li2S4,优化结构小分子跑散了,该如何处理

[复制链接 Copy URL]

7

帖子

0

威望

79

eV
积分
86

Level 2 能力者

Ni基MOF 表面吸附Li2S4模型,用multiwfn生成的输入文件。文献查得4.5K以下MOF为铁磁性,故设置Ni的MAGNETIZATION 为2.0,截取4X3X0.5的晶胞,体系中Ni的个数为48,所以设置MULTIPLICITY为97. 做结构优化,现在还没结束,但是截取最后一帧来看,Li2S4的结构已经跑散了如图,感觉不合理。网上说这种原子跑散了 一般都是输入文件设置有问题,但不知道问题在哪?还是再等等 它还能跑回去?

另外还有一个问题,做表面切片时候要考虑上表面的断键的合理性,那下表面需要吗

屏幕截图 2023-12-22 094042.jpg (85.12 KB, 下载次数 Times of downloads: 34)

屏幕截图 2023-12-22 094042.jpg

mof-9-1.inp

23.71 KB, 下载次数 Times of downloads: 5

7

帖子

0

威望

79

eV
积分
86

Level 2 能力者

13#
 楼主 Author| 发表于 Post on 2024-7-2 10:38:39 | 只看该作者 Only view this author
哪个参数

1

帖子

0

威望

19

eV
积分
20

Level 1 能力者

12#
发表于 Post on 2024-5-21 17:18:58 | 只看该作者 Only view this author
37010405老师您好,请问您还记得您做这个MOF结构优化的参数吗?本小白最近也在做这个结构的优化,但是切完平面和添加真空层后一直无法收敛。

7

帖子

0

威望

79

eV
积分
86

Level 2 能力者

11#
 楼主 Author| 发表于 Post on 2023-12-28 14:49:21 | 只看该作者 Only view this author
喝杯热水 发表于 2023-12-28 10:49
我最近算过的体系。一个就是完整的,另外一个就是不完整的。基底如果跟S强结合,那不完整很正常。

明白,谢谢

7

帖子

0

威望

79

eV
积分
86

Level 2 能力者

10#
 楼主 Author| 发表于 Post on 2023-12-28 14:47:35 | 只看该作者 Only view this author
sobereva 发表于 2023-12-28 10:14
体系不一样没有可比性
当前基底的原子和S成键作用如果太强,自然就会破坏Li2S4结构。你也可以先重复一下 ...

嗯,我对比一下看看,非常感谢

85

帖子

0

威望

735

eV
积分
820

Level 4 (黑子)

9#
发表于 Post on 2023-12-28 10:49:36 | 只看该作者 Only view this author
37010405 发表于 2023-12-27 10:48
我看文献上算出来的锂硫化物结构都是完整的,有办法固定锂硫化物吗

我最近算过的体系。一个就是完整的,另外一个就是不完整的。基底如果跟S强结合,那不完整很正常。

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

8#
发表于 Post on 2023-12-28 10:14:23 | 只看该作者 Only view this author
37010405 发表于 2023-12-27 10:48
我看文献上算出来的锂硫化物结构都是完整的,有办法固定锂硫化物吗

体系不一样没有可比性
当前基底的原子和S成键作用如果太强,自然就会破坏Li2S4结构。你也可以先重复一下已有文章里的计算,看是当前体系的实际特点还是纯粹计算问题
北京科音自然科学研究中心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!

7

帖子

0

威望

79

eV
积分
86

Level 2 能力者

7#
 楼主 Author| 发表于 Post on 2023-12-27 10:48:50 | 只看该作者 Only view this author
喝杯热水 发表于 2023-12-27 10:29
Li2S4这类结构本身稳定性就不是那么好。。

我看文献上算出来的锂硫化物结构都是完整的,有办法固定锂硫化物吗

85

帖子

0

威望

735

eV
积分
820

Level 4 (黑子)

6#
发表于 Post on 2023-12-27 10:29:45 | 只看该作者 Only view this author
Li2S4这类结构本身稳定性就不是那么好。。

7

帖子

0

威望

79

eV
积分
86

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2023-12-27 09:52:22 | 只看该作者 Only view this author
sobereva 发表于 2023-12-22 22:01
用中性状态计算应当更合理

老师 我用中性表面,把距离拉开到3  跑来还是散的

屏幕截图 2023-12-27 094748.jpg (61.03 KB, 下载次数 Times of downloads: 32)

屏幕截图 2023-12-27 094748.jpg

mof-9-5-1.inp

24.3 KB, 下载次数 Times of downloads: 0

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

4#
发表于 Post on 2023-12-22 22:01:16 | 只看该作者 Only view this author
37010405 发表于 2023-12-22 20:52
感谢卢老师回复。静电荷为-2是因为图中这两个位置,原来应该是-OH,考虑到Li应该跟O作用,为了把氧暴露出 ...

用中性状态计算应当更合理
北京科音自然科学研究中心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!

7

帖子

0

威望

79

eV
积分
86

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2023-12-22 20:52:08 | 只看该作者 Only view this author
sobereva 发表于 2023-12-22 17:11
净电荷为什么是-2?

Li2S4一开始摆远一点,看优化成什么

感谢卢老师回复。静电荷为-2是因为图中这两个位置,原来应该是-OH,考虑到Li应该跟O作用,为了把氧暴露出来,所以建模的时候这两个位置删掉了两个H,因此将静电荷设成-2.
另外,我在建模的时候对于周期性体系切表面后的静电荷设置还有一个疑惑,也想请教您。比如我这个体系中,与体相内部的Ni相比,表面的Ni少一个羧基的配位,那设置体系静电荷时候需要考虑每个表面Ni多出来的一个正电荷吗?

好的老师,我试一下把距离提高到3看看

屏幕截图 2023-12-22 203721.png (134.67 KB, 下载次数 Times of downloads: 33)

屏幕截图 2023-12-22 203721.png

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

2#
发表于 Post on 2023-12-22 17:11:20 | 只看该作者 Only view this author
净电荷为什么是-2?

Li2S4一开始摆远一点,看优化成什么
北京科音自然科学研究中心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 09:04 , Processed in 0.199072 second(s), 31 queries , Gzip On.

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