SNEWS 2.0: A Next-Generation SuperNova Early Warning System for Multi-messenger Astronomy - arXiv.org

Page created by Juanita Mcdaniel
 
CONTINUE READING
SNEWS 2.0: A Next-Generation SuperNova Early Warning System for Multi-messenger Astronomy - arXiv.org
SNEWS 2.0: A Next-Generation SuperNova Early
                                              Warning System for Multi-messenger Astronomy
arXiv:2011.00035v2 [astro-ph.HE] 8 Jan 2021

                                                     S. Al Kharusi1 , S. Y. BenZvi2 , J. S. Bobowski3 , W. Bonivento4 ,
                                                     V. Brdar5,6,7 , T. Brunner1,8 , E. Caden9 , M. Clark10 ,
                                                     A. Coleiro11 , M. Colomer-Molla11,12 , J. I. Crespo-Anadón13 ,
                                                     A. Depoian10 , D. Dornic14 , V. Fischer15 , D. Franco11 ,
                                                     W. Fulgione16 , A. Gallo Rosso17 , M. Geske18 , S. Griswold2 ,
                                                     M. Gromov19,20 , D. Haggard1 , A. Habig21 , O. Halim22 ,
                                                     A. Higuera23 , R. Hill17 , S. Horiuchi24 , K. Ishidoshrio25 ,
                                                     C. Kato26 , E. Katsavounidis27 , D. Khaitan2 , J. P. Kneller28 ,
                                                     A. Kopec10 , V. Kulikovskiy29 , M. Lai30,31 , M. Lamoureux32 ,
                                                     R. F. Lang10 , H. L. Li33 , M. Lincetto14 , C. Lunardini34 ,
                                                     J. Migenda35 , D. Milisavljevic10 , M. E. McCarthy2 ,
                                                     E. O’Connor36 , E. O’Sullivan37 , G. Pagliaroli38 , D. Patel39 ,
                                                     R. Peres40 , B. W. Pointon41,8 , J. Qin10 , N. Raj8 , A. Renshaw42 ,
                                                     A. Roeth43 , J. Rumleskie17 , K. Scholberg43 , A. Sheshukov20 ,
                                                     T. Sonley9 , M. Strait44 , V. Takhistov45 , I. Tamborra46 ,
                                                     J. Tseng47 , C.D. Tunnell23 , J. Vasel48 , C. F. Vigorito49 ,
                                                     B. Viren50 , C. J. Virtue17 , J. S. Wang47 , L. J. Wen33 ,
                                                     L. Winslow27 , F. L. H. Wolfs2 , X. J. Xu7 and Y. Xu23
                                                     1
                                                       Department of Physics, McGill University, Montréal, QC H3A 2T8, Canada
                                                     2
                                                       Department of Physics and Astronomy, University of Rochester, Rochester, NY
                                                     14627, USA
                                                     3
                                                       Physics, University of British Columbia, Kelowna, BC V1V 1V7, Canada
                                                     4
                                                       INFN sezione di Cagliari Istituto Nazionale, Complesso Universitario di Monserrato
                                                     - S.P. per Sestu Km 0.700, I-09042 Monserrato (Cagliari), Italy
                                                     5
                                                       Fermi National Accelerator Laboratory, Batavia, IL, 60510, USA
                                                     6
                                                       Northwestern University, Department of Physics & Astronomy, 2145 Sheridan
                                                     Road, Evanston, IL 60208, USA
                                                     7
                                                       Max-Planck-Institut für Kernphysik, Postfach 103980, D-69029 Heidelberg,
                                                     Germany
                                                     8
                                                       TRIUMF, Vancouver, BC V6T 2A3, Canada
                                                     9
                                                       SNOLAB, Creighton Mine #9, 1039 Regional Road 24, Sudbury ON P3Y 1N2,
                                                     Canada
                                                     10
                                                        Department of Physics and Astronomy, Purdue University, West Lafayette, IN
                                                     47907, USA
                                                     11
                                                        Université de Paris, CNRS, AstroParticule et Cosmologie, F-75013, Paris, France
                                                     12
                                                        Instituto de Fı́sica Corpuscular (CSIC - Universitat de València) c/ Catedrático
                                                     José Beltrán, 2 E-46980 Paterna, Valencia, Spain
                                                     13
                                                        Department of Physics, Columbia University, New York, NY 10027, USA
SNEWS 2.0: A Next-Generation SuperNova Early Warning System for Multi-messenger Astronomy - arXiv.org
SNEWS 2.0: A Next-Generation SuperNova Early Warning System for Multi-messenger Astronomy2
            14
               Aix Marseille Univ, CNRS/IN2P3, CPPM, Marseille, France
            15
               Department of Physics, University of California at Davis, Davis, CA 95616, U.S.A.
            16
               OATo Torino & INFN sezione dei Laboratori Nazionali del Gran Sasso (LNGS),
            Assergi, Italy
            17
               Department of Physics, Laurentian University, Sudbury ON P3E 2C6, Canada
            18
               Department of Physics, Gonzaga University, Spokane, WA 99258, USA
            19
               Lomonosov Moscow State University Skobeltsyn Institute of Nuclear Physics,
            119234 Moscow, Russia
            20
               Joint Institute for Nuclear Research, 141980 Dubna, Russia
            21
               Department of Physics and Astronomy, University of Minnesota Duluth, Duluth,
            MN, 55812, USA
            22
               Dipartimento di Fisica, Universit‘a di Trieste, & INFN sezione di Trieste, I-34127
            Trieste, Italy
            23
               Departments of Physics, Astronomy, and Computer Science, Rice University, 6100
            Main St, Houston, TX, 77005, USA
            24
               Center for Neutrino Physics, Department of Physics, Virginia Tech, Blacksburg,
            VA 24061, USA
            25
               Research Centre for Neutrino Science, Tohoku University, Sendai 980-8578, Japan
            26
               Department of Aerospace Engineering, Tohoku University, Sendai 980-8579, Japan
            27
               Massachusetts Institute of Technology, Cambridge, MA 02139, USA
            28
               Department of Physics, NC State University, Raleigh, NC 27695, USA
            29
               INFN Sezione di Genova, Via Dodecaneso 33, Genova, 16146 Italy
            30
               Department of Physics, Cagliari University, Cagliari, CA 09127, Italy
            31
               Istituto Nazionale di Fisica Nucleare INFN
            32
               INFN Sezione di Padova & Università di Padova, Dipartimento di Fisica, Padova,
            Italy
            33
               Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100049,
            China
            34
               Department of Physics, Arizona State University, Tempe, AZ 85287-1504, USA
            35
               Department of Physics, King’s College London, London WC2R 2LS, United
            Kingdom
            36
               Department of Astronomy and The Oskar Klein Centre, Stockholm University,
            AlbaNova, 109 61, Stockholm, Sweden
            37
               Deptartment of Physics and Astronomy, Uppsala University, Box 516, S-75120
            Uppsala, Sweden
            38
               Gran Sasso Science Institute (GSSI) & INFN sezione di Laboratori Nazionali del
            Gran Sasso (LNGS), Assergi, Italy
            39
               Department of Physics, University of Regina, Regina, SK S4S 0A2, Canada
            40
               Physik-Institut, Universität Zürich, Zürich, Switzerland
            41
               Department of Physics, British Columbia Institute of Technology, Burnaby, BC,
            V5G 3H2, Canada
            42
               Department of Physics, University of Houston, Houston, TX 77204, USA
            43
               Department of Physics, Duke University, Durham, NC 27708, USA
            44
               School of Physics and Astronomy, University of Minnesota Twin Cities,
            Minneapolis, Minnesota 55455, USA
            45
               Department of Physics, University of California Los Angeles, Los Angeles, U.S.A.
            46
               Niels Bohr International Academy and DARK, Niels Bohr Institute, Blegdamsvej
            17, 2100 Copenhagen, Denmark
            47
               Oxford University, Denys Wilkinson Building, Keble Road, Oxford OX1 3RH, UK
            48
               Department of Physics, Indiana University, Bloomington, IN 47405, USA
