N-body simulations of dark matter-baryon interactions
Abstract
Dark matter (DM) particles can interact with particles of the Standard Model. Although there exist constraints from direct and indirect detection experiments, the dynamical evolution of astrophysical objects could provide a promising probe for these interactions. Obtaining astrophysical predictions is challenging and limited by our ability to simulate scatterings between DM and baryonic particles within N-body and hydrodynamics simulations. We develop a novel scheme that allows simulating these interacting dark matter (IDM) models and accurately accounts for their angular and velocity dependence, as well as the mass ratio between the DM and baryonic scattering partners. To describe DM-baryon interactions, we use an N-body code together with its implementation of smoothed-particle hydrodynamics and meshless finite mass. The interaction itself is realised in a pairwise fashion by creating a virtual scattering partner from the baryonic particle and letting it interact with a DM particle using a scattering routine initially developed for self-interacting dark matter. After the interaction, the virtual particle is rejoined with the baryonic particle, fulfilling energy and momentum conservation. Through several test problems, we demonstrate that we can reproduce their analytic solutions with our IDM scheme. We comment on various numerical aspects and challenges as well as describe the limitations of our numerical scheme. Furthermore, we study the impact of IDM on halo formation with a collapsing overdensity. Overall, it is possible to accurately model IDM within N-body and hydrodynamics simulations, commonly used in astrophysics. In consequence, our scheme allows for making novel predictions and obtaining new constraints of DM-baryon scattering.