Search for ultra-high energy photons and neutrinos at the Pierre Auger Observatory Mathieu Tartare 1 on behalf of the Pierre Auger Collaboration 2 1 Laboratoire de Physique Subatomique et de Cosmologie (LPSC) - Grenoble, France UMR UJF/INPG/CNRS 2 Observatorio Pierre Auger, Av. San Martín Norte 304, (5613) Malargüe, Argentina Full author list: http://www.auger.org/archive/authors 2013 02.html
Context The UHECR puzzle Composition? Origin? GZK cutoff or maximum acceleration? If GZK cutoff : photons and neutrinos expected Constraints on production models (top-down models) Neutrinos point directly to sources Limits on EeV photon fraction reduce systematic in other analyses The Pierre Auger Observatory Surface Detector (SD) 1600 Cherenkov water tanks with a 1.5 km spacing covering 3000 km 2 61 tanks in 750 m grid ( infill low energies) 0 % duty cycle Fluorescence Detector (FD) 4 sites + 1 HEAT (low energies) 6 fluorescence telescopes per site Field of view : 30 30 14 % duty cycle Mathieu Tartare [tartare@lpsc.in2p3.fr] Rencontres de Moriond, 03/12/2013 1
Ultra-high energy neutrinos Mathieu Tartare [tartare@lpsc.in2p3.fr] Rencontres de Moriond, 03/12/2013 2
Identifying neutrinos in data Identify neutrino showers from nucleonic showers background Regular nucleonic showers : Interact high in the atmosphere At ground : mainly muons, flat shower front Neutrinos May interact closer to the ground level At ground : EM component, curved shower front Discrimination power enhanced at high zenith-angle Selection criteria : Inclined young (=deep) shower Mathieu Tartare [tartare@lpsc.in2p3.fr] Rencontres de Moriond, 03/12/2013 3
Identification criteria : Inclined Showers Elongated footprint Apparent speed along major axis L Vertical Shower T ij 0 V c Horizontal Shower V c Earth Skimming Elongated event : L/W 5 Down Going Elongated event : L/W 3 Apparent speed : 0.29 m ns V 0.31 m ns SD(V ) < 0.08 m.ns 1 Apparent speed : V 0.313 m ns SD(V ) V < 8% Rec. zenith angle θ 75 Mathieu Tartare [tartare@lpsc.in2p3.fr] Rencontres de Moriond, 03/12/2013 4
Identification criteria : Inclined Showers Elongated footprint Apparent speed along major axis L Vertical Shower T ij 0 V c Horizontal Shower V c Events 3 2 Events 3 2 Events 3 2 1 1 1 0 5 15 20 25 30 L/ W 0.2 0.3 0.4 0.5 0.6 <V> (m ns -1 ) 0 0.1 0.2 0.3 v (m ns-1 ) Mathieu Tartare [tartare@lpsc.in2p3.fr] Rencontres de Moriond, 03/12/2013 4
Identification criteria : Young Showers Young/Deep Showers = Broad Signal Signal extented in time Time over Threshold trigger Large Area over Peak value Station trigger Area over Peak ToT station ToT fraction > 0.6 Thr station ToT station ToT fraction < 0.6 Thr station Earth Skimming ToT fraction Down Going Fisher discriminant based on AoP Mathieu Tartare [tartare@lpsc.in2p3.fr] Rencontres de Moriond, 03/12/2013 5
Identification criteria : Young Showers Young/Deep Showers = Broad Signal Signal extented in time Time over Threshold trigger Large Area over Peak value Events 3 2 1 Training Data simulation 0 0.2 0.4 0.6 0.8 1 ToT stations Events 3 2 1-1 - -5 0 5 Fisher discriminant value Mathieu Tartare [tartare@lpsc.in2p3.fr] Rencontres de Moriond, 03/12/2013 5
Exposure Based on detailed dedicated Monte Carlo simulations 17 16 Earth-skimming (3.5 yr of full Auger) s sr] 2 Exposure [cm 15 14 13 Down-going (2 yr of full Auger) Total CC e CCµ CCτ NC x CCτ Mountains 17 18 ν energy (ev) 19 20 Mathieu Tartare [tartare@lpsc.in2p3.fr] Rencontres de Moriond, 03/12/2013 6
Limits to diffuse fluxes After unblinding : 0 candidates survive the cuts Assuming a flux φ(e) = k E 2 k < 3.2 8 GeV cm 2 s 1 sr 1 in 1.6 17 ev< E < 2 19 ev k < 1.7 7 GeV cm 2 s 1 sr 1 in 1 17 ev< E < 1 20 ev -5 Single flavour neutrino limits (90% CL) s -1 sr -1 ] -2 ) [ GeV cm ν Φ(E k = E ν 2-6 -7-8 -9 - -11 ν limits Auger downward-going Auger Earth-skimming IceCube-40 (333.5 days) Anita-II (28.5 days) 17 18 19 Cosmogenic models Ahlers 20 Kotera 20 20 E ν [ev] Mathieu Tartare [tartare@lpsc.in2p3.fr] Rencontres de Moriond, 03/12/2013 7 21 k = N up Emax φ(e)ede E min Search sample : Earth skimming : Jan. 04 - May Down going : Nov. 07 - May
Limits to point sources Search for point-like sources of UHEν over a broad range of declinations : [north -65, south +55 ] Single flavour neutrino limits (90% CL) s -1 ] -5 Auger downward-going Auger Earth-skimming -2 ) [GeV cm -6 F(E ν = E k PS 2 ν -7-80 -60-40 -20 0 20 40 60 80 Source declinationδ [deg] Mathieu Tartare [tartare@lpsc.in2p3.fr] Rencontres de Moriond, 03/12/2013 8
Ultra-high energy photons Mathieu Tartare [tartare@lpsc.in2p3.fr] Rencontres de Moriond, 03/12/2013 9
Identifying UHE photons Fluorescence Detector Deeper development of the air showers. larger X max Surface Detector Smaller signal at a given distance & same energy Fewer triggered stations S b = ( ) 4 Ri S i 00 i smaller S b Mathieu Tartare [tartare@lpsc.in2p3.fr] Rencontres de Moriond, 03/12/2013
Identifying UHE photons FD : Deeper X max SD : Smaller S b ) -2 (g cm X max 1200 00 Monte Carlo Simulations photon 18 < log (E /ev) < 18.5 Photon-like events 800 600 proton -3-2 -1 0 1 2 3 log (S ) b Mathieu Tartare [tartare@lpsc.in2p3.fr] Rencontres de Moriond, 03/12/2013 11
Photon selection Data selection Reconstruction Level : Good geometry and longitudinal profile Zenith angle < 60 X max within the field of view Quality cuts : Time periods with clouds rejected Active station within 2 km from the shower axis 4 Fisher analysis 3 separate energy intervals (1-3 EeV, 3- EeV, > EeV) Background : proton QGSJET-II-03 Event are tagged as photon candidates for X > X cut Proton background 1% Photon selection efficiency = 50% Mathieu Tartare [tartare@lpsc.in2p3.fr] Rencontres de Moriond, 03/12/2013 12
Photon candidates Using hybrid data from Jan. 2005 to Sep. 20 : 6, 0, 0, 0 and 0 candidates above 1, 2, 3, 5 and EeV Compatible with the expected nuclear background Example of a selected candidate )] -2 de/dx [PeV/(g cm 4 E [EeV] = 1.18 ± 0.09 Xmax [gcm 3.5-2 ] = 23 ± log(sb) = 1.16 ± 0.27 3 2.5 2 1.5 1 0.5 2 /Ndf= 93.8/95 700 800 900 00 10 1200 1300-2 slant depth [g cm ] ) -2 (g cm X max 1200 10 00 900 800 700 proton simulations ~ 2% of protons are marked as candidates photon-like events selected candidate 600-1 -0.5 0 0.5 1 1.5 2 2.5 log (S Dedicated proton simulations : same energy, arrival direction, core position and detector configuration of this selected candidate ) b Mathieu Tartare [tartare@lpsc.in2p3.fr] Rencontres de Moriond, 03/12/2013 13
Hybrid exposure Upper Limit to the Integral Photon Flux : Φ 95CL γ = N 95CL γ (E γ > E 0 ) E γ,min sr yr] 3 Nγ 95CL : Number of photon candidates (at 95% C.L.) with energy E γ above the threshold E 0. Nuclear background not subtracted (conservative approach) E γ,min : Minimum hybrid exposure for photons. Time dependent simulations 2 Hybrid Exposure for photons [km 2 1-1 Photon candidate level Realistic and time dependent simulations J anuary January 2005 2005 - September - 20 level Photon candidate 17.5 18 18.5 19 19.5 log (Energy/eV) (method according to Astrop. Phys. 34, 2011) Mathieu Tartare [tartare@lpsc.in2p3.fr] Rencontres de Moriond, 03/12/2013 14
Limits on photon flux sr -1 y -1 ] -2 [km 0 Integral Flux E>E 1-1 -2-3 Y Auger Hyb 2011 Auger Hyb 2009 Y upper limits 95% C.L. A TA 2011 Auger SD A SHDM SHDM TD Z-burst GZK Systematic uncertainties Exposure, X max, S b, Energy scale, hadronic interaction model and mass composition assumptions : +20% 64% (E 0 = 1 EeV) +15% 36% (E 0 > 1 EeV) 18 19 20 Energy[eV] Upper limit on the integral photon fraction assuming the Auger spectrum : 0.4%, 0.5%, 1.0%, 2.6% and 8.9% at E > 1, 2, 3, 5 and EeV Mathieu Tartare [tartare@lpsc.in2p3.fr] Rencontres de Moriond, 03/12/2013 15
Conclusions Neutrino search using surface detector Earth skimming (sensitive to ν τ ) : k < 3.2 8 GeV cm 2 s 1 sr 1 in 1.6 17 ev< E < 2 19 ev Down going (sensitive to all flavors) : k < 1.7 7 GeV cm 2 s 1 sr 1 in 1 17 ev< E < 1 20 ev The surface detector of the Pierre Auger Observatory is sensitive to potential point sources of UHE neutrinos in a broad range in declination. Photon search using hybrid data 6, 0, 0, 0 and 0 candidates above 1, 2, 3, 5 and EeV Upper limit to the integral photon flux : 8.2 2 km 2 sr 1 year 1 above 1 EeV 2.0 2 km 2 sr 1 year 1 above 2, 3, 5 and EeV Corresponding to photon fractions limit of : 0.4%, 0.5%, 1.0%, 2.6% and 8.9% at E > 1, 2, 3, 5 and EeV Mathieu Tartare [tartare@lpsc.in2p3.fr] Rencontres de Moriond, 03/12/2013 16
Perspective & Outlook GZK within reach in the next few years Photon & neutrino limits provide severe constraints on top-down models. Photon limits allow to reduce systematic uncertainties on mass composition, energy spectrum & cross section measurements Earth skimming and Downgoing combination in progress. Directional search for UHE photons in progress Search for UHE photons and neutrinos Mathieu Tartare [tartare@lpsc.in2p3.fr] Rencontres de Moriond, 03/12/2013 17
Perspective & Outlook GZK within reach in the next few years Photon & neutrino limits provide severe constraints on top-down models. Photon limits allow to reduce systematic uncertainties on mass composition, energy spectrum & cross section measurements Earth skimming and Downgoing combination in progress. Directional search for UHE photons in progress Search for UHE photons and neutrinos a postdoctoral position Mathieu Tartare [tartare@lpsc.in2p3.fr] Rencontres de Moriond, 03/12/2013 17
Backup slides Mathieu Tartare [tartare@lpsc.in2p3.fr] Rencontres de Moriond, 03/12/2013 18
General neutrino search strategy Mathieu Tartare [tartare@lpsc.in2p3.fr] Rencontres de Moriond, 03/12/2013 19
Limits to point-like sources s -1 ] -5-6 Single flavour neutrino limits (90% CL) Down-going (2yr) Earth-skimming (3.5yr) Fixed energy range 17 ev < E < 20 ev 1.6 17 ev < E < 2 19 ev g(e) [GeV cm -2 2 E -7-8 Cen A IceCube-40 (375.5 days) Energy range depends on declination Below 17 ev -9-80 -60-40 -20 0 20 40 60 80 Source declination [deg] Mathieu Tartare [tartare@lpsc.in2p3.fr] Rencontres de Moriond, 03/12/2013 20
Limits to CenA -4 Centaurus A - Single flavour neutrino limits (90% CL) s -1 ] -2-5 Auger Downward-going LUNASKA 2008 ) [GeV cm ν F(E -6-7 IceCube 2011b Auger Earth-skimming Cuoco 2008 = E 2 ν -8 k PS -9-15 16 17 18 E ν Kachelriess 2009 [ev] Mathieu Tartare [tartare@lpsc.in2p3.fr] Rencontres de Moriond, 03/12/2013 21 19 20 21 22 23
Point source as seen by Auger 1 Source atδ= -43 o 1 Source atδ= -80 o 0.5 0.5 θ = 75 o θ = 75 o cosθ 0 θ = 90 o cosθ 0 θ = 90 o -0.5-0.5-1 0 0.2 0.4 0.6 0.8 1 Time t [sidereal days] -1 0 0.2 0.4 0.6 0.8 1 Time [sidereal days] Mathieu Tartare [tartare@lpsc.in2p3.fr] Rencontres de Moriond, 03/12/2013 22
safe cut on Fisher value such that expected background < 1 event / 20 yearsof Auger data 11 Mathieu Tartare [tartare@lpsc.in2p3.fr] Rencontres de Moriond, 03/12/2013 23
Down- going: estimation of background Assume an exponential shape for the tail of background distribution of F extrapolation to find the value of F_cut corresponding to 1 background event in a given time Mathieu Tartare [tartare@lpsc.in2p3.fr] Rencontres de Moriond, 03/12/2013 24
Expected number of events Model& reference Earth-skimming Downward-going Cosmogenic(Fermi)[40] 0.6 0.1 BBR(AGNs)[41] 5.1 0.8 Exotic(SH relics)[42] 3.0 1.0 Table 4: Number of expected events for several theoretical models of UHE neutrino production, given the exposure of the surface detector of the Pierre Auger Observatory to Earth-skimming and downward-going neutrinos(table3). Integrated limit The limits indicate the level of the diffuse neutrino flux needed to detect Nup events with a Poisson probability of 90% given the exposure accumulated during the 3.5 years for Earth-skimming (2.0 years for downward-going) of equivalent time of a full SD Mathieu Tartare [tartare@lpsc.in2p3.fr] Rencontres de Moriond, 03/12/2013 25
Elongation rate 1200 <X max > (g cm -2 ) 10 00 900 800 700 Fly s Eye HiRes-MIA HiRes 2004 Yakutsk 2001 Yakutsk 2005 CASA-BLANCA HEGRA-AIROBICC SPASE-VULCAN DICE TUNKA photon proton photon with preshower 600 500 400 iron QGSJET 01 QGSJET II SIBYLL 2.1 14 15 16 17 18 19 20 21 E lab (ev) Mathieu Tartare [tartare@lpsc.in2p3.fr] Rencontres de Moriond, 03/12/2013 26
Systematics Mathieu Tartare [tartare@lpsc.in2p3.fr] Rencontres de Moriond, 03/12/2013 27
Photon candidates Hybrid data J an 2005 - Sep 20 (~ a factor 2 more events compared to previous analysis) 6, 0, 0, 0 and 0 candidates above 1, 2, 3, 5 and EeV compatible with the expected nuclear background E (EeV) Xmax (gcm -2 ) log(sb) 1.18 ± 0.09 23 ± 1.16 ± 0.27 1.59 ± 0.19 981 ± 16 1.35 ± 0.09 1.20 ± 0.09 952 ± 17 1.17 ± 0.07 1.25 ± 0.16 958 ± 27 1.08 ± 0.24 1.13 ± 0.17 1119 ± 29 1.56 ± 0.21 1.42 ± 0.19 959 ± 22 1.08 ± 0.55 ) -2 Xmax (g cm 1200 1150 10 50 00 950 Photon candidates Jan 2005 - Sept 20 1.18 EeV 1.59 EeV 1.25 EeV 1.20 EeV 1.42 EeV 1.13 EeV 900 850 0.8 1 1.2 1.4 1.6 1.8 log (S ) b Mathieu Tartare [tartare@lpsc.in2p3.fr] Rencontres de Moriond, 03/12/2013 28