|
本帖最后由 ganlanchina 于 2024-11-27 08:47 编辑
在优化金属-乙烯配合物过程中,设置了金属与乙烯碳原子初始距离为2.5A左右。进行结构优化。
顺利优化完成,无虚频。输出文件为附件中1-ethene。从结构上看貌似是目标结构。
然而读取log文件,仔细思考,如果确实发生了乙烯与金属的配位,乙烯双键应该被削弱,键长应该变长。
但在输出文件中,”配位乙烯”双键键长(1.35547A)与游离乙烯双键键长(1.35520A)基本无差。很明显,gaussian是没有认为乙烯和金属中心形成配位键来优化的。
于是,我缩短金属与乙烯碳原子距离为2.1A(这个距离文献有报道),再进行优化,仍在计算中。
但我调用还未计算结束的log文件(附件中1-modified),仍然可以发现乙烯正在远离金属中心,目前已经到2.2A了,照这个趋势很显然最后还是会优化成与之前类似的结果。
因为刚刚自学计算,不太理解高斯读取输入文件的逻辑。
我如果希望高斯朝着金属-乙烯配合物方向优化,是应该继续减小初始结构中的金属与乙烯碳原子的距离吗(小于2.1A)?
从更高角度去讨论,高斯究竟是怎么理解输入文件并选择优化的方向的呢?
|
|