[Wien] question on :CHAXXX
whoming
whoming at sohu.com
Thu Nov 6 11:23:37 CET 2003
Dear WIEN2K users,
In my recently calculation, :CHAXXX of one atom
in case.scf file varites between about 7.4 and 5.1 alternately as the following last four scf-cycle and at last
convergy to 5.141914:
:CHA035: TOTAL CHARGE INSIDE SPHERE 35 = 7.468536
:PCS035: PARTIAL CHARGES SPHERE = 35 S,P,D,F,PX,PY,PZ,D-Z2,D-X2Y2,D-XY,D-XZ,D-YZ
:QTL035: 0.0695 3.0379 4.3507 0.0090 1.0166 1.0047 1.0165 0.8659 0.8518 0.8743 0.8884 0.8700
:CHA035: TOTAL CHARGE INSIDE SPHERE 35 = 5.137844
:PCS035: PARTIAL CHARGES SPHERE = 35 S,P,D,F,PX,PY,PZ,D-Z2,D-X2Y2,D-XY,D-XZ,D-YZ
:QTL035: 0.0627 3.0285 2.0364 0.0081 1.0139 1.0007 1.0139 0.2548 0.3065 0.5783 0.3024 0.5944
:CHA035: TOTAL CHARGE INSIDE SPHERE 35 = 7.467941
:PCS035: PARTIAL CHARGES SPHERE = 35 S,P,D,F,PX,PY,PZ,D-Z2,D-X2Y2,D-XY,D-XZ,D-YZ
:QTL035: 0.0695 3.0380 4.3500 0.0090 1.0166 1.0047 1.0165 0.8654 0.8518 0.8743 0.8884 0.8704
:CHA035: TOTAL CHARGE INSIDE SPHERE 35 = 5.141914
:PCS035: PARTIAL CHARGES SPHERE = 35 S,P,D,F,PX,PY,PZ,D-Z2,D-X2Y2,D-XY,D-XZ,D-YZ
:QTL035: 0.0627 3.0285 2.0404 0.0081 1.0139 1.0007 1.0139 0.2567 0.3078 0.5775 0.3019 0.5967
The main difference originates from the d-charges, varing between 4.35 and 2.04.
And :ENE and :DIS from the last four-SCF file are:
:ENE : ********** TOTAL ENERGY IN Ry = -32986.463168 ------------\
:ENE : ********** TOTAL ENERGY IN Ry = -32986.464245 ------------/ difference is about -0.001077Ry
:ENE : ********** TOTAL ENERGY IN Ry = -32986.464216 ------------\ difference is about -0.000029Ry
:ENE : ********** TOTAL ENERGY IN Ry = -32986.464175 ------------/ difference is about -0.000041Ry
:DIS : CHARGE DISTANCE 0.0106337
:DIS : CHARGE DISTANCE 0.0037342
:DIS : CHARGE DISTANCE 0.0026811
:DIS : CHARGE DISTANCE 0.0017367
The 035th atom in case.inst is:
Co 1
Ar 3 5
3, 2,2.0 N
3, 2,2.0 N
3,-3,3.0 N
3,-3,0.0 N
4,-1,1.0 N
4,-1,1.0 N
So Co has 9 electrons in the valance band, and according to the :CHAR035 value,
the Co will have about +1.6(9-7.467941) or +3.9(9-5.141914) charges , corresponding
to Co^+2 or Co^+4 ionic state.
My question is whether the variation of :CHA035 means that it is still not converged?
How to get to the "real" ground state?
Any comments will be appreciated!
Thank you!
Yours,
Whoming
More information about the Wien
mailing list