Gaussian程序算MOF体系的几何优化任务中结构是否合理和SCF不收敛问题

筱婷397

看文献中关于MOFs材料的计算大多是从结构中提取金属节点，我想研究配体上各种取代基团对金属节点上吸附的影响，所以需要带着配体计算，也就是切出一个结构带着周围的配体进行几何优化（包含150个原子），现在有以下两个问题向各位老师请教：

1.金属节点带着配体去完全放开优化，配体部分结构扭曲程度较大，所以打算用sob老师讲过的簇模型冻结原子去处理，不知道冻结图中用蓝色圈出的碳原子是否合理？

2.结构中有Ce，用了PBE0结合SDD赝势进行几何优化，为了帮助SCF收敛，带上了以下关键词：SCF=(vshift=350,Fermi,conver=6,maxcyc=200)，也试了sob老师博文中解决SCF不收敛的方法，SCF还是不收敛，请问是否有其他帮助收敛的解决办法？

计算小白，有说错的地方还请各位老师指正。

wzkchem5

1. 冻结碳原子或冻结羧酸氧原子都有一定道理，可能后者稍好些？这个也看一下其他人的意见
2. 这些Ce有几个Ce(III)，几个Ce(IV)？如果有不止一个Ce(III)，且彼此呈反铁磁耦合的话，需要准备片段初猜

sobereva

1 取决于具体是什么取代基，加在哪里。取代基的不同会明显结构的区域不应当被冻结

2 下文没有可补充的地方
解决SCF不收敛问题的方法
http://sobereva.com/61

对含镧系体系更具体的参考
使用Gaussian做镧系金属配合物的量子化学计算
http://sobereva.com/581（http://bbs.keinsci.com/thread-21274-1-1.html）