MIMAC
MIcro-tpc MAtrix of Chambers
A Large TPC for Directional Dark Matter detection
Daniel Santos
Laboratoire de Physique Subatomique et de Cosmologie (LPSC-Grenoble) (Université Grenoble-Alpes -CNRS/IN2P3)
MIMAC MIcro-tpc MAtrix of Chambers A Large TPC for Directional Dark - - PowerPoint PPT Presentation
MIMAC MIcro-tpc MAtrix of Chambers A Large TPC for Directional Dark Matter detection Daniel Santos Laboratoire de Physique Subatomique et de Cosmologie (LPSC-Grenoble) (Universit Grenoble-Alpes -CNRS/IN2P3) MIMAC (MIcro-tpc MAtrix of
Laboratoire de Physique Subatomique et de Cosmologie (LPSC-Grenoble) (Université Grenoble-Alpes -CNRS/IN2P3)
LPSC (Grenoble) : D. Santos, F.Naraghi C.Couturier (post-doc), N. Sauzet
L.Gallin-Martel, F. Rarbi
J-F. Muraz IRFU (Saclay): P. Colas, E. Ferrer-Ribas, I. Giomataris CCPM (Marseille): J. Busto, D. Fouchez, C. Tao Tsinghua University (Beijing-China): C. Tao, I. Moric, Y. Tao XAO (Xinjiang-China): Chung-Lin Shan Neutron facility (AMANDE) : IRSN (Cadarache): T. Vinchon, B. Tampon (Ph. D.)
DM-France, Paris, Dec 1st 2016
<Vrot> ~ 220 km/s Cygnus
DM-France, Paris, Dec 1st 2016
DM-France, Paris, Dec 1st 2016
l b XG YG ZG Solar System Galactic Center VWIMP VSS Solar System’s orbit Dark matter Halo = gaz of WIMPs
Galactic coordinates
Cygnus Constellation (l = 90°,b = 0°)
After collision WIMP signal expected
DM-France, Paris, Dec 1st 2016
ZW XG YG ZG lW bW YW XW ϕR θR Recoil
108 Events with ER = [5,50] keV
Map of recoils in galactic coordinates (HealPix)
DM-France, Paris, Dec 1st 2016
DM-France, Paris, Dec 1st 2016
19F recoils (Ekin= 1-110 keV)
Angular distribution in the laboratory (with respect to the neutron direction) Produced by neutrons of 565 keV Validated experimentally at Cadarache !!
19F recoils (Ekin= 1- 40 keV)
Angular distribution in the laboratory Produced by neutrons of 200 keV Geant4 simulations ( N. Sauzet, DS. (2016)) The same kind of distributions for C !!
Wimp recoils Background 100 WIMP evts + 100 Background evts
DM-France, Paris, Dec 1st 2016
Proof of discovery: Signal pointing toward the Cygnus constellation Blind likelihood analysis in order to establish the galactic origin of the signal
Latitude galactique Longitude galactique
Signature angulaire
100 WIMP + 100 BKG
DM-France, Paris, Dec 1st 2016
Mass – cross section Galactic Halo shape Dark Matter signature
8 parameters simultaneouly constrained by only one 3D experiment
Mass Cross section σx σy σz l b
DM-France, Paris, Dec 1st 2016
Running in an Underground Laboratory
Evolution of the collected charges on the anode Scheme of a MIMAC µTPC
Measurement of the ionization energy: Charge integrator connected to the mesh coupled to a FADC sampled at 50 MHz
E~ 200 V/cm E~ 30 kV/cm
Drifting properties: V ≈ 20 µm/ns
DM-France, Paris, Dec 1st 2016
A low pressure multi-chamber detector
Strategy:
Bi-chamber module 2 x (10.8x 10.8x 25 cm3)
DM-France, Paris, Dec 1st 2016
MIMAC Target:
Magboltz Simulation
Too fast
DM-France, Paris, Dec 1st 2016
Time [ns]
1000 2000
A [ADC-Channels]
1000 2000
τ
Energy
Time [ns]
1000 2000
[ADC-Channels/ns] dA dt
5 10 peak
µ σr σl Apeak
preamplifier signal + FADC: Energy 3D - track
Dedicated fast electronics (self-triggered) Based on the MIMAC chip (64 channels)
DM-France, Paris, Dec 1st 2016
DM-France, Paris, Dec 1st 2016
19F ( 3 keV) in CF4 (50 mbar)
DM-France, Paris, Dec 1st 2016
Electrons of 7 keV In a gas detector the IQF depends strongly on the quality of the gas. The IQF needs to be measured periodically (in-situ) in a long term run experiment.
DM-France, Paris, Dec 1st 2016
0,000 ¡ 0,100 ¡ 0,200 ¡ 0,300 ¡ 0,400 ¡ 0,500 ¡ 0,600 ¡ 0,700 ¡ 0,800 ¡ 0,900 ¡ 0 ¡ 10 ¡ 20 ¡ 30 ¡ 40 ¡ 50 ¡ 60 ¡
Quenching ¡Factor ¡ Recoil ¡Energy ¡(keV) ¡
Fluorine ¡in ¡ CF4 ¡at ¡50 ¡ mbar ¡ He ¡in ¡He ¡+ ¡ 5% ¡C4H10 ¡at ¡ 350 ¡mbar ¡
DM-France, Paris, Dec 1st 2016
10 20 30 40 50 10 20 30 40 50 60 70 80 Total kinetic Energy (keV) Quenching Factor (%)
Pressure : 350 mbar 700 mbar 1000 mbar 1300 mbar
– Neutrons with a well defined energy from resonances of 7Li by a (p,n) reaction
DM-France, Paris, Dec 1st 2016
2 Eneutron cos2θ
Protons beam
7Li target
MIMAC one chamber proto
55Fe (5.9 keV) and 109Cd (3.1 keV)
sources
Neutrons F , C, H, nuclear recoils Electrons
7Li (p,n (565 keV)) nuclear reaction
22 observables built using the MIMAC readout…. and more … (Q. Riffard et al. arXiv: 1602.01738 (2016))
NR + e- Only e- NR + e- Only e-
With fast neutrons
IDM2016, 7/21/16
Focusing on the “Fluorine Endpoint”:
energies above 50 keV
rad max ~ 550 keV
IDM2016, 7/21/16
27
Method Neutron kinetic energy (mean ± 1σ) [keV] “Simple” = Joining barycenters of the extreme timeslices 542.8 ± 25.6 Fit of the centroids 541.3 ± 23.8 Fit of every (x,y,z) coincidence 545.8 ± 23.4 Principal Component Analysis 545.7 ± 23.5
/!\ Extrapolating the quenching factor above the measured range Focusing on the “Fluorine Endpoint”:
energies above 50 keV
rad max ~ 545 keV
IDM2016, 7/21/16
Theta_simple (radians)
MIMAC (bi-chamber module)at Modane Underground Laboratory (France) since June 22nd 2012. Upgraded in June 2013, and in June 2014.
(CF4+28% CHF3 + 2% C4H10)
and commanded
Many thanks to LSM staff
DM-France, Paris, Dec 1st 2016
500 600
ADC Channel
100 200 300 400
Energy [keV]
2 4 6 8 10 12 14
Cd Cr Fe Cu Cu
Date
02/07/13 31/08/13 31/10/13
a [keV/ADC-Channel]
0.005 0.01 0.015
Calibration: (once a week)
X-ray generator producing fluorescence photons from Cd, Fe, Cu foils. Threshold ~ 1 keV
Circulation system:
Excelent Gain stability in time
10 12
Eioni [keV]
2 4 6 8
Count
100 200 300 400 500 Cd Cr Fe Cu Cu
DM-France, Paris, Dec 1st 2016
DM-France, Paris, Dec 1st 2016
X-Y (anode) X-Z(t) Y-Z(t)
DM-France, Paris, Dec 1st 2016
DM-France, Paris, Dec 1st 2016
DM-France, Paris, Dec 1st 2016
222Rn chain:
222Rn
(3.8 days)
218Po
(3.1 min)
214Pb
(26.8 min)
214Bi
(19.8 min)
214Po
(0.2 ms)
210Pb
(22 years)
210Bi
(5.01 days)
210Po
(138 days)
206Pb
(stable)
Simulation (SRIM)
DM-France, Paris, Dec 1st 2016
Chamber coincidences:
3D tracks from nuclear recoil
218Po 214
DM-France, Paris, Dec 1st 2016
Nuclear recoil spectra
Mesure:{
RPR events occur at different positions in the detector…
« Anode » event
DM-France, Paris, Dec 1st 2016
Energy [keV]
5 10 15 20 25 30
MPD
1 2 3 4 5 6
10 1 10
(C. Couturier,Q. Riffard, N. Sauzet et al. in preparation ) Measurement in a MIMAC chamber of an alpha passing through the active volume parallel to the cathode at 10 cm distance.
(C. Couturier, Q. Riffard, N. Sauzet et al. 2016)
DM-France, Paris, Dec 1st 2016
Distance from grid [mm] 50 100 150 200 250 s] µ t [ ∆ 5 10 15 20 25 30 35 40
/ ndf
2
χ 7.701 / 10
retard
T 0.5266 ± 23.72
drift
v 1.139 ± 21.32 / ndf
2
χ 7.701 / 10
retard
T 0.5266 ± 23.72
drift
v 1.139 ± 21.32
DM-France, Paris, Dec 1st 2016
For more info on COMIMAC: (Muraz et al. NIM A, 2016)
X [stripe #]
120 122 124 126 128 130 132 134 136
Y [stripe #]
120 122 124 126 128 130 132 134 136
hAnode
Entries 320 Mean x 128.8 Mean y 129.2 RMS x 1.625 RMS y 1.486
2 4 6 8 10 12 hAnode
Entries 320 Mean x 128.8 Mean y 129.2 RMS x 1.625 RMS y 1.486
Anode projection
Time [timeslice #]
4 6 8 10 12 14 16 18 20 22
X [stripe #]
120 122 124 126 128 130 132 134 136
hZX
Entries 320 Mean x 14.43 Mean y 128.8 RMS x 3.057 RMS y 1.623
1 2 3 4 5 6 hZX
Entries 320 Mean x 14.43 Mean y 128.8 RMS x 3.057 RMS y 1.623
X/Time projection
Time [timeslice #]
4 6 8 10 12 14 16 18 20 22
Y [stripe #]
120 122 124 126 128 130 132 134 136
hZY
Entries 320 Mean x 14.35 Mean y 129.2 RMS x 3.014 RMS y 1.486
1 2 3 4 5 6 hZY
Entries 320 Mean x 14.35 Mean y 129.2 RMS x 3.014 RMS y 1.486
Y/Time projection
DM-France, Paris, Dec 1st 2016
# timeslices 2 4 6 8 10 12 14 16 18 20 22 24 Count 500 1000 1500 2000 2500 3000 3500 hTimeslice Entries 12227 Mean 10.04 RMS 3.962
hTimeslice
DM-France, Paris, Dec 1st 2016
Couturier et al. (in preparation)
Ionization energy (keV) 1 2 3 4 5 6 7 8 9 Angular resolution (deg) 12 14 16 18 20 22 24 26 28
Ionization energy (keV) 1 2 3 4 5 6 7 8 9 Track length (mm) 1 2 3 4 5 6
PRELIMINARY PRELIMINARY
DM-France, Paris, Dec 1st 2016
target = C (low masses), Ar, Br, Kr (high masses)
target = O (low masses), Zn, W (high masses)
target = F
X [stripe #] 120 122 124 126 128 130 132 134 136 Y [stripe #] 120 122 124 126 128 130 132 134 136 hAnode Entries 320 Mean x 128.8 Mean y 129.2 RMS x 1.625 RMS y 1.486 2 4 6 8 10 12 hAnode Entries 320 Mean x 128.8 Mean y 129.2 RMS x 1.625 RMS y 1.486Anode projection
Time [timeslice #] 4 6 8 10 12 14 16 18 20 22 X [stripe #] 120 122 124 126 128 130 132 134 136 hZX Entries 320 Mean x 14.43 Mean y 128.8 RMS x 3.057 RMS y 1.623 1 2 3 4 5 6 hZX Entries 320 Mean x 14.43 Mean y 128.8 RMS x 3.057 RMS y 1.623X/Time projection
Time [timeslice #] 4 6 8 10 12 14 16 18 20 22 Y [stripe #] 120 122 124 126 128 130 132 134 136 hZY Entries 320 Mean x 14.35 Mean y 129.2 RMS x 3.014 RMS y 1.486 1 2 3 4 5 6 hZY Entries 320 Mean x 14.35 Mean y 129.2 RMS x 3.014 RMS y 1.486Y/Time projection
No tracks ; only statistical distributions (!) D'Ambrosio et al. 2014 Capella et al. 2013
DM-France, Paris, Dec 1st 2016
(SRIM simulations)
DM-France, Paris, Dec 1st 2016
DM-France, Paris, Dec 1st 2016
DM-France, Paris, Dec 1st 2016
DM-France, Paris, Dec 1st 2016
For more information on the comparison: Couturier et al. (JCAP 2016) Initial direction of the recoil 𝜄 ↓𝑗 Direction at collision i collision i
Directionality D (preservation of the direction)
Crystal Emulsion TPC DM-France, Paris, Dec 1st 2016
DM-France, Paris, Dec 1st 2016
i) New technology anode 35cmx35cm ii) Stretched thin (12 um) grid at 512um. iii) New electronic board (1920 channels) iv) Only one big chamber New 20cmx20cm pixellized anode (1024 channels)
DM-France, Paris, Dec 1st 2016
Left: Top view of the new detector design using kapton and plexiglass instead of PCB. Right: Bottom view, showing the ASICs distribution to minimize the length of the connections.
DM-France, Paris, Dec 1st 2016
DM-France, Paris, Dec 1st 2016
30 kg.yr, 90% CL lower limits
A: 5 keV (threshold) no background 3D track with head-tail angular resolution 20o B: 20 keV background= 10evt/kg yr angular resolution 50o 3D with no head-tail
DM-France, Paris, Dec 1st 2016
DM-France, Paris, Dec 1st 2016
i) A new directional detector of nuclear recoils at low energies has been developed giving a lot of flexibility on targets, pressure, energy range… ii) Ionization quenching factor measurements have been determined experimentally and they can be checked in-situ. iii) Phenomenology studies performed by the MIMAC team show the impact of this kind of detector. iv) MIMAC bi-chamber module has been installed at Modane Underground Laboratory in June 2012. An upgraded versions in June 2013 and June 2014 and it shows an excellent gain stability. v) For the first time the 3D nuclear recoil tracks from Rn progeny have been observed. vi) New degrees of freedom are available to discriminate electrons from nuclear recoils to improve the DM search for. vii) Angular resolution and directional studies of 3D tracks are now possible with COMIMAC. viii) The 1 m3 will be the validation of a new generation of a large DM high definition detector including directionality (a needed signature for DM discovery)
DM-France, Paris, Dec 1st 2016
Discovery at 3σ With BKG (300) Without BKG à Even with a large number of background events, discovery is still possible à Only low number of WIMP events are required at low masses à A discovery (>3σ @90%CL) with BKG is possible down to 10-3-10-4 pb
Estimation of the discovery potential
MIMAC characteristics
MSSM NMSSM
DM-France, Paris, Dec 1st 2016
Up to 10-4 pb
0 WIMP, 300 bkg
Up to 10-6 pb
100 WIMP, 100 bkg
Simulated data
DM-France, Paris, Dec 1st 2016
19F : contenu en spin selon
les auteurs
DM-France, Paris, Dec 1st 2016
DRIFT MIMAC NEWAGE DMTPC Boulby Modane Kamioka SNOLAB Gas mix 73%CS2 +25%CF4 +2%O2 70%CF4 +28%CHF3 +2%C4H10 CF4 CF4 Current volume 800 L 6 L 37 L 1000 L Drift ion, 50 cm e−, 25 cm e−, 41 cm e−, 27 cm Threshold (keVee) 20 1 50 20 Readout Multi-Wire Proportional Counters Micromegas micro-pixel chamber +GEM CCD Adapted from Mayet et al. [arXiv:1602.03781]
DM-France, Paris, Dec 1st 2016
Non-baryonic matter is 6 times more important than baryonic one…
(1E0657-558) Z= 0.296 Total mass profiles Baryonic Matter
DM-France, Paris, Dec 1st 2016
Détection directe : mesure de l’énergie déposée lors de la diffusion élastique WIMP-noyau
Noyau cible
détecteur basse énergie
En tenant compte de la distribution de vitesse f(v), du facteur de forme F(q) :
nucléaire Astro SUSY ou Autres…