arXiv:1212.5241v2 [hep-ph] 11 Jan 2013 oe MS)wt conserved with (MSSM) Model the matter, be dark light. would cold too that of being candidates neutrinos properties the any with contain compatible not Standard the does Unfortunately, Model interacting weakly (WIMP). a particle is matter massive dark for candidate leading equation h ihetnurln,dntd˜ denoted neutralino, lightest the otsuiddr atrcniae h ieevolution density the time number The far its candidate. of by matter is dark phe- studied neutralino and most the theoretical advantages, other nomenological many offers supersymmetry l xaso rate expansion ble t1 at inadsproadt 1 osri h akmatter of value dark precise the very constrain the to [1] oscilla- WMAP density data acoustic relic the supernova baryonic of and with measurements Uni- tion combination recent the in in most satellite component The (CDM) siz- Matter verse. a Dark of evidence Cold convincing able today provide scales length ¶ ∗ § ‡ † etaioso oanhlto noeetoekguea gauge electroweak into co-annihilation Neutralino-stop [email protected] [email protected] [email protected] [email protected] [email protected] ncnrs,teMnmlSprymti Standard Supersymmetric Minimal the contrast, In ayatohsclosrain vrawd ag of range wide a over observations astrophysical Many σ d d ofiec ee,where level, confidence n t χ 2 AT,Uiested aoe/CR,9Cei eBleu,B. Bellevue, de Chemin 9 CNRS, / Savoie Universit´e de LAPTh, = 4 ASnmes 12.38.Bx,12.60.Jv,95.30.Cq,95.35.+d unc numbers: experimental PACS current on the corrections than the larger ne is including predicted space of parameter the impact calculati and the our sections that subtle detail demonstrate cross kinematical in co-annihilation the present the and We for treatment, LHC. infrared the the at scheme, GeV observ the 126 of light about the of in phenomenologi particular are in channels Super MSSM, nomenological Minimal annihilation the these in that show bosons We Higgs and gauge electroweak into nttt o hoeia hsc,KrsueIsiueof Institute Karlsruhe Physics, Theoretical for Institute Ω nvri´ oehFuirCR-NP/NG 3RedsMa des Rue 53 Fourier/CNRS-IN2P3/INPG, Universit´e Joseph − ecmuetefl O( full the compute We CDM 1 .INTRODUCTION I. 3 etce lkrnnSnhorn(EY,Ntetae85 Notkestraße (DESY), Elektronen-Synchrotron Deutsches .Harz, J. Hn H h 0 2 χ 3 nuiso 0 ms km 100 of units in n ntttfu hoeicePyi,Westf¨alische Wilhelm Physik, f¨ur Theoretische Institut h − 0 = χ 1, sdsrbdb h Boltzmann the by described is σ ∗ ann . 1126 .Herrmann, B. h v eoe h rsn Hub- present the denotes 5 i χ ihl-lm-tae9 -84 ¨ntr Germany M¨unster, D-48149 9, Wilhelm-Klemm-Straße aoaor ePyiu uaoiu td Cosmologie, de et Subatomique Physique de Laboratoire 1 0 R n ± tbeWM.Since WIMP. stable a , prt otiswith contains -parity χ 2 − 0 α . s 06(1.1) 0036 ( uesmercQDcretosfrnurln-tpco-a neutralino-stop for corrections QCD supersymmetric ) n − χ eq 1 2, ) Mpc 2 † EY1-0,LPH011,LS-231 ST-21,KA- MS-TP-12-17, LPSC-12-341, LAPTH-051/12, 12-205, DESY Dtd aur 5 2013) 15, January (Dated: , .Klasen, M. − 1 The . (1.2) 3, ‡ where ytmo otmn qain o ubrdensities number for ac- solve equations to and Boltzmann modified particles of be all system to between a matter interactions has dark the (1.2) the for Eq. count of case, density this relic In the particle. suf- affect for to Universe time the ficient in survive particles supersymmetric ido order of kind where niiaino particle of annihilation hr h rttr ntergthn-iecnann the containing right-hand-side parameter the Hubble on term first the where o ahsriigpril species, particle surviving each for notedr atrpril,terlvn uniyi the is quantity relevant density number the total particle, matter decay dark eventually the particles into heavier all As velocity. relative ilsb hi niiain(rain ncliin with term collisions The in (creation) par- particles. matter annihilation other dark their remaining of by The number ticles the Universe. the (increase) of reduce expansion terms the to due ter vrg fteanhlto rs eto fteneutralino the of section cross annihilation the of average iia oE.(.)frtettlnme est ihthe with equation density number Boltzmann total section the cross single annihilation for (1.2) a Eq. solving to similar into problem the rs section cross = h .Kovaˇr´ık,K. σ eew ilcnie h aewe eve,unstable heavier, when case the consider will we Here ann 8 d d ehooy -62 alrh,Germany Karlsruhe, D-76128 Technology, m n t taiorlcdniyaepeetd We presented. are density relic utralino v K i σ to faHgslk atcewt mass a with particle Higgs-like a of ation 4 ist eadesd ueia results Numerical addressed. be to ties i ranyfo MPdata. WMAP from ertainty K T al eeatwti h ocle phe- so-called the within relevant cally ij i h omlgclypeerdrgo of region preferred cosmologically the = = .10 -44 neyl-iu,France Annecy-le-Vieux, F-74941 110, P. ymti tnadMdl(MSSM). Model Standard symmetric r h oie eslfntoso h second the of functions Bessel modified the are ≡ 1 2 R 2 -20 abr,Germany Hamburg, D-22607 , n nldn h renormalization the including on, − tr,F306Geol,France Grenoble, F-38026 rtyrs, h 4, R s-Universit¨at M¨unster, i ve σv σ m/T σ 3 § frdtisseRf [2]). Ref. see details (for ( eff e χ Hn − dHgsbsn toeloop one at bosons Higgs nd n .L Boulc’h Le Q. and v i E − χ H i 1 i E oepeiey hscosscinis section cross this precisely, More . j /T Z h − n 1 tnsfrtedlto fdr mat- dark of dilution the for stands 4 /T → χ e m ∞ i − σ = 2 e E ihparticle with − ij X h σ σ 2 P E v /T stecosscino the of section cross the is ) h ann ij s 2 σ /T i − d i ann n v 3 d i i p 4 n n ecnreformulate can we and , 3 v 1 m elcdb neffective an by replaced i p nnihilation i d n 1 2 3 j d eoe h thermal the denotes 5, p 3 − √ 2 p ¶ j K s 2 n and , i eq 1 n j eq √ TP-42-2012 v s/T ij , stheir is (1.4) (1.3) s, ds n 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. the very specific cases of co-annihilation of a bino-like The second source of uncertainty in identifying (dis-) neutralino with a right-handed stop into a top quark and favored regions of the MSSM parameter space is con- a gluon as well as into a bottom quark and a W -boson. nected to the calculation of the essential parameters such However, depending on the considered region of param- as masses and couplings of supersymmetric particles. As eter space, many other final states, including those with the relic density is very sensitive to the mass of the neu- other electroweak gauge and Higgs bosons, can become tralino (see Eq. (1.7)), any uncertainty in the calcula- dominant. Moreover, in realistic supersymmetric sce- tion of its mass directly translates into an uncertainty on narios, helicity mixing in the stop sector is usually non- the calculated relic density. Moreover, the relic density negligible, as is the mixing of bino, wino, and higgsino also strongly depends on the (co-)annihilation cross sec- components in the lightest neutralino, which strongly tion, which in turn crucially depends on the masses of influences its couplings and preferred (co-)annihilation 3
q’ q’ χ0 χ0 0 q’ ˜ n ˜ n χ˜ n , q˜ j χ˜ k q φ φ φ
q˜ i q˜ i q˜ i q’ q’ χ˜ 0 χ˜ 0 χ˜ 0 q’ n n n q˜ , χ˜ q j k V V V q˜ i q˜ i q˜ i
FIG. 1. Leading-order Feynman diagrams for neutralino-squark co-annihilation into a quark and a Higgs boson (φ = h0,H0,A0,H±) or an electroweak gauge boson (V = γ, Z0,W ±). The u-channel is absent for a photon in the final state. channels. Therefore, we extend in this paper the analysis an NLSP is the scalar top, whose chirality eigenstates of QCD and SUSY-QCD corrections to co-annihilation of can mix significantly, e.g. when the trilinear coupling At a neutralino with a stop by computing the general case of becomes large, and which can then have a lower mass neutralino-stop co-annihilation into a quark and a Higgs eigenstate that is almost mass-degenerate with the light- or an electroweak vector boson. The paper is organized est neutralino [20]. as follows: In Sec. II, we first discuss the phenomenol- There is ample motivation for a light scalar top. First, ogy of neutralino-stop co-annihilation in the MSSM. We a light stop is a necessary ingredient to achieve elec- then describe in detail the calculation of the radiative troweak baryogenesis in the MSSM [21]. Second, “natu- corrections to the relevant processes in Sec. III. Numeri- ral” SUSY models [22, 23] require a light third genera- cal results for annihilation cross sections and dark matter tion of sfermions in order to reduce fine-tuning and stay relic densities in typical MSSM benchmark scenarios are compatible with experimental constraints at the same presented in Sec. IV, and our conclusions are given in time. This is due to the fact that the mass degeneracy Sec. V. between the lightest neutralino and NLSP weakens the LHC exclusion potential on the third-generation squark masses, since this degeneracy results in events with soft II. PHENOMENOLOGY OF jets [24, 25]. Third, interpreting the new boson with a NEUTRALINO-STOP CO-ANNIHILATION mass of about 126 GeV observed recently at the LHC [26– 28] as a light CP-even Higgs boson (h0) implies within the As discussed in Sec. I, the co-annihilation of the next- MSSM a particular choice of parameters in the stop and to-lightest supersymmetric particle (NLSP) with the sbottom sector [29]. The reason is that in the MSSM the lightest neutralino can in certain regions of the MSSM lightest Higgs boson mass receives a large contribution parameter space become dominant and lead to a relic from a loop containing scalar tops. The leading contri- density that is compatible with the observational limit bution to the mass coming from this loop together with of Eq. (1.1). A particularly important example of such the tree-level contribution can be expressed as [30, 31]
2 4 2 2 2 2 2 2 3g mt MSUSY Xt Xt m 0 = m cos 2β + log + 1 , (2.1) h Z 8π2m2 m2 M 2 − 12 M 2 W t SUSY SUSY
with Xt = At µ/ tan β and MSUSY = √mt˜1 mt˜2 . The man diagrams are depicted in Fig. 1. These processes maximal contribution− from stop mixing is then obtained compete with all other possible (co-)annihilation chan- for X √6MSUSY, which favors a sizable trilinear cou- nels of the lightest neutralino and in certain cases also | t|∼ pling At and consequently a rather light stop. with stop pair annihilation. At tree level, the co-annihilation of a neutralino and In order to quantify the relative importance of the pro- a stop into final states containing a quark and an elec- cesses in Fig. 1, we have performed a random scan in the troweak gauge or Higgs boson is mediated either by an phenomenological MSSM. In the following we describe s-channel quark, a t-channel squark, or a u-channel neu- the settings and discuss in detail the results of our scan. tralino or chargino exchange. The corresponding Feyn- According to the SPA convention [32] the soft-breaking 4
FIG. 2. Relative contributions of the neutralino-stop co-annihilation channels for the generated parameter points as a function of the input parameters M1, Mq˜3 , Tt, and tan β before (top) and after (bottom) applying the selection cuts of Eq. (2.3). Shown are the contributions from th0 (red), tg (green), tZ0 (blue), tH0 (yellow), bW + (cyan), tA0 (brown), bH+ (pink), and tγ (gray)
final states. The parameters M1, Mq˜3 , and Tt are given in GeV. parameters are defined at the scale Q = 1 TeV. We ter points within the following ranges for the eight input have made a few simplifying assumptions, which bring parameters: the number of parameters down to eight. In the squark 500 GeV Mq˜1,2 4000 GeV, sector, we use a common mass parameter Mq˜1,2 for the ≤ ≤ 100 GeV M˜3 2500 GeV, squarks of the first and second generation, leaving the ≤ q ≤ common mass parameter Mq˜3 for the left- and right- 500 GeV Mℓ˜ 4000 GeV, handed squarks of the third generation independent. In ≤ ≤ Tt 5000 GeV, contrast, the slepton sector is characterized by a single | | ≤ 200 GeV M1 1000 GeV, (2.2) mass parameter Mℓ˜ for all three generations. All trilin- ≤ ≤ ear couplings are set to zero except for the A in the stop 100 GeV mA 2000 GeV, t ≤ ≤ sector, which enters our calculations through the relation µ 3000 GeV, | | ≤ Tt = YtAt with the top Yukawa coupling Yt. All gaugino 2 tan β 50. masses are defined through the bino mass parameter M1. ≤ ≤ The wino and gluino masses are then fixed by the rela- For each set of parameters, the physical mass spectrum tion 2M1 = M2 = M3/3, which is deduced from gaugino and the related mixing matrices have been obtained us- mass unification at the GUT scale. Finally, the Higgs sec- ing SPheno [7] (version 3.2.1). The neutralino relic den- 2 tor is specified by the pole mass mA of the pseudoscalar sity Ωχh as well as the contributions from the individ- Higgs boson, the higgsino mass parameter µ, as well as ual (co-)annihilation channels have been computed using the ratio tan β of the two vacuum expectation values of micrOMEGAs (version 2.4.1). For the numerical values the Higgs doublets. In order to explore the parameter of the Standard Model parameters we refer the reader to space, we have randomly generated 1.2 million parame- Ref. [33]. For a substantial number of these scenarios, co-annihilation of the lightest neutralino with a scalar 5
M1 M˜ M˜ M˜ T m µ β m 0 m˜ m 0 m 0 q1,2 q3 ℓ t A tan χ˜1 t1 h H I 306.9 2037.7 709.7 1499.3 1806.5 1495.6 2616.1 9.0 307.1 350.0 124.43 1530.72 II 470.6 1261.2 905.3 1963.2 1514.8 1343.1 725.9 18.3 467.3 509.4 124.06 1342.77 III 314.4 2870.5 763.6 2417.7 1877.5 386.0 2301.5 10.3 316.5 371.9 123.43 367.45
TABLE I. Three characteristic scenarios chosen in the pMSSM, which will be considered in this study. Given are the input m 0 m˜ parameters as described in the text, the lightest neutralino mass χ˜1 , the lightest stop mass t1 , and the masses of the light and heavy CP-even Higgs bosons mh0 and mH0 . All values except for tan β are given in GeV.
2 0 0 0 0 0 0 0 + Ωχh χ˜1t˜1 → th χ˜1t˜1 → tH χ˜1t˜1 → tZ χ˜1t˜1 → bW Sum I 0.114 38.5% – 3.4% 5.9% 47.8% II 0.116 24.6% – 10.7% 3.4% 38.7% III 0.111 14.2% 20.7% 1.2% 2.1% 38.2%
TABLE II. Neutralino relic density and relative contributions of neutralino-stop co-annihilation into a quark and a Higgs or electroweak gauge boson for the characteristic scenarios of Tab. I. The last column gives the sum of the listed contributions. top plays an important role. This can be seen in the ers is that the lightest neutralino and the lightest scalar upper part of Fig. 2, where we show the relative contri- top are almost mass degenerate. This is reflected in the butions of the different final states channels to the total left and left-center columns of Fig. 2 where we can ob- (co-)annihilation cross section as a function on the phe- serve a strikingly similar dependence of fraction of co- nomenologically most relevant input parameters. annihilation processes on the gaugino mass parameter
Experimentally viable scenarios have to satisfy a num- M1 and the third-generation squark mass parameter Mq˜3 , ber of additional constraints. We therefore impose the which are largely responsible for the masses of neutrali- following cuts on the neutralino relic density, the mass nos and squarks of the third generation. For large values of the lightest Higgs boson, and the inclusive branching of both parameters, co-annihilations cease to to be im- ratio of the most sensitive B-meson decay, b sγ: portant and annihilations of stops take their place as the → dominating contribution of the total cross section. 0.0946 Ω h2 0.1306, ≤ χ ≤ In the right-center panel of Fig. 2 one can notice the in- 120 GeV m 0 130 GeV, (2.3) teresting feature that after applying the cut on the light- ≤ h ≤ 2.77 10−4 BR(b sγ) 4.33 10−4. est Higgs boson mass, large values of Tt are preferred and · ≤ → ≤ · the initially rather important percentage of top-gluon fi- The first cut selects the points which match the observed nal states is reduced. This is driven by the fact that relic density of Eq. (1.1) within a 5σ confidence inter- in contrast to the top-gluon final state, the Higgs bo- val. The second limit corresponds to a very conserva- son mass prefers a sizable trilinear coupling. Moreover, tive mass range for the new boson observed at the LHC positive values for Tt are slightly preferred, since they al- [26, 27]. Note that the theoretical uncertainty on the cal- low a better maximization of the Higgs boson mass [29]. culation of the lightest Higgs boson mass within SPheno This is well visible in the center-right column of Fig. 2, is estimated to be about 3 GeV [7]. Finally, the lim- where after applying the cuts two distinct ranges for the its on the branching ratio of b sγ correspond to a trilinear coupling parameter Tt can be observed. These 3σ interval around the observed value→ of BR(b sγ)= large values also enhance the Higgs-squark-squark cou- (3.55 0.26) 10−4 [34]. The points selected in→ this way pling, which is present in the t-channel of the th0 final are depicted± in· the lower part of Fig. 2, where we show state. Accordingly, this changes the relative importance again the relative contribution of the different neutralino- of the squark exchange with respect to the two other di- stop co-annihilation channels. Applying the experimen- agrams (quark or neutralino exchange, see Fig. 1). The tal cuts described above reduces the density of the points, t-channel enhancement also leads to an almost univer- but does not significantly change the shape of the distri- sal dominance of co-annihilation into Higgs final states butions. As can be seen, the statistically most impor- in the scenarios considered here. In other words, the tant final state is a top quark together with a light Higgs same mechanism which drives the mass of the lightest boson, followed by top quark and a gluon, a heavy CP- Higgs boson to the observed value through important even Higgs boson, or a Z-boson. Comparable in size to stop-loop contributions is responsible for the increase of the latter channel is the co-annihilation into a bottom neutralino-stop co-annihilation into the lightest Higgs bo- quark and a W -boson, whereas final states including a son together with a top quark. pseudoscalar Higgs boson, a charged Higgs boson, or a A similar connection between parameters that we men- photon are less important. tioned above for M1 and Mq˜3 is found for Tt and the Higgs One viable option how to satisfy the relic density parameter µ. After the cuts, large and positive values of bound and respect current exclusion limits from collid- µ are preferred, which also enhances the Higgs-sfermion- 6
2.5 ̃ ̃ → (Scenario I) ̃ ̃ → (Scenario III) 3.0 2.0 ) ) − − 1.5 2.0