Primary Neutron Spectra in Ion Vlasov-Fokker-Planck Simulations

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.