Sensitivity of multi-PMT Optical Modules in Antarctic Ice to Supernova Neutrinos of MeV energy

New optical sensors with a segmented photosensitive area are being developed for the next generation of neutrino telescopes at the South Pole. In addition to increasing sensitivity to high-energy astrophysical neutrinos, we show that this will also lead to a significant improvement in sensitivity to MeV neutrinos, such as those produced in core-collapse supernovae (CCSN). These low-energy neutrinos can provide a detailed picture of the events after stellar core collapse, testing our understanding of these violent explosions. We present studies on the event-based detection of MeV neutrinos with a segmented sensor and, for the first time, the potential of a corresponding detector in the deep ice at the South Pole for the detection of extra-galactic CCSN. We find that exploiting temporal coincidences between signals in different photocathode segments, a $27\ \mathrm{M}_{\odot}$ progenitor mass CCSN can be detected up to a distance of 269 kpc with a false detection rate of $0.01$ year$^{-1}$ with a detector consisting of 10000 sensors. Increasing the number of sensors to 20000 and reducing the optical background by a factor of 140 expands the range such that a CCSN detection rate of $0.08$ per year is achieved, while keeping the false detection rate at $0.01$ year$^{-1}$.


Introduction
A core-collapse supernova (CCSN) explosion is the final stage in the evolution of stars featuring masses greater than ∼ 8 solar masses (M ) [1]. Neutrinos play a crucial a email: c.lozano@wwu.de (corresponding author) role during these explosions, carrying away most of the energy that is radiated during the collapse in a short burst of ∼ 10 s [2,3]. In 1987, the first and so far only neutrinos from a CCSN were detected. This supernova, named SN 1987A, exploded in the Large Magellanic Cloud, a satellite galaxy of the Milky Way at a distance of 51.4 kpc from Earth, emitting a burst of neutrinos. During this episode, 25 neutrino events with estimated energies of ∼15 MeV were detected in temporal coincidence by three different neutrino observatories: twelve by the Kamiokande neutrino detector [4], eight by the IMB detector [5] and five by the Baksan scintillation telescope [6]. Even though this detection confirmed the general picture of CCSN explosions, more than 30 years later a detailed picture of the physics of the core collapse is still missing.
Neutrino telescopes that use abundances of natural water or ice as detection medium such as Antares [7] in the Mediterranean Sea, Baikal [8] in Lake Baikal or IceCube [9] at the South Pole have studied neutrinos in recent decades. The deep ice at the South Pole, with its high transparency and low radioactive contamination, has proven to be an excellent site for the detection of high-energy atmospheric and astrophysical neutrinos (see e.g. [10,11]). The IceCube Neutrino Observatory instruments one cubic kilometer of South Pole ice in depths between 1450 m and 2450 m with 5160 Digital Optical Modules (DOMs). DOMs are spherical glass pressure vessels equipped with one 10-inch photomultiplier tube (PMT) facing downwards and associated read-out electronics [12]. The modules are mounted on strings with a horizontal spacing of 125 m and a vertical distance of 17 m between modules. Neutrinos are detected indirectly via secondary charged particles which are produced after a neutrino interacts in the ice or bedrock below the detector. During their passage, these particles induce Cherenkov light emission in the ice, which is detected by the photomultipliers.
The current string layout of IceCube is optimized to reconstruct the energy and direction of neutrinos with energies above ∼100 GeV. The detection of light from MeV neutrinos by more than one module, however, the prerequisite for event-based reconstruction, is very unlikely, since the produced low-energy secondary particles travel only a few centimeters in the ice. Therefore, IceCube does not have the capability to detect MeV supernova neutrinos individually. Nevertheless, a nearby CCSN can be detected by the observation of an increased counting rate in all DOMs in a time window of ∼ 10 seconds, corresponding to the supernova neutrino burst [13].
Based on the success of IceCube, new extensions to the original detector are being planned that will further improve its performance. The IceCube upgrade [14] is scheduled for installation in the Australian summer of 2022/23 and consists of about 700 additional optical modules distributed across seven strings. This extension will significantly expand IceCube's sensitivity to neutrino oscillation parameters using atmospheric neutrinos. The construction of a large detector for high-energy astrophysical neutrinos, IceCube-Gen2 [15], is planned to begin in the second half of the decade and is designed to encompass a large volume with about 10 000 optical sensors on 120 strings. Beyond a mere increase in the number of sensors, novel module types with segmented photosensitive areas are being developed for IceCube-Gen2. It is expected that such modules will provide increased sensitivity not only for high-energy neutrinos, but also for MeV neutrinos. The segmentation will enable a new approach to study low-energy events using temporal coincidences between photons within the same module (local coincidences), improving the sensitivity to CCSN detection.
While closely spaced (unsegmented) sensors can also be used to detect CCSN neutrinos [16], this scenario is not applicable to very sparsely instrumented neutrino telescopes like IceCube-Gen2, which are optimized for high-energy astrophysical neutrinos. The KM3NeT Collaboration [17] is developing a local-coincidence method for the segmented sensors of their sparse detectors in the Mediterranean Sea [18]. In this work, we investigate CCSN detection via local coincidences for a detector in the South Pole ice equipped with 10 000 multi-PMT digital Optical Modules (mDOMs) [19] which has been developed for the IceCube Upgrade [20].
The paper is organized as follows: in Section 2 an introduction to the mDOM concept and its features is provided. In Section 3 we describe the simulation of su- pernova neutrinos used for the sensitivity studies. The final Section 4 describes a set of trigger conditions for MeV neutrino identification and background suppression, which is subsequently used to estimate the sensitivity of a detector equipped with mDOMs for the detection of extra-galactic CCSNe.

