micrOMEGAs is a code to compute dark matter observables in generic extensions of the standard model. This version of micrOMEGAs includes a generalization of the Boltzmann equations to take into account the possibility of two dark matter candidates. The modification of the relic density calculation to include interactions between the two DM sectors as well as semi-annihilation is presented. Both DM signals in direct and indirect detection are computed as well. An extension of the standard model with two scalar doublets and a singlet is used as an example.

With present and upcoming SUSY searches both directly, indirectly and at accelerators, the need for accurate calculations is large. We will here go through some of the tools available both from a dark matter point of view and at accelerators. For natural reasons, we will focus on public tools, even though there are some rather sophisticated private tools as well.

We calculate the relic density of dark matter in the MSSM with CP violation. We analyse various scenarios of neutralino annihilation: the cases of a bino, bino-wino and bino-Higgsino LSP, annihilation through Higgs, as well as sfermion coannihilation scenarios. Large phase effects are found, on the one hand due to shifts in the masses, on the other hand due to modifications of the couplings. Taking special care to disentangle the effects in masses and couplings, we demonstrate that the presence of CP phases can have a significant influence on the neutralino relic abundance. Typical variations in \Omega h^2 solely from modifications in the couplings are O(10%-100%), but can reach an order of magnitude in some cases.

We revisit the one-loop calculation of the annihilation of a pair of the lightest neutralinos into a pair of photons, a pair of gluons and also a Z photon final state. For the latter we have identified a new contribution that may not always be negligible. For all three processes we have conducted a tuned comparison with previous calculations for some characteristic scenarios. The approach to the very heavy higgsino and wino is studied and we argue how the full one-loop calculation should be matched into a more complete treatment that was presented recently for these extreme regimes. We also give a short description of the code that we exploited for the automatic calculation of one-loop cross sections in the MSSM that could apply both for observables at the colliders and for astrophysics or relic density calculations. In particular the automatic treatment of zero Gram determinants which appear in the latter applications is outlined. We also point out how generalised non-linear gauge fixing constraints can be exploited.

We present a code that calculates the relic density of the LSP in the MSSM. All tree-level processes for the annihilation of the LSP are included as well as all possible coannihilation processes with neutralinos, charginos, sleptons,squarks and gluinos. In all we have included over 2800 processes not counting charged conjugate states.The cross-sections extracted from CompHEP are calculated exactly. Relativistic formulae for the thermal average are used and care is taken to handle poles and thresholds by adopting specific integration routines. The input parameters can be either the soft SUSY parameters in a general MSSM or the five parameters of a SUGRA model. A link with ISASUGRA/ISAJET allows to calculate the parameters of the MSSM at the weak scale for an input at the GUT scale. The Higgs masses are calculated with FeynHiggsFast. Routines calculating g-2 and bsgamma are also included.