Blazar TXS Research Articles

Investigating the Blazar TXS 0506+056 through Sharp Multiwavelength Eyes During 2017–2019

Abstract The blazar TXS 0506+056 got into the spotlight of the astrophysical community in 2017 September, when a high-energy neutrino detected by IceCube (IceCube-170922A) was associated at the 3σ level with a γ-ray flare from this source. This multi-messenger photon-neutrino association remains, as per today, the most significant association ever observed. TXS 0506+056 was a poorly studied object before the IceCube-170922A event. To better characterize its broadband emission, we organized a multiwavelength campaign lasting 16 months (2017 November to 2019 February), covering the radio band (Metsähovi, OVRO), the optical/UV (ASAS-SN, KVA, REM, Swift/UVOT), the X-rays (Swift/XRT, NuSTAR), the high-energy γ rays (Fermi/LAT), and the very high-energy (VHE) γ rays (MAGIC). In γ rays, the behavior of the source was significantly different from the behavior in 2017: MAGIC observations show the presence of flaring activity during 2018 December, while the source only shows an excess at the 4σ level during the rest of the campaign (74 hr of accumulated exposure); Fermi/LAT observations show several short (on a timescale of days to a week) flares, different from the long-term brightening of 2017. No significant flares are detected at lower energies. The radio light curve shows an increasing flux trend that is not seen in other wavelengths. We model the multiwavelength spectral energy distributions in a lepto-hadronic scenario, in which the hadronic emission emerges as Bethe-Heitler and pion-decay cascade in the X-rays and VHE γ rays. According to the model presented here, the 2018 December γ-ray flare was connected to a neutrino emission that was too brief and not bright enough to be detected by current neutrino instruments.

The spectra of IceCube neutrino (SIN) candidate sources – II. Source characterization

ABSTRACT Eight years after the first detection of high-energy astrophysical neutrinos by IceCube, we are still almost clueless as regards to their origin, although the case for blazars being neutrino sources is getting stronger. After the first significant association at the $3\!-\!3.5\, \sigma$ level in time and space with IceCube neutrinos, i.e. the blazar TXS 0506+056 at z = 0.3365, some of us have in fact selected a unique sample of 47 blazars, out of which ∼16 could be associated with individual neutrino track events detected by IceCube. Building upon our recent spectroscopy work on these objects, here we characterize them to determine their real nature and check if they are different from the rest of the blazar population. For the first time we also present a systematic study of the frequency of masquerading BL Lacs, i.e. flat-spectrum radio quasars with their broad lines swamped by non-thermal jet emission, in a γ-ray- and IceCube-selected sample, finding a fraction >24 per cent and possibly as high as 80 per cent. In terms of their broad-band properties, our sources appear to be indistinguishable from the rest of the blazar population. We also discuss two theoretical scenarios for neutrino emission, one in which neutrinos are produced in interactions of protons with jet photons and one in which the target photons are from the broad-line region. Both scenarios can equally account for the neutrino–blazar correlation observed by some of us. Future observations with neutrino telescopes and X-ray satellites will test them out.

Hadronuclear interpretation of the possible neutrino emission from PKS B1424-418, GB6 J1040+0617 and PKS 1502+106

In addition to neutrino event IceCube-170922A which is observed to be associated with a γ-ray flare from blazar TXS 0506+056, there are also several neutrino events that may be associated with blazars. Among them, PKS B1424-418, GB6 J1040+0617 and PKS 1502+106 are low synchrotron peaked sources, which are usually believed to have the broad line region in the vicinity of the central black hole. They are considered as counterparts of IceCube event 35, IceCube-141209A and IceCube-190730A, respectively. By considering the proton-proton (pp) interactions between the dense gas clouds in the broad line region and the relativistic protons in the jet, we show that the pp model that is applied in this work can not only reproduce the multi-waveband spectral energy distribution but also suggest a considerable annual neutrino detection rate. We also discuss the emission from the photopion production and Bethe-Heitler pair production with a sub-Eddington jet power that is suggested in our model and find that it has little effect on the spectrum of total emission for all of three sources.

Multimessenger Astronomy of Transient Point Sources at the Pierre Auger Observatory