A multi-PMT optical module for the South Pole ice
The mDOM [19], depicted in Fig. 1, features 24 PMTs with 80 mm photocathode diameter housed inside a pressure vessel. PMTs and support structure are coupled to the pressure vessel with optical gel. In contrast to the single large PMT in the current IceCube DOM, the PMTs in the mDOM cover the whole solid angle with high homogeneity. They are surrounded by reflector cones to further increase the PMT's effective area. This design results in several advantages: -Increased sensitive area per module by more than a factor two; -A near uniform 4π coverage; -PMT orientations provide intrinsic information about the photon direction, thereby improving event reconstruction; -Simultaneous photon detection in several PMTs of the same module (local coincidence) which can be utilized to identify neutrino interactions over background.
In particular, the last point has the potential to distinguish photons from interacting MeV neutrinos from detector background. In this context, each module can be treated as an individual detector which performs an event-by-event detection. Thus, the sensitivity scales with the number of sensors and the results of our study are independent of the inter-module spacing.

Simulation
A GEANT4 [21] Monte Carlo simulation of the mDOM [22] has been adapted and is used to study the sensitivity of the module to MeV supernova neutrinos. The simulation takes into account the interactions of MeV neutrinos in the ice and the mDOM response, as well as the optical background produced by radioactive decays in the glass of the pressure vessel.

Detector simulation
The mDOM simulation includes all relevant mechanical components of the module: pressure vessel, PMTs, reflector cones, support structure and optical gel. The PMT is modeled as a solid glass body with the photocathode on the inside. Each component has its own material properties. The support structure of the PMTs is simulated as a total absorbing massive object. When a photon reaches the surface of the photocathode it is removed and saved. Its detection probability is given by the quantum efficiency (QE) of the PMT. A further simulation of the PMT or front-end electronics is not included. The medium around the module was modeled after the properties of the South Pole ice [23]. These vary with depth and have a layered structure reflecting climatological changes or volcanic activities in the past [24]. In the simulation, absorption and scattering are parameterized as a function of the mDOM depth, however, on the scale of the simulated volume around the module, the ice properties are assumed to be constant. The simulations are performed using a single mDOM and afterwards are scaled to a detector with 10 000 modules. The varying ice properties in depths between 1400 to 2490 m 1 are taken into account by using the effective per-module volume (see section 3.3). With a vertical module spacing of approximately 13 m in the IceCube-Gen2 detector, shadowing effects can be neglected. 1 IceCube-Gen2 strings will probably be longer and extended into greater depths, increasing the number of modules in cleaner ice and thus the detector performance. However, in this work, the module distribution is limited to depths with measured ice properties by IceCube's LED calibration system [23].

