计算化学公社

标题: 求助铁和有机酸的IRC跑完之后铁与有机为什么没有分开？ [打印本页]

作者
Author: studyday 时间: 2023-6-5 20:05
标题: 求助铁和有机酸的IRC跑完之后铁与有机为什么没有分开？
各位老师好，采用B3LYP/6-31+g(d,p)对结构进行优化、频率分析以及IRC分析，相关输出log文件已附上。
几何优化参数： # opt=(calcfc,rfo,ts,noeigen) freq rb3lyp/6-31+g(d,p) nosymm geom=connectivity scf=noincfock
IRC参数： # irc=(calcfc,maxpoints=80,lqa) ub3lyp/6-31+g(d,p) geom=connectivity scf=noincfock
结构优化使用的初始结构为MS-Dmol3找到的过渡态，然后用Gaussian09进行OPT+Fr计算，收敛之后，虚频数为1，然后进行IRC，IRC结果如下，能量好像稳定了，但是Fe3+ 和有机酸并没有分开。
想请问各位老师，计算参数设置是否可行？出现这个结果是错的吗？是结构优化的初始结构不对吗？
感谢各位老师。

(, 下载次数 Times of downloads: 11) (, 下载次数 Times of downloads: 12) (, 下载次数 Times of downloads: 11)
(, 下载次数 Times of downloads: 11) (, 下载次数 Times of downloads: 10) (, 下载次数 Times of downloads: 11)

作者
Author: wzkchem5 时间: 2023-6-5 20:11
计算设置严重不合理，你现在的计算结果是没法用的。
（1）这种有机酸和Fe3+配位是需要脱质子的
（2）你的电荷和自旋多重度都设错了，建议把金属有机、配位化学的教材翻出来，重新学一下怎么数电子，怎么判断一个配合物是高自旋还是低自旋
（3）必须加隐式溶剂模型

作者
Author: studyday 时间: 2023-6-5 20:32
本帖最后由 studyday 于 2023-6-5 20:45 编辑

感谢老师的回复。好的老师，我下去再学习学习。再次感谢老师！

作者
Author: xexlalalan 时间: 2023-6-6 03:56
没看你输入文件。单从给的结构截图上来看，你IRC跑重了，前后是往一个方向跑的，都是跑到反应物去了。按楼上说的改好并把结构按正常收敛限收敛之后，加个calcall和tight收敛收敛，再跑下irc试试

作者
Author: studyday 时间: 2023-6-6 20:01

xexlalalan 发表于 2023-6-6 03:56
没看你输入文件。单从给的结构截图上来看，你IRC跑重了，前后是往一个方向跑的，都是跑到反应物去了。按楼 ...

好的老师，感谢！

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