计算化学公社

标题: 用ORCA进行Local Energy Decomposition(LED)分析怎么读frozen state冻结项？ [打印本页]

作者
Author: liruixue 时间: 2025-10-16 16:44
标题: 用ORCA进行Local Energy Decomposition(LED)分析怎么读frozen state冻结项？
本帖最后由 liruixue 于 2025-10-16 17:00 编辑

各位老师好，我在用DLPNO-CCSD(T)/def2-TZVP进行LED分析时，可以得出各个能量项，但是进行具体的共价程度讨论时，从下面表中的能量项中貌似分析不了共价程度。在读LED相关文献时看到静电和交换项等还可以接着分解，想着是不是可以从更深一层的能量分解项的轨道弛豫校正项进行解释呢？我在手册里看到需要gbw提供波函数文件，便用高斯优化(M06l/def2-TZVP)出的二体复合物的fch文件转为gbw文件用于ORCA计算，但我在orca计算得出的out文件中不会读frozen state冻结项，求大神帮忙

作者
Author: wzkchem5 时间: 2025-10-16 18:00
frozen electrostatics就是你截图的手册第二段做完后的out文件里的electrostatics能量。正常计算的electrostatics，减去frozen state计算的electrostatics，等于induced electrostatics。frozen electronic preparation, frozen exchange同理。

作者
Author: liruixue 时间: 2025-10-16 19:20

wzkchem5 发表于 2025-10-16 18:00
frozen electrostatics就是你截图的手册第二段做完后的out文件里的electrostatics能量。正常计算的electros ...

老师，比如frozen state计算的electrostatics是-0.362771255051a.u.，正常计算的electrostatics是-0.107945737390a.u.,得到的induced electrostatics是正值0.254825518a.u.，这种情况是正常的吗

作者
Author: wzkchem5 时间: 2025-10-16 20:14

liruixue 发表于 2025-10-16 19:20
老师，比如frozen state计算的electrostatics是-0.362771255051a.u.，正常计算的electrostatics是-0.1079 ...

感觉induced electrostatics可能偏大，但不是完全不可能

作者
Author: liruixue 时间: 2025-10-16 21:18
本帖最后由 liruixue 于 2025-10-16 21:53 编辑

wzkchem5 发表于 2025-10-16 20:14
感觉induced electrostatics可能偏大，但不是完全不可能

老师我还有问题，为什么我算了ClAg(etu)3和I2之间的LED分析得到的I2的Intra REF. energy差值过大呢？LED正常算的I2的Intra REF. energy是-593.257，但frozen state计算的I2的Intra REF. energy是-168.463。frozen electronic preparation项我不太会看，正常的electronic preparation项我是看的手册里的，分别读取LED的out文件两个单体的Intra REF. energy分别减去复合物中单体单点的能量再相加得到的。老师麻烦您再看一下我上面算冻结项的输入文件有没有问题呢

欢迎光临计算化学公社 (http://bbs.keinsci.com/)