To determine the non-spherical wavefunctions the potential has to be lm1 and lm2 dependent. the potential is stored only as lm dependent , therefore a transformation in the following way has to be done : where are the gaunt coeffients. (see notes by B. Drittler)
Note
attention : The gaunt coeffients are stored in an index array only for lm1.gt.lm2 (see subroutine gaunt) R. Zeller Sep. 2000: modified
Warning
In the old version of KKRimp the following is only done for lm1<=lmmaxd && lm2<=lmmaxd && lm3>1
Warning
The following part is commented out in the old version of KKRimp routine!
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=dp), | intent(out), | dimension (lmmaxd, lmmaxd, irmind:irm) | :: | vnspll | ||
real(kind=dp), | intent(in), | dimension (irmind:irm, lmpot) | :: | vins |
Non-spherical part of the potential |
|
real(kind=dp), | intent(in), | dimension (ncleb, 2) | :: | cleb |
GAUNT coefficients (GAUNT) |
|
integer, | intent(in), | dimension (ncleb, 4) | :: | icleb |
Pointer array |
|
integer, | intent(in) | :: | iend | |||
integer, | intent(in) | :: | irm |
Maximum number of radial points |
||
integer, | intent(in) | :: | ncleb |
Number of Clebsch-Gordon coefficients |
||
integer, | intent(in) | :: | lmpot |
(LPOT+1)**2 |
||
integer, | intent(in) | :: | irmind |
IRM-IRNSD |
||
integer, | intent(in) | :: | lmmaxd |
(KREL+KORBIT+1)(LMAX+1)^2 |