From Variability to SED Modeling: A Multiwavelength Study of the Neutrino Blazar TXS 0506+056

Abstract: The blazar TXS 0506+056 is the first source that was reported to be associated with high-energy extragalactic neutrino events and is one of the major targets for multi-messenger studies. We carried out multi-wavelength optical monitoring of this object on 24 nights in the period from 2018 to 2023. The overall light curves exhibit a dimming trend superposed by some small-amplitude fluctuations, and intraday variability was detected on four nights. Bluer-when-brighter behaviors were observed on both intraday and long timescales and were more pronounced on long timescales, while a weak redder-when-brighter trend was detected on one night. No significant time lags were found between variations at different optical wavelengths. We also retrieved the multi-broadband data from some monitoring programs. The data reveal complex, asynchronous flaring in different wavebands. A cross-correlation analysis shows that the high-energy emission (optical to gamma-ray) is co-spatial and leads the radio emission by a substantial time of about 800 to 900 days, suggesting that the radio emission originates from a downstream region of the jet. We performed time-dependent lepto-hadronic modeling of the spectral energy distributions for three representative epochs, the 2017 neutrino-associated flare, a post-flare phase, and a deep quiescent state, revealing an evolution in the radiative properties of the emission regions. The modeling results provide a phenomenological framework for interpreting the long-term multiwavelength behavior of TXS 0506+056 in a multi-messenger context.

• The paper combines five years of optical monitoring with roughly 15 years of multiwavelength data, detecting intraday variability on four nights with amplitudes up to 29.90% and a long-term bluer-when-brighter trend.
• The analysis finds no significant optical interband lags but measures approximately 800–900-day delays between high-energy and radio emission, supporting compact upstream and extended downstream jet zones.
• The paper’s SED models indicate that the neutrino-associated flare arose from a compact, strongly magnetized, proton-rich region, whereas later states were larger, lepton-dominated, and produced weaker predicted neutrino fluxes.

The blazar TXS 0506+056 ($z=0.3365$) is the first extragalactic source associated with high-energy neutrino events, and it remains a central target for multi-messenger studies. This paper by Du, Xia, Wu, and Fang presents a dedicated optical monitoring campaign of the source together with a long-term multi-broadband analysis and time-dependent lepto-hadronic spectral energy distribution (SED) modeling, aiming to connect the source's variability behavior to the physical conditions of its jet across distinct activity states (2608.17526).

Observations and data

The optical monitoring used the 85 cm telescope at Xinglong Station of NAOC, collecting roughly 3,000 data points on 24 nights between October 2018 and November 2023 in the $R$ and $I$ bands (2018) and subsequently in $B$, $V$, and $R$. Differential photometry relative to three comparison stars was performed with an aperture of $2\times$FWHM, and the full light curves are released in machine-readable form. These data were supplemented with archival photometry from AAVSO, ZTF, ASAS-SN, and KAIT spanning MJD 55800–60600; $\gamma$-ray light curves from the Fermi-LAT Lightcurve Repository (3-day bins, MJD 55000–60600) plus unbinned Pass 8 likelihood analysis in 0.1–300 GeV; Swift-XRT spectra (0.3–10 keV) fitted with an absorbed power law; and radio light curves at 1 mm (SMA), 37 GHz (Metsähovi), and 15 GHz (MOJAVE).

Optical variability

Intraday variability (IDV) was assessed conservatively: a night was counted as variable only when the $\chi^2$, $F$, and one-way ANOVA tests all rejected constancy at the 99% confidence level. IDV was confirmed on four nights (MJD 58420, 58515, 58821, and 60275), with variability amplitudes between 3.88% ($R$0 band, MJD 58515) and 29.90% ($R$1 band, MJD 60275); most other detections were $R$2. The 29.90% $R$3-band amplitude is by far the largest reported for this source and stands out against the otherwise modest intraday fluctuations, although the $R$4 and $R$5 data on that night were of insufficient quality to confirm IDV in those bands.

Color behavior was examined through color–magnitude diagrams using BCES regression and Spearman statistics. A significant bluer-when-brighter (BWB) trend was found on MJD 58821, while MJD 58420 showed a weak redder-when-brighter (RWB) tendency with $R$6 but a large $R$7-value, so the authors treat it cautiously. The BWB trend is interpreted within the shock-in-jet framework as rapid electron acceleration relative to cooling; the RWB tendency may reflect the source's masquerading-BL-Lac/FSRQ nature, where a less-variable quasar-like component dilutes the color response. Notably, no color reversal or spectral hysteresis loops were detected, which the authors take to indicate a stable ratio of acceleration to cooling timescales during the observed events. On long timescales the source dimmed overall, from $R$8 mag on MJD 58418 to $R$9 mag on MJD 60274, a total amplitude of approximately 145.6%, with a clear long-term BWB trend that the authors attribute to intrinsic energetic processes such as synchrotron-peak shifts rather than pure geometric effects.

