计算化学公社

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

[新手求助] 关于分子之间弱相互作用能计算中的疑问

[复制链接 Copy URL]

54

帖子

0

威望

286

eV
积分
340

Level 3 能力者

各位老师,大家好,请教大家一个问题!

我们在计算两个分子之间的弱相互作用能(Binding Energy) 时,通常用的以下的公式:


Ebind= EA–B - (EA + EB) +BSSE

公式中在计算EA之前,我需要重新对A分子进行结构优化吗?还是直接用A分子在优化后的A-B复合物中的结构来计算能量就行呢?


还请各位老师不吝赐教,感谢!

1665

帖子

5

威望

4788

eV
积分
6553

Level 6 (一方通行)

喵星人

2#
发表于 Post on 2021-1-28 20:17:23 | 只看该作者 Only view this author
那要看你算得是哪种相互作用能

54

帖子

0

威望

286

eV
积分
340

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2021-1-28 20:33:18 | 只看该作者 Only view this author
喵星大佬 发表于 2021-1-28 20:17
那要看你算得是哪种相互作用能

老师您好!  我是想比较一个分子的 羧基位点 是与丙酮分子的作用力更强  还是与乙腈分子的作用力更强

6万

帖子

99

威望

6万

eV
积分
125197

管理员

公社社长

4#
发表于 Post on 2021-1-28 20:48:43 | 只看该作者 Only view this author
需要对A和B都优化。
另外,对于你的实际目的,用结合自由能来讨论更严格。
北京科音自然科学研究中心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

54

帖子

0

威望

286

eV
积分
340

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2021-1-28 21:31:19 | 只看该作者 Only view this author
sobereva 发表于 2021-1-28 20:48
需要对A和B都优化。
另外,对于你的实际目的,用结合自由能来讨论更严格。

sobereva老师,非常感谢您的解答!

我还有一个疑问想请教你一下,我注意到您在许多帖子中提到的一篇您的论文 Exploring Nature and Predicting Strength of Hydrogen Bonds: A Correlation Analysis Between Atoms-in-Molecules Descriptors, Binding Energies, and Energy Components of Symmetry-Adapted Perturbation Theory.

论文里您也算了 binding energy (BE),我注意到您在 Computational Details 中描述到The deformation energy, which is related to the monomer geometry distortion due to the formation of complexes, was not taken into account in the BE evaluations, because BE is best to exhibit intrinsic binding strength of the monomers in a complex without this component. Note that all the monomers considered in present HB set are rather rigid; hence the deformation energy is essentially negligible.

我理解的是这篇论文里的情况是不是没有对 A分子和B分子重新进行结构优化呢?


1665

帖子

5

威望

4788

eV
积分
6553

Level 6 (一方通行)

喵星人

6#
发表于 Post on 2021-1-29 02:47:41 | 只看该作者 Only view this author
Jerry_Wang 发表于 2021-1-28 21:31
sobereva老师,非常感谢您的解答!

我还有一个疑问想请教你一下,我注意到您在许多帖子中提到的一篇您 ...

你可以先看下这篇文章了解下不同的结合能是什么含义

Phys. Chem. Chem. Phys. 2017, 19, 5944

6万

帖子

99

威望

6万

eV
积分
125197

管理员

公社社长

7#
发表于 Post on 2021-1-29 19:17:44 | 只看该作者 Only view this author
Jerry_Wang 发表于 2021-1-28 21:31
sobereva老师,非常感谢您的解答!

我还有一个疑问想请教你一下,我注意到您在许多帖子中提到的一篇您 ...

文中这些体系都是刚性很强的小分子,优不优化都一样,形成二聚体造成的变形能可以忽略不计
北京科音自然科学研究中心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

54

帖子

0

威望

286

eV
积分
340

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2021-1-30 15:30:04 | 只看该作者 Only view this author
sobereva 发表于 2021-1-29 19:17
文中这些体系都是刚性很强的小分子,优不优化都一样,形成二聚体造成的变形能可以忽略不计

明白了! 谢谢 sobereva老师!

sobereva老师,我还有一个问题想请教一下您。

我在优化A-B复合物结构时,一开始用 B97-D3/6-31G(d,p) 来优化,即使用了 opt=tight 发现也会存在虚频,而换成 B3LYP-D3/6-31G(d,p) 来优化,则不会有虚频。我是不是可以理解为,泛函的选择也会影响体系虚频的出现呢?有时候如果出现虚频,我们也可以尝试改变泛函呢?

54

帖子

0

威望

286

eV
积分
340

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2021-1-30 15:30:41 | 只看该作者 Only view this author
喵星大佬 发表于 2021-1-29 02:47
你可以先看下这篇文章了解下不同的结合能是什么含义

Phys. Chem. Chem. Phys. 2017, 19, 5944

谢谢老师,我好好研究一下!感谢!

1万

帖子

0

威望

9895

eV
积分
22145

Level 6 (一方通行)

10#
发表于 Post on 2021-1-30 16:04:06 | 只看该作者 Only view this author
Jerry_Wang 发表于 2021-1-30 15:30
明白了! 谢谢 sobereva老师!

sobereva老师,我还有一个问题想请教一下您。

B97可能对格点比较敏感一些,可以加大格点试试
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
本团队长期招收研究生,有意者可私信联系

6万

帖子

99

威望

6万

eV
积分
125197

管理员

公社社长

11#
发表于 Post on 2021-1-31 12:31:53 | 只看该作者 Only view this author
Jerry_Wang 发表于 2021-1-30 15:30
明白了! 谢谢 sobereva老师!

sobereva老师,我还有一个问题想请教一下您。

很多情况都是数值层面的巧合,不要因为个例就认为在方法层面上有什么明显差别
B97-D3算弱相互作用明显不及B3LYP-D3(BJ)
北京科音自然科学研究中心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

本版积分规则 Credits rule

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

GMT+8, 2026-2-26 18:23 , Processed in 0.171838 second(s), 23 queries , Gzip On.

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