Abstract
Supernova detection is a major objective of the Super-Kamiokande (SK) experiment. In the next stage of SK (SK-Gd), gadolinium (Gd) sulfate will be added to the detector, which will improve the ability of the detector to identify neutrons. A core-collapse supernova (CCSN) will be preceded by an increasing flux of neutrinos and antineutrinos, from thermal and weak nuclear processes in the star, over a timescale of hours; some of which may be detected at SK-Gd. This could provide an early warning of an imminent CCSN, hours earlier than the detection of the neutrinos from core collapse. Electron antineutrino detection will rely on inverse beta decay events below the usual analysis energy threshold of SK, so Gd loading is vital to reduce backgrounds while maximizing detection efficiency. Assuming normal neutrino mass ordering, more than 200 events could be detected in the final 12 hr before core collapse for a 15–25 solar mass star at around 200 pc, which is representative of the nearest red supergiant to Earth,
Export citation and abstract BibTeX RIS
1. Introduction
A core-collapse supernova (CCSN) produces an ∼10 s burst of neutrinos at tens of MeV, which is large enough that it can be detected by Super-Kamiokande (SK) and other neutrino detectors (Abe et al. 2016a) if it is in or near the Milky Way. Much of what is known about galactic supernova explosions (SNe) comes from the detection of supernova (SN) neutrinos in 1987 (Alekseev et al. 1988; Bratton et al. 1988; Hirata et al. 1988). In the case of another SN in or near our galaxy, the current generation of neutrino detectors would be capable of collecting a much larger sample of SN neutrinos, improving our understanding and resolving outstanding questions about SNe. Neutrinos arrive before the electromagnetic radiation produced by an SN, so can generate a warning enabling astronomers to start observing as early as possible. Alert systems already exist for this purpose, linking together many detectors for maximum effect (Antonioli et al. 2004). SK has the currently unique ability of determining the direction of an SN from elastically scattered electrons, which is useful for guiding optical instruments (Abe et al. 2016a); SN neutrino detection is a major goal for SK.
Prior to collapse, as a star approaches the end of its life, the temperature and density increases, causing neutrinos from thermal and weak nuclear processes to become the main source of cooling as their emission is highly temperature dependent. For a nearby star, neutrino emission increases over a relatively short timescale to detectable levels (Odrzywolek & Heger 2010), which could give advanced warning before core collapse occurs. This earlier alert to the astronomy community could aid in observing the early light from an SN. Advanced warning could prevent SK (and other SN neutrino detectors) from missing a nearby core collapse due to planned down-time. SK has run for over 20 yr so far, and has on average >90% up-time, making it especially useful for the detection of SNe, which are rare and could happen at any time. KamLAND, another neutrino detection experiment, has a functioning pre-SN alert system (Asakura et al. 2016).
Furthermore, astrophysical neutrinos with known sources have only ever been detected from the Sun (for a review, see Kirsten 1999), from SN1987A, and blazar TXS 0506+056 (The IceCube Collaboration et al. 2018). The only direct observation of SN neutrinos so far is SN1987A, and neutrino-cooled stars have never been observed, so detection would contribute to our understanding of late stage stellar evolution (Odrzywolek et al. 2007; Kato et al. 2015; Yoshida et al. 2016; Patton et al. 2017b).
In the next stage of the Super-Kamiokande experiment (SK-Gd), gadolinium sulfate will be added to the detector, which will improve the ability of the detector to identify neutrons, and therefore low energy through inverse beta decay (IBD). Previous estimates (Odrzywolek et al. 2004; Patton et al. 2017b; Fernandez Menendez 2017) of SK's ability to detect pre-SN neutrinos have assumed the energy threshold of the SK solar analysis (3.5 MeV positron kinetic energy) (Abe et al. 2016b). In fact, detection will be possible at lower energies, albeit with reduced efficiency.
In this article, the sensitivity of the SK-Gd detector to this pre-SN neutrino flux is assessed. The article is structured as follows. Section 2 describes the pre-SN neutrino emission processes, Section 3 gives necessary information concerning SK-Gd, and theoretical pre-SN flux estimates are discussed in Section 4. The detection efficiency and backgrounds for low energy IBD are assessed in Section 5. Results of the study are presented in Section 6, followed by a conclusion in Section 7.
2. Pre-supernova Neutrinos
In a massive star at the end of its life, the fusion of hydrogen (H) and helium (He) nuclei is insufficient to stabilize the temperature and density of the star. That is, cooling from radiation is greater than heating from fusion, leading to contraction under gravity. The higher density, and hence temperature, then enables the fusion of heavier nuclei, initially at the core of the star, and then in shells propagating outwards. Neutrino emission is strongly temperature dependent. The increased temperature leads to an increased rate of cooling due to neutrino emission, which leads to further contraction and heating, and the fusion of even heavier nuclei (see Woosley et al. 2002 for a review of late stellar evolution leading up to SNe). This proceeds in stages driven by carbon (C), neon (Ne), oxygen (O), and silicon (Si) burning. Si burning creates an iron core, which can lead to a CCSN; these stars are sometimes called pre-supernova stars (pre-SN). Following the ignition of C burning, neutrino emission is the greatest source of cooling (Odrzywolek et al. 2007); so these stars are also referred to as neutrino-cooled stars.
The neutrino-cooled stage of a massive star's life is remarkably short compared to usual astrophysical timescales, with the C burning stage lasting hundreds of years, Ne and O stages lasting under a year, and Si burning lasting under two weeks (Odrzywolek et al. 2004). While these processes occur in the interior of the star, it may be that nothing changes in the star's electromagnetic emissions or at the outer surface of the star, so such a state would not be observed by electromagnetic astronomy (Odrzywolek et al. 2007).
Stars capable of CCSN are usually specified as having a zero age main sequence (ZAMS) mass >8 solar masses (M⊙), i.e., those capable of CCSN. However, lighter stars in this category might not enter the Si burning stage, instead undergoing core collapse with an O–Ne–Mg core (Odrzywolek et al. 2007; Kato et al. 2015). Very massive stars >30 M⊙ may collapse directly to a black holes as failed supernovae, but could still produce an increasing and detectable neutrino flux in late stages prior to collapse (Patton et al. 2017a). Aside from ZAMS mass, stellar evolution also depends on metallicity and rotation; the models of pre-SN neutrino emission considered in this study assume solar metallicity and no rotation effects.
Thermal processes, as well as emission. At these high temperatures a large equilibrium population of positrons exists, leading to neutrino emission by pair annihilation
, which is important for detection prospects due to the high flux and relatively high average energy of the
(Odrzywolek et al. 2007; Patton et al. 2017b).
SK-Gd has a chance of detecting a pre-SN star following the ignition of Si burning, as the rate of emission (Figure 1(a)), and crucially the
average energy (Figure 1(b)), increase as the star approaches core collapse (Odrzywolek & Heger 2010). As shown by Figure 2, pre-SN
emission is at much lower energies than those of an SN. Furthermore, pre-SN emission is over a very long timescale compared to SN emission. The mean
energy is below the IBD threshold, meaning only the tail could possibly be detected through IBD; the increasing average energy will mean that the proportion above threshold will also increase. For a sufficiently nearby star, SK-Gd would see a rapid increase in the rate of low energy IBD candidate events. This would presage a CCSN by hours.
Figure 1. Comparison of L mean energy, of emitted
against time to core collapse in models by Odrzywolek & Heger (2010) and Patton et al. (2017b). Spikes in both quantities are caused by ignition of fusion of heavier isotopes in the core or shell of the star.
Download figure:
Standard image High-resolution imageFigure 2. Total emission vs. energy for from pre-SN and from an SN burst, integrated over the full time range of the respective models. The dotted line shows the reaction threshold of IBD. Pre-SN spectrum shown is a 30 M⊙ model from Patton et al. (2017b), and SN spectrum is 30 M⊙ from Nakazato et al. (2013). 30 M⊙ was chosen so that the two would have the same mass. Note that the pre-SN neutrinos are emitted over a much longer timescale (1000 hr vs. 20 s), and that the IBD cross section is strongly energy dependent.
Download figure:
Standard image High-resolution imageFrom Nakamura et al. (2016; their Table 2), there are 41 red supergiant (RSG) stars with distance estimates within 1 kpc, 16 within 0.5 kpc, and 5 within 0.2 kpc. Famous nearby RSGs include Betelgeuse Peg. Wolf–Rayet stars are also possible supernova progenitors, e.g.,
3. Super-Kamiokande with Gadolinium
SK is a large water Cerenkov detector, and is well described elsewhere (Fukuda et al. 2003; Abe et al. 2014). It consists mainly of a tank of 50 kilotons (kt) of ultra-pure water. The inner detector (ID) is 32 kt, and the fiducial volume (FV) is usually given as 22.5 kt, although in practice a smaller or larger FV is used by different analyses as appropriate. In this paper quoted efficiencies assume the full ID volume. The ID is instrumented with around 11,000 50 cm photomultiplier tubes (PMTs). Charged particles are detected through their emission of Cerenkov radiation.
Super-Kamiokande with Gadolinium (SK-Gd, formerly GADZOOKS!) is the next phase of the SK experiment. The main aim of SK-Gd is to detect the supernova relic neutrino signal within a few years of adding gadolinium (Gd; Beacom & Vagins 2004).
SK-Gd will add gadolinium sulfate () to SK's pure water. Naturally abundant isotopes of Gd have some of the highest thermal neutron capture (TNC) cross sections. TNC on Gd is followed by a
The main channel for detection of at low energy (roughly <10 MeV) in SK is IBD on hydrogen (
), as its cross section is relatively high. The neutron takes a short time to thermalize in water and capture, and travels only a short distance, meaning that the positron and TNC form a delayed coincidence (DC), in which two events are reconstructed within a short time and distance of each other. This method of detection is made possible by the upgrade to QBee electronics described in Yamada et al. (2010) and Nishino et al. (2009). The probability of uncorrelated events producing this signature is low, so neutron tagging allows electron antineutrino events to be distinguished from background events, including neutrino events. The high TNC cross section of Gd makes the time between the prompt and delayed parts of the event shorter than with H (∼20
Note that in low energy IBD, the direction of the incoming cannot be reconstructed from the direction of the emitted positron (Vogel & Beacom 1999), and the number of elastic scattering events will be small for a pre-SN, so this technique will not have any SN pointing ability. The direction of neutron travel cannot be resolved in SK, so a positron-neutron vector cannot be used to infer the
direction either.
It is planned that SK's ultra-pure water will be loaded with gadolinium sulfate in two steps, first to 0.02% by mass, then to 0.2%; leading to 50% and 90% of neutrons capturing on Gd respectively, with the rest mainly capturing on H. This paper assumes 0.2% Gd loading, so it should be noted that SK-Gd will begin with a period of reduced sensitivity.
Research and development of the required technologies for SK-Gd has been undertaken by the EGADS experiment, which has successfully operated a Gd-loaded water Cerenkov detector for over two years at 0.2% gadolinium sulfate loading (Ikeda et al. 2019).
4. Electron Antineutrino Flux
Neutrino emissions are calculated from stellar models. Although there are several sets of predictions published for the flux of from a pre-SN, this study primarily uses the data sets of Odrzywolek & Heger (2010; data downloaded from Odrzywolek, Web) and Patton et al. (2017b; data downloaded from Patton et al. 2019). Patton et al. (2017b) predict similar
total emission rates to Odrzywolek & Heger (2010) for a 15 M⊙ star, as shown in Figures 1 and 3, though the time and energy dependent emission rates differ. The flux estimates of Odrzywolek & Heger (2010) are calculated by post-processing the output of an already existing stellar model, and isotopic composition is calculated by assuming nuclear statistical equilibrium (Odrzywolek 2009). Patton et al. (2017b, 2017a) use a more modern stellar evolution code, which fully couples the isotopic composition to the stellar evolution, and tracks the rates of a larger number of isotopes individually. Isotopic composition especially affects the rate of neutrinos produced by weak nuclear processes.
Figure 3. Comparison of the total number of IBDs in a 12 hr window against time to core collapse, predicted in models by Odrzywolek & Heger (2010) and Patton et al. (2017b). Detection efficiency is not taken into account. A distance of 200 pc is assumed. Note that a reduced time range is used compared to Figure 1. Interpolation has been used to produce equally spaced time and energy points.
Download figure:
Standard image High-resolution imageTime t = 0 is taken to be the moment at which the stellar simulation is terminated (when the infall velocity exceeds some threshold), which can be taken as the beginning of core collapse. Figure 3 contrasts the time dependent predicted IBD rates and the positron true energy spectra from the models considered.
The electron flavor ratio of the neutrinos is affected as it passes through the dense stellar medium. Following Patton et al. (2017b) and Kato et al. (2017), an adiabatic transition is assumed, and the flavor is changed by the Mikheyev–Smirnov–Wolfenstein high resonance, dependent on the neutrino mass ordering (MO). The assumed transition probability is in the normal (inverted) mass ordering case, and
. In the data of Patton et al. (2017b), the
flux is provided. For the data of Odrzywolek & Heger (2010), the initial
ratio assumed to be 0.19 following Asakura et al. (2016). This is assumed to not be energy or time dependent, and includes the effect of vacuum oscillations. Earth matter effects have not been included. Only the electron antineutrino flavor will interact through IBD, the rest of the flux is assumed to be invisible in this energy range at SK.
From the flux as a function of neutrino true energy and time, the expected rate of IBD reactions at SK-Gd are calculated. The effect of distance is simply a factor of 1/R2. An accurate approximation of the energy dependent cross section is used (Strumia & Vissani 2003, Equation (25)). The number of targets for IBD is the number of hydrogen nuclei in the SK ID. Detection efficiency is dependent on the energy of the positron, and energy resolution is taken into account with a smearing matrix calculated from detector simulation.
5. Detection Efficiency and Backgrounds
5.1. Energy Thresholds at SK-Gd
The most fundamental energy thresholds for SK-Gd are E is the nucleon mass difference. Positrons above the Cerenkov threshold may not be reconstructed if there are not enough photons detected. Cerenkov photons detection inefficiencies arise due to attenuation in the water, the photocathode coverage of the detector, and the quantum and collection efficiencies of the PMTs (Fukuda et al. 2003).
In SK, a PMT is considered a hit if the charge collected by that PMT exceeds a threshold. In order to reduce the data rate from dark noise and radioactivity, a threshold is usually applied to the number of hits in a time window (Abe et al. 2016b). The Wide-band Intelligent Trigger (WIT) is an independent trigger system at SK that uses parallel computing to reconstruct very low energy events which do not meet the usual thresholds. It is close to 100% efficient for electrons around 4 MeV total energy within the FV (Carminati 2015). Below this, the efficiency drops, as shown by Figure 4.
Figure 4. WIT efficiency against electron true total energy, evaluated with simulation. Note that true energy is shown, not reconstructed energy. 5000 MC events were generated for each 0.2 MeV true energy bin. The efficiency for events generated in the whole ID is compared to those generated only inside the FV. The Cerenkov threshold and SK-IV solar neutrino analysis threshold (4 MeV total electron energy Abe et al. 2016b) are shown for reference.
Download figure:
Standard image High-resolution imageThere is no way to reliably distinguish positrons at very low energy (<4 MeV total) from the intrinsic radioactive background of SK, which increases by an order of magnitude for each MeV the energy threshold is lowered. Requiring DC with a TNC on Gd can make background rates to IBD manageable, even at the lowest energies SK can reconstruct. Furthermore, even if the IBD positron is not detected and successfully reconstructed, the
5.2. Signal Event Model
Low energy positrons in 0.2 MeV true total energy bins from 0.8 to 7 MeV, as well as
The spectrum of
Real data from the SK detector were used as background noise, into which were injected the PMT hits produced by the simulated events. The hybrid data were then subjected to the same triggering algorithms used by the WIT system. It was assumed that all events had to be independently reconstructed and selected by WIT. In the future, requirements on reconstruction quality could be loosened in the WIT for events in coincidence, or prompt and delayed vertices could be fit simultaneously, slightly improving triggering efficiencies for low energy prompt events.
The important features of a DC are the temporal and spatial separation of its two parts, shown in Figure 5. The distribution of spatial separation is dependent on the energy of the positron through the vertex reconstruction resolution, and is modeled by pairing positrons and neutrons created at the same true locations in the detector simulation. Neutron transport was not modeled, i.e., the neutrons were assumed to capture at the same point at which they were produced. This is a reasonable assumption as, at these low energies, the distance traveled by neutrons before capturing (order 10 cm) is much smaller than the position reconstruction resolution (order 1 m). The distribution of DC times was assumed to have the form , where
Figure 5. Distribution of the reconstructed distance (a) and reconstructed time (b) between simulated event candidates (with 3 MeV prompt positrons) before selection. Signal MC is shown in black and accidental background from data in red, area normalized to one.
Download figure:
Standard image High-resolution image5.3. Background Model
The background to the single neutron channel is mostly events from dark noise and radioactive decays with similar characteristics to TNCs ("fake neutrons"), and to a lesser extent events which are real TNCs not from pre-SN neutrinos. Backgrounds to DC type events include accidental DC between unrelated background events, real DC for certain radioactive decays, and real DC from reactor neutrinos and geoneutrinos.
It is anticipated that rates for intrinsic backgrounds will be known quite precisely once in situ measurements become possible, but estimates have been made during the planning and development of SK-Gd. In order to estimate the rate of fake neutrons intrinsic to SK, events recorded by the WIT system during 6000 hr of normal pure water data from the fourth run period of SK (SK-IV) were used. An initial event selection for neutrons was designed to reduce the rate of fake neutrons to a reasonable level, while retaining as much efficiency for simulated TNC events as possible. These data were also used to estimate the accidental DC rate by searching for pairs of events in DC.
A large part of the intrinsic background of the SK detector comes from radon emanation into the water from materials inside the detector and radioactive decays in the detector materials; this is mainly concentrated at the edges and bottom of the detector (Takeuchi et al. 1999; Nakano 2017), and is reduced by fiducial volume and energy threshold cuts.
Any radioactive impurities left in the gadolinium sulfate will be distributed throughout the detector volume by loading. Much of the work preparing for SK-Gd has been in quantifying these, including with low background counting using high-purity germanium detectors, and working with chemical manufacturers to reduce the level of contamination such that it does not detrimentally impact other SK analyses (Ikeda et al. 2019). In this study, to allow for this additional contamination, backgrounds calculated from SK-IV data, and assumed backgrounds from (
The process , and its equivalent with
, produce neutrons. There was some concern that this could produce a high rate of neutrons, especially if there were a high contamination of
series isotopes that are
decays by neutron emission, each reaction of this kind produces two neutrons, so it would be possible for one of the TNCs to be mistaken for a positron, creating a DC background. Rates were calculated using the SOURCES code Wilson et al. (2002). At expected levels of contamination from the
(<5 mBq kg−1),
(<3 mBq kg−1), and
(<0.05 mBq kg−1) chains (Ikeda et al. 2019), we estimate these processes will contribute 6–12 pairs of neutrons per day before detection efficiency. SF of
also produces on average more than one neutron per fission (Verbeke et al. 2007), and
contamination of the gadolinium sulfate of 5 mBq kg−1, the SF rate is calculated to be 0.6 per hour in the ID before efficiencies. The
Reactor and geo electron antineutrinos are an irreducible background, being the same particles in the same energy range as the signal. The reactor background depends strongly on the number of Japanese nuclear reactors that are running. To account for this, the reactor and geo neutrino flux was calculated using the geoneutrinos.org web app (Barna & Dye 2015), which combines reference models for reactor neutrinos (Baldoncini et al. 2015) and geoneutrinos (Huang et al. 2013). Reactor activity is assumed as the mean values given by IAEA's PRIS database for the years 2010 and 2017 (IAEA 2018). During 2010 most Japanese nuclear reactors were running as normal, however, many were switched off in 2011. Some reactors began to be returned to operation after 2010, so 2017 is taken to represent the lowest flux which is likely in the future, and 2010 the highest. Two reactors at Takahama were restarted in 2017 May and June, so the average July–December was used for 2017. If Japanese nuclear reactor activity returns to 2010 levels, then reactor neutrinos will be an important background to the DC channel of this analysis.
Fast neutrons from cosmic-ray muons are not a concern. Neutrons produced by muons not passing through the detector should be few in number; fast neutrons will mostly not penetrate to the FV before capturing as they would need to pass through >4.5 m of water to do so. Cosmic-ray muons passing through the detector are detected very efficiently; products of spallation are already rejected at SK by vetoing events associated in time and space with a muon track (Zhang et al. 2016). Backgrounds from fast neutrons could be controlled by simply vetoing for 1 ms after each muon, which would introduce negligible dead time.
A fraction of muons create unstable daughter nuclei through spallation, and those that do can be efficiently identified through their light emission profiles (Li & Beacom 2015a, 2015b). Some unstable isotopes produce
Reactor and geo neutrino IBD, spallation daughters decaying by
5.4. Event Selection
The rate at which WIT recorded data in SK-IV was typically on the order of 108 event candidates per day. Most low energy triggers have a reconstructed position near the detector wall, and follow an exponentially falling energy and hit count distribution.
Initial cuts on the reconstructed vertex location and number of hit PMTs are used to select neutron candidates for the single neutron channel, reducing the background rate by a factor of 104. These cuts are based on the number of on-time hits, quality of reconstruction, and reconstructed vertex location. This is 47% efficient for simulated TNC on Gd candidates within the ID. The initial selection was based on the reconstructed vertex location and number of selected PMT hits, and the distance from the reconstructed vertex to the detector wall in the reconstructed direction of the event.
DC events are selected by searching for events close in time that have been reconstructed within the FV. The time between the neutron and positron, and the distance between their reconstructed positions, is shown in Figure 5 for DC signal MC with 3 MeV positrons, and accidental background. Other than the standard FV cut, selection of positron candidates is kept very inclusive, in order to achieve the lowest possible neutrino energy threshold. The efficiencies of these cuts are somewhat energy dependent due to the deterioration of vertex resolution at low energy. A cut is also applied to delayed event candidates based on the distance from the reconstructed vertex to the detector wall, and the number of selected PMT hits. A boosted decision tree (BDT) classifier based on the angular distribution of hits, distance from the detector wall, reconstructed energy, and reconstruction quality is used to select good neutron candidates for both the DC and single neutron channels. The BDT was trained on a random 5% subsample of the WIT background data, with performance evaluated on the whole data set. Although the details of the Gd
BDTs are a multivariate method capable of combining multiple discriminating variables into a single score. Cut criteria (decision trees) are automatically generated to separate signal and background training samples. A weight is assigned to each tree according to the accuracy of the tree. Events which are poorly classified are given boosted weight in subsequent trees. The overall score comes from the weighted sum of the response of all decision trees in the ensemble. The TMVA (Hoecker et al. 2007) software package was used. For a more detailed, rigorous explanation, see Hoecker et al. (2007) and Freund & Schapire (1997).
Cuts were optimized to give the greatest range for efficient detection of a 25 M⊙ pre-SN star using the model of Odrzywolek & Heger (2010) 0.1 hr before core collapse, assuming the more optimistic background and efficiency cases. The distribution of the variables used in the final selection of DC events are shown in Figure 6 for the two largest backgrounds in the most important variables. While the accidental backgrounds are well controlled by the combination of coincidence variables and the neutron BDT, the background from pairs of neutrons is harder to reduce. Selection of single neutrons was based mostly on the BDT classifier, with distributions in the final selection shown in Figure 7. The cuts shown in Figure 6 are delayed candidate BDT score >0, coincidence distance <300 cm, and coincidence time <100
Figure 6. Distribution of DC channel variables in final selection for signal and largest backgrounds. Solid lines show before selection, dotted after. Background is normalized to 12 hr. Signal is normalized to the flux in the final 12 hr before CCSN for a 25 M⊙ star at 500 pc.
Download figure:
Standard image High-resolution imageFigure 7. Distribution of single neutron channel Gd BDT score in final selection for signal and largest background. Solid lines show the distribution before selection, the dotted line shows the cut point. Background is normalized to 12 hr. Signal is normalized to the flux in the final 12 hr before CCSN for a 25 M⊙ star at 500 pc. The distribution of the Gd BDT score is different from that shown in Figure 6 as single neutron events had a higher hit threshold applied in pre-selection.
Download figure:
Standard image High-resolution imageTotal trigger and selection efficiency depends on the energy distribution of the flux. In the final 12 hr before core collapse, the proportion of all IBD events which are triggered and selected is 4.3%–6.7% for DC events, and 9.5%–10% for neutron singles. Using the alternative MC for the
Figure 8. Trigger and selection efficiency for DC events against prompt positron true total energy. Above 4 MeV positrons are efficiently triggered and reconstructed, and so the efficiency is determined by the triggering and selection of the delayed part of the event. The result with the alternative Gd
Download figure:
Standard image High-resolution imageThe selection requirements on the DC time, and the positron quality requirements, are very efficient for signal. DC signal efficiency is lower at low energies due to the lower trigger efficiencies, and worse vertex reconstruction leading to increased DC distance. It should not be alarming that this efficiency is lower and these backgrounds higher than are typically quoted for other SK analyses (such as the supernova relic neutrino analysis), which are at higher energy. At higher energy, reconstruction is much better, and backgrounds much lower. Almost all prompt events in this analysis reconstruct below the energy threshold of any other analysis at SK. If the TNCs were only on H, rather than Gd, there would be no sensitivity at all in this analysis.
Remaining backgrounds for single neutrons, dominated by fake neutrons from the intrinsic radioactivity of the detector, total 66–140 per 12 hr window. For DC events, remaining backgrounds are 5–11 accidental, 0.2–0.4 from SF
5.5. Detection Strategy
The simplest way of searching for a rapid excursion in the candidate event rate is to define some signal time window (e.g., 12 hr), and background window (e.g., 30 days), then perform a hypothesis test. The null hypothesis is that the observed rate in the signal window is consistent with the observed rate in the background window, taking into account Poissonian fluctuations. The alternative hypothesis is simply that the rate in the signal window is higher than that in the background window. A Poisson likelihood for the null hypothesis is calculated for the total detected event rate, combining both the DC and single neutron channels.
The time from the start of a data block for reconstruction, event selection, and hypothesis testing will be around 10 minutes. This sets an estimated typical latency for an alarm.
The length of the signal window should be similar to the timescale over which the event rate would change if a pre-SN star was detected. The choice of signal window size can have a dramatic effect on the sensitivity of the analysis, as it affects the statistical fluctuations, the background level, and the trial factor. A range from 1 to 72 hr were tested and the 12 hr signal window performed best. For the model used, the flux multiplied by the IBD cross section, integrated over the final 12 hr before collapse is summarized in Table 1.
Table 1. Antielectron Neutrino Flux Multiplied by IBD Cross Section for Models Considered
MO | Mass | Model |
![]() |
---|---|---|---|
(M⊙) | (![]() |
||
NO | 15 | Odrzywolek & Heger (2010) | 1.6 |
15 | Patton et al. (2017b) | 1.9 | |
25 | Odrzywolek & Heger (2010) | 3.3 | |
30 | Patton et al. (2017b) | 3.8 | |
IO | 15 | Odrzywolek & Heger (2010) | 0.44 |
15 | Patton et al. (2017b) | 0.53 | |
25 | Odrzywolek & Heger (2010) | 0.93 | |
30 | Patton et al. (2017b) | 1.2 |
Download table as: ASCIITypeset image
A longer background window would always be better due to reduced uncertainty in the background rate; however, the background rate may change slowly over time. A gradual change might be expected, due to the gradual increase in PMT gain, seasonal variations in the radon concentration in the mine air (Nakano et al. 2017; Pronost et al. 2018), changes in water flow affecting radon activity in the fiducial volume, or changes in nearby nuclear power station activity. In this study it is assumed that the background level is known precisely when the signal is detected.
These methods are chosen for the purpose of benchmarking performance. In practice, greater sensitivity can be achieved by properly accounting for the likelihood of the event rate across separate bins, and by assuming a more complicated alternative hypothesis, for example, by calculating the rate of increase of the candidate event rate. Furthermore, a model including background time variation, measurement uncertainty, and correlation between bins, should be included when a large enough background sample has been collected during SK-Gd.
In attempting to get an SN early warning from the detection of pre-SN neutrinos, there are four variables which together describe a detector's performance.
- 1.Alarm efficiency, i.e., the probability of correctly detecting a true pre-SN.
- 2.False positive rate (FPR).
- 3.Expected time of early warning that the detector would provide.
- 4.Expected distance to which the warning would be efficient.
Tolerating a higher FPR would allow for greater range and more warning, but would reduce trust in the warning system; a choice needs to be made on what FPR is acceptable. Note that an SN within 1 kpc is a rare occurrence—on the order of 1 in 10,000 yr based on historical data (Adams et al. 2013). FPR levels shown in this paper are 1 per year and 1 per century (cy.), and are assumed to be set by Poisson fluctuations only. Range is defined as the point at which alarm efficiency is 50%. By formulating the problem in terms of FPR, a trials factor is incorporated.
6. Results
Figure 9 shows the distance to the pre-SN star at which the null hypothesis would be rejected before core collapse. Figure 10 shows the largest amount of time before core collapse at which a pre-SN is expected to be detected. The width of the bands shows uncertainty due to levels of background and the difference between models of
Figure 9. Expected maximum range of detection in the final 12 hr before collapse, for 15 M⊙ and 25 M⊙ stars. The width of the bands shows uncertainty due to levels of background in SK-Gd, and the difference between models of
Download figure:
Standard image High-resolution imageFigure 10. Expected early warning against range of detection, for 15 M⊙ and 25 M⊙ stars. It is assumed the range is the largest distance at which probability of detection is greater than 50%.
Download figure:
Standard image High-resolution imageTable 2. Maximum Range of Detection under Various Assumptions
Max. Range (pc) | ||||
---|---|---|---|---|
Mass | given FPR | |||
(M⊙) | Model | 1 yr−1 | 1/cy. | |
NO | 15 | Odrzywolek | 300–400 | 250–300 |
15 | Patton | 420–600 | 360–500 | |
25 | Odrzywolek | 330–400 | 280–400 | |
30 | Patton | 480–600 | 410–500 | |
IO | 15 | Odrzywolek | 160–200 | 130–200 |
15 | Patton | 220–300 | 190–200 | |
25 | Odrzywolek | 180–200 | 150–200 | |
30 | Patton | 270–400 | 230–300 |
Note. Range is assumed to be the highest distance at which alarm efficiency is greater than 50%. Uncertainty comes from the TNC
Download table as: ASCIITypeset image
Results depend on the neutrino mass ordering (normal NO or inverted IO), the ZAMS mass of the star, the distance to the star, and the background level in SK-Gd. Questions of detector model uncertainty will be resolved by in situ measurements once SK-Gd is loaded. An inverted neutrino mass ordering is detrimental to this analysis as it reduces the fraction of the pre-SN flux.
Discussions of pre-SN stars often focus on
Figure 11. Expected signal events in a 12 hr time window, after detection efficiencies are taken into account as per Section 5. A distance of 200 pc and NO is assumed. Dotted lines show the high and low background assumptions.
Download figure:
Standard image High-resolution imageTable 3. Expected Numbers of Events at 200 kpc in the Final 12 hr before Collapse at SK-Gd
Mass | Single | |||
---|---|---|---|---|
Model | (M⊙) | Neutrons | DC | |
NO | Odrzywolek | 15 | 55–71 | 33–36 |
Patton | 15 | 65–84 | 45–50 | |
Odrzywolek | 25 | 120–160 | 59–65 | |
Patton | 30 | 130–170 | 100–110 | |
IO | Odrzywolek | 15 | 16–20 | 9–10 |
Patton | 15 | 18–23 | 13–15 | |
Odrzywolek | 25 | 34–44 | 17–18 | |
Patton | 30 | 40–52 | 34–37 |
Note. Uncertainty comes from TNC
Download table as: ASCIITypeset image
Table 4. Time at Which Expected Signal Exceeds Threshold, for Some Assumptions Chosen to Represent Betelgeuse
MO | Assumed | Warning (hr) | |||
---|---|---|---|---|---|
Mass | Distance | Given FPR | |||
(M⊙) | Model | (pc) | 1 yr−1 | 1/cy. | |
NO | 15 | 150 | Odrzywolek | 5.3–8.4 | 3.4–6.3 |
15 | 150 | Patton | 7.1–14.1 | 5.1–9.8 | |
25 | 250 | Odrzywolek | 4.7–7.4 | 3.3–5.7 | |
30 | 250 | Patton | 1.0–1.6 | 0.7–1.1 | |
IO | 15 | 150 | Odrzywolek | 0.1–2.0 | 0.0–0.8 |
15 | 150 | Patton | 0.3–4.1 | 0.0–2.2 | |
25 | 250 | Odrzywolek | 0.0–0.6 | 0.0–0.0 | |
30 | 250 | Patton | 0.1–0.4 | 0.0–0.1 |
Note. 30 M⊙ at 250 pc is far from the range of mass estimates for Betelgeuse, but is provided anyway for comparison. Uncertainty comes from the TNC
Download table as: ASCIITypeset image
The KamLAND collaboration published an analysis of their own sensitivity to pre-SN neutrinos (Asakura et al. 2016), which is compared to the expected sensitivity of SK-Gd. KamLAND is a liquid scintillator detector based in the same mine as SK. It has lower energy thresholds than SK, and so would detect IBD events from a pre-SN more efficiently, and has lower background rates. However, the mass of SK is more than 20 times larger than that of KamLAND, so more events are seen in total.
The nominal performance of KamLAND is taken from Asakura et al. (2016). KamLAND background rates are 0.071–0.355 events per day depending on Japanese nuclear reactor power. Events are integrated over a 48 hr sliding window every 15 minutes. Assumed signal rates in the final 48 hr at 200 pc are 25.7(7.28) in the 25 M⊙ case, 12.0(3.38) in the 15 M⊙ case, for the NO(IO) case. The pre-SN models used were those of Odrzywolek & Heger (2010). Not enough information is provided in Asakura et al. (2016) to directly and fairly compare warning times.
Figure 12 shows the probability of detection before core collapse (t = 0) against distance to the pre-SN star. The estimated range for KamLAND is also shown. KamLAND has a latency of 25 minutes, which is not taken into account. The FPR is set to match that of the 3 per 48 hr for 3
per 48 hr for 5
Figure 12. As in Figure 9, but with reference numbers for KamLAND included. The FPR is fixed at the equivalent of 5
Download figure:
Standard image High-resolution imageNext generation liquid scintillator and Gd-loaded water Cerenkov detectors could provide earlier warning to longer distances due to their large target masses. A future dark-matter direct-detection experiment could also detect significant numbers of pre-SN neutrinos through coherent scattering, with the advantage of being sensitive to all flavors (Raj et al. 2019).
7. Conclusion
Electron antineutrinos from a pre-SN star precede those from a CCSN by hours or days, increasing in flux and energy rapidly over a period of hours: this has never been detected. In the next stage of SK, gadolinium loading will enable efficient identification of neutrons, enabling the reduction in the energy threshold for the detection of .
The background rates and signal efficiencies for an SK-Gd low energy analysis capable of detecting pre-SN have been quantified. This requires detection of events below the usual energy thresholds of SK, for which trigger efficiency and reconstruction are poorer, and backgrounds higher. Gadolinium loading is essential to detecting these events. Through a rapid increase in the number of event candidates, additional warning of a very nearby SN can be achieved, and useful information provided about late stellar burning processes that lead up to a supernova.
Based on this and the predicted fluxes of Odrzywolek & Heger (2010) and Patton et al. (2017b), estimates were produced of the distance at which a pre-SN star could be observed, and the amount of additional early warning that could be expected. Uncertainty in the future capabilities of the detector arises mainly from the future internal contamination of the SK detector, which is the main source of backgrounds at low energy. This uncertainty will be reduced once in situ measurements of background become available. An inverted neutrino mass ordering would have a detrimental effect on the range of this technique by reducing the fraction of the
flux leaving the star.
The nearest RSG star to Earth is
A pre-SN alert could be provided by SK to the astrophysics community following gadolinium loading. Future large neutrino detectors will improve the potential range of detection, especially if they have a sufficiently small low energy threshold and the ability to tag neutrons from IBD. It could also be possible to use multiple detectors in combination to provide a pre-SN alert with higher confidence.
We gratefully acknowledge the cooperation of the Kamioka Mining and Smelting Company. The SuperKamiokande experiment has been built and operated from funding by the Japanese Ministry of Education, Culture, Sports, Science and Technology, the U.S. Department of Energy, and the U.S. National Science Foundation. Some of us have been supported by funds from the National Research Foundation of Korea NRF20090083526 (KNRC) funded by the Ministry of Science, ICT, and Future Planning and the Ministry of Education (2018R1D1A3B07050696, 2018R1D1A1B07049158), the Japan Society for the Promotion of Science, the National Natural Science Foundation of China under grants No. 11235006, the National Science and Engineering Research Council (NSERC) of Canada, the Scinet and Westgrid consortia of Compute Canada, the National Science Centre, Poland (2015/18/E/ST2/00758), the Science and Technology Facilities Council (STFC) and GridPP, UK, and the European Unions H2020-MSCA-RISE-2018 JENNIFER2 grant agreement no.822070.