24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1235  |  回复: 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
回复此楼
开心工作,快乐生活
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangmt

至尊木虫 (著名写手)

我叫MT

引用回帖:
Originally posted by yuguotianqing at 2009-11-13 20:53:
下面的这个我试了下,是可以的。
但前面讲的scan 那个要是加上优化结构opt,在高斯03中还有方法运行么

用我后面的方法啊,那个一样可以优化的(每个点都优化),为什么非要用scan关键字呢。g03中scan和opt两个关键字不同时使用的。
一群自以为正义凛然的年轻人将一切不能以科学解释的事情定性为封建迷信并大刀阔斧地进行消灭,其实这是修养不足学识浅薄的一种体现,也是可恶的偏执和愚蠢的自以
12楼2009-11-13 22:11:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答

zhangmt

至尊木虫 (著名写手)

我叫MT

★ ★ ★ ★ ★
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuguotianqing

木虫 (正式写手)

我现在计算的时候需要优化的。你加了OPT后的输入文件在高斯中不能运行。
还请指点。

[ Last edited by yuguotianqing on 2009-11-13 at 17:19 ]
开心工作,快乐生活
3楼2009-11-13 17:17:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

forestwolf9291

木虫 (著名写手)


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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见