SN neutrino burst
The gravitational collapse of type II supernovae ends in a neutron star or, if the mass of the progenitor is larger than ∼ 10 M , a stellar-mass black hole might be the outcome [25]. During this collapse, electron capture by protons produces a temporally sharp electron neutrino burst lasting ∼ 10 ms, which is called neutronization peak. After that, the main reaction that produces neutrinos and antineutrinos of all flavors is pair production. The neutrino energy spectra are parameterized as [26,27]: where E ν and E ν are the neutrino energy and its mean, respectively, t is the time after the collapse, and α(t) can be calculated from the moments of the neutrino energy spectrum, In this work, we employ fluxes from two different CC-SNe models [28], which are based on the LS220 (Lattimer -Swesty) Equation of State (EoS) [29]: The luminosity, mean energy, and mean squared energy as a function of time from the models are shown in Fig. 2.
Upon reaching Earth, some neutrinos will interact within the instrumented volume. The most important reactions for a detector in the ice are elastic neutrinoelectron scattering (ENES) for electron neutrinos originating from the neutronization peak, and inverse beta decay (IBD) which has the largest cross section of all MeV neutrino interactions in ice [13,30]. The generated e ± usually exceed the Cherenkov threshold and therefore produce visible light that propagates through the ice and can be detected by the optical modules. Note that we do not take into account interactions with oxygen nuclei or the scattering of other neutrino flavors, which reduces the event rate by about 5% (see Table 1    energy neutrinos produce higher energy secondary particles in the ice, which in turn produce more photons, increasing the probability of being detected in one or more PMTs. Thus, although the mean energy of the CCSN neutrinos range mainly from about 5 to 15 MeV (see Fig. 2), the larger effective cross section of the higher energy neutrinos and their higher Cherenkov photon yield lead to a mean energy of about 25 MeV for the detected neutrinos. Results shown in this work do not include any oscillation scenario, which is a conservative choice. In both normal and inverted mass ordering scenarios a larger number of detected events is expected than in the nooscillation case, with the inverted scenario providing the largest signal because energeticν x would oscillate intoν e [13].

Effective volume
The effective detection volume of a mDOM for MeV neutrinos is calculated for different depths with the simulation described in Section 3.1. Electrons are injected homogeneously into a sufficiently large spherical volume of South Pole ice surrounding the module. As they travel through the ice, they induce Cherenkov radiation that can be detected by the module's PMTs. The effective volume is given by where N det and N gen are the number of detected and generated electrons, respectively, and V gen is the generation volume. Electrons producing at least one photon registered by a PMT are considered detected. A sufficiently large generation volume means that an increase in the generation volume does not increase the effective volume. To reduce computational effort, the mDOM is assumed to be spherically symmetric, requiring particles to be simulated only from one direction. This approximation is sufficient for the re-scaling of the effective volume used later in the analysis. The same light production is assumed for electrons and positrons.
Simulations at different depths have been performed for an electron energy of 25 MeV, where each depth has a corresponding absorption and scattering length taken from [24]. The effective volume as a function of the module depth is shown in Fig. 4. Table 1 presents a comparison to the effective volumes of the IceCube detector, which was calculated in [13]. For this result, it is assumed that the 10 000 Table 1 Effective volume and effective mass for IceCube and a detector equipped with 10 000 mDOMs. Quantities have been obtained for the detection of e − with an energy of 25 MeV. The results for the mDOM detector include statistical uncertainties only, while the IceCube results also include systematic uncertainties. IceCube data taken from [13].
IceCube with DOMs 725 ± 95 (3.8 ± 0.5) × 10 6 3.5 ± 0.5 Detector with 10 000 mDOMs 1776 ± 6 (1.705 ± 0.005) × 10 7 15.71 ± 0.05 mDOMs are evenly distributed over the same depths as the IceCube modules. An increase of about a factor 2.4 for the effective volume per module is obtained for the mDOM (in agreement with the larger effective photocathode area) resulting in a factor ∼ 4.5 larger effective volume of a detector equipped with 10 000 mDOMs compared to the current IceCube detector.

CCSN simulation
In order to perform the simulation efficiently, positrons from IBD and electrons from ENES are generated in the ice at random positions inside the simulated volume following the model described in section 3.2. It is assumed that the neutrino flux comes from the zenith. To simulate the events, a time t of the flux is sampled according to the flux distribution from Fig. 2. Using the associated mean energy and mean squared energy, the energy spectrum at each point in time is calculated using Eq. 1. Afterwards, a single energy value from the spectrum is sampled according to its probability. Next, the corresponding angular cross section for this energy and the interaction that is being simulated is calculated, and a direction for the generated particle is sampled. Detected events are weighted according to their interaction probability and the total flux of (anti) neutrinos through the simulated volume. The weight W is given by The first element W int (E) represents the interaction probability of the neutrino traveling through the simulated volume. The generation volume surrounding the module is simulated as a cylinder of South Pole ice with a radius of 20 m and a length of 40 m with the mDOM in its center. The cylinder base is facing the SN, which simplifies the weight calculations. Thus W int is given by where σ(E) is the total cross section of the simulated interaction, n target is the number of targets per volume in the ice for such an interaction, and l = 40 m is the length of the simulated volume. Note that W int (E) is different for each simulated particle since it depends on the energy.
The flux weight W flux is given by where N gen is the number of particles that have been generated, r = 20 m is the radius of the cylindrical cap of ice facing the SN, and d is the distance at which the SN is being simulated. L(t) and E(t) are the fluxes and mean energies from the CCSN models. Finally, W eff accounts for the fact that the simulation is performed for a single mDOM at a particular depth z sim . Since the ice properties vary with depth, the results are scaled using the mean effective volume of the whole detector: where N mDOM is the number of mDOMs that are being simulated and V eff (m, z sim ) is the effective volume for multiplicity m at the depth where the simulation is performed. We define multiplicity as the number of PMTs per module that have detected photons from the same event. The effective volume does not scale in the same way for events detected in a single PMT as it does for events detected in multiple PMTs. This is because coincident events are generally closer to the module than events detected in a single PMT, and thus are less affected by a variation in absorption length. Therefore, events are assigned a different effective volume weight V eff (m) depending on the multiplicity m, with each one calculated using the method of Section 3.3 with the appropriate multiplicity condition. The limited size of the simulated volume, chosen to speed up the simulations, does not cover the whole sensitive volume for a detection threshold of a single photon registered in a PMT. Therefore, the following calculations slightly underestimate the sensitivities. However, this does not impact the conclusions in this work, since our main results are based on events detected in coincidence in several PMTs. The distribution of event multiplicities for a CCSN simulation at 10 kpc distance are shown in Fig. 5. For the high-mass CCSN we obtain more than 2.4×10 6 neutrinos with at least a single detected photon and ∼ 1.4 × 10 6 for the low-mass progenitor. In both cases, more than 12% of the detected events are in local coincidence, i.e. at least two photons are registered in different PMTs, while the value is slightly lower for the low-mass progenitor. Figure 6 shows the distribution of time differences between the first and last registered photon for events detected in coincidence. It demonstrates that ∼ 85% of the local coincidences occur in less than 1 ns between the first and last registered photon, and ∼ 99% within 10 ns (note that the simulation does not include PMT pulse generation and readout). Hence, photons from the same event arrive almost simultaneously at different PMTs. This property can be exploited to suppress background and improve the identification of individual MeV neutrinos.

Identification of MeV neutrinos from CCSNe over background
In this section, a study of local and global coincidence conditions is presented with the aim to optimize the signal-to-background ratio. To this end, we define trigger conditions for identifying single neutrino interactions based on fast local coincidences, as well as a minimum number of such interactions for identifying a supernova outburst. The trigger conditions are defined as follows: a) Once a PMT detects a photon, a time window of length ∆t coin is opened. If within ∆t coin at least m−1 different PMTs in the same module detect one or more additional photons the trigger condition is fulfilled and a neutrino event is detected; b) After condition a) is satisfied, a second time window opens with ∆T SN = 10 s. If N ν − 1 further neutrino events according to a) from any mDOM in the detector fall into this time window, a supernova detection is claimed.
The trigger conditions are designed to be strict enough to effectively suppress backgrounds and keep the false detection rate at an acceptable level. In the following sections, the main background sources are discussed and an estimate of the sensitivity of the method to extra-galactic CCSNe is presented.

