Double Cascade Signature in Neutrino Detectors

• Double cascade signature is defined by two sequential cascades from tau-neutrino interactions that enable clear flavor discrimination in neutrino telescopes.
• Reconstruction methods, including Taupede fitting and waveform double-pulse analysis, enhance event resolution and improve tau neutrino identification.
• Detector geometry and energy thresholds, as seen in IceCube and KM3NeT, play a critical role in resolving cascade separations above ~0.2 PeV.

“Double cascade signature” is a term used in several research literatures for observables generated by sequential two-stage processes. In current arXiv usage, its dominant meaning is the topology produced by a charged-current $\nu_\tau$ interaction in a large-volume optical Cherenkov detector: a first compact cascade at the interaction vertex, followed—after propagation of the $\tau$ lepton—by a second cascade from $\tau$ decay. This topology is central to $\nu_\tau$ identification in IceCube, Baikal-GVD, and KM3NeT/ARCA, where it helps break the $\nu_e$–$\nu_\tau$ degeneracy of single-cascade samples and constrains the astrophysical flavor composition at Earth (Chen et al., 11 Jul 2025, Stachurska, 2019, Eeden et al., 2022). The same phrase also appears in distinct contexts, including moiré interlayer excitons, collider cascades in the type-II seesaw model, and coupled contagion processes in financial networks (Tan et al., 2022, Bai et al., 2021, Hurd et al., 2013).

1. Physical origin in neutrino telescopes

In neutrino astronomy, the double cascade is the hallmark of a charged-current tau-neutrino interaction. The underlying chain is

$\nu_\tau + N \to \tau + X,$

followed by $\tau$ decay. The first interaction produces a nearly point-like hadronic or electromagnetic shower at the interaction vertex. If the $\tau$ decays hadronically or electronically, together about $83\%$ of the time, the decay yields a second compact cascade; if it decays muonically, about $\tau$0 of the time, the event contains a short track rather than a second cascade and is not targeted by cascade-based selections (Chen et al., 11 Jul 2025).

IceCube descriptions emphasize the topology as two time-ordered energy depositions connected by a faint, short $\tau$1 track. In that terminology, the first cascade is hadronic at the interaction vertex, while the second arises from $\tau$2 decay products. This is the direct event class that gives IceCube sensitivity to $\tau$3 at high energy and breaks the flavor degeneracy that persists when only single cascades are observed (Stachurska, 2019).

The spatial separation is controlled by the $\tau$4 decay length,

$\tau$5

with $\tau$6 and $\tau$7. Numerically,

$\tau$8

or, equivalently in the KM3NeT formulation, $\tau$9 (Stachurska, 2019, Eeden et al., 2022). Representative separations are therefore about $\tau$0 at $\tau$1, about $\tau$2 at $\tau$3, and about $\tau$4 at $\tau$5 (Chen et al., 11 Jul 2025).

2. Resolvability and detector-scale constraints

Whether a double cascade is observable depends on the relation between $\tau$6 and detector granularity, containment, and optical transport. In IceCube, the Digital Optical Modules are spaced by about $\tau$7 vertically and about $\tau$8 horizontally. The 11-year cascade analysis states that the topology becomes resolvable once $\tau$9 exceeds roughly $\nu_\tau$0, which sets an effective energy scale of order a few $\nu_\tau$1 for clearly separated double cascades, although identification can sometimes occur at somewhat lower energies (Chen et al., 11 Jul 2025).

Baikal-GVD presents a similar geometric threshold. With $\nu_\tau$2 vertical optical-module spacing and about $\nu_\tau$3 central-to-peripheral string spacing within a cluster, practical minimum separations are also stated to be of order $\nu_\tau$4. At lower separations, Baikal-GVD supplements geometric separation with a waveform-level double-pulse method. In synthetic Gumbel-modeled waveforms, a local minimum between two peaks requires a time delay of about $\nu_\tau$5, corresponding to a light-path difference of order $\nu_\tau$6–$\nu_\tau$7 at the sensor (Allakhverdyan et al., 2021).