One of the key challenges in astroparticle physics is the identification of the sources of cosmic rays at the highest energies (above 1018 eV). In this context, the search for neutral messenger particles in the ultra-high energy (UHE) regime is of high interest. The sources of the gravitational waves (GWs) that can be observed with the current generation of GW detectors provide extreme astrophysical environments that are most likely to be unique in the universe. Another extraordinary source candidate is the anomalous blazar TXS 0506+056 which has been found to be coincident with two periods of enhanced high-energy neutrino flux in 2014/15 and 2017 as reported by the IceCube Collaboration. Due to their distance and transient nature, the capabilities of these sources to produce UHE radiation can only be studied through neutral messengers like photons and neutrinos. The Pierre Auger Observatory near Malargue, Argentina, is the largest air shower experiment for the detection of ultra-high energy cosmic rays. With its surface detector, consisting of a grid of 1660 water-Cherenkov detectors covering an area of 3000 km2, it has a unique exposure to UHE photons and neutrinos and has published first constraints on these particles from GW sources and TXS 0506+056.

The relativistic parsec-scale jets of the blazars TXS 0506+056 and PKS 0502+049 and their possible association with gamma-ray flares and neutrino production

ABSTRACT The physical nature of the mechanism responsible for the emission of neutrinos in active galactic nuclei (AGNs) has been matter of debate in the literature, with relativistic jets of radio-loud AGNs as possible candidates to be the sources of high-energy neutrinos. The most prominent candidate so far is the blazar TXS 0506+056, which is found to be associated with the neutrino event IceCube-170922A. Furthermore, the IceCube reported an excess of neutrinos towards TXS 0506+056 between September 2014 and March 2015, even though this association needs additional investigation, considering the presence of a nearby gamma-ray source, the quasar PKS 0502+049. Motivated by this, we studied the parsec-scale structures of TXS 0506+056 and PKS 0502+049 through radio interferometry at 8 and 15 GHz. We identified twelve jet components in TXS 0506+056 and seven components in PKS 0502+049. The most reliable jet components show superluminal speeds ranging from 9.5c to 66c in the case of TXS 0506+056, and from 14.3c to 59c for PKS 0502+049, which were used to estimate a lower (upper) limit for the Lorentz factor (jet viewing angle) for both sources. A novel approach using simultaneously the brightness temperature of the core region and the apparent speeds of the jet components allowed us to infer basic jet parameters for TXS 0506+056 at distinct epochs. We also found that the emergence of new jet components coincides with the occurrence of gamma-ray flares. Interestingly, two of these coincidences in the case of PKS 0502+049 and one for TXS 0506+056 seems to be correlated with neutrino events detected by the IceCube Observatory.

Time-dependent lepto-hadronic modelling of the emission from blazar jets with SOPRANO: the case of TXS 0506 + 056, 3HSP J095507.9 + 355101, and 3C 279

