计算化学公社

 找回密码 Forget password
 注册 Register
Views: 7100|回复 Reply: 10
打印 Print 上一主题 Last thread 下一主题 Next thread

[xtb] 如何处理文献给出的构型xtb却找不到过渡态

[复制链接 Copy URL]

31

帖子

0

威望

1052

eV
积分
1083

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 绫波レイ 于 2020-5-11 14:25 编辑

测试文献中构型,用文献提供的b3lyp/6-31+G(d),过渡态金属用lanl08(f)能找到过渡态,但是这个构型放到xtb里却找不到过渡态,类似的情况应该怎么处理呢
是因为xtb默认的方法和基组的原因吗,还是这个体系用xtb不太好使的缘故?

3-ts.gjf

1.19 KB, 下载次数 Times of downloads: 14

输入文件

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

11#
发表于 Post on 2020-5-13 07:00:04 | 只看该作者 Only view this author
绫波レイ 发表于 2020-5-11 16:59
优化构型不是在气相而是#p opt=(GDIIS,calcfc,ts,noeigen,recalcfc=8) b3lyp/genecp nosymm scrf=(smd,so ...

差这么多通常要么过渡态没找对,要么是对机理理解有问题,要么数据弄错了,要么忽略了什么因素
什么可能导致催化看具体反应和环境,诸如怀疑溶剂产生催化效应就找过渡态的时候加入显式溶剂
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

31

帖子

0

威望

1052

eV
积分
1083

Level 4 (黑子)

10#
 楼主 Author| 发表于 Post on 2020-5-11 16:59:40 | 只看该作者 Only view this author
sobereva 发表于 2020-5-11 15:08
作为预优化用,小体系也可以用xtb。当前就甭用xtb就完了

没听说过DCE溶剂模型

优化构型不是在气相而是#p opt=(GDIIS,calcfc,ts,noeigen,recalcfc=8) b3lyp/genecp nosymm scrf=(smd,solvent=1,2-Dichloroethane) em=gd3bj temperature=333.15,水也是在这样下优化的,催化剂用的Cu(OTf)2,脱掉了Cu(OTf)2的反应物经过H转移,最后脱水生成最后产物,就是最后这步计算溶质自由能相差80kcal/mol,在溶质自由能文献中看很多过渡态和产物最多也就相差50kcal,所以纠结是不是这样子投稿会被审稿人直接拒了(说的有点啰嗦,见谅),请问如何考虑催化效应能影响呢,谢谢

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

9#
发表于 Post on 2020-5-11 15:08:41 | 只看该作者 Only view this author
绫波レイ 发表于 2020-5-11 15:04
xtb至少得体系50原子以上才有用武之地吗?使用DFT找到的过渡态优化确实结构逐渐跑远了;

另外想请问DC ...

作为预优化用,小体系也可以用xtb。当前就甭用xtb就完了

没听说过DCE溶剂模型
没更具体信息没法判断。注意恰当考虑催化效应
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

31

帖子

0

威望

1052

eV
积分
1083

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 2020-5-11 15:04:57 | 只看该作者 Only view this author
sobereva 发表于 2020-5-11 14:42
这个体系这么小,本来也没必要用xtb
如果怀疑GFN-xTB描述过渡态的能力,用你之前在DFT下找到的过渡态结构 ...

xtb至少得体系50原子以上才有用武之地吗?使用DFT找到的过渡态优化确实结构逐渐跑远了;

另外想请问DCE溶剂模型下优化结构来计算溶质的自由能,在一个脱水反应中过渡态和产物相差80kcal/mol,这是不是说明猜的这个反应错了,因为能量差太大了?谢谢

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

7#
发表于 Post on 2020-5-11 14:42:17 | 只看该作者 Only view this author
这个体系这么小,本来也没必要用xtb
如果怀疑GFN-xTB描述过渡态的能力,用你之前在DFT下找到的过渡态结构当初始结构再试。如果优化过程中逐渐结构跑走了,说明GFN2-xTB不适合描述当前体系的这种过渡态

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

31

帖子

0

威望

1052

eV
积分
1083

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 2020-5-11 14:26:13 | 只看该作者 Only view this author
sobereva 发表于 2020-5-9 02:52
本来xtb就没有找过渡态的功能,说结果之前先说清楚怎么实现的
最好的方式就是3L里提到的做法
所有细节都 ...

感谢卢老师的回复,输入文件已放在1#描述中,谢谢

31

帖子

0

威望

1052

eV
积分
1083

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2020-5-11 14:22:33 | 只看该作者 Only view this author
sobereva 发表于 2020-5-9 02:52
本来xtb就没有找过渡态的功能,说结果之前先说清楚怎么实现的
最好的方式就是3L里提到的做法
所有细节都 ...

在xtb.sh里添加--gfn 2使用GFN2-xTB理论找过渡态,还是没找到过渡态,难道必须用文献的方法和基组才能找到吗?

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

4#
发表于 Post on 2020-5-9 02:52:14 | 只看该作者 Only view this author
本来xtb就没有找过渡态的功能,说结果之前先说清楚怎么实现的
最好的方式就是3L里提到的做法
所有细节都要交代清楚,xtb里具体用的什么理论方法你也没说,GFN0-xTB和GFN2-xTB精度差得十万八千里。
分清楚理论方法和程序,不要把xtb和GFN-xTB理论划等号。
你问的问题和Gaussian一点关系都没有,为什么要选Gaussian/gview分类?已经给你改了,以后注意。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

1102

帖子

18

威望

6643

eV
积分
8105

Level 6 (一方通行)

計算化学の社畜

3#
发表于 Post on 2020-5-8 15:53:54 | 只看该作者 Only view this author
本帖最后由 冰释之川 于 2020-5-9 08:42 编辑

xtb找过渡态??   是用GSM方法么?  你可以试试与高斯联用:将Gaussian与Grimme的xtb程序联用搜索过渡态、产生IRC、做振动分析
http://sobereva.com/421http://bbs.keinsci.com/thread-10106-1-1.html
Stand on the shoulders of giants

831

帖子

1

威望

7185

eV
积分
8036

Level 6 (一方通行)

2#
发表于 Post on 2020-5-8 13:18:25 | 只看该作者 Only view this author
这很正常。dft能算得动就不需要xtb了。

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

GMT+8, 2024-11-24 07:44 , Processed in 0.229137 second(s), 25 queries , Gzip On.

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