New Signs Pointing Toward a Correlation Between Astrophysical Neutrinos and Radio Flares

arXiv:2511.21627 · astro-ph.GA, astro-ph.HE · Submitted 2025-11-26 · Read on arXiv

Listen

Transcript

Introduction to the show: ident: Astrophysics Radio. Generated commentary on the latest astrophysics papers.

Vera: Today's paper: "New Signs Pointing Toward a Correlation Between Astrophysical Neutrinos and Radio Flares".

Jocelyn: The gist: We find an excess number of associations between flaring radio sources and neutrinos that were detected between the first and second VLASS observations at > 2σ confidence,

Vera: First, who's behind it and why it matters.

Title and authors: Vera: To unpack that a bit more, Jocelyn, what does this excess mean in simple terms for us here listening right now? Is it just a statistical fluke or something more concrete?

Jocelyn: Well, they are comparing the actual count of these associations to what you’d expect if you just randomly matched two different populations together. The real number of associations they found is higher than that random background suggests, which is what gives them that two sigma confidence.

Subrahmanyan: From a theoretical side, that excess suggests that radio flares might be contributing about thirteen percent of all the astrophysical neutrinos we've seen so far in the IceCube data <ref:2511.21627#pg1>. That’s a specific contribution number they are working with.

Vera: Thirteen percent of the total neutrino count is significant when you consider how rare these high-energy neutrino events are in general. It gives us a concrete fraction to work with when modeling where these neutrinos might be coming from.

Jocelyn: And they also found that those radio flares that actually match up with neutrinos are usually located closer on the sky to the neutrino event than if they were randomly placed, meaning there’s a spatial bias.

Subrahmanyan: That spatial bias is an important piece of information because it suggests the source isn't just a random coincidence; it implies some kind of physical connection between the flare and the neutrino production site.

The paper's summary: Vera: So, to summarize what they actually did with this paper, they took data from VLASS and IceCube to search for this link between radio variability and neutrinos, and their main finding is that there is an excess of these matches at over two sigma confidence.

Jocelyn: And beyond that excess, the paper looks at the properties of those associated radio flares. They found that out of the sixty-six flaring radio sources they studied, only nine, which is about fourteen percent, had an X-ray counterpart in MORX.

Subrahmanyan: That lack of a strong high-energy electromagnetic counterpart is something they discuss. They suggest it might be because the obscuration is too thick or perhaps because their radio observations are sampling a larger volume than where the gamma-ray and X-ray data are coming from.

Vera: That makes sense in terms of how different wavelengths interact with the same source, but they also looked at the infrared colors of these radio flares, and those didn't look much different from what you’d expect for any random radio variable source.

Jocelyn: And interestingly, when they look at blazars specifically within their samples of radio flares associated with neutrinos, about fifty percent of them are blazars. That’s pretty similar to other findings that suggest a small fraction of astrophysical neutrinos originate from blazars in the first place.

Subrahmanyan: So, the summary is really tying together the detection rate, the spatial bias, and then discussing why we don't always see an immediate high-energy electromagnetic signal when we look at these neutrino events.

The paper's improvements: Vera: Now they discuss how this work could be improved or what future steps should be taken based on their results. They suggest that incorporating time-domain information into the search for electromagnetic counterparts can really cut down on the uncertainty when matching individual neutrino events to potential sources.

Jocelyn: That’s smart because it means they can specifically look for things like, "was the neutrino detected in a time window defined by the observation dates of the first and second VLASS epochs?" that tightens up the search significantly.

Subrahmanyan: They also talk about developing a better way to classify these sources based on their multiwavelength properties. They point out that since over eighty percent of the radio sources associated with IceCat-one events don't show up in gamma-rays or X-rays, we need to use radio surveys more heavily for finding these counterparts <ref:2511.21627#pg1>.

Vera: And they predict that if their results hold up, the significance of this excess will get better by the time VLASS Epoch four is done, and they expect it to go past three sigma confidence.

Jocelyn: If we get improved uncertainties in the sky position of those neutrino events with IceCat-two data, they think that could push that significance even higher, possibly to over five sigma confidence when you look at the excess number of associations between radio flares and high energy astrophysical neutrinos <ref:2511.21627#pg1,excess number of associations between>.

Conclusion: Vera: So to wrap up this discussion on "New Signs Pointing Toward a Correlation Between Astrophysical Neutrinos and Radio Flares," the main thing is that they’ve found an excess number of associations between flaring radio sources and neutrinos that are both spatially and temporally associated at over two sigma confidence.