KM3NeT/ARCA frames the same issue through reconstruction performance rather than a hard geometric threshold. For selected double-cascade events with reconstructed energy above $\nu_\tau$8 and a contained vertex, the single-cascade median angular deviation remains at about $\nu_\tau$9, whereas the double-cascade reconstruction drops below $\nu_e$0 for $\nu_e$1 lengths larger than $\nu_e$2 and reaches $\nu_e$3 at $\nu_e$4 (Eeden et al., 2022).

Detector	Resolvable regime	Representative metric
IceCube	Minimum reconstructed decay length $\nu_e$5	Double cascades resolvable once $\nu_e$6 exceeds roughly $\nu_e$7
Baikal-GVD	Geometric separation $\nu_e$8; double-pulse for smaller $\nu_e$9	Minimal double-pulse split about $\nu_\tau$0
KM3NeT/ARCA	Selected events above $\nu_\tau$1; performance improves with $\nu_\tau$2 length	Sub-degree angular deviation for lengths $\nu_\tau$3

These thresholds are not identical selection rules across experiments, but they establish a common operational regime: double-cascade observability improves rapidly once the two light sources are separated at the scale of several to several tens of meters.

3. Reconstruction methodologies and discriminating observables

IceCube’s current cascade-based search reconstructs events under single- and double-cascade hypotheses with Taupede, which fits the positions and energies of two spatially separated showers. Two topological observables are emphasized: the reconstructed $\nu_\tau$4 decay length and the energy asymmetry

$\nu_\tau$5

Preselection requires reconstructed energy above $\nu_\tau$6, reconstructed decay length $\nu_\tau$7, and containment of both cascades within the instrumented volume and outside the optically adverse dust layer. Two Boosted Decision Trees are then applied: the first suppresses single-cascade backgrounds from $\nu_\tau$8 charged-current and neutral-current interactions of all flavors, and the second removes starting-track backgrounds from $\nu_\tau$9 charged-current events. The most important inputs include the single- and double-cascade log-likelihoods, reconstructed decay length, energy asymmetry, individual cascade energies, and total recorded charge. With optimized selections, the sample reaches a signal-to-background ratio of about $\nu_\tau + N \to \tau + X,$0, $\nu_\tau + N \to \tau + X,$1-neutrino purity of about $\nu_\tau + N \to \tau + X,$2, a weighted mean reconstruction error on the $\nu_\tau + N \to \tau + X,$3 decay length of about $\nu_\tau + N \to \tau + X,$4, and an expected selected rate of about $\nu_\tau + N \to \tau + X,$5 $\nu_\tau + N \to \tau + X,$6 double-cascade events per year (Chen et al., 11 Jul 2025).

IceCube’s earlier topology-based analysis reconstructed all events under single-cascade, track, and double-cascade hypotheses and used three observables for classification: double-cascade length, energy asymmetry, and energy confinement. Acceptance of a double-cascade reconstruction required convergence, at least $\nu_\tau + N \to \tau + X,$7 reconstructed energy in each cascade, an opening angle within $\nu_\tau + N \to \tau + X,$8 between the best-fit double-cascade and track hypotheses, and soft containment such that neither cascade extended more than $\nu_\tau + N \to \tau + X,$9 outside the instrumented volume. The explicit selection cuts were length $\tau$0, energy asymmetry $\tau$1, and energy confinement $\tau$2–$\tau$3 (Stachurska, 2019).

Baikal-GVD has developed two complementary approaches. The waveform-level method identifies double pulses by a sign change in the first derivative, then fits each pulse with a Gumbel function,

$\tau$4

and suppresses pedestal-induced fakes with a ROOT TMVA BDT trained on 12 waveform parameters. At BDT cut $\tau$5, the reported signal efficiency is $\tau$6 and the background efficiency is $\tau$7 (Allakhverdyan et al., 2021).

