计算化学公社

标题: 溶剂化结构的NBO分析以及离子与溶剂分子之间的二阶微扰能的计算问题 [打印本页]

作者
Author: Dreamer1 时间: 2025-5-5 21:50
标题: 溶剂化结构的NBO分析以及离子与溶剂分子之间的二阶微扰能的计算问题
本帖最后由 Dreamer1 于 2025-5-5 21:52 编辑

老师好，我在分子动力学模拟后得到了电解液中最有可能存在的溶剂化结构：1Li（阳离子）—2FSI(阴离子）—2DME(溶剂分子）。现在我想分析这个结构的NBO，参照用NBO计算原子轨道杂化后的能量变化 - 思想家公社的门口：量子化学·分子模拟·二次元的帖子，我在结构优化的输入文件中写入了““pop=nboread”以及“$NBO plot file=solvation $END”字样，得到了fchk以及.37文件。我的下一步就是想分别详细地计算Li与FSI之间的二阶微扰能、Li与DME之间的二阶微扰能等（如下图参考文献所示）。现在的问题是：（1）这个二阶微扰能是如何通过Multiwfn分析得到的？还是说在别的地方有展示这个数据？（2）文献中每个二阶微扰能对应的轨道怎么确定呢？因为我也要做每个能量所对应的图片，我不知道该如何确定对应的轨道是哪个？  下面附有我优化结构的输入文件的开头和结尾的截图以及参考文献的截图。

   (, 下载次数 Times of downloads: 50)

(, 下载次数 Times of downloads: 50)

    (, 下载次数 Times of downloads: 57)

作者
Author: sobereva 时间: 2025-5-5 22:26
E2值是NBO给出的，此处Multiwfn是用来图形化展示轨道重叠

作者
Author: Dreamer1 时间: 2025-5-5 22:31

sobereva 发表于 2025-5-5 22:26
E2值是NBO给出的，此处Multiwfn是用来图形化展示轨道重叠

好的，老师，这点明白了，那么怎么才能得到E2值呢？这个我不太明白

作者
Author: sobereva 时间: 2025-5-5 22:41

Dreamer1 发表于 2025-5-5 22:31
好的，老师，这点明白了，那么怎么才能得到E2值呢？这个我不太明白

参考量子化学波函数分析与Multiwfn程序培训班（http://www.keinsci.com/WFN）专门讲NBO部分的一页幻灯片：

(, 下载次数 Times of downloads: 54)

作者
Author: Dreamer1 时间: 2025-5-5 22:46

sobereva 发表于 2025-5-5 22:41
参考量子化学波函数分析与Multiwfn程序培训班（http://www.keinsci.com/WFN）专门讲NBO部分的一页幻灯 ...

谢谢老师的指点，我明白了

作者
Author: Dreamer1 时间: 2025-5-8 01:14

sobereva 发表于 2025-5-5 22:41
参考量子化学波函数分析与Multiwfn程序培训班（http://www.keinsci.com/WFN）专门讲NBO部分的一页幻灯 ...

老师，我想再多问一句，您给的这张图里面，每一行的序号是对应的轨道的编号吗？比如第一行的3和13是不是分别对应的轨道编号，那么我要画这个二阶微扰能对应的轨道图，是不是就是画3和13号轨道？

作者
Author: sobereva 时间: 2025-5-8 03:03

Dreamer1 发表于 2025-5-8 01:14
老师，我想再多问一句，您给的这张图里面，每一行的序号是对应的轨道的编号吗？比如第一行的3和13是不是 ...

是。注意NBO plot文件里才对应这个序号。fch文件里的序号会被重排，下文说了这点
使用Multiwfn绘制NBO及相关轨道
http://sobereva.com/134

作者
Author: Dreamer1 时间: 2025-5-8 11:28

sobereva 发表于 2025-5-8 03:03
是。注意NBO plot文件里才对应这个序号。fch文件里的序号会被重排，下文说了这点
使用Multiwfn绘制NBO及 ...

好的，老师，我明白了，感谢您的耐心指导

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