[Wien] Question about wave function and case.almblm file
gutowska at agh.edu.pl
gutowska at agh.edu.pl
Tue Apr 13 18:14:39 CEST 2021
Dear Professor Blaha and users of Wien2k,
I am trying to get wave function from *.radwf and *.almblm files (in
order to calculate Hopfield parameter). However, I have a problem with
local orbitals.
In User Guide there is an information that local orbitals are
k-independent, so I thought, that in order to get local orbital I should
multiply the corresponding basis functions from *radwf by A_lm, B_lm and
C_lm of freely chosen k-point from *almblm file.
However, in *.almblm file the A_lm coefficients of local orbital (8th
and 9th columns) seem to be k-dependent (they are different for
different k-points), e.g.
K-POINT: 0.0208333333 0.0208333333 0.0208333333 83 13
1
1 1 5 jatom,nemin,nemax
1 ATOM
1 1.446759259259260E-004 NUM, weight
0 0 1 0.11018906E-01 0.00000000E+00 -0.34384483E-01
0.00000000E+00 -0.99978407E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00
1 -1 2 0.10910331E-02 -0.10910331E-02 -0.28815122E-02
0.28815122E-02 -0.16313988E-02 0.16313988E-02 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00
1 0 3 0.00000000E+00 -0.15429538E-02 0.00000000E+00
0.40750736E-02 0.00000000E+00 0.23071463E-02 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00
1 1 4 0.10910331E-02 0.10910331E-02 -0.28815122E-02
-0.28815122E-02 -0.16313988E-02 -0.16313988E-02 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00
...
K-POINT: 0.0416666667 0.0416666667 0.0208333333 83 13
2
1 1 5 jatom,nemin,nemax
1 ATOM
1 4.340277777777782E-004 NUM, weight
0 0 1 0.10715529E-01 0.00000000E+00 -0.33458040E-01
0.00000000E+00 -0.99972925E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00
1 -1 2 0.21642513E-02 -0.10795369E-02 -0.57119968E-02
0.28476341E-02 -0.32356411E-02 0.16136986E-02 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00
1 0 3 0.00000000E+00 -0.30607136E-02 0.00000000E+00
0.80779833E-02 0.00000000E+00 0.45758876E-02 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00
1 1 4 0.21642513E-02 0.10795369E-02 -0.57119968E-02
-0.28476341E-02 -0.32356411E-02 -0.16136986E-02 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00
How should I handle that?
I would appreciate any answer or hint.
Best regards,
Sylwia Gutowska
PhD student, AGH UST in Krakow, Poland
More information about the Wien
mailing list