SNEWS 2.0: A Next-Generation SuperNova Early Warning System for Multi-messenger Astronomy3
            49
               Department of Physics, University of Torino & INFN, via Pietro Giuria 1, 10125
            Torino, Italy
            50
               Physics Department, Brookhaven National Laboratory, Upton, NY 11973, USA
SNEWS 2.0: A Next-Generation SuperNova Early Warning System for Multi-messenger Astronomy4

             Abstract. The next core-collapse supernova in the Milky Way or its satellites will
             represent a once-in-a-generation opportunity to obtain detailed information about the
             explosion of a star and provide significant scientific insight for a variety of fields because
             of the extreme conditions found within. Supernovae in our galaxy are not only rare
             on a human timescale but also happen at unscheduled times, so it is crucial to be
             ready and use all available instruments to capture all possible information from the
             event. The first indication of a potential stellar explosion will be the arrival of a
             bright burst of neutrinos. Its observation by multiple detectors worldwide can provide
             an early warning for the subsequent electromagnetic fireworks, as well as signal to
             other detectors with significant backgrounds so they can store their recent data. The
             Supernova Early Warning System (SNEWS) has been operating as a simple coincidence
             between neutrino experiments in automated mode since 2005. In the current era
             of multi-messenger astronomy there are new opportunities for SNEWS to optimize
             sensitivity to science from the next Galactic supernova beyond the simple early alert.
             This document is the product of a workshop in June 2019 towards design of SNEWS 2.0,
             an upgraded SNEWS with enhanced capabilities exploiting the unique advantages of
             prompt neutrino detection to maximize the science gained from such a valuable event.

Submitted to: New J. Phys.
CONTENTS                                                                                                                        5

Contents

1 Introduction                                                                      7
  1.1 Current Configuration (SNEWS 1.0) . . . . . . . . . . . . . . . . . . . . 8
  1.2 SNEWS 2.0 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
  1.3 Lowering the Threshold . . . . . . . . . . . . . . . . . . . . . . . . . . . 10

2 Stellar Core Collapse Signals                                                                                                11
  2.1 Neutrinos from Core Collapse Supernovae . . .                    .   .   .   .   .   .   .   .   .   .   .   .   .   .   12
  2.2 Gravitational Wave Signals from Core Collapse                    .   .   .   .   .   .   .   .   .   .   .   .   .   .   13
  2.3 Electromagnetic Signals . . . . . . . . . . . . .                .   .   .   .   .   .   .   .   .   .   .   .   .   .   14
  2.4 Other Transients . . . . . . . . . . . . . . . . .               .   .   .   .   .   .   .   .   .   .   .   .   .   .   14
       2.4.1 Type Ia Supernovae . . . . . . . . . . . .                .   .   .   .   .   .   .   .   .   .   .   .   .   .   14
       2.4.2 Pair Instability Supernovae . . . . . . .                 .   .   .   .   .   .   .   .   .   .   .   .   .   .   15
       2.4.3 Compact Object Mergers . . . . . . . . .                  .   .   .   .   .   .   .   .   .   .   .   .   .   .   15

3 Pointing to the Supernova with Neutrinos                                                                                     16
  3.1 Anisotropic Interactions . . . . . . . . . .         .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   17
      3.1.1 Water Cherenkov . . . . . . . . . .            .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   17
      3.1.2 Liquid Argon . . . . . . . . . . . .           .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   19
      3.1.3 Liquid Scintillator . . . . . . . . . .        .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   19
  3.2 Triangulation . . . . . . . . . . . . . . . .        .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   20

4 Presupernova Neutrinos                                                                                                       24

5 The   SNEWS Alert and Followup                                                                                               26
  5.1   Real-Time Algorithmics . . . . . . . . . . . . . . . . .                       .   .   .   .   .   .   .   .   .   .   27
  5.2   Multimessenger Follow-Up . . . . . . . . . . . . . . . .                       .   .   .   .   .   .   .   .   .   .   27
  5.3   Alert Broadcasting and Optimized Observing Strategies                          .   .   .   .   .   .   .   .   .   .   29
  5.4   Latency . . . . . . . . . . . . . . . . . . . . . . . . . .                    .   .   .   .   .   .   .   .   .   .   30
  5.5   Data Sharing . . . . . . . . . . . . . . . . . . . . . . .                     .   .   .   .   .   .   .   .   .   .   31
  5.6   Walkthrough Example . . . . . . . . . . . . . . . . . .                        .   .   .   .   .   .   .   .   .   .   33

6 Supernova-Neutrino Sensitive         Detectors                                                                               34
  6.1 Water Cherenkov Detectors        . . . . . . .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   34
      6.1.1 Super-Kamiokande .         . . . . . . .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   34
      6.1.2 Hyper-Kamiokande .         . . . . . . .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   35
      6.1.3 IceCube . . . . . . .      . . . . . . .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   36
      6.1.4 KM3NeT . . . . . .         . . . . . . .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   37
  6.2 Scintillator Detectors . . . .   . . . . . . .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   38
      6.2.1 Baksan . . . . . . . .     . . . . . . .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   38
      6.2.2 LVD . . . . . . . . .      . . . . . . .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   38
      6.2.3 Borexino . . . . . . .     . . . . . . .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   39
CONTENTS                                                                                                 6

         6.2.4 KamLAND . . . . . . . . . . . . . . . . . . . . . . .        .   .   .   .   .   .   .   40
         6.2.5 JUNO . . . . . . . . . . . . . . . . . . . . . . . . . .     .   .   .   .   .   .   .   41
         6.2.6 SNO+ . . . . . . . . . . . . . . . . . . . . . . . . . .     .   .   .   .   .   .   .   42
         6.2.7 NOvA . . . . . . . . . . . . . . . . . . . . . . . . . .     .   .   .   .   .   .   .   43
   6.3   Lead-Based Detectors . . . . . . . . . . . . . . . . . . . . . .   .   .   .   .   .   .   .   44
         6.3.1 HALO . . . . . . . . . . . . . . . . . . . . . . . . . .     .   .   .   .   .   .   .   44
         6.3.2 Future Lead-Based Detectors . . . . . . . . . . . . .        .   .   .   .   .   .   .   45
   6.4   Liquid Noble Dark Matter Detectors . . . . . . . . . . . . .       .   .   .   .   .   .   .   45
         6.4.1 Global Argon Dark Matter Collaboration . . . . . . .         .   .   .   .   .   .   .   46
         6.4.2 Xenon . . . . . . . . . . . . . . . . . . . . . . . . . .    .   .   .   .   .   .   .   46
   6.5   Liquid Argon Time Projection Chamber Neutrino Detectors            .   .   .   .   .   .   .   47
         6.5.1 DUNE . . . . . . . . . . . . . . . . . . . . . . . . . .     .   .   .   .   .   .   .   47
   6.6   Detection in Other Low-background Detectors . . . . . . . .        .   .   .   .   .   .   .   48
         6.6.1 The nEXO Experiment . . . . . . . . . . . . . . . . .        .   .   .   .   .   .   .   48
   6.7   Detection Estimates . . . . . . . . . . . . . . . . . . . . . .    .   .   .   .   .   .   .   49

