24小时热门版块排行榜    

查看: 2718  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 能否退出参与的面上项目解除限项 +20 koalala 2026-08-24 22/1100 2026-08-24 18:01 by Nano-Firefly
[基金申请] 没有任何消息-是不是就凉了 +6 图啦图啦 2026-08-24 6/300 2026-08-24 17:44 by 小小红帆
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +5 scalable 2026-08-24 5/250 2026-08-24 17:03 by spiderzhuzhi
[基金申请] 2026国自然函评费到账 +17 羊腰板 2026-08-21 19/950 2026-08-24 16:52 by iaeyuan
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 今日不放榜?网传国自然预计 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
[基金申请] 我面上完蛋了 +8 且听虎啸 2026-08-20 9/450 2026-08-24 13:52 by zzuzxg
[基金申请] 范进中举一文的中心思想 +6 炎黄贵胄 2026-08-22 7/350 2026-08-24 11:58 by 6543yes
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 明天应该可查了!? +4 chengyan1220 2026-08-23 4/200 2026-08-24 08:47 by gloomy6159
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 看来今天不会放榜了? +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 紫杉醇
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
请填处理意见