计算化学公社

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

[新手求助] 离子液体催化体系的弥散函数该怎么加?

[复制链接 Copy URL]

206

帖子

0

威望

1860

eV
积分
2066

Level 5 (御坂)

跳转到指定楼层 Go to specific reply
楼主
请教各位老师,我计算的是离子液体催化中性有机分子的反应,一共有62个原子,计算用自建的集群(并行8核),用Gaussian M062X/6-31g** em=gd3优化了4天左右才算完(计算速度太慢)。然后我用ORCA B97-3c预优化一下,然后再用Gaussian M062X/6-31g** em=gd3优化,大概1天能优化完。但是之前没有加弥散函数,而加弥散函数会大大增加计算时间。
我需要计算的是反应过程能量变化以及离子液体与有机分子之间的相互作用,结构优化的时候基组改用加弥散的2-zeta 6-31+g*好还是不加弥散的3-zeta 6-311g**?还是说现在用的6-31g**也还能接受? 另外单点能计算,我是打算用M062X/ma-def2-TZVP。

请各位老师赐教,谢谢!


1万

帖子

0

威望

9888

eV
积分
22133

Level 6 (一方通行)

11#
发表于 Post on 2022-8-4 15:12:24 | 只看该作者 Only view this author
一条君 发表于 2022-8-4 04:13
老师,中性/阳离子A不加弥散,阴离子B加弥散,请问组合AB是组合所有原子加弥散(用一个基组),还是用混合 ...

一般来说,最好的是对B里面带显著负电的原子加弥散,次好的是对B里的所有原子加弥散,再次是对A和B的所有原子都加弥散。但三种做法都是对的,至少如果你没有遇到基组线性相关问题的话是对的。
例外情况是假如阳离子上有一个原子在后续反应里会显著带负电,那需要加弥散。比如离子液体的阳离子上有一个羟基,后续反应里会脱质子变成氧负离子,那这个羟基氧需要加弥散。

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
一条君 + 5 谢谢

查看全部评分 View all ratings

Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

191

帖子

0

威望

1098

eV
积分
1289

Level 4 (黑子)

10#
发表于 Post on 2022-8-4 11:13:15 | 只看该作者 Only view this author
本帖最后由 一条君 于 2022-8-4 11:30 编辑
sobereva 发表于 2021-6-20 18:13
你这体系已经算挺小的了,阴离子部分全都加弥散用6-31+G**,阳离子部分6-31G**
opt freq用B3LYP-D3(BJ), ...

老师,中性/阳离子A不加弥散,阴离子B加弥散,请问组合AB是组合所有原子加弥散(用一个基组),还是用混合基组仅对B加弥散呀,还是两种都可以,谢谢老师(之前忘记混合基组了,都是按AB用一个基组加弥散处理的)
科音学员

206

帖子

0

威望

1860

eV
积分
2066

Level 5 (御坂)

9#
 楼主 Author| 发表于 Post on 2021-6-20 18:38:10 | 只看该作者 Only view this author
sobereva 发表于 2021-6-20 18:32
再正常不过
量子化学研究中切换泛函应当注意的问题
http://sobereva.com/415(http://bbs.keinsci.com/ ...

好的,谢谢您!

6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

8#
发表于 Post on 2021-6-20 18:32:16 | 只看该作者 Only view this author
yaoyuan0711 发表于 2021-6-20 18:30
可能是机器比较老,算的很慢。优化用B3LYP的话,单点能用M062X,这样换泛函好吗?

再正常不过
量子化学研究中切换泛函应当注意的问题
http://sobereva.com/415http://bbs.keinsci.com/thread-9878-1-1.html

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
yaoyuan0711 + 5 谢谢

查看全部评分 View all ratings

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

206

帖子

0

威望

1860

eV
积分
2066

Level 5 (御坂)

7#
 楼主 Author| 发表于 Post on 2021-6-20 18:30:20 | 只看该作者 Only view this author
sobereva 发表于 2021-6-20 18:13
你这体系已经算挺小的了,阴离子部分全都加弥散用6-31+G**,阳离子部分6-31G**
opt freq用B3LYP-D3(BJ), ...

可能是机器比较老,算的很慢。优化用B3LYP的话,单点能用M062X,这样换泛函好吗?

6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

6#
发表于 Post on 2021-6-20 18:13:20 | 只看该作者 Only view this author
你这体系已经算挺小的了,阴离子部分全都加弥散用6-31+G**,阳离子部分6-31G**
opt freq用B3LYP-D3(BJ),比M06-2X快不少,而且可以用更便宜的积分格点
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

206

帖子

0

威望

1860

eV
积分
2066

Level 5 (御坂)

5#
 楼主 Author| 发表于 Post on 2021-6-20 18:13:13 | 只看该作者 Only view this author
wzkchem5 发表于 2021-6-20 17:58
N的负电荷可以共振到O上。不过毕竟共振的比较少,所以O上不加弥散可能也是可以的。如果是我的话,保险起 ...

好的,谢谢您了

1万

帖子

0

威望

9888

eV
积分
22133

Level 6 (一方通行)

4#
发表于 Post on 2021-6-20 17:58:17 | 只看该作者 Only view this author
yaoyuan0711 发表于 2021-6-20 10:35
谢谢指导。还有个问题,NTF2-中不应该是N显著带负电荷吗,为什么O也要加上弥散?

N的负电荷可以共振到O上。不过毕竟共振的比较少,所以O上不加弥散可能也是可以的。如果是我的话,保险起见O上我也会加

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
北大-陶豫 + 5 我很赞同

查看全部评分 View all ratings

Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

206

帖子

0

威望

1860

eV
积分
2066

Level 5 (御坂)

3#
 楼主 Author| 发表于 Post on 2021-6-20 17:35:54 | 只看该作者 Only view this author
wzkchem5 发表于 2021-6-20 15:14
只对那些带显著负电荷的原子加弥散就行了,比如如果你的离子液体阴离子是BF4-,对F加弥散;如果阴离子是NTf ...

谢谢指导。还有个问题,NTF2-中不应该是N显著带负电荷吗,为什么O也要加上弥散?

1万

帖子

0

威望

9888

eV
积分
22133

Level 6 (一方通行)

2#
发表于 Post on 2021-6-20 15:14:47 | 只看该作者 Only view this author
只对那些带显著负电荷的原子加弥散就行了,比如如果你的离子液体阴离子是BF4-,对F加弥散;如果阴离子是NTf2-,对N、O加弥散。这样你加弥散的原子只占原子总数的很小一部分,用3-zeta也可以算得动

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
yaoyuan0711 + 5 谢谢

查看全部评分 View all ratings

Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

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

GMT+8, 2026-2-24 15:42 , Processed in 0.162953 second(s), 22 queries , Gzip On.

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