计算化学公社

标题: 过渡金属催化过渡态算不出来，求帮助 [打印本页]

作者
Author: typ_gauss 时间: 2019-12-23 16:34
标题: 过渡金属催化过渡态算不出来，求帮助
算一个过渡金属催化的过渡态，153个原子：
# opt=(modredundant,noeigen) b3lyp/genecp nosymm 5d scf=conver=6

B 1 2 F
B 2 75 F

C H O N S P 0
6-31G**
****
Ru 0
lanl2dz
****

# opt=(ts,calcfc,nofreeze,noeigen) b3lyp 5d scf=conver=6 nosymm genecp geom=check guess=read freq=noraman

C H O N S P 0
6-31G**
****
Ru 0
lanl2dz
****

冻结的是反应核心的断键和成键，几何优化收敛了，但是过渡态一直算不出来了，看了下out文件是有虚频的，不知道为就是不能收敛，力一般都能收敛，最大位移和根均方位移总是不能收敛，跑了500多个构象，感觉没太大希望，就把进程杀了。请问这种情况是初始构象不对，还是头文件需要修改，谢谢

作者
Author: sobereva 时间: 2019-12-23 18:26
图也没有坐标也没有输出文件收敛曲线也没有，这问题完全unanswerable

关键词就有一堆毛病
写5d完全多余，genecp的时候默认就是5d
conver=6会降低几何优化精度，没事甭用
nosymm完全多余
noraman完全多余，仔细看
常见的多余的和被滥用的Gaussian关键词
http://sobereva.com/331（http://bbs.keinsci.com/thread-3460-1-1.html）

用了genecp连赝势都没指定，结果能有意义就怪了
Gaussian计算根本就没有什么“头文件”的概念。

作者
Author: typ_gauss 时间: 2019-12-23 20:15

sobereva 发表于 2019-12-23 18:26
图也没有坐标也没有输出文件收敛曲线也没有，这问题完全unanswerable

关键词就有一堆毛病

多谢sob老师的指点，本来是想传out文件的，后来发现太大了，上传不上来，我截几个图放在下面。
5d, 和 conver = 6， noraman确实不知道，谢谢老师的提醒。
nosymm是因为这个过渡态我个人感觉是没有什么对称性的，所以就写了这个，不过看您博文即便这样也不用写，我下次把这个关键词去掉。
赝势我是指定了的，我的gjf没贴完整，抱歉。
我把我的gjf上传一下，out文件太大。

谢谢sob老师了。

作者
Author: sobereva 时间: 2019-12-24 00:06

typ_gauss 发表于 2019-12-23 20:15
多谢sob老师的指点，本来是想传out文件的，后来发现太大了，上传不上来，我截几个图放在下面。
5d, 和 c ...

较大的文本文件在上传前先压缩，这在公社首页公告栏里明确提及了

光从你的图什么也看不出来。看实际结构的话估计连过渡态的苗头都没有，8成是初猜严重不合理。优化过程没有往期望的过渡态走的话就不要继续跑，否则纯粹是浪费时间。像你图中的情况跑50步就足够判断肯定没戏了，毫无收敛趋势。

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