{"id":"07011ea6-6b45-4b60-a901-dfba948519b9","arxiv_id":"1909.00886","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"The authors added an IceCube detector model to the SNOwGLoBES supernova neutrino event-rate software using an effective volume approach.","lead":"This paper adds a simulation of the IceCube neutrino detector to the SNOwGLoBES software package, so researchers can compute how many supernova neutrinos IceCube would detect. It gives the supernova community a standard way to compare IceCube with other neutrino observatories in the same tool.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No benchmark: the 'working' claim rests on an unvalidated effective-volume mapping; the constant 0.95 deadtime factor in Eq. 3.2 likely biases rates by tens of percent.","rationale":"The paper is a modest software methods contribution: adding IceCube to SNOwGLoBES. I read the central claim as 'a working implementation' (Section 4), evidenced by submitted code. What must be true for this claim: the effective-volume renormalization in Eq. 3.2 must produce IceCube supernova rates consistent with the collaboration's own published response, within the intended accuracy. The internal algebra is plausible: setting M*T_I = rho_ice*V_eff,tot(E_k) maps the SNOwGLoBES mass/efficiency product onto the IceCube rate formula Eq. 2.1, and the target-weight table in Section 3.1 handles the water composition. The submitted pull request gives some confidence that the code runs. However, the paper provides no external benchmark. Figure 1 shows only the implementation's own predictions, not a comparison to [3] or [4]. The constant deadtime factor is flagged as a simplification in Section 3.2, but its magnitude is not quantified; from [3] the efficiency varies with instantaneous supernova rate, so the fixed 0.95 can overestimate the peak rate by tens of percent. This is a correctness risk for the main use case, not merely a stylistic gap. I do not see an internal inconsistency sufficient for rejection; the approximations are transparent and the code is available for testing. The appropriate remedy is a benchmark comparison, which matches the reader's CONDITIONAL verdict. I mark agreement as 'partial' because my emphasis is less on the abstract validity of Eq. 3.2 and more on the absence of any quantitative check of it against IceCube's published rates, with the deadtime approximation as the most concrete source of expected bias.","tokens_in":5120,"tokens_out":11375,"duration_ms":122990,"concrete_test":"Using the submitted SNOwGLoBES files, reproduce the IceCube supernova rate curve from Abbasi et al. 2011 ([3]) under identical inputs (e.g., an 8.8 solar-mass LS220 model at 10 kpc, same IBD and electron-scattering channels); compare the time-dependent rate and total counts to the published curve, then rerun replacing epsilon_noise = 0.95 with epsilon_noise(t) = 0.87/(1 + 250 microseconds * r_SN(t)). If peak rates differ by more than 20% or the time profile is distorted, the mapping in Eq. 3.2 is not validated and the implementation should be labeled approximate pending revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that Eq. 3.2 lets SNOwGLoBES compute IceCube supernova rates by setting the detector mass to M = 51600 kton and absorbing rho*V_eff,tot(E_k)/M into the post-smearing efficiency. Mathematically this is a clean renormalization, and the code has been submitted as PR #8, so mechanical operation is credible. The load-bearing gap is that nowhere is this mapping compared with IceCube's own published simulation or data, e.g. [3]. The paper's only figure is its own prediction (Fig. 1), not a benchmark. This matters because the mapping contains approximations with known sign and magnitude: the constant deadtime factor epsilon_noise = 0.95 differs from the time-dependent epsilon_noise ~ 0.87/(1 + tau*r_SN(t)) used in [3], which at peak per-DOM supernova rates of several hundred Hz to 1 kHz reduces the efficiency to roughly 0.7-0.85. The constant factor therefore likely overpredicts the signal by 10-35% during the accretion phase. Without a benchmark, a user cannot tell whether a disagreement with IceCube comes from the effective-volume parameterization, the deadtime simplification, or a coding error. Thus the paper's 'working implementation' claim is conditional on validation that is absent.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents an implementation of IceCube in the SNOwGLoBES package for computing supernova neutrino detection rates. To work within SNOwGLoBES's fixed-mass framework, the authors replace IceCube's position-dependent sensitivity with an effective volume V_eff,tot and absorb the product rho_ice*V_eff,tot into the post-smearing efficiency T_I(E_k) via Eq. (3.2), using a constant fiducial mass M=51600 kton, a constant deadtime efficiency epsilon_noise=0.95, and an assumed 98% DOM livetime. The implementation includes inverse beta decay, neutrino-electron scattering, and oxygen channels, with cross-section choices listed in Table 1. A sample calculation for 20 supernova progenitor masses is shown in Figure 1. The code has been submitted as pull request #8 to the SNOwGLoBES repository.","tokens_in":5376,"tokens_out":9646,"duration_ms":92949,"significance":"If validated, this implementation would fill a genuine gap in SNOwGLoBES by enabling direct, fast comparisons between IceCube and other supernova neutrino detectors within a widely used community tool. The effective-volume-to-efficiency mapping is a clean workaround for SNOwGLoBES's inability to vary detector mass with energy or channel, and the approach is potentially generalizable to other large water Cherenkov detectors such as KM3Net. The submitted code as a pull request is a concrete reproducibility asset. However, the paper's central claim of a 'working implementation' currently rests on an unvalidated effective-volume mapping and an unquantified deadtime approximation, so the practical significance is conditional on the benchmarks and sensitivity studies requested below.","major_comments":[{"comment":"The claim that this is a 'working implementation' (Section 4) is not supported by any benchmark against IceCube's published supernova rates or full detector simulation, e.g., reference [3]. The only quantitative output, Figure 1, is a self-prediction rather than a validation. Because Eq. (3.2) absorbs the effective volume into the post-smearing efficiency, any discrepancy between the SNOwGLoBES output and IceCube's published rates could originate from the effective-volume parametrization, the deadtime simplification, or an implementation error. The manuscript should compare its output with the rates in [3] (or an equivalent reference simulation) for a common supernova model and show agreement within the expected systematic uncertainties.","section":"Section 4 and Section 3.3 (Eq. 3.2)"},{"comment":"The constant deadtime efficiency epsilon_noise = 0.95 replaces the time-dependent expression epsilon_noise ~ 0.87/(1 + tau*r_SN(t)) from [3]. During the accretion phase, where the per-DOM supernova rate can be several hundred Hz to 1 kHz, the time-dependent factor is expected to drop to roughly 0.7-0.85, so the constant factor likely overpredicts the signal by 10-35%. The paper should quantify this bias for the models shown in Figure 1, and should either provide a mechanism for users to supply a time-dependent deadtime factor or explicitly restrict the validity of the default constant value.","section":"Section 3.2 and Eq. (2.1)"},{"comment":"The generation of the post-smearing efficiency functions T_I(E_k) for the inverse beta decay and neutrino-electron scattering channels is not described in sufficient detail for reproducibility. The text says these functions were 'generated using the sources listed in Table 1' but does not specify the exact functional form that converts V_eff,tot(E_k) into T_I(E_k) for each channel, nor how the relationship between the incoming neutrino energy E_j and the produced electron/positron energy E_k is handled in the summation of Eq. (2.4). The authors should provide the explicit algorithm or code so that a user can verify the implementation independently.","section":"Section 3.1 and Section 3.3 (Table 1, Eq. 3.2)"}],"minor_comments":[{"comment":"The fiducial mass M = 51600 kton is an arbitrary normalization scale, not the physical mass of IceCube; the text should state more prominently that the choice of M cancels in the effective rate through Eq. (3.2), to prevent users from misinterpreting this value as IceCube's detector mass.","section":"Section 3.3"},{"comment":"The relationship between N_gamma(E') in Eq. (2.1) and <V_eff,pm>(E') in Eq. (2.3) is not explicit; a reader may be confused about whether N_gamma(E') is given by theta(E'-E_ch)*(E'-E_ch)*C_pm*dN_gamma/dx or is an additional factor. Please clarify the definition of N_gamma(E').","section":"Section 2.2"},{"comment":"The caption should specify the assumed supernova distance, the luminosity normalization, and whether the shown rates include the constant deadtime factor epsilon_noise = 0.95, so that the predictions can be compared with other calculations.","section":"Figure 1"},{"comment":"The sentence 'recent advances such as the development of HitSpooling is expected to allow' contains a subject-verb agreement error; it should read 'are expected to allow.'","section":"Section 3.2"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is written as a concise proceedings contribution, and the authors may argue that the code itself is the validation. However, if this is to be considered a journal publication, the absence of any benchmark against published IceCube rates is a serious gap that should be addressed before the paper can be regarded as establishing a reliable tool. The requested fixes (a benchmark comparison and a quantitative sensitivity study of the deadtime approximation) are well within the manuscript's scope and should not require a fundamentally new approach."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Good to see this. The paper fills a real gap: SNOwGLoBES is the standard quick-comparison tool for supernova neutrino detectors, and IceCube was missing. The design choice in Eq. (3.2), folding the effective volume into the post-smearing efficiency with a fixed fiducial mass, is clean and sensible given SNOwGLoBES's fixed-mass constraint. The channel list and cross-sections lean on established work, and the code is in a pull request, so there are concrete artifacts. Credit where due: it's a useful, honest methods note.\n\nThe soft spot is exactly the one you'd expect: no benchmark. The paper claims a 'working implementation' but shows no comparison with IceCube's own published simulation or data, e.g. the rate curves from [3]. The single figure is just the implementation's own prediction. That matters because the mapping has known approximations. The constant deadtime factor epsilon_noise = 0.95 is a simplification; the time-dependent factor in [3] is epsilon_noise ~ 0.87/(1 + tau*r_SN(t)), which during the accretion peak can drop to ~0.7-0.85. So the implementation likely overpredicts signal by 10-35% in that phase. Also the assignment of effective volume to the post-smearing function is plausible but not cross-checked. These aren't fatal for a software paper—users can rescale and the structure is right—but the 'working' claim is stronger than the evidence.\n\nThe absence of a benchmark also makes it hard to debug. If a user's rate disagrees with IceCube's official sensitivity, they can't tell whether it's the effective volume parameterization, the deadtime simplification, or a coding mistake. A simple reproduction of the results from [3] with this SNOwGLoBES module would fix that. I'd suggest the authors add that, or at least clearly state the expected bias from the constant deadtime factor.\n\nBottom line: modest but genuinely useful contribution for the supernova neutrino community. It deserves proper review. I'd send it to a referee, with the request that the validation gap be addressed. The paper is a solid starting point for adding IceCube to SNOwGLoBES; right now it's a conditional acceptance.","headline":"A useful, clearly described SNOwGLoBES module for IceCube that fills a real gap, but the 'working implementation' claim needs a benchmark against IceCube's own simulation or data before users can trust the predicted rates.","tokens_in":5854,"tokens_out":2034,"would_cite":true,"duration_ms":19341,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"IceCube has been implemented in SNOwGLoBES, so researchers can compute its expected supernova neutrino event rates with the standard detector-comparison package.","keywords":["IceCube","SNOwGLoBES","supernova neutrinos","effective volume","post-smearing efficiency","inverse beta decay","core-collapse supernova","detector simulation"],"falsifier":"Take the same supernova flux used in the paper's demonstration and run a full Monte Carlo simulation of IceCube's DOM response (or use the published IceCube supernova analysis) to produce a hit-rate time series; compare it to the SNOwGLoBES calculation for the same input. Agreement within the expected statistical and systematic errors would support the implementation; a systematic divergence, especially at low neutrino energies or early times, would show that the constant-mass, effective-volume mapping in Eq. (3.2) is not adequate.","tokens_in":4934,"feed_emoji":"🧊","tokens_out":7430,"duration_ms":65966,"temperature":0.7,"pith_summary":"The paper reports a working implementation of IceCube in SNOwGLoBES, the software package used to estimate neutrino detection rates from core-collapse supernovae. The motivation is practical: SNOwGLoBES is a standard tool for comparing how different detectors would respond to a galactic supernova, but it had no IceCube simulation files. The implementation works by giving IceCube a fixed fiducial mass of 51,600 kilotonnes and folding the detector's energy-dependent effective volume into the post-smearing efficiency function. The authors demonstrate the result by computing IceCube signal time profiles for 20 supernova progenitor models under two equations of state. If the mapping is faithful, IceCube can now be included in routine SNOwGLoBES sensitivity comparisons without running a full detector simulation.","feed_headline":"IceCube supernova rates now computable in SNOwGLoBES","feed_subtitle":"The world's largest neutrino detector joins the standard package for comparing supernova responses.","key_machinery":"The load-bearing object is the mapping between IceCube's effective volume and SNOwGLoBES's post-smearing efficiency function. In the SNOwGLoBES rate formula $R_I(t) = M_{\\rm detector} n_{{\\rm weight},I} (\\Delta E)^2 \\sum_{j,k} F_a(E_j) \\sigma_I(E_j) k_I(E_k,E_j) T_I(E_k)$, the factor $T_I(E_k)$ normally describes the probability that an interaction product of energy $E_k$ triggers a detection. Here that factor is set equal to $\\rho_{\\rm ice} V_{\\rm eff,tot}(E_k)/M$ with fixed $M = 51600$ kton, so the energy dependence of IceCube's detection efficiency enters through $T_I$ rather than through a varying detector mass. This lets the implementation reuse SNOwGLoBES's existing channel and cross-section machinery while still capturing the detector's energy-dependent response.","core_discovery":"The central claim is that IceCube has been successfully added to SNOwGLoBES, so that the expected supernova neutrino event rate in IceCube can be calculated with the package's standard tools. The technical workaround that makes this possible is Eq. (3.2): since SNOwGLoBES requires a single constant detector mass, the implementation sets $M = 51600$ kton and embeds all effective-volume information in the post-smearing efficiency $T_I(E_k) = \\rho_{\\rm ice} V_{\\rm eff,tot}(E_k) / M$. The effective volume itself follows the IceCube supernova literature: $V_{\\rm eff,tot}(E') = N_{\\rm DOM} \\varepsilon_{\\rm noise} \\theta(E' - E_{\\rm ch}) (E' - E_{\\rm ch}) C_{\\pm} (dN_\\gamma/dx) \\langle V_{\\rm eff,\\gamma}\\rangle$, with $N_{\\rm DOM} = 0.98 \\times 5160$ and $\\varepsilon_{\\rm noise} = 0.95$. Using SNOwGLoBES's standard water-Cherenkov channels, the authors demonstrate the implementation by calculating IceCube's expected signal for 20 supernova progenitor masses under normal and inverted mass hierarchies and two equations of state.","pith_inferences":["Inference beyond the paper: The paper does not test Eq. (3.2) against a full IceCube event simulation, so the most useful next step would be a direct comparison of SNOwGLoBES output with published IceCube supernova rates under identical fluxes and deadtime assumptions.","Inference beyond the paper: Because $\\varepsilon_{\\rm noise} = 0.95$ is constant while real deadtime losses depend on the instantaneous DOM hit rate, the implementation likely overestimates the signal during the early, high-rate phase of a supernova; a time-dependent noise model would require extending SNOwGLoBES itself.","Inference beyond the paper: If merged into the standard SNOwGLoBES distribution, IceCube would become a default reference point in future supernova neutrino phenomenology, making detector-complementarity comparisons easier but also inheriting the systematic uncertainty of the effective-volume approximation."],"forward_implications":["SNOwGLoBES users can now compute IceCube supernova neutrino event rates for arbitrary input flux spectra, not only precomputed models.","IceCube can be included alongside other detectors in SNOwGLoBES-based sensitivity comparisons for core-collapse supernovae.","A user-supplied noise-screening factor can be recovered by dividing the computed rates by the built-in constant 0.95.","The same effective-volume-to-efficiency procedure applies to other large water Cherenkov neutrino detectors, given an estimate of their effective volume."],"supporting_citations":[{"why":"Supplies the SNOwGLoBES package and its pre-packaged water Cherenkov channel and cross-section files that the implementation reuses.","marker":"[1]"},{"why":"Provides the IceCube supernova event-rate formula, effective-volume plots, deadtime treatment, and the 98% DOM uptime assumption.","marker":"[3]"},{"why":"Provides IceCube DOM noise rates and supernova sensitivity context used for the noise treatment.","marker":"[4]"},{"why":"Gives the effective-volume parameterization used in Eq. (2.3), including the Cherenkov threshold, $C_{\\pm}$, $dN_\\gamma/dx$, and $\\langle V_{\\rm eff,\\gamma}\\rangle$.","marker":"[7]"},{"why":"Supplies the inverse-beta-decay cross section used for the main IBD channel.","marker":"[9]"},{"why":"Supplies the neutrino-electron scattering cross sections used for the scattering channels.","marker":"[10]"}],"fun_headline_variants":["IceCube joins SNOwGLoBES for supernova neutrino rates","Supernova neutrino rates for IceCube now in SNOwGLoBES","SNOwGLoBES gets IceCube simulation files","IceCube supernova sensitivity now in standard package","SNOwGLoBES implementation computes IceCube supernova signals"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that IceCube's position- and energy-dependent response to supernova neutrinos can be represented in SNOwGLoBES by a fixed mass of 51,600 kilotonnes, an energy-dependent post-smearing efficiency derived from the effective volume, and a constant 0.95 deadtime factor, without biasing the calculated event rates.","fun_headline_variants_meta":{"raw":{"variants":["IceCube joins SNOwGLoBES for supernova neutrino rates","Supernova neutrino rates for IceCube now in SNOwGLoBES","SNOwGLoBES gets IceCube simulation files","IceCube supernova sensitivity now in standard package","SNOwGLoBES implementation computes IceCube supernova signals"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000159,"raw_usage":{"total_tokens":1191,"prompt_tokens":872,"completion_tokens":319,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":488,"completion_tokens_details":{"reasoning_tokens":229}},"tokens_in":488,"tokens_out":319,"duration_ms":3490,"temperature":1.0,"reasoning_tokens":229,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T05:32:40.206038+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the same supernova flux used in the paper's demonstration and run a full Monte Carlo simulation of IceCube's DOM response (or use the published IceCube supernova analysis) to produce a hit-rate time series; compare it to the SNOwGLoBES calculation for the same input. Agreement within the expected statistical and systematic errors would support the implementation; a systematic divergence, especially at low neutrino energies or early times, would show that the constant-mass, effective-volume mapping in Eq. (3.2) is not adequate.","supporting_citations":[{"cited_title":"Albert et al., 2018","cited_arxiv_id":null,"evidence_quote":"Supplies the SNOwGLoBES package and its pre-packaged water Cherenkov channel and cross-section files that the implementation reuses."},{"cited_title":"Abbasi et al","cited_arxiv_id":null,"evidence_quote":"Provides the IceCube supernova event-rate formula, effective-volume plots, deadtime treatment, and the 98% DOM uptime assumption."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides IceCube DOM noise rates and supernova sensitivity context used for the noise treatment."},{"cited_title":"Kowarik, 2010","cited_arxiv_id":null,"evidence_quote":"Gives the effective-volume parameterization used in Eq. (2.3), including the Cherenkov threshold, $C_{\\pm}$, $dN_\\gamma/dx$, and $\\langle V_{\\rm eff,\\gamma}\\rangle$."},{"cited_title":"Strumia & F","cited_arxiv_id":null,"evidence_quote":"Supplies the inverse-beta-decay cross section used for the main IBD channel."},{"cited_title":"Marciano & Z","cited_arxiv_id":null,"evidence_quote":"Supplies the neutrino-electron scattering cross sections used for the scattering channels."}],"review_version":1}