24小时热门版块排行榜    

查看: 3418  |  回复: 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的回帖

枪下游魂

木虫 (著名写手)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +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的回帖
查看全部 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的回帖

豆豆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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +7 scalable 2026-08-24 7/350 2026-08-24 23:01 by anata1209
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 0XLacIJUOj8D 2026-08-24 3/150 2026-08-24 21:59 by BZKMTicpDhFj
[基金申请] 我面上完蛋了 +9 且听虎啸 2026-08-20 10/500 2026-08-24 20:08 by 苦寒来香
[基金申请] filecode,4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤?
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 21/1050 2026-08-24 14:21 by refreshing11
[基金申请] 估计是周四 +4 archvillain 2026-08-18 4/200 2026-08-24 13:53 by zzuzxg
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 什么时候开奖? +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 by 开心的小狮子
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见