Background sources
Due to the low radioactivity of the deep ice at the South Pole, background events originate almost exclusively from the intrinsic noise of the optical module and from the interaction of neutrinos of similar energy in the ice from other sources. These backgrounds are described in the following.

PMT dark rate
PMTs produce a measurable signal (dark rate) even if operated in complete darkness [31]. The dark rate of the PMT model in the mDOM at −30°C, a typical operating temperature in the deep ice, was measured with ∼ 30 s −1 for a threshold of 0.3 pe [32]. The exact value varies from PMT to PMT and ranges up to ∼ 50 s −1 . We assume a conservative dark rate of f d = 50 s −1 for the purpose of this study. Following a similar calculation performed in [16], the probability that the PMT produces a dark rate signal within a time window ∆t coin is given by the complementary probability of not registering any photon: Once a photon is detected, at least m − 1 other PMTs in the module are required to register dark counts in the time window ∆t coin to satisfy trigger condition a). The probability for this can be calculated again using the complementary probability where N PMT = 24 is the number of PMTs in a mDOM and B cum (m|n, p) = m k=0 n k p k (1 − p) n−k is the cumulative binomial distribution for m successes out of n tries when the probability for success is p.
The rate at which trigger condition a) is satisfied due to the PMT dark rate in a module anywhere in the detector is given by with N tot being the total number of mDOMs in the detector.

