24小时热门版块排行榜    

查看: 1360  |  回复: 7
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

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,不知那个同行有何高见。谢谢。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫


★ ★
zhou2009(金币+2): 2012-03-14 16:50:10
引用回帖:
6楼: Originally posted by kathy2008 at 2012-03-13 20:14:31:
即使对N原子的LP部分进行指定,运行后,出错信息如下
NATURAL BOND ORBITAL ANALYSIS:
Reading $CHOOSE input list.

Error in input of bond orbitals:
Keyword for orbital type is not LONE, BOND, ...

有些NBO版本在CHOOSE这部分是有bug的,即便是闭壳层也必须写ALPHA...END段落,否则不认。
如果你的情况不写ALPHA...END就已经可以识别而写的话反倒不识别,就没必要写它。
LP可以不明确地写出来,会自动搜索。
7楼2012-03-13 21:49:17
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 招08考数学 +6 laoshidan 2026-03-20 14/700 2026-03-23 14:37 by 15614359529
[考研] 工科材料085601 279求调剂 +8 困于星晨 2026-03-17 10/500 2026-03-23 13:05 by 醉在风里
[考研] 333求调剂 +6 87639 2026-03-21 10/500 2026-03-23 10:41 by Iveryant
[考研] 070300化学求调剂 +5 苑豆豆 2026-03-20 5/250 2026-03-23 09:39 by guoliang1816
[考研] 一志愿武理材料工程348求调剂 +5  ̄^ ̄゜汗 2026-03-19 7/350 2026-03-22 19:44 by 公瑾逍遥
[考研] 08工科 320总分 求调剂 +11 梨花珞晚风 2026-03-17 11/550 2026-03-22 17:42 by luoyongfeng
[考研] 324求调剂 +6 lucky呀呀呀鸭 2026-03-20 6/300 2026-03-22 16:01 by ColorlessPI
[考研] 269专硕求调剂 +6 金恩贝 2026-03-21 6/300 2026-03-22 14:31 by ColorlessPI
[考研] 一志愿华中科技大学071000,求调剂 +4 沿岸有贝壳6 2026-03-21 4/200 2026-03-22 07:21 by ilovexiaobin
[考研] 材料求调剂 +5 @taotao 2026-03-21 5/250 2026-03-21 20:55 by lbsjt
[考研] 0703化学297求调剂 +3 Daisy☆ 2026-03-20 3/150 2026-03-21 17:45 by ColorlessPI
[考研] 一志愿重庆大学085700资源与环境总分308求调剂 +7 墨墨漠 2026-03-20 7/350 2026-03-21 16:36 by barlinike
[考研] 化学求调剂 +4 临泽境llllll 2026-03-17 5/250 2026-03-21 02:23 by JourneyLucky
[考研] 一志愿武汉理工材料工程专硕调剂 +9 Doleres 2026-03-19 9/450 2026-03-20 22:36 by JourneyLucky
[考研] 求调剂 +3 @taotao 2026-03-20 3/150 2026-03-20 19:35 by JourneyLucky
[考研] 085410人工智能专硕317求调剂(0854都可以) +4 xbxudjdn 2026-03-18 4/200 2026-03-20 09:07 by 不168
[考研] 085601材料工程专硕求调剂 +10 慕寒mio 2026-03-16 10/500 2026-03-19 15:26 by 丁丁*
[考研] 085600材料与化工求调剂 +6 绪幸与子 2026-03-17 6/300 2026-03-19 13:27 by houyaoxu
[考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 283求调剂 +3 听风就是雨; 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
信息提示
请填处理意见