darksusy is hosted by Hepforge, IPPP Durham

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 freeze-out 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 two-body annihilation vs. allowing for off-shell 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 freeze-out 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 gS 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

Freeze-out beyond kinetic equilibrium

  • Description
    Dark matter annihilation via an s-channel resonance is one of the examples where the usual Boltzmann equation may be incorrect because kinetic equilibrium is not maintained during the entire freeze-out 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

Self-interactions 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 Sommerfeld-enhanced dark matter annihilation into mediator pairs). The blue band indicates the resulting dark matter self-interaction 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 model-dependent and spans about 8 orders of magnitude.

  • Journal Ref
    New J. Phys. 11 (2009) 105027 [arXiv:0903.0189]

Cosmic-ray 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 sub-GeV dark matter masses. This specific plot assumes a constant spin-independent 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 tree-level 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 ]