Radioactive decays in the pressure vessel glass
Scintillation and Cherenkov radiation from radioactive decays in the glass of the pressure vessel are the main background source in IceCube DOMs [9] and the mDOM. The radiation levels of the borosilicate glass used in the mDOM have been measured 2 and are listed in Table 2 assuming secular equilibrium between isotopes within the natural decay chains. Radioactive decays in other parts of the module, e.g. the gel, are found to be negligible. A GEANT4 simulation is performed to estimate the noise contribution of these decays in the mDOM. Within a simulation run, the number of isotopes decaying inside a time window of 20 min is drawn randomly from a Poisson distribution. All decay chains are simulated independently and the photons reaching a PMT are saved after taking its QE into account. Subsequently, the photons are mixed within the time window. Temporal correlations between mother-daughter isotopes are preserved if the decay time is within the time window. Otherwise, the decay time is changed to a random time within the window so that the number of decays remains constant (following the secular equilibrium argument). Finally, the time-ordered hits are checked for coincidences.

Solar neutrinos
Solar neutrinos with an energy above the Cherenkov threshold in ice originate mainly from β + -decays of 8   similar to that of supernova neutrinos. Their spectrum is simulated in GEANT4 with a total flux of Φ νe = 1.7 × 10 6 cm −2 s −1 and Φ νµ,τ = 3.3 × 10 6 cm −2 s −1 using the data from [34]. We consider only the elastic scattering process with electrons in the ice, since it dominates the interaction rate by two orders of magnitude [35].

Further background sources
The following background sources are not included in this study: -Cosmic-ray induced atmospheric muons: these highenergy muons can penetrate into the instrumented volume, generate photons, and produce a detection pattern in a mDOM similar to that of MeV neutrinos. It is expected that due to their extended signature in the detector, the majority of these muons can be identified and rejected, resulting in a short dead time (not accounted for in this work). However, some of them might escape the identification algorithm and thus become part of the background for this analysis. In fact, for KM3NeT, atmospheric muons cause the majority of background events at multiplicities above six [18]. On the other hand, the scattering length in ice is much shorter than in water. Hence, photons produced simultaneously by a muon arrive spread out in time at a module. Therefore, for our case, we expect a smaller contribution of these muons to the background. Corresponding studies require a full detector simulation, though, which goes beyond the scope of this work. -Muons or electromagnetic showers from low-energy atmospheric neutrinos: these neutrinos have energies in the range of 10 MeV -1 GeV, but their interaction rate is only ∼1 s −1 in a Gton volume of ice [36].
The expected background event rates for different multiplicities in a time window of ∆t coin = 20 ns in a detector with 10 000 mDOMs are depicted in Fig. 7. The largest contribution to the background comes from radioactive decays in the pressure vessel glass. Scintillation is the main background source for the overall detection rate, producing ∼ 98% of the detected light after a radioactive decay. Nevertheless, since this light Table 3 Number of signal and background events and effective volume for different multiplicities (see trigger conditions in Section 4; ∆t coin = 20 ns, ∆T SN = 10 s), for a detector with 10 000 mDOMs. is emitted over a long period of time, the probability for high multiplicity coincidences is low compared to that from Cherenkov emission. For m ≥ 6, this allows the use of simulations without scintillation, which significantly reduces computational power and, as a consequence, increases statistics by a factor of ten. A total of 4500 days of radioactive background in a single mDOM is simulated for the no-scintillation case. Note that the cutoff at high multiplicities for radioactive and solar background in Fig. 7 is due to limitation in statistics, since these events are very rare. Due to the lack of noise events at higher multiplicities, only trigger conditions up to a threshold of m = 8 are considered.

