Computer simulations of particle flux and energy peaking

A computer simulation of the channeling of 300 keV deuterons in a niobium lattice has been performed. The model calculations included: binary Thomas-Fermi interaction between projectiles and lattice atoms, thermal displacements of the lattice atoms in the transverse plane of particle motion, electron scattering of the deuterons dependent on the local electron density, and electronic energy loss dependent on the local electron density. Flux distributions and mean particle energy distributions of the deuterons were calculated in their development with penetration depth. Flux peaking and energy peaking was found to depend sensitively on the local variation of the electron density in the crystal.