We present a comprehensive study of a model where the dark matter is composed of a singlet real scalar that couples to the Standard Model predominantly via a Yukawa interaction with a light quark and a colored vector-like fermion. A distinctive feature of this scenario is that thermal freeze-out in the early universe may be driven by annihilation both into gluon pairs at one-loop ($gg$) and by virtual internal Bremsstrahlung of a gluon ($q \bar{q} g$). Such a dark matter candidate may also be tested through direct and indirect detection and at the LHC; viable candidates have either a mass nearly degenerate with that of the fermionic mediator or a mass above about 2 TeV.

Gamma ray spectral features are of interest for indirect searches of dark matter (DM). Following Barger et al, we consider 3 simple scenarios of DM that annihilates into Standard Model (SM) fermion pairs. Scenario 1 is a Majorana DM candidate coupled to a charged scalar, scenario 2 is a Majorana DM coupled to a charged gauge boson and scenario 3 is a real scalar DM coupled a charged vector-like fermion. As shown by Barger et al, these 3 scenarios share precisely the same internal Bremsstrahlung spectral signature into gamma rays. Their phenomenology is however distinct. In particular for annihilation into light SM fermions, in the chiral limit, the 2-body annihilation cross section is p-wave suppressed for the Majorana candidates while it is d-wave suppressed for the real scalar. In the present work we study the annihilation into 2 gammas, showing that these three scenarios have distinct, and so potentially distinguishable, spectral signatures into gamma rays. In the case of the real scalar candidate we provide a new calculation of the amplitude for annihilation into 2 gammas.