Anisotropic target materials are promising candidates for dark matter direct detection experiments, providing a directional sensitivity that can be used to distinguish a dark matter (DM) signal from the various Standard Model backgrounds. In this paper we introduce the Julia and Python implementations of Vector Spaces for Dark Matter (VSDM), which handle the difficult scattering rate computation for these rotating, three-dimensional response functions by calculating a partial rate matrix for every combination of DM velocity distribution, material response function, and particle DM properties.
Lillard et al. (Tue,) studied this question.