Results From Core-Collapse Simulations with Multi-Dimensional, Multi-Angle Neutrino Transport

We present new results from the only 2D multi-group, multi-angle calculations of core-collapse supernova evolution. The first set of results from these calculations was published in Ott et al. (2008). We have followed a nonrotating and a rapidly rotating 20 solar mass model for ~400 ms after bounce. We show that the radiation fields vary much less with angle than the matter quantities in the region of net neutrino heating. This obtains because most neutrinos are emitted from inner radiative regions and because the specific intensity is an integral over sources from many angles at depth. The latter effect can only be captured by multi-angle transport. We then compute the phase relationship between dipolar oscillations in the shock radius and in matter and radiation quantities throughout the postshock region. We demonstrate a connection between variations in neutrino flux and the hydrodynamical shock oscillations, and use a variant of the Rayleigh test to estimate the detectability of these neutrino fluctuations in IceCube and Super-K. Neglecting flavor oscillations, fluctuations in our nonrotating model would be detectable to ~10 kpc in IceCube, and a detailed power spectrum could be measured out to ~5 kpc. These distances are considerably lower in our rapidly rotating model or with significant flavor oscillations. Finally, we measure the impact of rapid rotation on detectable neutrino signals. Our rapidly rotating model has strong, species-dependent asymmetries in both its peak neutrino flux and its light curves. The peak flux and decline rate show pole-equator ratios of up to ~3 and ~2, respectively.

💡 Research Summary

This paper presents the most advanced two‑dimensional core‑collapse supernova simulations to date that incorporate full multi‑group, multi‑angle neutrino transport. Building on the initial results reported by Ott et al. (2008), the authors follow two 20‑solar‑mass progenitors—a non‑rotating model and a rapidly rotating model—for approximately 400 ms after core bounce. The key methodological advance is the solution of the Boltzmann transport equation for neutrinos with explicit angular resolution, allowing the radiation field to be constructed as an integral over many deep sources rather than relying on diffusion‑type approximations such as multi‑group flux‑limited diffusion (MGFLD).

The first major finding is that, in the net heating region behind the stalled shock, the neutrino radiation field varies far less with angle than the hydrodynamic quantities (density, temperature, velocity). This occurs because the bulk of the neutrino luminosity originates deep inside the proto‑neutron star where conditions are nearly spherical, and the specific intensity at a given point is an angular average over many interior emission directions. Such “multi‑source averaging” can only be captured with true multi‑angle transport; simpler schemes systematically underestimate angular uniformity.

The second result concerns the phase relationship between dipolar oscillations of the shock radius and fluctuations in matter and radiation variables throughout the post‑shock region. The authors show that while density and temperature oscillations are nearly in phase with the shock motion, the neutrino flux exhibits a roughly 90‑degree phase lag: the flux peaks when the shock is at its minimum radius and reaches a trough when the shock expands maximally. This lag reflects the finite response time of the radiation field to changes in the underlying matter distribution and is quantified using a variant of the Rayleigh statistical test.

The third component of the study translates these intrinsic flux variations into observable signals for existing detectors. Assuming no flavor conversion, the non‑rotating model’s flux fluctuations would be detectable by IceCube out to roughly 10 kpc and by Super‑Kamiokande out to about 5 kpc, with a detailed power‑spectral analysis possible within the nearer range. In contrast, the rapidly rotating model displays strong, species‑dependent asymmetries: the neutrino flux is up to three times larger along the rotation axis than in the equatorial plane, and the decline rates differ by a factor of two. Consequently, the effective detection horizon for the rotating case is substantially reduced, and any realistic inclusion of flavor oscillations would further diminish the observable signal.

Finally, the authors quantify how rapid rotation imprints itself on the neutrino signal. The rotating core develops an oblate shape, leading to pole‑equator differences in the peak ν_e and ν̄_e fluxes of up to a factor of three, and similar anisotropies in the subsequent light curves. These anisotropies imply that an observer’s viewing angle could dramatically alter the inferred neutrino luminosity and temporal evolution, complicating the extraction of intrinsic explosion properties from a single detector.

In summary, the paper demonstrates that multi‑angle neutrino transport is essential for accurately capturing the subtle coupling between hydrodynamic shock oscillations and neutrino emission in core‑collapse supernovae. It also shows that rapid rotation introduces pronounced directional dependence in the neutrino signal, thereby affecting the prospects for detection with current facilities. The work provides a critical benchmark for future theoretical models and for the interpretation of the next Galactic supernova neutrino burst.