Neutrino telescope searches for dark matter in the Sun

Pat Scott Fundamental Physics Section, Department of Physics, Imperial College London, Blackett Laboratory, Prince Consort Road, London SW7 2AZ, UK

I give a brief review of a few recent developments and future directions in the search for dark matter using high-energy neutrinos from the Sun. This includes the ability to recast neutrino telescope limits on nuclear scattering of dark matter to arbitrary new theories, and new calculations of the solar atmospheric background relevant to such searches. I also touch on applications to global searches for new physics, and prospects for improving searches for asymmetric dark matter in the Sun.

1 Current status

High-energy neutrinos from the Sun provide one of the cleanest potential discovery channels for weakly-interacting dark matter (DM). Weakly-interacting DM particles passing through the Sun are expected to scatter on solar nuclei. Some of these collisions reduce the kinetic energy of the DM particle enough for it to become gravitationally bound to the Sun, causing it to return on a bound orbit and undergo subsequent scattering, eventually thermalising and settling down to the solar core. If DM is able to annihilate, either with itself of with anti-DM captured in a similar manner, high-energy SM particles will be produced in the solar core. Even if neutrinos are not amongst those particles produced in the annihilation hard process, they will still be generated with quite high energies in the decay and subsequent interaction of other SM particles with nuclei in the Sun. Unlike the other SM particles, these GeV-scale neutrinos are then able to travel unhindered from the centre of the Sun to the surface, and across space to Earth, where they may be detected with terrestrial experiments. The directionality of the signal is the primary means by which it can be distinguished from the atmospheric neutrino background, caused by cosmic ray interactions with the Earth’s at- mosphere. The only known background to the signal is therefore the analogous production of high-energy neutrinos in the atmosphere of the Sun, due to interactions of cosmic rays with solar arXiv:1710.05190v1 [hep-ph] 14 Oct 2017 nuclei. The capture of dark matter by the Sun typically becomes the rate-limiting step in the pro- duction of any signal, rather than the annihilation. Searches for high-energy neutrinos from the Sun are therefore most useful for constraining the interaction cross-section of dark matter with nuclei. Spin-dependent interactions are particularly relevant, as the Sun consists mostly of hydrogen, which possesses nuclear spin. Current limits from neutrino telescope and direct searches for dark matter scattering are shown in Fig.1. The IceCube neutrino telescope presently provides the leading sensitivity to spin-dependent scattering with protons at high DM masses1, whereas Super-Kamiokande2 and PICO-603 have the leading sensitivity at low masses. ANTARES4 also provides complementary constraints. Direct searches lead the way for spin-independent interactions, and spin-dependent interactions with neutrons. -36 10 10 0 IceCube (2011-2014) -39 -3 Super-K (1996-2012) 10 10 Antares (2007-2012) IceCube (2011-2014) Super-K (1996-2012) -37 -1 PICO-2L PICO-60 10 -40 -4 10 (2015) (2013-2014) 10 Antares (2007-2012) 10

-41 10 10-5 -38 -2 ] 10 10 2

] -42 -6 2 10 10 [cm [pb] [cm p -39 -3 [pb] 10 p

10 − -43 -7

SI 10 χ

SD 10 σ

-44 -8 -40 10 10 10 10 -4 LUX (2013,2014-2016) -45 Xenon100 (2010-2014) 10 bb¯ 10-9 bb¯ PandaX-II (2016) + − + W W -41 W W -5 + − 10 + 10 τ τ -46 10 10-10 101 102 103 104 10 1 10 2 10 3 10 4 m [GeV] m [GeV] χ 0 36 1 37 10 10− 10− 10−

DAMA ] ] 2

2 2 38 1 XENON100 37 10− 10− 10− 10− LUX WS2013

