计算化学公社

 找回密码 Forget password
 注册 Register

用xtb得到了虚频正确过渡态,再用gaussian反而算不出来

查看数: 7893 | 评论数: 19 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2021-11-1 14:42

正文摘要:

首先用的xtb,得到了TS结构,freq和IRC都是合理的 再把TS结构用GAUSSIAN计算,输入文件1是这样的: # opt=(calcfc,ts,noeigen,cartesian) freq b3lyp def2svp em=gd3bj 输出文件1显示结构优化不收敛L9999 于 ...

回复 Reply

sobereva 发表于 Post on 2021-11-4 14:44:43
魔卡少女晗 发表于 2021-11-3 20:35
老师,我用B3LYP/def2svp进行了柔性扫描,扫描范围包括了反应物到生成物的对应键长,这样的结果是不是说 ...

结合轨迹看
如果扫描轨迹确实能体现反应过程主要结构变化特征,而且计算级别合理,可以这么说
sobereva 发表于 Post on 2021-11-4 14:44:01
魔卡少女晗 发表于 2021-11-3 20:32
老师,键长从1.44-2.44,已经超过了反应产物的对应键长了,这样是不是可说明这个反应没有过渡态呢

贴图方式不对,其他人看不到。重新编辑帖子并仔细看置顶的新社员必读贴第4节了解怎么正确贴图,此问题在这里还特意强调了:http://bbs.keinsci.com/thread-18961-1-1.html
魔卡少女晗 发表于 Post on 2021-11-3 20:32:52
本帖最后由 魔卡少女晗 于 2021-11-4 15:13 编辑
wzkchem5 发表于 2021-11-1 16:19
做一下粗略看看倒是可以,但是要想发表的话,柔性扫描证明没有过渡态所用的理论级别必须和同一篇文章的结 ...


sobereva 发表于 Post on 2021-11-2 12:13:26
魔卡少女晗 发表于 2021-11-1 16:07
用BP86/def2-svp这种级别可以吗老师?体系里面一共36个原子,CONH这种的,最重的只有一个I原子。

在def2-SVP这么小基组下,BP86在Gaussian里开了RI也比B3LYP快不了多少。算有机体系用BP86完全莫名其妙,本来这在算有机体系方面就是主流泛函里垫底的。
魔卡少女晗 发表于 Post on 2021-11-2 10:14:11
fatpig 发表于 2021-11-1 17:35
guess做了stable=opt吗?有可能gaussian的guess和xtb不一样。

是说XTB过渡态的波函数稳定性分析吗
RAL 发表于 Post on 2021-11-1 19:00:31
魔卡少女晗 发表于 2021-11-1 18:55
您指的反应中心的原子,是主要发生变化的原子吗

所有发生成、断键的原子
魔卡少女晗 发表于 Post on 2021-11-1 18:55:37
RAL 发表于 2021-11-1 17:29
xtb的过渡态结构和DFT相比通常有系统误差,且方向不固定,直接用DFT优化有一定概率会失败。我一般会先把反 ...

您指的反应中心的原子,是主要发生变化的原子吗
魔卡少女晗 发表于 Post on 2021-11-1 18:54:41
snljty 发表于 2021-11-1 17:02
用Gaussian,BP86应该不会比B3LYP明显快很多,Gaussian的RI做的又不好,纯泛函没那么大优势的。不过什么 ...

同时用的人太多,还得排队
魔卡少女晗 发表于 Post on 2021-11-1 18:54:17
wzkchem5 发表于 2021-11-1 16:19
做一下粗略看看倒是可以,但是要想发表的话,柔性扫描证明没有过渡态所用的理论级别必须和同一篇文章的结 ...

谢谢老师,我回去试一下
fatpig 发表于 Post on 2021-11-1 17:35:32
guess做了stable=opt吗?有可能gaussian的guess和xtb不一样。
RAL 发表于 Post on 2021-11-1 17:29:39
xtb的过渡态结构和DFT相比通常有系统误差,且方向不固定,直接用DFT优化有一定概率会失败。我一般会先把反应中心的原子固定住,用小一点的基组预优化其余的部分,然后放开限制,改回原来的基组优化过渡态
snljty 发表于 Post on 2021-11-1 17:02:41
魔卡少女晗 发表于 2021-11-1 16:07
用BP86/def2-svp这种级别可以吗老师?体系里面一共36个原子,CONH这种的,最重的只有一个I原子。

用Gaussian,BP86应该不会比B3LYP明显快很多,Gaussian的RI做的又不好,纯泛函没那么大优势的。不过什么服务器36个原子都算不动?
wzkchem5 发表于 Post on 2021-11-1 16:19:03
魔卡少女晗 发表于 2021-11-1 09:07
用BP86/def2-svp这种级别可以吗老师?体系里面一共36个原子,CONH这种的,最重的只有一个I原子。

做一下粗略看看倒是可以,但是要想发表的话,柔性扫描证明没有过渡态所用的理论级别必须和同一篇文章的结构优化理论级别严格一致

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-23 09:28 , Processed in 0.220428 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list