Jocelyn: And they estimate that this excess suggests radio flares might be contributing about thirteen percent of the astrophysical neutrinos observed by IceCube, which aligns with some previous studies <ref:2511.21627#pg1,of the astrophysical neutrinos observed by IceCube>. They also found a spatial bias where the real associations have a median angular separation of one hundred forty arcminutes compared to random associations at one hundred eighty arcminutes.

Subrahmanyan: The implication for us is that we need to keep using radio time-domain surveys because they seem crucial for finding electromagnetic counterparts to astrophysical neutrinos. It suggests that the production mechanism might involve radio flares as a significant part of how these high-energy neutrinos are generated.

Vera: So, the paper "New Signs Pointing Toward a Correlation Between Astrophysical Neutrinos and Radio Flares" shows that there is evidence linking radio variability and neutrino detections through spatial and temporal coincidence. It’s a solid piece of observational evidence we can build on as we look deeper into these high-energy phenomena.

Jocelyn: We’ll keep an eye on those future predictions they made about the significance improving with more data from IceCat-two to see if that thirteen percent contribution becomes clearer <ref:2511.21627#pg1>.

Subrahmanyan: It’s a compelling piece of work that pushes the connection between high-energy neutrinos and observable radio phenomena further into the realm of possibility.

Department of Physics, University of Wisconsin-Madison · Department of Astronomy, University of Washington

astro-ph.GA, astro-ph.HE

Submitted: 2025-11-26

Updated: 2026-10-08

Comments: 15 pages, 11 figures, 5 tables. Published in OJAp. Supplementary data available at https://zenodo.org/records/23223067

License: http://creativecommons.org/licenses/by/4.0/

Importance score: 76/100

The gist: The gist: We find an excess number of associations between flaring radio sources and neutrinos that were detected between the first and second VLASS observations at > 2σ confidence, which is

Key concepts

Astrophysical Neutrinos: These are subatomic particles produced in extreme cosmic events, like those from active galactic nuclei. They are unique because they interact very weakly with matter, allowing them to travel across vast cosmic distances and penetrate dense environments that block light.
VLASS Survey: The Very Large Array Sky Survey is a multi-epoch radio survey mapping the sky north of $\delta = -40^°$ at 3 GHz. It tracks variable radio sources over time, providing data on how these sources change their brightness across different observation periods.
Association Analysis: This method compares the actual number of detected neutrino events to the expected random background using Monte Carlo simulations. The study specifically looked for both spatial proximity and temporal coincidence between radio flares and neutrino detections to establish a correlation.
HE Neutrinos: High-energy neutrinos are those with energies greater than 1 TeV. These particles are crucial because their detection confirms the existence of high-energy astrophysical processes in the Universe, such as those occurring near blazars or other powerful cosmic accelerators.

Terminology

AGN — 11× in this paper · explained in 1 episode
Very Large Array Sky Survey — 5× in this paper
Active Galactic Nuclei (AGN) — 2× in this paper · explained in 1 episode
Adam — 2× in this paper · explained in 2 episodes
Clustering — 1× in this paper · explained in 2 episodes
Selection Bias — 1× in this paper · explained in 4 episodes
Sparsity — 1× in this paper · explained in 4 episodes
Trace — 1× in this paper · explained in 2 episodes

Summary

The gist: We find an excess number of associations between flaring radio sources and neutrinos that were detected between the first and second VLASS observations at > 2σ confidence, which is consistent with radio flares contributing ∼ 13 % of the astrophysical neutrinos observed by IceCube.

Introduction and Motivation

