## Firefly and PC GAMESS-related discussion club

Learn how to ask questions correctly

** Re: SCF problems **
Alex Granovsky

gran@classic.chem.msu.su

Dear Pedro,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...

Tue Sep 28 '10 7:39pm

This message read **753** times