The VERITAS Dark Matter Program
In the cosmological paradigm, cold dark matter (DM) dominates the mass content of the Universe and is present at every scale. Candidates for DM include many extensions of the Standard Model with weakly interacting massive particles (WIMPs) in the mass range from 10 GeV to greater than 10 TeV. The self-annihilation or decay of WIMPs in astrophysical regions of high DM density can produce secondary particles including very-high-energy (VHE; E¿100 GeV) gamma rays. VERITAS, an array of atmospheric Cherenkov telescopes, sensitive to VHE gamma rays in the 85 GeV-30 TeV energy range, has been utilized for indirect DM searches. Astrophysical candidates for indirect DM detection by VERITAS are dwarf spheroidal galaxies (dSphs) of the Local Group and the Galactic Center, among others. This is a report on the current status of VERITAS observations of dark matter targets and the results of a combined search of four dSphs.
VERITAS Dark Matter Program
\FullConference35th International Cosmic Ray Conference â ICRC217-
10-20 July, 2017
Bexco, Busan, Korea
The search for standard model particles resulting from the annihilation of dark matter provides an important complement to the efforts of direct searches for dark matter interactions and searches for dark matter production at particle accelerators. Among the theoretical candidates for the dark matter particle above a few GeV, Weakly Interacting Massive Particles (WIMPs) are well motivated [1, 2] as they naturally provide the measured present day cold dark matter density [3, 4, 5, 7, 6]. In such models, the WIMPs either decay or annihilate into standard model particles that produce mono-energetic gamma-ray lines and/or a continuum of gamma rays with energies up to the dark matter particle mass.
Good targets for searches with gamma-ray instruments are those with a large inferred dark matter (DM) density and low astrophysical backgrounds; i.e. those that would not produce gamma rays through conventional processes. Dwarf Spheroidal Galaxies (dSphs) fit these criteria since they are close by (tens to hundreds of kpc away) and with no theoretical basis for emitting very-high-energy (VHE) gamma rays. The center of our galaxy is also a tantalizing target since it is even closer, about 8 kpc away. However the center of the galaxy has many astrophysical backgrounds including diffuse emission along the galactic plane ; the radio source at the center of our galaxy, SgrA*, is also a very bright VHE emitter. The metric for a good DM target is the -factor (or factor in the case of decay) defined as:
where is the DM density and is the direction along the line-of-sight. This work will focus on DM annihilations. We will also refer to the -profile which is equal to -factor per solid angle. The -factor is directly proportional to the expected gamma-ray flux from DM annihilations,
where is the WIMP particle mass, is the energy spectrum of a single annihilation with branching ratio , and is the velocity-averaged annihilation cross section.
VERITAS (Very Energetic Radiation Imaging Telescope Array System) is an array of four imaging atmospheric Cherenkov telescopes (IACTs), each 12 m in diameter, located at the Fred Lawrence Whipple Observatory in southern Arizona, USA. VERITAS is sensitive to gamma rays from approximately 85 GeV (after camera upgrade) to greater than 30 TeV with a typical energy resolution of and an angular resolution (68% containment) of 0.1 degrees per event. The flux sensitivity of the standard analysis is such that a source with a flux of order of 1% of the Crab Nebula flux can be detected in approximately 25 hours of observation.
2 Combined Analysis of Dwarf Galaxies from 2007 to 2013
From the beginning of four-telescope operations in 2007 to the summer of 2013, five dSphs in the northern hemisphere have been observed: Segue 1, Ursa Minor, Draco, Boötes and Willman 1. No significant gamma-ray excess was observed in the direction of any of the five dSphs. A dark matter search was performed with 216 hours of VERITAS data, including 92 hours on the dSph Segue 1. Table 1 summarizes the analysis results of each of the five dSphs.
The 95% confidence level limits on the thermally-averaged annihilation cross section are shown in Figure 1. The search for DM annihilation and the limits on the cross section were derived from the Event Weighting method presented in Geringer Sameth et. al. . Willman 1 was excluded from the combined search and limits since the -factor can not be calculated reliably due to irregular kinematics . Additionally, Bonnivard et. al.  have pointed out the possibility of contamination of the stellar samples used to perform the Jeans analysis which could bias -profiles towards larger annihilation signals, possibly by orders of magnitude. Therefore, two different limits are shown in each panel in Figure 1: including Segue 1 and without. Each panel in Figure 1 represents a 100% branching ratio into a given standard model particle. The limits are presented as a band with the representing the 1 systematic uncertainty of the -profile.
3 A Decade of Dwarf Galaxy Observations: 2007 to 2017
The indirect detection of dark matter is one of the collaboration’s key science projects, and dwarf galaxies are an integral part of the program. Because VERITAS is not a survey instrument, we have to prioritize the DM targets within our observing budget. DSphs observed by VERITAS are subdivided into two categories: ‘Deep Exposure’ dSphs and ‘Survey’ dSphs. ‘Deep Exposure’ dSphs are those with the best -factors in the literature which get the majority of the observing time. This time is typically split between the ultra-faint and classical dSphs. This is done as a way to ensure that the overall DM program is not severely penalized in case one particular dSph falls out of favor as a good indirect DM target. Most of the dSphs that do not have deep exposures are in the survey for much shorter exposures, typically a few hours per year. The survey provides DM sensitivity to nearly all dSphs in the northern hemisphere as a form of insurance in case misunderstanding of velocity dispersion data leads to a dSph having a much larger -factor than otherwise indicated. The survey and the deep exposure dSphs can be also be combined or ‘stacked’ into a more sensitive limit than any one individual dSph as shown in the previous section.
Table 2 lists all dSph observations since the beginning of four-telescope operations in 2007 to mid-2017. The exposure accumulated on for each dSph are divided into observing epochs: the original array configuration (v4), after the move of Telescope 1 (v5) and after the camera upgrade (v6). See  for details on the performance of the different epochs. The observation times in Table 2 do not include time intervals with bad weather or hardware malfunctions, but rather the total observing time so these numbers should be considered an upper limit to the possible exposure of each dSph. The analysis of many of these dSph data sets is still preliminary.
3.1 Ursa Major II
As shown in Table 2, a large amount of the total dSph observing time after the camera upgrade has been devoted to Ursa Major II. Ursa Major II was discovered in SDSS data with follow-up observations from Subaru . It is one of the faintest known dSphs, with velocity dispersion measurements made with only 20 stars , but with a large mass-to-luminosity ratio. -factors for Ursa Major II typically indicate that it is one of the highest, often similar to or even exceeding Segue 1 (e.g. [10, 11].
VERITAS observed Ursa Major II between late 2013 and spring 2017. After dead-time corrections and removal of time intervals with bad weather or hardware performance, a live time of 145 hours remain. Observation runs with fewer than all four telescopes were also removed from this analysis. The analysis of the data follows the same procedure as in . The hFit algorithm is used for event reconstruction  and the crescent background method is used to subtract the remaining cosmic-ray background . An excess of 1 above the background was found at the Ursa Major II position. We report a flux upper limit of 1.410ms above 200 GeV at a 95% confidence level using the method of Rolke  assuming a powerlaw with index of -2.4.
Figure 2 shows preliminary annihilation cross section limits for Ursa Major II. A version of the full maximum likelihood introduced by Aleksić et. al.  was used for the search for DM annihilations and for determining the limit to the velocity-averaged cross section. The -profile is convolved with the VERITAS PSF, generated from monte carlo simulations has described in . We do this because the DM profile is wider than the PSF, particularly at lower energies. The median and 1 -profile was provided in the supplemental material of Geringer-Sameth et al. . The annihilation spectrum is the one used from PPPC 4 DM ID . As shown in Figure 2, the VERITAS limit for Ursa Major II is more constraining on the cross section than the 216 hrs combined dSph  limit at every DM mass by a factor of 30% to more than a factor of 2, depending on DM mass.
|Dwarf||Exposure v4||Exposure v5||Exposure v6||Total Exposure|
|[GeV cm]||[GeV cm]||[GeV cm]||[min]||[min]||[min]||[min]|
|Ursa Major II||19.4||19.9||18.4||0||0||10869||10869|
|Canes Ver. II||17.7||18.5||17.0||0||0||864||864|
|Canes Ver. I||17.4||17.5||17.6||0||0||850||850|
|Ursa Major I||17.9||18.7||17.6||0||0||482||482|
4 Galactic Center
The center of our galaxy harbors a 410 solar mass black hole, coincident with the radio source SgrA*. The region around the center of our galaxy is very complex in gamma rays. This region includes bright sources such as the supernova remnant G0.6+0.1 , diffuse emission along the plane as observed by HESS  and Fermi-LAT as well as the extended Fermi bubbles . The existence of several possible astrophysical backgrounds makes finding a DM signal from WIMP annihilation a complicated but not impossible prospect. Because VERITAS is a Northern Hemisphere instrument, the Galactic Center never transits above 30 elevation. Observing at these elevations has the effect of raising the energy threshold and widening the PSF, but increasing sensitivity for higher energy showers. The degradation of the PSF is mostly offset by the displacement method of shower reconstruction . VERITAS is therefore the most sensitive instrument for observing the Galactic Center at energies greater than 2 TeV, with a PSF of 0.12 (68% containment).
5 Dark Matter Subhalos
Cosmological N-body simulations  indicate that DM halos contain significant substructures . Because of tidal disruption near the Galactic disk, most of the sub-halos are thought to survive at high Galactic latitude. The lack of material in these regions prevents the DM over-densities from attracting enough baryonic matter to trigger star formation. DM clumps would therefore be invisible to most astronomical observations from radio to X-rays. DM structures residing in the the Milky Way halo could potentially be nearby the Sun and therefore have a bright gamma-ray annihilation signal. These clumps would therefore only be visible at gamma-ray energies and not been in any other wavelength.
The Second Fermi-LAT Catalog (2FGL) contains 576 high-energy gamma-ray sources detected by the LAT after the first 24 months of observations without any clear associations in any other wavelength. To identify the unassociated Fermi objects (UFOs) most likely to be DM clump candidates we adapt the selection criteria from , requiring that they are located outside of the Galactic Plane, that they be non-variable, they exhibit power law spectra and that they have no known counterparts. The two best candidates were selected for VERITAS observations: 2FGL J0545.6+6018 and 2FGL J1115.0-0701. It should be noted that at the time these objects were proposed for observation, that the 2FGL catalog was the most up-to-date Fermi-LAT catalog available. 2FGL J1115.0-0701 was observed for 13.8 hours and 2FGL J0545.6+6018 was observed for 8.5 hours. No emission was detected by VERITAS, giving flux upper limits 0.5% and 0.6% of the Crab nebula flux above 250 GeV, respectively. Later analysis of 2FGL J1115.0-0701 in Fermi-LAT showed that the flux was variable, excluding the possibility of it being a DM subhalo. Updated Fermi-LAT analysis of 2FGL J0546.6+6018 shows that a curved power law is actually favored over a simple power law, but DM annihilation to quarks or bosons also provide a good spectral fit. Optical follow-up observations were performed within the 95% error region of 2FGL J0546.6+6018 by the 1.3m McGraw-Hill telescope located at Kitt Peak, AZ, USA. A previously undiscovered bright UV source was observed within the error circle, the analysis and interpretation of this is still ongoing. Swift-XRT observations show no new associations .
6 Summary and Future Work
VERITAS has an active Dark Matter and Astroparticle group that leads this key science project of indirect DM detection for the collaboration. Continued observations of dSphs and the Galactic Center region are a crucial part of the VERITAS long term plan . Unfortunately, no potential subhalo candidates have been found in the 3FGL catalog. The Galactic Center region could potentially be the best indirect -ray source for dark matter discovery or stringent limits on annihilation cross sections at the highest masses. Future work for dwarf galaxies involve a new combined analysis with not only new data, but with other -ray instruments, namely Fermi-LAT and HAWC, which is still in the preliminary stages. New analysis techniques such as boosted decision trees (BDTs)  and template reconstruction  have potential for improving the overall dark matter sensitivity of the experiment. Additionally many dSphs are extended for IACTs, so an extended source analysis could give a boost to DM sensitivity by as much as a factor of two in some cases.
This research is supported by grants from the U.S. Department of Energy Office of Science, the U.S. National Science Foundation and the Smithsonian Institution, and by NSERC in Canada. We acknowledge the excellent work of the technical support staff at the Fred Lawrence Whipple Observatory and at the collaborating institutions in the construction and operation of the instrument. The VERITAS Collaboration is grateful to Trevor Weekes for his seminal contributions and leadership in the field of VHE gamma-ray astrophysics, which made this study possible.
- Jungman, G., Kamionkowski, M., & Griest, K. 1996, physrep, 267, 195
- Servant, G., & Tait, T. M. P. 2003, Nuclear Physics B, 650, 391
- Planck Collaboration, Ade, P. A. R., Aghanim, N., et al. 2014, aap, 571, A16
- Steigman, G. 1979, Annual Review of Nuclear and Particle Science, 29, 313
- Zeldovic, Y. B., Okun, L. B., & Pikelner, S. B. 1965, PhL, 17, 164
- Chiu, H.-Y. 1966, Physical Review Letters, 17, 712
- Zel’dovich, Y. B. 1965, Advances in Astronomy and Astrophysics, ed. Z. Kopal, Vol. 3 (Academic Press), 242
- Aharonian, F., Akhperjanian, A. G., Bazer-Bachi, A. R., et al. 2006, Nature, 439, 695
- Geringer-Sameth, A., Koushiappas, S. M., & Walker, M. G. 2015, Phys. Rev. D, 91, 083535
- Geringer-Sameth, A., Koushiappas, S. M., & Walker, M. 2015, apj, 801, 74
- Bonnivard, V., Combet, C., Maurin, D., & Walker, M. G. 2015, MNRAS, 446, 3002
- Li, T.-P., & Ma, Y.-Q. 1983, apj, 272, 317
- Rolke, W. A., López, A. M., & Conrad, J. 2005, Nuclear Instruments and Methods in Physics Research A, 551, 493
- Park, N., & for the VERITAS Collaboration 2015, arXiv:1508.07070
- Zucker, D. B., Belokurov, V., Evans, N. W., et al. 2006, apjl, 650, L41
- McConnachie, A. W. 2012, aj, 144, 4
- Christiansen, J., & VERITAS Collaboration. 2012, in American Institute of Physics Conference Series, Vol. 1505, American Institute of Physics
- Aleksić, J., Ansoldi, S., Antonelli, L. A., et al. 2014, jcap, 2, 008
- Archambault et al. (The VERITAS Collaboration) Phys. Rev. D 95, 082001
- Cirelli, M., Corcella, G., Hektor, A., et al. 2011, jcap, 3, 051
- Archer, A., Benbow, W., Bird, R., et al. 2016, apj, 821, 129
- Su, M., Slatyer, T. R., & Finkbeiner, D. P. 2010, apj, 724, 1044
- Springel, V., et al. 2008, MNRAS, 391, 1685
- Diemand, J., et al. 2008, Nature, 454, 735
- Nieto, D., et al., 2011, 32nd ICRC, arXiV:1109.5935
- Stroh, M. C., & Falcone, A. D. 2013, apjs, 207, 28
- Staszak, D., & VERITAS Collaboration 2015, 34th International Cosmic Ray Conference (ICRC2015), 34, 868
- Krause, M., Pueschel, E., & Maier, G. 2017, Astroparticle Physics, 89, 1
- Christiansen, J. & the VERITAS Collaboration, these proceedings