jiahui 发表于 2024-12-12 21:03
卢老师，我又重新建模用gaussian 优化了一次，参考了http://sobereva.com/406
Gaussian输出文件如zip所示 ...

这输入文件里包围有机物的水分子看起来是随机朝向的，没形成氢键网络，这样Gaussian优化自然困难。

社长那篇博文有点老了，现在更推荐的做法是用Molclus（http://www.keinsci.com/research/molclus.html）结合genmer生成水分子包围有机物的团簇初始构型，再调用xtb在便宜的计算级别下跑动力学模拟、批量优化并去重，待团簇只剩下少数比较合理的候选结构再上Gaussian做DFT计算。

编辑：上述操作处理的孤立团簇尺寸足够的话相当于有机物处于水溶液环境的情况。至于有机物放在冰面的情况，可以考虑用genice生成冰晶体的表面并建模，然后用CP2K等做周期性体系计算，或者截取足够大的冰晶体基底结合适当边缘冻结来按簇模型计算。

Gaussian比Dmol3干这种事强太多了。Dmol3里你选的TS色散校正早就过时了
为什么我非常不建议购买和使用D摩3量子化学程序
http://sobereva.com/508