Antideuterons from dark matter annihilations and hadronization model dependence

We calculate the antideuteron yield in dark matter annihilations on an event-by-event basis using the HERWIG++ Monte Carlo generator. We present the resulting antideuteron fluxes for quark and gauge boson final states. As deuteron production in the coalescence model depends on momentum differences between nucleons that are small compared to Lambda(QCD), it is potentially very sensitive to the hadronization model employed. We therefore compare our antideuteron yields to earlier results based on PYTHIA, thereby estimating their uncertainties. We also briefly discuss the importance of n > 2 final states for annihilations of heavy dark matter particles.