## Firefly and PC GAMESS-related discussion club

Learn how to ask questions correctly

** SCMCF convergence problems **
Andrey

andrei_ipatov@mail.ru

Sorry, dear collegues, in my calculations of the O2 excitation spectra I see some strange (to my up to now not very experienced point of view) dependence of the results on the basis size.! ================================================================

! Input file for FireFly

! ================================================================

$CONTRL SCFTYP=MCSCF RUNTYP=ENERGY EXETYP=RUN $END

$SYSTEM MWORDS=12 $END

$CONTRL COORD=UNIQUE UNITS=ANGS $END

$MCSCF CISTEP=ALDET SOSCF=1 $END

$DET GROUP=C1 NSTATE=20 NCORE=2 NELS=12 NACT=8 $END

$DATA

Oxygen molecule

Dnh 4

O 8.000000 0.000000 0.000000 -0.6035

! here was the basis set data !

$END

$END

$GUESS GUESS=MOREAD NORB=40 $END

$VEC

! here were the guess vectors !

$END

I made a set of calculations with different basises,

the basis from the internal GAMESS databese 6-311G** (the smallest one among the tested ones), then with the basis sets from the EMSL web page :

PVTZ-Gamess, aug-cc-pVDZ, aug-cc-pVDZ-DK, aug-cc-pVTZ and finally the totally uncontracted aug-cc-pVTZ. I expected that the ground state (triplet X³\Sigma⁻_g ) energy and the energies of the lowest excited states ( singlets a¹\Delta_g and b¹\Sigma⁺_g ) will lower with increase of the basis size, but instead I got the minimal values for the PVTZ-Gamess case and, on the contrary, increase of the total energies for aug-cc-pVDZ and aug-cc-pVDZ-DK.

Moreover, for the larger basis sets, aug-cc-pVTZ, in the converged situation I see, that the double degeneracy of the \Delta state (and also for the other higher states with the angular momentum projection > 0, that must be doubly degenerated in the case of O2 molecule) has disappeared, i.e. these levels arise to be splitted ! In contrast, for the smaller basis sets the degeneracy was practically perfect.

What can be the reason of such strange unphysical dependence ? What I am doing in a wrong way ?

Thanks in advance for any help.

Sincerely

Andrey

Sat Mar 10 '12 5:34pm

This message read **1315** times