7 Amateur Astronomer Engagement                                                     49
  7.1 Thrusts of Amateur Astronomer Engagement . . . . . . . . . . . . . . . 50
  7.2 Assessing Observational Readiness . . . . . . . . . . . . . . . . . . . . . . 51

8 Summary                                                                                               52

9 Acknowledgements                                                                                      52
CONTENTS                                                                                 7

1. Introduction

The explosion of a star within the Milky Way Galaxy will provide us with a front row
seat of physics under conditions that could never be produced in a terrestrial experiment.
While the remnants of the explosion will be observable for many thousands of years, the
information about what occurred in the core of the star to cause the explosion will be
most easily found in the first tens of seconds. It is therefore imperative that we be able
to detect the supernova as soon as it begins — if not sooner.
     The familiar blast of light visible across the universe is not really when a supernova
begins: that electromagnetic radiation starts when the shockwave from the stellar core’s
collapse reaches the surface and breaks out. Neutrinos are produced at the start of the
core collapse process, and escape a supernova explosion well before the photon emission is
visible and thus provide the earliest opportunity to anticipate the imminent appearance
of a galactic supernova in time to alert observatories. The Supernova Early Warning
System (SNEWS) is an open, public alert system that has provided the capability for
such an early warning since 2005 by combining the detection capabilities of a variety
of neutrino detectors worldwide (Antonioli et al., 2004). If several detectors report
a potential supernova within a small time window, SNEWS will issue an alert to its
subscribers which include astronomical observatories, neutrino detectors, and amateur
astronomers and citizen scientists. SNEWS is one of the few successful examples of
cyberinfrastructure spanning major neutrino experiments.
     Since the SNEWS network was first established over a decade ago, the particle
astrophysics landscape has evolved considerably. The detection of gravitational waves by
LIGO/Virgo along with electromagnetic observations of a neutron star merger (Abbott
et al., 2017a), and the subsequent possible observation of neutrinos from an active
blazar by IceCube (Aartsen et al., 2018) have ushered in a new era of multi-messenger
astrophysics. At the same time, neutrino detector technologies and data analysis
techniques have progressed in recent years, and the ability of detectors to detect and
analyze neutrinos from galactic supernovae in real time has improved substantially.
     In its current form, SNEWS is designed to send a prompt alert based on a simple
coincidence, but its functionality can be extended to take advantage of these recent
advances. The overarching aim is to enhance the overall science obtained from the next
galactic core-collapse supernova (CCSN). These are expected to occur rarely enough
(1.63 ± 0.46/century in the Milky Way (Rozwadowska et al., 2021)) that we need to
extract all the information possible from the next such CCSN to happen. Specifically,
the goals of SNEWS 2.0 are to:
  • reduce the threshold for generating alerts in order to gain sensitivity;
  • reduce alert latency;
  • combine pointing information from individual experiments and enhance it via
    timing triangulation;
  • implement a pre-supernova alert based on the rising neutrino flux which precedes
CONTENTS                                                                               8

    core-collapse;
  • develop a follow-up observing strategy to prepare the astronomical community for
    the next galactic supernova; and
  • engage amateur astronomer and citizen science communities through alert
    dissemination and outreach.
     In this paper, we describe enhancements to SNEWS to exploit new opportunities
in the era of multi-messenger astrophysics as a means of realizing these goals. Section 1
introduces the existing network and the overall plan. Section 2 provides background on
the types of transient events that are of interest to SNEWS and the characteristics
of their signals at Earth. Section 3 describes how pointing information can be
extracted from a neutrino signal via anisotropic interactions and signal triangulation
between multiple detectors. Section 4 explores the possibility of producing an earlier
warning by measuring the pre-supernova neutrino flux from stars during silicon burning,
which directly precedes core-collapse. Section 5 is a review of the SNEWS design,
how SNEWS 2.0 alerts will be disseminated, and how follow-up observations can be
incorporated, while Section 6 describes the experiments involved. Finally, Section 7
discusses public outreach, including how SNEWS 2.0 will interface with amateur
astronomers and citizen scientists.

1.1. Current Configuration (SNEWS 1.0)
The SNEWS 1.0 system was designed to give a supernova neutrino alert that was
  • Prompt, providing an alert within minutes and followup within hours;
  • Positive, with less than one false alarm per century; and
  • Pointing, providing a sky location if and when possible by passing along
    experiments’ estimates.
