应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » NBO/AIM » 怎么样用choose把CH3NH2分成两个unit
51/1返回列表
查看: 1403  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

kathy2008

木虫 (正式写手)

[求助] 怎么样用choose把CH3NH2分成两个unit

在Gamess下编译NBO后，怎么样用choose把CH3NH2分成两个unit。其中一个单位是CH3，另一个单位是NH2，我的输入文件如下。
$CONTRL SCFTYP=RHF RUNTYP=ENERGY MULT=1 ICHARG=0 MAXIT=200
  COORD=UNIQUE UNITS=ANGS NZVAR=0 NPRINT=7 $END
$SYSTEM MEMORY=100000000 $END
$SCF CONV=1.0d-8  $END
$BASIS  GBASIS=N21 NGAUSS=3 $END
$NBO $END
$DATA
ch3nh2
C1
C  6                0.00000000    0.00000000    0.00000000
N  7                0.00000000    0.00000000    1.47000000
H  1                1.02766199    0.00000000   -0.36333296
H  1               -0.51383099   -0.88998139   -0.36333296
H  1               -0.51383099    0.88998139   -0.36333296
H  1                0.47611863    0.82466165    1.80666632
H  1                0.47611863   -0.82466165    1.80666632
$END
$CHOOSE
  LONE 2 1  END
  BOND S 1 2
       S 1 3
       S 1 4
       S 1 5
       S 2 6
       S 2 7 END
$END
输出的结果总是
NATURAL BOND ORBITALS (Summary):

                                                     Principal Delocalizations
           NBO                 Occupancy    Energy   (geminal,vicinal,remote)
===============================================================================
Molecular unit  1  (CH5N)
   1. BD ( 1) C  1- N  2        1.99858    -0.89908
   2. BD ( 1) C  1- H  3        1.99860    -0.69181  14(v)
   3. BD ( 1) C  1- H  4        1.99399    -0.68892  27(v),26(g)
   4. BD ( 1) C  1- H  5        1.99399    -0.68892  28(v),25(g)
      Donor NBO (i)              Acceptor NBO (j)       kcal/mol   a.u.    a.u.
===============================================================================

within unit  1
   2. BD ( 1) C 1- H 3       14. RY*( 1) N 2              0.84    2.18    0.038
   3. BD ( 1) C 1- H 4       26. BD*( 1) C 1- H 5         0.52    1.39    0.024
   3. BD ( 1) C 1- H 4       27. BD*( 1) N 2- H 6         3.03    1.37    0.057
   4. BD ( 1) C 1- H 5       25. BD*( 1) C 1- H 4         0.52    1.39    0.024
   4. BD ( 1) C 1- H 5       28. BD*( 1) N 2- H 7         3.03    1.37    0.057
   5. BD ( 1) N 2- H 6       10. RY*( 1) C 1              0.56    1.78    0.028
   5. BD ( 1) N 2- H 6       25. BD*( 1) C 1- H 4         2.85    1.51    0.059
   6. BD ( 1) N 2- H 7       10. RY*( 1) C 1              0.56    1.78    0.028
   6. BD ( 1) N 2- H 7       26. BD*( 1) C 1- H 5         2.85    1.51    0.059
   7. CR ( 1) C 1            16. RY*( 3) N 2              0.61   13.11    0.080
   7. CR ( 1) C 1            18. RY*( 1) H 3              1.40   11.99    0.116
   7. CR ( 1) C 1            19. RY*( 1) H 4              1.55   11.99    0.122
   7. CR ( 1) C 1            20. RY*( 1) H 5              1.55   11.99    0.122
   8. CR ( 1) N 2            10. RY*( 1) C 1              1.51   16.23    0.140
   8. CR ( 1) N 2            12. RY*( 3) C 1              0.84   16.77    0.106
   8. CR ( 1) N 2            21. RY*( 1) H 6              0.61   16.26    0.089
   8. CR ( 1) N 2            22. RY*( 1) H 7              0.61   16.26    0.089
   9. LP ( 1) N 2            24. BD*( 1) C 1- H 3         8.13    1.13    0.086
   9. LP ( 1) N 2            25. BD*( 1) C 1- H 4         1.46    1.14    0.037
   9. LP ( 1) N 2            26. BD*( 1) C 1- H 5         1.46    1.14    0.037
也就是说没有分成两个unit，不知那个同行有何高见。谢谢。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-03-13 13:02:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kathy2008

木虫 (正式写手)
引用回帖:
楼: Originally posted by sobereva at 2012-03-13 21:49:17:
有些NBO版本在CHOOSE这部分是有bug的，即便是闭壳层也必须写ALPHA...END段落，否则不认。
如果你的情况不写ALPHA...END就已经可以识别而写的话反倒不识别，就没必要写它。
LP可以不明确地写出来，会自动搜索。

