An analytic n-body potential for bcc Iron

We have developed an analytic n-body phenomenological potential for bcc iron made of two electron-density functionals representing repulsion via the Thomas–Fermi free-electron gas kinetic energy term and attraction via a square root functional similar to the second moment approximation of the tight-binding scheme. Electron-density is given by radial, hydrogen-like orbitals with effective charges taken as adjustable parameters fitted on experimental and ab-initio data. Although the set of adjustable parameters is small, prediction of static and dynamical properties of iron is in excellent agreement with the experiments. Advantages and shortcomings of this model are discussed with reference to published works.