Firefly and PC GAMESS-related discussion club

Learn how to ask questions correctly  
We are NATO-free zone

basic question


Given a basis set, need expansion coefficients of exact AOs of isolated open shell atoms and ions (atomic MOs).
The fixed potential is spherical (around the unique nucleus), but depending on the symmetry fixed in $data, calculated occupied orbitals do not satisfy the 3D spanning by proper degeneracy: the only formally respected degeneracy is that of the group Irreps, The R3 group (with S, P, D,...degeneracies) is not available, but the spatial degeneracy can however be recovered by approximating R3 by polyhedral groups like Oh (but the results obtained in C1 symmetry are restored by adding NOSYM=1 in $contrl).
The explanation is likely simple, but I'd like to know why.
Moreover, not only the total energy of the atom/ion depends on the fixed symmetry, but also the lowest value may not correspond to symmetry closest to R3.

Example with Cu2+:

C1 ->TOTAL ENERGY =   -1639.4442846853
$contrl scftyp=uhf DFTTYP=B3PW91 runtyp=energy nprint=7 icharg=2
mult=2 $end
$system mwords=128 timlim=1200 $end
$BASIS gbasis=TZV ndfunc=1 npfunc=1 $END
Cu2+ B3PW91/TZV energy C1
Cu     29      0.000     0.0000          0.0000

D6h -> TOTAL ENERGY =   -1639.4118631252

OH -> TOTAL ENERGY =   -1639.4013245610  (with 3D spanning by degenerated OAs)      

Thanks for any comment.

[ Previous ] [ Next ] [ Index ]           Thu Mar 31 '16 2:41pm
[ Reply ] [ Edit ] [ Delete ]           This message read 535 times