Firefly and PC GAMESS-related discussion club

Learn how to ask questions correctly

Alex Granovsky
gran@classic.chem.msu.su

what is the symmetry of the system?
The 6 MB attachment should not be the problem and would be indeed very helpful...

Cheers,
Alex


On Tue Sep 28 '10 2:58pm, Pedro Silva wrote
-------------------------------------------
>Dear FF community,

>I am reasonably acquainted with Firefly, but I have just met a nasty problem I cannot solve:
>in one of my RODFT numerical hessian computations, I get an irrealistically high negative frequency (-55400 cm-1).

>The numerical computation of the hessian involves the computation of energy of the structure upon small sequential displacements of each atom along the x, y and z axis, and these values are then used by FF to compute the hessian matrix.
>I found that although most of these energies in my molecule are around -1860.68, several SCF cycles yield final results around -1860.40. this is obviously the reason for the wrong computation of frequencies, but i have not been able to make the "wrong" SCF cycles to converge into the proper energies around -1860.68. I have used very small VIBSIZ (0.001) and I keep getting trouble in only some SCF cycles, but those are enough to throw the computation to the garbage bin.
>My system is MULT=2, the basis set is 6-31+G* (plus the SBKJC pseudopotential in an iron atom), I include initial guess orbitals from a previous computation, $SCF DIIS=.t. $END

>Any suggestions?

>Thank you!
>Pedro S.
>
>
>PS: I can provide the input file upon request. I did not attach it because of its 6 Mb size...