计算化学公社

标题: muitiwfn使用和NBO分析 [打印本页]

作者
Author: 风玲儿 时间: 2017-5-28 09:22
标题: muitiwfn使用和NBO分析
我计算了从基态到激发态的Electron density difference ，得到图一的结果，查文献有得到图二的结果，感觉自己的别扭，应如何修改成图二的视图呢？

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

(, 下载次数 Times of downloads: 20) (, 下载次数 Times of downloads: 21)
图一图二

作者
Author: 冰释之川 时间: 2017-5-28 11:26
本帖最后由冰释之川于 2017-5-28 11:27 编辑

通过计算基态和垂直激发态下的片段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

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

作者
Author: 风玲儿 时间: 2017-5-29 20:04

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

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

作者
Author: 冰释之川 时间: 2017-5-29 20:37

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

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

作者
Author: 风玲儿 时间: 2017-5-29 20:39

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

哦哦，明白了，谢谢老师

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