MOF团簇对H2的吸附的优化的问题

筱婷397 · 发表于 Post on 2023-10-11 17:29:29

本帖最后由筱婷397 于 2023-11-22 20:13 编辑

各位老师好，我做的是MOF对H2的吸附，阅读了sob老师簇模型的博文，截取了中间含有金属和氧的团簇加一个配体，边缘用H或CH3饱和，固定所有边缘的C原子。随后构造一个配体缺陷暴露活性位点，保留一个OH。

但是计算的时候发现：①用H饱和的结构，SCF不收敛
②用CH3饱和的结构，几何优化会震荡，而且OH的位置会发生很大变化，与我参考的文献差别较大。
请问各位老师，我这样截取结构进行计算是否合理，可以采用哪些方法解决遇到的问题。（关键词为“#p opt(maxstep=5,NoTrust) pbe1pbe/genecp SCF=(vshift=400,Fermi,conver=6,maxcyc=200)”，附上两种情况的输入输出文件）

输入输出文件.rar (5.67 MB, 下载次数 Times of downloads: 6)

wzkchem5 · 发表于 Post on 2023-10-11 17:51:00

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

sobereva · 发表于 Post on 2023-10-11 17:56:21

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

筱婷397 · 发表于 Post on 2023-10-11 18:13:14

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

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

筱婷397 · 发表于 Post on 2023-10-11 18:16:06

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

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

Childhood · 发表于 Post on 2023-10-11 18:45:08

“随后构造一个配体缺陷暴露活性位点，保留一个OH。”

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

筱婷397 · 发表于 Post on 2023-10-16 10:31:35

本帖最后由筱婷397 于 2023-10-16 10:47 编辑

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

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

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

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

[Gaussian/gview] MOF团簇对H2的吸附的优化的问题