24小时热门版块排行榜    

查看: 1735  |  回复: 9

陈艳cy

铜虫 (小有名气)

[交流] 【求助】高斯03出错 已有3人参与

我的输入文件是
%chk=C:\G03W\C5H4O2.chk
%mem=10MW
%nprocshared=1
#p opt(maxcyc=200) freq=raman ub3lyp/genecp

C5H4O2bAu

0 2
C                  1.59743700   -0.89657200    0.00005100
C                  1.98525300    0.41728100   -0.00000200
C                  0.77983900    1.17535900   -0.00009500
C                 -0.25565800    0.26755700    0.00008200
O                  0.24957800   -1.00509200   -0.00002600
H                  2.14514500   -1.82690400    0.00017400
H                  2.99989200    0.78819300   -0.00015100
H                  0.67883300    2.25193400   -0.00013500
O                 -2.53711600   -0.41423000   -0.00010000
C                 -1.69661500    0.46968500    0.00006000
H                 -1.98510800    1.54149500    0.00055000
Au                -1.70605158   -1.73744284    0.24816673

C H O 0
6-31+G(d)
****
Au 0
Lanl2dz
****

Au 0
Lanl2dz
输出文件的最后几行是
No NMR shielding tensors so no spin-rotation constants.
Leave Link  601 at Sun Oct 24 01:34:32 2010, MaxMem=   10485760 cpu:       1.0
(Enter C:\G03W\l9999.exe)


THAT'S WHAT MAKES US A GREAT COUNTRY.
THE LITTLE THINGS ARE SERIOUS AND THE BIG ONES ARE NOT.
           WILL ROGERS
Error termination request processed by link 9999.
Error termination via Lnk1e in C:\G03W\l9999.exe at Sun Oct 24 01:34:38 2010.
Job cpu time:  0 days 11 hours 24 minutes 36.0 seconds.
File lengths (MBytes):  RWF=     34 Int=      0 D2E=      0 Chk=      6 Scr=      1
麻烦大家帮我看下出了什么问题。谢谢
回复此楼

» 收录本帖的淘帖专辑推荐

高斯使用注意

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ykwang

金虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
这种情况表示优化没成功,多半是因为优化循环的次数达到最大值了。你检查一下。
Nothing_Is_Impossible!
2楼2010-10-24 11:03:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

陈艳cy

铜虫 (小有名气)

引用回帖:
Originally posted by ykwang at 2010-10-24 11:03:39:
这种情况表示优化没成功,多半是因为优化循环的次数达到最大值了。你检查一下。

谢谢。那我应该怎样修改输入文件呢?
3楼2010-10-24 12:27:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ykwang

金虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by 陈艳cy at 2010-10-24 12:27:30:

谢谢。那我应该怎样修改输入文件呢?

你应该列出较详细的情况(比如四个优化指标中有几个“YES”等)才好判断并提出具体的建议。如果仅是优化步骤达到了最大值,那么只需取其中能量最小的构型作为输入几何重新优化即可,因为你的输入文件中所设的参数MaxCyc=200已经不小了。
Nothing_Is_Impossible!
4楼2010-10-24 16:36:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

陈艳cy

铜虫 (小有名气)

引用回帖:
Originally posted by ykwang at 2010-10-24 16:36:00:

你应该列出较详细的情况(比如四个优化指标中有几个“YES”等)才好判断并提出具体的建议。如果仅是优化步骤达到了最大值,那么只需取其中能量最小的构型作为输入几何重新优化即可,因为你的输入文件中所设的参 ...

这是后面一部分输出文件

       Item               Value     Threshold  Converged?
Maximum Force            0.000016     0.000450     YES
RMS     Force            0.000004     0.000300     YES
Maximum Displacement     0.027442     0.001800     NO
RMS     Displacement     0.006013     0.001200     NO
Predicted change in Energy=-6.954421D-07
Optimization stopped.
    -- Number of steps exceeded,  NStep= 100
    -- Flag reset to prevent archiving.
                           ----------------------------
                           ! Non-Optimized Parameters !
                           ! (Angstroms and Degrees)  !
