计算化学公社

标题: HOMO分子轨道图的组成与pop=orbital关键词得到的输出文件中的不同，为什么？ [打印本页]

作者
Author: chunhuaqiushi 时间: 2015-10-12 19:31
标题: HOMO分子轨道图的组成与pop=orbital关键词得到的输出文件中的不同，为什么？
优化了一个Li（cl）2小分子，无论是利用gaussview还是MUltiwfn，得到的HOMO轨道图都显示是由一个原子提供，但是根据pop=orbital关键词得到的HOMO轨道式两个分子共同组成，这中间有什么问题呢？从对称性及电荷分布看，都应该是两个cl原子共同提供HOMO，但是HOMO图只显示一个呢？有没有人能帮我解释及怎么解决问题呢？是计算出错还是显示有问题？多谢了！
输出文件是这样的：19号轨道式HOMO轨道
Atomic contributions to Alpha molecular orbitals:
Alpha occ 8 OE=-7.780 is Cl3-p=1.00
Alpha occ 9 OE=-7.780 is Cl1-p=1.00
Alpha occ 10 OE=-7.780 is Cl1-p=1.00
Alpha occ 11 OE=-2.249 is Li2-s=1.00
Alpha occ 12 OE=-0.829 is Cl3-s=0.54 Cl1-s=0.51
Alpha occ 13 OE=-0.827 is Cl1-s=0.55 Cl3-s=0.52
Alpha occ 14 OE=-0.261 is Cl3-p=0.47 Cl1-p=0.46
Alpha occ 15 OE=-0.258 is Cl1-p=0.46 Cl3-p=0.45
Alpha occ 16 OE=-0.241 is Cl3-p=0.48 Cl1-p=0.47
Alpha occ 17 OE=-0.241 is Cl3-p=0.48 Cl1-p=0.47
Alpha occ 18 OE=-0.236 is Cl1-p=0.50 Cl3-p=0.50
Alpha occ 19 OE=-0.236 is Cl1-p=0.50 Cl3-p=0.50
Alpha vir 20 OE=0.090 is Li2-p=0.99
Alpha vir 21 OE=0.090 is Li2-p=0.99
Alpha vir 22 OE=0.091 is Li2-s=0.62 Cl3-s=0.23 Cl1-s=0.17
Alpha vir 23 OE=0.091 is Li2-p=0.87 Cl1-s=0.16
Alpha vir 24 OE=0.168 is Li2-s=4.96 Cl3-s=-1.52 Cl1-s=-1.52 Li2-d=-0.53 Cl1-p=-0.19 Cl3-p=-0.19
Alpha vir 25 OE=0.187 is Li2-p=1.38 Cl3-p=-0.21 Cl1-p=-0.21
Alpha vir 26 OE=0.187 is Li2-p=1.38 Cl3-p=-0.21 Cl1-p=-0.21
Alpha vir 27 OE=0.192 is Li2-p=3.57 Cl1-s=-1.60 Cl3-s=-1.60 Cl3-p=0.32 Cl1-p=0.32
Alpha vir 28 OE=0.256 is Cl1-s=1.23 Cl3-s=1.23 Li2-s=-1.12 Li2-d=-0.61 Cl3-p=0.12 Cl1-p=0.12
Alpha vir 29 OE=0.272 is Cl3-p=0.75 Cl1-p=0.75 Li2-d=-0.49
Alpha vir 30 OE=0.272 is Cl3-p=0.75 Cl1-p=0.75 Li2-d=-0.49
Alpha vir 31 OE=0.292 is Li2-s=-4.53 Cl3-s=1.41 Cl1-s=1.41 Cl3-p=0.91 Cl1-p=0.91 Li2-d=0.89
Alpha vir 32 OE=0.292 is Li2-d=1.00
Alpha vir 33 OE=0.292 is Li2-d=1.00
Alpha vir 34 OE=0.294 is Cl1-p=0.83 Cl3-p=0.83 Li2-p=-0.44 Cl3-s=-0.12 Cl1-s=-0.11
Alpha vir 35 OE=0.312 is Cl1-p=0.65 Cl3-p=0.65 Li2-p=-0.31
Alpha vir 36 OE=0.312 is Cl1-p=0.65 Cl3-p=0.65 Li2-p=-0.31
Alpha vir 37 OE=0.332 is Li2-p=0.99
Alpha vir 38 OE=0.332 is Li2-p=0.99

作者
Author: sobereva 时间: 2015-10-12 21:18
你既然都用Multiwfn了，何必还用高斯那个很渣的pop=orbital关键词，直接用Multiwfn做轨道成份分析多好
谈谈轨道成份的计算方法
http://sobereva.com/131

确保轨道编号没搞错。如果用Multiwfn算的轨道成份还是和图对不上，把fch文件发上来。

作者
Author: chunhuaqiushi 时间: 2015-10-13 13:51
本帖最后由 chunhuaqiushi 于 2015-10-13 13:54 编辑

因为想知道HOMO轨道的具体组成，无论从电荷分布还是结构对称性上来看，两个Cl原子都应该是等价的，而且从输出文件的NBO轨道系数上看（直接计算NBO轨道，没用orbital），HOMO轨道中两个Cl原子的p轨道的轨道系数也是一样的，只是符号不同，但HOMO轨道图却只是一个原子，有点对不上。
我用MUtiwfn分析轨道成分，看起来好像只有一个原子提供HOMO轨道。那要是这样那就是轨道没错，只是和电荷及轨道系数不太相符。

作者
Author: sobereva 时间: 2015-10-13 14:13
你分析的到底是NBO还是MO？完全不一样。你如果分析MO就不要牵扯NBO。

作者
Author: chunhuaqiushi 时间: 2015-10-13 17:21
本帖最后由 chunhuaqiushi 于 2015-10-13 17:22 编辑

不好意思，说错了，是 Molecular Orbital Coefficients，这个系数显示HOMO轨道里两个Cl原子的p轨道都是组成部分，

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

作者
Author: sobereva 时间: 2015-10-13 17:29
既然这样，Multiwfn做轨道成分分析肯定两个Cl都有贡献
你得到那样只有一边的结果，要么是轨道编号选错了，要么是你用saveNBO把NBO保存到了fch里，导致分析的不是MO而是NBO。你在主功能0里一看轨道就知道了，Multiwfn分析的结果肯定主功能0看到的轨道分布是一样的。

作者
Author: chunhuaqiushi 时间: 2015-10-13 22:24
多谢了，我是采用的saveNBO，我再试试不用这个关键词的

作者
Author: chunhuaqiushi 时间: 2015-10-13 22:28
非常感谢，的确如此，在不使用saveNBO关键词的情况下，得到的HOMO轨道是两个Cl原子都提供了，非常感谢

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