
Examples
Here we showcase some selected physics applications that illustrate
results you can obtain with DarkSUSY. Many of those are based on example programs
located in exampels/aux. Have you obtained interesting results with DarkSUSY
that you want us to advertise here? Let us know!
Thermal annihilation cross section



Description
Thermally averaged annihilation rate during freezeout that is
needed to obtain the observed dark matter relic density. Often
used for benchmarking purposes, in particular in the context of
indirect searches for dark matter. The inset shows the impact of
a hard kinematic cutoff for twobody annihilation vs. allowing
for offshell final states.
 Code
examples/aux/oh2_generic_wimp.f
 Journal Ref
JCAP 1807 (2018) 033
[arXiv:1802.03399]

Dark sector relic density



Description
Thermally averaged annihilation rate for dark matter freezeout
in a secluded dark sector, fully taking into account the evolution
of the temperature ratio between the two sectors. This specific plot
assumes a constant thermally averaged annihilation rate (with
g_{S} indicating additional dark sector degrees of freedom).
The same example program can be used for arbitrary dark sector models
featuring 2→2 annihilations.
 Code
examples/aux/oh2_dark_sector.f
 Journal Ref
arXiv:2007.03696

Freezeout beyond kinetic equilibrium



Description
Dark matter annihilation via an schannel resonance is one of the
examples where the usual Boltzmann equation
may be incorrect because kinetic equilibrium is not
maintained during the entire freezeout process. The plot illustrates
the size of this effect for the Scalar Singlet model.
(The couplings are here chosen as indicated in the
bottom panel; for the standard  in this case incorrect 
calculation this would result in a relic density matching the measured one).
 Code
examples/aux/oh2_cBE_ScalarSinglet.f
 Journal Ref
Phys. Rev. D 96 (2017) 11
[arXiv:1706.07433]
arXiv:2103.01944

Selfinteractions and late kinetic decoupling



Description
A simple dark sector model with a scalar mediator, where the coupling is fixed by the relic density for the purpose of this plot (taking into account Sommerfeldenhanced dark matter annihilation into mediator pairs). The blue band indicates the resulting dark matter selfinteraction strength in dwarf galaxies, while the green band shows the cutoff mass in the matter power spectrum due to late kinetic decoupling.
 Code
examples/aux/oh2_vdSIDM.f
 Journal Ref
JCAP 1807 (2018) 033
[arXiv:1802.03399]

Kinetic decoupling in the MSSM



Description
Neutralino dark matter kinetically decouples much earlier than the dark sector example above. The resulting cutoff in the power spectrum (aka the smallest protohalo mass) is strongly modeldependent and spans about 8 orders of magnitude.
 Journal Ref
New J. Phys. 11 (2009) 105027
[arXiv:0903.0189]

Cosmicray accelerated dark matter



Description
Cosmic rays necessarily accelerate a subdominant part of the Galactic dark matter population to relativistic velocities. This allows conventional direct detection experiments (but also neutrino detectors) to probe otherwise inaccessible subGeV dark matter masses. This specific plot assumes a constant spinindependent scattering rate; the same example program can be used for scattering rates with arbitrary dependence on energy and/or momentum transfer.
 Code
examples/aux/DDCR_limits.f
 Journal Ref
Phys. Rev. Lett. 122 (2019) 171801
[arXiv:1810.10543]
JHEP 03 (2020) 118
[arXiv:1909.08632]

Particle yields with U(1), SU(2) and SU(3) corrections



Description
A crucial input for indirect dark matter searches is the particle yield resulting from dark matter decay or annihilation. Radiative corrections can lead to significant modifications of the results from event generators based on treelevel rates. For the MSSM module, all leading correctiosn are fully implemented, stemming from final states with a fermion pair and an additional photon, gluon, electroweak gauge boson or Higgs boson.
 Journal Ref
JHEP 09 (2017) 041
[arXiv:1705.03466 ]