--------------------------                            --------------------------
! Name  Definition              Value          Derivative Info.                !
--------------------------------------------------------------------------------
! R1    R(1,2)                  1.3714         -DE/DX =    0.0                 !
! R2    R(1,5)                  1.3504         -DE/DX =    0.0                 !
! R3    R(1,6)                  1.0803         -DE/DX =    0.0                 !
! R4    R(2,3)                  1.4226         -DE/DX =    0.0                 !
! R5    R(2,7)                  1.0811         -DE/DX =    0.0                 !
! R6    R(3,4)                  1.3792         -DE/DX =    0.0                 !
! R7    R(3,8)                  1.0822         -DE/DX =    0.0                 !
! R8    R(4,5)                  1.3691         -DE/DX =    0.0                 !
! R9    R(4,10)                 1.4486         -DE/DX =    0.0                 !
! R10   R(5,12)                 3.5038         -DE/DX =    0.0                 !
! R11   R(9,10)                 1.2256         -DE/DX =    0.0                 !
! R12   R(9,12)                 2.7699         -DE/DX =    0.0                 !
! R13   R(10,11)                1.1078         -DE/DX =    0.0                 !
! R14   R(10,12)                3.6708         -DE/DX =    0.0                 !
! A1    A(2,1,5)              111.0645         -DE/DX =    0.0                 !
! A2    A(2,1,6)              133.2139         -DE/DX =    0.0                 !
! A3    A(5,1,6)              115.7217         -DE/DX =    0.0                 !
! A4    A(1,2,3)              105.7161         -DE/DX =    0.0                 !
! A5    A(1,2,7)              126.5696         -DE/DX =    0.0                 !
! A6    A(3,2,7)              127.7142         -DE/DX =    0.0                 !
! A7    A(2,3,4)              106.5839         -DE/DX =    0.0                 !
! A8    A(2,3,8)              127.4656         -DE/DX =    0.0                 !
! A9    A(4,3,8)              125.9504         -DE/DX =    0.0                 !
! A10   A(3,4,5)              109.5241         -DE/DX =    0.0                 !
! A11   A(3,4,10)             130.6213         -DE/DX =    0.0                 !
! A12   A(5,4,10)             119.8546         -DE/DX =    0.0                 !
! A13   A(1,5,4)              107.1113         -DE/DX =    0.0                 !
! A14   A(1,5,12)             144.2005         -DE/DX =    0.0                 !
! A15   A(4,5,12)             103.1884         -DE/DX =    0.0                 !
! A16   A(4,10,9)             126.2159         -DE/DX =    0.0                 !
! A17   A(4,10,11)            113.4927         -DE/DX =    0.0                 !
! A18   A(4,10,12)             94.1337         -DE/DX =    0.0                 !
! A19   A(9,10,11)            120.2914         -DE/DX =    0.0                 !
! A20   A(11,10,12)           147.7723         -DE/DX =    0.0                 !
! A21   A(5,12,9)              52.7            -DE/DX =    0.0                 !
! A22   A(5,12,10)             39.657          -DE/DX =    0.0                 !
! D1    D(5,1,2,3)             -0.0225         -DE/DX =    0.0                 !
! D2    D(5,1,2,7)            179.949          -DE/DX =    0.0                 !
! D3    D(6,1,2,3)           -179.9932         -DE/DX =    0.0                 !
! D4    D(6,1,2,7)             -0.0217         -DE/DX =    0.0                 !
! D5    D(2,1,5,4)             -0.0016         -DE/DX =    0.0                 !
! D6    D(2,1,5,12)          -146.6138         -DE/DX =    0.0                 !
! D7    D(6,1,5,4)            179.9747         -DE/DX =    0.0                 !
! D8    D(6,1,5,12)            33.3625         -DE/DX =    0.0                 !
! D9    D(1,2,3,4)              0.0375         -DE/DX =    0.0                 !
! D10   D(1,2,3,8)            179.916          -DE/DX =    0.0                 !
! D11   D(7,2,3,4)           -179.9336         -DE/DX =    0.0                 !
! D12   D(7,2,3,8)             -0.0551         -DE/DX =    0.0                 !
! D13   D(2,3,4,5)             -0.0398         -DE/DX =    0.0                 !
! D14   D(2,3,4,10)           179.9628         -DE/DX =    0.0                 !
! D15   D(8,3,4,5)           -179.9207         -DE/DX =    0.0                 !
! D16   D(8,3,4,10)             0.0819         -DE/DX =    0.0                 !
! D17   D(3,4,5,1)              0.0263         -DE/DX =    0.0                 !
! D18   D(3,4,5,12)           160.7199         -DE/DX =    0.0                 !
! D19   D(10,4,5,1)          -179.976          -DE/DX =    0.0                 !
! D20   D(10,4,5,12)          -19.2824         -DE/DX =    0.0                 !
! D21   D(3,4,10,9)          -179.1585         -DE/DX =    0.0                 !
! D22   D(3,4,10,11)            0.8985         -DE/DX =    0.0                 !
! D23   D(3,4,10,12)         -162.0831         -DE/DX =    0.0                 !
! D24   D(5,4,10,9)             0.8443         -DE/DX =    0.0                 !
! D25   D(5,4,10,11)         -179.0986         -DE/DX =    0.0                 !
! D26   D(5,4,10,12)           17.9198         -DE/DX =    0.0                 !
! D27   D(1,5,12,9)           167.7867         -DE/DX =    0.0                 !
! D28   D(1,5,12,10)          157.5039         -DE/DX =    0.0                 !
! D29   D(4,5,12,9)            20.4833         -DE/DX =    0.0                 !
! D30   D(4,5,12,10)           10.2006         -DE/DX =    0.0                 !
! D31   D(4,10,12,5)           -9.4037         -DE/DX =    0.0                 !
! D32   D(11,10,12,5)        -159.1831         -DE/DX =    0.0                 !
--------------------------------------------------------------------------------
GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad

