Dmol3计算轨道能量和布居分析问题

柚子不开花

老师您好，这半年刚学习MS，还不会分析outmol结果文件，看到老师的答疑帖子，故冒昧打扰想请教老师两个MS里Dmol3版块小问题。
第一个是，计算HOMO或LUMO轨道能量时，在outmol文件里结果经常显示有UP和down值，有时又直接给出了结果，如下两种结果。我想知道直接给出的结果就是前线轨道能量吗？如果给出的是up和down值，那又该怎么计算前线轨道能量呢？
Energy of Highest Occupied Orbital up:   -0.155397 Ha     -4.229 eV  down:  -0.162915 Ha     -4.433 eV
Energy of Lowest Unoccupied Orbital up:  -0.134810 Ha     -3.668 eV  down:  -0.195532 Ha     -5.321 eV
Energy of Highest Occupied Molecular Orbital  -0.18196Ha    -4.951eV
Energy of Lowest Unoccupied Molecular Orbital  -0.11032Ha    -3.002eV
第二个问题时，mulliken电荷布局分析在outmol文件中结果如下：计算的体系是钴卟啉，请问H（1）为什么有5个都是什么意思？结果中的10、20、21、32、都代表什么意思？即下面结果每列都有什么含义？是用来分析什么的？
Mulliken Population analysis

Population analysis for representation  1    a   
              spin up    down    charge    spin
H (  1) 10 0   0.420   0.420     0.841   0.000
H (  1) 10 0   0.019   0.019     0.038   0.000
H (  1) 21-1   0.006   0.006     0.012   0.000
H (  1) 21 0   0.004   0.004     0.009   0.000
H (  1) 21 1   0.004   0.004     0.009   0.000
C (  2) 10 0   1.000   1.000     2.000   0.000
C (  2) 20 0   0.630   0.630     1.260   0.000
C (  2) 20 0  -0.007  -0.007    -0.015   0.000
C (  2) 21-1   0.447   0.447     0.894   0.000
C (  2) 21-1   0.002   0.002     0.004   0.000
C (  2) 21 0   0.491   0.491     0.982   0.000
C (  2) 21 0  -0.004  -0.004    -0.007   0.000
C (  2) 21 1   0.469   0.469     0.939   0.000
C (  2) 21 1   0.003   0.003     0.007   0.000
C (  2) 32-2   0.007   0.007     0.014   0.000
C (  2) 32-1   0.005   0.005     0.010   0.000
C (  2) 32 0   0.004   0.004     0.008   0.000
C (  2) 32 1   0.004   0.004     0.007   0.000
C (  2) 32 2   0.007   0.007     0.013   0.000
C (  3) 10 0   1.000   1.000     2.000   0.000
C (  3) 20 0   0.630   0.630     1.260   0.000
C (  3) 20 0  -0.007  -0.007    -0.015   0.000
C (  3) 21-1   0.460   0.459     0.919   0.000
期待老师回帖，万分感谢。

sobereva

帖子主题请体现帖子内容，光写个MS别人不知道问的什么，给你改了

非限制性开壳层计算必然需要单独讨论两套轨道，HOMO、LUMO、gap都得对alpha和beta分开讨论


第二个问题这里有详述
谈谈自旋密度、自旋布居以及在Multiwfn中的绘制和计算
http://sobereva.com/353

柚子不开花

sobereva 发表于 2017-1-6 20:24
帖子主题请体现帖子内容，光写个MS别人不知道问的什么，给你改了

非限制性开壳层计算必然需要单独讨论两 ...

明白了，老师，非常感谢。

卡开发发

10、20、21、32表示一个数值轨道，举其中一例，21这个指标实际上是2 1，表示轨道的主量子数n=2 角量子数l=1，后续还有0 1 -1表示p型轨道的伸展方向。

在你的计算中显然使用了极化函数（DNP），如氢原子1s轨道有两个指标一模一样，碳原子的2s和2p也分别有两个指标一模一样的，这是因为采用了双zeta，此时价层电子的轨道采用两个数值轨道描述一个原子轨道。之后对于碳原子多出了3d型的函数，而氢原子多出了2p型的函数，这两个都是角动量和主量子数高出价电子本身的极化基。对于DND而言，则少去氢原子上的2p型数值轨道。

柚子不开花

卡开发发 发表于 2017-1-7 03:11
10、20、21、32表示一个数值轨道，举其中一例，21这个指标实际上是2 1，表示轨道的主量子数n=2 角量子数l=1 ...

哦哦，原来是 这样，之前有帖子说10、 20、21那个代表1S、 2S、 2P、还把我给唬住了。非常感谢您的回答。上面老师回答了前线轨道的问题，我还是不知道up和down值是怎么回事？基础薄弱，老师能不能再深入解释一下，谢谢老师。

卡开发发

柚子不开花 发表于 2017-1-7 09:51
哦哦，原来是 这样，之前有帖子说10、 20、21那个代表1S、 2S、 2P、还把我给唬住了。非常感谢您的回答。 ...

sob前辈在2#其实已经给出了很详尽的解释，但估计你没有仔细看。自旋极化（非限制自旋）计算的时候，alpha电子（即spin up）和beta电子（即spin down）无论是轨道还是密度矩阵等都是分开的，布居分析的时候这些也是分开的，对应的spin up和spin down就是这个，而charge指的是总原子电荷=spin_up+spin_down，spin则是spin_up-spin_down。

你再把2#的帖子好好看看吧。

柚子不开花

谢谢两个前辈，您这么一解释再结合SOb社长的图就懂了，非常谢谢两个前辈。很好的平台，很好的前辈，谢谢。

CocuLiu

您好，看到您也用dmol3 计算了布局分析，看到sob老师的博文中只提到怎样分析高斯计算下的结果文件。
那么请问您对于dmol3计算下的.outmol文件中的结果数据是用什么方法进行分析布局分析的呢？图形是用什么方法绘制的？