The design was primarily driven by the Positive requirement. Only extremely high
quality coincidences could automatically trigger an alert. All other coincidences would
require human intervention before an alert could be triggered.
     SNEWS (https://snews.bnl.gov) currently involves an international collabora-
tion of supernova neutrino detectors: Super-Kamiokande, LVD, IceCube, Borexino,
KamLAND, HALO, and Daya Bay (with NOvA, KM3NET, and Baksan testing their
connections to join soon). SNEWS has been operational since 1998 and has been run-
ning in a fully-automated mode since 2005 with near-100% up-time. The main idea
of SNEWS is to provide prompt, high-confidence alerts of nearby CCSN by requiring a
burst coincidence between detectors; this allows alarms from individual detectors to go
out promptly without needing a human check.
     SNEWS operates two “coincidence servers”: a primary server at Brookhaven
National Laboratory and a backup at the University of Bologna. The participating
experiments each run their own online supernova monitors, and run client code provided
by SNEWS to send datagrams to the servers if supernova-like bursts are observed. The
CONTENTS                                                                                  9

minimal information provided in the client datagram is the experiment, time of the first
event of the burst, and a burst-quality parameter. Experiments may choose to provide
directional and burst size information if it is available promptly. “Gold” alerts, for which
input datagrams must satisfy several quality criteria, are sent out automatically by the
server to a mailing list if a coincidence within 10 s is found. Email alerts are provided
first to the other experiments and to “express-line” subscribers (LIGO, ANTARES, and
the Gamma-Ray Coordination Network (GCN)(Barthelmy et al., 1995)), and they are
also available by direct socket connection (NOvA and XENON1T). “Silver” alerts are
sent to the experiments only.
      The current SNEWS requirement for accidental-coincidence alerts is that they
must occur less than once per century. If individual input alarm rates become too
high, or there are other low-quality indicators in the input, coincidence output is
demoted to “silver”. Single experiments may also send datagrams to SNEWS with
sufficiently well-vetted alarms to be propagated automatically as “individual” alerts.
Detailed information on coincidence criteria for the configuration of SNEWS can be
found in (Antonioli et al., 2004; Scholberg, 2008). A weekly test alert is sent via GCN
every Tuesday at noon Eastern.

1.2. SNEWS 2.0
SNEWS 2.0 will be an upgrade of the SNEWS system for the age of multi-messenger
astronomy. In this environment, false alarms are acceptable, low probability events
should be reported, and SNEWS will be one of many multi-messenger alert systems.
Nevertheless, SNEWS remains unique by combining data from different neutrino
observatories, and providing clear summaries of neutrino data for astronomers. This
section outlines a number of improvements which will be made using more detectors
and different techniques than the original SNEWS.
     Recent additions to the suite of potential detectors for the next Galactic supernova
are large dark matter detectors. Capitalizing on the recently-discovered coherent elastic
neutrino-nucleus scattering (CEνNS) (Akimov et al., 2017), those detectors rely on
the coherent enhancement of the neutrino cross-section for supernova burst detection
that can be probed thanks to low (< keV) energy thresholds. Those detectors provide
a flavor-insensitive detection channel and thus a total-flux measurement of the total
energy going into neutrinos, independent of, for example, uncertainties from neutrino
oscillations (Lang et al., 2016; Chakraborty et al., 2014). Furthermore, the combination
of CEνNS and inelastic interaction channels will help to disentangle oscillation effects.
Currently, these detectors need improved understanding of backgrounds at the lowest
energies that are relevant here (Aprile et al., 2014; Sorensen, 2017; Sorensen and Kamdin,
2018). So far, XENON1T has a dedicated trigger following a SNEWS alert. SNEWS 2.0
aims to go a step further to enable the integration of dark-matter detectors as inputs to
SNEWS.
     SNEWS 2.0 also intends to develop and provide a true pre-supernova alert—a
CONTENTS                                                                                10

pre-core-collapse alert. This is based on the predicted uptick in neutrino production
that accompanies the final burning stages of a doomed star (Odrzywolek et al.,
2004a; Kato et al., 2017; Patton et al., 2017). This alert has been implemented in
KamLAND (Asakura et al., 2016), and provides a 3 σ detection 48 hours prior to the
explosion of a 25 M star at 690 pc. Extending this alert to the network should expand
the sensitivity to a larger fraction of the galaxy.
     SNEWS 2.0 could also be used to communicate and organize planned shutdowns
or downtime in each detector to ensure that the overall supernova detection livetime is
not affected, since other participating experiments can cover.

1.3. Lowering the Threshold
One general benefit of combining detectors’ data real-time would be the lowering of the
effective threshold for observing a signal. An astrophysical neutrino signal would be
observable in many detectors at once, but might be not strong enough to be significant
in any one detector. Since most detectors are large enough to be sensitive to a CCSN
somewhere in our galaxy, the most obvious benefit of being able to see farther against
the inverse-square law flux suppression with distance is not as useful as it might seem
at first, with the exception of a supernova in the Magellanic clouds, where the distance
is substantial and the flux is borderline for most detectors. However, a combination
of detector signals allowing for more sensitivity to low flux would enhance the world’s
ability to notice any unusually low flux events. Four important examples of this are the
pre-supernova neutrinos introduced in the previous section (elaborated on in Sec. 4), the
neutrinos emitted at the latest times in the burst (Weishi Li et al., 2020), and neutrinos
from three other types of potential transient bursts described in Sec. 2.4. As the current
range of such detection is only hundreds of parsecs, increasing sensitivity via comparing
sub-threshold signals in different experiments will increase the number of progenitors
under observation by a factor of distance cubed.
      These pre-supernova neutrinos have a lower energy (∼ a few MeV) and a much
lower flux than core-collapse supernova neutrinos. Being low energy, there are also
more potential background events to confuse with a potential signal. In the near
future, lowering the energy thresholds and better background rejection are in the plans
for both running and planned detectors, expanding the experiments able to do so
(currently, only KAMLAND has this capability). For example, delayed coincidence
with gadolinium is being implemented Super-Kamiokande, which effectively lowers the
energy threshold by the confirmation of inverse-β decay events (Beacom and Vagins,
2004). At liquid-scintillator detectors, the low energy threshold (1.8 MeV) is achieved via
good scintillation light production. The detection of keV-neutrinos will be practical via
CEνNS at large dark-matter detectors. A supernova alert with pre-supernova neutrinos
has been investigated in several recent works (Kato et al., 2017; Asakura et al., 2016;
Raj et al., 2020; Simpson et al., 2019; Li et al., 2020).
      In addition to the developments in individual neutrino detectors, their combination
CONTENTS                                                                                11

via SNEWS 2.0 reduces the uncertainty of the supernova alert and effectively lowers the
threshold for the alert issue. Because the energy of pre-supernova neutrinos increases as
the pre-supernova stellar core evolves, an earlier alert is possible by combining different
experiments’ sub-threshold data to provide plenty of preparation time for the detection
of other observables. This earlier alert will maximize the information to be gained
from multi-messenger astronomy, yielding information of supernovae from a different
perspective. Moreover, pre-supernova neutrinos will be one of the useful tools to prove
the theory of stellar evolution (Kato et al., 2015; Yoshida et al., 2016) and long term
detection over several stellar evolutionary phases and experiments will improve the
results.

2. Stellar Core Collapse Signals

The dominant source of supernova neutrino bursts are Core Collapse SuperNovae
(CCSN). These type of supernova occur when a massive (more than ∼8–10 M ) star,
after successively burning elements from hydrogen to silicon, forms an inert, but growing,
iron core. This core soon reaches the effective Chandrasekhar mass and collapses due
to the unmatchable strength of gravity. The collapse continues until the density reaches
nuclear densities where the equation of state stiffens and the nuclear force is able to
stabilize the core against gravity. The formation of the protoneutron star also leads to
the formation of a shock wave. It is this shock wave that for successful supernovae
will traverse the star over the course of minutes to hours and unbind all but the
innermost material. Core-collapse supernovae emit signals in three cosmic messengers—
electromagnetic emission, neutrinos, and gravitational waves—and also cosmic rays at
later times. Figure 1 shows an example of the expected time sequence for these signals,
starting from before (left panel) and after (right panel) core bounce.
     While each cosmic messenger is valuable by itself, when analyzed together, they
provide a comprehensive understanding that is impossible to achieve from any single one
of them alone. Multi-messenger astrophysics had two foundational discoveries in 2017:
a binary neutron star merger that produced gravitational waves and electromagnetic
radiation (Abbott et al., 2017a), and the coincident detection of high-energy neutrinos
and electromagnetic emission from a blazar (Aartsen et al., 2018). We expect multi-
messenger observations of the next Galactic core-collapse supernova offer similar
synthesis opportunities. Coordinating timely follow-up observations that enable a
true multi-messenger analysis demands rapid identification and characterization of the
neutrino signal, along with prompt broadcasting to ensure that transient emission only
detectable on short time scales is recorded. Enabling this coordination is the SNEWS2.0
raison-d’être.
CONTENTS                                                                                                 12

                                             54
                                                                                    νe
                                             52                                     νe
                                                                                    νx

                Log (luminosity [erg s-1])

                                                                       ALERT	
  
                                             50                                    GW
                                                                                   EM
                                             48   pre-SN νe
                                             46                                      SBO
                                             44
                                                                                         plateau
                                             42
                                             40       progenitor
                                             38
                                                  9     6 3 0 -2            0   2    4    6   8
                                                        Log (time relative to bounce [s])

              Figure 1. Time sequence for multi-messenger signals pre- (left panel) and post-
              (right panel) core collapse of a non-rotating 17 M progenitor star. Neutrinos (νe ,
              ν̄e , and νx are shown by red, thick red, and magenta lines, respectively, where νx
              represents heavy-lepton neutrinos: νµ , ντ , ν̄µ , and ν̄τ ), gravitational waves (blue line),
              and electromagnetic signals (black line) are shown. Solid lines are predictions from
              a hydro-dynamical simulation with axis-symmetric radiation, while dashed lines are
              approximate predictions. Neutrino emission prior to collapse arises from the last
              moments of stellar evolution, but is quickly overtaken during collapse by the neutrino
              burst. The electromagnetic signal exhibits the shock breakout (SBO), plateau, and
              decay components. Note that the height of the curves does not reflect the energy
              output in each messenger; the total energy emitted after the bounce in the form of ν̄e ,
              photons, and gravitational waves are ∼6 × 1052 erg, ∼4 × 1049 erg, and ∼7 × 1046 erg,
              respectively. The focus of SNEWS 2.0 is to establish the neutrino burst as an alert for
              gravitational waves and electromagnetic followup, as shown by arrows. Adapted from
              (Nakamura et al., 2016).

2.1. Neutrinos from Core Collapse Supernovae
The neutrino emission from a core collapse supernova in our Galaxy cannot be hidden in
any way. The neutrinos are not obscured by dust as electromagnetic signals may be, nor
would failure of the explosion mean the supernova would evade our detection: a large
burst of neutrinos would still be emitted prior the formation of a black hole. Finally,
the present detection horizon for neutrinos reaches out beyond the edge of the Milky
Way. For all these reasons, neutrinos are a unique messenger to provide a compelling
trigger for an alert. Coupled with gravitational waves (whose detection will also be
enhanced by the precise timing information provided by neutrinos) and electromagnetic
observations, the neutrinos will allow us to piece together a comprehensive picture of the
supernova from the moment of core collapse to supernova shock breakout and beyond.
Expected features in the neutrino signal will permit us to probe a long list of topics,
including: key aspects of the supernova explosion mechanism (e.g., fluid instabilities vs.
CONTENTS                                                                               13

core rotation), the nuclear equation of state, the stellar radius and interior structure,
explosive nucleosynthesis, the nature of the remnant core (neutron star vs. black hole),
as well as answer questions about the fundamental properties of neutrinos, and even test
Beyond-Standard-Model physics (Horiuchi and Kneller, 2018). To fully develop these
prospects, it is essential the supernova be detected to the latest times possible with
good flavor and spectrum information (Weishi Li et al., 2020; Nakazato and Suzuki,
2020). The multi-messenger nature of the supernova signal greatly helps in extracting
this information from the neutrinos. For example, it has been shown by (Warren et al.,
2020) that the neutrino emission is correlated with the gravitational wave signal which
would aid in disentangling the neutrino oscillation effects.

2.2. Gravitational Wave Signals from Core Collapse
Together with neutrinos, gravitational waves provide a unique probe of the core collapse
in realtime. The emission of gravitational waves is strongly dependent on the asymmetry
of the collapsing core and the nuclear equation of state, opening a view of the collapsing
core complementary to neutrinos (Janka, 2017; Kotake, 2013; Morozova et al., 2018). By
combining gravitational waves with neutrinos and electromagnetic waves, key aspects of
the collapse, from the spin of the collapsed core to the supernova explosion mechanism
and black hole formation, become be more robustly probed. At present, even for
a conservative prediction of the emitted gravitational wave signal, detectors such as
Advanced LIGO, Advanced Virgo, and KAGRA are able to detect CCSN gravitational
waves out to a few kpc from the Earth, while future detectors such as the Einstein
Telescope can reach the entire Milky Way. The detection horizon of the circular
polarization can be significantly larger than the gravitational wave amplitude, and can
also help reveal inner dynamics (Hayama et al., 2018).
     Several mechanisms can generate gravitational waves during a CCSN, for a recent
review see (Abdikamalov et al., 2020) and references therein. The majority of these
signals have the common feature of being short and “burst like”, i.e. impulsive signals
lasting less than a second and very difficult to model. These characteristics make
detection more challenging. The identification of a temporal window in which to look
for the signal significantly increases the detection efficiency. Neutrinos can provide the
best temporal trigger for this gravitational wave search; indeed the neutrino signal for
a Galactic CCSN allows the time of the core “bounce” to be identified within a window
of ∼ 10 ms or less (Pagliaroli et al., 2009a; Halzen and Raffelt, 2009). The use of this
information improves the background reduction of gravitational wave detectors, with
consequent increases of the detection capability (Nakamura et al., 2016). In the case
of long-lasting GW emission due to neutron star oscillations (Radice et al., 2019) the
identification of the time of the bounce through neutrinos could provide a reference
point to start the search.
CONTENTS                                                                              14

2.3. Electromagnetic Signals
EM radiation in the first hours to days after core collapse explosion provides critical
information about the progenitor star and the overall energy budget and dynamics of the
core collapse explosion. A few hours to days after the core collapse, the supernova shock
breaks out of the progenitor surface, suddenly releasing the photons behind the shock
in a flash bright in UV and X-rays, known as shock breakout (SBO) emission. SBO has
been observed on rare occasions in extragalactic systems (Soderberg et al., 2008; Gezari
et al., 2010; Bersten et al., 2018). The SBO signal provides important information about
the supernova, such as the radius, mass, and structure of the progenitor star, and the
kinematic energy associated with the rapidly expanding ejecta. Initial observations of
the gamma flux from the first moments of a SN will be also be important to help
constrain the terrestrial effects of gamma rays from historical SNe on atmospheric
chemistry and climate science (Brakenridge, 2020; Jull et al., 2018). Knowledge of where
and when to anticipate the signal will ensure that the peak luminosity and duration of
the SBO (strongest at UV and soft X-ray wavelengths) is not lost. Even including
SN1987A, the precise time between onset of core collapse and shock break out has never
been measured (Arnett et al., 1989; Ensman and Burrows, 1992). Prompt alert and
coordinated follow up with SNEWS 2.0 will make this possible.

2.4. Other Transients
While core collapse supernovae are expected to be the dominant type of supernovae in
the Milky Way, they are not the only astrophysical sources of neutrino bursts. Bursts are
also expected from Type Ia supernovae (SNIae), pair-instability supernovae (PISNe),
compact object mergers, and possibly others yet unknown. There are a number of
questions, many fundamental, about these other neutrino transients so that a neutrino
signal from any one of them would represent as rich an opportunity to advance our
knowledge as the signal from a core collapse.

2.4.1. Type Ia Supernovae The progenitor systems of Type Ia supernovae and their
associated explosion mechanisms remain debated. The possible progenitors of SNIae —
and the observational constraints upon the various scenarios — are discussed extensively
in Maoz, Mannucci and Nelemans (Maoz et al., 2014) and Ruiz-Lapuente (Ruiz-
Lapuente, 2014). Even if we accept the canonical model of a SNIa as the disruption of
a Chandrasekhar mass (1.4 M ) carbon-oxygen white dwarf, many different scenarios
for how the explosion proceeds can be found in the literature (Khokhlov, 1991; Plewa
et al., 2004; Bravo and Garcı́a-Senz, 2006). We refer the reader to Hillebrandt et al.
(Hillebrandt et al., 2013) for a review.
     The neutrino emission from a limited number of SNIa simulations has been
computed (Odrzywolek and Plewa, 2011; Wright et al., 2016, 2017). Wright et al.
considered the most optimistic case (known as the DDT) and a more general case (their
GCD case). The number of events they expect in a 374 kt water-Cherenkov detector
CONTENTS                                                                                 15

from a SNIa at a distance of 10 kpc is of order 1 for the DDT case and 0.01 for the less
optimistic GCD. A SNIa would have to be within a few kpc in order to detect tens of
events but the probability the next Galactic supernova is within 5 kpc is only of order
10% according to Adams et al. (Adams et al., 2013).

2.4.2. Pair Instability Supernovae Very massive stars can explode as a PISN if they
form a carbon-oxygen core in the range of 64 M < MCO < 133 M (Heger and Woosley,
2002). The temperatures in these cores are sufficiently high and the electron degeneracy
sufficiently low that electron-positron pairs are created. The formation of the pairs
softens the equation of state causing a contraction of the core triggering explosive
burning of the oxygen (Barkat et al., 1967; Rakavy and Shaviv, 1967; Fraley, 1968).
The energy released is enough to unbind the entire star leaving behind no remnant.
Some models of PISNe produce very large amounts of 56 Ni and PISN are candidates for
some superluminous supernovae (Smith et al., 2007; Gal-Yam et al., 2009; Cooke et al.,
2012; Lunnan et al., 2016).
     The long-standing expectation of theorists is that only metal-free stars could remain
sufficiently massive to explode as PISN (Heger and Woosley, 2002). However this
expectation has been challenged in recent years. (Langer et al., 2007) found PISN can
occur in stars with metallicities as large as Z /3 while (Georgy et al., 2017) obtained
the conditions for a PISN at near solar metallicities if they included surface magnetic
fields. Thus, theoretically at least, a PISN in the Milky Way or one of its satellites
cannot be ruled out.
     The rate of PISNe is uncertain because a) observationally we lack an unambiguous
