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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » NBO/AIM » 怎么样用choose把CH3NH2分成两个unit
51/1返回列表
查看: 1363  |  回复: 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的回帖

vigaryang

木虫 (正式写手)
请教sobereva老师：不需要在LP部分对N的孤对电子进行指定么？
赞一下 回复此楼
3楼2012-03-13 19:43:40
已阅   回复此楼   关注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的回帖

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的回帖

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的回帖
查看全部 8 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0854 考研调剂 招生了！AI 方向 +3 pk3725069 2026-03-19 14/700 2026-03-23 14:44 by Jack_333
[考研] 336求调剂 +3 收到VS 2026-03-20 3/150 2026-03-23 13:23 by Iveryant
[考研] 306求0703调剂一志愿华中师范 +6 纸鱼ly 2026-03-21 6/300 2026-03-23 12:08 by Iveryant
[考研] 一志愿武理材料工程348求调剂 +5 ￣^￣゜汗 2026-03-19 7/350 2026-03-22 19:44 by 公瑾逍遥
[考研] 求调剂 +6 十三加油 2026-03-21 6/300 2026-03-22 17:00 by i_cooler
[考研] 环境学硕288求调剂 +6 皮皮皮123456 2026-03-22 6/300 2026-03-22 16:52 by i_cooler
[考研] 298求调剂一志愿211 +3 上岸6666@ 2026-03-20 3/150 2026-03-22 15:50 by ColorlessPI
[考研] 285求调剂 +6 ytter 2026-03-22 6/300 2026-03-22 12:09 by 星空星月
[考研] 一志愿华中科技大学071000，求调剂 +4 沿岸有贝壳6 2026-03-21 4/200 2026-03-22 07:21 by ilovexiaobin
[考研] 一志愿东华大学控制学硕320求调剂 +3 Grand777 2026-03-21 3/150 2026-03-21 19:23 by 简之-
[考研] 332求调剂 +3 凤凰院丁真 2026-03-20 3/150 2026-03-21 10:27 by luoyongfeng
[考研] 一志愿中海洋材料工程专硕330分求调剂 +8 小材化本科 2026-03-18 8/400 2026-03-20 23:16 by JourneyLucky
[考研] 求调剂，一志愿:南京航空航天大学大学 ，080500材料科学与工程学硕，总分289分 +4 @taotao 2026-03-19 4/200 2026-03-20 22:14 by JourneyLucky
[考研] 求调剂一志愿南京航空航天大学289分 +3 @taotao 2026-03-19 3/150 2026-03-20 21:34 by JourneyLucky
[考研] 求调剂 +3 @taotao 2026-03-20 3/150 2026-03-20 19:35 by JourneyLucky
[考研] 0703化学调剂 +5 pupcoco 2026-03-17 8/400 2026-03-19 13:58 by houyaoxu
[考研] 085600材料与化工求调剂 +6 绪幸与子 2026-03-17 6/300 2026-03-19 13:27 by houyaoxu
[考研] 材料考研调剂 +3 xwt。 2026-03-19 3/150 2026-03-19 11:22 by w沐阳w
[考研] 材料，纺织，生物（0856、0710），化学招生啦 +3 Eember. 2026-03-17 9/450 2026-03-18 10:28 by Eember.
[考研] 考研调剂 +3 淇ya_~ 2026-03-17 5/250 2026-03-17 09:25 by Winj1e
信息提示
关闭
请填处理意见
关闭