24小时热门版块排行榜    

查看: 2393  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 科学岛 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +37 医学老男孩 2026-08-13 86/4300 2026-08-18 00:18 by longsantaizi
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 时间戳变了,能看出什么问题? +4 基诺咪客 2026-08-17 4/200 2026-08-17 16:10 by Vivilian
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +12 Ldrop2023 2026-08-13 15/750 2026-08-17 15:02 by 小豌豆_发芽
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
信息提示
请填处理意见