 code to determine the particle spectrum and couplings, including the running Yukawa couplings, and compute the thermally averaged cross section for neutralino annihilation, including co-annihilation and solve the equation describing the evolution of the number density for the DM candidate. While the assumptions of the cMSSM are quite helpful for defining a set of benchmark points, the cMSSM is not representative of the generic MSSM, since it implies several mass relations, and its assumptions have no strong physics justification. Therefore, in studying the accuracy on $\Omega_{\chi} h^2$, the full set of MSSM parameters must be scanned in an uncorrelated way and the mass spectrum evaluated for each parameter set. A detailed study has recently been performed [@Baltz:2006fm]. I summarise here some of the findings, Table \[tab:dmsummary\] gives results for the neutralino relic density estimates in MSSM for the LHC, the ILC at 0.5 TeV and the ILC at 1 TeV.

The LCC1 point is in the bulk region and the model contains light sleptons, with masses just above that of the lightest neutralino. The most important annihilation reactions are those with t-channel slepton exchange. At the LHC, many of the SUSY spectrum parameters can be determined from kinematic constraints. At the ILC masses can be determined both by the two-body decay kinematics of the pair-produced SUSY particles and by dedicated threshold scans. Let us consider the two body decay of a scalar quark $\tilde q \to q \chi^0_1$. If the scalar quarks are pair produced $e^+e^- \to \tilde q \tilde q$, $E_{\tilde q} = E_{beam}$ and the $\chi^0_1$ escapes undetected, only the $q$ (and the $\bar q$) are observed in the detector. In a 1994 paper, J. Feng and D. Finnell [@Feng:1993sd] pointed out that the minimum and maximum energy of production for the quark can be related to the mass difference between the scalar quark $\tilde q$ and the $\chi^0_1$:

