计算化学公社

标题: 求助sobeda方法能量分解的计算数值偏高 [打印本页]

作者
Author: 假如. 时间: 2025-1-8 15:34
标题: 求助sobeda方法能量分解的计算数值偏高
使用Sobereva老师的sobeda方法进行结构的分解，安装包里的例子能够顺利算完且结果一致，但是针对我自己的结构分解后数值高达-633.51 kcal/mol，感觉格外的高，并且类似结构稳定在同一水平的高，我自己核对了分片段的顺序，没有发现问题，希望各位能帮忙看一下这么高的数值是否正常，或者知道我的问题所在，谢谢（sobEDA.sge.o16512为结果文件）

作者
Author: sobereva 时间: 2025-1-8 20:37
不要自己在标题里手写【求助】这种标签，http://bbs.keinsci.com/thread-9348-1-1.html里明确说了。这次给你改了，以后注意

作者
Author: sobereva 时间: 2025-1-8 20:39
同时形成那么多配位键，而且还是真空下阴阳离子结合，再加上还同时形成了片段间的三个显著的氢键，总相互作用能不很高才怪
如果死活就是怀疑结果，先自己手动在同级别下计算总相互作用能并确保每个计算的波函数的稳定性，如果也是sobEDA.sh直接给出的那个数，就更体现出sobEDA的结果毫无问题

作者
Author: 假如. 时间: 2025-1-9 09:51
本帖最后由假如. 于 2025-1-15 17:24 编辑

sobereva 发表于 2025-1-8 20:39
同时形成那么多配位键，而且还是真空下阴阳离子结合，再加上还同时形成了片段间的三个显著的氢键，总相互作 ...

好的好的，老师我还有个新的问题，我想对相同配体通过不同配位方式的一系列结构进行片段间的总弱相互作用能（包括氢键以及片段间互斥等方面）的对比，来说明某些配位方式的优越性，并且与ΔG方面的得到的结论来相互印证，所以为了排除作用力较强的配位键，我通过不同的片段分法想通过手动减去配位键的作用能以得到其他弱相互作用能的数值，但是得到的二片段分法总作用能更强。
1.由于是锕系原子，是否是我选择的方法基组造成的数据反常？
2.这结果是合理且值得分析的么？如果是合理的从什么方面进行分析较为合适？
3.有没有其他方式能够在排除铀形成的配位键的前提下进行其他部分作用能的分解分析？（因为其实二片段分法也忽略了其他配体有可能与铀酰氧分子形成的氢键作用）

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