| 查看: 1638 | 回复: 12 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
yuguotianqing木虫 (正式写手)
|
[交流]
【求助】用Scan来做一个多原子的分子的势能面。输入文件怎么建。
|
||
|
我在算一个CH3Ti(O)Cl这个分子的势能面，改变下C-Ti的键长，其它参数都放开。在高斯03中，请问下输入文件怎么建。谢谢。 下面是这个分子的内坐标。 0 1 C H 1 B1 H 1 B2 2 A1 H 1 B3 3 A2 2 D1 Ti 1 B4 4 A3 3 D2 Cl 5 B5 1 A4 4 D3 O 5 B6 1 A5 4 D4 B1 1.09674443 B2 1.09670621 B3 1.09582506 B4 2.08376622 B5 2.25473385 B6 1.60266480 A1 109.87425335 A2 107.34028553 A3 109.47198097 A4 115.58707478 A5 107.70460198 D1 -117.89720170 D2 120.47969812 D3 71.10197770 D4 -62.17865020 |
回复此楼» 猜你喜欢
调剂招生
已经有0人回复
-
调剂招生
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有165人回复
-
调剂招生
已经有1人回复
-
实际气体内压强的微观机制或将改写--------从一个基本的力学佯谬说起
已经有2人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有5人回复
-
招收0703开头的学硕调剂生 线上复试
已经有5人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有37人回复
-
【接测试】TPD-TPR、TPO、脉冲吸附、原位红外，电化学原位红外，站内联系
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
zhangmt
至尊木虫 (著名写手)
我叫MT
- QC强帖: 5
- 应助: 99 (初中生)
- 金币: 6961.8
- 散金: 10406
- 红花: 49
- 帖子: 1761
- 在线: 763.2小时
- 虫号: 880392
- 注册: 2009-10-22
- 性别: GG
- 专业: 理论和计算化学
★ ★ ★ ★ ★
yuguotianqing(金币+4,VIP+0): 11-13 20:54
yjcmwgk(金币+1,VIP+0): 11-13 21:00
yuguotianqing(金币+4,VIP+0): 11-13 20:54
yjcmwgk(金币+1,VIP+0): 11-13 21:00
|
抱歉，忘了一般版本在scan的时候不给优化的（内部版据说有哦） 可以用另一种方法实现，试试这个： windows下跑了4分半，跑了18个点。图形就不贴了。 #p Opt(ModRedundant) hf/sto-3g Title 0 1 C H 1 B1 H 1 B2 2 A1 H 1 B3 3 A2 2 D1 Ti 1 B4 4 A3 3 D2 Cl 5 B5 1 A4 4 D3 O 5 B6 1 A5 4 D4 B1 1.09674443 B2 1.09670621 B3 1.09582506 B4 2.08376622 B5 2.25473385 B6 1.60266480 A1 109.87425335 A2 107.34028553 A3 109.47198097 A4 115.58707478 A5 107.70460198 D1 -117.89720170 D2 120.47969812 D3 71.10197770 D4 -62.17865020 1 5 1.5 S 20 0.1 [ Last edited by zhangmt on 2009-11-13 at 17:57 ] |
6楼2009-11-13 17:41:01
zhangmt
至尊木虫 (著名写手)
我叫MT
- QC强帖: 5
- 应助: 99 (初中生)
- 金币: 6961.8
- 散金: 10406
- 红花: 49
- 帖子: 1761
- 在线: 763.2小时
- 虫号: 880392
- 注册: 2009-10-22
- 性别: GG
- 专业: 理论和计算化学
★ ★ ★ ★ ★
yuguotianqing(金币+5,VIP+0):谢谢 ，但加opt后运行不了的 11-13 17:20
yuguotianqing(金币+5,VIP+0):谢谢 ，但加opt后运行不了的 11-13 17:20
|
前面自己定义chk路径，内存大小，cpu个数，然后自己改算法和基组。 扫描的地方可以写成如下： #p opt hf/sto-3g scan Title 0 1 C H 1 B1 H 1 B2 2 A1 H 1 B3 3 A2 2 D1 Ti 1 B4 4 A3 3 D2 Cl 5 B5 1 A4 4 D3 O 5 B6 1 A5 4 D4 B1 1.09674443 B2 1.09670621 B3 1.09582506 B5 2.25473385 B6 1.60266480 A1 109.87425335 A2 107.34028553 A3 109.47198097 A4 115.58707478 A5 107.70460198 D1 -117.89720170 D2 120.47969812 D3 71.10197770 D4 -62.17865020 B4 1.5 20 0.1 我没使用opt，直接扫了一个莫斯势的曲线，windows下40秒扫完，图形挺好看的。 [ Last edited by zhangmt on 2009-11-13 at 15:16 ] |
2楼2009-11-13 15:12:29
yuguotianqing
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1402.5
- 散金: 228
- 帖子: 371
- 在线: 80.6小时
- 虫号: 509471
- 注册: 2008-02-23
- 性别: GG
- 专业: 有机合成
3楼2009-11-13 17:17:46
forestwolf9291
木虫 (著名写手)
- 应助: 7 (幼儿园)
- 金币: 2702.3
- 散金: 464
- 红花: 4
- 沙发: 1
- 帖子: 1048
- 在线: 289.7小时
- 虫号: 695237
- 注册: 2009-02-02
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
★
yjcmwgk(金币+1,VIP+0): 11-16 12:56
yjcmwgk(金币+1,VIP+0): 11-16 12:56
|
以前做过，好像可以这样。opt也可以，你试试。 #p hf/sto-3g popt Title 0 1 C H 1 B1 H 1 B2 2 A1 H 1 B3 3 A2 2 D1 Ti 1 B4 4 A3 3 D2 Cl 5 B5 1 A4 4 D3 O 5 B6 1 A5 4 D4 B1 1.09674443 B2 1.09670621 B3 1.09582506 B5 2.25473385 B6 1.60266480 A1 109.87425335 A2 107.34028553 A3 109.47198097 A4 115.58707478 A5 107.70460198 D1 -117.89720170 D2 120.47969812 D3 71.10197770 D4 -62.17865020 B4 1.5 S 20 0.1 [ Last edited by forestwolf9291 on 2009-11-13 at 22:17 ] |
4楼2009-11-13 17:27:09
20