24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2139  |  回复: 7

科学岛

铁虫 (初入文坛)

[求助] 用gaussian计算负一价离子的NBO时,输出文件中电荷显示为何为0

Sample Text
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lishijunzong

铁杆木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
贴出你的输入与输出文件
2楼2013-07-03 08:38:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

科学岛

铁虫 (初入文坛)

引用回帖:
2楼: Originally posted by lishijunzong at 2013-07-03 08:38:13
贴出你的输入与输出文件

输出文件中共有3处Summary of Natural Population Analysis:,前两处显示total=-1.000,但最后一处却显示total=0.不知是和原因,上网查过,加上(7/76=8)后仍然是同样的错误。简单的输入输出如下:
输入:%nproc=4     
%mem=500MB
# pbe1pbe/genecp Pop=NBORead

输出:前两处中的一处, Summary of Natural Population Analysis:                  
                                                         
                                       Natural Population
                Natural  -----------------------------------------------
    Atom  No    Charge         Core      Valence    Rydberg      Total
-----------------------------------------------------------------------
     Au    1   -0.06858     67.99088    11.06733    0.01038    79.06858
     Au    2   -0.06858     67.99088    11.06733    0.01038    79.06858
      C    3   -0.86284      1.99974     4.80188    0.06123     6.86284
=======================================================================
   * Total *   -1.00000    137.98149    26.93653    0.08198   165.00000
Summary of Natural Population Analysis:                  
                                                         
                                       Natural Population
                Natural  -----------------------------------------------
    Atom  No    Charge         Core      Valence    Rydberg      Total
-----------------------------------------------------------------------
     Au    1   -0.02329     33.99503     5.52360    0.00466    39.52329
     Au    2   -0.02329     33.99503     5.52360    0.00466    39.52329
      C    3    0.04659      0.99986     1.92479    0.02877     2.95341
=======================================================================
   * Total *    0.00000     68.98993    12.97199    0.03808    82.00000

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : reopt_reopt_Anion_Geom_31_.out
  • 2013-07-03 08:53:07, 204.33 K
  • 附件 2 : reopt_reopt_Anion_Geom_31_.gjf
  • 2013-07-03 08:53:22, 8.36 K
3楼2013-07-03 08:53:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

科学岛

铁虫 (初入文坛)

引用回帖:
3楼: Originally posted by 科学岛 at 2013-07-03 08:53:39
输出文件中共有3处Summary of Natural Population Analysis:,前两处显示total=-1.000,但最后一处却显示total=0.不知是和原因,上网查过,加上(7/76=8)后仍然是同样的错误。简单的输入输出如下:
输入:%nproc ...

请大牛们指点一下啊,万分感谢!
4楼2013-07-03 10:54:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫


【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
gkf高: 金币+3 2013-07-03 15:06:11
科学岛: 金币+5, ★★★★★最佳答案 2013-07-03 19:17:15
第一遍输出的是对体系的总密度分析,所以加和为-1。第二遍是对alpha电子进行分析,第三遍是对beta电子进行分析。由于你的情况是自旋多重度为2的开壳层,所以alpha电子比beta电子多1,而NBO在计算alpha和beta部分的原子电荷都时将原子核电荷除以2,故beta电子和核电荷正好抵消了,而alpha电子比核电荷多1,故alpha的净电荷为-1。
5楼2013-07-03 14:49:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

科学岛

铁虫 (初入文坛)

引用回帖:
5楼: Originally posted by sobereva at 2013-07-03 14:49:27
第一遍输出的是对体系的总密度分析,所以加和为-1。第二遍是对alpha电子进行分析,第三遍是对beta电子进行分析。由于你的情况是自旋多重度为2的开壳层,所以alpha电子比beta电子多1,而NBO在计算alpha和beta部分的原 ...

那比如说分析Au2C_阴离子团簇时,能否将Au2_看做一个整体,而分析C和Au2_之间的电荷转移情况。
                                      Natural Population                 
               Natural  -----------------------------------------------  
   Atom  No    Charge         Core      Valence    Rydberg      Total   
-----------------------------------------------------------------------  
    Au    1   -0.06858     67.99088    11.06733    0.01038    79.06858   
    Au    2   -0.06858     67.99088    11.06733    0.01038    79.06858   
     C    3   -0.86284      1.99974     4.80188    0.06123     6.86284   
=======================================================================  
  * Total *   -1.00000    137.98149    26.93653    0.08198   165.00000   
在上表中,是不是就是有0.86284个电荷由Au2_转移到C上。
有一篇文献中有这样一句话:The calculated nonorthogonal natural atomic overlap populations between Au and Au are larger than those between Au and Zn. 这是不是就要看输出文件中的体系总密度分析后的
Atom-atom overlap-weighted NAO bond order:                        
                                                                  
     Atom    1       2       3                                    
     ---- ------  ------  ------                                   
   1. Au  0.0000  0.0851  0.8473                                   
   2. Au  0.0851  0.0000  0.8473                                   
   3.  C  0.8473  0.8473  0.0000                                   
  部分。                                                                 
                                                                  
