Firefly and PC GAMESS-related discussion club

Learn how to ask questions correctly

Gordon
gordon.driver@abo.fi

I hadn't thought of it in this way previously, thanks for the "input"!

Cheers,
Gordon

On Mon May 17 '10 3:03pm, Vyacheslav wrote
------------------------------------------
>Hi!
>Sometimes it is necessary to overthrow the starting geometry in other valley - what leads to a local minimum instead of to a saddle point. An example of calculation for your structure is applied. However in this run symmetry of structure has strongly decreased.
>Best regards,
>          Vyacheslav
>--------------------------------------------------------------
>
>
>
>On Mon May 17 '10 8:10am, Gordon wrote
>--------------------------------------
>>Dear forum followers,

>>I have been optimising geometries to calculate gas phase proton affinities with FF, using MP2 6-311+G(d,p) which will generally yield calc'd values within +/- 10 kJ/mol of experimental values. Now I am trying the same business using $BASIS NDFUNC=2 NFFUNC=1 NPFUNC=2 DIFFSP=.TRUE. DIFFS=.TRUE. $END. With this I can optimise the geometry, BUT, I also observe a single imaginary frequency. To try to fix this, I then used the same input and included $CONTRL RUNTYP=Optimize ICUT=30 ITOL=30 INTTYP=HONDO $END and $SCF DIRSCF=.TRUE. NCONV=10 FDIFF=.FALSE. $END.  (I have attached the two outputs for the acetate anion here).  I have also tried with ICUT=25, but I obtain the same result.

>>I am not sure what I could do now to get further down the PE well!  Perhaps directed adjustment of these ICUT, ITOL, NCONV values are in order?  After 3 tries I was hoping to receive a suggestion from those who are more knowledgeable.  Perhaps there is another trick to try?

>>I would be most grateful for any suggestions.

>>Cheers,
>>Gordon