Astrophysical neutrinos present a unique way to study the cosmos because their low interaction cross section allows them to pass through highly dense environments that are opaque to electromagnetic radiation (particularly at shorter wavelengths) Moreover, the extremely low probability of interaction means that neutrinos can reach Earth from the farthest and earliest reaches of the Universe—unlike light, neutrinos have no effective horizon. The IceCube Neutrino Observatory has confirmed the existence of high energy (HE; Eν ≳ 1 TeV) astrophysical neutrinos (IceCube Collaboration 2013). Some HE neutrinos are diffuse emission from the plane of the Milky Way (Icecube Collaboration et al. 2023), potentially resulting from cosmic rays in the Milky Way’s interstellar medium. Two extragalactic point sources of HE neutrinos have been identified so far: the blazar TXS 0506 + 065, and the Seyfert galaxy NGC 1068 (IceCube Collaboration et al. 2018, 2022). Active galactic Nuclei (AGN) such as TXS 0506 + 065 and NGC 1068 are intriguing candidates as neutrino sources, with their central engine’s potentially hosting a number of mechanisms that should in theory produce HE neutrinos, e.g., magnetic reconnection events in the accretion disk halo, or shocks in the jet or accretion disk (Mannheim 1995; Bednarek & Protheroe 1999; Khiali & de Gouveia Dal Pino 2016; Blandford et al. 2019; Murase & Stecker 2023). Identifying electromagnetic counterparts to neutrino detections has proven difficult because the typically large positional uncertainties for neutrino detections encompass many potential sources of origin.

Data and Methodology

The study utilizes data from the first two epochs of the Very Large Array Sky Survey (VLASS) and the IceCube Neutrino Observatory. VLASS is a multi-epoch survey of the entire sky north of δ = −40◦ at ν ∼ 3 GHz that began in 2017, with a synthesized beam size of ∼ 3′′ and a typical rms noise of 130 µJy beam−1 in each epoch. The analysis focuses on identifying associations between variable radio sources detected in VLASS and HE neutrinos detected by IceCube. The researchers define three samples of VLASS variables with which to search for a correlation with HE neutrinos. One sample is the full sample of 1, 928 variables, for which they search for a spatial association only. A second sample is the subset of 1, 093 radio variables that brightened between VLASS Epoch 1 and Epoch 2. The third sample consists of the same 1, 093 radio sources as in the flaring sample, but for the temporal correlation they shift the time window in which they allow the neutrino to be associated by adding 150 days to the neutrino’s event time in IceCat-1.

Cross Association Analysis

The researchers match all three samples of VLASS variables to the IceCat-1 data using positional error ellipses defined by the uncertainties of both sources. They define a VLASS source as temporally associated with a neutrino if the neutrino was detected in a time window defined by the observation dates of the first and second VLASS epochs. The analysis compares observed numbers of associations to expected random background contributions using 1,000 iteration Monte Carlo simulations for each test sample. For example, in the flaring sample, the real count of 64 associations has a p-value of 0.034 based on the simulated distribution. The significance of the excess number of associations is clearly driven by the number of radio sources rather than the neutrino counts.

Multiwavelength Properties and Interpretation

The study investigates multiwavelength properties to see if radio flares are associated with high energy electromagnetic counterparts. Of the 66 flaring radio sources, just 9 (14 %) have an X-ray counterpart in MORX. The authors suggest that the lack of high energy electromagnetic counterparts might be due to obscuration or that radio observations are probing a larger cosmic volume than available γ-ray and X-ray data. Furthermore, the IR colors for radio flares associated with IceCat-1 events are not significantly different from what would be expected for any radio variable source. The percentage of blazars in their samples of radio flares associated with neutrinos is ∼ 50 %, which is similar to findings that suggest ∼ 6 − 7 % of astrophysical neutrinos originate from blazars.

Conclusion and Future Prospects

The key findings are that they find an excess number of associations between flaring radio sources and neutrinos that are both spatially and temporally associated when compared to the expected number of random associations, at > 2σ confidence. The majority (> 80 %) of radio sources associated with IceCat-1 events are not detected at γ-rays or X-rays, highlighting the importance of radio timedomain surveys in identifying candidate electromagnetic counterparts to astrophysical neutrinos. Assuming their results are representative, they predict that the significance of the excess number of associations between radio flares and neutrinos will exceed 3σ confidence by the end of VLASS Epoch 4. Should improved uncertainties in the sky position of neutrino events in IceCat-2 also lower the shot noise as expected, they anticipate > 5σ confidence in the excess number of associations between radio flares and HE astrophysical neutrinos. The analysis suggests that a contribution of ∼ 13 % to the astrophysical neutrino counts by flaring radio sources seems realistic to first order.

Summary and Conclusions

The study investigates whether variable radio sources contribute to the astrophysical neutrinos observed by IceCube. Their key findings are summarized below.

"We find an excess number of associations between flaring radio sources and neutrinos that are both spatially and temporally associated when compared to the expected number of random associations, at > 2σ confidence"