For geometrically resolved events, Baikal-GVD uses causality-based hit selection, many five-hit space-time seeds, splitting of hits into two cascade subsets, a timing $\tau$8 fit for positions and times, and a Poisson likelihood for charges including non-detections. In the 2023 formulation, the algorithm reconstructs positions and times of both cascade vertices, the direction defined by the line joining them, and the two cascade energies. On simulated $\tau$9 events with true vertex separation above $\tau$0 and $\tau$1 energy above $\tau$2, the reported preliminary precisions are mean position errors of $\tau$3 for cascade A and $\tau$4 for cascade B, and mean $\tau$5 (Aynutdinov et al., 2023).

KM3NeT/ARCA uses a three-stage maximum-likelihood chain: a single-cascade prefit with Aashowerfit, a $\tau$6-length prefit along the prefit direction, and a full two-cascade fit over first hits on each PMT. The fitted parameter vector includes the interaction vertex, direction, $\tau$7 length, and the energy-asymmetry parameter. The paper reports a median $\tau$8-length error of $\tau$9 with $83\%$0 quantiles of $83\%$1, while the summary quotes a $83\%$2-length resolution of $83\%$3. The visible-energy reconstruction has median error $83\%$4 with $83\%$5 quantiles of $83\%$6, while the summary quotes an energy resolution of $83\%$7 (Eeden et al., 2022).

4. Scientific role in IceCube: candidate events, flux measurement, and flavor composition

The first explicit IceCube double-cascade candidates were reported in the 7.5-year High-Energy Starting Events sample. Above $83\%$8 reconstructed deposited energy, that sample contained 60 events: 42 single cascades, 16 tracks, and 2 double cascades. The two candidates, “Big Bird” and “Double Double,” had reconstructed lengths of $83\%$9 and $\tau$00, respectively. Their a-posteriori $\tau$01 charged-current probabilities were about $\tau$02 for Big Bird and at least $\tau$03 for Double Double. In the same topology-separated likelihood fit, the best-fit Earthly flavor composition was $\tau$04, consistent with the standard $\tau$05 expectation and also statistically consistent with zero astrophysical $\tau$06 within uncertainties (Stachurska, 2019).

The 11-year IceCube cascade analysis generalizes the role of the double-cascade signature from candidate identification to joint flux-and-flavor inference. The cascade sample used for diffuse-flux measurement contains over 14,000 events across 2010–2020. For a single power law, the per-flavor flux is parameterized as

$\tau$07

with best-fit values $\tau$08 and $\tau$09. A broken power law yields $\tau$10, $\tau$11, $\tau$12, and $\tau$13. The corresponding binned goodness-of-fit values are $\tau$14 and $\tau$15, indicating a markedly better fit for the broken-power-law hypothesis in this cascade channel (Chen et al., 11 Jul 2025).

Flavor composition is then constrained with a joint binned maximum-likelihood fit combining three statistically independent samples: the 11-year $\tau$16-enriched double-cascade subset, the remaining 11-year single-cascade events, and a 9.5-year northern-sky track sample dominated by $\tau$17 charged-current interactions. The flavor fractions at Earth are parameterized as $\tau$18 with $\tau$19. Sensitivity is evaluated with an Asimov data set under $\tau$20, and the inclusion of the $\tau$21-enriched sample significantly tightens sensitivity to $\tau$22 relative to fits using only cascades and tracks; the reported $\tau$23 and $\tau$24 contours, under both single- and broken-power-law assumptions, encompass the canonical $\tau$25 point (Chen et al., 11 Jul 2025).

5. Backgrounds, uncertainties, and limitations

The principal challenge is that true $\tau$26 double cascades are rare and can be mimicked by other topologies. In the 11-year IceCube search, the dominant residual backgrounds are single cascades from $\tau$27 charged-current and neutral-current interactions of all flavors, and starting $\tau$28 charged-current tracks with large stochastic energy losses. The multistage BDT chain is designed to suppress these while retaining efficiency for $\tau$29 signal events (Chen et al., 11 Jul 2025).