method for discriminating these kind of supernovae from the others and b) theorists
have not reached a consensus on which masses at a given metallicity produce these
kinds of events. The estimate by Langer et al. is for a rate of 10−4 yr−1 but that could
be larger by an order of magnitude if the recent revisions to the progenitor are correct.
     The neutrino signals from two PISN simulations have been computed by (Wright
et al., 2017). The two models they considered were a low-mass and high-mass case so
that the computed signals spanned the range of possibilities. The flux at Earth from
a ‘small’ PISN at 10 kpc was similar to the most optimistic SNIa case i.e. around 1–2
events, but for a ‘large’ PISN at the same distance the flux was much larger, between
50–100 events depending upon the equation of state and the neutrino mass ordering.

2.4.3. Compact Object Mergers The neutrino emission from merging neutron stars has
been computed by (Rosswog and Liebendörfer, 2003) and the neutrino emission from a
black-hole-neutron star merger simulation was computed by (Caballero et al., 2009). In
both cases the neutrino emission is similar to that from a core-collapse supernova i.e., the
neutrino luminosities and mean energies are within a factor of a a few of those found in
core-collapse simulations), and therefore give similar event rates in detectors. However
there are differences: in a core-collapse there are more neutrinos than antineutrinos
emitted and the duration of the burst is of order 10 s. In a neutron star merger the
CONTENTS                                                                                  16

