[molpro-user] puzzle about the reaction path following options

MHT mht at mrdlab.icas.ac.cn
Wed Sep 6 09:00:33 BST 2006


Hi,

I did a reaction path calculation/IRC with idir=1, which means the first 
optimization step will be along the transition vector from TS.
Normally, it will go to "an asymptote".
But sometimes this asymptote is not so "smooth", there are one or two or 
more points which are a little higher than this asymptote. 
Although, this job could be normally terminated, it is quite strange. How 
can the energy of these strange points go higher than asymptote.
See the following output:

 END OF GEOMETRY OPTIMIZATION.    TOTAL CPU:      3942.9 SEC
                                                                                                                               
 ITER.   ENERGY(OLD)    ENERGY(NEW)      DE          GRADMAX     GRADNORM    
GRADRMS     STEPMAX     STEPLEN     STEPRMS
   1  -152.90707395  -152.90707561    -0.00000165  0.00000117  0.00000143  
0.00000010  0.01158242  1.00000000  0.06666667
   2  -152.90707561  -152.90713573    -0.00006013  0.00010077  0.00018109  
0.00001207  0.07969603  0.13160862  0.00877391
   3  -152.90713573  -152.90721859    -0.00008286  0.00226480  0.00269265  
0.00017951  0.04601674  0.07766024  0.00517735
   4  -152.90721859  -152.90729683    -0.00007823  0.00073856  0.00127870  
0.00008525  0.04450214  0.07562299  0.00504153
   5  -152.90729683  -152.90738954    -0.00009272  0.00158575  0.00221004  
0.00014734  0.04626073  0.07690365  0.00512691
   6  -152.90738954  -152.90748236    -0.00009282  0.00162790  0.00214846  
0.00014323  0.04388172  0.07636533  0.00509102
   7  -152.90748236  -152.90758151    -0.00009915  0.00190875  0.00260383  
0.00017359  0.04739069  0.07687329  0.00512489
   8  -152.90758151  -152.90767629    -0.00009478  0.00183900  0.00242181  
0.00016145  0.04508402  0.07641860  0.00509457
   9  -152.90767629  -152.90777065    -0.00009435  0.00183494  0.00252928  
0.00016862  0.04810958  0.07647086  0.00509806
  10  -152.90777065  -152.90785786    -0.00008722  0.00145822  0.00204495  
0.00013633  0.04665480  0.07566115  0.00504408
  11  -152.90785786  -152.90793198    -0.00007412  0.00066289  0.00106386  
0.00007092  0.04768685  0.07596757  0.00506450
  12  -152.90793198  -152.90799646    -0.00006448  0.00055258  0.00089460  
0.00005964  0.05031503  0.08149318  0.00543288
  13  -152.90799646  -152.90803970    -0.00004324  0.00042039  0.00068701  
0.00004580  0.04668051  0.07634139  0.00508943
  14  -152.90803970  -152.90806354    -0.00002384  0.00026877  0.00045157  
0.00003010  0.04571338  0.07621047  0.00508070
  15  -152.90806354  -152.90806761    -0.00000406  0.00010776  0.00019717  
0.00001314  0.02897178  0.04941772  0.00329451
  16  -152.90806761  -152.90806766    -0.00000005  0.00004433  0.00006608  
0.00000441  0.00248125  0.00380770  0.00025385
  17  -152.90806766  -152.90806766     0.00000000  0.00000817  0.00001507  
0.00000100  0.00013310  0.00020669  0.00001378
  18  -152.90806766  -152.90806766     0.00000000  0.00000240  0.00000329  
0.00000022  0.00009225  0.00014318  0.00000955
  19  -152.90806766  -152.90806766     0.00000000  0.00000026  0.00000055  
0.00000004  0.00000895  0.00000925  0.00000062
  20  -152.90806766  -152.90790849     0.00015917  0.00000040  0.00000054  
0.00000004  0.13188970  0.19228404  0.01281894
  21  -152.90790849  -152.90806766    -0.00015917  0.00522778  0.00692745  
0.00046183  0.13189403  0.19229438  0.01281963
  22  -152.90806766  -152.90806766     0.00000000  0.00000014  0.00000020  
0.00000001  0.00000188  0.00000288  0.00000019
  23  -152.90806766  -152.90758378     0.00048388  0.00000036  0.00000045  
0.00000003  0.11612722  0.16830718  0.01122048
  24  -152.90758378  -152.90806766    -0.00048388  0.01702191  0.01971109  
0.00131407  0.11612209  0.16830180  0.01122012
  25  -152.90806766  -152.90806766     0.00000000  0.00000006  0.00000011  
0.00000001  0.00000449  0.00000483  0.00000032


yours sincerely
MHT




More information about the Molpro-user mailing list