Prospects for dark matter detection with inelastic transitions of - - PowerPoint PPT Presentation

prospects for dark matter detection with inelastic
SMART_READER_LITE
LIVE PREVIEW

Prospects for dark matter detection with inelastic transitions of - - PowerPoint PPT Presentation

Prospects for dark matter detection with inelastic transitions of xenon Christopher M c Cabe preliminary results work in progress TeVPA, Tokyo, Japan - 27th October 2015 An old idea The original direct detection paper: Christopher


slide-1
SLIDE 1

Prospects for dark matter detection with inelastic transitions of xenon

TeVPA, Tokyo, Japan - 27th October 2015

Christopher McCabe

preliminary results —work in progress—

slide-2
SLIDE 2

An old idea…

Christopher McCabe GRAPPA - University of Amsterdam

  • The original direct detection paper:
slide-3
SLIDE 3

An old idea… Inelastic scattering

Christopher McCabe GRAPPA - University of Amsterdam

  • What is it?
  • Why is it interesting?
  • Why consider it now?

Can it ever be detected?

slide-4
SLIDE 4

What is it?

Christopher McCabe GRAPPA - University of Amsterdam

N N r e c

  • i

l DM DM N N* r e c

  • i

l DM DM N

γ

elastic scattering: inelastic scattering: measure: N’s recoil energy measure: N’s recoil energy + photon energy

slide-5
SLIDE 5

What is a good target?

Christopher McCabe GRAPPA - University of Amsterdam

XENON

slide-6
SLIDE 6

Why Xenon?

Christopher McCabe GRAPPA - University of Amsterdam

Inelastic scattering is not A2 enhanced

★ Only accessible for spin-dependent interactions ➡ Elastic and inelastic scattering rates comparable ★ Ideal target should have

  • i. good spin-dependent sensitivity
  • ii. a low lying excitation

Vietze et al arXiv:1412.6091

(. EDM−kinetic ≈ 100 keV)

slide-7
SLIDE 7

Exp

129Xe

3/2

+

11/2

  • 1/2
+

(9/2)

  • 3/2
+

5/2

+

1/2

+

(5/2)

+

(5/2)

+

7/2

+

1/2

+

1/2

+

5/2

+

5/2

+

5/2

+

7/2

+

3/2

+

3/2

+

11/2

  • 9/2
  • 11/2

100 200 300 400 500 600 700 800 900

Excitation energy (keV)

Theory Exp

131Xe

3/2

+

1/2

+

11/2

  • 9/2
  • 5/2
+

3/2

+

7/2

+

7/2

  • 3/2
+

(1/2,3/2)

+

Why Xenon?

Christopher McCabe GRAPPA - University of Amsterdam

  • 47.6% of xenon sensitive to spin-dependent interactions:

129Xe

Natural abundance: 26.4% Lowest excitation: 39.6 keV Lifetime: 0.97 ns

131Xe

Natural abundance: 21.2% Lowest excitation: 80.2 keV Lifetime: 0.48 ns

slide-8
SLIDE 8

Previous studies

Christopher McCabe GRAPPA - University of Amsterdam

  • Previous searches with single phase-detectors
  • No limits or studies for two-phase detectors (LUX, XENON)
slide-9
SLIDE 9

Why is it interesting?

Christopher McCabe GRAPPA - University of Amsterdam

Inferring properties of dark matter is difficult! We should search for all signals that provide information

  • A detection should:
  • give independent evidence for dark matter scattering
  • point strongly to a spin-dependent interaction
  • help with mass reconstruction (because of different kinematics)
slide-10
SLIDE 10

Why now?

Christopher McCabe GRAPPA - University of Amsterdam

We can accurately quantify the signal and background

  • Structure functions known (needed for cross-section)
  • Backgrounds are more-or-less known
  • Future detector properties are more-or-less known
slide-11
SLIDE 11

An old idea… Inelastic scattering

Christopher McCabe GRAPPA - University of Amsterdam

Can it ever be detected?

slide-12
SLIDE 12

Scattering rate

Christopher McCabe GRAPPA - University of Amsterdam

100 200 300 400 500 600 700 800

vmin (km/s)

10

  • 4

0.001 0.01 0.1 1

g(vmin)/g(0)

Standard Halo Model Double Power Law Tsallis Model Inelastic

131Xe

Inelastic

129Xe

Elastic

Baudis et al 1309.0825

dR dER ∝ g(vmin) = Z

vmin

d3v f(v) v

  • Rate depends on the DM velocity distribution:
  • vmin is higher for inelastic

(DM kinetic energy must also excite the nucleus)

  • This suppresses

the inelastic rate by factor ~10

slide-13
SLIDE 13

Structure functions

Christopher McCabe GRAPPA - University of Amsterdam Baudis et al 1309.0825

  • (+)
  • (+)
  • ER [keV]

SA

n (ER)

dR dER / dσ dER / Sn

A =

  • hXe∗| ¯

ψqγµγ5ψq|Xei

  • 2
  • Known for axial-vector interaction:
  • Rate depends on the structure functions
  • Smaller for inelastic

(Small ER most relevant)

  • This suppresses

the inelastic rate by factor ~10

slide-14
SLIDE 14

The rate

Christopher McCabe GRAPPA - University of Amsterdam

  • Rate as a function recoil energy (not directly measured)
  • Inelastic rate smaller by factor ~100

➡ Always see an elastic signal first

= σ

=-
  • []