Largest change from initial coordinates is atom   11       0.618 Angstoms.
Leave Link  103 at Sun Oct 24 01:34:30 2010, MaxMem=   10485760 cpu:       0.0
(Enter C:\G03W\l202.exe)
                          Input orientation:                          
---------------------------------------------------------------------
Center     Atomic     Atomic              Coordinates (Angstroms)
Number     Number      Type              X           Y           Z
---------------------------------------------------------------------
    1          6             0        1.649640   -0.807202   -0.086449
    2          6             0        2.052469    0.499521    0.018302
    3          6             0        0.857411    1.271140    0.033728
    4          6             0       -0.189163    0.377935   -0.062100
    5          8             0        0.303149   -0.897454   -0.135872
    6          1             0        2.184443   -1.744593   -0.134562
    7          1             0        3.071172    0.856794    0.076867
    8          1             0        0.769572    2.347195    0.107662
    9          8             0       -2.486483   -0.262203   -0.168087
   10          6             0       -1.620137    0.601381   -0.092706
   11          1             0       -1.901149    1.671774   -0.041573
   12         79             0       -2.435497   -2.883023    0.733366
---------------------------------------------------------------------
                    Distance matrix (angstroms):
                    1          2          3          4          5
     1  C    0.000000
     2  C    1.371411   0.000000
     3  C    2.227459   1.422602   0.000000
     4  C    2.187770   2.246366   1.379245   0.000000
     5  O    1.350418   2.243976   2.244720   1.369097   0.000000
     6  H    1.080292   2.253182   3.299087   3.185024   2.063229
     7  H    2.194609   1.081125   2.252617   3.298242   3.283991
     8  H    3.280613   2.251158   1.082163   2.196811   3.287036
     9  O    4.172674   4.606197   3.684223   2.387193   2.861229
    10  C    3.560279   3.675695   2.569593   1.448638   2.438728
    11  H    4.330756   4.124179   2.788518   2.146005   3.386555
    12  Au   4.655048   5.665223   5.346943   4.038889   3.492598
                    6          7          8          9         10
     6  H    0.000000
     7  H    2.756484   0.000000
     8  H    4.336273   2.742191   0.000000
     9  O    4.900628   5.674477   4.181733   0.000000
    10  C    4.469919   4.701316   2.966264   1.225567   0.000000
    11  H    5.326563   5.040059   2.758843   2.024572   1.107847
    12  Au   4.836648   6.688843   6.165967   2.771987   3.672639
                   11         12
    11  H    0.000000
    12  Au   4.651047   0.000000
Stoichiometry    C5H4AuO2(2)
Framework group  C1[X(C5H4AuO2)]
Deg. of freedom    30
Full point group                 C1      NOp   1
Largest Abelian subgroup         C1      NOp   1
Largest concise Abelian subgroup C1      NOp   1
                         Standard orientation:                        
