The Hartree-Fock potential

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


  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:

1 Add a tag, such as **oxygen, with '**' in columns 1,2.
2 one line with Ndata, Z2
3 Ndata lines with radius(A), density(electrons/A**3)
4 To check the data, run program "hfshow" (or run 'make test'), and compare the graph to the one in the ZBL book.
5 Also check if the total charge printed has the correct value.
6 and last but not least:

Mail a copy to


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.