Interband time lags

Using the interpolated cross-correlation function (ICCF) with 10,000 Monte Carlo realizations and alias removal, no significant optical interband lags were found on intraday timescales, implying co-spatial $I$0-, $I$1-, and $I$2-band emission regions within the observational precision. On multi-year timescales the picture is different: high-energy emission (optical to $I$3-ray) is co-spatial and contemporaneous, but the radio emission lags by approximately 800–900 days. Specifically, the $I$4-ray versus 1 mm lag is $I$5 days, versus 37 GHz $I$6 days, and versus 15 GHz $I$7 days, forming a coherent frequency-stratified sequence (1 mm leads 15 GHz by $I$8 days; 37 GHz leads 15 GHz by $I$9 days). The authors contrast this result with the ~100-day optical-leads-radio lag reported in earlier work based on OVRO data that did not cover the post-2017 outburst, and they note that adiabatic expansion alone cannot explain the 2–3 year delay, favoring instead a re-dissipation event at parsec scales. The practical implication is a two-zone jet geometry: a compact upstream region producing high-energy flares and neutrinos, and a larger downstream region responsible for the delayed radio outbursts.

SED modeling

Three epochs were modeled with the time-dependent, one-zone lepto-hadronic AM$B$0 code: Epoch 1 (MJD 58027–58036), the $B$1-ray flare contemporaneous with IceCube-170922A; Epoch 2 (MJD 58830–58850), a post-flare phase coinciding with the radio outburst peak; and Epoch 3 (MJD 59830–59864), a quiescent state. The high-energy spectral analysis shows the $B$2-ray flux dropping from $B$3 (Epoch 1) to $B$4 erg cm$B$5 s$B$6 (Epoch 2), while X-ray spectra exhibit a "softer-when-brighter" behavior, with the photon index increasing from 1.77 (Epoch 2) to 2.62 (Epoch 1) as the flux roughly doubled—interpreted as injection of fresh, lower-energy electrons.

The best-fit parameters reveal a pronounced evolution in the emission region. Epoch 1 requires a compact region ($B$7 cm) with a relatively strong magnetic field ($B$8 G), high bulk Lorentz factor ($B$9), and a large injected proton luminosity ($V$0 erg s$V$1), under which the X-ray emission is hadronically dominated and a significant neutrino flux is predicted—consistent with the neutrino detection. Epochs 2 and 3 are instead lepton-dominated, with much larger radii ($V$2 and $V$3 cm), weaker fields (0.0077 and 0.032 G), and proton luminosities lower by two to nearly three orders of magnitude ($V$4 and $V$5 erg s$V$6), yielding lower predicted neutrino fluxes. A notable feature is that the Epoch 2 and 3 models reproduce the radio data reasonably well, whereas the Epoch 1 compact-zone model underestimates the radio flux—consistent with radio emission arising from extended, optically thin downstream regions.

Limitations and open questions

Several caveats bear directly on these results. The IDV detection on MJD 60275 rests on the $V$7 band alone because of poor $V$8- and $V$9-band data quality. The 1 mm light curve has a gap over MJD 58200–59300, weakening the constraints on lags involving that band, and the $R$0-ray–X-ray and $R$1-ray–optical lags were not consistently recovered across data segments or search ranges. The SED modeling assumes a single spherical blob with fixed parameters ($R$2, $R$3, $R$4, escape rate), and radio data were excluded from the fits on the standard assumption that they arise from a separate zone—an assumption that the two-zone interpretation of the time lags supports but does not prove. The paper also leaves open whether the 800–900 day radio delay is produced by re-dissipation at parsec scales or by another mechanism, and whether a single model can simultaneously account for both the 2017 flare and the 2014–2015 archival neutrino excess.

Conclusion

This work combines five years of dedicated optical monitoring with a roughly 15-year multi-broadband dataset and multi-epoch lepto-hadronic modeling of TXS 0506+056. The principal results are the detection of IDV on four nights with amplitudes up to 29.90%, a dominant BWB color trend on long timescales, no optical interband lags, a robust ~800–900 day radio lag behind the high-energy emission, and an inferred evolution from a compact, strongly magnetized, proton-loaded emission region during the neutrino-associated flare to larger, lepton-dominated, downstream regions in the post-flare and quiescent states. Together these findings provide a coherent, if phenomenological, framework linking the source's electromagnetic variability to its multi-messenger behavior.