24小时热门版块排行榜    

查看: 1364  |  回复: 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的回帖

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 个回答

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了。

ALPHA在NBO说明书上,针对的是open-shell。也许对这样的中性分子不适合。
5楼2012-03-13 20:08:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 急发核心期刊论文 +3 贤达问津 2026-03-23 5/250 2026-03-23 17:13 by 妹子不好惹
[考研] 一志愿重庆大学085700资源与环境,总分308求调剂 +6 墨墨漠 2026-03-23 7/350 2026-03-23 16:09 by 墨墨漠
[考研] 081700 调剂 267分 +5 迷人的哈哈 2026-03-23 5/250 2026-03-23 14:59 by lbsjt
[考研] 0854电子信息求调剂 324 +3 Promise-jyl 2026-03-23 3/150 2026-03-23 13:43 by wangkm
[考研] 求老师收我 +3 zzh16938784 2026-03-23 3/150 2026-03-23 12:56 by ztnimte
[考研] 北科281学硕材料求调剂 +8 tcxiaoxx 2026-03-20 9/450 2026-03-23 12:16 by tcxiaoxx
[考研] 306求0703调剂一志愿华中师范 +6 纸鱼ly 2026-03-21 6/300 2026-03-23 12:08 by Iveryant
[考研] 一志愿北京化工大学070300 学硕336求调剂 +5 vv迷 2026-03-21 8/400 2026-03-22 14:20 by ColorlessPI
[考研] 0703化学调剂 +4 妮妮ninicgb 2026-03-21 4/200 2026-03-21 18:39 by 学员8dgXkO
[考研] 求调剂 +6 Mqqqqqq 2026-03-19 6/300 2026-03-21 08:04 by JourneyLucky
[考研] 二本跨考郑大材料306英一数二 +3 z1z2z3879 2026-03-17 3/150 2026-03-21 02:29 by JourneyLucky
[考研] 317求调剂 +5 申子申申 2026-03-19 9/450 2026-03-20 22:26 by JourneyLucky
[考研] 材料与化工 322求调剂 +4 然11 2026-03-19 4/200 2026-03-20 22:12 by luoyongfeng
[考研] 一志愿苏州大学材料求调剂,总分315(英一) +5 sbdksD 2026-03-19 5/250 2026-03-20 22:10 by luoyongfeng
[考研] 261求B区调剂,科研经历丰富 +3 牛奶很忙 2026-03-20 4/200 2026-03-20 19:34 by JourneyLucky
[考研] 材料学硕318求调剂 +5 February_Feb 2026-03-19 5/250 2026-03-19 23:51 by 23Postgrad
[考研] 材料考研调剂 +3 xwt。 2026-03-19 3/150 2026-03-19 11:22 by w沐阳w
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
[硕博家园] 湖北工业大学 生命科学与健康学院-课题组招收2026级食品/生物方向硕士 +3 1喜春8 2026-03-17 5/250 2026-03-17 17:18 by ber川cool子
[考研] 有没有道铁/土木的想调剂南林,给自己招师弟中~ +3 TqlXswl 2026-03-16 7/350 2026-03-17 15:23 by TqlXswl
信息提示
请填处理意见