计算化学公社

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

[综合交流] 离子液体体系氢键强度计算问题

[复制链接 Copy URL]

26

帖子

0

威望

657

eV
积分
683

Level 4 (黑子)

(1) 计算几种阳离子A+,B+,C+与H2O之间的结合能(根据结合能定义),采用的基组是B3LYP/6-311+G(d,p), em=d3bj
得到的结合能顺序是:A-H2O<B-H2O<C-H2O
(2)利用sob老师文章J. Comput. Chem., 40, 2868-2881 (2019) DOI: 10.1002/jcc.26068中,在上述得到的fchk文件,通过Multiwfn软件,采用E_HB = -332.34*ρ(BCP)-1.0661计算出来的阳离子-H2O之间的氢键结合能顺序如下:
A-H2O>B-H2O>C-H2O
请问,两个顺序相反,以哪一个为准? 如何判断阳离子与H2O之间相互作用的强弱呢?

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

2#
发表于 Post on 2019-12-5 13:35:19 | 只看该作者 Only view this author
求助帖不要发到“资源分享”版,给你移动了,以后注意

没有具体细节没法判断是否适合用那个预测公式。当前用DFT算的结合能是没有问题的,只要过程正确。
北京科音自然科学研究中心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

26

帖子

0

威望

657

eV
积分
683

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2019-12-5 13:56:43 | 只看该作者 Only view this author
如图为阳离子-H2O优化的几何结构,结构中看到的确有氢键形成,并且做了CVB指数(顺序:A-H2O<B-H2O<C-H2O)

1.jpg (95.55 KB, 下载次数 Times of downloads: 28)

1.jpg

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

4#
发表于 Post on 2019-12-5 16:49:31 | 只看该作者 Only view this author
yuclark 发表于 2019-12-5 13:56
如图为阳离子-H2O优化的几何结构,结构中看到的确有氢键形成,并且做了CVB指数(顺序:A-H2O

这种情况作用能相差肯定非常小,必须用很可靠的计算级别才能充分区分

B3LYP/6-311+G(d,p), em=d3bj用于优化够了,建议拿双杂化结合更大的基组算结合能

至于CVB那些方法,不容易区分出这三个体系作用能的微小差异,甚至可能由于噪音导致顺序是错的
北京科音自然科学研究中心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

95

帖子

0

威望

3396

eV
积分
3491

Level 5 (御坂)

5#
发表于 Post on 2019-12-7 14:08:55 | 只看该作者 Only view this author
sobereva 发表于 2019-12-5 16:49
这种情况作用能相差肯定非常小,必须用很可靠的计算级别才能充分区分

B3LYP/6-311+G(d,p), em=d3bj用 ...

老师,像这种离子液体,只是烷基链的长度不一样,使得相互作用产生微小的差异,应该从哪些方面分析这些差异的原因呢?

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

6#
发表于 Post on 2019-12-7 19:24:03 | 只看该作者 Only view this author
风起时想你 发表于 2019-12-7 14:08
老师,像这种离子液体,只是烷基链的长度不一样,使得相互作用产生微小的差异,应该从哪些方面分析这些差 ...

表面静电势极值点分析是一定要做的,这和作用强度关系最直接,参考
使用Multiwfn+VMD快速地绘制静电势着色的分子范德华表面图和分子间穿透图(含视频演示)
http://sobereva.com/443http://bbs.keinsci.com/thread-11080-1-1.html
使用Multiwfn的定量分子表面分析功能预测反应位点、分析分子间相互作用
http://sobereva.com/159

也可以计算氢的原子电荷予以说明
北京科音自然科学研究中心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

95

帖子

0

威望

3396

eV
积分
3491

Level 5 (御坂)

7#
发表于 Post on 2019-12-8 11:17:27 | 只看该作者 Only view this author
sobereva 发表于 2019-12-7 19:24
表面静电势极值点分析是一定要做的,这和作用强度关系最直接,参考
使用Multiwfn+VMD快速地绘制静电势着 ...

感谢老师的解答。

15

帖子

0

威望

614

eV
积分
629

Level 4 (黑子)

8#
发表于 Post on 2021-4-15 22:28:37 | 只看该作者 Only view this author
本帖最后由 danyj 于 2021-4-15 22:38 编辑
sobereva 发表于 2019-12-7 19:24
表面静电势极值点分析是一定要做的,这和作用强度关系最直接,参考
使用Multiwfn+VMD快速地绘制静电势着 ...

卢老师,那请问一下这种通过计算氢的原子电荷来说明氢键强度的差异选用哪种原子电荷合理呢?看了老师那篇原子电荷计算方法的对比,是否该选用ADCH电荷?在文献中看到有用CHELPG电荷的,这种电荷在老师您的文章中没有提及,又是否可以采用呢?

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

9#
发表于 Post on 2021-4-16 09:12:59 | 只看该作者 Only view this author
danyj 发表于 2021-4-15 22:28
卢老师,那请问一下这种通过计算氢的原子电荷来说明氢键强度的差异选用哪种原子电荷合理呢?看了老师那篇 ...

ADCH、CHELPG都可以
北京科音自然科学研究中心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

15

帖子

0

威望

614

eV
积分
629

Level 4 (黑子)

10#
发表于 Post on 2021-4-16 16:32:54 | 只看该作者 Only view this author
sobereva 发表于 2021-4-16 09:12
ADCH、CHELPG都可以

好的 感谢老师解惑

本版积分规则 Credits rule

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

GMT+8, 2026-2-21 02:42 , Processed in 0.195066 second(s), 23 queries , Gzip On.

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