Multiplicity
In particular, the 232 Th chain produces, by far, the most important contribution to the background due to 208 Tl, which decays into an excited state of 208 Pb that produces a γ of ∼ 2.6 MeV. This γ can either interact with the vessel glass or the surrounding ice. In both cases, some of its energy is transferred to an electron, which then generates Cherenkov radiation that produces a detection pattern in the module similar to MeV neutrinos.
Note that ∆t coin = 20 ns is still conservative with respect to the time resolution of the module, which we expect to be around 10 ns. As the coincidence time of most photons from supernova neutrino interactions is lower than 1 ns (Section 4.1.4), reducing the time window to a few nanoseconds would further suppress the background. Table 3 shows the number of detected events for signal and background for ∆t coin = 20 ns and ∆T SN = 10 s as well as effective detector volumes for a detector with 10 000 mDOMs. Theses conditions will also be used for the sensitivity estimation in the next chapter.

CCSN sensitivity
The (13) Here we neglect triggers from the combination of background photons from different sources as their contribution is small. With this background, the expected number of false supernova detections N fSN within an observation time δt is given by (14) where P cdf is the cumulative Poisson distribution.
On the other hand, the probability for a supernova at distance d to produce N ν neutrinos that meet trigger conditions a) is (note that these neutrinos are registered within ∆T SN by construction) where µ 0 is the expected number of triggered neutrino events from the simulation of a CCSN at 10 kpc for the chosen multiplicity m condition. Table 4 presents trigger conditions for which the false CCSN detection rate per year (δt = 1 year) is ≤ 1 and ≤ 0.01, respectively. The distance at which a CCSN is still detectable with a probability of 50% for such conditions is also listed. Due to the high rate of low multiplicities from background sources, the best results are obtained for high multiplicities, with the largest distance reached at m ≥ 7. The trigger condition (m ≥ 7 , N ν ≥ 7) can be used to send supernova alerts with very high confidence (about one false detection per century), and identify CCSN at a distance of 269 kpc with 50% probability. With a relaxed set of conditions of (m ≥ 7 , N ν ≥ 6), SNe up to 291 kpc can be detected with less than one false CCSN detection per year. A detailed discussion of acceptable false detection rates goes beyond the scope of this paper but the assumed rates are rather conservative according to [37]. Figure 8 shows the probability from Eq. 15 as a function of distance that a 27 M CCSNe is detected for the aforementioned cases. These results clearly demonstrate the potential of the method to increase the sensitivity to distant CCSNe compared to the current Ice-Cube detector, which can detect CCSNe up to 50 kpc with a false detection rate of 0.1 year −1 [13].
In the case that other telescopes detect a CCSN, we can also analyse archival data in search for a signal. Let us assume that the time at which the neutrinos should have arrived at the detector is known with an accuracy of one hour [38]. Then Eq. 14 gives the expected number of false supernova detections N fSN generated by the background within δt = 1 h. The probability Q that the background did not produce a false supernova-like signal within this interval is Q = P (0, N fSN (δt = 1 h)).
(16) Figure 9 shows Q in terms of one-sided Gaussian standard deviations for the detection of a 27 M CCSN as a function of the minimum number of neutrino events N ν together with the distances at which such a CCSN can be detected. Only the case m ≥ 7 is shown since it provides better results than other multiplicity conditions. Q drastically increases with the increase of N ν while the range for identifying distant supernovae decreases only moderately. For example, for a number of   Table 4.
detected events N ν = 5 a background origin can be excluded at 3.2 σ, while at least a corresponding number of events will be detected in 50% of cases from a 27 M CCSNe at a distance of 322 kpc. If N ν = 7 events with m ≥ 7 are detected we obtain a 4.9 σ confidence that such signal was not produced by background with a 50% detection probability at 269 kpc distance. The KM3NeT collaboration estimates a 5 σ discovery range of 30 kpc for a CCSN with a 27 M progenitor star [18] if the arrival time of the burst is exactly known 3 . We can make the same assumption and, following their approach, obtain the sensitivity for such a detection from Z = 2((s + b)ln(1 + s/b − s). Here, s is the expected number of signal and b the expected number of background events. The 5σ discovery horizon in this scenario reaches 315 kpc for a 27 M CCSN using m ≥ 7, and 234 kpc for the 9.6 M model. It should be noted that the KM3NeT result is based on a full detector simulation with atmospheric muon background rejection while this study does not. As stated before, we expect the contribution of the muon background in ice to be less significant than in water at high coincidence multiplicities. The trigger conditions applied in this work can also be further optimized, e.g., by reducing ∆t coin or by shortening the ∆T SN window. The Q (σ) Fig. 9 Detection prospects for a CCSN whose time is known to within 1 h (m ≥ 7 , ∆t coin = 20 ns , ∆T SN = 10 s). Left axis: probability in σ that the signal is not produced by background fluctuations. Right axis: distance at which a 27 M progenitor mass CCSN is detected with 10% (upper boundary), 50% (middle mark) and 90% probability (lower boundary), when at least Nν detected events are required.
latter is based on the fact that the flux is expected to be most intense in the first few seconds of the burst.
In this work, we have used the detection horizon as a performance metric. The ranges achieved reach up to several hundred kpc. However, most CCSNe candidates within this range are located within ∼ 50 kpc in the Milky Way or the Small and Large Magellanic Clouds. Beyond, CCSN progenitors are expected again in the nearest galaxy, M31, which is ∼ 800 kpc away, and in more distant galaxies. Using the estimated CCSNe population based on recent observations and scaled to the star formation rate [16], we obtain a CCSNe detection rate of 0.046 year −1 , i.e. one SN about every 20 years (m ≥ 7 , N ν ≥ 7). One way to increase the detection horizon and thus the SN detection rate is to use a pressure vessel with lower radioactive contamination. In fact, this is not an unrealistic scenario, as pressure vessels for IceCube-Gen2 are currently under investigation with radioactive contamination levels about a factor of 100 lower than those assumed in this work. Figure 10 shows the expected SN detection rate for hypothetical detectors with 10 000 and 20 000 mDOMs as a function of the false SN rate and the noise reduction factor. In order to double the CCSN detection rate, the noise level must be significantly reduced and the number of modules increased. Under the challenging scenario that the radioactive noise of the vessel glass could be reduced by a factor of ∼ 140, a doubling of the number of installed modules would allow the false SN detection rate to be kept below 0.01 year −1 while doubling the expected CCSN detection to one every ∼12 years. The same would be achieved with a reduction factor of ∼100 and a relaxation of the trigger conditions allowing a false detection rate of 0.1 year −1 .

Conclusions
In this paper, we have shown that exploiting temporal coincidences between detected photons within a segmented photosensor significantly increases the sensitivity of sparsely instrumented neutrino telescopes to MeV neutrinos and thus distant CCSNe. Due to its negligible optical background, the deep ice at the South Pole is particularly well suited for this purpose. For a detector equipped with 10 000 sensors consisting of 24 3-inch photomultipliers, we find that CCSNe up to a distance of 269 kpc can be identified with 50% probability with 0.01 false SN detection per year.
Increasing the number of installed modules to 20 000 and using pressure vessels with significantly reduced optical background could extend the range such that one CCSN every ∼12 years can be observed. If the arrival time of CCSN neutrinos is known from an independent observation with δt = 1 h, a 27 M CCSN at [322, 269] kpc can be detected in 50% of cases and with a [3.2, 4.9] σ certainty that the signal was not produced by background. We note that our studies, which are based on the simulation of a single sensor, do not account for background from atmospheric muons. However, we ex-pect that this background can be effectively suppressed with tailored selection and reconstruction algorithms.
Since each sensor represents a self-contained detector for MeV neutrino, the sensitivities depend only on the number of photosensors but not on the inter-module spacing. Hence, future sparsely instrumented neutrino telescopes like IceCube-Gen2, optimized for TeV-PeV neutrinos, will not only increase the sensitivity to the high energy universe but, if equipped with segmented photosensors, also significantly increase our reach for the detection of CCSNe with MeV neutrinos beyond the galactic range.