[Wien] final states and case.inc

Michael Gurnett michael.gurnett at kau.se
Sat Oct 2 21:03:08 CEST 2004


I have removed half an electron in case.inc, so that for atom 7 I have

 4 0.00     NUMBER OF ORBITALS (EXCLUDING SPIN), SHIFT
1,-1,2               ( N,KAPPA,OCCUP)
2,-1,2               ( N,KAPPA,OCCUP)
2, 1,1.5             ( N,KAPPA,OCCUP) <== changed from 2 to 1.5
2,-2,4               ( N,KAPPA,OCCUP)

I have also changed case.inm

BROYD  -0.5   YES  (BROYD/PRATT, extra charge (+1 for additional e), norm)
0.10            mixing FACTOR for BROYD/PRATT scheme
0.50  1.50      PW and CLM-scaling factors


However, once the scf cycle has been performed a look at the 2p energies 
show that atom 7 has gone from -6.5 ryd ( unaltered) to -7.22 ryd (when half 
an electron is removed). This change seems rather large.

Is it possible that I have forgotten to do something. I would really 
appreciate if someone could help me with this

Thanks

Michael 





More information about the Wien mailing list