计算化学公社

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

[综合交流] 请教三聚体相互作用能的计算问题

[复制链接 Copy URL]

627

帖子

11

威望

1万

eV
积分
12459

Level 6 (一方通行)

跳转到指定楼层 Go to specific reply
#
请问大神,用超分子方法计算三聚体的相互作用能,有下面两个方法(都不考虑BSSE的问题),

(1)   E(ABC) - E(A) - E(B) - E(C)
(2)   E(ABC) - E(AB) - E(BC) - E(AC) + E(A) + E(B) + E(C)

由于方法(2)考虑了耦合项,是不是要更合理一些?
谢谢指导!


计算化学与分子模拟

627

帖子

11

威望

1万

eV
积分
12459

Level 6 (一方通行)

11#
 楼主 Author| 发表于 Post on 2020-9-19 12:01:16 | 只看该作者 Only view this author
sobereva 发表于 2020-9-19 11:37
(1)没问题
(2)如果B-C彼此间的极化作用对于A-B、A-C相互作用能的影响可以忽略的话,可以这么算

好的,明白了, 谢谢大神指导!
计算化学与分子模拟

6万

帖子

99

威望

6万

eV
积分
125181

管理员

公社社长

10#
发表于 Post on 2020-9-19 11:37:34 | 只看该作者 Only view this author
captain 发表于 2020-9-16 18:48
好的,谢谢,基本明白了!
另外如果只是简单的定性的研究,(2)式是否也可以用来说明一下问题呢?

(1)没问题
(2)如果B-C彼此间的极化作用对于A-B、A-C相互作用能的影响可以忽略的话,可以这么算
北京科音自然科学研究中心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

627

帖子

11

威望

1万

eV
积分
12459

Level 6 (一方通行)

9#
 楼主 Author| 发表于 Post on 2020-9-16 18:48:55 | 只看该作者 Only view this author
本帖最后由 captain 于 2020-9-16 18:52 编辑
wzkchem5 发表于 2020-9-16 18:15
对,Int(A-BC)应该以A为一个单体、BC复合物为另一个单体进行计算

好的,谢谢,基本明白了!
另外如果只是简单的定性的研究,(2)式是否也可以用来说明一下问题呢?

还想再听一听社长的高见!
计算化学与分子模拟

1万

帖子

0

威望

9892

eV
积分
22140

Level 6 (一方通行)

8#
发表于 Post on 2020-9-16 18:15:31 | 只看该作者 Only view this author
captain 发表于 2020-9-16 17:04
谢谢热心解答,我没有表达清楚,又编辑了一下帖子。
如果deltaE是基于超分子方法计算的,那么(1)式确实 ...

对,Int(A-BC)应该以A为一个单体、BC复合物为另一个单体进行计算
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
本团队长期招收研究生,有意者可私信联系

627

帖子

11

威望

1万

eV
积分
12459

Level 6 (一方通行)

7#
 楼主 Author| 发表于 Post on 2020-9-16 17:04:31 | 只看该作者 Only view this author
本帖最后由 captain 于 2020-9-16 17:18 编辑
wzkchem5 发表于 2020-9-16 14:02
极化和多体效应都考虑了啊,都在DeltaE(A-BC)里。
把所有的DeltaE都展开,可以得到E(ABC)-E(A)-E(BC)-E( ...

谢谢热心解答,我没有表达清楚,又编辑了一下帖子。
如果deltaE是基于超分子方法计算的,那么(1)式确实是严格的。
如果deltaE是基于SAPT方法计算的,那么(2)式就不严格了吧!
计算化学与分子模拟

1万

帖子

0

威望

9892

eV
积分
22140

Level 6 (一方通行)

6#
发表于 Post on 2020-9-16 14:02:29 | 只看该作者 Only view this author
captain 发表于 2020-9-16 09:38
还想请教大神一下,这篇文献里(J. Chem. Theory Comput. 2006, 2, 1444-1452)提到,
The interaction  ...

极化和多体效应都考虑了啊,都在DeltaE(A-BC)里。
把所有的DeltaE都展开,可以得到E(ABC)-E(A)-E(BC)-E(AC)+E(A)+E(C)=E(ABC)-E(BC)-E(AC)+E(C),其中E(ABC)显然包含了所有的多体效应。
如果改成DeltaEB(AC),那就成为了一个鸡生蛋蛋生鸡的问题——你怎么计算DeltaEB(AC)?是不是又要套DeltaEC(AB)的公式,得到DeltaE(A-BC)-DeltaEB(AC)?跟原式一联立,DeltaEB(AC)和DeltaEC(AB)都消掉了,你啥也没得到。
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
本团队长期招收研究生,有意者可私信联系

627

帖子

11

威望

1万

eV
积分
12459

Level 6 (一方通行)

5#
 楼主 Author| 发表于 Post on 2020-9-16 09:38:01 | 只看该作者 Only view this author
本帖最后由 captain 于 2020-9-16 17:00 编辑

还想请教大神一下,这篇文献里(J. Chem. Theory Comput. 2006, 2, 1444-1452)提到,

The interaction energy between A and B in the presence of C, DeltaEC(A-B),
was calculated as DeltaE(A-BC) - DeltaE(AC)    (1)
感觉这种处理方式没有考虑单体间的极化和多体效应,应该是一种不太严格的计算方法吧?
另外DeltaE(AC)是不是改为DeltaEB(AC),即B存在时A-C相互作用能更合理一些?(但是这样也不好进行实际的计算)

如果Int是由SAPT计算的双体相互作用能,那么下面这个式子没有考虑单体间的极化和多体效应,应该是不严格的吧?
Int(A-BC)=Int(A-B)+Int(A-C)      (2)

计算化学与分子模拟

627

帖子

11

威望

1万

eV
积分
12459

Level 6 (一方通行)

4#
 楼主 Author| 发表于 Post on 2020-6-21 08:17:31 | 只看该作者 Only view this author

谢谢大神指导!
计算化学与分子模拟

6万

帖子

99

威望

6万

eV
积分
125181

管理员

公社社长

3#
发表于 Post on 2020-6-21 02:57:33 | 只看该作者 Only view this author
captain 发表于 2020-6-18 20:24
大神,我看的就是那篇文献。
(2)被称为 the three-body energy,看来意思是三体耦合作用能。
那么用 ...

北京科音自然科学研究中心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

627

帖子

11

威望

1万

eV
积分
12459

Level 6 (一方通行)

2#
 楼主 Author| 发表于 Post on 2020-6-18 20:24:08 | 只看该作者 Only view this author
本帖最后由 captain 于 2020-6-18 20:29 编辑
sobereva 发表于 2020-6-18 20:10
第二个式子不是计算三体结合能的,而是计算三体耦合作用的(即三体结合能中无法被双体作用描述的那部分), ...

大神,我看的就是那篇文献。
(2)被称为 the three-body energy,看来意思是三体耦合作用能。
那么用(1)来计算应该是合理的吧!

计算化学与分子模拟

6万

帖子

99

威望

6万

eV
积分
125181

管理员

公社社长

楼主
发表于 Post on 2020-6-18 20:10:51 | 只看该作者 Only view this author
第二个式子不是计算三体结合能的,而是计算三体耦合作用的(即三体结合能中无法被双体作用描述的那部分),明显二者不是一码事。看DOI: 10.1021/acs.jctc.5b00281的2.3节的描述
北京科音自然科学研究中心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

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

GMT+8, 2026-2-26 02:43 , Processed in 0.163168 second(s), 21 queries , Gzip On.

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