muitiwfn使用和NBO分析

风玲儿

我计算 了从基态到激发态的Electron density difference  ，得到图一的结果，查文献有得到图二的结果，感觉自己的别扭，应如何修改成图二的视图呢？

另外，图二中的定量的确定了电子转移数，说是从NBO中得到的NPA，应如何得到的呢？







        图一                                                                      图二

冰释之川

通过计算基态和垂直激发态下的片段NPA电荷，最后求差得到转移量。
高斯里单点加上 pop=NPA关键词可以获得NPA电荷。
获取激发态NPA电荷 td=(root=1) density pop=NPA。（这里假设获取第一激发态的NPA电荷）
最后的输出文件里你会看到类似如下信息：

Summary of Natural Population Analysis:                  
                                                         
                                       Natural Population
                      Natural  -----------------------------------------------
    Atom  No    Charge         Core      Valence    Rydberg      Total
     -----------------------------------------------------------------------
      C    1   -0.55206      1.99932     4.54336    0.00939     6.55206
      C    2   -0.36114      1.99925     4.34905    0.01284     6.36114
      C    3   -0.38205      1.99921     4.36850    0.01433     6.38205
      C    4   -0.15959      1.99918     4.14114    0.01926     6.15959
      N    5   -0.49285      1.99939     5.47688    0.01658     7.49285
      C    6   -0.24634      1.99905     4.22825    0.01904     6.24634
      C    7   -0.00107      1.99868     3.98868    0.01370     6.00107
      C    8   -0.22222      1.99845     4.21861    0.00516     6.22222
      C    9   -0.24754      1.99906     4.22951    0.01898     6.24754
      C   10   -0.01928      1.99870     4.00694    0.01363     6.01928
      C   11   -0.20169      1.99845     4.19777    0.00547     6.20169
      H   12    0.19303      0.00000     0.80572    0.00125     0.80697
      H   13    0.18685      0.00000     0.81140    0.00175     0.81315
      H   14    0.18505      0.00000     0.81320    0.00176     0.81495
      H   15    0.18230      0.00000     0.81513    0.00258     0.81770
      H   16    0.18247      0.00000     0.81481    0.00272     0.81753
      H   17    0.19722      0.00000     0.79866    0.00412     0.80278
      H   18    0.18039      0.00000     0.81609    0.00352     0.81961
      H   19    0.15060      0.00000     0.84324    0.00616     0.84940
      H   20    0.18538      0.00000     0.81185    0.00277     0.81462
      H   21    0.21377      0.00000     0.78403    0.00220     0.78623
      H   22    0.19472      0.00000     0.80072    0.00455     0.80528
      H   23    0.21725      0.00000     0.78202    0.00073     0.78275
      H   24    0.21848      0.00000     0.77969    0.00183     0.78152
      H   25    0.17897      0.00000     0.81525    0.00579     0.82103
      H   26    0.21934      0.00000     0.77989    0.00077     0.78066
=======================================================================
   * Total *    0.00000     21.98875    59.82038    0.19087    82.00000

sobereva

1 你做了多余的步骤，你把密度差给转化成了高斯函数分布描述。不打算以这种方式讨论就不要做这个额外步骤。
普通的密度差图做法仔细效仿此文
使用Multiwfn作电子密度差图
http://sobereva.com/113

风玲儿

冰释之川 发表于 2017-5-28 11:26
通过计算基态和垂直激发态下的片段NPA电荷，最后求差得到转移量。
高斯里单点加上 pop=NPA关键词可以获得N ...

谢谢老师，获得基态和激发态的ＮＰＡ电荷输入文件用相同的分子坐标吗，还是个自优化后取最不同的坐标？

冰释之川

风玲儿 发表于 2017-5-29 20:04
谢谢老师，获得基态和激发态的ＮＰＡ电荷输入文件用相同的分子坐标吗，还是个自优化后取最不同的坐标？

必然是相同坐标的，都是垂直激发的，不然你做不出EDD

风玲儿

冰释之川 发表于 2017-5-29 20:37
必然是相同坐标的，都是垂直激发的，不然你做不出EDD

哦哦，明白了，谢谢老师