---------------------------------------------------------------------
Center     Atomic     Atomic              Coordinates (Angstroms)
Number     Number      Type              X           Y           Z
---------------------------------------------------------------------
    1          6             0       -2.907997   -1.519623    0.210558
    2          6             0       -4.077499   -0.912907   -0.170136
    3          6             0       -3.768963    0.469449   -0.303160
    4          6             0       -2.431683    0.605819    0.005710
    5          8             0       -1.907674   -0.619160    0.320780
    6          1             0       -2.643373   -2.543637    0.430559
    7          1             0       -5.030682   -1.396419   -0.332889
    8          1             0       -4.439267    1.268764   -0.591022
    9          8             0       -0.392299    1.807817    0.313616
   10          6             0       -1.585858    1.781479    0.036596
   11          1             0       -2.120310    2.716917   -0.221556
   12         79             0        1.535005   -0.153167   -0.038451
---------------------------------------------------------------------
Rotational constants (GHZ):      2.9838513      0.4425714      0.3890560
Leave Link  202 at Sun Oct 24 01:34:30 2010, MaxMem=   10485760 cpu:       0.0
(Enter C:\G03W\l601.exe)
Copying SCF densities to generalized density rwf, ISCF=1 IROHF=0.








No NMR shielding tensors so no spin-rotation constants.
Leave Link  601 at Sun Oct 24 01:34:32 2010, MaxMem=   10485760 cpu:       1.0
(Enter C:\G03W\l9999.exe)


THAT'S WHAT MAKES US A GREAT COUNTRY.
THE LITTLE THINGS ARE SERIOUS AND THE BIG ONES ARE NOT.
           WILL ROGERS
Error termination request processed by link 9999.
Error termination via Lnk1e in C:\G03W\l9999.exe at Sun Oct 24 01:34:38 2010.
Job cpu time:  0 days 11 hours 24 minutes 36.0 seconds.
File lengths (MBytes):  RWF=     34 Int=      0 D2E=      0 Chk=      6 Scr=      1
5楼2010-10-24 17:26:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ykwang

金虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
注意其中的
    -- Number of steps exceeded,  NStep= 100
表示优化达到了最大次数100了(100是缺省的优化次数)。要改变这一缺省次数为200,可在输入文件添加关键字MaxStep=200,并在已有的结构基础上继续优化即可。
Nothing_Is_Impossible!
6楼2010-10-24 17:35:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

陈艳cy

铜虫 (小有名气)

引用回帖:
Originally posted by ykwang at 2010-10-24 17:35:12:
注意其中的
    -- Number of steps exceeded,  NStep= 100
表示优化达到了最大次数100了(100是缺省的优化次数)。要改变这一缺省次数为200,可在输入文件添加关键字MaxStep=200,并在已有的结构基础上继续优化 ...

Nstep不是步长吗?优化最大次数不是用maxcyc来定义的吗
7楼2010-10-24 18:13:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ykwang

金虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by 陈艳cy at 2010-10-24 18:13:59:

Nstep不是步长吗?优化最大次数不是用maxcyc来定义的吗

MaxCycle是用来定义SCF即自洽迭代的最大次数的,MaxStep才是用来定义几何优化的最大次数的。详见Gaussian03/09手册。
Nothing_Is_Impossible!
8楼2010-10-24 18:37:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

陈艳cy

铜虫 (小有名气)

引用回帖:
Originally posted by ykwang at 2010-10-24 18:37:29:

MaxCycle是用来定义SCF即自洽迭代的最大次数的,MaxStep才是用来定义几何优化的最大次数的。详见Gaussian03/09手册。

呵呵,谢谢!
9楼2010-10-25 21:29:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xxuusq

新虫 (初入文坛)

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
小红豆(金币+3): 感谢交流 2011-07-11 19:00:44
Originally posted by ykwang at 2010-10-24 18:37:29:

MaxCycle是用来定义SCF即自洽迭代的最大次数的,MaxStep才是用来定义几何优化的最大次数的。详见Gaussian03/09手册。

-----------------------------------------------
这个不对,在opt关键字中,maxstep是定义的最大步长,maxcycle才是定义几何优化的最大次数的。scf关键字中MaxCycle是用来定义自洽迭代的最大次数。
10楼2011-07-11 00:37:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 陈艳cy 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见