Full text of publication follows: The slow down density of recoils produced by 14 MeV neutron in LiD is investigated. Of particular interest are the deuterium (D) recoils in the energy range of 0.1 to 12.444 MeV to bombard the tritium in LiD to induce the fusion reaction. The Deuteron recoil density is the key issue to evaluate whether the fusion reaction rate reaches saturation at the tritium (T) concentration {approx}0.5x10{sup 22} Tcm{sup -3} in LiD with certain tritium production rate. The slow down density by internal irradiation and the thermal displacement mode in LiD is calculated with the Marlowe code. LiD is considered to have a body centre cubic lattice. This procedure allows evaluating the total kinetic energy of moving atoms and their energy spectra at various times during the development of the cascade. Fast Exponential-Sum Screened-Coulomb Potential is taken for the interatomic potential and inelastic energy loss is taken into account as non-local and local inelastic energy losses. Results are fitted with the following formulae: q1(E, E{sub r})={xi}2E/E{sub r} (E > E{sub r}), q1(E, E{sub r})= 0 (E < E{sub r}), where q1(E, E{sub r}) is the number of recoils with energy E{sub r} produced by a PKA with energy E. {xi} is the fitting coefficient reflecting the lattice, interatomic potentials effects and inelastic energy loss in the LiD lattice. Results are presented here.(authors)