DESY 12-205, LAPTH-051/12, LPSC-12-341, MS-TP-12-17, KA-TP-42-2012 Neutralino-stop co-annihilation into electroweak gauge and Higgs bosons at one loop J. Harz,1, ∗ B. Herrmann,2, † M. Klasen,3, ‡ K. Kovaˇr´ık,4, § and Q. Le Boulc’h5, ¶ 1Deutsches Elektronen-Synchrotron (DESY), Notkestraße 85, D-22607 Hamburg, Germany 2LAPTh, Universit´ede Savoie / CNRS, 9 Chemin de Bellevue, B.P. 110, F-74941 Annecy-le-Vieux, France 3Institut f¨ur Theoretische Physik, Westf¨alische Wilhelms-Universit¨at M¨unster, Wilhelm-Klemm-Straße 9, D-48149 M¨unster, Germany 4Institute for Theoretical Physics, Karlsruhe Institute of Technology, D-76128 Karlsruhe, Germany 5Laboratoire de Physique Subatomique et de Cosmologie, Universit´eJoseph Fourier/CNRS-IN2P3/INPG, 53 Rue des Martyrs, F-38026 Grenoble, France (Dated: January 15, 2013) We compute the full O(αs) supersymmetric QCD corrections for neutralino-stop co-annihilation into electroweak gauge and Higgs bosons in the Minimal Supersymmetric Standard Model (MSSM). We show that these annihilation channels are phenomenologically relevant within the so-called phe- nomenological MSSM, in particular in the light of the observation of a Higgs-like particle with a mass of about 126 GeV at the LHC. We present in detail our calculation, including the renormalization scheme, the infrared treatment, and the kinematical subtleties to be addressed. Numerical results for the co-annihilation cross sections and the predicted neutralino relic density are presented. We demonstrate that the impact of including the corrections on the cosmologically preferred region of parameter space is larger than the current experimental uncertainty from WMAP data. PACS numbers: 12.38.Bx,12.60.Jv,95.30.Cq,95.35.+d I. INTRODUCTION where the first term on the right-hand-side containing the Hubble parameter H stands for the dilution of dark mat- Many astrophysical observations over a wide range of ter due to the expansion of the Universe. The remaining length scales provide today convincing evidence of a siz- terms reduce (increase) the number of dark matter par- able Cold Dark Matter (CDM) component in the Uni- ticles by their annihilation (creation) in collisions with other particles. The term σannv denotes the thermal verse. The most recent measurements of the WMAP h i satellite in combination with baryonic acoustic oscilla- average of the annihilation cross section of the neutralino tion and supernova data [1] constrain the dark matter −E1/T −E2/T 3 3 relic density to the very precise value of σv e e d p1d p2 σannv = (1.3) −E1/T −E2/T 3 3 2 h i e e d p1d p2 ΩCDMh = 0.1126 0.0036 (1.1) R ± 1 ∞ R 2 √ √ = 4 2 σ s 4m sK1 s/T ds , at 1σ confidence level, where h denotes the present Hub- 8m TK m/T 4 2 − −1 −1 2 Z m ble expansion rate H0 in units of 100 km s Mpc . The leading candidate for dark matter is a weakly interacting where Ki are the modified Bessel functions of the second massive particle (WIMP). Unfortunately, the Standard kind of order i (for details see Ref. [2]). Model does not contain any candidates that would be Here we will consider the case when heavier, unstable compatible with the properties of cold dark matter, the supersymmetric particles survive in the Universe for suf- neutrinos being too light. ficient time to affect the relic density of the dark matter In contrast, the Minimal Supersymmetric Standard arXiv:1212.5241v2 [hep-ph] 11 Jan 2013 particle. In this case, Eq. (1.2) has to be modified to ac- Model (MSSM) with conserved R-parity contains with count for the interactions between all particles and solve the lightest neutralino, denotedχ ˜0, a stable WIMP. Since 1 a system of Boltzmann equations for number densities ni supersymmetry offers many other theoretical and phe- for each surviving particle species, nomenological advantages, the neutralino is by far the most studied dark matter candidate. The time evolution dni eq eq of its number density n is described by the Boltzmann = 3Hni σij vij ninj n n , (1.4) χ dt − − h i − i j equation where σ σ(χ χ X) is the cross section of the dnχ 2 eq 2 ij i j = 3Hn σannv n (n ) , (1.2) ≡ → dt − χ − h i χ − χ annihilation of particle i with particle j, and vij is their relative velocity. As all heavier particles eventually decay into the dark matter particle, the relevant quantity is the ∗ [email protected] total number density nχ = i ni, and we can reformulate † [email protected] the problem into solving a single Boltzmann equation ‡ [email protected] similar to Eq. (1.2) for theP total number density with the § [email protected] annihilation cross section σannv replaced by an effective ¶ [email protected] h i cross section σeff v . More precisely, this cross section is h i 2 given by the remaining particles and their couplings to the neu- tralino. In the MSSM, the mass and the couplings of the eq neq ni j neutralino as well as any other relevant couplings and σeff v = σij vij eq eq , (1.5) h i nχ nχ masses are typically obtained using a dedicated spectrum i,j X calculator (see, e.g., Ref. [7]), which evolves all parame- where the sum runs over all MSSM particles i and j (for ters down from a grand unification scale and calculates a detailed discussion see Refs. [3, 4]). The ratio of their the masses and couplings for all particles at the weak eq respective number density in thermal equilibrium, ni , scale. Different treatments of the radiative corrections eq and the number density of the neutralino, nχ , at the for masses and couplings as well as different implemen- temperature T is Boltzmann suppressed, tations of the renormalization group equations in various MSSM spectrum calculations can lead to differences in neq m m i exp i − χ , (1.6) the predictions for the relic density and thus in the pre- neq ∼ − T χ ferred/excluded regions of the MSSM parameter space (for details see, e.g., Ref. [8]). so that only particles, whose masses are almost degener- ate with the one of lightest neutralino, can give sizable The uncertainty which we will address in this paper contributions. In the MSSM, typical examples of relevant does not fall into either of the above mentioned cate- co-annihilations are those of the neutralino with the light- gories, but concerns the precision, with which the (co-) est slepton or squark, or with another gaugino. More- annihilation cross sections in Eq. (1.5) are computed. over, pair annihilations of the next-to-lightest superpart- The cross sections in public dark matter tools such as ner can be non-negligible. Having solved the Boltzmann DarkSUSY [9] or micrOMEGAs [10] are implemented us- equation numerically, the relic density is finally obtained ing only an effective tree-level calculation. It is, how- through ever, well known that higher-order corrections, particu- larly those involving the strong coupling constant, can 2 mχnχ Ωχh = , (1.7) have a sizable impact on such processes. The impact of ρcrit next-to-leading order corrections to neutralino annihila- tion on the neutralino relic density has been discussed where nχ is the current number density of the neutralino in several previous analyses. SUSY-QCD corrections to and ρcrit is the critical density of the Universe. Com- paring the predicted value obtained by solving the Boltz- neutralino pair annihilation into quark-antiquark pairs mann equation to the observational limits in Eq. (1.1) al- have been studied in Refs. [11–13], while the correspond- lows one to identify cosmologically (dis-)favored regions ing electroweak corrections have been evaluated in Refs. of the MSSM parameter space and thus to obtain im- [14–16]. The authors of Refs. [15, 16] have also discussed portant information that is complementary to collider the case of co-annihilation of a neutralino with another searches and precision measurements. gaugino. Further studies rely on effective coupling ap- The procedure described above is unfortunately sub- proaches in order to capture certain classes of corrections ject to several uncertainties. The first source of uncer- to neutralino pair annihilation and co-annihilation with tainty lies in the extraction of the relic density of CDM a tau slepton [17, 18]. All these analyses show that radia- from cosmological data as given in Eq. (1.1). The ex- tive corrections are not negligible in the context of relic traction is based on a simple cosmological model, the density calculations, the impact of the corrections being ΛCDM model, which uses a minimal set of six parame- larger than the experimental uncertainty from WMAP ters to fit the available cosmological data and bases its in many regions of parameter space. With the Planck conclusions on the Standard Model of cosmology [1]. It satellite data providing more precise cosmological mea- has been shown that changing either the number of free surements in the very near future, it becomes even more parameters of the model used to fit the cosmological data pressing that theoretical predictions match the experi- [5] or modifying the assumptions contained in the Stan- mental precision. dard Model of cosmology (e.g., altering the expansion The important case of SUSY-QCD corrections to co- rate in the primordial Universe or later, but still before annihilation of a neutralino with a scalar top has so far Big Bang Nucleosynthesis [6]), may change the extracted only been considered in Ref. [19]. This study concerns 2 central value of ΩCDMh along with the confidence levels.