24小时热门版块排行榜    

查看: 2348  |  回复: 12

freefly2016

木虫 (正式写手)

[求助] Gaussian结构优化报错

各位老师好,我用Gaussian 09W优化一个物质的稳定构型,分子式为(C15H16O2)
手绘的化学结构式见附件:

参数为:Job Type: Optimization to a Minimum
Method: Ground State  DFT b3lyp/6-31G (d, p)
Charge/Mult: 0 1

但是运行不了,每次都error,
提示如下:
请问我哪里设置错了吗? 还是我画的结构式有问题呢?谢谢大家。


Gaussian 09:  IA32W-G09RevA.02 11-Jun-2009
                04-May-2016
******************************************
%mem=2000MW
%rwf=BPA-.rwf
%nosave
%chk=C:\Users\X240\Desktop\Gaussiantest\BPA-.chk
----------------------------------------
# opt b3lyp/6-31g(d,p) geom=connectivity
----------------------------------------
1/14=-1,18=20,19=15,26=3,38=1,57=2/1,3;
2/9=110,12=2,17=6,18=5,40=1/2;
3/5=1,6=6,7=101,11=2,16=1,25=1,30=1,71=1,74=-5/1,2,3;
4//1;
5/5=2,38=5/2;
6/7=2,8=2,9=2,10=2,28=1/1;
7//1,2,3,16;
1/14=-1,18=20,19=15/3(2);
2/9=110/2;
99//99;
2/9=110/2;
3/5=1,6=6,7=101,11=2,16=1,25=1,30=1,71=1,74=-5/1,2,3;
4/5=5,16=3/1;
5/5=2,38=5/2;
7//1,2,3,16;
1/14=-1,18=20,19=15/3(-5);
2/9=110/2;
6/7=2,8=2,9=2,10=2,19=2,28=1/1;
99/9=1/99;
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

小范范1989

木虫 (著名写手)

1:不要把chk放到桌面上
2:添加#p 看输出细节错误
3:删除geom=connectivity 和你输入文件后面的 那些信息
It doesn't matter how slow you are, as long as you're determined to get there, you'll get there.
4楼2016-05-04 18:51:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

xieruobing

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
freefly2016: 金币+5 2016-05-04 20:24:46
有个建议,chk尽量不放C盘,另外遇到问题首先把输入输出贴出来,这样才可以看到症结
2楼2016-05-04 16:31:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

freefly2016

木虫 (正式写手)

引用回帖:
2楼: Originally posted by xieruobing at 2016-05-04 16:31:00
有个建议,chk尽量不放C盘,另外遇到问题首先把输入输出贴出来,这样才可以看到症结

谢谢您的回答,输入的结构式我放附件里,谢谢您。

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : BPA.gjf
  • 2016-05-04 16:45:22, 2.47 K
3楼2016-05-04 16:45:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

freefly2016

木虫 (正式写手)

引用回帖:
4楼: Originally posted by 小范范1989 at 2016-05-04 18:51:36
1:不要把chk放到桌面上
2:添加#p 看输出细节错误
3:删除geom=connectivity 和你输入文件后面的 那些信息

谢谢您的回答,按照您的建议1和3改了,还是link die,请问 #p  要在哪里添加呢?  谢谢。

输入和输出文件如下,谢谢您。
Gaussian结构优化报错
g1.jpg

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : BPA.gjf
  • 2016-05-04 20:21:36, 2.51 K
  • 附件 2 : BPA.out
  • 2016-05-04 20:21:41, 4.13 K
5楼2016-05-04 20:23:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

freefly2016

木虫 (正式写手)

引用回帖:
4楼: Originally posted by 小范范1989 at 2016-05-04 18:51:36
1:不要把chk放到桌面上
2:添加#p 看输出细节错误
3:删除geom=connectivity 和你输入文件后面的 那些信息

想给您发金币,发不了哦,您不是应助回贴哦。
6楼2016-05-04 20:25:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xieruobing

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
freefly2016: 金币+15, 谢谢 2016-05-06 09:19:49
freefly2016: 金币+20 2016-05-10 22:50:14
我看了你的输入结构,没有多大问题,另外 你的原子比较少,又是在计算机上,%mem=2000MW 约1.6G应该不合理,你设个60-100MW试试,还有,rwf磁盘分区你并没有真正理解,那就去掉吧,
你可以试试这样写

