计算化学公社

标题: Fe(H2O)2Cl4-结构优化中出现的化学键断裂 [打印本页]

作者
Author: Qing 时间: 2020-6-8 11:11
标题: Fe(H2O)2Cl4-结构优化中出现的化学键断裂
卢老师您好，我在计算FeCl4-和两个水分子配位生成Fe(H2O)2Cl4-的结构的时候，发现最后输出的out文件放入Gaussview中，出现了铁氯和氧氢化学键断裂的情况，可是正常情况下应该不是这样，是铁和两个水分子形成配位键，，您能看一下是怎么回事吗

作者
Author: Qing 时间: 2020-6-8 11:22
老师，输入文件那里打错了，Cl用的是6-311g+**基组，O和H用的是6-311g**基组，其它地方都是没有错误的，计算结果也是按这个基组计算的

作者
Author: thanhtam 时间: 2020-6-8 11:33
本帖最后由 thanhtam 于 2020-6-8 11:35 编辑

先看：谈谈原子间是否成键的判断问题
http://sobereva.com/414（http://bbs.keinsci.com/thread-9840-1-1.html）

其次，6-311g+**的写法是错的，应该是6-311+g**

作者
Author: Qing 时间: 2020-6-8 11:35
好的，谢谢老师

作者
Author: sobereva 时间: 2020-6-8 12:48
lanl2TZ+根本就不是内置基组，你这输入文件根本不可能跑
基组名也写错了。O还定义了两遍，最后还少四个*
你还没给Fe定义赝势

你的输入文件简直就是反面教材
把此文仔细看了
详解Gaussian中混合基组、自定义基组和赝势基组的输入
http://sobereva.com/60

先把基组和赝势修正了再说计算结果。
键的判断问题前面的人提到的我的博文里已经写得很明白了，不要在意gview怎么显示的。

作者
Author: Qing 时间: 2020-6-30 19:28
老师，我在四氢化萘溶剂中对三氯化铁进行几何优化时出现了两个较大的虚频（-103和-89），采用的初始构型为气相中的优化构型，我看了您的帖子后，分别在关键词中添加了opt=tight与opt=calcall的关键词，发现仍然不能够消除虚频。之后，我查看了两个虚频的振动方向，按照振动的方向重新调整了初始构型，但发现此时几何优化报错，出现振荡，我分别使用了采用gdiis优化方法，使用笛卡尔坐标系与内坐标系，使用IOp(1/8=5)调整步长上限，采用PBE0理论方法的方法，结果仍然不能够解决问题，所以想要请教sob老师，有没有什么建议
我采用了b3lypd3理论方法，铁采用lanl2tz基组，氯采用6-311g**基组

作者
Author: sobereva 时间: 2020-7-7 05:17

Qing 发表于 2020-6-30 19:28
老师，我在四氢化萘溶剂中对三氯化铁进行几何优化时出现了两个较大的虚频（-103和-89），采用的初始构型为 ...

没具体输入文件没法判断
如果你用的SMD，改成默认的IEFPCM

欢迎光临计算化学公社 (http://bbs.keinsci.com/)