计算化学公社

 找回密码 Forget password
 注册 Register
楼主 Author: sobereva
打印 Print 上一主题 Last thread 下一主题 Next thread

[量化理论] 谈谈量子化学中基组的选择

  [复制链接 Copy URL]

38

帖子

0

威望

265

eV
积分
303

Level 3 能力者

31#
发表于 Post on 2021-3-9 15:30:20 | 只看该作者 Only view this author
请问sob老师,我在做关于C6H6分子双电离后的反应路径,在寻找第一步开环的过渡态是先用便宜的b3lyp/6-31g(d) 成功的找到了过渡态,但是在此基础上换成m062x/def2tzvp后始终无法收敛,请问老师是基组的问题吗,谢谢老师!

487

帖子

1

威望

4023

eV
积分
4530

Level 6 (一方通行)

32#
发表于 Post on 2021-3-9 15:38:47 | 只看该作者 Only view this author
lyt0226 发表于 2021-3-9 15:30
请问sob老师,我在做关于C6H6分子双电离后的反应路径,在寻找第一步开环的过渡态是先用便宜的b3lyp/6-31g(d ...

什么是双电离,是指电离掉两个电子的苯分子?那样的话电离掉的应该是π电子,σ键键级应该不受影响,为什么会开环?
现代化学以狄拉克的一句“一切化学问题业已解决”为嚆矢。滥觞于经验主义传统的期望正失去它们的借鉴意义。但面对看似不可达的通往天堂之阶梯,我想循伍德沃德“最好的模型是你底物的对映异构体”的信仰好过过早地振翮。
我们怀揣热忱的灵魂天然被赋予对第一性的追求,不屑于单一坐标的约束,钟情于势能面彼端的芬芳。但

38

帖子

0

威望

265

eV
积分
303

Level 3 能力者

33#
发表于 Post on 2021-3-9 17:11:37 | 只看该作者 Only view this author
本帖最后由 lyt0226 于 2021-3-9 17:17 编辑
北大-陶豫 发表于 2021-3-9 15:38
什么是双电离,是指电离掉两个电子的苯分子?那样的话电离掉的应该是π电子,σ键键级应该不受影响,为什 ...

我目前所理解的双电离就只是将初始输入文件中苯的电荷态设置为2,不知道是不是理解有误?(因为想找的是二价苯碎裂成两个一价离子的反应路径:类似于C6H6 2+  ——> CH3+ +C5H3+)
我是看目前的文献中关于该路径的研究都是以开环作为第一个过渡态,我想把第一步开环和质子迁移的过渡态都算出来看哪一个势垒更低。

1万

帖子

0

威望

7390

eV
积分
18143

Level 6 (一方通行)

34#
发表于 Post on 2021-3-9 19:55:33 | 只看该作者 Only view this author
lyt0226 发表于 2021-3-9 15:30
请问sob老师,我在做关于C6H6分子双电离后的反应路径,在寻找第一步开环的过渡态是先用便宜的b3lyp/6-31g(d ...

这种体系还是多参考态方法靠谱,建议用CASPT2
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

538

帖子

2

威望

2476

eV
积分
3054

Level 5 (御坂)

35#
发表于 Post on 2021-3-9 21:57:15 | 只看该作者 Only view this author
sobereva 发表于 2020-10-2 12:42
这种体系用PBE0就足够描述得合理(而且在普通泛函里基本没有明显更好的),完全没必要用CAM-B3LYP。用M06 ...

老师,您好,请问是否可以写个关于使用CP2K计算时,方法和基组选择的帖子呢?

5万

帖子

99

威望

5万

eV
积分
112349

管理员

公社社长

36#
 楼主 Author| 发表于 Post on 2021-3-10 04:15:38 | 只看该作者 Only view this author
风飞 发表于 2021-3-9 21:57
老师,您好,请问是否可以写个关于使用CP2K计算时,方法和基组选择的帖子呢?

目前没打算写
北京科音自然科学研究中心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!

5万

帖子

99

威望

5万

eV
积分
112349

管理员

公社社长

37#
 楼主 Author| 发表于 Post on 2021-3-10 04:17:36 | 只看该作者 Only view this author
lyt0226 发表于 2021-3-9 15:30
请问sob老师,我在做关于C6H6分子双电离后的反应路径,在寻找第一步开环的过渡态是先用便宜的b3lyp/6-31g(d ...

跟基组没有必然联系,跟理论方法关系通常更大
先改用M06-2X/6-31G*再试
北京科音自然科学研究中心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!

38

帖子

0

威望

265

eV
积分
303

Level 3 能力者

38#
发表于 Post on 2021-3-10 11:05:32 | 只看该作者 Only view this author
sobereva 发表于 2021-3-10 04:17
跟基组没有必然联系,跟理论方法关系通常更大
先改用M06-2X/6-31G*再试

谢谢老师的回答,

我在以b3lyp/6-31g(d) 成功找到的开环过渡态结构作为初参条件下,用M06-2X/6-31G*以及m062x/def2tvzp等去寻找开环的过渡态,最终产物除了不收敛以外还总是跑向一个质子发生迁移的状态,这种情况是不是可以说明是找不到这个开环过渡态的

11

帖子

0

威望

506

eV
积分
517

Level 4 (黑子)

39#
发表于 Post on 2021-3-10 11:11:42 | 只看该作者 Only view this author
精彩

5万

帖子

99

威望

5万

eV
积分
112349

管理员

公社社长

40#
 楼主 Author| 发表于 Post on 2021-3-11 01:14:34 | 只看该作者 Only view this author
lyt0226 发表于 2021-3-10 11:05
谢谢老师的回答,

我在以b3lyp/6-31g(d) 成功找到的开环过渡态结构作为初参条件下,用M06-2X/6-31G*以 ...

你也可以试试wB97XD。若能算得动双杂化的话以双杂化为准
另外也可以在B3LYP优化的TS结构下用M06-2X做个振动分析看看虚频情况了解M06-2X在此结构附近的势能面特征
也可以考虑做个柔性扫描,这样便于全面了解不同计算级别下势能面的差异,更容易搞清楚是怎么回事,以及避免漏掉可能的TS
北京科音自然科学研究中心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!

38

帖子

0

威望

265

eV
积分
303

Level 3 能力者

41#
发表于 Post on 2021-3-11 08:39:53 | 只看该作者 Only view this author
sobereva 发表于 2021-3-11 01:14
你也可以试试wB97XD。若能算得动双杂化的话以双杂化为准
另外也可以在B3LYP优化的TS结构下用M06-2X做个 ...

谢谢老师的回答!我尝试一下!

149

帖子

0

威望

1873

eV
积分
2022

Level 5 (御坂)

42#
发表于 Post on 2022-1-4 15:43:56 | 只看该作者 Only view this author
凡是文献里和本文内容有冲突的,一律以本文为准。

22

帖子

0

威望

139

eV
积分
161

Level 3 能力者

43#
发表于 Post on 2022-1-6 13:16:44 | 只看该作者 Only view this author
Sob老师,我采用了B3LYP和MP2等好几种方法和def2tzvp基组对CuAr2进行几何结构优化,得到的基态结构具有C2V对称性,但是其他文献中给的是Dinfh对称的,而我计算得到的Dinfh对称性的结构有一个虚频,请问是我的方法和基组选的有问题吗?

1632

帖子

4

威望

4088

eV
积分
5800

Level 6 (一方通行)

喵星人

44#
发表于 Post on 2022-1-6 13:40:22 | 只看该作者 Only view this author
danni511 发表于 2022-1-6 13:16
Sob老师,我采用了B3LYP和MP2等好几种方法和def2tzvp基组对CuAr2进行几何结构优化,得到的基态结构具有C2V ...

自旋设对了么

22

帖子

0

威望

139

eV
积分
161

Level 3 能力者

45#
发表于 Post on 2022-1-6 13:58:24 | 只看该作者 Only view this author

我算的Cu2+Ar2,电荷和自旋多重度设的是2,2,您看对吗,谢谢

本版积分规则 Credits rule

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

GMT+8, 2024-11-23 05:07 , Processed in 0.179721 second(s), 21 queries , Gzip On.

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