opposite matter-anti-matter ratio is expected and the signal lasts for 1 second unless
the supermassive neutron star can be prevented from forming a black hole. The rate
of black-hole - neutron star mergers in not known precisely but the rate of neutron
star-neutron star mergers can be better estimated because there exist a number of such
systems in the Milky Way. (Abadie et al., 2010) calculate the likely event rate to be
10−4 yr−1 while (Kalogera et al., 2004a,b) give the plausible range to be from 10−6 yr−1
to 10−3 yr−1 .
     To detect neutrinos from black-hole-neutron star and neutron star-neutron star
mergers is challenging. A promising strategy is to search for neutrinos in time-
coincidence with detections of mergers in gravitational waves, using a time window
of, e.g., 1 s after each merger (Kyutoku and Kashiyama, 2018; Lin and Lunardini,
2020). This strategy reduces the backgrounds very effectively so that, in fortunate
circumstances, even the detection of a single, time-coincident neutrino can be
statistically significant. If alerts from Advanced LIGO (Abadie et al., 2010) (which
has a distance of sensitivity to mergers of about 200 Mpc) were used, a megaton water
Cherenkov detector could record about 1 neutrino detection per century (Kyutoku and
Kashiyama, 2018). When operating in synergy with third-generation gravitational wave
observatories, like the proposed Einstein Telescope (Punturo et al., 2010), and Cosmic
Explorer (Abbott et al., 2017b) (sensitivity up to redshift z ∼ 2), the same detector
could identify of up to a few neutrinos from mergers per decade, and start placing
constraints on the parameters space already after a decade or so of operation (Lin and
Lunardini, 2020).

3. Pointing to the Supernova with Neutrinos

While supernovae are optically highly luminous, a large fraction are anticipated to be
heavily attenuated by dust along the line of sight, typically within the disk of the Milky
Way. The supernova optical signal as observed at Earth has been estimated adopting a
standard intrinsic supernova luminosity distribution, Galactic distribution for supernova
occurrence, and a simple model of Galactic dust extinction. According to (Nakamura
et al., 2016; Adams et al., 2013), the dominant fraction (some 50%) will be observable
with 1–2 meter class telescopes in the optical band. An additional 10% of supernova
can be observed by larger 4–8 meter class telescopes, while the faintest 25% will be as
faint as 25–26 magnitudes and require dedicated observations by the largest available
telescopes (however, dust attenuation uncertainties are large in this regime, being at
least a few magnitudes). This is illustrated in the top panel of Figure 2. The fields
of view of the relevant telescopes typically do not cover more than several degrees in
a single pointing, as shown by the rectangular boxes on the bottom panel of Figure 2.
This highlights the quantitative demands for a combined rapid (in time) and accurate
(in direction pointing) alert.
     In order to be an effective trigger, the neutrino alert needs to be sent faster than the
delay between neutrinos and the first electromagnetic signal, and also alert the pointing
CONTENTS                                                                                                                                   17
                                               20
                                                          9.8%       15.6%       15.5%   13.5%    11.5%           9.4%        7.4%

                             Probability (%)
                                               15
                                               10
                                                5
                                                      1.2%                                                                    16.1%
                                                0
                                                                             Naked eye
                                                                                            Evryscope      1-2m 4m >8m
                                               100

                       FOV diameter (deg)
                                               10                                                ASAS-SN
                                                                                                            ZTF

                                                                                                        Pan-STARRS     LSST

                                                                                                                  Blanco
                                                                                                                           Subaru

                                                1                                                                  CFHT

                                               0.1
                                                     -5          0           5        10     15             20             25         30
                                                                                  Optical magnitude

              Figure 2. Optical follow-up requirements for the next Galactic supernova. The top
              panel shows a histogram of the apparent magnitude probability distribution for the
              shock breakout signal of a Galactic supernova (uncertainties not shown; see text). The
              bottom panel shows the magnitude sensitivity range and fields of view (FOV) of optical
              telescopes: ASAS-SN, Blanco, CFHT, Evryscope, LSST, Pan-STARRS, Subaru, and
              ZTF. When the optical magnitude is brighter than ∼ 15 mag, early detection is feasible
              thanks to the wide FOV of small-aperture telescopes. However, fainter cases are
              more challenging since there are no > 1 m telescope with a FOV larger than ∼ 6◦
              diameter. Therefore, this sets the target accuracy for triangulation by SNEWS 2.0.
              Note that for the brightest supernovae, telescopes will need high-quality filters for
              accurate photometry, shown by the fading color. Adapted from (Nakamura et al.,
              2016).

in the sky to within a few degrees. The SNEWS 2.0 alert will be automated and sent
electronically to meet the timing demands (see Section 5). Triangulation-based pointing
will be implemented in SNEWS 2.0. Such pointing techniques were originally explored
in (Beacom and Vogel, 1999) and further developed by (Mühlbeier et al., 2013; Brdar
et al., 2018; Linzer and Scholberg, 2019).
     An important consideration is that an alert may be associated with an event that
is challenging to observe at EM wavelengths. Possible scenarios include a distant
supernova associated with large extinction due to dust and/or formation of a black
hole with a weak explosion. In these cases follow up strategies informed by the neutrino
alert are needed. For example, in extreme scenarios the Vera Rubin Observatory’s field
of view and photometric depth are uniquely suited (Walter et al., 2019).

