计算化学公社

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

[新手求助] 体系里有阴离子反应、中性分子反应,弥散是都要加上吗?

[复制链接 Copy URL]

146

帖子

0

威望

435

eV
积分
581

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
各位老师,我想请问一下在一整个反应过程中,有的反应是有阴离子参与,有的反应就都是中性分子反应,那所有反应需要都加上弥散函数吗?是可以选择有阴离子的加,没有阴离子的不加?

谢谢各位老师。

3806

帖子

4

威望

7996

eV
积分
11882

Level 6 (一方通行)

MOKIT开发者

2#
发表于 Post on 2020-3-16 16:48:01 | 只看该作者 Only view this author
同一条反应路径里计算精度必须保持精度一致,比如你给过渡态加了弥散,那么反应物和过渡态也要加弥散。算的动或者懒就全加弥撒,如6-31G(d,p)换成6-31+G(d,p)省事。算不动或者想算得经济一些,只需给涉及阴离子的原子加即可,例如羧酸根离子的两个氧,需使用混合/自定义基组。

当然,一般的有机反应,有时候有阴离子、有时候没有,你不加弥散别人也不会说什么,毕竟只是个结构优化而已。算高精度单点时可以考虑加弥散。
自动做多参考态计算的程序MOKIT

146

帖子

0

威望

435

eV
积分
581

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2020-3-16 17:18:21 | 只看该作者 Only view this author
本帖最后由 wwp 于 2020-3-16 17:53 编辑
zjxitcc 发表于 2020-3-16 16:48
同一条反应路径里计算精度必须保持精度一致,比如你给过渡态加了弥散,那么反应物和过渡态也要加弥散。算的 ...

嗯好的,谢谢老师
那老师,那我可以用b3lyp/6-31+G**优化带阴离子的体系,用b3lyp/6-31G**优化不带阴离子的体系吗?在一整个反应路线中。
谢谢老师了。

3806

帖子

4

威望

7996

eV
积分
11882

Level 6 (一方通行)

MOKIT开发者

4#
发表于 Post on 2020-3-17 00:22:11 | 只看该作者 Only view this author
本帖最后由 zjxitcc 于 2020-3-17 00:24 编辑
wwp 发表于 2020-3-16 17:18
嗯好的,谢谢老师
那老师,那我可以用b3lyp/6-31+G**优化带阴离子的体系,用b3lyp/6-31G**优化不带阴离 ...

不可以,都说了精度要保持一致。

我上面说的“只需给涉及阴离子的原子加即可”不是你这么理解的。举例:羧酸根,假设在反应物里是COO-,在产物里是COOH,那么反应物、过渡态和产物你要么全给氧加弥散,要么全不加。不是指你可以反应物加弥散、产物不加弥散。
自动做多参考态计算的程序MOKIT

146

帖子

0

威望

435

eV
积分
581

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2020-3-17 07:14:39 | 只看该作者 Only view this author
zjxitcc 发表于 2020-3-17 00:22
不可以,都说了精度要保持一致。

我上面说的“只需给涉及阴离子的原子加即可”不是你这么理解的。举例 ...

好的,谢谢老师.

156

帖子

0

威望

554

eV
积分
710

Level 4 (黑子)

6#
发表于 Post on 2021-9-12 20:22:41 | 只看该作者 Only view this author
您好,我想问一下,您的阴离子反应需要在溶剂模型下计算高精度单点能吗?如果需要的话,可以直接使用SMD溶剂模型计算吗?

5万

帖子

99

威望

5万

eV
积分
112384

管理员

公社社长

7#
发表于 Post on 2021-9-29 16:47:10 | 只看该作者 Only view this author
化小白 发表于 2021-9-12 20:22
您好,我想问一下,您的阴离子反应需要在溶剂模型下计算高精度单点能吗?如果需要的话,可以直接使用SMD溶 ...

显然需要
可以用uESE算溶解自由能,加到高精度量化算的气相自由能上。仔细看
比SMD算溶解自由能更好的溶剂模型uESE的使用
http://sobereva.com/593http://bbs.keinsci.com/thread-22879-1-1.html
谈谈隐式溶剂模型下溶解自由能和体系自由能的计算
http://sobereva.com/327http://bbs.keinsci.com/thread-3345-1-1.html
北京科音自然科学研究中心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!

37

帖子

0

威望

316

eV
积分
353

Level 3 能力者

8#
发表于 Post on 2022-6-16 18:43:47 | 只看该作者 Only view this author
zjxitcc 发表于 2020-3-17 00:22
不可以,都说了精度要保持一致。

我上面说的“只需给涉及阴离子的原子加即可”不是你这么理解的。举例 ...

