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

24小时热门版块排行榜    

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

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

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

kathy2008

木虫 (正式写手)
引用回帖:
3楼: Originally posted by vigaryang at 2012-03-13 19:43:40:
请教sobereva老师：不需要在LP部分对N的孤对电子进行指定么？

即使对N原子的LP部分进行指定，运行后，出错信息如下
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 ')
80000  WORDS OF DYNAMIC MEMORY USED
EXECUTION OF GAMESS TERMINATED NORMALLY Tue Mar 13 21:12:02 2012
问题的关键应该是ALPHA针对的是open-shell。
赞一下 回复此楼
6楼2012-03-13 20:14:31
已阅   回复此楼   关注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可以不明确地写出来，会自动搜索。
赞一下(2人) 回复此楼
7楼2012-03-13 21:49:17
已阅   回复此楼   关注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的回帖
相关版块跳转 我要订阅楼主 kathy2008 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 306求调剂 +9 chuanzhu川烛 2026-03-18 9/450 2026-03-23 13:17 by luoyongfeng
[考研] 工科0856求调剂 +4 沐析汀汀 2026-03-21 4/200 2026-03-23 13:08 by 醉在风里
[考研] 工科材料085601 279求调剂 +8 困于星晨 2026-03-17 10/500 2026-03-23 13:05 by 醉在风里
[考研] 276求调剂。有半年电池和半年高分子实习经历 +8 材料学257求调剂 2026-03-23 9/450 2026-03-23 13:01 by ztnimte
[考研] 一志愿东华大学化学070300，求调剂 +7 2117205181 2026-03-21 8/400 2026-03-22 22:55 by chixmc
[考研] 303求调剂 +5 安忆灵 2026-03-22 6/300 2026-03-22 12:46 by 素颜倾城1988
[考研] 285求调剂 +6 ytter 2026-03-22 6/300 2026-03-22 12:09 by 星空星月
[考研] 材料工程专硕 348分求调剂 +3 冬辞. 2026-03-17 5/250 2026-03-21 18:47 by 学员8dgXkO
[考研] 266求调剂 +3 哇呼哼呼哼 2026-03-20 3/150 2026-03-21 16:46 by barlinike
[考研] 316求调剂 +6 梁茜雯 2026-03-19 6/300 2026-03-21 06:32 by Ecowxq666！
[考研] 299求调剂 +6 △小透明* 2026-03-17 6/300 2026-03-21 02:42 by JourneyLucky
[考研] 求调剂 +3 Ma_xt 2026-03-17 3/150 2026-03-21 02:05 by JourneyLucky
[考研] 085700资源与环境308求调剂 +12 墨墨漠 2026-03-18 13/650 2026-03-21 01:42 by JourneyLucky
[考研] 274求调剂 +10 S.H1 2026-03-18 10/500 2026-03-20 23:51 by JourneyLucky
[考研] 295求调剂 +4 一志愿京区211 2026-03-18 6/300 2026-03-20 23:41 by JourneyLucky
[考研] 一志愿南理工085701环境302求调剂院校 +3 葵梓卫队 2026-03-20 3/150 2026-03-20 19:28 by zhukairuo
[考研] 生物学调剂招人！！！ +3 山海天岚 2026-03-17 4/200 2026-03-19 21:34 by 怎么释怀
[考研] 0703化学调剂 +5 pupcoco 2026-03-17 8/400 2026-03-19 13:58 by houyaoxu
[考研] 一志愿福大288有机化学，求调剂 +3 小木虫200408204 2026-03-18 3/150 2026-03-19 13:31 by houyaoxu
[考研] 一志愿南京大学，080500材料科学与工程，调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
信息提示
关闭
请填处理意见
关闭