$$E_{max,~min} = \frac{E_{beam}}{2} \big( 1 \pm \sqrt{1 - \frac{m_{\tilde q}^2}{E_{beam}^2}}
\big) \big( 1 - \frac{
code to determine the particle spectrum and couplings, including the run Yukawa coupling, and compute the thermally averaged cross department for neutralino annihilation, including co - annihilation and clear the equality describing the evolution of the number density for the DM candidate. While the assumption of the cMSSM are quite helpful for defining a set of benchmark period, the cMSSM is not representative of the generic MSSM, since it implies several mass relations, and its assumptions have no potent physics justification. Therefore, in studying the accuracy on $ \Omega_{\chi } h^2 $, the wide set of MSSM parameters must be scanned in an uncorrelated manner and the mass spectrum evaluated for each parameter set. A detailed study has recently been perform   [ @Baltz:2006fm ]. I summarise here some of the findings, Table   \[tab: dmsummary\ ] gives results for the neutralino relic density estimates in MSSM for the LHC, the ILC at 0.5   TeV and the ILC at 1   TeV. 

 The LCC1 point is in the bulk region and the model contains light sleptons, with masses just above that of the lightest neutralino. The most important annihilation chemical reaction are those with t - channel slepton substitution. At the LHC, many of the SUSY spectrum parameters can be determined from kinematic constraints. At the ILC masses can be settle both by the two - body decay kinematics of the pair - produced SUSY particles and by consecrated threshold scans. Let us consider the two body decay of a scalar quark $ \tilde q \to q \chi^0_1$. If the scalar quarks are pair produced $ e^+e^- \to \tilde q \tilde q$, $ E_{\tilde q } = E_{beam}$ and the $ \chi^0_1 $ escapes undetected, only the $ q$ (and the $ \bar q$) are observed in the detector. In a 1994 newspaper, J.   Feng and D.   Finnell   [ @Feng:1993sd ] pointed out that the minimal and maximal department of energy of production for the quark can be related to the mass difference between the scalar quark $ \tilde q$ and the $ \chi^0_1 $: 

 $ $ E_{max,~min } = \frac{E_{beam}}{2 } \big (1 \pm \sqrt{1 - \frac{m_{\tilde q}^2}{E_{beam}^2 } } 
 \big) \big (1 - \frac {
 cofe to determine the partlcle spectrum anb couplmngs, indluding ghe running Yukawa couplings, abd conpute the thermally avdraged crlss sectuon hor neutralino aihihilatljn, ihgludiug co-annihilatiok and solve the equation gercxibing the evolution of the number dqnsity gog the DM candiqate. Rhils the assumptions of the cMSSM are quite ielpful for defoning a set of benchmark plintd, the cMSSM is not representajjve if the generkc MSSM, since it impligs several mass relations, and itr assbmptions hace no vtrong physmcs juftification. Bnerefose, in syudying the acguracb on $\Omega_{\chi} h^2$, the full vet of MSSM paramejers must te scanned in an unxoerelajed wdy ave tfe jaxs spectgum evaluated for each pqrameter set. A detaolqe study has rscentlr feen performed [@Baltz:2006fm]. I summarise here voms of the findings, Table \[rab:dmsummary\] gives redults for the neutralino relic density estimates in MSSM fmr thx UHC, tme IUX wt 0.5 TeV and the ILC at 1 TeV.

The LCC1 point is in ege bllk region and thc model contains logjt fleptons, with masses jhst above that of hhe ligrtest neutraligo. Tne most important annihilatuon reactionf are those with t-ckannel sleptun ecchanbe. At the LHC, many of tke SUSG spectrum oarametera can be determindd nrok kinemaukc constraints. At the ILC nasszs can bd dejermineq both by hhe two-body decay kinemahics lf the pair-pgoduced SUSY particles and by dxvicated thresnond vcans. Lec us cpnsider the tro body decay pf a scclar qjark $\tilde q \to q \chi^0_1$. If the scalar quarkv are pair prmduced $e^+q^- \to \tilee q \tiuae q$, $E_{\tilde q} = E_{beam}$ anb the $\chi^0_1$ escapes undetccted, knly the $q$ (and cke $\bar q$) are obsetvea ig nhe dqdector. In a 1994 papdr, G. Geng xnd D. Finnelo [@Weng:1993xd] pointed out that dhe jinimum and maximuk cnergy of productyon for the qiark can be relatef to vhe mars digfetence between the scalar quark $\filde q$ ajd bhe $\chi^0_1$:

$$E_{max,~mig} = \fvac{E_{feam}}{2} \big( 1 \'m \sqrt{1 - \frac{m_{\tilde q}^2}{E_{beam}^2}}
\big) \big( 1 - \frac{
code to determine the particle spectrum and the Yukawa couplings, compute the thermally annihilation, co-annihilation and solve equation describing the of the number density for the candidate. While the assumptions of the cMSSM are quite helpful for defining a of benchmark points, the cMSSM is not representative of the generic MSSM, since implies mass and assumptions have no strong physics justification. Therefore, in studying the accuracy on $\Omega_{\chi} h^2$, the full of MSSM parameters must be scanned in an way and the mass evaluated for each parameter set. detailed has recently performed I here some of findings, Table \[tab:dmsummary\] gives results for the neutralino relic density estimates in MSSM for the LHC, the at 0.5 the ILC 1 The point is in region and the model contains light just above that of the lightest neutralino. The important annihilation are those with t-channel slepton exchange. the LHC, many of the SUSY spectrum parameters be determined from kinematic constraints. At the ILC masses can be determined both by the kinematics of the pair-produced particles and by threshold Let consider two body of a scalar quark $\tilde q \to q \chi^0_1$. If the quarks are pair produced $e^+e^- \to \tilde q \tilde q$, = and the $\chi^0_1$ undetected, only the $q$ the q$) are observed in In 1994 and Finnell pointed out that the and maximum energy of production the quark can be between the scalar quark $\tilde q$ and the $$E_{max,~min} = \frac{E_{beam}}{2} \big( 1 \pm \sqrt{1 \frac{m_{\tilde q}^2}{E_{beam}^2}} \big) \big( 1 - \frac{
 code to determine the particlE spectrum aNd couPliNgs, InCludIng tHe running YukawA CoupLings, and compute the therMally AvERageD CrOss seCtion foR NeUTRalInO aNniHiLAtIon, inCluDing co-aNnihilatioN anD sOlve the equatIOn Describing The Evolution of tHe nUmber dEnSitY For thE DM CandiDate. WhILe the aSsumptionS oF The cMSsm are quiTE HeLpfuL for defining a set oF BeNChmark points, thE cMSSM Is NOt REPreSenTative of thE gEneriC mSSM, sinCE iT IMPliES several mass rElations, and ITs aSsumptIoNs hAVe no stRong pHySIcs JustificatiOn. ThErefore, in StudyiNG the accURacy on $\OMega_{\chI} h^2$, tHe fUll sET oF MsSM PaRAmeTErS muST be Scanned iN aN uNcorrElatED WAY and The Mass SpectRum evaluated fOr eAch pARamEter sEt. A deTailEd Study Has recEntly BeEn performed [@BaltZ:2006fm]. I Summarise HerE sOme Of The fiNDings, TAblE \[taB:dmsummAry\] giveS ResUlTS FOr The neutralino relic DeNSItY estimatEs in MSsm fOr THe LHC, the iLc at 0.5 teV aND The ILc at 1 TEv.

THe LCC1 poiNt is in THe BuLk regioN aNd the mOdEl cOntAins lIGht sLeptonS, with masSes juST above that of thE Lightest neutrALiNO. thE Most ImpOrtant annihIlatIOn reActiONs Are THose wIth t-cHaNNeL Slepton exchange. At thE LhC, many Of the sUSY spectrum pArameters cAN BE determiNed fROm KInematic constrAints. at the ILC maSSes can be DeterMined botH by the two-BODy decay kIneMatIcs Of tHE PaIr-produced SUSy PArtiClEs and by DedIcated tHreShoLd sCanS. LEt us consiDer the twO bOdY dEcAy oF a scaLAr quark $\tIlDe q \To Q \chI^0_1$. If thE Scalar QuarkS are PaIr PRodUced $e^+e^- \tO \TiLDE q \tiLdE q$, $e_{\tilDe q} = e_{bEam}$ anD the $\CHi^0_1$ eScapes uNdetected, OnlY The $q$ (AnD tHe $\bar q$) aRe observed in tHe Detector. In A 1994 pApeR, J. Feng AND D. FinnelL [@Feng:1993sd] pointed out that thE Minimum And MaximUm enErgy of proDucTion foR thE Quark cAn be reLated To The MASs difFEReNce BeTween the scALAr qUark $\tIlDe q$ aNd the $\chI^0_1$:

$$E_{max,~min} = \frac{E_{beam}}{2} \bIG( 1 \pm \Sqrt{1 - \frac{m_{\tildE q}^2}{E_{Beam}^2}}
\BIG) \bIg( 1 - \fRAc{
 code to determine the par ticle spec trumand  co up ling s, i ncluding the r u nnin g Yukawa couplings, an d com pu t e th e  t herma lly ave r ag e d  cr os ssec ti o nfor n eut ralinoannihilati on,  i ncluding co- a nn ihilationand  solve the e qua tion d es cri b ing t heevolu tion o f  the n umber den si t y fort he DM c a n di date . While the assum p ti o ns of the cMSS M arequ i te h elp ful  for defin in g a s e t of be n ch m a r k p o ints, the cMS SM is not r e pre sentat iv e o f  the g eneri cM SSM , since itimpl ies sever al mas s  relati o ns, and  its a ssu mpt ions ha ve  no  s t ron g  p hys i csjustific at io n. Th eref o r e ,  instu dyin g the  accuracy on$\O mega _ {\c hi} h ^2$,thefu ll se t of M SSM p ar ameters must be  sca nned in a n u nc orr el atedw ay and  th e m ass spe ctrum e v alu at e d fo r each parameter s et . Adetailed  study ha sr ecentlybe enperf o r med [ @Bal t z: 2006fm].  I sum m ar is e hereso me ofth e f ind ings, Tabl e \[ta b:dmsumm ary\] gives resultsf or the neutra l in o re l ic d ens ity estimat es i n  MSS M fo r  t heL HC, t he IL Ca t0 .5 TeV and the ILCat  1 TeV .

Th e LCC1 pointis in theb u l k region  and th e  model contain s lig ht slepton s , with m asses  just ab ove thato f  the lig hte stneu tra l i no . The most im p o rtan tannihil ati on reac tio nsare  th os e with t- channelsl ep to nexc hange .  At theLH C,ma nyof th e  SUSYspect rumpa ra m ete rs canb ed e term in ed  fro m k in emati c co n str aints.At the IL C m a sses  c an  be det ermined bothby  the two-b od y d ecay k i n ematicsof the pair-produced SU S Y parti cle s and  bydedicated  th reshol d s c ans. L et usconsi de r t h e  twob o dy  de ca y of a sca l a r q uark$\ tild e q \to  q \chi^0_1$. If t h e s calar quarksare  pai r pr odu c ed $e^ +e ^ - \ t o  \tilde q \tild e q$, $E_{ \t i ld e q} = E_{ b eam }$  and th e $\chi ^0_1$ escapes  undetect ed, onlyth e $q $ (an d the $\ba r q$) ar e observe d  in t h edetec tor . In a  1 994  pape r, J.F eng  andD. Fin ne ll [@F eng:1 99 3sd] poi nted out that the minim um and  maxi mum  energy o f p r odu ction for  the  quark can  be  re lated  to the m assd if fer e nce b etwe e n the sca l ar  qu a r k$\tilde q$a n d  th e $\c hi^ 0 _1$:
$$E_ {max,~min} = \fra c {E_{beam}}{2}\big ( 1 \ pm\ sqrt {1  - \frac{m_{\t ild eq } ^2}{E_{b ea m}^2}}
\big ) \big(1-  \fra c{
 code_to determine_the particle spectrum and_couplings, including_the_running Yukawa_couplings,_and compute the_thermally averaged cross_section for neutralino annihilation,_including co-annihilation and_solve_the equation describing the evolution of the number density for the DM candidate. While_the_assumptions of_the_cMSSM_are quite helpful for defining_a set of benchmark points,_the cMSSM_is not representative of the generic MSSM, since_it_implies several mass_relations, and its assumptions have no strong physics justification._Therefore, in studying the accuracy on_$\Omega_{\chi} h^2$, the_full_set_of MSSM parameters must_be scanned in an uncorrelated way_and the mass spectrum evaluated for_each parameter set. A detailed study has_recently been performed [@Baltz:2006fm]. I summarise here_some of the findings, Table \[tab:dmsummary\]_gives results_for the neutralino relic density_estimates in MSSM_for the_LHC, the ILC_at 0.5 TeV and the ILC at_1 TeV.

The LCC1 point_is in the bulk region and_the_model contains light_sleptons,_with_masses just_above that of_the_lightest neutralino._The_most important annihilation reactions are those_with_t-channel slepton exchange. At the LHC, many_of the SUSY spectrum_parameters_can be determined from_kinematic constraints. At the ILC_masses can be determined both by_the two-body_decay kinematics_of the pair-produced SUSY particles and by dedicated threshold scans. Let_us consider the two body decay_of a scalar quark_$\tilde q_\to_q \chi^0_1$. If_the_scalar quarks_are pair produced $e^+e^- \to \tilde q_\tilde q$,_$E_{\tilde q} = E_{beam}$ and the_$\chi^0_1$ escapes undetected, only_the_$q$ (and the $\bar q$) are_observed in the detector. In a_1994 paper, J. Feng and D. Finnell [@Feng:1993sd]_pointed_out_that the minimum and maximum_energy of production for the quark_can be related_to the mass difference between the scalar_quark_$\tilde q$ and the $\chi^0_1$:

$$E_{max,~min} =_\frac{E_{beam}}{2}_\big( 1 \pm \sqrt{1 - \frac{m_{\tilde_q}^2}{E_{beam}^2}}
\big)_\big(_1 - \frac{