The 2019 IceCube analysis described the same problem in event-topology terms. True single cascades, especially $\tau$30 charged-current events, can fluctuate into two nearby depositions with high energy confinement. $\tau$31 charged-current events with a large localized stochastic energy loss can resemble two “blobs,” and atmospheric muons or bundles with complex light patterns can be misreconstructed as contained double cascades. The analysis therefore used containment, fit-quality, asymmetry, and energy-confinement cuts, with a length resolution of about $\tau$32 (Stachurska, 2019).

Detector and medium modeling are a second limiting factor. In IceCube, the 11-year cascade analysis uses updated ice modeling, specifically SPICE-3.2.1 with a unified two-parameter hole-ice description, plus pass-2 reprocessing with recalibrated DOM efficiencies. The analysis reports improved data–MC agreement, with a binned $\tau$33 p-value of $\tau$34 for the inherited cascade selection. Systematic uncertainties are dominated by bulk ice and refrozen hole ice, DOM efficiency and calibration, atmospheric self-veto modeling, and Monte Carlo statistics propagated through an effective likelihood. Relative to a pure Poisson likelihood, the effective likelihood broadens the $\tau$35 contours on normalization and spectral index by about $\tau$36 without shifting the best-fit point (Chen et al., 11 Jul 2025).

Geometry and containment impose hard phenomenological limits. IceCube explicitly imposes $\tau$37, which pushes the selected sample toward higher energies; the analysis notes that typical resolvable separations correspond to $\tau$38. Very long separations become difficult because containment efficiency falls once one or both cascades approach or exceed the instrumented volume, with an effective length cutoff around $\tau$39–$\tau$40 in the earlier HESE analysis (Stachurska, 2019). The dust layer in IceCube is excluded to protect reconstruction quality, at the cost of some acceptance (Chen et al., 11 Jul 2025).

Baikal-GVD and KM3NeT/ARCA are at earlier stages in this respect. Baikal-GVD explicitly notes that detailed background rates and systematics are not yet provided in the relevant reconstruction papers, and that extension to separations below $\tau$41 is under development (Aynutdinov et al., 2023). KM3NeT/ARCA likewise states that event selection and background response are not yet covered and will be the topic of future efforts (Eeden et al., 2022).

6. Other domain-specific meanings of the term

In moiré interlayer excitons, “double cascade signature” refers not to spatially separated particle showers but to sequential interlevel transitions among moiré-quantized exciton levels. Tan et al. describe a two-step cascade $\tau$42 identified by time-ordered photoluminescence dynamics: the delayed rise of $\tau$43 matches the fast decay constant of $\tau$44, and the delayed rise of $\tau$45 matches the fast decay constant of $\tau$46. The fitted seven-level rate model requires finite transition rates $\tau$47 and $\tau$48 (Tan et al., 2022).

In collider phenomenology, the phrase denotes a two-step decay chain within the non-degenerate Higgs Triplet Model. The relevant sequence is

$\tau$49

followed by $\tau$50 or $\tau$51. In the cited study, same-sign tetraleptons arise when both doubly charged Higgs bosons decay through these channels, and the rate is controlled by the near-degeneracy of $\tau$52 and $\tau$53, the parameter $\tau$54, and the mass-splitting parameter $\tau$55 (Bai et al., 2021).

In financial-network theory, “double cascade” denotes a coupled mapping between default contagion and funding-liquidity stress contagion. The model iterates these two channels to a fixed point and predicts, in the absence of asset fire sales, a negative relation between eventual defaults and the intensity of liquidity hoarding: stronger hoarding induces more stress but fewer defaults (Hurd et al., 2013).

These usages are not interchangeable. In neutrino telescopes, the term names a detector topology linked to $\tau$56 charged-current interactions; in excitonics it denotes sequential interlevel relaxation; in collider studies it labels a decay chain through intermediate resonances; and in systemic-risk modeling it refers to coupled propagation mechanisms. The shared feature is a resolved two-stage cascade process, but the underlying observables, reconstruction methods, and scientific objectives are domain-specific.