对mofs结构优化的求助

lchperuk

各位老师同学好，我最近在尝试重复一篇文章中对一个mof结构的计算。（文章名：Cerium Metal−Organic Framework for Photocatalysis）
其中对uio-66-ce的cluster modle进行了结构优化，我依照文章中的说法，先用VASP并依照所给的条件优化了它的周期模型，然后从中截取了我们想要的部分并加氢修饰。但是在优化增加的氢原子时遇到了SCF不收敛的问题，我按照sob老师文章的建议加了一些关键词并尝试从较低的基组开始算起，仍然没有解决问题。我联系了这篇文章的作者，他说并没有用什么特殊的收敛方法。因为从来没有在gaussian上算过这么大的结构，所以我想请问使用gaussian对mof的cluster modle进行优化还有结合能计算是否合理，还是说这只是单纯的计算的出现问题，仍然需要按照社长文章中的方法继续尝试使其收敛。谢谢各位。

sobereva

完全是胡算
Ce都没设赝势
用3-21G是搞笑
完全没基本Gaussian和量子化学计算常识

甭用PBE算这个，改成PBE0

把这些都仔细看了
使用Gaussian做镧系金属配合物的量子化学计算
http://sobereva.com/581（http://bbs.keinsci.com/thread-21274-1-1.html）
谈谈量子化学中基组的选择
http://sobereva.com/336（http://bbs.keinsci.com/thread-3545-1-1.html）
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272（http://bbs.keinsci.com/thread-536-1-1.html）
详解Gaussian中混合基组、自定义基组和赝势基组的输入
http://sobereva.com/60

lchperuk

sobereva 发表于 2022-5-31 06:53
完全是胡算
Ce都没设赝势
用3-21G是搞笑

谢谢老师，确实是第一次接触量子化学计算。用3-21g基组是看了您的文章，在scf不收敛的情况下可以考虑从小基组开始逐步提升，实际上原文中用的是6-311+（d）。

sobereva

lchperuk 发表于 2022-5-31 22:02
谢谢老师，确实是第一次接触量子化学计算。用3-21g基组是看了您的文章，在scf不收敛的情况下可以考虑从小 ...

完全没必要带弥散，不要效仿原文

当前体系没必要从3-21G这么小的基组开始去试。直接6-31G*就完了。
PS：第一次接触量子化学计算的人不建议研究这种难搞的体系。

cccxq

老师，请问您在高斯优化MOF结构，选的是不对称单元，还是以金属原子为中心，生长周围配体

Mortal

你好，请问你用的是Materials studio软件做的吗，这个结构怎么截取的