3.1. Anisotropic Interactions
3.1.1. Water Cherenkov Large water Cherenkov detectors (WCDs), such as Super-K
and the future Hyper-K, have potentially good supernova pointing capability through
the anisotropic neutrino-electron elastic scattering (ES) interaction (Beacom and Vogel,
CONTENTS                                                                                           18

1999; Tomas et al., 2003). The majority of supernova neutrinos interact in WCDs
through inverse beta decay (IBD). In IBD the direction of the outgoing positron is
nearly random with respect to the direction of the incoming supernova neutrino. The
angular distribution of measured positron directions from a large number of IBD events
is necessary to detect (and possibly use) the small anisotropy in the direction of the
neutrino flux.
     Fortunately, a few percent of the interactions are due to elastic scatter of the
supernova neutrinos from electrons. The outgoing electrons are preferentially forward-
scattered and the reconstructed direction of the scattered electron is correlated to the
direction of the incoming neutrino. The angular distribution of measured electron
directions from ES shows a strong correlation with the direction of the neutrino flux.
The magnitude of the anisotropy varies with neutrino energy and flavor. However, since
the ES cross sections are small, a large detector mass is required to measure enough
ES interactions for direction finding. At present WCDs cannot accurately differentiate
between IBD and ES interactions, although adding increasing amounts of Gadolinium
to Super-K will improve this. Thus, the high ratio of IBD to ES interactions reduces
the signal-to-noise ratio of the direction signal, as shown in Figure 3.

              Figure 3. Reconstructed skymap of MC simulated Super-K supernova event direction
              vectors. Red points: anisotropic elastic scatter events, blue points: IBD and other
              nearly isotropic events, star point: direction vector from supernova (Abe et al., 2016).

    In a large WCD the direction of the neutrino flux (and therefore the supernova
CONTENTS                                                                                 19

direction) may be found by analyzing the energies and direction vectors of the
electrons/positrons, reconstructed from the Cherenkov ring for each event. For example,
the current Super-K real-time supernova burst monitor performs SN direction finding
using a maximum likelihood method (Abe et al., 2016). The direction and energy
of each event is used to calculate the likelihood of a given supernova direction and
event reaction channel based on a probability density function (PDF) determined from
supernova neutrino flux models and Super-K MC simulations. The supernova direction
angles, and other parameters, are varied until the total likelihood is maximized. The
pointing accuracy for a supernova at 10 kpc is estimated using a modern supernova
model to be 4.3–5.9◦ (68.2% C.L.) covering all combinations of neutrino oscillations and
mass orderings. This will improve to 3.3–4.1◦ for the fully doped SK-Gd.
     The accuracy of supernova direction finding based on the anisotropy of ES events
depends on the number of events. This varies with detector volume, supernova distance,
neutrino oscillations and supernova mass and neutrino emission mechanisms. The
coming Mtonne-scale water Cherenkov detectors, such as Hyper-Kamiokande, will also
have improved direction finding due to increased statistics in the ES channel. The
pointing accuracy for Hyper-K is expected to be 1–1.3◦ (Abe et al., 2018).
     The introduction of small amounts of gadolinium into the water volume of WCDs
will allow accurate tagging of individual IBD events. Thus, direction finding routines
could de-weight or exclude the IBD events, potentially increasing the speed and pointing
accuracy. Super-K will be implementing such an upgrade in the near future.

3.1.2. Liquid Argon Liquid argon time projection chambers (LArTPC) detectors have
the ability to do fine-grained tracking of the final-state particles, and, like water
Cherenkov detectors, can exploit the intrinsic directionality of anisotropic interactions in
the detector. Ability to tag different interaction channels also helps. Elastic scattering
interactions on electrons with well-known energy dependence can be used, as well as the
charged-current νe absorption interactions on 40 Ar. The latter have a relatively weak
anisotropy but large statistics.
      Unlike water Cherenkov signals, the tracked electrons have a head-tail ambiguity
that results in about half of them with a fake reconstructed backwards direction. This
ambiguity can be resolved statistically using sophisticated reconstruction techniques.
Improvement by using the directionality of bremsstrahlung gammas, which are emitted
preferentially in the electron travel direction, has been demonstrated in DUNE. Using a
likelihood technique with the ensemble of electron scattering and νe CC events, DUNE
has demonstrated about 5◦ pointing for a 10 kpc supernova signal (Abi et al., 2020a).

3.1.3. Liquid Scintillator Liquid scintillator and water Cherenkov detectors alike are
mostly sensitive to IBD interactions – the major difference between the two being that
such interactions are considered a background for supernova pointing in the latter while
they are considered a signal in the former. Indeed, while considered isotropic at first
order, the positron and neutron emitted after an IBD interaction both possess a slight
CONTENTS                                                                                   20

energy-dependent anisotropy (Strumia and Vissani, 2003). At the energies of interest
for supernova neutrino detection, the positron quickly deposits its energy and therefore
most of the anisotropy is carried away by the neutron, always emitted in the forward
direction. Although this appears to be similar to the forward emission of an electron in
elastic scattering, detecting a neutron direction is arduous in large scintillator detectors.
In the vast majority of cases, only the position of the neutron capture vertex after
thermalization and diffusion can be determined. For each IBD interaction, a direction
vector, starting at the reconstructed positron vertex and ending at the reconstructed
neutron capture vertex, can be defined. Due to the smearing caused by the neutron
transport after its creation, a single IBD vector is not sufficient to efficiently reconstruct
its neutrino incoming direction. However,the analysis of thousands or more of IBD
interactions can help reconstruct the statistical direction of an incoming neutrino flux,
as demonstrated by the CHOOZ collaboration with about 2,700 events (Apollonio et al.,
2000).
      Such an analysis can be performed to determine the expected direction of a
supernova-induced neutrino flux, as shown in (Fischer et al., 2015). In this study, the
statistical nature of the supernova direction reconstruction through IBD anisotropy was
exploited by combining the direction vectors of all IBD interactions from several liquid
scintillator-based detectors, existing or proposed. While the pointing capabilities of
individual existing detectors, shown in Figure 4, are no match for the accuracy of Super-
K, their combination, as well as the introduction of JUNO in a near future, provides
non-negligible pointing information. With the addition of JUNO to the existing large
liquid scintillator detectors, supernova pointing accuracy through IBD interactions could
reach 12 degrees (68% C.L) for a supernova located 10 kpc away.
      It is worth noting that efforts are underway to extract directional information
from the small amount of Cherenkov light which leads the largely isotropic scintillation
light (Aberle et al., 2014). The CHESS experiment found that a time resolution of
338 ± 12 ps (FWHM) was required for reasonable efficiency in separating the two
light components in a mixture of LAB with 2g/L PPO (Caravaca et al., 2017). An
alternative to fast PMT’s is to slow the emission of scintillation light (e.g., (Wang and
Chen, 2020; Biller et al., 2020)), which is possible with different scintillators and fluors.
Finally, the different spectra of the two components may be exploited using dichroic
filters (Kaptanoglu et al., 2019, 2020). Such ideas may be exploited in future liquid
scintillator detectors, or upgrades of current ones.

3.2. Triangulation
In triangulation, the time delays of neutrino events observed between detectors at
different geographical locations are used to infer the direction of the supernova. For
a pair of detectors i and j, the delay between them ∆tij is defined as follows:
               ∆tij = d~ij · ~n/c,                                                        (1)
