计算化学公社

标题: Dmol3计算轨道能量和布居分析问题 [打印本页]

作者
Author: 柚子不开花 时间: 2017-1-6 15:12
标题: 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
期待老师回帖，万分感谢。

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

非限制性开壳层计算必然需要单独讨论两套轨道，HOMO、LUMO、gap都得对alpha和beta分开讨论
(, 下载次数 Times of downloads: 122)

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

作者
Author: 柚子不开花 时间: 2017-1-6 21:27

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

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

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

作者
Author: 卡开发发 时间: 2017-1-7 03:11
10、20、21、32表示一个数值轨道，举其中一例，21这个指标实际上是2 1，表示轨道的主量子数n=2 角量子数l=1，后续还有0 1 -1表示p型轨道的伸展方向。

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

作者
Author: 柚子不开花 时间: 2017-1-7 09:51

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

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

作者
Author: 卡开发发 时间: 2017-1-7 10:12

柚子不开花发表于 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#的帖子好好看看吧。

作者
Author: 柚子不开花 时间: 2017-1-7 11:06
谢谢两个前辈，您这么一解释再结合SOb社长的图就懂了，非常谢谢两个前辈。很好的平台，很好的前辈，谢谢。

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

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