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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】用Scan来做一个多原子的分子的势能面。输入文件怎么建。
131/2返回列表12下一页
查看: 1681  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 WqtG0Aem8Ufl 2026-10-24 5/250 2026-10-26 02:55 by cLOojpOxkX0G
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 3xN9SualEUPt 2026-10-24 5/250 2026-10-26 02:38 by cLOojpOxkX0G
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 3xN9SualEUPt 2026-10-24 5/250 2026-10-26 02:25 by cLOojpOxkX0G
[论文投稿] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 3xN9SualEUPt 2026-10-24 5/250 2026-10-26 02:18 by cLOojpOxkX0G
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 3xN9SualEUPt 2026-10-24 4/200 2026-10-26 02:16 by cLOojpOxkX0G
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 3xN9SualEUPt 2026-10-24 5/250 2026-10-26 02:12 by cLOojpOxkX0G
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +4 3xN9SualEUPt 2026-10-24 5/250 2026-10-26 02:05 by cLOojpOxkX0G
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 3xN9SualEUPt 2026-10-24 4/200 2026-10-26 02:03 by cLOojpOxkX0G
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 DI1qFHFK0KbE 2026-10-24 4/200 2026-10-26 01:50 by cLOojpOxkX0G
[硕博家园] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 DI1qFHFK0KbE 2026-10-24 3/150 2026-10-26 01:37 by cLOojpOxkX0G
[考博] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 DI1qFHFK0KbE 2026-10-24 5/250 2026-10-26 01:33 by cLOojpOxkX0G
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +5 7xHjyB3bAiwg 2026-10-23 12/600 2026-10-26 00:54 by cLOojpOxkX0G
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 7xHjyB3bAiwg 2026-10-23 12/600 2026-10-26 00:45 by cLOojpOxkX0G
[考研] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 WqtG0Aem8Ufl 2026-10-24 3/150 2026-10-25 20:11 by cLOojpOxkX0G
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 3xN9SualEUPt 2026-10-24 3/150 2026-10-25 14:23 by XiYP04It3Evb
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +5 7xHjyB3bAiwg 2026-10-23 10/500 2026-10-25 12:51 by XiYP04It3Evb
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 7xHjyB3bAiwg 2026-10-24 4/200 2026-10-25 01:22 by XcOYG53d3x4p
[文学芳草园] 亲密关系 +7 宋昙 2026-10-21 13/650 2026-10-23 17:07 by 小牧11
[考博] 27秋申博自荐 +4 且听风吟时 2026-10-21 4/200 2026-10-23 13:26 by 大大大肥虫
[考博] 郑州大学申博 +5 stars~ 2026-10-19 8/400 2026-10-22 15:06 by stars~
信息提示
关闭
请填处理意见
关闭