| 查看: 1411 | 回复: 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，不知那个同行有何高见。谢谢。 |
回复此楼» 猜你喜欢
调剂招生
已经有0人回复
-
调剂招生
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有250人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有1人回复
-
实际气体内压强的微观机制或将改写--------从一个基本的力学佯谬说起
已经有1人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有5人回复
-
招收0703开头的学硕调剂生 线上复试
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于两个offer分两批申报的问题~~不要拍砖
已经有25人回复
- 怎么用柱色谱分开这两个东西。
已经有27人回复
- NiO/Al2O3的H2-TPR峰高于NiO，MoO3/Al2O3的H2-TPR峰低于MoO3，求原因
已经有6人回复
- CeO2 H2-TPR问题
已经有14人回复
- 如何得到Ar和H2的混合气体
已经有7人回复
- 【求助】H2 Binding Energy@SIESTA
已经有13人回复
- 【求助】一甲基胺盐酸盐如何制备成一甲胺CH3NH2的THF溶液或者甲醇溶液
已经有16人回复
- 【求助】请教两个问题，一个关于推荐信，一个关于套磁的人选
已经有4人回复
- 【求助】请教C6H2、C8H2是什么物质？
已经有4人回复
kathy2008
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2161.5
- 红花: 2
- 帖子: 328
- 在线: 213.2小时
- 虫号: 558502
- 注册: 2008-05-15
- 性别: GG
- 专业: 材料物理化学
|
即使对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
sobereva
至尊木虫 (著名写手)
本人已永久离开小木虫
- QC强帖: 17
- 应助: 226 (大学生)
- 贵宾: 5.4
- 金币: 2180.5
- 红花: 217
- 帖子: 1241
- 在线: 366.8小时
- 虫号: 281271
- 注册: 2006-09-23
- 专业: 宇宙学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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
感谢参与，应助指数 +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
vigaryang
木虫 (正式写手)
- 应助: 16 (小学生)
- 金币: 3466.3
- 散金: 677
- 红花: 6
- 帖子: 471
- 在线: 495.7小时
- 虫号: 227097
- 注册: 2006-03-23
- 性别: GG
- 专业: 应用化学
3楼2012-03-13 19:43:40
kathy2008
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2161.5
- 红花: 2
- 帖子: 328
- 在线: 213.2小时
- 虫号: 558502
- 注册: 2008-05-15
- 性别: GG
- 专业: 材料物理化学
4楼2012-03-13 19:58:15
