Dear all,(adsbygoogle = window.adsbygoogle || []).push({});

Are you familiar with the Elk-code software? It is a software about FLAPW (full potential linear augumented plane wave).Recently,I am reading the source code,but I meet many difficulties.

For example,I've a question related to charge density.In the Elk source code, the charge density is divided two parts,i.e.,rhomt(:,:,:) in the spherical coordinates and rhoir(:) in the Cartesian coordinates.But I don't find the angle (theta,phi) for the rhomt(:,:,:) and the position(x,y,z) for the rhoir(:) in the "rhomag.f90".

In addition,I want to learn about the calculating of the exchange-correlation energy density and potential,but I still can't find the (theta,phi) and (x,y,z),taking the xc_xalpha.f90 as an exmaple as follows:

! loop over density points

do i=1,n

r=rho(i)

if (r.gt.1.d-12) then

rs=(3.d0/(4.d0pir))(1.d0/3.d0)

rsm1=1.0d0/rs

! compute energy density

exc(i)=-alphaefacrsm1

! compute potential

vxc(i)=-alphavfacrsm1

end if

end do

Here,the "i=1,n" stands for the value of the radial meshgrid in the muffin-tin sphere and the value of G-point in the interstitial region.As my understanding,the energy density or potential should be for every point in the real space,so I don't clear why the (theta,phi) and (x,y,z) can't be found in those subroutines.

If anyone give me a hint or a suggestion,I will be thankful.

Regards,

**Physics Forums - The Fusion of Science and Community**

Join Physics Forums Today!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# [Fortran] Is anyone familiar with the Elk-code software?

Loading...

Similar Threads for Fortran anyone familiar |
---|

Fortran Successive Over Relaxation Method in FORTRAN |

If statement for Monte Carlo integration |

C/++/# How to use C++ in studying calculus |

Fortran How to buy Fortran Programming Software Commersial License |

Fortran Fortran for Computational Fluid Dynamics Problems |

**Physics Forums - The Fusion of Science and Community**