%nprocshared=12
%mem=60MW
%chk=BPA.chk
# opt b3lyp/6-31g(d,p) geom=connectivity          (PS:在本地计算机上这个无所谓)
7楼2016-05-05 13:40:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

freefly2016

木虫 (正式写手)

引用回帖:
7楼: Originally posted by xieruobing at 2016-05-05 13:40:26
我看了你的输入结构,没有多大问题,另外 你的原子比较少,又是在计算机上,%mem=2000MW 约1.6G应该不合理,你设个60-100MW试试,还有,rwf磁盘分区你并没有真正理解,那就去掉吧,
你可以试试这样写

%nprocsha ...

谢谢您的回答,我试试看
8楼2016-05-06 09:20:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

freefly2016

木虫 (正式写手)

引用回帖:
7楼: Originally posted by xieruobing at 2016-05-05 13:40:26
我看了你的输入结构,没有多大问题,另外 你的原子比较少,又是在计算机上,%mem=2000MW 约1.6G应该不合理,你设个60-100MW试试,还有,rwf磁盘分区你并没有真正理解,那就去掉吧,
你可以试试这样写

%nprocsha ...

谢谢您的回答,这回可以运行了,可是问题是运行了一天都没有完成,卡住不动了。请问又是哪里出了问题呢?谢谢您。


Harris functional with IExCor=  402 diagonalized for initial guess.
ExpMin= 1.61D-01 ExpMax= 5.48D+03 ExpMxC= 8.25D+02 IAcc=1 IRadAn=         1 AccDes= 0.00D+00
HarFok:  IExCor=  402 AccDes= 0.00D+00 IRadAn=         1 IDoV= 1
ScaDFX=  1.000000  1.000000  1.000000  1.000000
FoFCou: FMM=F IPFlag=           0 FMFlag=      100000 FMFlg1=           0
         NFxFlg=           0 DoJE=T BraDBF=F KetDBF=T FulRan=T
         Omega=  0.000000  0.000000  1.000000  0.000000  0.000000 ICntrl=     500 IOpCl=  0
         NMat0=    1 NMatS0=    1 NMatT0=    0 NMatD0=    1 NMtDS0=    0 NMtDT0=    0
         I1Cent=           4 NGrid=           0.
Petite list used in FoFCou.
Requested convergence on RMS density matrix=1.00D-08 within 128 cycles.
Requested convergence on MAX density matrix=1.00D-06.
Requested convergence on             energy=1.00D-06.
No special actions if energy rises.
Integral accuracy reduced to 1.0D-05 until final iterations.
Initial convergence to 1.0D-05 achieved.  Increase integral accuracy.
Gaussian结构优化报错-1
0507.jpg

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : BPA.gjf
  • 2016-05-07 10:12:07, 2.45 K
  • 附件 2 : BPA.out
  • 2016-05-07 10:12:11, 827.75 K
9楼2016-05-07 10:13:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xieruobing

木虫 (著名写手)

引用回帖:
9楼: Originally posted by freefly2016 at 2016-05-07 10:13:06
谢谢您的回答,这回可以运行了,可是问题是运行了一天都没有完成,卡住不动了。请问又是哪里出了问题呢?谢谢您。


Harris functional with IExCor=  402 diagonalized for initial guess.
ExpMin= 1.61D- ...

我看了一下,结构并没有收敛,
         Item               Value     Threshold  Converged?
Maximum Force            0.002724     0.000450     NO
RMS     Force            0.000398     0.000300     NO
Maximum Displacement     0.285639     0.001800     NO
RMS     Displacement     0.064916     0.001200     NO
Predicted change in Energy=-2.622864D-04
还在继续计算,用的内存太小,所以很慢(因为60MW连0.5G都不到,所以。。。。因为不清楚你的电脑性能,就给了一个通常都承受得起的),
只要不报错,就给他继续算,(PS你也要观察优化趋势是否趋于你想要的结构)

如果你电脑可以,就尽量加大看看(可以在优化曲线总提取能量比较低的结构作为初猜进行计算,这样比较省事省时),
10楼2016-05-08 09:54:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 freefly2016 的主题更新
信息提示
请填处理意见