计算化学公社

标题: MOF团簇对H2的吸附的优化的问题 [打印本页]

作者
Author: 筱婷397 时间: 2023-10-11 17:29
标题: MOF团簇对H2的吸附的优化的问题
本帖最后由筱婷397 于 2023-11-22 20:13 编辑

各位老师好，我做的是MOF对H2的吸附，阅读了sob老师簇模型的博文，截取了中间含有金属和氧的团簇加一个配体，边缘用H或CH3饱和，固定所有边缘的C原子。随后构造一个配体缺陷暴露活性位点，保留一个OH。
(, 下载次数 Times of downloads: 4)
但是计算的时候发现：①用H饱和的结构，SCF不收敛
②用CH3饱和的结构，几何优化会震荡，而且OH的位置会发生很大变化，与我参考的文献差别较大。
请问各位老师，我这样截取结构进行计算是否合理，可以采用哪些方法解决遇到的问题。（关键词为“#p opt(maxstep=5,NoTrust) pbe1pbe/genecp SCF=(vshift=400,Fermi,conver=6,maxcyc=200)”，附上两种情况的输入输出文件）
(, 下载次数 Times of downloads: 6)

作者
Author: wzkchem5 时间: 2023-10-11 17:51
确定氢原子数正确吗？如果氢原子数过多或过少，会导致金属价态不对，进而导致难收敛、且即使收敛了结果也不能用。
另外一个小问题：RCOOH里面的R和H一般呈反式比较稳定，你画成顺式了

作者
Author: sobereva 时间: 2023-10-11 17:56
conver=6太糙
优化目的至少用=7。顶多是初期用=6，末尾还得改成=7重新优化

SCF不收敛该怎么解决怎么解决
解决SCF不收敛问题的方法
http://sobereva.com/61
优化震荡该怎么解决怎么解决
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164

镧系配合物本来就难搞，参考
使用Gaussian做镧系金属配合物的量子化学计算
http://sobereva.com/581（http://bbs.keinsci.com/thread-21274-1-1.html）

考虑反铁磁性状态的可能，设置方式参考
谈谈片段组合波函数与自旋极化单重态
http://sobereva.com/82

作者
Author: 筱婷397 时间: 2023-10-11 18:13

wzkchem5 发表于 2023-10-11 17:51
确定氢原子数正确吗？如果氢原子数过多或过少，会导致金属价态不对，进而导致难收敛、且即使收敛了结果也不 ...

氢原子数是正确的，除了我饱和边缘手动加的H，其他的H都是参考了别人的文献中的结构，确定没有错误。我改成反式再试试，感谢老师的回复

作者
Author: 筱婷397 时间: 2023-10-11 18:16

sobereva 发表于 2023-10-11 17:56
conver=6太糙
优化目的至少用=7。顶多是初期用=6，末尾还得改成=7重新优化

谢谢老师的指导，我再好好学习一下这几篇博文

作者
Author: Childhood 时间: 2023-10-11 18:45
“随后构造一个配体缺陷暴露活性位点，保留一个OH。”

UiO-66的缺陷应该不会出现OMS，应该用水去饱和配位。

作者
Author: 筱婷397 时间: 2023-10-16 10:31
本帖最后由筱婷397 于 2023-10-16 10:47 编辑

Childhood 发表于 2023-10-11 18:45
“随后构造一个配体缺陷暴露活性位点，保留一个OH。”

UiO-66的缺陷应该不会出现OMS，应该用水去饱和配 ...

我是参考了一篇文献的结论，构造的这种结构，文献中说a图是最稳定的缺陷构型。我的想法是实验中温度较高之后水会脱掉，暴露金属位点去裂氢。我也还在学习中，有不对的地方还请多指教。
(, 下载次数 Times of downloads: 2)

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