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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】用Scan来做一个多原子的分子的势能面。输入文件怎么建。
131/2返回列表12下一页
查看: 1340  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

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
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
开心工作,快乐生活
1楼 2009-11-13 14:27:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

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 ]
赞一下(3人) 回复此楼
一群自以为正义凛然的年轻人将一切不能以科学解释的事情定性为封建迷信并大刀阔斧地进行消灭，其实这是修养不足学识浅薄的一种体现，也是可恶的偏执和愚蠢的自以
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 ]
赞一下(10人) 回复此楼
4楼2009-11-13 17:27:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuguotianqing

木虫 (正式写手)
这个文件在高斯中运行几步也出错了，还请你帮我看下
赞一下 回复此楼
开心工作,快乐生活
5楼2009-11-13 17:40:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangmt

至尊木虫 (著名写手)

我叫MT
★ ★ ★ ★ ★
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 ]
赞一下(11人) 回复此楼
一群自以为正义凛然的年轻人将一切不能以科学解释的事情定性为封建迷信并大刀阔斧地进行消灭，其实这是修养不足学识浅薄的一种体现，也是可恶的偏执和愚蠢的自以
6楼2009-11-13 17:41:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuy069

铁杆木虫 (正式写手)
★
yuguotianqing(金币+1,VIP+0):谢谢参与 11-13 20:55
引用回帖:
Originally posted by zhangmt at 2009-11-13 17:41:
抱歉，忘了一般版本在scan的时候不给优化的（内部版据说有哦）
可以用另一种方法实现，试试这个：

windows下跑了4分半，跑了18个点。图形就不贴了。

#p Opt(ModRedundant) hf/sto-3g

Title

0 1
C
...

照着你的试了下，不行啊 连L1都没过去，输入文件无法识别。兄弟，请教一下是怎么回事？
赞一下 回复此楼
7楼2009-11-13 18:16:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangmt

至尊木虫 (著名写手)

我叫MT
★
yjcmwgk(金币+1,VIP+0): 11-16 12:57
让我猜猜看，如果你是在windows下做的，你可能漏掉了文件结尾的空行。空两行更保险。如果你是在linux下做的，需要先执行dos2unix来转换成linux下的文本格式。
不然，难道你用的是g94的版本？
还有，下次说出错的时候，要把最后那几行出错信息贴出来。别让大家猜。。。
这个我已经跑过了的程序，不会有错啊。
赞一下(10人) 回复此楼
一群自以为正义凛然的年轻人将一切不能以科学解释的事情定性为封建迷信并大刀阔斧地进行消灭，其实这是修养不足学识浅薄的一种体现，也是可恶的偏执和愚蠢的自以
8楼2009-11-13 20:25:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuguotianqing

木虫 (正式写手)
下面的这个我试了下，是可以的。
但前面讲的scan 那个要是加上优化结构opt，在高斯03中还有方法运行么
赞一下 回复此楼
开心工作,快乐生活
9楼2009-11-13 20:53:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuy069

铁杆木虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by zhangmt at 2009-11-13 20:25:
让我猜猜看，如果你是在windows下做的，你可能漏掉了文件结尾的空行。空两行更保险。如果你是在linux下做的，需要先执行dos2unix来转换成linux下的文本格式。
不然，难道你用的是g94的版本？
还有，下次说出错的 ...

还是不行啊 下面出错
******************************************
Gaussian 03:  IA32W-G03RevC.02 12-Jun-2004
                13-Nov-2009
******************************************
%mem=512Mb
%nproc=1
Will use up to    1 processors via shared memory.
-----------------------------------
#p opt(ModRedundant) hf/sto-3g scan
-----------------------------------
CompJT:  unrecognized IType=10.
Error termination via Lnk1e in D:\G03W\l1.exe at Fri Nov 13 21:30:58 2009.
Job cpu time:  0 days  0 hours  0 minutes  2.0 seconds.
File lengths (MBytes):  RWF=      7 Int=      0 D2E=      0 Chk=      1 Scr=      1
输入文件就用你最后一次给的
最后一行15什么意思啊  我把15改成B4也不行
Linux下G03 G09和G03w我都试了，都不行 在Linux下当然用了dos2unix，输入文件最后也空了两行
请指教下
赞一下(1人) 回复此楼
10楼2009-11-13 21:39:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yuguotianqing 的主题更新
131/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 211本，11408一志愿中科院277分，曾在中科院自动化所实习 +3 Losir 2026-03-12 4/200 2026-03-16 21:52 by Losir
[考研] 材料专硕326求调剂 +5 墨煜姒莘 2026-03-15 5/250 2026-03-16 21:30 by 木瓜膏
[考研] 一志愿苏州大学材料工程（085601）专硕有科研经历三项国奖两个实用型专利一项省级立项 +3 大火山小火山 2026-03-16 5/250 2026-03-16 16:54 by barlinike
[考研] 0703一志愿211 285分求调剂 +5 ly3471z 2026-03-13 5/250 2026-03-16 16:16 by 哦哦123
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[考研] 326求调剂 +4 上岸的小葡 2026-03-15 5/250 2026-03-16 08:39 by Linda Hu
[考研] 【0703化学调剂】-一志愿华中师范大学-六级475 +5 Becho359 2026-03-11 5/250 2026-03-14 11:35 by 哦哦123
[考研] 学硕285求调剂 +13 Wisjxn 2026-03-12 46/2300 2026-03-14 10:33 by JourneyLucky
[考研] 271求调剂 +10 生如夏花… 2026-03-11 10/500 2026-03-14 00:35 by 卖报员小雨
[考研] 材料工程,326分，求调剂 +6 KRSLSR 2026-03-10 6/300 2026-03-13 23:47 by JourneyLucky
[考研] 材料工程调剂 +4 咪咪空空 2026-03-11 4/200 2026-03-13 19:57 by JourneyLucky
[考研] 求调剂 +5 一定有学上- 2026-03-12 5/250 2026-03-13 18:31 by ms629
[考研] 材料专硕350 求调剂 +4 王金科 2026-03-12 4/200 2026-03-13 16:02 by ruiyingmiao
[考研] 化工学硕306求调剂 +9 42838695 2026-03-12 9/450 2026-03-13 10:16 by houyaoxu
[考研] 工科0856专硕化学工程269能调剂吗 +10 我想读研11 2026-03-10 10/500 2026-03-13 10:14 by Yuyi.
[考研] 270求调剂 085600材料与化工专硕 +3 YXCT 2026-03-11 3/150 2026-03-13 10:13 by houyaoxu
[考研] 298求调剂 +3 Vv呀！ 2026-03-10 3/150 2026-03-10 22:40 by 剑诗杜康
[考研] 求调剂材料专硕293 +6 段_(:з」∠)_ 2026-03-10 6/300 2026-03-10 18:22 by ms629
[考研] 0703化学调剂 +3 三dd. 2026-03-10 3/150 2026-03-10 15:45 by peike
[硕博家园] 木虫好像不热闹了，是不是？ +4 偏振片 2026-03-10 4/200 2026-03-10 09:51 by longwave
信息提示
关闭
请填处理意见
关闭