ABSTRACT The observation of a very high-energy neutrino by IceCube (IceCube-170922A) and its association with the flaring blazar TXS 0506 + 056 provided the first multimessenger observations of blazar jets, demonstrating the important role of protons in their dynamics and emission. In this paper, we present SOPRANO (https://www.amsdc.am/soprano), a new conservative implicit kinetic code that follows the time evolution of the isotropic distribution functions of protons, neutrons, and the secondaries produced in photo-pion and photo-pair interactions, alongside with the evolution of photon and electron/positron distribution functions. SOPRANO is designed to study leptonic and hadronic processes in relativistic sources such as blazars and gamma-ray bursts. Here, we use SOPRANO to model the broadband spectrum of TXS 0506 + 056 and 3HSP J095507.9 + 355101, which are associated with neutrino events, and of the extreme flaring blazar 3C 279. The SEDs are interpreted within the guise of both a hadronic and a hybrid model. We discuss the implications of our assumptions in terms of jet power and neutrino flux.

A Test of the Hadronic Origin of γ-Rays from Blazars with Follow-up up to a Month Later of IceCube Alerts with Imaging Air Cherenkov Telescopes

Abstract The sources of IceCube neutrinos are as yet unknown. The multimessenger observation of their emission in γ-rays can be a guide to their identification, as exemplified by the case of TXS 0506+056. We suggest a new method of searching for γ-rays with Imaging Air Cherenkov Telescopes from sources coinciding with possible astrophysical neutrinos. We propose that searches of γ-rays be extended from the current practice of only a few days to up to one month from a neutrino alert. We test this strategy on simulated sources modeled after the blazar TXS 0506+056, emitting neutrinos and γ-rays via photohadronic interactions: the γ-rays are subsequently reprocessed in the VHE range. Using MAGIC as a benchmark example, we show that current Cherenkov Telescopes should be able to detect γ-ray counterparts to neutrino alerts with a rate of approximately one per year. It has been proposed that the high-energy diffuse neutrino flux can be explained by ∼5% of all blazars flaring in neutrinos once every 10 yr, with a neutrino luminosity similar to that of TXS 0506+056 during the 2014–2015 neutrino flare. The implementation of our strategy could lead, over a timescale of one or a few years, either to the detection of this subclass of blazars contributing to the diffuse neutrino flux or to a constraint on this model.

First detection of VHE gamma-ray emission from TXS 1515–273, study of its X-ray variability and spectral energy distribution

ABSTRACT We report here on the first multiwavelength (MWL) campaign on the blazar TXS 1515–273, undertaken in 2019 and extending from radio to very-high-energy gamma-rays (VHE). Up until now, this blazar had not been the subject of any detailed MWL observations. It has a rather hard photon index at GeV energies and was considered a candidate extreme high-synchrotron-peaked source. MAGIC observations resulted in the first-time detection of the source in VHE with a statistical significance of 7.6σ. The average integral VHE flux of the source is 6 ± 1 per cent of the Crab nebula flux above 400 GeV. X-ray coverage was provided by Swift-XRT, XMM–Newton, and NuSTAR. The long continuous X-ray observations were separated by ∼9 h, both showing clear hour scale flares. In the XMM–Newton data, both the rise and decay time-scales are longer in the soft X-ray than in the hard X-ray band, indicating the presence of a particle cooling regime. The X-ray variability time-scales were used to constrain the size of the emission region and the strength of the magnetic field. The data allowed us to determine the synchrotron peak frequency and classify the source as a flaring high, but not extreme synchrotron-peaked object. Considering the constraints and variability patterns from the X-ray data, we model the broad-band spectral energy distribution. We applied a simple one-zone model, which could not reproduce the radio emission and the shape of the optical emission, and a two-component leptonic model with two interacting components, enabling us to reproduce the emission from radio to VHE band.

Cosmic Neutrinos from Temporarily Gamma-suppressed Blazars

Abstract Despite the uncovered association of a high-energy neutrino with the apparent flaring state of blazar TXS 0506+056 in 2017, the mechanisms leading to astrophysical particle acceleration and neutrino production are still uncertain. Recent studies found that when transparent to γ-rays, γ-flaring blazars do not have the opacity for protons to produce neutrinos. Here we present observational evidence for an alternative explanation, in which γ-ray emission is suppressed during efficient neutrino production. A large proton and target photon density helps produce neutrinos while temporarily suppressing the observable γ-emission due to a large γ γ opacity. We show that the Fermi-LAT γ-flux of blazar PKS 1502+106 was at a local minimum when IceCube recorded the coincident high-energy neutrino IC-190730A. Using data from the OVRO 40 m Telescope, we find that radio emission from PKS 1502+106 at the time period of the coincident neutrino IC-190730A was in a high state, in contrast to earlier time periods when radio and γ fluxes are correlated for both low and high states. This points to an active outflow that is γ-suppressed at the time of neutrino production. We find similar local γ-suppression in other blazars, including in MAGIC’s TeV flux of TXS 0506+056 and Fermi-LAT’s flux of blazar PKS B1424-418 at the time of coincident IceCube neutrino detections. Using temporary γ-suppression, neutrino–blazar coincidence searches could be substantially more sensitive than previously assumed, enabling the identification of the origin of IceCube’s diffuse neutrino flux possibly with already existing data.

High-Energy Neutrino Astronomy—Baikal-GVD Neutrino Telescope in Lake Baikal

High-energy neutrino astronomy is a fascinating new field of research, rapidly developing over recent years. It opens a new observation window on the most violent processes in the universe, fitting very well to the concept of multi-messenger astronomy. This may be exemplified by the recent discovery of the high-energy neutrino emissions from the γ-ray loud blazar TXS 0506+056. Constraining astrophysical neutrino fluxes can also help to understand the long-standing mystery of the origin of the ultra-high energy cosmic rays. Astronomical studies of high-energy neutrinos are carried out by large-scale next-generation neutrino telescopes located in different regions of the world, forming a global network of complementary detectors. The Baikal-GVD, being currently the largest neutrino telescope in the Northern Hemisphere and still growing up, is an important constituent of this network. This paper briefly reviews working principles, analysis methods, and some selected results of the Baikal-GVD neutrino telescope.

Hadronic X-Ray Flares from Blazars

Abstract The detection of a high-energy neutrino from the flaring blazar TXS 0506+056 and the subsequent discovery of a neutrino excess from the same direction have strengthened the hypothesis that blazars are cosmic neutrino sources. The lack, however, of γ-ray flaring activity during the latter period challenges the standard scenario of correlated γ-ray and high-energy neutrino emission in blazars. We propose instead that TeV–PeV neutrinos are produced in coincidence with X-ray flares that are powered by proton synchrotron radiation. In this case, neutrinos are produced by photomeson interactions of protons with their own synchrotron radiation, while MeV to GeV γ-rays are the result of synchrotron-dominated electromagnetic cascades developed in the source. Using a time-dependent approach, we find that this “pure hadronic flaring” hypothesis has several interesting consequences. The X-ray flux is a good proxy for the all-flavor neutrino flux, while certain neutrino-rich X-ray flares may be dark in GeV–TeV γ-rays. Lastly, hadronic X-ray flares are accompanied by an equally bright MeV component that is detectable by proposed missions like e-ASTROGAM and AMEGO. We applied this scenario to the extreme blazar 3HSP J095507.9+355101 which has been associated with IceCube-200107A while undergoing an X-ray flare. We showed that the number of muon and anti-muon neutrinos above 100 TeV during hadronic flares can be up to ∼3–10 times higher than the expected number in standard leptohadronic models. Still, frequent hadronic flaring activity is necessary for explaining the detected neutrino event IceCube-200107A.

A Two-zone Blazar Radiation Model for “Orphan” Neutrino Flares

Abstract In this work, we investigate the 2014–2015 neutrino flare associated with the blazar TXS 0506+056 and a recently discovered muon neutrino event IceCube-200107A in spatial coincidence with the blazar 4FGL J0955.1+3551, under the framework of a two-zone radiation model of blazars where an inner/outer blob close to/far from the supermassive black hole is invoked. An interesting feature that the two sources have in common is that no evidence of GeV gamma-ray activity is found during the neutrino detection period, probably implying a large opacity for GeV gamma rays in the neutrino production region. In our model, continuous particle acceleration/injection takes place in the inner blob at the jet base, where the hot X-ray corona of the supermassive black hole provides target photon fields for efficient neutrino production and strong GeV gamma-ray absorption. We show that this model can self-consistently interpret the neutrino emission from both blazars in a large parameter space. In the meantime, the dissipation processes in outer blob are responsible for the simultaneous multiwavelength emission of both sources. In agreement with previous studies of TXS 0506+056, an intense MeV emission from the induced electromagnetic cascade in the inner blob is robustly expected to accompany the neutrino flare in our model and could be used to test the model using the next-generation MeV gamma-ray detector in the future.

Recurrent Neutrino Emission from Supermassive Black Hole Mergers

Abstract The recent detection of possible neutrino emission from the blazar TXS 0506+056 was the first high-energy neutrino associated with an astrophysical source, making this special type of active galaxies a promising neutrino emitter. The fact that two distinct episodes of neutrino emission were detected with a separation of around 3 yr raises the possibility that emission could be periodic. Periodic emission is expected from supermassive binary black hole systems due to jet precession close to the binary’s merger. Here, we show that, if TXS 0506+056 is a binary source, then the next neutrino flare could already have occurred, possibly still hidden in IceCube’s not-yet-analyzed data. We derive the binary properties that would lead to the detection of gravitational waves from this system by the Laser Interferometer Space Antenna (LISA) over the next decade. Our results for the first time quantify the timescale of these correlations for the example of TXS 0506+056, connecting the possible neutrino and gravitational-wave signatures of such sources.

The intrinsic gamma-ray spectrum of TXS 0506+056: intergalactic propagation effects

ABSTRACT The recent observation of high-energy neutrinos from the 2017 flare of the blazar TXS 0506+056, together with counterparts across the whole electromagnetic spectrum, opens up new possibilities for investigating the properties of this class of objects as well as the traversed medium. Propagation effects such as the attenuation of the very-high-energy gamma-ray component by the extragalactic background light are well known, and usually taken into account when fitting spectral energy distributions of objects. Other effects such as those of intergalactic magnetic fields are, however, often neglected. In this work, we present a comprehensive study of the influence of these fields and the extragalactic background light on the determination of the intrinsic gamma-ray spectrum of this blazar.

On the Detection Potential of Blazar Flares for Current Neutrino Telescopes

Abstract Blazar jets are extreme environments, in which relativistic proton interactions with an ultraviolet photon field could give rise to photopion production. High-confidence associations of individual high-energy neutrinos with blazar flares could be achieved via spatially and temporally coincident detections. In 2017, the track-like, extremely high-energy neutrino event IC 170922A was found to coincide with increased γ-ray emission from the blazar TXS 0506+056, leading to the identification of the most promising neutrino point-source candidate so far. We calculate the expected number of neutrino events that can be detected with IceCube, based on a broadband parameterization of bright short-term blazar flares that were observed in the first 6.5 yr of Fermi/Large Area Telescope observations. We find that the integrated keV-to-GeV fluence of most individual blazar flares is far too small to yield a substantial Poisson probability for the detection of one or more neutrinos with IceCube. We show that the sample of potentially detectable high-energy neutrinos from individual blazar flares is rather small. We further show that the blazars 3C 279 and PKS 1510−089 dominate the all-sky neutrino prediction from bright and short-term blazar flares. In the end, we discuss strategies to search for more significant associations in future data unblindings of IceCube and KM3NeT.

A Search for Ultra-high-energy Neutrinos from TXS 0506+056 Using the Pierre Auger Observatory

Abstract Results of a search for ultra-high-energy neutrinos with the Pierre Auger Observatory from the direction of the blazar TXS 0506+056 are presented. They were obtained as part of the follow-up that stemmed from the detection of high-energy neutrinos and gamma rays with IceCube, Fermi-LAT, MAGIC, and other detectors of electromagnetic radiation in several bands. The Pierre Auger Observatory is sensitive to neutrinos in the energy range from 100 PeV to 100 EeV and in the zenith-angle range from θ = 60° to θ = 95°, where the zenith angle is measured from the vertical direction. No neutrinos from the direction of TXS 0506+056 have been found. The results were analyzed in three periods: one of 6 months around the detection of IceCube-170922 A, coinciding with a flare period of TXS 0506+056, a second one of 110 days during which the IceCube collaboration found an excess of 13 neutrinos from a direction compatible with TXS 0506+056, and a third one from 2004 January 1 up to 2018 August 31, over which the Pierre Auger Observatory has been taking data. The sensitivity of the Observatory is addressed for different spectral indices by considering the fluxes that would induce a single expected event during the observation period. For indices compatible with those measured by the IceCube collaboration the expected number of neutrinos at the Observatory is well below one. Spectral indices as hard as 1.5 would have to apply in this energy range to expect a single event to have been detected.

Multifrequency Observations of the Candidate Neutrino-emitting Blazar BZB J0955+3551

Abstract The recent spatial and temporal coincidence of the blazar TXS 0506+056 with the IceCube-detected neutrino event IC-170922A has opened up a realm of multimessenger astronomy with blazar jets as a plausible site of cosmic-ray acceleration. After TXS 0506+056, a second blazar, BZB J0955+3551, was recently found to be spatially coincident with the IceCube-detected neutrino event IC-200107A and undergoing its brightest X-ray flare measured so far. Here we present the results of our multifrequency campaign to study this peculiar event that includes observations with the NuSTAR, Swift, Neutron star Interior Composition Explorer (NICER), and 10.4 m Gran Telescopio Canarias (GTC). The optical spectroscopic observation from GTC secured its redshift as and the central black hole mass as . Both NuSTAR and NICER data reveal a rapid flux variability, albeit at low significance ( ). We explore the origin of the target photon field needed for the photopion production using analytical calculations and considering the observed optical-to-X-ray flux level. We conclude that seed photons may originate from outside the jet, similar to that reported for TXS 0506+056, although a scenario invoking a comoving target photon field (e.g., electron synchrotron) cannot be ruled out. The electromagnetic output from the neutrino-producing photohadronic processes are likely to make only a subdominant contribution to the observed spectral energy distribution, suggesting that the X-ray flaring event may not be directly connected with IC-200107A.

Multimessenger Constraints on Intergalactic Magnetic Fields from the Flare of TXS 0506+056

Abstract The origin of magnetic fields in the universe is an open problem. Seed magnetic fields possibly produced in early times may have survived up to the present day close to their original form, providing an untapped window to the primeval universe. The recent observations of high-energy neutrinos from the blazar TXS 0506+056 in association with an electromagnetic counterpart in a broad range of wavelengths can be used to probe intergalactic magnetic fields via the time delay between the neutrinos and gamma-rays as well as the time dependence of the gamma-ray fluxes. Using extensive three-dimensional Monte Carlo simulations, we constrain both the magnetic-field strength and, for the first time, its coherence length, considering six orders of magnitude for each.

Limiting superluminal neutrino velocity and Lorentz invariance violation by neutrino emission from the blazar TXS 0506+056

The detection of high-energy neutrino coincident with the blazar TXS 0506+056 provides a unique opportunity to test Lorentz invariance violation (LIV) in the neutrino sector. Thanks to the precisely measured redshift, i.e., $z=0.3365$, the comoving distance of the neutrino source is determined. In this work, we obtain and discuss the constraints on the superluminal neutrino velocity $\delta_\nu$ and the LIV by considering the energy loss of superluminal neutrino during propagation. Given superluminal electron velocity ($\delta_e \ge 0$), a very stringent constraint on superluminal neutrino velocity can be reached, i.e., $\delta_\nu \lesssim 1.3\times 10^{-18}$, corresponding to the quantum gravity (QG) scale $M_{\rm QG,1} \gtrsim 5.7 \times 10^{3} M_{\rm Pl}$ and $M_{\rm QG,2} \gtrsim 9.3 \times 10^{-6} M_{\rm Pl}$ for linear (quadratic) LIV, which are $\sim 12$ orders of magnitude tighter for linear LIV and $\sim 9$ orders tighter for quadratic LIV compared to the time-of-flight constraint from MeV neutrinos of SN 1987A. While given the subluminal electron velocity, a weaker constraint on the superluminal neutrino velocity is obtained, i.e., $\delta_\nu \lesssim 8 \times 10^{-17}$, which is consistent with the conclusions of previous works. We also study the neutrino detection probability due to the distortion of neutrino spectral shape during propagation, which gives slightly weaker constraints than above by a factor of $\sim2$.

Electron–positron pair plasma in TXS 0506+056 and the ‘neutrino flare’ in 2014–2015

ABSTRACT The detection of a prolonged flaring activity from blazar TXS 0506+056 in temporal and spatial coincidence with the energetic neutrino IceCube-170922A provided evidence about the photohadronic interactions in this source. However, analysis of the archival neutrino and multiwavelength data from the direction of this blazar between 2014 September and 2015 March revealed a ‘neutrino flare’ without observing quasi-simultaneous activity in the gamma-ray bands, posing challenges to established models. Electron–positron (e±) pairs generated from the accretion discs have been amply proposed as a mechanism of bulk acceleration of sub-relativistic and relativistic jets. These pairs annihilate inside the source producing a line around the electron mass, which is blueshifted in the observed frame (on Earth) and redshifted in the frame of the dissipation region of the jet. The redshifted photons in the dissipation region interact with accelerated protons, producing high-energy neutrinos that contribute significantly to the diffuse neutrino flux in the ∼10–20 TeV energy range in connection with gamma-rays from the photopion process, which can be detected by future MeV orbiting satellites. Based on this phenomenological model, we can explain the ‘neutrino flare’ reported in 2014–1015.