The implementation of the Hartee-Fock potential, used in flux, is based on the solid-state electron densities as tabulated in the book:
J.F. Ziegler, J.P. Biersack and U. Littmark,
"The Stopping and Range of Ions in Solids",
(Pergamon, New York, 1985)
This implementation is already described in the original flux paper.
P.J.M. Smulders and D.O Boerma
Nucl. Instr. and Meth. B29 (1987) 471
The use of a Hartree Fock potential instead of the ZBL potential was found, at least for He in Si, to be superior long ago:
P.J.M. Smulders, A. Dygo and D.O. Boerma
Giant focusing peak and potential dependence observed in a transition
from axial to planar channeling in Si
Nucl. Instr. and Meth. in Phys. Res. B67 (1992) 185-188
Louis Selen at the Technical University of Eindhoven, and I have found a case where this difference is even more striking
L.J.M. Selen, L.J. van IJzendoorn, P.J.M. Smulders, M.J.A. de Voigt,
Planar MeV ion channeling on strained buried nanofilms,
Nucl. Instr. & Meth. in Physics Research, B190 (2002) 570
and
L.J.M. Selen,
Ion-Channeling on Nanostructured Semiconductors,
Thesis, Technische Universiteit Eindhoven, 2001, ISBN 90-386-1779-8
The routines in hf.f describe the potential in terms of electron densities.
Subroutine HFZBL(Z2) reads the relevant data from the file of electron densities and does some initializing.
Functions POT(Z) and DPOT(Z) calculate the continuum string potential and its derivative, as a function of distance Z. They are actually dependent on Z2 as well! The parameter Z2 (IZ2 in the source), is passed via common block COMPOTZ
The routines ff.f and screening.f contain a few other potential- dependent bits and pieces.
The location of the data file is set in file $FLUX7/FLUXLIB/FILE.LOC.
This data file contains subfiles for various values of Z.
To add one or more elements to this file:
Mail a copy to p.j.m.smulders@home.nl
This program asks for a value of Z and then plots the total charge density 4*pi*r**2 * rho(r), by table lookup in file "HFCHARGE.DAT". The resulting graph should be similar to the one given in the ZBL book.
As a free bonus, the program planpot has been thrown in. Program planpot has nothing to do with Monte Carlo simulations, but, to the contrary, with the continuum planar model, a useful tool to describe planar channeling. The program calculates and prints within this model the potential and force, and the planar oscillation wavelength as a function of amplitude, for various potentials (ZBL, Moliere, Lindhard, HF). Only for the first two of these the modifications due to thermal vibrations are incorporated. The program is limited to mono- atomic crystals.