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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 含双键二面角，柔性扫描势能曲线，计算到一半就出错（gaussian）
91/1返回列表
查看: 2812  |  回复: 8
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

豆豆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
1楼2014-04-02 08:51:02
已阅   回复此楼   关注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的错误大多数是来自于分子坐标和结构问题上。
首先你先看下在扫描碳链的时候相应的其他原子是否跟着一起旋转;
然后你要把需要跟着一起动的原子的相对位置固定；
另外可以尝试下改成笛卡尔坐标。
赞一下(1人) 回复此楼
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 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