"Radio flares that are spatially and temporally associated with neutrinos are typically located closer on the sky to the neutrino event than simulated random associations. Real associations have a median angular separation of 140′ compared to 180′ for random associations."

"The majority (> 80 %) of radio sources associated with IceCat-1 events are not detected at γ-rays or X-rays, highlighting the importance of radio timedomain surveys in identifying candidate electromagnetic counterparts to astrophysical neutrinos."

The excess number of neutrinos identified by these associations is consistent with flaring radio sources contributing ∼ 13 % of the astrophysical neutrinos observed by IceCube, in agreement with previous studies.

"The number of associations between variable radio sources and neutrinos is consistent with the expected background when no time correlation is required, i.e. there is no evidence of a spatial-only correlation between radio variables and neutrinos."

This paper was built using the Open Journal of Astrophysics LATEX template. The OJA is a journal which provides fast and easy peer review for new papers in the astro-ph section of the arXiv, making the reviewing process simpler for authors and referees alike. The VLA is operated by NRAO, a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. WISE is a joint project of the University of California, Los Angeles, and the Jet Propulsion Laboratory/California Institute of Technology, and NEOWISE, which is a project of the Jet Propulsion Laboratory/California Institute of Technology. The KM3NeT Collaboration, MessMapp Group, Fermi-LAT Collaboration, et al. 2025. The VLA is operated by NRAO, a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. WISE is a joint project of the University of California, Los Angeles, and the Jet Propulsion Laboratory/California Institute of Technology, and NEOWISE, which is a project of the Jet Propulsion Laboratory/California Institute of Technology The KM3NeT Collaboration, MessMapp Group, Fermi-LAT Collaboration, et al. 2025 WISE and NEOWISE are funded by the National Aeronautics and Space Administration. The VLA is operated by NRAO, a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc.

Improvements for AI systems

A model capable of performing cross-association stacking experiments could be improved by incorporating time-domain constraints, as incorporating time-domain information into searches for the electromagnetic counterparts to HE neutrinos can reduce uncertainty in the associations between individual neutrino events and potential sources of origin. This allows the system to search for correlations between radio variability and neutrinos, specifically identifying candidates where the neutrino was detected in a time window defined by the observation dates of the first and second VLASS epochs.
A system could be enhanced to predict future detection significance based on updated instrumentation, as "We can estimate the significance for a given number of associations between neutrinos and radio flares, Nassoc = 62 compared to an expected λ = 30 using the IceCat-2 data, an observation that would have p ∼ 10−7, or > 5σ confidence." This enables proactive assessment of detection capabilities as new catalogs like IceCat-2 are released.
A sophisticated classification module could be developed to assess source type likelihood based on multiwavelength properties, leveraging the finding that "the majority (> 80 %) of radio sources associated with IceCat-1 events are not detected at γ-rays or X-rays, highlighting the importance of radio observations for identifying potential electromagnetic counterparts to astrophysical neutrinos." This allows the system to weigh the reliability of a potential counterpart based on its absence or presence in high-energy emission channels.
A search algorithm could be refined to prioritize spatial and temporal coincidence simultaneously, as the paper demonstrates that Real associations have a median angular separation of 140′ compared to 180′ for random associations, suggesting a bias toward closer associations. This refinement would improve the identification of true physical counterparts by filtering out purely spatial matches.

Abstract

We use data from the first two epochs of the Very Large Array Sky Survey (VLASS) and the IceCube Neutrino Observatory to search for evidence of a correlation between radio variability and the detection of astrophysical neutrinos. Our work probes the radio population down to S 3, GHz a few mJy, more than an order of magnitude fainter than previous analyses that aim to cross correlate radio variability with neutrino detections. We find an excess number of associations between flaring radio sources and neutrinos that were detected between the first and second VLASS observations at >2σ confidence. This excess is consistent with radio flares contributing about13,% of the astrophysical neutrinos observed by IceCube. Notably, about90,% of the radio flares associated with neutrinos are not detected at either γ-ray or X-ray wavelengths, highlighting the importance of deep radio observations for identifying potential electromagnetic counterparts to astrophysical neutrinos. No excess in the number of associations between the wider radio-variable population and the IceCube neutrinos is seen when no time constraint is placed on the neutrino detection. We predict that data from future VLASS epochs will see an excess number of associations between radio flares and neutrinos at the >3σ level, and expected improvements to the positional constraints on the neutrinos may increase that confidence to >5σ, should our results be representative.