Use of event-level IceCube data in SUSY scans Pat Scott Department of Physics, McGill University July 3, 202 Based on: PS, Chris Savage, Joakim Edsjö and The IceCube Collaboration (esp. Matthias Danninger and Klas Hultqvist), arxiv:207.soon Slides available from www.physics.mcgill.ca/ patscott
Outline How to find DM with neutrino telescopes How to find DM with neutrino telescopes 2
Introduction The short version:
Introduction The short version: Halo WIMPs crash into the Sun
Introduction The short version: Halo WIMPs crash into the Sun 2 Some lose enough energy in the scatter to be gravitationally bound
Introduction The short version: Halo WIMPs crash into the Sun 2 Some lose enough energy in the scatter to be gravitationally bound 3 Scatter some more, sink to the core
Introduction The short version: Halo WIMPs crash into the Sun 2 Some lose enough energy in the scatter to be gravitationally bound 3 Scatter some more, sink to the core 4 Annihilate with each other, producing neutrinos
Introduction The short version: Halo WIMPs crash into the Sun 2 Some lose enough energy in the scatter to be gravitationally bound 3 Scatter some more, sink to the core 4 Annihilate with each other, producing neutrinos 5 Propagate+oscillate their way to the Earth, convert into muons in ice/water
Introduction The short version: Halo WIMPs crash into the Sun 2 Some lose enough energy in the scatter to be gravitationally bound 3 Scatter some more, sink to the core 4 Annihilate with each other, producing neutrinos 5 Propagate+oscillate their way to the Earth, convert into muons in ice/water 6 Look for Čerenkov radiation from the muons in IceCube, ANTARES, etc
What can the muon signal tell me? Roughly: Number how much annihilation is going on in the Sun = info on σ SD, σ SI and σv Spectrum sensitive to WIMP mass m χ and branching fractions BF into different annihilation channels X Direction how likely it is that they come from the Sun In model-independent analyses a lot of this information is either discarded or not given with final limits Goal: Use as much of this information on σ SD, σ SI, σv, m χ and BF (X) as possible to directly constrain specific points and regions in WIMP model parameter spaces (+LHC+DD+... )
What can the muon signal tell me? The focus here is supersymmetry (SUSY) but this talk is on a framework, applicable to any model. Goal: Use as much of this information on σ SD, σ SI, σv, m χ and BF (X) as possible to directly constrain specific points and regions in WIMP model parameter spaces (+LHC+DD+... )
What can the muon signal tell me? The focus here is supersymmetry (SUSY) but this talk is on a framework, applicable to any model. All the methods discussed here are available in DarkSUSY 5.: www.darksusy.org Goal: Use as much of this information on σ SD, σ SI, σv, m χ and BF (X) as possible to directly constrain specific points and regions in WIMP model parameter spaces (+LHC+DD+... )
What can the muon signal tell me? The focus here is supersymmetry (SUSY) but this talk is on a framework, applicable to any model. All the methods discussed here are available in DarkSUSY 5.: www.darksusy.org All IceCube data used are available at http://icecube.wisc.edu/science/data/ic22-solar-wimp (and in DarkSUSY, for convenience) Goal: Use as much of this information on σ SD, σ SI, σv, m χ and BF (X) as possible to directly constrain specific points and regions in WIMP model parameter spaces (+LHC+DD+... )
Outline How to find DM with neutrino telescopes How to find DM with neutrino telescopes 2
SUSY Scanning with IceCube Simple Likelihood Simplest way to do anything is to make it a counting problem... Compare observed number of events n and predicted number θ for each model, taking into account error σ ɛ on acceptance: L num(n θ BG + θ sig ) = 2πσɛ 0 (θ BG + ɛθ sig ) n e (θ BG+ɛθ sig ) n! [ ɛ exp 2 ( ) ] ln ɛ 2 dɛ. σ ɛ () Nuisance parameter ɛ takes into account systematic errors on effective area, from theory, etc. σ ɛ 20% for IceCube.
SUSY Scanning with IceCube Full Likelihood Full unbinned likelihood with number (L num ), spectral (L spec ) and angular (L ang ) parts with L spec,i (N i, Ξ) = and L ang,i ( cos φ i ) = L = L num(n ) n L spec,i L ang,i (2) i= θbg dp BG dn i (N i )+ θsignal θbg 0 E disp(n i E i ) dpsignal de i (E i, Ξ) de i (3) dp BG ( cos φ i ) + θsignal PSF( cos φ i ) (4) d cos φ i
SUSY Scanning with IceCube Full Likelihood Full unbinned likelihood with number (L num ), spectral (L spec ) and angular (L ang ) parts L = L num(n ) with Number of lit channels (energy estimator) L spec,i (N i, Ξ) = θbg dp BG (N i )+ θsignal E disp(n i E dn i and L ang,i ( cos φ i ) = θbg n L spec,i L ang,i (2) i= 0 i ) dpsignal de i (E i, Ξ) de i (3) dp BG ( cos φ i ) + θsignal PSF( cos φ i ) (4) d cos φ i
SUSY Scanning with IceCube Full Likelihood Full unbinned likelihood with number (L num ), spectral (L spec ) and angular (L ang ) parts L = L num(n ) with Number of lit channels (energy estimator) L spec,i (N i, Ξ) = θbg dp BG (N i )+ θsignal E disp(n i E dn i and L ang,i ( cos φ i ) = θbg n L spec,i L ang,i (2) i= 0 i ) dpsignal de SUSY parameters i (E i, Ξ) de i (3) dp BG ( cos φ i ) + θsignal PSF( cos φ i ) (4) d cos φ i
SUSY Scanning with IceCube Full Likelihood Full unbinned likelihood with number (L num ), spectral (L spec ) and angular (L ang ) parts with L spec,i (N i, Ξ) = and L ang,i ( cos φ i ) = L = L num(n ) n L spec,i L ang,i (2) i= θbg dp BG dn i (N i )+ θsignal θbg Predicted signal spectrum (from theory) 0 E disp(n i E i ) dpsignal de i (E i, Ξ) de i (3) dp BG ( cos φ i ) + θsignal PSF( cos φ i ) (4) d cos φ i
SUSY Scanning with IceCube Full Likelihood Full unbinned likelihood with number (L num ), spectral (L spec ) and angular (L ang ) parts with L spec,i (N i, Ξ) = and L ang,i ( cos φ i ) = L = L num(n ) n L spec,i L ang,i (2) i= θbg dp BG dn i (N i )+ θsignal θbg Predicted signal spectrum (from theory) E disp(n i E (E de i i, Ξ) de i (3) Instrument response function 0 i ) dpsignal dp BG ( cos φ i ) + θsignal PSF( cos φ i ) (4) d cos φ i
SUSY Scanning with IceCube Full Likelihood Full unbinned likelihood with number (L num ), spectral (L spec ) and angular (L ang ) parts with L spec,i (N i, Ξ) = and L ang,i ( cos φ i ) = L = L num(n ) n L spec,i L ang,i (2) i= θbg dp BG dn i (N i )+ θsignal Observed BG distribution θbg Predicted signal spectrum (from theory) 0 E disp(n i E (E de i i, Ξ) de i (3) Instrument response function i ) dpsignal dp BG ( cos φ i ) + θsignal PSF( cos φ i ) (4) d cos φ i
SUSY Scanning with IceCube Full Likelihood Full unbinned likelihood with number (L num ), spectral (L spec ) and angular (L ang ) parts with L spec,i (N i, Ξ) = and L ang,i ( cos φ i ) = L = L num(n ) n L spec,i L ang,i (2) i= θbg dp BG dn i (N i )+ θsignal θbg Event arrival angle 0 E disp(n i E i ) dpsignal de i (E i, Ξ) de i (3) dp BG ( cos φ i ) + θsignal PSF( cos φ i ) (4) d cos φ i
SUSY Scanning with IceCube Full Likelihood Full unbinned likelihood with number (L num ), spectral (L spec ) and angular (L ang ) parts with L spec,i (N i, Ξ) = and L ang,i ( cos φ i ) = L = L num(n ) n L spec,i L ang,i (2) i= θbg dp BG dn i (N i )+ θsignal θbg 0 E disp(n i E i ) dpsignal de i (E i, Ξ) de i (3) dp BG ( cos φ i ) + θsignal PSF( cos φ i ) (4) d cos φ i Predicted signal direction (δ function at Sun)
SUSY Scanning with IceCube Full Likelihood Full unbinned likelihood with number (L num ), spectral (L spec ) and angular (L ang ) parts with L spec,i (N i, Ξ) = and L ang,i ( cos φ i ) = L = L num(n ) n L spec,i L ang,i (2) i= θbg dp BG dn i (N i )+ θsignal θbg 0 E disp(n i E (E de i i, Ξ) de i (3) Instrument response function i ) dpsignal dp BG ( cos φ i ) + θsignal PSF( cos φ i ) (4) d cos φ i Predicted signal direction (δ function at Sun)
SUSY Scanning with IceCube Full Likelihood Full unbinned likelihood with number (L num ), spectral (L spec ) and angular (L ang ) parts with L spec,i (N i, Ξ) = and L ang,i ( cos φ i ) = L = L num(n ) n L spec,i L ang,i (2) i= θbg dp BG dn i (N i )+ θsignal Observed BG distribution θbg 0 E disp(n i E (E de i i, Ξ) de i (3) Instrument response function i ) dpsignal dp BG ( cos φ i ) + θsignal PSF( cos φ i ) (4) d cos φ i Predicted signal direction (δ function at Sun)
SUSY Scanning with IceCube Statistics 0 Why simple IN/OUT analyses are not enough:... Only partial goodness of fit, no measure of convergence, no idea how to generalise to regions or whole space. Frequency/density of models in IN/OUT scans means essentially nothing. More information comes from a global statistical fit. = parameter estimation exercise Composite likelihood made up of observations from all over: dark matter relic density from WMAP precision electroweak tests at LEP LEP limits on sparticle masses B-factory data (rare decays, b sγ) muon anomalous magnetic moment LHC searches, direct detection (only roughly implemented for now)
Example: SUSY Scanning with IceCube Global Fits CMSSM, IceCube-22 events m 0 m /2 and m χ 0 nuclear scattering cross-sections Scott, Danninger, Savage, Edsjö & Hultqvist 202 42 Scott, Danninger, Savage, Edsjö & Hultqvist 202 38 Scott, Danninger, Savage, Edsjö & Hultqvist 202 IC22 global fit IC22 global fit m0 (TeV) 3 2 Mean IC22 global fit log 0 ( σsi,p/cm 2 ) 43 44 45 46 log 0 ( σsd,p/cm 2 ) 39 40 4 42 43 Mean Mean 0.0 0.5.5 2.0 2.5 m (TeV) 2 0.0 m χ 0 (TeV) 0.0 m χ 0 (TeV) Contours indicate σ and 2σ credible regions Grey contours correspond to fit without IceCube data Shading+contours indicate relative probability only, not overall goodness of fit
Example: SUSY Scanning with IceCube Global Fits Base Observables Scott, Danninger, Savage, Edsjö, Hultqvist & The IceCube Collab. (202) 42 Scott, Danninger, Savage, Edsjö, Hultqvist & The IceCube Collab. (202) 38 Scott, Danninger, Savage, Edsjö, Hultqvist & The IceCube Collab. (202) Old data Old data m0 (TeV) 3 2 Mean Old data log 0 ( σsi,p/cm 2 ) 43 44 45 46 log 0 ( σsd,p/cm 2 ) 39 40 4 42 43 Mean Mean 0.0 0.5.5 2.0 2.5 m (TeV) 2 0.0 m χ 0 (TeV) 0.0 m χ 0 (TeV)
Example: SUSY Scanning with IceCube Global Fits Base Observables + XENON-00 Grey contours correspond to Base Observables only Scott, Danninger, Savage, Edsjö, Hultqvist & The IceCube Collab. (202) 42 Scott, Danninger, Savage, Edsjö, Hultqvist & The IceCube Collab. (202) 38 Scott, Danninger, Savage, Edsjö, Hultqvist & The IceCube Collab. (202) XENON00 XENON00 m0 (TeV) 3 2 Mean XENON00 log 0 ( σsi,p/cm 2 ) 43 44 45 46 log 0 ( σsd,p/cm 2 ) 39 40 4 42 43 Mean Mean 0.0 0.5.5 2.0 2.5 m (TeV) 2 0.0 m χ 0 (TeV) 0.0 m χ 0 (TeV)
Example: SUSY Scanning with IceCube Global Fits Base Observables + XENON-00 + CMS 5 fb Grey contours correspond to Base Observables only m0 (TeV) 3 2 Scott, Danninger, Savage, Edsjö, Hultqvist & The IceCube Collab. (202) Mean XENON00 +CMS 5 fb log 0 ( σsi,p/cm 2 ) 42 43 44 45 46 Scott, Danninger, Savage, Edsjö, Hultqvist & The IceCube Collab. (202) XENON00 +CMS 5 fb Mean log 0 ( σsd,p/cm 2 ) 38 39 40 4 42 43 Scott, Danninger, Savage, Edsjö, Hultqvist & The IceCube Collab. (202) XENON00 +CMS 5 fb Mean 0.0 0.5.5 2.0 2.5 m (TeV) 2 0.0 m χ 0 (TeV) 0.0 m χ 0 (TeV)
Example: SUSY Scanning with IceCube Global Fits Base Observables + XENON-00 + CMS 5 fb + IC22 00 Grey contours correspond to Base Observables only m0 (TeV) 3 2 Scott, Danninger, Savage, Edsjö, Hultqvist & The IceCube Collab. (202) Mean IC22 00 +XENON00 +CMS 5 fb log 0 ( σsi,p/cm 2 ) 42 43 44 45 46 Scott, Danninger, Savage, Edsjö, Hultqvist & The IceCube Collab. (202) IC22 00 +XENON00 +CMS 5 fb Mean log 0 ( σsd,p/cm 2 ) 38 39 40 4 42 43 Scott, Danninger, Savage, Edsjö, Hultqvist & The IceCube Collab. (202) Mean IC22 00 +XENON00 +CMS 5 fb 0.0 0.5.5 2.0 2.5 m (TeV) 2 0.0 m χ 0 (TeV) 0.0 m χ 0 (TeV) CMSSM, IceCube-22 with 00 boosted effective area (kinda like IceCube-86+DeepCore)
Example: Model Recovery m0 (TeV) 3 2 Scott, Danninger, Savage, Edsjö & Hultqvist 202 Mean IC22 00 Predicted signal events 40 20 00 80 60 40 20 IC22 00 Scott, Danninger, Savage, Edsjö & Hultqvist 202 Mean True value CMSSM, IceCube-22 00 signal reconstruction 60 signal events, 500 GeV, χχ W + W 0.0 0.5.5 2.0 m (TeV) 2 0 200 400 600 800 m χ 0 (GeV) 42 Scott, Danninger, Savage, Edsjö & Hultqvist 202 38 Scott, Danninger, Savage, Edsjö & Hultqvist 202 IC22 00 IC22 00 log 0 ( σsi,p/cm 2 ) 43 44 45 Mean log 0 ( σsd,p/cm 2 ) 39 40 4 42 43 Mean Grey contours correspond to reconstruction without energy information 0 200 400 600 800 m χ 0 (GeV) 0 200 400 600 800 m χ 0 (GeV)
Closing remarks A framework for directly comparing event-level neutrino telescope data to individual points in theory parameter spaces is in place These iclike extensions are available in DarkSUSY 5. IceCube event data has been released in a form digestible by DarkSUSY 5. Direct SUSY analyses of IC79 data are on the way Many models exist that only IC86 will be sensitive to The extensions can be used equally well for non-susy BSM scenarios too
Outline Backup Slides 3 Backup Slides
Advertisement Backup Slides If you liked any of the plots in this talk... Flat priors pippi v Best fit Mean True value 2.5 3.0 ) 3.5 log 0 (m /TeV 2 log 0 (m0/tev) 3.5 3.0 2.5 2.0 Best fit Mean True value pippi v Flat priors 2.5 3.0 ) 3.5 log 0 (m /TeV 2 A0 (TeV) 4 2 0 2 Flat priors Prof. likelihood Best fit Mean pippi v 0 20 30 40 50 60 tan β Profile likelihood ratio Λ = L/Lmax Pippi parse it, plot it arxiv:206.2245 http://github.com/patscott/pippi Generic pdflatex sample parser, post-processor & plotter
CMS 5 fb analyses Backup Slides [GeV] m /2 800 700 600 500 400 = LSP τ CMS Preliminary m(q ~ ) = 500 m(q ~ ) = 000 SS Dilepton Razor L int m(q ~ ) = 2000 - = 4.98 fb, m(g ~ ) = 500 MT2 m( ~ q) = 2500 s = 7 TeV tan(β)=0 A 0 = 0 GeV µ > 0 m t = 73.2 GeV ~ ± LEP2 l ± LEP2 χ m(g ~ ) = 000 300 200 00 m(g ~ ) = 500 OS Dilepton Multi-Lepton Lepton Jets+MHT Non-Convergent RGE's 500 000 500 2000 2500 3000 [GeV] m 0 No EWSB
Backup Slides XENON-00 00-day analysis WIMP-Nucleon Cross Section [cm 2 ] 0 0 0 0 0 0-39 -40-4 -42-43 -44 DAMA/Na CoGeNT DAMA/I CDMS (20) CDMS (200) XENON0 (S2 only, 20) EDELWEISS (20) XENON00 (20) observed limit (90% CL) Expected limit of this run: ± σ expected ± 2 σ expected XENON00 (200) Trotta et al. 0-45 Buchmueller et al. 6 7 8 90 20 30 40 50 00 200 300 400 000 2 WIMP Mass [GeV/c ]