SUSY Dark Matter Loops and Precision from Particle Physics Fawzi BOUDJEMA LAPTH, CNRS, France in collaboration with Andrei Semenov and David Temes in parts with Ben Allanach, Genevi` eve B´ elanger and Sacha Pukhov RADCOR05 F. BOUDJEMA, SUSY Dark Matter: Loops and Precision from Particle Physics – p. 1/3
35
Embed
SUSY Dark Matter Loops and Precision from Particle Physics
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
SUSY Dark Matter
Loops and Precision
from Particle Physics
Fawzi BOUDJEMA
LAPTH, CNRS, France
in collaboration with Andrei Semenov and David Temes
in parts with Ben Allanach, Genevieve Belanger and Sacha Pukhov
RADCOR05 F. BOUDJEMA, SUSY Dark Matter: Loops and Precision from Particle Physics – p. 1/35
Plan
Cosmology in the era of precision measurements
− Dark Matter is New Physics
− Evidence for and Precision on the matter content of the universe
− New Paradigm: can LHC/ILC match the precision of the upcoming
cosmo/astro experiments and indirectly probe the history of the early universe?
Entropy non-conservation, e.g., through decays(entropy increase will reduce the relic
abundance for example)
RADCOR05 F. BOUDJEMA, SUSY Dark Matter: Loops and Precision from Particle Physics – p. 11/35
Thermal average
must calculate all annihilation, co-annihilation processes. Each annihilation can consist of
tens of cross sections...
χ0i χ
0j → XSMYSM , χ0
1f1 → XSMYSM ,...
< σv >=
P
i,jgigj
R
(mi + mj)2ds
√sK1(
√s/T ) p2
ijσij(s)
2T`
P
igim2
i K2(mi/T )´2
,
pij is the momentum of the incoming particles in their center-of-mass frame.
pij =1
2
»
(s − (mi + mj)2)(s − (mi − mj)
2)
s
–
12
→ v
v = 0v × σv
RADCOR05 F. BOUDJEMA, SUSY Dark Matter: Loops and Precision from Particle Physics – p. 12/35
SUSY as an example
Will concentrate on supersymmetry, SUSY in particular the MSSM (no CP violation)
sometimes assumes the mSUGRA scenario (bring down the number of parameters to
4 + 1/2, but relies on RGE
RADCOR05 F. BOUDJEMA, SUSY Dark Matter: Loops and Precision from Particle Physics – p. 13/35
micrOMEGAs G. Belanger, FB, A. Semenov, A. Pukhov (CPC 02, CPC 05)
given any set of parameters it can identify LSP, NLSP, generate and calculate ω
Model defined in Lanhep (more later)
Fed into CompHEP ....tree-level, some 3000 processes could be needed.
Higgs sector: improved Higgs masses/mixings (read from FeynHiggs, for example) but
interpreted in terms of an effective scalar potential (GI), following FB and A.
Semenov (PRD 02)
Effective Lagrangian also includes important RC (Higgs couplings, ∆mb effects,..)
Interfaced with Isajet, Suspect, SoftSUSY parameters at high scale run down to the ew scale
(g − 2)µ, b → sγ, Bs → µ−µ+
NMSSM done (with C. Hugonie), JHEP 05
CP violation in progress
“open source”: procedure to define your own model, soon
RADCOR05 F. BOUDJEMA, SUSY Dark Matter: Loops and Precision from Particle Physics – p. 14/35
The mSUGRA inspired regions
100 200 300 400 500 600 700 800 900 10000
100
200
300
400
500
600
700
800
100 200 300 400 500 600 700 800 900 10000
100
200
300
400
500
600
700
800
mh = 114 GeV
m0 (
GeV
)
m1/2 (GeV)
tan β = 10 , µ < 0
τ1/χ01 co-annihilation strip
mτ1< m
χ01
WMAP
b→
sγ
(exc.)
Pre-WMAP (bulk)
Bulk region: bino LSP, lR exchange,
(small m0, M1/2)
τ1 co-annihilation: NLSP thermally
accessible, ratio of the two populations
exp(−∆M/Tf ) small m0,
M1/2 : 350 − 900GeV
Higgs Funnel: Large tan β,
χ01χ0
1 → A → bb, (τ τ),
M1/2 : 250 − 1100GeV,
m0 : 450 − 1000GeV
Focus region: small µ ∼ M1, important
higgsino component, requires very large
TeV m0
RADCOR05 F. BOUDJEMA, SUSY Dark Matter: Loops and Precision from Particle Physics – p. 15/35
Relic density around Higgs Pole with and without RC
RADCOR05 F. BOUDJEMA, SUSY Dark Matter: Loops and Precision from Particle Physics – p. 16/35
The Approach
Use micrOMEGAs+SOFTSUSY →...Softmicro..fix A0, tan β, sgn(µ) but scan on M1/2
WMAP strips imply m0 = f(M1/2): slopes
RGE also needs SM input parameters!
scale dependence of relic: default MSUSY =√
mt1mt2
: scale of EWSB conditions
theoretical uncertainty: effect of different refinements in RGE and threshold
corrections
derive accuracy within mSUGRA, relying completely on mSUGRA. accuracies on high
scale parameters and SM inputs
model independent approach: find out most relevant parameters and extract accuracy
on these (weak scale parameters)
accuracies derived in an iterative procedure and refer to the 10% WMAP precision
RADCOR05 F. BOUDJEMA, SUSY Dark Matter: Loops and Precision from Particle Physics – p. 17/35
τ1 co-annihilation region: Theory uncertainty
0.11
0.115
0.12
0.125
0.13
0.135
0.14
0.145
100 150 200 250 300 350 400
Ωh
2
Mτ~ (GeV)
Ωh2
Scale variation
0.02
0.04
0.06
0.08
0.1
0.12
0.14
0.16
0.18
100 150 200 250 300 350 400 450Ω
h2
Mτ~ (GeV)
-1 loop tau Yukawafull calculation
-1 loop neutralino-2 loop gaugino RGEs
Scale variation: From 5% (at small mτ1) to 20% at large mτ1)
2-loop gaugino RGE’s ESSENTIAL as is 1-loop threshold correction to χ01
RADCOR05 F. BOUDJEMA, SUSY Dark Matter: Loops and Precision from Particle Physics – p. 18/35
τ1 co-annihilation region: accuracy within mSUGRA
0
0.02
0.04
0.06
0.08
0.1
0.12
0.14
0.16
100 150 200 250 300 350 400 45060
80
100
120
140
160
180
a
∆A
0- (
GeV
)
mτ~ (GeV)
a(M1/2)a(m0)
∆A0-
a(tanβ)
accuracy on m0, M1/2 demanding:
3% : 1% may be achievable at LHC
for tan β require 10percent.
RADCOR05 F. BOUDJEMA, SUSY Dark Matter: Loops and Precision from Particle Physics – p. 19/35
τ1 co-annihilation region: Model Independent
0.05
0.1
0.15
0.2
0.25
0.3
0.35
0.4
100 150 200 250 300 350 400 4500.6
0.7
0.8
0.9
1
1.1
1.2
1.3
∆co
s2
θ τ- r
eq
uir
ed
, a(M
)
δM r
eq
uir
ed
(G
eV
)
m~τ (GeV)
∆cos2θτ- required
δM requireda(Mχ)
∆M must be measured to less than 1GeV
mixing angle accuracy should be feasible at ILC
accuracy on LSP mass not demanding but this is
because we have constrained ∆M .
other slepton masses need also be measured
in terms of physical parameters residual µ tan β
accuracies not demanding
Preliminary studies indicate these accuracies will
be met for the lowest mχ01
-0.1
-0.05
0
0.05
100 200 300 400 500
∆Ωh2 /Ω
h2
m~τ (GeV)
∆tanβ=10∆µ=µ
2
3
4
5
6
7
150 200 250 300 350 400 450
GeV
me~ (GeV)
δ-(∆M)
RADCOR05 F. BOUDJEMA, SUSY Dark Matter: Loops and Precision from Particle Physics – p. 20/35
direct and indirectDirect and Indirect Searches
p, e+ , γ, ν, . . .
χ01
CDMS, Edelweiss, DAMA, Genuis, ..
χ
ν
χχ → νν
Amanda, Antares, Icecube, ..
RADCOR05 F. BOUDJEMA, SUSY Dark Matter: Loops and Precision from Particle Physics – p. 21/35
Underground direct detection
W−
χ0 χ− χ0
h0, H0
(a)
χ0 χ− χ0
W− W−
q q′ q
(b)
Loops for direct detection
χm100 200 300 400 500 600 700 800 900
(pb)
S.I.
σ
-1110
-1010
-910
-810
-710
-610
-510 Edelweiss II
Zeplin II
Zeplin IV
Xenon
Genius
χm100 200 300 400 500 600 700 800 900
(pb)
S.I.
σ
-1110
-1010
-910
-810
-710
-610
-510 Edelweiss II
Zeplin II
Zeplin IV
Xenon
Genius
tan β = 10 tan β = 50
• within WMAP
RADCOR05 F. BOUDJEMA, SUSY Dark Matter: Loops and Precision from Particle Physics – p. 22/35
Annihilation into photons
dΦγ
dΩdEγ=
∑
i
dN iγ
dEγσiv
1
4πm2χ
︸ ︷︷ ︸
Physique des Particules
∫
ρ2dl︸ ︷︷ ︸
Astro
γ′s: Point to the source, independent of propagation model(s)
• continuum spectrum from χ01χ0
1 → ff , . . ., hadronisa-
tion/fragmentation (→ π0 → γ ) done through isajet/herwig
• Loop induced mono energetic photons,γγ, Zγ final states
ACT: HESS,
Magic, VERITAS,
Cangoroo, ...
Space-based:
AMS, GLAST,
Egret,...RADCOR05 F. BOUDJEMA, SUSY Dark Matter: Loops and Precision from Particle Physics – p. 23/35
FB, A. Semenov, D. Temes
S LOOP S
Need for an automatic tool for susy calculations
handles large numbers of diagrams both for tree-level
and loop level
able to compute loop diagrams at v = 0 : dark matter, LSP, move at galactic
velocities, v = 10−3
ability to check results: UV and IR finiteness but also gauge parameter independence
for example
ability to include different models easily and switch between different renormalisation
schemes
RADCOR05 F. BOUDJEMA, SUSY Dark Matter: Loops and Precision from Particle Physics – p. 24/35
Non-linear gauge implementation
LGF = − 1
ξW|(∂µ − ieαAµ − igcW βZµ)Wµ + + ξW
g
2(v + δhh + δHH + iκχ3)χ+|2
− 1
2ξZ(∂.Z + ξZ
g
2cW(v + ǫhh + ǫHH)χ3)
2 − 1
2ξγ(∂.A)2
• quite a handful of gauge parameters, but with ξi = 1, no “unphysical threshold”
• more important: no need for higher (than the minimal set)for higher rank tensors and
tedious algebraic manipulations
RADCOR05 F. BOUDJEMA, SUSY Dark Matter: Loops and Precision from Particle Physics – p. 25/35
Strategy: Exploiting and interfa ing modulesfrom dierent odes
Lagrangian of the modeldened in LanHEP- parti le ontent- intera tion terms- shifts in elds and parameters- ghost terms onstru ted by BRST↓ ↓Generi Model Classes Model-kinemati al stru tures -Feynman rules, in luding CT
⇓Evaluation viaFeynArts-FormCal LoopTools modied!!tensor redu tion inappropriate for small relative velo ities(Zero Gram determinants)
⇑Renormalisation s heme- denition of renorm. onst. in the lasses modelNon-Linear gauge-xing onstraints, gauge parameter dependen e he ks
RADCOR05 F. BOUDJEMA, SUSY Dark Matter: Loops and Precision from Particle Physics – p. 26/35
m0 = 113GeV, m1/2 = 375 GeV, A = 0, tan β = 20, µ > 0
γ lines could be distinguished from diffuse background
RADCOR05 F. BOUDJEMA, SUSY Dark Matter: Loops and Precision from Particle Physics – p. 31/35
Summary
Precise calculations for the relic density are necessary with the foreseen precision on
some cosmological parameters
Rates for direct and direct detection also need some loop calculations, here
measurements most probably will constrain the astrophysical parameters
A general code is being developed for minimal susy, but could be extended, to provide
annihilation rates for cosmo/astro and corrected cross sections for colliders
but still much to be done and to be improved
RADCOR05 F. BOUDJEMA, SUSY Dark Matter: Loops and Precision from Particle Physics – p. 32/35
Extra Expansion of Universe, Einstein Equations
Einstein Rµν − 1
2Rgµν = 8πG
„
Tµν − Λ
8πG
«
Isotropic and Homogeneous
ds2 = −dt2 + a2(t)
»
dr2
1 − kr2+ r2(dθ2 + sin2 θdφ2)
–
conservation H2 =
„
a
a
«2
=8πG
3
X
i
ρi −k
a2
→X
M
ΩM + ΩΛ + Ωk = 1 ΩM =ρM
ρcρc =
3H2
8πG
Acceleration
„
a
a
«
= −4πG
3
X
i
(ρi + 3pi) p = wρ
ρ(a) ∝ 1
a(t)3(1+w)wrad = 1/3 wM = 0 wΛ = −1
RADCOR05 F. BOUDJEMA, SUSY Dark Matter: Loops and Precision from Particle Physics – p. 33/35
Application to χ01χ0
1 → γγ, Zγ, ggHiggsino and Wino limits
• σv vs. Mχ01
(v = 0)
tan β = 10, A = 0, mA = 100 TeV, mf = 4105 TeV, M1 = 2105 TeV
M2 = 2M1 (higgsino), µ = M1 (wino)
• σv vs. % higgsino/wino (v = 0)
M2 = 2M1, µ = 50 TeV (higgsino), µ = M1, M2 = 10 TeV (wino)
δM = Mχ+1− Mχ0
1Asymptoti value for large Mχ01, σv ∼ 1/M2
WLargest ross se tion for Zγ in wino aseSmooth behaviour in higgsino ase, δM ∼ m2z/M1Constant value after transition in wino ase, δM ∼ m4
z/M31Numeri al results reprodu e analyti al behaviour
RADCOR05 F. BOUDJEMA, SUSY Dark Matter: Loops and Precision from Particle Physics – p. 34/35
Mat hing with non-perturbative omputationHisano et al.'04In the extreme higgsino and wino limits,
• one-loop treatment breaks unitarity
• non-perturbative non-relativisti approa h...Non-perturbative omputation and one-loop resultshiggsino ase (v = 0)
mf = mA = 100 TeV, M2 = 2M1 = 50 TeV, tan β = 10, A = 0, δM = 0.1 GeVResonan es an enhan e result several orders of magnitudeMat hing will take pla e around 400 − 500 GeV RADCOR05 F. BOUDJEMA, SUSY Dark Matter: Loops and Precision from Particle Physics – p. 35/35