计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 第一性原理 (First-principles) › 计算MOF-H2的吸附能输入文件咨询

计算MOF-H2的吸附能输入文件咨询

查看数: 1351 | 评论数: 9 | 收藏 Add to favorites 1

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

honest

发布时间: 2024-12-11 21:41

正文摘要:

各位老师好，我用CP2K计算MOF对氢气的吸附能。然后我准备先几何优化，优化后用restart文件进行能量计算。对这个过程可能理解不全面，所以想请教一下这样是否合理。几何优化的设置如下，PBE（这里有没有必要使用PBE ...

回复 Reply

sobereva 发表于 Post on 2025-1-1 01:12:33

honest 发表于 2024-12-30 10:38
卢老师，在您的建议下，我切换成了CG，最后用了28核，跑了三天，任务结束了。但是energy看起来很震荡，最 ...

从最后部分的收敛情况来看没问题。如果几何结构看起来也合理就可以跑后面的任务

sobereva 发表于 Post on 2024-12-24 13:31:14

honest 发表于 2024-12-24 09:21
卢老师，按照您的之前的回复。
我又重新摆放了H2的位置，进行优化，试了几次还是没有收敛。在这次优化的 ...

对于CP2K的BFGS，几何优化60多步没达到收敛限很正常，只要看到结构还在明显变化就继续跑
用CG往往收敛得更快

sobereva 发表于 Post on 2024-12-13 01:44:15

honest 发表于 2024-12-12 11:26
卢老师，再次请教一下。在您的指导下。我构建的MOF-H2模型几何优化跑完了，但是结构基本没有什么变化，请 ...

被吸附分子位置没怎么变化通常都是收敛限太松，或者一开始摆得太远导致受力太小、过早地满足了收敛限

CP2K的.wfn根本就不是一般意义的.wfn文件，别光看后缀想当然
CP2K的波函数给Multiwfn进行分析必须用molden文件，下文明确说了
详谈使用CP2K产生给Multiwfn用的molden格式的波函数文件
http://sobereva.com/651（http://bbs.keinsci.com/thread-31827-1-1.html）

sobereva 发表于 Post on 2024-12-11 23:18:45

honest 发表于 2024-12-11 22:12
谢谢卢老师的回复。

卢老师，在算能量的时候使用TZ档次的基组可以吧。

优化时候也用更好

不需要

honest 发表于 Post on 2024-12-11 22:12:37

sobereva 发表于 2024-12-11 21:46
没必要用LBFGS，体系又不算很大

给最关键的H可以用TZ档次的基组。其它没问题

谢谢卢老师的回复。

卢老师，在算能量的时候使用TZ档次的基组可以吧。

再就是我是根据GCMC模拟后获得的吸附密度分布图，根据密度分布初步判断的吸附位点，然后我将H2分子放在了对应的位置附近，这个时候需要GFN1-xTB先跑一下吗？然后再进一步进行几何优化？（这个我是参考您沸石吸附甲苯的做法，但是不知道是否适用于我的这个工况）

sobereva 发表于 Post on 2024-12-11 21:46:43

没必要用LBFGS，体系又不算很大

给最关键的H可以用TZ档次的基组。其它没问题

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register