mod_vintras Module

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)


Subroutines

public subroutine vintras(cmom, cminst, lmax, nspin, nstart, nend, rho2ns, v, r, drdi, irws, ircut, ipan, kshape, ntcell, ilm_map, ifunm, imaxsh, gsh, thetas, lmsp, lmpot, natyp)

Author
B. Drittler
License
Creative Commons License
Category
potential, KKRhost

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 . )

Read more…

Arguments

Type IntentOptional 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