24小时热门版块排行榜    

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

枪下游魂

木虫 (著名写手)

出错原因是什么?
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的回帖

豆豆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的回帖

豆豆sk8

银虫 (小有名气)

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

改了坐标,固定了一个键角...又出错了...也还是之前的问题,优化到最后就是结构错误。。。
Berny optimization.
Using GEDIIS/GDIIS optimizer.
Bend failed for angle     1 -     2 -     3
Tors failed for dihedral     9 -     1 -     2 -     3
Tors failed for dihedral    10 -     1 -     2 -     3
Tors failed for dihedral    11 -     1 -     2 -     3
Tors failed for dihedral     1 -     2 -     3 -     4
Tors failed for dihedral     1 -     2 -     3 -     6
FormBX had a problem.
Error termination via Lnk1e in /home/yhm/g09/l103.exe at Wed Apr  2 04:53:20 2014.
Job cpu time:       0 days  2 hours  1 minutes  4.1 seconds.
File lengths (MBytes):  RWF=     10 Int=      0 D2E=      0 Chk=      3 Scr=      1

换了基组也不行...
怎么都不行...
准备换个简单类似结构算算....唉
Justskate
6楼2014-04-03 09:17:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

枪下游魂

木虫 (著名写手)

★ ★
gkf高: 金币+2 2014-04-03 11:15:10
引用回帖:
6楼: Originally posted by 豆豆sk8 at 2014-04-03 09:17:35
改了坐标,固定了一个键角...又出错了...也还是之前的问题,优化到最后就是结构错误。。。
Berny optimization.
Using GEDIIS/GDIIS optimizer.
Bend failed for angle     1 -     2 -     3
Tors failed  ...

看下是转到几度的时候停掉的
尝试把范围设成-180~180
7楼2014-04-03 11:06:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

豆豆sk8

银虫 (小有名气)

引用回帖:
7楼: Originally posted by 枪下游魂 at 2014-04-03 11:06:15
看下是转到几度的时候停掉的
尝试把范围设成-180~180...

出错的时候大概是在20~40度之间的时候,我是从180开始往下减小计算的,还么有计算到-180就出错了....
Justskate
8楼2014-04-08 09:29:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ficlerlin

木虫 (小有名气)

引用回帖:
8楼: Originally posted by 豆豆sk8 at 2014-04-08 09:29:00
出错的时候大概是在20~40度之间的时候,我是从180开始往下减小计算的,还么有计算到-180就出错了.......

我也遇到了同样的问题,期待专家解答
坚持就是成功
9楼2014-09-11 16:42:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 豆豆sk8 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 315分,诚求调剂,材料与化工085600 +3 13756423260 2026-03-22 3/150 2026-03-22 20:11 by edmund7
[考研] 284求调剂 +5 Zhao anqi 2026-03-22 5/250 2026-03-22 17:38 by barlinike
[考研] 一志愿中南大学化学学硕0703总分337求调剂 +3 niko- 2026-03-22 3/150 2026-03-22 15:15 by 杨杨杨紫
[考研] 311求调剂 +3 26研0 2026-03-20 3/150 2026-03-22 14:46 by ColorlessPI
[考研] 303求调剂 +5 安忆灵 2026-03-22 6/300 2026-03-22 12:46 by 素颜倾城1988
[考研] 一志愿华中科技大学071000,求调剂 +4 沿岸有贝壳6 2026-03-21 4/200 2026-03-22 07:21 by ilovexiaobin
[考研] 268求调剂 +9 简单点0 2026-03-17 9/450 2026-03-21 15:37 by lature00
[考研] 22 350 本科985求调剂,求老登收留 +3 李轶男003 2026-03-20 3/150 2026-03-21 13:28 by 搏击518
[考研] 279分求调剂 一志愿211 +14 chaojifeixia 2026-03-19 15/750 2026-03-21 13:24 by zhukairuo
[考研] 一志愿中国石油大学(华东) 本科齐鲁工业大学 +3 石能伟 2026-03-17 3/150 2026-03-21 02:22 by JourneyLucky
[考研] 材料 336 求调剂 +3 An@. 2026-03-18 4/200 2026-03-21 01:39 by JourneyLucky
[考研] 271材料工程求调剂 +8 .6lL 2026-03-18 8/400 2026-03-21 00:58 by JourneyLucky
[考研] 一志愿西南交大,求调剂 +5 材化逐梦人 2026-03-18 5/250 2026-03-21 00:26 by JourneyLucky
[考研] 304求调剂 +7 司空. 2026-03-18 7/350 2026-03-20 23:08 by JourneyLucky
[考研] 求调剂,一志愿:南京航空航天大学大学 ,080500材料科学与工程学硕,总分289分 +4 @taotao 2026-03-19 4/200 2026-03-20 22:14 by JourneyLucky
[考研] 085600材料与化工 +8 安全上岸! 2026-03-16 8/400 2026-03-20 22:13 by luoyongfeng
[考研] 290求调剂 +7 ^O^乜 2026-03-19 7/350 2026-03-20 21:43 by JourneyLucky
[考研] 材料学硕297已过四六级求调剂推荐 +11 adaie 2026-03-19 11/550 2026-03-20 21:30 by laoshidan
[考研] 招收调剂硕士 +4 lidianxing 2026-03-19 12/600 2026-03-20 12:25 by lidianxing
[考研] 301求调剂 +4 A_JiXing 2026-03-16 4/200 2026-03-17 17:32 by ruiyingmiao
信息提示
请填处理意见