还是下面这部分
Atom-atom overlap-weighted NAO bond order, Totals by atom:        
                                                                  
     Atom    1                                                     
     ---- ------                                                   
   1. Au  0.9325                                                   
   2. Au  0.9325                                                   
   3.  C  1.6947                                                   
请指点一下啊,谢谢!
6楼2013-07-03 19:40:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫


你的理解正确。

"Atom-atom overlap-weighted NAO bond order, Totals by atom"和那句话无关。"Atom-atom overlap-weighted NAO bond order"其实也不是那句话所精确对应的,只能说是相似(结论估计能够定性一致),因为其中用的是NAO的密度矩阵而非PNAO的密度矩阵。关于NBO输出的这个东西的含义和计算方式见《谈谈NBO中用bndidx关键词输出的三种键级》
http://hi.baidu.com/sobereva/item/1fb69cd491e5e3ba32db90c3)。要想精确算那个量就用DMPNAO打印出PNAO下的密度矩阵,然后和相应的PNAO重叠矩阵元相乘来得到。

至于那文献作者怎么讨论的这个量这不好说,很有可能作者自己不懂十分清楚理论,误用Atom-atom overlap-weighted NAO bond order来讨论这个量。
7楼2013-07-03 21:17:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

科学岛

铁虫 (初入文坛)

引用回帖:
7楼: Originally posted by sobereva at 2013-07-03 21:17:08
你的理解正确。

"Atom-atom overlap-weighted NAO bond order, Totals by atom"和那句话无关。"Atom-atom overlap-weighted NAO bond order"其实也不是那句话所精确对应的,只能说是相似(结 ...

非常感谢您。
8楼2013-07-04 08:08:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 科学岛 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +4 电气小神童 2026-04-04 6/300 2026-04-07 00:14 by guanxin1001
[考研] 304求调剂 +7 素年祭语 2026-04-06 10/500 2026-04-06 23:31 by 素年祭语
[考研] 305求调剂 +3 77Qi 2026-04-06 3/150 2026-04-06 15:36 by going home
[考研] 调剂 一志愿吉林大学357分 +5 .Starry. 2026-04-04 5/250 2026-04-06 09:28 by cql1109
[考研] 一志愿211,化学学硕,310分,本科重点双非,求调剂 +13 努力奋斗112 2026-04-04 13/650 2026-04-06 07:13 by jj987
[考研] 求调剂求调剂 +8 121. 2026-04-02 8/400 2026-04-05 20:15 by lys0704
[考研] 282求调剂 +3 aaa车辆 2026-04-01 3/150 2026-04-05 17:03 by yulian1987
[考研] 0832食品科学与工程学硕282调剂 +6 鱼在水中游a 2026-04-02 9/450 2026-04-05 11:45 by flysky1234
[考研] 材料求调剂 +10 呢呢妮妮 2026-04-01 10/500 2026-04-04 23:12 by 无际的草原
[考研] 278求调剂 +14 范婷娜 2026-04-04 15/750 2026-04-04 22:15 by lqwchd
[考研] 材料调剂 +18 一样YWY 2026-04-02 19/950 2026-04-04 22:14 by hemengdong
[考研] 一志愿东北大学085901土木专硕345求调剂 +3 zxt11111 2026-04-04 3/150 2026-04-04 14:21 by 土木硕士招生
[考研] 一志愿重庆大学085404,总分314分,求调剂 +4 zf83hn 2026-04-03 4/200 2026-04-03 21:25 by 啵啵啵0119
[考研] 295求调剂 +3 尚偌呀 2026-04-03 4/200 2026-04-03 21:23 by zhq0425
[考研] 081200-11408-276学硕求调剂 +5 崔wj 2026-04-03 5/250 2026-04-03 15:06 by arrow8852
[考研] 一志愿深大085601材料工程专业(专硕)300分可以调剂去哪 +8 10160315 2026-04-02 8/400 2026-04-03 09:36 by hypershenger
[考研] 26考研调剂 +4 Wnz.20030617 2026-04-01 5/250 2026-04-02 16:11 by 1939136013狗壮
[考研] 0710生物学求调剂 +9 manman511 2026-04-01 9/450 2026-04-02 10:00 by zxl830724
[考研] 303分 0807学硕求调剂 +3 TYC3632 2026-04-01 3/150 2026-04-01 19:24 by lwk2004
[考研] 一志愿 南京航空航天大学 ,080500材料科学与工程学硕 +10 @taotao 2026-03-31 11/550 2026-04-01 09:43 by xiayizhi
信息提示
请填处理意见