Arxiv:1303.1406V1 [Astro-Ph.HE] 6 Mar 2013 the flux of Gamma-Rays from Dark Matter Annihilation Takes a Surprisingly Simple Form 2
Total Page:16
File Type:pdf, Size:1020Kb
4th Fermi Symposium : Monterey, CA : 28 Oct-2 Nov 2012 1 The VERITAS Dark Matter Program Alex Geringer-Sameth∗ for the VERITAS Collaboration Department of Physics, Brown University, 182 Hope St., Providence, RI 02912 The VERITAS array of Cherenkov telescopes, designed for the detection of gamma-rays in the 100 GeV-10 TeV energy range, performs dark matter searches over a wide variety of targets. VERITAS continues to carry out focused observations of dwarf spheroidal galaxies in the Local Group, of the Milky Way galactic center, and of Fermi-LAT unidentified sources. This report presents our extensive observations of these targets, new statistical techniques, and current constraints on dark matter particle physics derived from these observations. 1. Introduction Earth's atmosphere as a target for high-energy cosmic particles. An incoming gamma-ray may interact in the The characterization of dark matter beyond its Earth's atmosphere, initiating a shower of secondary gravitational interactions is currently a central task particles that travel at speeds greater than the local of modern particle physics. A generic and well- (in air) speed of light. This entails the emission of motivated dark matter candidate is a weakly in- ultraviolet Cherenkov radiation. The four telescopes teracting massive particle (WIMP). Such particles of the VERITAS array capture images of the shower have masses in the GeV-TeV range and may inter- using this Cherenkov light. The images are analyzed act with the Standard Model through the weak force. to reconstruct the direction of the original particle as Searches for WIMPs are performed at particle accel- well as its energy. erators [Fairbairn et al. 2007], where dark matter may VERITAS is sensitive to showers initiated by be produced in high-energy collisions, and by low- gamma-rays with energies from around 100 GeV to background direct detection experiments which look more than 30 TeV. Because of the steep energy spec- for the scattering of dark matter particles off nuclei tra of typical sources, fluxes in the VERITAS en- [Gaitskell 2004]. In astrophysics, the problem is be- ergy band are much smaller than those seen with ing addressed through the field of indirect detection Fermi. However, the difficulty of measuring such small [Jungman et al. 1996, Bergstr¨om2000]. Here the goal fluxes is mitigated by the enormous effective areas of is the detection of the end-products of dark matter ACTs. VERITAS is sensitive to gamma-ray show- annihilation which can take place throughout the Uni- ers incident over an (energy-dependent) area of ap- verse. proximately 109 cm2 (compare with ∼ 104 cm2 for the Dark matter annihilation into Standard Model par- LAT). At 1 TeV the energy resolution of VERITAS is ticles generically gives rise to high-energy gamma- approximately 15%. The array has an angular resolu- rays. There are several classes of experiments cur- tion of about 0:14◦ at 200 GeV and 0:1◦ at 1 TeV. rently searching for hints of such gamma-ray emis- sion. The Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope constantly surveys the entire sky at energies of about 100 MeV up to a 3. Dark matter targets few hundred GeV. At higher energies, ground-based atmospheric Cherenkov telescopes perform searches for emission from specific targets. In this contribu- The particle physics governing the annihilation of tion I review the dark matter searches currently being cold dark matter is the same everywhere. Therefore, carried out by the Very Energetic Radiation Imaging any location in the Universe that hosts a suitably high Telescope Array System (VERITAS). dark matter density is a target for an indirect search. arXiv:1303.1406v1 [astro-ph.HE] 6 Mar 2013 The flux of gamma-rays from dark matter annihilation takes a surprisingly simple form 2. VERITAS dF hσvi dNγ (E) dJ(θ) (E; θ) = 2 : (1) The VERITAS array [Weekes et al. 2002, Holder dEdΩ 8πMχ dE dΩ et al. 2006] consists of four 12-meter imaging atmo- spheric Cherenkov telescopes (ACTs) located at the Here, dF=dEdΩ is the flux of gamma-rays per energy base of Mount Hopkins at the Fred Lawrence Whip- per solid angle, hσvi is the velocity-averaged cross ple Observatory in Arizona. Such telescopes use the section for dark matter self-annihilation and Mχ is the mass of the dark matter particle. The quantity dNγ =dE is the energy spectrum of gamma-rays arising from a single annihilation. The last factor quantifies ∗Electronic address: [email protected] the distribution of dark matter along the line of sight. eConf C121028 2 4th Fermi Symposium : Monterey, CA : 28 Oct-2 Nov 2012 1178 ACCIARI ET AL. Vol. 720 10−10The so-called J value is defined as 10 −20 Draco Willman I Z Bootes I 10 −21 dJ 2 Draco = dlρ ; Willman I (2) χ −22 Ursa Minor ) dΩ 10 1 −11 Ursa Minor − 10 Bootes I −23 s 10 2 where ρ is the dark matter density and the integral 1 χ − s −24 3 is taken along a particular line of sight. − 10 Note that the dependence on the environment is −25 10−12 10 completely captured in the J value; the other quanti- v> cm −26 σ 10 ties are universal properties of the dark matter parti- < 10 −27 dN/dE (ergs cm cle. In particular, the energy spectrum of annihilation 2 −13 E 10 −28 is governed by the annihilation channel (e.g annihila- 10 tion into quarks, leptons, photons, etc). Only at cos- −29 10 mological distances is this spectrum expected to be −30 −14 10 2 3 10 attenuated−1 by absorption or pair-production. 10 10 10 1 10 Mχ (GeV) Unlike charged particles,E (TeV) gamma-rays travel unde- flected from their point of emission. The excellent2 1 Figure 2. Exclusion regions in the (Mχ , σv )parameterspacebasedonthe Figure 1.angular95% CL upper resolution limits on theof spectral VERITAS energy densitytherefore (erg cm makes− s− ) it resultsFigure of 1: the observations. Upper limits It is on computed the dark# according$ matter to annihilation Equation (4)usinga as a function of gamma-ray energy for the four dSphs. compositecross section neutralino derived spectrum from (see VERITAS Wood et al. 2008 observations)andthevaluesof of fourJ from (A colora version powerful of this figure instrument is available because in the online it journal.) can target specific TableMilky1.BlackasterisksrepresentpointsfromMSSMmodelsthatfallwithin Way dwarf galaxies. Asterisks correspond to a locations with a high dark matter density (large J model3 standard scan deviations of MSSM of the parameter relic density space measured subject in the three-yearto a relic WMAP value). The Collaboration is currently focused on four data± set (Spergel et al. 2007). the dwarf galaxy center, ρ is the DM mass density, and λ is abundance constraint. Figure taken from Acciari et al. classes of dark matter targets: nearby dwarf galaxies, the line-of-sight distance to the differential integration volume. [2010]. the Milky Way galactic center, nearby galaxy clusters, The astrophysical contribution to the flux can be expressed by hour, significantly constraining upper limits on σ v could be E. ALIU et al. 4 # $ PHYSICAL REVIEW D 85, 062001 (2012) the dimensionlessand unidentified factor J sources in the Fermi catalog. established if J is on the order of 10 . Because the factor, J, is proportional to the DM6 density by Rolke et al. [85] is used in our analysis. As discussed 1 squared, it is subject to considerable uncertainty in its exper- J ( ) d ρ2(λ, ) dλ, (3) in the following sections, these ULs will serve for the 3.1.∆Ω Dwarf galaxies2 Ω Ω = ρ RH imental determination. For example, the mass of a DM4 halo is computation of integrated flux ULs and for constraining ! c " #∆Ω # 17 Milky Way dwarf spheroidal galaxies are dark- determined by the interaction of a galaxy with its neighbors and some dark matter models. To make the computation of which has been normalized to the product of the square of the is concentrated in the outer regions of the galaxy. Unlike the matter dominated systems30 that orbit3 inside the grav- 2 integrated flux ULs robust, we define a minimum energy, critical density, ρc 9.74 10− gcm− and the Hubble DM halo mass, the neutralino annihilation flux is determined] itational potential= of× the Galaxy [Walker 2012]. They [deg] above which the energy reconstruction bias is less than 5%. radius, RH 4.16 Gpc following Wood et al. (2008). by the inner regions of the galaxy where the density is highest. are= excellent targets for indirect detection because The energy reconstruction bias as a function of the recon- Based on Equation (2), the upper limits on the gamma- ForJ2000 these regions, the stellar kinematic data do not fully con- they are nearby, have large dark matter densities, and 0 ray rate, Rγ (95% CL), constrain the WIMP parameter space strain16 the DM density profile due to limited statistics. Various structed -ray energy has been studied with Monte Carlo perhaps most importantly, contain no known gamma- (mχ , σ v )accordingto parameterizations of the DM mass density profile have been put simulations. The -ray selection cuts used in this analysis # ray$ sources. Therefore, any detection of gamma-rays -2 forward (Navarro et al. 1997;Kazantzidisetal.2004;deBloksignificance [ set the minimum reconstructed energy to Emin Rfromγ (95% this CL class) of objectsJ becomesσ intriguingv evidence declination ¼ > # $ et al. 2001;Burkert1995)basedonempiricalfitsandstudies 300 GeV. The ULs on the number of -rays computed forhr the1 annihilation1.09 of10 dark4 3 matter.10 26 cm3 s 1 of simulated cold dark matter (CDM) halos. We adopt the as- with the Rolke prescription are displayed in Table I, along − ! − − " -4 In the absence of× a signal, constraints× 2 can be placed sumption15 of the NFW profile (Navarro et al.