求助关于使用TrEsp进行激子耦合计算

yangjunfang · 发表于 Post on 2025-12-1 15:21:16

本帖最后由 yangjunfang 于 2025-12-1 15:22 编辑

老师好，我想请教一下关于舒勇TrEsp进行激子耦合计算，我的计算步骤如下：
先将二聚体的两个单体根据该关键词进行计算，#p td IOp(9/40=4) cam-b3lyp/6-311+g**，得到激发态，计算S0->S1的TDM。
然后将fchk输入multiwfn中，根据下面步骤产生TDM.fch
18->9->1->log文件->2->y->y产生两个单体的TDM.fch

然后将TDM.fch输入multiwfn产生每个单体的chg
7->12->5->3->1这样就产生了TrEsp电荷

然后将两个chg文件合并，继续载入multiwfn，通过步骤7->-2得到激子耦合结果如下图，这个很不合理，请问我是哪里有问题嘛

sobereva · 发表于 Post on 2025-12-2 06:57:13

舒勇语义不明

chg文件里的坐标明显离谱，载入Multiwfn进主功能0一看便知

记得恰当使用nosymm避免调整到标准朝向下
谈谈Gaussian中的对称性与nosymm关键词的使用
http://sobereva.com/297（http://bbs.keinsci.com/thread-1355-1-1.html）

yangjunfang · 发表于 Post on 2025-12-2 10:12:10

sobereva 发表于 2025-12-2 06:57
舒勇语义不明

chg文件里的坐标明显离谱，载入Multiwfn进主功能0一看便知

老师，您好，这个坐标确实和初始的坐标不一致，然后我换了和初始坐标一致的坐标，计算的结果 Electrostatic interaction energy: 13466044.64 kJ/mol 3218464.17 kcal/mol，感觉也是不合理的
因为我这个需要考虑对称性，所以没使用nosymm

wal · 发表于 Post on 2025-12-2 16:22:41

你这两个单体互相连一点偏移都没有，优化出来的无虚频结构真是这样吗？可以上传一下结构优化的log么

yangjunfang · 发表于 Post on 2025-12-3 15:50:12

wal 发表于 2025-12-2 16:22
你这两个单体互相连一点偏移都没有，优化出来的无虚频结构真是这样吗？可以上传一下结构优化的log么

我这个分子没有优化

wal · 发表于 Post on 2025-12-3 16:14:52

本帖最后由 wal 于 2025-12-3 16:37 编辑

yangjunfang 发表于 2025-12-3 15:50
我这个分子没有优化

没优化那你算出什么结果都在意料之内

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Multiwfn使用咨询] 求助关于使用TrEsp进行激子耦合计算

浏览过的版块