计算化学公社

标题: 铜簇800+原子几何优化结果离谱 [打印本页]

作者
Author: 不是兔纸不是猫 时间: 2022-2-15 18:34
标题: 铜簇800+原子几何优化结果离谱
本帖最后由不是兔纸不是猫于 2022-3-24 14:21 编辑

想请教社长和大家，我用CP2K跑几何优化，铜簇800+原子，初始结构是已经通过单晶xrd得到，跑后结果非常离谱，竟然连苯环都跑乱了，不知道是不是我输入文件哪里设置有问题，我将输入文件贴上来，烦请社长及大家赐教，非常感激！
(, 下载次数 Times of downloads: 41) (, 下载次数 Times of downloads: 45) (, 下载次数 Times of downloads: 52)

作者
Author: Eudaimonia 时间: 2022-2-15 19:11
初始结构直接上传 .inp 文件啊，你这个直接看看不出来

Cu边缘做饱和了吗

作者
Author: 不是兔纸不是猫 时间: 2022-2-15 19:13

Eudaimonia 发表于 2022-2-15 19:11
初始结构直接上传 .inp 文件啊，你这个直接看看不出来

Cu边缘做饱和了吗

我只是删除了坐标部分，其他都截图了，弱弱的问句Cu边缘做饱和是什么意思，应该是没有

作者
Author: Eudaimonia 时间: 2022-2-15 19:43
使用量子化学程序基于簇模型计算金属表面吸附问题
用簇模型的弊端在于：
• 需要合理考虑边界效应。边界效应此处指的是使用有限的簇模型作为周期性体系的近似导致的对被研究问题的计算精度的影响。边界效应取决于两点：(1)簇的尺寸 (2)对边界的处理。簇的尺寸越大、让边界原子所处的环境越像体相，则边界效应越小。
簇的大小选取有很大任意性。固体部分截的原子太少的话边界效应太强，结果不准；而截得大的话，则计算耗时太高。除非团簇截得极大，否则都要考虑怎么恰当处理边界。比如有机体系（如石墨烯），边界往往用氢饱和，见比如《18碳环（cyclo[18]carbon）与石墨烯的相互作用：基于簇模型的研究一例》（http://sobereva.com/615）；如果是无机晶体，有用赝氢、capped ECP、用大量背景电荷表现更长程区域原子的静电势等做法；如果是纯金属的话倒是不需要做特别的考虑。

作者
Author: 不是兔纸不是猫 时间: 2022-2-15 20:00

Eudaimonia 发表于 2022-2-15 19:43
使用量子化学程序基于簇模型计算金属表面吸附问题
用簇模型的弊端在于：
• 需要合理考虑边界效应。 ...

感谢回复，但是我的这个感觉好像跟社长文章中描述还不是一类体系，可能我比较菜，理解有问题，我的铜簇是总共800多个原子，只有几十个铜原子在中心部分，周围一圈全是小分子配体，所以您说的这个给铜做饱和的问题处理方式跟社长这篇文章相似么？是类似体系么？

作者
Author: Eudaimonia 时间: 2022-2-15 20:08

不是兔纸不是猫发表于 2022-2-15 20:00
感谢回复，但是我的这个感觉好像跟社长文章中描述还不是一类体系，可能我比较菜，理解有问题，我的铜簇是 ...

你这个看起来比较像 MOF/COF 一类的东西，可能和这种用簇模型模拟表面不太一样；
你可以看看相关帖子有没有类似的讨论；

作者
Author: 不是兔纸不是猫 时间: 2022-2-15 20:19

Eudaimonia 发表于 2022-2-15 20:08
你这个看起来比较像 MOF/COF 一类的东西，可能和这种用簇模型模拟表面不太一样；
你可以看看相关帖子有 ...

好的同样非常感谢

作者
Author: wzkchem5 时间: 2022-2-15 21:50
检查一下是不是初始结构没有去除无序原子、没有加氢之类的。如果这方面没问题，把一些不重要的取代基用氢代替，跑一下试试

作者
Author: sobereva 时间: 2022-2-16 10:10
直接上传个输入文件，你自己省得解释，别人也省得猜

仔细阅读
在网上求助计算化学问题的时候必须把问题描述得详细、具体、准确、清楚
http://sobereva.com/620（http://bbs.keinsci.com/thread-25787-1-1.html）

作者
Author: 不是兔纸不是猫 时间: 2022-2-16 14:57

sobereva 发表于 2022-2-16 10:10
直接上传个输入文件，你自己省得解释，别人也省得猜

仔细阅读

感谢社长批评，已在一楼上传inp文件，烦请社长百忙之中指点一二，非常感谢！

作者
Author: 不是兔纸不是猫 时间: 2022-2-16 14:58

wzkchem5 发表于 2022-2-15 21:50
检查一下是不是初始结构没有去除无序原子、没有加氢之类的。如果这方面没问题，把一些不重要的取代基用氢代 ...

非常感谢回复，已在一楼上传inp文件，如果可能烦请百忙之中指点一二，感激不尽！

作者
Author: wzkchem5 时间: 2022-2-16 16:49

不是兔纸不是猫发表于 2022-2-16 07:58
非常感谢回复，已在一楼上传inp文件，如果可能烦请百忙之中指点一二，感激不尽！

是所有苯环都跑乱了吗？还是只有一部分跑乱了？
如果是一部分，具体是哪些苯环跑乱了？

作者
Author: chrinide 时间: 2022-2-16 16:54
建议先用xtb预优化一下，再拿来DFT搞

作者
Author: 不是兔纸不是猫 时间: 2022-2-16 18:19

wzkchem5 发表于 2022-2-16 16:49
是所有苯环都跑乱了吗？还是只有一部分跑乱了？
如果是一部分，具体是哪些苯环跑乱了？

全乱了，所以我感觉应该就是我哪里设置的不对

作者
Author: 晓来雨过 时间: 2022-2-16 18:40
如果结构是单晶XRD得到的，那是不是可以固定非氢原子再做优化呢？

作者
Author: 不是兔纸不是猫 时间: 2022-2-16 19:45

晓来雨过发表于 2022-2-16 18:40
如果结构是单晶XRD得到的，那是不是可以固定非氢原子再做优化呢？

感谢回复，那也就是说我可以只优化氢，其他都不需要动是么？不需要全局优化是么？

作者
Author: sobereva 时间: 2022-2-17 08:26

不是兔纸不是猫发表于 2022-2-16 19:45
感谢回复，那也就是说我可以只优化氢，其他都不需要动是么？不需要全局优化是么？

如果晶体解析度够高，只优化氢就行了
你当前体系净电荷设为了+2，注意考虑这个的合理性。其它地方没显著问题
说什么“跑乱了”必须把截图贴出来，否则别人没法猜是什么情况

作者
Author: 晓来雨过 时间: 2022-2-17 12:20

不是兔纸不是猫发表于 2022-2-16 19:45
感谢回复，那也就是说我可以只优化氢，其他都不需要动是么？不需要全局优化是么？

我觉得是的，首先要保证你的晶体结构完整度和GooF值符合要求就好。

作者
Author: 不是兔纸不是猫 时间: 2022-2-17 14:26

sobereva 发表于 2022-2-17 08:26
如果晶体解析度够高，只优化氢就行了
你当前体系净电荷设为了+2，注意考虑这个的合理性。其它地方没显著 ...

好的，感谢社长指点

欢迎光临计算化学公社 (http://bbs.keinsci.com/)