/ [///]

slide-15
SLIDE 15

Two-phase xenon detectors

Christopher McCabe GRAPPA - University of Amsterdam

  • Express the signal in terms of measured quantities:

g1, g2 and drift field are the crucial parameters

E

field

Particle

e-

γ

S 1 S 2

52 phe 4540 phe

S1 = g1nγ S2 = g2ne

slide-16
SLIDE 16

Mock detectors

Christopher McCabe GRAPPA - University of Amsterdam

  • I’ll consider two benchmark scenarios:
  • Number of photons & electrons modelled with NEST

γ

XenonA200 g1=0.07 PE/ g2=12.5 PE/e

(50% extraction efficiency)

drift field=200 V/cm XenonB1000 g1=0.12 PE/ g2=50 PE/e

(100% extraction efficiency)

drift field=1000 V/cm

γ

Szydagis et al 1106.1613

slide-17
SLIDE 17
  • Include detector and recombination fluctuations
  • For same energy, electronic recoils produce a much

larger S1 and S2

Mock signals

Christopher McCabe GRAPPA - University of Amsterdam

γ

+

γ

+

⨯ ⨯

[] []

γ

+

γ

+

⨯ ⨯ ⨯ ⨯

[] []

slide-18
SLIDE 18

Mock signals

  • Looks like real data…

Christopher McCabe GRAPPA - University of Amsterdam

S1 [PE] 100 200 300 400 500 600 S2 [PE]

10000 20000 30000 40000 50000 60000

1 10

2

10

3

10

+NR ee 80 keV +NR ee 40 keV NR

Data from PandaX-I arXiv:1505.00771

😄

γ

+

γ

+

⨯ ⨯

[] []

slide-19
SLIDE 19

Background

Christopher McCabe GRAPPA - University of Amsterdam

  • Background spectra expected in LZ/XENONnT:
  • 2-neutrino — 2-beta decay of 136Xe dominates above 20 keV

νββ (±%) (±%) (±%) (±%) (±%) (±%)

  • []

/ [///]

136Xe

LZ Design: 1509.02910

129Xe 131Xe

slide-20
SLIDE 20

Reminder: Usual signal plane

Christopher McCabe GRAPPA - University of Amsterdam

LUX arXiv:1310.8214

electronic recoil band nuclear recoil band signal region S1 < 30 PE

slide-21
SLIDE 21
  • Signal region at higher values of S1
  • Large backgrounds…but some signal-to-background discrimination
  • Better discrimination for higher drift fields

Background versus signal

Christopher McCabe GRAPPA - University of Amsterdam

= σ

=-
  • []

(/)

= σ

=-
  • []

(/)

slide-22
SLIDE 22
  • Quantify the sensitivity of future experiments with a

‘discovery limit’ The smallest cross-section at which 90% of experiments can make a 3σ detection of the signal

  • Profile likelihood ratio:
  • Include background uncertainties

Discovery limit

Christopher McCabe GRAPPA - University of Amsterdam Billard et al 1110.6079

(0) = L(0

n = 0,

ˆ ˆ ~ ABG) L( ˆ

n, ˆ

~ ABG)

slide-23
SLIDE 23
  • Compare discovery limit with current/future (elastic) constraints
  • Detectable if XENON1T make discovery in next run

Discovery limit

Christopher McCabe GRAPPA - University of Amsterdam

() ()

  • mDM [GeV]

σn

0 [cm2]

() ()

  • mDM [GeV]

σn

0 [cm2]

slide-24
SLIDE 24

Summary

Christopher McCabe GRAPPA - University of Amsterdam

  • Dark matter can excite the 129Xe and 131Xe isotopes

➡ signal is nuclear recoil + photon

  • Signal is always smaller than elastic rate

➡ Can it be detected?

Yes! …need an (elastic) discovery signal in the next run of XENON1T

slide-25
SLIDE 25

Thank you

Christopher McCabe GRAPPA - University of Amsterdam

slide-26
SLIDE 26

Backup

Christopher McCabe GRAPPA - University of Amsterdam

slide-27
SLIDE 27

Christopher McCabe GRAPPA - University of Amsterdam

ne = ni − rni n = nex + rni

Gammas have shorter tracks, more recombination (r bigger) so ne smaller, ngamma bigger

slide-28
SLIDE 28

Christopher McCabe GRAPPA - University of Amsterdam

L(0

n, ~

ABG) = ⇣ µDM + P6

j=1 µBGj

⌘N N! exp @−µDM +

6

X

j=1

µBGj 1 A ·

6

Y

m=1

Lm(ABGm) ·

N

Y

i=1

" µDM µDM + P6

k=1 µBGk

fDM(S1i, log10(S2b/S1)i) +

6

X

j=1

µBGj µDM + P6

k=1 µBGk

fBGj(S1i, log10(S2b/S1)i) # ,

(0) = L(0

n = 0,

ˆ ˆ ~ ABG) L( ˆ

n, ˆ

~ ABG)

slide-29
SLIDE 29

Single-phase experiments

Christopher McCabe GRAPPA - University of Amsterdam

  • Detecting this signal could be difficult…

…impossible for single phase (S1-only)?

νββ

  • = σ
=-
  • []

/ [///]

νββ

  • = σ
=-
  • []

/ [///]

slide-30
SLIDE 30

Improvements?

Christopher McCabe GRAPPA - University of Amsterdam

  • Could have a larger exposure

➡ background dominated so only scales with the

square root

  • Could reduce backgrounds
  • Largest: 2-beta—2-neutrino decay of 136Xe

➡ Remove the 136Xe isotope

  • Try to search for displaced the S2 signal from the

recoil and photon?