where d~ij is the vector connecting two detector sites and ~n is the unit vector defining
CONTENTS                                                                                                                  21

                                                                                                   KamLAND (1kt)
                                                                                                   SNO+ (0.8kt)
                               80                                                                  Borexino (0.3kt)
                                                                                                   DayaBay (0.3kt)
                                                                                                   DoubleChooz (0.05kt)
                               70                                                                  RENO (0.1kt)
                                                                                                   MiniBoone (0.7kt)
                                                                                                   JUNO (20kt)
                               60
                                                                                                   LENA (50kt)
           Angular error [°]

                               50

                               40

                               30

                               20

                               10

                                    2    4    6         8       10      12     14   16   18   20
                                                  Distance to the Supernova [kpc]

                                Figure 4. Angular uncertainty (68% C.L.) as a function of a galactic supernova
                                distance for different existing and proposed detectors along, with their associated
                                masses (Fischer et al., 2015).

                                Figure 5. Left: Sky area determined at 1σ by combining IceCube timing information
                                with Super-K, assuming normal hierarchy and the (Hudepohl, 2014) model for a
                                supernova at 10 kpc. The true direction is shown with a black dot. Right: Sky
                                area determined by combining IceCube, DUNE, JUNO, and Hyper-K (Linzer and
                                Scholberg, 2019).

the direction of the CCSN. The vector ~n is calculated from the right ascension, α,
declination, δ, of the source in the geographic horizontal coordinate system, and the
event Greenwich mean sidereal time, γ, expressed as an angle (Coleiro et al., 2020):
                                ~n = (− cos(α − γ) cos δ, − sin(α − γ) cos δ, − sin δ).                                   (2)
For simplicity, in recent studies (Brdar et al., 2018; Linzer and Scholberg, 2019; Coleiro
et al., 2020) γ has been fixed to 0◦ ; we note that the results are expected to be
qualitatively insensitive on the choice of such parameter, especially for supernovae in
the galactic center (where α is large). For a known ∆tij and d~ij , Eq. 1 defines a cone
that has a thickness 2δ(cos θij ) due to the uncertainty δ(∆tij ). Typically, ∆tij ≈ 30 ms
for pairs of neutrino detectors since the Earth diameter corresponds to a time delay of
CONTENTS                                                                                                                                                     22

                                 90°                                                                           90°
                                              JU δ=60°                                                                             δ=60°
                          erK
                                                NO                                                    erK                     JU
                       p                           -S
                                                      up                                         Sup                            NO
                   -Su                                   erK
                                                               30°                          v A-                                  -S
                                                                                                                                          up      30°
                vA                                                                       NO                                                 erK
              NO                                                                                                                vA
                                                                                                                               O
   -180°   -120°     -60°        α=0° OvA 60°     120°          0° 180°   -180°       -120°      -60°         α=0°        IC-N
                                                                                                                             60°          120°     0° 180°
                                       N
                                   IC-                                                                       IC-
                                IC                                                                              Su
                                  -S            LVD                                                                  pe
                                     up
                                       erK  IC-                                                                        rK          VD
                                                               -30°
                                                                                                JUN                            IC-L               -30°
                   JUN                                                                                O-
                         O-I                                                                            IC
                            C                   -60°                                                                               -60°
                                -90°                                                                          -90°

                   Figure 6. Regions constrained at 1σ CL by two-detector combinations, adopted from
                   Ref. (Brdar et al., 2018). The left and right panels show scenarios of the supernova
                   core-collapse into a neutron star or a black hole, respectively. All regions expectedly
                   overlap at the supernova location (black dot) which are set in the Galactic center.

∼40 ms. The uncertainty δ(∆tij ) can be evaluated for each detector pair as in (Linzer
and Scholberg, 2019; Coleiro et al., 2020); or as δ(∆tij ) = Max(δti , δtj ), where δti , δtj
are each detector uncertainties defined independently as in (Brdar et al., 2018). The
probability that a test position in the sky (α, δ) coincides with the equatorial coordinates
of the CCSN can be evaluated with the following χ2 function:
                                                                          !2
                                             ∆tij (α, δ) − ∆tdata
                                                             ij
                   χ2ij (α, δ) =                                                  ,                                                                          (3)
                                                   δ(∆tij )
the minimum of the function gives the best estimate for the angles (α, δ) for the searched
CCSN location in the sky.
    Different detector pairs can be combined into a total χ2 by summing each
contribution:
                                       i
CONTENTS                                                                                          23

        Experiments      major process                    target         δt   δt (BH)
         Super-K       ν e + p → e+ + n                32 kt H2 O      0.9 ms 0.14 ms
           JUNO        ν e + p → e+ + n               20 kt Cn Hm      1.2 ms 0.19 ms
           DUNE     νe + 40 Ar → e− + 40 K∗            40 kt LAr       1.5 ms 0.18 ms
           NOνA        ν e + p → e+ + n               14 kt Cn Hm      1.4 ms 0.24 ms
           CJPL        ν e + p → e+ + n                 3 kt H2 O      3.8 ms 0.97 ms
          IceCube         noise excess                     H2 O         1 ms  0.16 ms
        ANTARES           noise excess                     H2 O        100 ms 32 ms
         Borexino      ν e + p → e+ + n               0.3 kt Cn Hm     16 ms   5.5 ms
            LVD        ν e + p → e+ + n                1 kt Cn Hm      7.5 ms 2.4 ms
        XENON1T       coherent scattering                 2 t Xe       27 ms   10 ms
         DARWIN       coherent scattering                40 t Xe       1.3 ms 0.7 ms
              Table 1. A summary of supernova neutrino arrival time uncertainties (δt) estimated
              in Ref. (Brdar et al., 2018). In the second and third columns, the main detection
              channel as well as the target are shown for each of the experiments listed in the first
              column. The next (last) two columns show the δt values for the galactic supernova
              core-collapse into a neutron star (black hole).

for the event rates were assumed. See also (Hansen et al., 2020) for a recent timing
analysis.
     Fig. 6 shows the 1σ regions of supernova directions constrained by several two-
detector combinations. The left and right panel correspond to the case of core-collapse
into a neutron star and black hole, respectively. Table 1 summarizes arrival time
uncertainties for a number of present and future neutrino detectors, for neutron star
and black hole final states. The advantage of this method is most evident in cases with
rapid temporal variation, in particular the sharp cut-off in the flux arising from the
formation of a black hole. Namely, it was found in (Brdar et al., 2018) that a small
fraction of events around the cut-off chiefly determines the timing uncertainty in this
scenario (while the events around the onset were also considered, their effect turned out
to be marginal for obtaining δt in the performed statistical analysis). The disadvantage
of the applied fit is its dependence on the theoretical prediction of the flux.
     Ideally, one would want to exploit advantages of both first-event method (Linzer
and Scholberg, 2019) and the χ2 fit of the full-spectrum (Brdar et al., 2018). For instance,
including the first couple of events in the fit (or last few events in case of the black hole
scenario), would be less model-dependent than the latter and statistically more robust
than the former.
     In order to reduce model dependency for the χ2 fit of the full light-curve, direct
matching of the detected neutrino light-curves has been explored in (Coleiro et al., 2020)
using two different techniques to evaluate the signal arrival time and its uncertainty: χ2
and normalized cross-correlation. The results reproduced in Fig. 7 and Table 2 show that
an uncertainty area of ∼70 deg2 (at 1σ level) in the sky can be achieved when combining
four current and near-future detectors sensitive to IBD (IceCube, KM3NeT-ARCA,
You can also read