计算化学公社

 找回密码 Forget password
 注册 Register

基于TrEsp计算的激子耦合能中出现的计算问题以及对说明书中公式的理解问题

查看数: 6208 | 评论数: 16 | 收藏 Add to favorites 1
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2021-7-31 10:31

正文摘要:

本帖最后由 量化新生 于 2021-7-31 15:04 编辑 3.8版本的Multiwfn说明书中4.A.9 Calculate TrEsp 3. Evaluating excitonic coupling energy based on TrEsp 补充了基于TrEsp计算的激子耦合能的计算方法以及公式, ...

回复 Reply

895098540wkw 发表于 Post on 6 day ago
你好,请问你算对了吗?这第二步里是计算每个单体计算S0和S1之间的tresp吗
zhangwei2713 发表于 Post on 2025-6-13 08:58:52
sobereva 发表于 2025-6-13 02:21
诸如这作为第一个片段。断C-C键的地方补氢

谢谢sob老师,我试了一下不补氢,结果得到的coupling是0
sobereva 发表于 Post on 2025-6-13 02:21:11
zhangwei2713 发表于 2025-6-12 10:02
谢谢sob老师的回复。比如向附件图片这种PDI二体,其通过一个C原子连接,将其拆分后每个单体都会缺失两个C ...

诸如这作为第一个片段。断C-C键的地方补氢

zhangwei2713 发表于 Post on 2025-6-12 10:02:28
sobereva 发表于 2025-6-12 05:22
通常单体是指一个分子。你的情况需要给出具体结构图

谢谢sob老师的回复。比如向附件图片这种PDI二体,其通过一个C原子连接,将其拆分后每个单体都会缺失两个C原子。
请问此时计算两个单体时,需要分别补足两个氢原子吗?
这篇文章也是采用TrEsp方法进行的coupling计算,但是没有写细节

Picture1.png (72.11 KB, 下载次数 Times of downloads: 116)

Picture1.png
sobereva 发表于 Post on 2025-6-12 05:22:50
zhangwei2713 发表于 2025-6-11 21:49
还有一个问题想请问一下sob老师。采用TrEsp方法计算同一个分子的两个片段(比如两个单体组成的dimer)之 ...

通常单体是指一个分子。你的情况需要给出具体结构图
zhangwei2713 发表于 Post on 2025-6-11 21:49:24
sobereva 发表于 2025-5-29 07:40
算进去

1-14和15-30

还有一个问题想请问一下sob老师。采用TrEsp方法计算同一个分子的两个片段(比如两个单体组成的dimer)之间的coupling时,根据手册(英文版,962-963页)中的介绍需要分别计算两个单体并得到chg文件。我有一个疑问,此时拆分得到的单体需要补足氢原子吗?因为直接将两个单体相连的共价键打断并生成两个单体的gjf文件,不补氢原子则会使自旋变为2。请问一下补和不补哪个才是对的?
sobereva 发表于 Post on 2025-5-29 07:40:58
zhangwei2713 发表于 2025-5-28 10:38
请问一下sob老师,这种计算两个单体之间的耦合时,如果两个单体上各连着一些取代基,计算时是否需要把取 ...

算进去

1-14和15-30
zhangwei2713 发表于 Post on 2025-5-28 10:38:04
sobereva 发表于 2021-8-1 10:34
分清楚3.8和3.8(dev),前者迄今没发布过


请问一下sob老师,这种计算两个单体之间的耦合时,如果两个单体上各连着一些取代基,计算时是否需要把取代基也算进去?比如单体原子序号分别为1-10,15-24,单体取代基分别是11-14,25-30。
我计算的时候写原子序号只写1-10和15-24还是1-14和15-30?因为我发现这两种写法最后结果不一样,不清楚哪个才是准确的结果。请sob老师指导一下。谢谢
sobereva 发表于 Post on 2024-11-14 08:16:27
chunlinxxx 发表于 2024-11-14 07:50
谢谢老大。那所得的值默认就是真空吗

评分 Rate

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

查看全部评分 View all ratings

chunlinxxx 发表于 Post on 2024-11-14 07:50:59
sobereva 发表于 2024-11-14 06:03
Multiwfn用的原子单位,此时没1/(4πε0) ...

谢谢老大。那所得的值默认就是真空吗
sobereva 发表于 Post on 2024-11-14 06:03:48
chunlinxxx 发表于 2024-11-13 18:09
请问老师们,如果只考虑晶体环境,用默认Multiwfn计算需要考虑特定的真空电容率值吗?我看文献公式前面还乘 ...

Multiwfn用的原子单位,此时没1/(4πε0)
chunlinxxx 发表于 Post on 2024-11-13 18:09:28
本帖最后由 chunlinxxx 于 2024-11-13 18:11 编辑

请问老师们,如果只考虑晶体环境,用默认Multiwfn计算需要考虑特定的真空电容率值吗?我看文献公式前面还乘以了一个1/(4πε0)项
sobereva 发表于 Post on 2021-8-1 10:34:37
分清楚3.8和3.8(dev),前者迄今没发布过


就是你说的含义。Multiwfn目前只能基于TrESP电荷近似计算,没法按照你截图里那个式子严格按照积分来计算。
怎么算在4.A.9第三部分已经明确描述了,先用Multiwfn算相应跃迁对应的TrESP电荷,再通过TrESP电荷让Multiwfn直接根据库仑公式计算


就是优化基态二聚体结构(如果你想计算其它有特殊意义的结构下的激子耦合也可以用其它结构)


明显不能直接用B3LYP优化,这是弱相互作用体系最基本的计算常识,仔细看
谈谈量子化学研究中什么时候用B3LYP泛函优化几何结构是适当的
http://sobereva.com/557http://bbs.keinsci.com/thread-17899-1-1.html

你贴的文件我没法正常载入
mol不可能不能上传,例: a.mol (533 Bytes, 下载次数 Times of downloads: 1)
重新用恰当的级别优化、算激发态,如果还有问题,直接上传结构文件我再看


你说的那篇文章我目前没时间看。你可以看看之前一些Multiwfn用户通过Multiwfn算TrESP研究库仑耦合的文章,试图予以重复:
J. Phys. Chem. Lett. 2020, 11, 7908−7913
Chemical Physics https://doi.org/10.1016/j.chemphys.2019.110439
J. Phys. Chem. B 2020, 124, 10653−10662
Phys. Chem. Chem. Phys. DOI: 10.1039/c9cp05753f
NATURE COMMUNICATIONS, (2019) 10, 2048

量化新生 发表于 Post on 2021-7-31 15:06:30
sobereva 发表于 2021-7-31 10:37
有别人回复之前若需要对帖子进行修改、补充,应直接编辑原帖,不要通过回帖进行补充,这点在置顶的新社员 ...

抱歉抱歉,我以为发完贴就编辑不了了。不好意思。主要还是望解答疑惑,萌新下次注意
sobereva 发表于 Post on 2021-7-31 10:37:10
量化新生 发表于 2021-7-31 10:33
mol文件有一丢丢问题,串行了两个,重新附上
DPA-dimer
  MatStudio         3D                         ...

有别人回复之前若需要对帖子进行修改、补充,应直接编辑原帖,不要通过回帖进行补充,这点在置顶的新社员必读贴里明确说了

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

GMT+8, 2025-8-12 18:17 , Processed in 0.281851 second(s), 32 queries , Gzip On.

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