The longitudinal and transverse density- and spin-density response functions in Pure Neutron Matter (PNM) are derived over a wide range of densities within the Time Dependent Local Spin Density Approximation (TDLSDA). The underlying density functional was derived from an Auxiliary Field Diffusion Monte Carlo (AFDMC) computation of the equation of state of unpolarized and fully spin polarized pure neutron matter. In order to assess the dependence of the results on the specific underlying nucleon-nucleon Hamiltonian, we used both the phenomenological Argonne AV8'+UIX force, and local chiral forces up to N2LO. The resulting response function has then been applied to the study of the neutrino mean free path in PNM.