  The energy spectra of unscattered neutrons produced by deuterium-deuterium
and deuterium-tritium fusion reactions are an important diagnostic in High
Energy Density Physics experiments as the spectra are sensitive to the
velocities of reacting ions. Methods exist for calculating these spectra in
radiation-hydrodynamic ("hydro") and Particle-in-Cell ("PiC") simulations. The
spectra are particularly sensitive to the high energy tail of ion velocity
distribution functions since reaction cross sections increase rapidly with the
kinetic energy of a reacting ion pair at the energies achieved in laboratory
experiments. This means both the hydro and PiC method may not be suitable in
certain plasma regimes. The hydro method assumes that the ion velocity
distribution is locally Maxwellian, while the PiC method is subject to
statistical noise that makes it challenging to accurately simulate finer
details of the spectra. In this work, we present a complementary approach: a
method for calculating the neutron spectra in ion Vlasov-Fokker-Planck
simulations in which the velocity distribution function is fully-resolved. The
method is implemented in the spherically-symmetric code iFP which is used to
simulate laser-driven Inertial Confinement Fusion experiments. The method is
computationally intensive as it requires a five-dimensional numerical integral,
but no approximations of the distribution functions or differential cross
sections are required. Results show that deviations of the ion distribution
functions from Maxwellian can have a noticeable effect on neutron spectra in
shock-driven ICF implosions. The method should facilitate more accurate
benchmarking of simulations and experiments.
