Initialize the potential with the electron-intracell-potentials the intracell-potential is expanded into spherical harmonics . the lm-term of the intracell-potential of the representive atom i is given by the lm contribution of the charge moment of the representive atom i is given by (see notes by b.drittler and u.klemradt) is muffin tin or Wigner-Seitz sphere radius, depending on kshape turned on or off
Note
Attention : rho2ns(...,1)
is the real charge density times
developed into spherical harmonics . (see deck rholm)
C a l c u l a t e
t h e
e l e c t r o n - i n t r a c e l l - p o t e n t i a l s
a n d
t h e
c h a r g e - m o m e n t s
o f
g i v e n
c h a r g e
d e n s i t i e s .
( F o r
e a c h
s p i n - d i r e c t i o n
t h e
p o t e n t i a l
i s
t h e
s a m e
i n
t h e
p o l a r i z e d
c a s e . )
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=dp), | intent(out), | dimension (lmpot, natyp) | :: | cmom |
LM moment of total charge |
|
real(kind=dp), | intent(out), | dimension (lmpot, natyp) | :: | cminst | ||
integer, | intent(in) | :: | lmax |
Maximum l component in wave function expansion |
||
integer, | intent(in) | :: | nspin |
Counter for spin directions |
||
integer, | intent(in) | :: | nstart | |||
integer, | intent(in) | :: | nend | |||
real(kind=dp), | intent(in), | dimension (irmd, lmpot, natyp, 2) | :: | rho2ns |
radial density |
|
real(kind=dp), | intent(out), | dimension (irmd, lmpot, npotd) | :: | v | ||
real(kind=dp), | intent(in), | dimension (irmd, natyp) | :: | r |
Radial mesh ( in units a Bohr) |
|
real(kind=dp), | intent(in), | dimension (irmd, natyp) | :: | drdi |
Derivative dr/di |
|
integer, | intent(in), | dimension (natyp) | :: | irws |
R point at WS radius |
|
integer, | intent(in), | dimension (0:ipand, natyp) | :: | ircut |
R points of panel borders |
|
integer, | intent(in), | dimension (natyp) | :: | ipan |
Number of panels in non-MT-region |
|
integer, | intent(in) | :: | kshape |
Exact treatment of WS cell |
||
integer, | intent(in), | dimension (natyp) | :: | ntcell |
Index for WS cell |
|
integer, | intent(in), | dimension (ngshd, 3) | :: | ilm_map | ||
integer, | intent(in), | dimension (natyp, lmxspd) | :: | ifunm | ||
integer, | intent(in), | dimension (0:lmpot) | :: | imaxsh | ||
real(kind=dp), | intent(in), | dimension (ngshd) | :: | gsh | ||
real(kind=dp), | intent(in), | dimension (irid, nfund, ncelld) | :: | thetas |
shape function THETA=0 outer space THETA=1 inside WS cell in spherical harmonics expansion |
|
integer, | intent(in), | dimension (natyp, lmxspd) | :: | lmsp |
0,1 : non/-vanishing lm=(l,m) component of non-spherical potential |
|
integer, | intent(in) | :: | lmpot |
(LPOT+1)**2 |
||
integer, | intent(in) | :: | natyp |
Number of kinds of atoms in unit cell |