非常感谢sobereva老师。去掉Lone描述后，成功了。
赞一下 回复此楼
8楼2012-03-13 22:19:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
gmy1990(金币+2): 2012-03-13 19:23:35
kathy2008(金币+5): ★有帮助 2012-03-13 20:00:09
kathy2008: 金币+15, ★★★很有帮助 2012-04-05 16:37:08
CHOOSE段落写成
$CHOOSE
ALPHA
BOND S 1 3
        S 1 4
        S 1 5
        S 2 6
        S 2 7
END
END
$END
这样由于1 2之间没有了BD，所以CH3和NH2就不会被认为是一个unit了。
赞一下 回复此楼
2楼2012-03-13 17:20:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vigaryang

木虫 (正式写手)
请教sobereva老师：不需要在LP部分对N的孤对电子进行指定么？
赞一下 回复此楼
3楼2012-03-13 19:43:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kathy2008

木虫 (正式写手)
引用回帖:
楼: Originally posted by sobereva at 2012-03-13 17:20:38:
CHOOSE段落写成
$CHOOSE
ALPHA
BOND S 1 3
        S 1 4
        S 1 5
        S 2 6
        S 2 7
END
END
$END
这样由于1 2之间没有了BD，所以CH3和NH2就不会被认为是一个unit了。

谢谢回复，但运行后，提示
NATURAL BOND ORBITAL ANALYSIS:
Reading $CHOOSE input list.

Error in input of bond orbitals:
Keyword for orbital type is not LONE, BOND, or 3CBOND (read `ALPHA ')
赞一下 回复此楼
4楼2012-03-13 19:58:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 302求调剂 +7 易！? 2026-04-13 7/350 2026-04-13 10:59 by Faiz5552
[考研] 0854调剂 +10 长弓傲 2026-04-11 11/550 2026-04-13 10:38 by wp06
[考研] 297求调剂 +18 ORCHID1 2026-04-10 19/950 2026-04-13 09:21 by lhj2009
[考研] 一志愿华工085600 331分 +7 天下ww 2026-04-09 7/350 2026-04-13 09:01 by lhj2009
[考研] 本科郑州大学，一志愿华东师范大学282求调剂 +24 熊哥xtk 2026-04-07 27/1350 2026-04-13 08:30 by lhj2009
[考研] 0854调剂 +12 长弓傲 2026-04-09 13/650 2026-04-12 09:56 by 逆水乘风
[找工作] 山东高校教师考核超级无底线，员工过不下去啦 +4 qut2026 2026-04-09 9/450 2026-04-12 00:54 by qut2026
[考研] 一志愿211，0703化学305分求调剂 +26 严西西戏 2026-04-06 33/1650 2026-04-11 23:01 by 314126402
[考研] 086003调剂求助 +21 苏弋万 2026-04-09 22/1100 2026-04-11 20:25 by dongdian1
[考研] 283求调剂 +22 那个噜子 2026-04-09 22/1100 2026-04-11 10:41 by 逆水乘风
[考研] 工科273调剂 +6 X1999 2026-04-09 7/350 2026-04-11 10:23 by zhq0425
[考研] 一志愿211，化学310分，本科重点双非，求调剂 +23 努力奋斗112 2026-04-08 23/1150 2026-04-10 23:29 by 314126402
[考研] 调剂 +19 小张ZA 2026-04-10 20/1000 2026-04-10 22:08 by 猪会飞
[考研] 265求调剂 +12 风说她早忘了 2026-04-10 13/650 2026-04-10 18:56 by chemisry
[考研] 材料专业344求调剂 +16 hualkop 2026-04-10 21/1050 2026-04-10 17:28 by laoshidan
[论文投稿] mdpi小修rvr时间四五天了 20+3 哈哈high 2026-04-08 5/250 2026-04-10 16:02 by 北京莱茵润色
[基金申请] 有爆料，一个青年教师卖房得400万，然后换了一个四青帽子 +9 babu2015 2026-04-08 9/450 2026-04-10 11:43 by 苏东坡二世
[考研] 材料调剂 +10 18815505510 2026-04-09 11/550 2026-04-09 17:07 by 544594351
[考研] 388求调剂 +6 四川王涛 2026-04-07 8/400 2026-04-08 00:17 by JourneyLucky
[考研] 331求调剂 +5 张元一 2026-04-07 6/300 2026-04-07 22:13 by hemengdong
信息提示
关闭
请填处理意见
关闭