24小时热门版块排行榜    

查看: 3433  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

豆豆sk8

银虫 (小有名气)

[求助] 含双键二面角,柔性扫描势能曲线,计算到一半就出错(gaussian) 已有1人参与

各位大神,如题:
用gaussian对丁烯醛C-C=C-C二面角进行柔性扫描。之前使用相同的命令对C=C-C=O成功进行了扫描,但是换成这个中间是双键的二面角就老是出错。在我这里是计算到一般出错,在师兄那里是直接就报错,求指点
%chk=C:\g09w\tests\butenal 2(1.15)\3-24 opt=modredundant trans\notstandard\build structure.chk
# opt=modredundant b3lyp/6-311g(d,p) nosymm scf(maxcycle=300,qc) geom=connectivity

Title Card Required

0 1
C              
H                  1    1.09054186
H                  1    1.09432804    2  107.93359226
H                  1    1.09432776    2  107.93350013    3  115.33687104    0
C                  1    1.50239580    2  111.80729538    4  122.33147321    0
H                  5    1.08929069    1  117.01196721    2   -0.00383048    0
C                  5    1.33840408    1  125.49243777    2  179.99643945    0
C                  7    1.46865845    5  121.10962089    1  179.99971109    0
H                  8    2.18241339    7   26.46832174    5 -179.99972750    0
H                  8    1.11380015    7  114.35094064    5    0.00000000    0
O                  8    1.21028815    7  124.66829891    5  179.99783613    0

1 2 1.0 3 1.0 4 1.0 5 1.0
2
3
4
5 6 1.0 7 2.0
6
7 9 1.0 8 1.0
8 10 1.0 11 2.0
9
10
11

D 1 5 7 8 S 90 -4.00000000
回复此楼
Justskate
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

豆豆sk8

银虫 (小有名气)

引用回帖:
4楼: Originally posted by 枪下游魂 at 2014-04-02 10:46:25
这应该不是方法和基组的问题,换linux也会报错。
l103的错误大多数是来自于分子坐标和结构问题上。
首先你先看下在扫描碳链的时候相应的其他原子是否跟着一起旋转;
然后你要把需要跟着一起动的原子的相对位置固 ...

其他原子是跟着转的,比如C -C =CH1-CHO中 H1和CHO就是同时旋转。当时是专门去编译了坐标让它俩同时旋转的,不知道有没有带来影响。是不是用柔性扫描时不需要让他们同时旋转?这么改是因为之前用SCAN扫描的时候俩分子会在旋转180度后重合。
采用了后两条建议,改成了笛卡尔并同时固定了一个键角,不知道能不能安全着陆....

非常感谢,有问题再向您请教~
Justskate
5楼2014-04-02 11:16:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

枪下游魂

木虫 (著名写手)

出错原因是什么?
2楼2014-04-02 09:17:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

豆豆sk8

银虫 (小有名气)

引用回帖:
2楼: Originally posted by 枪下游魂 at 2014-04-02 09:17:29
出错原因是什么?

刚才正在补充,结果断网了。。。
还想发个图片上来...
首先我看了一下这几次扫描的最后一个图,有一个C和H键是断开的。不知道是因为没有计算完所以断开了还是因为有断开的键才终止了计算。
另外,错误显示2070,最后几句话是:
Old curvilinear step not converged, using linear step
Error imposing constraints
Error termination via Lnk1e in C:\G09W\l103.exe at Thu Mar 27 22:47:20 2014.
Job cpu time:  0 days  2 hours 42 minutes 26.0 seconds.
File lengths (MBytes):  RWF=     10 Int=      0 D2E=      0 Chk=      9 Scr=      1
=======================================
现在能想到的方法就是:
1、师姐建议先换个计算方法和基组再算一下;
2、找个linux系统的gaussian计算一下;
========================================
Justskate
3楼2014-04-02 09:32:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

枪下游魂

木虫 (著名写手)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
gkf高: 金币+3 2014-04-03 11:15:04
引用回帖:
3楼: Originally posted by 豆豆sk8 at 2014-04-02 09:32:10
刚才正在补充,结果断网了。。。
还想发个图片上来...
首先我看了一下这几次扫描的最后一个图,有一个C和H键是断开的。不知道是因为没有计算完所以断开了还是因为有断开的键才终止了计算。
另外,错误显示2070, ...

这应该不是方法和基组的问题,换linux也会报错。
l103的错误大多数是来自于分子坐标和结构问题上。
首先你先看下在扫描碳链的时候相应的其他原子是否跟着一起旋转;
然后你要把需要跟着一起动的原子的相对位置固定;
另外可以尝试下改成笛卡尔坐标。
4楼2014-04-02 10:46:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +7 ASdOkHsho7FD 2026-08-28 10/500 2026-08-30 00:41 by ZPa0EcMwuECS
[基金申请] 投票:  有多少人是今天查系统知道结果的? +16 爱看书的可乐 2026-08-26 18/900 2026-08-30 00:33 by winsaint
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 12/600 2026-08-30 00:31 by ZPa0EcMwuECS
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 G6APbkg8SA6w 2026-08-29 4/200 2026-08-29 22:27 by ZPa0EcMwuECS
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 gy1nBQXYQJqL 2026-08-29 5/250 2026-08-29 21:33 by 4qO8zhfjr2AX
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 gy1nBQXYQJqL 2026-08-29 3/150 2026-08-29 21:32 by 4qO8zhfjr2AX
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-29 21:22 by 4qO8zhfjr2AX
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 为什么到现在没收到通知? +4 tannykie 2026-08-29 4/200 2026-08-29 17:14 by lmz0216
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 14:13 by jCd0dEvKHShX
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-29 09:26 by G6APbkg8SA6w
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
信息提示
请填处理意见