老师,请问一下如果我的反应中基本都是中性物质和阳离子,只有一个阴离子(碳酸根或氢氧根),计算单点能的过程中可以不加弥散函数吗?对计算反应自由能精度的影响可以接受吗?

3806

帖子

4

威望

7996

eV
积分
11882

Level 6 (一方通行)

MOKIT开发者

9#
发表于 Post on 2022-6-16 18:49:31 | 只看该作者 Only view this author
本帖最后由 zjxitcc 于 2024-3-5 14:03 编辑
Qing 发表于 2022-6-16 18:43
老师,请问一下如果我的反应中基本都是中性物质和阳离子,只有一个阴离子(碳酸根或氢氧根),计算单点能 ...

要加(只给阴离子加,比较划算)。如果你已经算完了,又想省时间/偷懒,可以找个决速步,加/不加弥散基函数 都算一下单点,对比一下,若能垒差别很小,说明不需要加,展示此数据给别人/审稿人看(注意,前提是“已经算完了”,如果你还没开始做计算,建议不要偷懒)。

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
一条君 + 2 GJ!

查看全部评分 View all ratings

自动做多参考态计算的程序MOKIT

37

帖子

0

威望

316

eV
积分
353

Level 3 能力者

10#
发表于 Post on 2022-6-16 19:27:29 | 只看该作者 Only view this author
zjxitcc 发表于 2022-6-16 18:49
要加(只给阴离子加,比较划算)。如果你已经算完了,又想省时间/偷懒,可以找个决速步,加/不加弥散基函 ...

好的,谢谢老师!

37

帖子

0

威望

316

eV
积分
353

Level 3 能力者

11#
发表于 Post on 2022-6-21 13:44:11 | 只看该作者 Only view this author
zjxitcc 发表于 2022-6-16 18:49
要加(只给阴离子加,比较划算)。如果你已经算完了,又想省时间/偷懒,可以找个决速步,加/不加弥散基函 ...

老师,还是之前的那个问题,我现在采用的做法是原有的中性分子和阳离子仍然采用原先没有使用弥散函数的单点能数据,对于反应物中的阴离子(氢氧根中的氧原子),以及涉及到该阴离子的产物(水分子中的氧原子),添加弥散函数,请问老师这样的话计算是否合理呢?

1万

帖子

0

威望

7396

eV
积分
18151

Level 6 (一方通行)

12#
发表于 Post on 2022-6-21 16:22:01 | 只看该作者 Only view this author
Qing 发表于 2022-6-21 06:44
老师,还是之前的那个问题,我现在采用的做法是原有的中性分子和阳离子仍然采用原先没有使用弥散函数的单 ...

一般情况下可以。例外情况是如果这样做会导致破坏体系对称性的话,最好在其他原子上也加一些弥散函数来恢复体系对称性。例如反应 HOCH2CH2Cl + OH- -> HOCH2CH2OH + Cl-,如果只给OH-以及对应的产物氧原子加弥散的话,会破坏HOCH2CH2OH的对称性,算后续反应就会发现两个羟基是不等价的。此时就需要给两个氧都加弥散。

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
chands + 2

查看全部评分 View all ratings

BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

37

帖子

0

威望

316

eV
积分
353

Level 3 能力者

13#
发表于 Post on 2022-6-21 16:32:34 | 只看该作者 Only view this author
wzkchem5 发表于 2022-6-21 16:22
一般情况下可以。例外情况是如果这样做会导致破坏体系对称性的话,最好在其他原子上也加一些弥散函数来恢 ...

好的,谢谢老师!

137

帖子

0

威望

3126

eV
积分
3263

Level 5 (御坂)

14#
发表于 Post on 2022-6-21 17:24:26 | 只看该作者 Only view this author
Qing 发表于 2022-6-21 16:32
好的,谢谢老师!

如果算单点能用的Pople基组,注意5d/6d

173

帖子

0

威望

949

eV
积分
1122

Level 4 (黑子)

15#
发表于 Post on 2022-8-3 23:27:26 | 只看该作者 Only view this author
本帖最后由 一条君 于 2022-8-3 23:33 编辑
wzkchem5 发表于 2022-6-21 16:22
一般情况下可以。例外情况是如果这样做会导致破坏体系对称性的话,最好在其他原子上也加一些弥散函数来恢 ...

楼主情况是能量计算中,中性分子和阳离子不加弥散,阴离子加弥散
请问老师,怎么给别人解释这样看似有其他变量的处理呢,(只说http://sobereva.com/119中,阴离子能量需要加,而中性会导致问题所以不加,感觉不能说服;还是告诉别人同一级别可以这样),谢谢老师
科音学员

本版积分规则 Credits rule

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

GMT+8, 2024-11-25 09:43 , Processed in 0.482898 second(s), 21 queries , Gzip On.

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