请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 196|回复: 5

[Multiwfn使用咨询] muitiwfn使用和NBO分析

[复制链接]

20

帖子

0

威望

41

eV
积分
61

Level 2 能力者

发表于 2017-5-28 09:22:07 | 显示全部楼层 |阅读模式
我计算 了从基态到激发态的Electron density difference  ,得到图一的结果,查文献有得到图二的结果,感觉自己的别扭,应如何修改成图二的视图呢?

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






图一

图一

图二

图二

        图一                                                                      图二

514

帖子

2

威望

1650

eV
积分
2204

Level 5 (御坂)

发表于 2017-5-28 11:26:53 | 显示全部楼层
本帖最后由 冰释之川 于 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


Stand on the shoulders of giants

9493

帖子

16

威望

1万

eV
积分
21642

管理员

公社社长

发表于 2017-5-28 12:37:51 | 显示全部楼层
1 你做了多余的步骤,你把密度差给转化成了高斯函数分布描述。不打算以这种方式讨论就不要做这个额外步骤。
普通的密度差图做法仔细效仿此文
使用Multiwfn作电子密度差图
http://sobereva.com/113
思想家公社的门口Blog:http://sobereva.com
北京科音自然科学研究中心:http://www.keinsci.com
Multiwfn主页:http://sobereva.com/multiwfn
计算化学公社:http://bbs.keinsci.com
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论内容相同,可加入任意其一但不可都加入,申请信息必须注明研究方向,否则一概不批。

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

20

帖子

0

威望

41

eV
积分
61

Level 2 能力者

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

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

514

帖子

2

威望

1650

eV
积分
2204

Level 5 (御坂)

发表于 2017-5-29 20:37:20 | 显示全部楼层
风玲儿 发表于 2017-5-29 20:04
谢谢老师,获得基态和激发态的NPA电荷输入文件用相同的分子坐标吗,还是个自优化后取最不同的坐标?

必然是相同坐标的,都是垂直激发的,不然你做不出EDD
Stand on the shoulders of giants

20

帖子

0

威望

41

eV
积分
61

Level 2 能力者

 楼主| 发表于 2017-5-29 20:39:35 | 显示全部楼层
冰释之川 发表于 2017-5-29 20:37
必然是相同坐标的,都是垂直激发的,不然你做不出EDD

哦哦,明白了,谢谢老师
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

Archiver|手机版|小黑屋|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2017-6-29 01:20 , Processed in 0.133170 second(s), 24 queries .

Powered by Discuz! X3.2

© 2001-2013 Comsenz Inc.

快速回复 返回顶部 返回列表