<div dir="ltr"><br><br><div class="gmail_quote"><div dir="ltr">---------- Forwarded message ---------<br>From: Cayo Emilio <<a href="mailto:cayoquimica@gmail.com">cayoquimica@gmail.com</a>><br>Date: ter, 16 de jun de 2015 às 15:02<br>Subject: [molpro-user] Non-convergence error in OPTG program<br>To:  <<a href="mailto:molpro-user@molpro.net">molpro-user@molpro.net</a>><br></div><br><br><div dir="ltr">Dear MOLPRO community.<div>I am trying to compute the energy of a 3-atom system using the program MULTI in different multiplicities. The singlet works just fine, but triplet is a headache</div><div>My input is:</div><div><br></div><div><div> memory,200,m</div><div> thresh,pspace=20</div><div> basis=AVQZ</div><div> SYMMETRY,X,Y</div><div> </div><div> R1=8.20 bohr</div><div> R2=[0.5,1.0,1.5,2.0,2.25,2.5,2.75,3.0,3.25,3.5,3.75,4.0,4.5,4.75,5.0,5.5,6.0]</div><div> theta1=90 deg</div><div> theta2=0 deg</div><div> theta3=-90 deg</div><div> R3=1.0 bohr</div><div> </div><div> geometry={Si1;</div><div>   q1,Si1,R1/2;</div><div>   q2,q1,R3,Si1,theta1;</div><div>   Si2,q1,R1/2,q2,theta3,si1,theta2;</div><div>   Si3,q1,RR,Si1,theta1,q2,theta2}</div><div> </div><div> DO i=1,#R2</div><div> RR=R2(i)</div><div> R(i)=RR</div><div> {multi;maxit,40;wf,42,1,2;state,2}</div><div> {optg,maxiter=60,VARIABLE=energy(1)</div><div>  active,R1}</div><div> RR1(i)=R1</div><div> theta(i)=theta1</div><div> ecas1(i)=energy(1)</div><div> ecas2(i)=energy(2)</div><div> table,R,ecas1,ecas2,RR1</div><div> ENDDO</div></div><div><br></div><div>I am trying to approximate one of the atoms point-by-point to the dimer, and usually after two or three points of the loop the MOLPRO crash with the error: "No convergence. The problem occurs in Multi.". The log says " WARNING: exit with wrong ci vector!!" for many interactions, and i don't understand what this means. Sometimes I try to reedit the input to start from the point (of the R2 vector) it had just crashed and it goes fine, but after several tries and small changes on my parameters values. </div><div>The problem occur only in triplet calculations, with all symmetries of C2v group. The calculations for all four symmetries of the singlet state, for the same parameters of the input above, end correctly. </div><div>If I do not give a higher thresh for p-space the program crashes with the "Non-zero p-space gradient" error, precisely in the same point.<br></div><div><br>I appreciate any help.</div><div>Thank you,</div><div>Cayo E. M. Gonçalves.</div></div></div></div>