计算化学公社

标题: 过渡态优化过程中震荡，只有RMS Force是YES [打印本页]

作者
Author: liyaohx 时间: 2021-10-25 16:08
标题: 过渡态优化过程中震荡，只有RMS Force是YES
体系是一个含有121原子的Fe的配合物，带了一个冠醚，关键词是
#p opt=(ts,modredundant,noeigentest,readfc,maxstep=5,notrust) geom=check guess=check tpssh/def2svp gfinput pop=no em=GD3BJ IOp(3/174=1000000) IOp(3/175=2238200) IOp(3/177=452900) IOp(3/178=4655000)
读取的力场数也是精确Hessian, 查了一下曲线，发现振荡的主要来源是Displacement(如下图)，请问各位老师该怎么办呢？

作者
Author: liyaohx 时间: 2021-10-25 16:22
后面的IOp主要是B3DJ的参数

作者
Author: wzkchem5 时间: 2021-10-25 16:50
可以考虑加大格点。meta-GGA普遍对格点要求比较高

作者
Author: 喵星大佬 时间: 2021-10-25 16:54
优化过渡态的时候recalc吧

作者
Author: sobereva 时间: 2021-10-25 17:51
上传压缩后的输入、输出文件。结构也不给，初始结构合不合理、有没有优化到无意义的结构上都没法判断
说清楚modredundant是干嘛的

作者
Author: liyaohx 时间: 2021-10-25 21:12
本帖最后由 liyaohx 于 2021-10-25 21:13 编辑

modredundant这个关键词是我复制师姐的，她是要在优化过渡态的时候把之前冻住的键放开，我再这里忘记删除了

。输入和输出文件上传了，多谢老师！
这个过渡态是C-H活化的，想要找到C(75)和H(2)成键的过程。

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