3 XENON100 39 LUX WS2013+WS2014–16 10− LUX WS2013 10− PandaX-II PandaX-II 2 38 10− 10− LUX WS2013+WS2014–16 ¯ ) bb 4 40 10− 10− PICO-2L IceCube ( 3 39 10 ) 10 MSSM (Strege et al, 2014) − − − W + 5 41 10− 10− W ¯ SuperK ( ττ) ) τ¯ 4 τ 40 ATLAS: 10− IceCube ( 10− = 1 6 gq = 0.25,gχ MSSM (GAMBIT, 2017) 42 WIMP–proton cross section [pb] 10 10

WIMP–neutron cross section [pb] − − WIMP–proton cross section [cm WIMP–neutron cross section [cm IceCube ( PICO-60 CMS: gq = 0.25,gχ = 1 5 41 7 43 10− 10− 10− 10− 101 102 103 104 105 101 102 103 104 105 WIMP Mass [GeV/c2] WIMP Mass [GeV/c2]

Figure 1 – Current limits on nuclear scattering of dark matter from neutrino telescopes and direct detection. The top row shows limits from the leading neutrino telescopes on both the spin-dependent scattering cross-section with protons, and the spin-independent cross-section (with any nucleon). The lower row compares limits on the spin-dependent interactions with protons and neutrons, illustrating the role of various direct detection and collider experiments, as well as the evolution of corresponding supersymmetric theory predictions over time. Figures from IceCube1 (top row) and LUX5 (bottom row). 10-38

-39 F 8 10 C 3 ] 0 6 τ¯ 2 - τ CO ) m PI 6 c -40 1

[ 10 0 Super-K(2015) τ τ¯ (2 D p C S χ I σ 10-41

10-42 ττ¯ floor ττ¯ floor Super-K IceCube 10-43 101 102 103 104 mχ [ GeV ]

Figure 2 – Current limits on spin-dependent nuclear scattering of DM, compared to older calculations of the neutrino floor for neutrino searches toward the Sun; more accurate calculations of the floor are now available9,10, but give broadly similar results. Figure from Ng et al.11

2 Improved background calculations

Previous predictions of the background rate of high-energy neutrinos from the Sun, due to interactions of cosmic rays with nuclei in the solar atmosphere, were computed more than a decade ago.6,7,8 However, two more recent recalculations have appeared.9,10 Compared to the older predictions, the new calculations make use of modern knowledge on neutrino oscillations, production and interaction cross-sections. One of these10 also makes use of up-to-date models of the solar composition and structure, and carries out extensive Monte Carlo simulations of neutrino production, interaction and oscillation. Both studies (and another at the same time, based on the old flux estimates11) show that the solar atmospheric background lies barely an order of magnitude below current sensitivity limits for some models (Fig.2). This suggests that future neutrino telescopes might be able to directly measure this irreducible ‘neutrino floor’, and that the improved calculations of the background rates should be included in future phenomenological studies of DM scattering and annihilation in the Sun.

3 Recasting efforts

One of the major difficulties in interpreting the results of neutrino searches for DM in the Sun is the model-dependence of most published limits. This arises from the fact that the signal pre- diction is highly model-dependent, as the capture rate, annihilation rate, annihilation branching fractions and resulting neutrino spectrum predicted at Earth all enter the calculation of the predicted signal in a significant way. Traditional presentations from experiments1,4,15 give limits on scattering cross-sections as a function of DM mass, under various limiting assumptions about the dominant annihilation channel. Recently, a more general and flexible method for presenting the results of such searches has been developed16,12, allowing existing searches to be easily recast to provide detailed and consistent constraints on alternative DM models. This allows existing results to be converted to limits on different annihilation channels than those assumed in the original analysis (Fig.3), and for them to be applied to much more complex models (Fig.4), including arbitrary combinations of different annihilation final states and nuclear interactions (spin-dependent, spin-independent, and even more general forms). Public software exists to 33 IceCube Collaboration 2016 ) 10− 2 + gg W W − 34 (cm 10− b¯b ZZ ,p + SD 35 hh τ τ − σ 10− tt¯ νν¯ 36 10−

37 10−

38 10−

39 10−

40 10−

41 10− IC79

Dark matter-proton cross-section 42 10− 101 102 103 104 Dark matter mass mχ (GeV)

Figure 3 – 90%CL upper limits on the DM-proton nuclear scattering cross-section derived by recasting the IceCube 79-string search for DM in the Sun. Different curves assume different annihilation final states. Figure from IceCube.12

35 IceCube Collaboration 2016

) 10− 2

36 b ¯ 10− b (cm Area where essentially

,p τ 37 + all MSSM models are 10− τ 0 SD − excluded by IC79 σ 38 10− MSSM excluded by 39 10− IC79 depending on channel 40 10− Well-tempered MSSM models not 41 neutralinos 10− excluded by IC79

42 10− χ± coannihilation -proton cross-section 0 1 43 ˜ χ 10− B q˜ Pure h˜ MSSM-25 benchmarks coann.− 44 excluded at >90% CL A funnel 10− tension (68% – 90% CL) 45 10− allowed Effective 101 102 103 104 Dark matter mass mχ (GeV)

Figure 4 – 25-parameter supersymmetric benchmark models from Silverwood et al.13 and Cahill-Rowley et al.14, colour-coded according to how strongly disfavoured they are by the recast limit from the 79-string IceCube search + for DM annihilation in the Sun. The borders of the grey regions show the recast limits for typical ‘hard’ (τ τ −) and ‘soft’ (b¯b) spectra seen in supersymmetric models. Figure from IceCube.12 perform the recast (nulike: http://nulike.hepforge.org). This recastable form also has the advantage of providing substantial additional information compared to the traditional presentation, as it provides event-level information on neutrino/muon arrival angles and energies, and a detailed approximation to the full likelihood form of the experiment. This is particularly important when including searches for high-energy neutrinos from the Sun in global analyses of new physics scenarios such as supersymmetry. IceCube searches have been used in this format in the most recent reference global analyses of supersymmetric DM17,18 and scalar singlet DM19, done in the context of the GAMBIT project.20,21,22,23,24,25

4 Prospects for combined analysis and application to new models

Whilst the recastable form12 of the 79-string IceCube data15 is far more flexible and generally useful to the phenomenological community, there is some amount of overhead required to reduce the data to the necessary form. For this reason, recastable versions of the 86-string IceCube search1, and the latest results from ANTARES4 and Super-Kamiokande2 are not yet available. It is hoped that this will soon change. One significant driver for such a development is the prospect that the data of all three neutrino telescopes could be seamlessly combined, to give a single unified and strengthened limit. Indeed, as soon as each of the individual datasets is available in recastable form, the combination would be extremely straightforward to perfom via the composite likelihood method. Including the combined constraint in global analyses of searches for new physics would be similarly straightforward. A related recasting application will be to rigorously apply neutrino telescope limits to models of asymmetric DM that exhibit both symmetric and asymmetric components26, allowing strong constraints to be placed on their asymmetry parameter r . Here the capture rates of such models need to be carefully determined via a low-energy effective∞ operator treatment27, and the full range of possible operators for both scattering and annihilation, along with their interferences, taken into account. The results of such an exercise will be especially interesting to compare to helioseismological and low-energy solar neutrino observables, given recent suggestions of a possible signal of DM from this sector. 28,29,30,31,32,33,a

Acknowledgments

I am supported by STFC (ST/K00414X/1, ST/P000762/1, ST/L00044X/1), and thank my co- authors on a number of the works discussed here.

1. M. G. Aartsen, M. Ackermann, et. al., Search for annihilating dark matter in the Sun with 3 years of IceCube data, Eur. Phys. J. C 77 (2017) 146, [arXiv:1612.05949]. 2. K. Choi, K. Abe, et. al., Search for Neutrinos from Annihilation of Captured Low-Mass Dark Matter Particles in the Sun by Super-Kamiokande, Phys. Rev. Lett. 114 (2015) 141301. 3. C. Amole, M. Ardid, et. al., Dark Matter Search Results from the PICO -60 C3F8 Bubble Chamber, Phys. Rev. Lett. 118 (2017) 251301, [arXiv:1702.07666]. 4. S. Adrián-Martínez, A. Albert, et. al., Limits on dark matter annihilation in the sun using the ANTARES neutrino telescope, Phys. Lett. B 759 (2016) 69–74, [arXiv:1603.02228]. 5. LUX Collaboration, D. S. Akerib, et. al., Limits on spin-dependent WIMP-nucleon cross section obtained from the complete LUX exposure, Phys. Rev. Lett. 118 (2017) 251302, [arXiv:1705.03380]. 6. G. Ingelman and M. Thunman, High energy neutrino production by cosmic ray interactions in the Sun, Phys. Rev. D 54 (1996) 4385–4392, [hep-ph/9604288].

aSee also A. Vincent, these proceedings. 7. C. Hettlage, K. Mannheim, and J. G. Learned, The sun as a high energy neutrino source, Astropart. Phys. 13 (2000) 45–50, [astro-ph/9910208]. 8. G. L. Fogli, E. Lisi, A. Mirizzi, D. Montanino, and P. D. Serpico, Oscillations of solar atmosphere neutrinos, Phys. Rev. D 74 (2006) 093004, [hep-ph/0608321]. 9. C. A. Argüelles, G. de Wasseige, A. Fedynitch, and B. J. P. Jones, Solar atmospheric neutrinos and the sensitivity floor for solar dark matter annihilation searches, JCAP 7 (2017) 024, [arXiv:1703.07798]. 10. J. Edsjö, J. Elevant, R. Enberg, and C. Niblaeus, Neutrinos from cosmic ray interactions in the Sun, JCAP 6 (2017) 033, [arXiv:1704.02892]. 11. K. C. Y. Ng, J. F. Beacom, A. H. G. Peter, and C. Rott, Solar Atmospheric Neutrinos: A New Neutrino Floor for Dark Matter Searches, arXiv:1703.10280. 12. IceCube Collaboration: M. G. Aartsen et. al., Improved limits on dark matter annihilation in the Sun with the 79-string IceCube detector and implications for supersymmetry, JCAP 04 (2016) 022, [arXiv:1601.00653]. 13. H. Silverwood, P. Scott, et. al., Sensitivity of IceCube-DeepCore to neutralino dark matter in the MSSM-25, JCAP 3 (2013) 27, [arXiv:1210.0844]. 14. M. W. Cahill-Rowley, J. L. Hewett, A. Ismail, M. E. Peskin, and T. G. Rizzo, pMSSM Benchmark Models for Snowmass 2013, arXiv:1305.2419. 15. IceCube Collaboration: M. G. Aartsen, R. Abbasi, et. al., Search for Dark Matter Annihilations in the Sun with the 79-String IceCube Detector, Phys. Rev. Lett. 110 (2013) 131302, [arXiv:1212.4097]. 16. P. Scott, C. Savage, J. Edsjö, and the IceCube Collaboration: R. Abbasi et al., Use of event-level neutrino telescope data in global fits for theories of new physics, JCAP 11 (2012) 57, [arXiv:1207.0810]. 17. GAMBIT Collaboration: P. Athron, C. Balázs, et. al., Global fits of GUT-scale SUSY models with GAMBIT, Eur. Phys. J. C in press (2017) [arXiv:1705.07935]. 18. GAMBIT Collaboration: P. Athron, C. Balázs, et. al., A global fit of the MSSM with GAMBIT, Eur. Phys. J. C in press (2017) [arXiv:1705.07917]. 19. GAMBIT Collaboration: P. Athron, C. Balázs, et. al., Status of the scalar singlet dark matter model, Eur. Phys. J. C 77 (2017) 568, [arXiv:1705.07931]. 20. GAMBIT Collaboration: P. Athron, C. Balazs, et. al., GAMBIT: The Global and Modular Beyond-the-Standard-Model Inference Tool, Eur. Phys. J. C in press (2017) [arXiv:1705.07908]. 21. GAMBIT Collider Workgroup: C. Balázs, A. Buckley, et. al., ColliderBit: a GAMBIT module for the calculation of high-energy collider observables and likelihoods, Eur. Phys. J. C in press (2017) [arXiv:1705.07919]. 22. GAMBIT Dark Matter Workgroup: T. Bringmann, J. Conrad, et. al., DarkBit: A GAMBIT module for computing dark matter observables and likelihoods, Eur. Phys. J. C in press (2017) [arXiv:1705.07920]. 23. GAMBIT Models Workgroup: P. Athron, C. Balázs, et. al., SpecBit, DecayBit and PrecisionBit: GAMBIT modules for computing mass spectra, particle decay rates and precision observables, submitted to Eur. Phys. J. C (2017) [arXiv:1705.07936]. 24. GAMBIT Flavour Workgroup: F. U. Bernlochner, M. Chrząszcz, et. al., FlavBit: A GAMBIT module for computing flavour observables and likelihoods, Eur. Phys. J. C in press (2017) [arXiv:1705.07933]. 25. GAMBIT Scanner Workgroup: G. D. Martinez, J. McKay, et. al., Comparison of statistical sampling methods with ScannerBit, the GAMBIT scanning module, Eur. Phys. J. C in press (2017) [arXiv:1705.07959]. 26. K. Murase and I. M. Shoemaker, Detecting asymmetric dark matter in the Sun with neutrinos, Phys. Rev. D 94 (2016) 063512, [arXiv:1606.03087]. 27. R. Catena and B. Schwabe, Form factors for dark matter capture by the Sun in effective theories, JCAP 4 (2015) 042, [arXiv:1501.03729]. 28. M. T. Frandsen and S. Sarkar, Asymmetric Dark Matter and the Sun, Phys. Rev. Lett. 105 (2010) 011301, [arXiv:1003.4505]. 29. M. Taoso, F. Iocco, G. Meynet, G. Bertone, and P. Eggenberger, Effect of low mass dark matter particles on the Sun, Phys. Rev. D 82 (2010) 083509, [arXiv:1005.5711]. 30. A. C. Vincent and P. Scott, Thermal conduction by dark matter with velocity and momentum-dependent cross-sections, JCAP 4 (2014) 19, [arXiv:1311.2074]. 31. A. C. Vincent, P. Scott, and A. Serenelli, Possible Indication of Momentum-Dependent Asymmetric Dark Matter in the Sun, Phys. Rev. Lett. 114 (2015) 081302, [arXiv:1411.6626]. 32. A. C. Vincent, A. Serenelli, and P. Scott, Generalised form factor dark matter in the Sun, JCAP 8 (2015) 40, [arXiv:1504.04378]. 33. A. C. Vincent, P. Scott, and A. Serenelli, Updated constraints on velocity and momentum-dependent asymmetric dark matter, JCAP 11 (2016) 007, [arXiv:1605.06502].