-
Association of the IceCube neutrinos with blazars in the CGRaBS sample
Authors:
Pouya M. Kouch,
Elina Lindfors,
Talvikki Hovatta,
Ioannis Liodakis,
Karri I. I. Koljonen,
Kari Nilsson,
Sebastian Kiehlmann,
Walter Max-Moerbeck,
Anthony C. S. Readhead,
Rodrigo A. Reeves,
Timothy J. Pearson,
Jenni Jormanainen,
Vandad Fallah Ramazani,
Matthew J. Graham
Abstract:
The origin of high-energy (HE) astrophysical neutrinos has remained an elusive hot topic in the field of HE astrophysics for the past decade. Apart from a handful of individual associations, the vast majority of HE neutrinos arise from unknown sources. While there are theoretically-motivated candidate populations, such as blazars -- a subclass of AGN with jets pointed towards our line-of-sight --…
▽ More
The origin of high-energy (HE) astrophysical neutrinos has remained an elusive hot topic in the field of HE astrophysics for the past decade. Apart from a handful of individual associations, the vast majority of HE neutrinos arise from unknown sources. While there are theoretically-motivated candidate populations, such as blazars -- a subclass of AGN with jets pointed towards our line-of-sight -- they have not yet been convincingly linked to HE neutrino production. Here, we perform a spatio-temporal association analysis between a sample of blazars (from CGRaBS catalog) in the radio and optical bands and the most up-to-date IceCube HE neutrino catalog. We find that if the IceCube error regions are enlarged by 1$^\circ$ in quadrature, to account for unknown systematic errors at maximal level, a spatio-temporal correlation between the multiwavelength light curves of the CGRaBS blazars and the IceCube HE neutrinos is hinted at least at a 2.17$σ$ significance level. On the other hand, when the IceCube error regions are taken as their published values, we do not find any significant correlations. A discrepancy in the blazar-neutrino correlation strengths, when using such minimal and enlarged error region scenarios, was also obtained in a recent study by the IceCube collaboration. In our study, this difference arises because several flaring blazars -- coinciding with a neutrino arrival time -- happen to narrowly miss the published 90\%-likelihood error region of the nearest neutrino event. For all of the associations driving our most significant correlations, the flaring blazar is much less than 1$^\circ$ away from the published error regions. Therefore, our results indicate that the question of the blazar-neutrino connection is highly sensitive to the reconstruction of the neutrino error regions, whose reliability is expected to improve with the next generation of neutrino observatories.
△ Less
Submitted 9 July, 2024;
originally announced July 2024.
-
The variability patterns of the TeV blazar PG 1553+113 from a decade of MAGIC and multi-band observations
Authors:
MAGIC Collaboration,
H. Abe,
S. Abe,
J. Abhir,
V. A. Acciari,
I. Agudo,
T. Aniello,
S. Ansoldi,
L. A. Antonelli,
A. Arbet Engels,
C. Arcaro,
M. Artero,
K. Asano,
D. Baack,
A. Babić,
A. Baquero,
U. Barres de Almeida,
I. Batković,
J. Baxter,
J. Becerra González,
E. Bernardini,
J. Bernete,
A. Berti,
J. Besenrieder,
C. Bigongiari
, et al. (242 additional authors not shown)
Abstract:
PG 1553+113 is one of the few blazars with a convincing quasi-periodic emission in the gamma-ray band. The source is also a very high-energy (VHE; >100 GeV) gamma-ray emitter. To better understand its properties and identify the underlying physical processes driving its variability, the MAGIC Collaboration initiated a multiyear, multiwavelength monitoring campaign in 2015 involving the OVRO 40-m a…
▽ More
PG 1553+113 is one of the few blazars with a convincing quasi-periodic emission in the gamma-ray band. The source is also a very high-energy (VHE; >100 GeV) gamma-ray emitter. To better understand its properties and identify the underlying physical processes driving its variability, the MAGIC Collaboration initiated a multiyear, multiwavelength monitoring campaign in 2015 involving the OVRO 40-m and Medicina radio telescopes, REM, KVA, and the MAGIC telescopes, Swift and Fermi satellites, and the WEBT network. The analysis presented in this paper uses data until 2017 and focuses on the characterization of the variability. The gamma-ray data show a (hint of a) periodic signal compatible with literature, but the X-ray and VHE gamma-ray data do not show statistical evidence for a periodic signal. In other bands, the data are compatible with the gamma-ray period, but with a relatively high p-value. The complex connection between the low and high-energy emission and the non-monochromatic modulation and changes in flux suggests that a simple one-zone model is unable to explain all the variability. Instead, a model including a periodic component along with multiple emission zones is required.
△ Less
Submitted 4 March, 2024;
originally announced March 2024.
-
Unprecedented extreme high-frequency radio variability in early-stage active galactic nuclei
Authors:
E. Järvelä,
T. Savolainen,
M. Berton,
A. Lähteenmäki,
S. Kiehlmann,
T. Hovatta,
I. Varglund,
A. C. S. Readhead,
M. Tornikoski,
W. Max-Moerback,
R. A. Reeves,
S. Suutarinen
Abstract:
We report on the discovery of one of the most extreme cases of high-frequency radio variability ever measured in active galactic nuclei (AGN), observed on timescales of days and exhibiting variability amplitudes of three to four orders of magnitude. These sources, all radio-weak narrow-line Seyfert 1 (NLS1) galaxies, were discovered some years ago at Aalto University Metsähovi Radio Observatory (M…
▽ More
We report on the discovery of one of the most extreme cases of high-frequency radio variability ever measured in active galactic nuclei (AGN), observed on timescales of days and exhibiting variability amplitudes of three to four orders of magnitude. These sources, all radio-weak narrow-line Seyfert 1 (NLS1) galaxies, were discovered some years ago at Aalto University Metsähovi Radio Observatory (MRO) based on recurring flaring at 37 GHz, strongly indicating the presence of relativistic jets. In subsequent observations with the Karl G. Jansky Very Large Array (JVLA) at 1.6, 5.2, and 9.0~GHz no signs of jets were seen. To determine the cause of their extraordinary behaviour, we observed them with the JVLA at 10, 15, 22, 33, and 45 GHz, and with the Very Long Baseline Array (VLBA) at 15 GHz. These observations were complemented with single-dish monitoring at 37 GHz at MRO, and at 15 GHz at Owens Valley Radio Observatory (OVRO). Intriguingly, all but one source either have a steep radio spectrum up to 45 GHz, or were not detected at all. Based on the 37 GHz data the timescales of the radio flares are a few days, and the derived variability brightness temperatures and variability Doppler factors comparable to those seen in blazars. We discuss alternative explanations for their extreme behaviour, but so far no definite conclusions can be made. These sources exhibit radio variability at a level rarely, if ever, seen in AGN. They might represent a new type of jetted AGN, or a new variability phenomenon, and thus deserve our continued attention.
△ Less
Submitted 4 December, 2023;
originally announced December 2023.
-
CLASS Observations of Atmospheric Cloud Polarization at Millimeter Wavelengths
Authors:
Yunyang Li,
John W. Appel,
Charles L. Bennett,
Ricardo Bustos,
David T. Chuss,
Joseph Cleary,
Jullianna Denes Couto,
Sumit Dahal,
Rahul Datta,
Rolando Dünner,
Joseph R. Eimer,
Thomas Essinger-Hileman,
Kathleen Harrington,
Jeffrey Iuliano,
Tobias A. Marriage,
Matthew A. Petroff,
Rodrigo A. Reeves,
Karwan Rostem,
Rui Shi,
Deniz A. N. Valle,
Duncan J. Watts,
Oliver F. Wolff,
Edward J. Wollack,
Zhilei Xu
Abstract:
The dynamic atmosphere imposes challenges to ground-based cosmic microwave background observation, especially for measurements on large angular scales. The hydrometeors in the atmosphere, mostly in the form of clouds, scatter the ambient thermal radiation and are known to be the main linearly polarized source in the atmosphere. This scattering-induced polarization is significantly enhanced for ice…
▽ More
The dynamic atmosphere imposes challenges to ground-based cosmic microwave background observation, especially for measurements on large angular scales. The hydrometeors in the atmosphere, mostly in the form of clouds, scatter the ambient thermal radiation and are known to be the main linearly polarized source in the atmosphere. This scattering-induced polarization is significantly enhanced for ice clouds due to the alignment of ice crystals under gravity, which are also the most common clouds seen at the millimeter-astronomy sites at high altitudes. This work presents a multifrequency study of cloud polarization observed by the Cosmology Large Angular Scale Surveyor (CLASS) experiment on Cerro Toco in the Atacama Desert of northern Chile, from 2016 to 2022, at the frequency bands centered around 40, 90, 150, and 220 GHz. Using a machine-learning-assisted cloud classifier, we made connections between the transient polarized emission found in all four frequencies with the clouds imaged by monitoring cameras at the observing site. The polarization angles of the cloud events are found to be mostly $90^\circ$ from the local meridian, which is consistent with the presence of horizontally aligned ice crystals. The 90 and 150 GHz polarization data are consistent with a power law with a spectral index of $3.90\pm0.06$, while an excess/deficit of polarization amplitude is found at 40/220 GHz compared with a Rayleigh scattering spectrum. These results are consistent with Rayleigh-scattering-dominated cloud polarization, with possible effects from supercooled water absorption and/or Mie scattering from a population of large cloud particles that contribute to the 220 GHz polarization.
△ Less
Submitted 13 September, 2023;
originally announced September 2023.
-
CLASS Angular Power Spectra and Map-Component Analysis for 40 GHz Observations through 2022
Authors:
Joseph R. Eimer,
Yunyang Li,
Michael K. Brewer,
Rui Shi,
Aamir Ali,
John W. Appel,
Charles L. Bennett,
Sarah Marie Bruno,
Ricardo Bustos,
David T. Chuss,
Joseph Cleary,
Sumit Dahal,
Rahul Datta,
Jullianna Denes Couto,
Kevin L. Denis,
Rolando Dünner,
Thomas Essinger-Hileman,
Pedro Fluxá,
Johannes Hubmayer,
Kathleen Harrington,
Jeffrey Iuliano,
John Karakla,
Tobias A. Marriage,
Carolina Núñez,
Lucas Parker
, et al. (9 additional authors not shown)
Abstract:
Measurement of the largest angular scale ($\ell < 30$) features of the cosmic microwave background (CMB) polarization is a powerful way to constrain the optical depth to reionization and search for the signature of inflation through the detection of primordial $B$-modes. We present an analysis of maps covering 73.6\% of the sky made from the $40\,\mathrm{GHz}$ channel of the Cosmology Large Angula…
▽ More
Measurement of the largest angular scale ($\ell < 30$) features of the cosmic microwave background (CMB) polarization is a powerful way to constrain the optical depth to reionization and search for the signature of inflation through the detection of primordial $B$-modes. We present an analysis of maps covering 73.6\% of the sky made from the $40\,\mathrm{GHz}$ channel of the Cosmology Large Angular Scale Surveyor (CLASS) from 2016 August to 2022 May. Taking advantage of the measurement stability enabled by front-end polarization modulation and excellent conditions from the Atacama Desert, we show this channel achieves higher sensitivity than the analogous frequencies from satellite measurements in the range $10 < \ell < 100$. Simulations show the CLASS linear (circular) polarization maps have a white noise level of $125 \,(130)\,\mathrm{μK\, arcmin}$. We measure the Galaxy-masked $EE$ and $BB$ spectra of diffuse synchrotron radiation and compare to space-based measurements at similar frequencies. In combination with external data, we expand measurements of the spatial variations of the synchrotron spectral energy density (SED) to include new sky regions and measure the diffuse SED in the harmonic domain. We place a new upper limit on a background of circular polarization in the range $5 < \ell < 125$ with the first bin showing $D_\ell < 0.023$ $\mathrm{μK^2_{CMB}}$ at 95\% confidence. These results establish a new standard for recovery of the largest-scale CMB polarization from the ground and signal exciting possibilities when the higher sensitivity and higher-frequency CLASS channels are included in the analysis.
△ Less
Submitted 14 February, 2024; v1 submitted 1 September, 2023;
originally announced September 2023.
-
Long-term multi-wavelength study of 1ES 0647+250
Authors:
MAGIC Collaboration,
V. A. Acciari,
T. Aniello,
S. Ansoldi,
L. A. Antonelli,
A. Arbet Engels,
C. Arcaro,
M. Artero,
K. Asano,
D. Baack,
A. Babić,
A. Baquero,
U. Barres de Almeida,
J. A. Barrio,
I. Batković,
J. Becerra González,
W. Bednarek,
E. Bernardini,
M. Bernardos,
A. Berti,
J. Besenrieder,
W. Bhattacharyya,
C. Bigongiari,
A. Biland,
O. Blanch
, et al. (195 additional authors not shown)
Abstract:
The BL Lac object 1ES 0647+250 is one of the few distant $γ$-ray emitting blazars detected at very high energies (VHE, $\gtrsim$100 GeV) during a non-flaring state. It was detected with the MAGIC telescopes during its low activity in the years 2009-2011, as well as during three flaring activities in the years 2014, 2019 and 2020, with the highest VHE flux in the latter epoch. An extensive multi-in…
▽ More
The BL Lac object 1ES 0647+250 is one of the few distant $γ$-ray emitting blazars detected at very high energies (VHE, $\gtrsim$100 GeV) during a non-flaring state. It was detected with the MAGIC telescopes during its low activity in the years 2009-2011, as well as during three flaring activities in the years 2014, 2019 and 2020, with the highest VHE flux in the latter epoch. An extensive multi-instrument data set was collected within several coordinated observing campaigns throughout these years. We aim to characterise the long-term multi-band flux variability of 1ES 0647+250, as well as its broadband spectral energy distribution (SED) during four distinct activity states selected in four different epochs, in order to constrain the physical parameters of the blazar emission region under certain assumptions. We evaluate the variability and correlation of the emission in the different energy bands with the fractional variability and the Z-transformed Discrete Correlation Function, as well as its spectral evolution in X-rays and $γ$ rays. Owing to the controversy in the redshift measurements of 1ES 0647+250 reported in the literature, we also estimate its distance in an indirect manner through the comparison of the GeV and TeV spectra from simultaneous observations with Fermi-LAT and MAGIC during the strongest flaring activity detected to date. Moreover, we interpret the SEDs from the four distinct activity states within the framework of one-component and two-component leptonic models, proposing specific scenarios that are able to reproduce the available multi-instrument data.
△ Less
Submitted 23 November, 2022;
originally announced November 2022.
-
Long-Timescale Stability in CMB Observations at Multiple Frequencies using Front-End Polarization Modulation
Authors:
Joseph Cleary,
Rahul Datta,
John W. Appel,
Charles L. Bennett,
David T. Chuss,
Jullianna Denes Couto,
Sumit Dahal,
Francisco Espinoza,
Thomas Essinger-Hileman,
Kathleen Harrington,
Jeffrey Iuliano,
Yunyang Li,
Tobias A. Marriage,
Carolina Nunez,
Matthew A. Petroff,
Rodrigo A. Reeves,
Rui Shi,
Duncan J. Watts,
Edward J. Wollack,
Zhilei Xu
Abstract:
The Cosmology Large Angular Scale Surveyor (CLASS) is a telescope array observing the Cosmic Microwave Background (CMB) at frequency bands centered near 40, 90, 150, and 220 GHz. CLASS measures the CMB polarization on the largest angular scales to constrain the inflationary tensor-to-scalar ratio and the optical depth due to reionization. To achieve the long time-scale stability necessary for this…
▽ More
The Cosmology Large Angular Scale Surveyor (CLASS) is a telescope array observing the Cosmic Microwave Background (CMB) at frequency bands centered near 40, 90, 150, and 220 GHz. CLASS measures the CMB polarization on the largest angular scales to constrain the inflationary tensor-to-scalar ratio and the optical depth due to reionization. To achieve the long time-scale stability necessary for this measurement from the ground, CLASS utilizes a front-end, variable-delay polarization modulator on each telescope. Here we report on the improvements in stability afforded by front-end modulation using data across all four CLASS frequencies. Across one month of modulated linear polarization data in 2021, CLASS achieved median knee frequencies of 9.1, 29.1, 20.4, and 36.4 mHz for the 40, 90, 150, and 220 GHz observing bands. The knee frequencies are approximately an order of magnitude lower than achieved via CLASS pair-differencing orthogonal detector pairs without modulation.
△ Less
Submitted 9 August, 2022;
originally announced August 2022.
-
Investigating the blazar TXS 0506+056 through sharp multi-wavelength eyes during 2017-2019
Authors:
MAGIC Collaboration,
V. A. Acciari,
T. Aniello,
S. Ansoldi,
L. A. Antonelli,
A. Arbet Engels,
M. Artero,
K. Asano,
D. Baack,
A. Babić,
A. Baquero,
U. Barres de Almeida,
J. A. Barrio,
I. Batković,
J. Becerra González,
W. Bednarek,
E. Bernardini,
M. Bernardos,
A. Berti,
J. Besenrieder,
W. Bhattacharyya,
C. Bigongiari,
A. Biland,
O. Blanch,
H. Bökenkamp
, et al. (192 additional authors not shown)
Abstract:
The blazar TXS 0506+056 got into the spotlight of the astrophysical community in September 2017, when a high-energy neutrino detected by IceCube (IceCube-170922A) was associated at the 3 $σ$ level to a $γ$-ray flare from this source. This multi-messenger photon-neutrino association remains, as per today, the most significant one ever observed. TXS 0506+056 was a poorly studied object before the Ic…
▽ More
The blazar TXS 0506+056 got into the spotlight of the astrophysical community in September 2017, when a high-energy neutrino detected by IceCube (IceCube-170922A) was associated at the 3 $σ$ level to a $γ$-ray flare from this source. This multi-messenger photon-neutrino association remains, as per today, the most significant one ever observed. TXS 0506+056 was a poorly studied object before the IceCube-170922A event. To better characterize its broad-band emission, we organized a multi-wavelength campaign lasting 16 months (November 2017 to February 2019), 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 behaviour of the source was significantly different from the 2017 one: MAGIC observations show the presence of flaring activity during December 2018, while the source only shows an excess at the 4$σ$ level during the rest of the campaign (74 hours of accumulated exposure); Fermi/LAT observations show several short (days-to-week timescale) 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, not seen in other wavelengths. We model the multi-wavelength 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 December 2018 $γ$-ray flare was connected to a neutrino emission that was too brief and not bright enough to be detected by current neutrino instruments.
△ Less
Submitted 1 May, 2022; v1 submitted 5 February, 2022;
originally announced February 2022.
-
New Tests of Millilensing in the Blazar PKS 1413+135
Authors:
A. L. Peirson,
I. Liodakis,
A. C. S. Readhead,
M. L. Lister,
E. S. Perlman,
M. F. Aller,
R. D. Blandford,
K. J. B. Grainge,
D. A. Green,
M. A. Gurwell,
M. W. Hodges,
T. Hovatta,
S. Kiehlmann,
A. Lähteenmäki,
W. Max-Moerbeck,
T. Mcaloone,
S. O'Neill,
V. Pavlidou,
T. J. Pearson,
V. Ravi,
R. A. Reeves,
P. F. Scott,
G. B. Taylor,
D. J. Titterington,
M. Tornikoski
, et al. (4 additional authors not shown)
Abstract:
Symmetric Achromatic Variability (SAV) is a rare form of radio variability in blazars that has been attributed to gravitational millilensing by a ~$10^2 - 10^5$ $M_\odot$ mass condensate. Four SAVs have been identified between 1980 and 2020 in the long-term radio monitoring data of the blazar PKS 1413+135. We show that all four can be fitted with the same, unchanging, gravitational lens model. If…
▽ More
Symmetric Achromatic Variability (SAV) is a rare form of radio variability in blazars that has been attributed to gravitational millilensing by a ~$10^2 - 10^5$ $M_\odot$ mass condensate. Four SAVs have been identified between 1980 and 2020 in the long-term radio monitoring data of the blazar PKS 1413+135. We show that all four can be fitted with the same, unchanging, gravitational lens model. If SAV is due to gravitational millilensing, PKS 1413+135 provides a unique system for studying active galactic nuclei with unprecedented microarcsecond resolution, as well as for studying the nature of the millilens itself. We discuss two possible candidates for the putative millilens: a giant molecular cloud hosted in the intervening edge-on spiral galaxy, and an undetected dwarf galaxy with a massive black hole. We find a significant dependence of SAV crossing time on frequency, which could indicate a fast shock moving in a slower underlying flow. We also find tentative evidence for a 989-day periodicity in the SAVs, which, if real, makes possible the prediction of future SAVs: the next three windows for possible SAVs begin in August 2022, May 2025, and February 2028.
△ Less
Submitted 8 January, 2022; v1 submitted 4 January, 2022;
originally announced January 2022.
-
The Unanticipated Phenomenology of the Blazar PKS~2131$-$021: A Unique Super-Massive Black Hole Binary Candidate
Authors:
S. O'Neill,
S. Kiehlmann,
A. C. S. Readhead,
M. F. Aller,
R. D. Blandford,
I. Liodakis,
M. L. Lister,
P. Mróz,
C. P. O'Dea,
T. J. Pearson,
V. Ravi,
M. Vallisneri,
K. A. Cleary,
M. J. Graham,
K. J. B. Grainge,
M. W. Hodges,
T. Hovatta,
A. Lähteenmäki,
J. W. Lamb,
T. J. W. Lazio,
W. Max-Moerbeck,
V. Pavlidou,
T. A. Prince,
R. A. Reeves,
M. Tornikoski
, et al. (2 additional authors not shown)
Abstract:
Most large galaxies host supermassive black holes in their nuclei and are subject to mergers, which can produce a supermassive black hole binary (SMBHB), and hence periodic signatures due to orbital motion. We report unique periodic radio flux density variations in the blazar PKS~2131$-$021, which strongly suggest an SMBHB with an orbital separation of $\sim 0.001-0.01$ pc. Our 45.1-year radio lig…
▽ More
Most large galaxies host supermassive black holes in their nuclei and are subject to mergers, which can produce a supermassive black hole binary (SMBHB), and hence periodic signatures due to orbital motion. We report unique periodic radio flux density variations in the blazar PKS~2131$-$021, which strongly suggest an SMBHB with an orbital separation of $\sim 0.001-0.01$ pc. Our 45.1-year radio light curve shows two epochs of strong sinusoidal variation with the same period and phase to within $<2\%$ and $\sim 10\%$, respectively, straddling a 20-year period when this variation was absent. Our simulated light curves accurately reproduce the ``red noise'' of this object, and Lomb-Scargle, weighted wavelet Z-transform, and least-squares sine wave analyses demonstrate conclusively, at the $4.6σ$ significance level, that the periodicity in this object is not due to random fluctuations in flux density. The observed period translates to $2.082\pm 0.003$ years in the rest frame at the $z=1.285$ redshift of PKS~2131$-$021. The periodic variation in PKS~2131$-$021 is remarkably sinusoidal. We present a model in which orbital motion, combined with the strong Doppler boosting of the approaching relativistic jet, produces a sine-wave modulation in the flux density which easily fits the observations. Given the rapidly-developing field of gravitational wave experiments with pulsar timing arrays, closer counterparts to PKS~2131$-$021 and searches using the techniques we have developed are strongly motivated. These results constitute a compelling demonstration that the phenomenology, not the theory, must provide the lead in this field.
△ Less
Submitted 24 January, 2022; v1 submitted 3 November, 2021;
originally announced November 2021.
-
The Relativistic Jet Orientation and Host Galaxy of the Peculiar Blazar PKS 1413+135
Authors:
A. C. S. Readhead,
V. Ravi,
I. Liodakis,
M. L. Lister,
V. Singh,
M. F. Aller,
R. D. Blandford,
I. W. A. Browne,
V. Gorjian,
K. J. B. Grainge,
M. A. Gurwell,
M. W. Hodges,
T. Hovatta,
S. Kiehlmann,
A. Lähteenmäki,
T. McAloone,
W. Max-Moerbeck,
V. Pavlidou,
T. J. Pearson,
A. L. Peirson,
E. S. Perlman,
R. A. Reeves,
B. T. Soifer,
G. B. Taylor,
M. Tornikoski
, et al. (4 additional authors not shown)
Abstract:
PKS 1413+135 is one of the most peculiar blazars known. Its strange properties led to the hypothesis almost four decades ago that it is gravitationally lensed by a mass concentration associated with an intervening galaxy. It exhibits symmetric achromatic variability, a rare form of variability that has been attributed to gravitational milli-lensing. It has been classified as a BL Lac object, and i…
▽ More
PKS 1413+135 is one of the most peculiar blazars known. Its strange properties led to the hypothesis almost four decades ago that it is gravitationally lensed by a mass concentration associated with an intervening galaxy. It exhibits symmetric achromatic variability, a rare form of variability that has been attributed to gravitational milli-lensing. It has been classified as a BL Lac object, and is one of the rare objects in this class with a visible counterjet. BL Lac objects have jet axes aligned close to the line of sight. It has also been classified as a compact symmetric object, which have jet axes not aligned close to the line of sight. Intensive efforts to understand this blazar have hitherto failed to resolve even the questions of the orientation of the relativistic jet, and the host galaxy. Answering these two questions is important as they challenge our understanding of jets in active galactic nuclei and the classification schemes we use to describe them. We show that the jet axis is aligned close to the line of sight and PKS 1413+135 is almost certainly not located in the apparent host galaxy, but is a background object in the redshift range $0.247 < z < 0.5$. The intervening spiral galaxy at $z = 0.247$ provides a natural host for the putative lens responsible for symmetric achromatic variability and is shown to be a Seyfert 2 galaxy. We also show that, as for the radio emission, a "multizone" model is needed to account for the high-energy emission.
△ Less
Submitted 7 December, 2020;
originally announced December 2020.
-
Filament rotation in the California L1482 cloud
Authors:
Rodrigo H. Álvarez-Gutiérrez,
Amelia M. Stutz,
Chi Y. Law,
Stefan Reissl,
Ralf S. Klessen,
Nathan W. C. Leigh,
Hong-Li Liu,
Rodrigo A. Reeves
Abstract:
We analyze the gas mass distribution, the gas kinematics, and the young stellar object (YSO) content of the California Molecular Cloud (CMC) L1482 filament. We derive a Gaia DR2 YSO distance of 511$^{+17}_{-16}$ pc. We derive scale-free power-laws for the mean gas line-mass (M/L) profiles; we calculate the gravitational potential and field profiles consistent with these. We present IRAM 30 m C…
▽ More
We analyze the gas mass distribution, the gas kinematics, and the young stellar object (YSO) content of the California Molecular Cloud (CMC) L1482 filament. We derive a Gaia DR2 YSO distance of 511$^{+17}_{-16}$ pc. We derive scale-free power-laws for the mean gas line-mass (M/L) profiles; we calculate the gravitational potential and field profiles consistent with these. We present IRAM 30 m C$^{18}$O (1-0) (and other tracers) position-velocity (PV) diagrams that exhibit complex velocity twisting and turning structures. We find a rotational profile in C$^{18}$O perpendicular to the southern filament ridgeline. The profile is regular, confined ($r\lesssim0.4$ pc), anti-symmetric, and to first order linear with a break at $r\sim0.25$ pc. The timescales of the inner (outer) gradients are $\sim$0.7 (6.0) Myr. We show that the centripetal force, compared to gravity, increases toward the break; when the ratio of forces approaches unity, the profile turns over, just before filament breakup is achieved. The timescales and relative roles of gravity to rotation indicate that the structure is stable, long lived ($\sim$ a few times 6 Myr), and undergoing outside-in evolution. Moreover, this filament has practically no star formation, a perpendicular Planck plane-of-the-sky (POS) magnetic field morphology, and POS "zig-zag" morphology, which together with the rotation profile lead to the suggestion that the 3D shape is a corkscrew filament with a helical magnetic field. These results, combined with results in Orion and G035.39-00.33, suggest evolution toward higher densities as rotating filaments shed angular momentum. Thus, magnetic fields may be an essential feature of high-mass (M $\sim10^5$ M$_{\odot}$) cloud filament evolution toward cluster formation.
△ Less
Submitted 17 December, 2020; v1 submitted 21 October, 2020;
originally announced October 2020.
-
Association of IceCube neutrinos with radio sources observed at Owens Valley and Metsähovi Radio Observatories
Authors:
T. Hovatta,
E. Lindfors,
S. Kiehlmann,
W. Max-Moerbeck,
M. Hodges,
I. Liodakis,
A. Lahteenmaki,
T. J. Pearson,
A. C. S. Readhead,
R. A. Reeves,
S. Suutarinen,
J. Tammi,
M. Tornikoski
Abstract:
Identifying the most likely sources for high-energy neutrino emission has been one of the main topics in high-energy astrophysics ever since the first observation of high-energy neutrinos by the IceCube Neutrino Observatory. Active galactic nuclei with relativistic jets, also known as blazars, have been considered to be one of the main candidates because of their ability to accelerate particles to…
▽ More
Identifying the most likely sources for high-energy neutrino emission has been one of the main topics in high-energy astrophysics ever since the first observation of high-energy neutrinos by the IceCube Neutrino Observatory. Active galactic nuclei with relativistic jets, also known as blazars, have been considered to be one of the main candidates because of their ability to accelerate particles to high energies. We study the connection between radio emission and IceCube neutrino events using data from the Owens Valley Radio Observatory (OVRO) and Metsähovi Radio Observatory blazar monitoring programs. We identify sources in our radio monitoring sample that are positionally consistent with IceCube high-energy neutrino events. We estimate their mean flux density and variability amplitudes around the neutrino arrival time, and compare these with values from random samples to establish the significance of our results. We find radio source associations within our samples with 15 high-energy neutrino events detected by IceCube. Nearly half of the associated sources are not detected in the $γ$-ray energies, but their radio variability properties and Doppler boosting factors are similar to the $γ$-ray detected objects in our sample, meaning that they could still be potential neutrino emitters. We find that the number of strongly flaring objects in our statistically complete OVRO samples is unlikely to be a random coincidence (at $2σ$ level). Based on our results, we conclude that although it is clear that not all neutrino events are associated with strong radio flaring blazars, observations of large-amplitude radio flares in a blazar at the same time as a neutrino event are unlikely to be a random coincidence.
△ Less
Submitted 15 April, 2021; v1 submitted 22 September, 2020;
originally announced September 2020.
-
Multi-wavelength characterization of the blazar S5~0716+714 during an unprecedented outburst phase
Authors:
MAGIC Collaboration,
M. L. Ahnen,
S. Ansoldi,
L. A. Antonelli,
C. Arcaro,
D. Baack,
A. Babić,
B. Banerjee,
P. Bangale,
U. Barres de Almeida,
J. A. Barrio,
J. Becerra González,
W. Bednarek,
E. Bernardini,
R. Ch. Berse,
A. Berti,
W. Bhattacharyya,
A. Biland,
O. Blanch,
G. Bonnoli,
R. Carosi,
A. Carosi,
G. Ceribella,
A. Chatterjee,
S. M. Colak
, et al. (165 additional authors not shown)
Abstract:
The BL Lac object S5~0716+714, a highly variable blazar, underwent an impressive outburst in January 2015 (Phase A), followed by minor activity in February (Phase B). The MAGIC observations were triggered by the optical flux observed in Phase A, corresponding to the brightest ever reported state of the source in the R-band. The comprehensive dataset collected is investigated in order to shed light…
▽ More
The BL Lac object S5~0716+714, a highly variable blazar, underwent an impressive outburst in January 2015 (Phase A), followed by minor activity in February (Phase B). The MAGIC observations were triggered by the optical flux observed in Phase A, corresponding to the brightest ever reported state of the source in the R-band. The comprehensive dataset collected is investigated in order to shed light on the mechanism of the broadband emission. Multi-wavelength light curves have been studied together with the broadband Spectral Energy Distributions (SEDs). The data set collected spans from radio, optical photometry and polarimetry, X-ray, high-energy (HE, 0.1 GeV < E < 100 GeV) with \textit{Fermi}-LAT to the very-high-energy (VHE, E>100 GeV) with MAGIC. The flaring state of Phase A was detected in all the energy bands, providing for the first time a multi-wavelength sample of simultaneous data from the radio band to the VHE. In the constructed SED the \textit{Swift}-XRT+\textit{NuSTAR} data constrain the transition between the synchrotron and inverse Compton components very accurately, while the second peak is constrained from 0.1~GeV to 600~GeV by \textit{Fermi}+MAGIC data. The broadband SED cannot be described with a one-zone synchrotron self-Compton model as it severely underestimates the optical flux in order to reproduce the X-ray to $γ$-ray data. Instead we use a two-zone model. The EVPA shows an unprecedented fast rotation. An estimation of the redshift of the source by combined HE and VHE data provides a value of $z = 0.31 \pm 0.02_{stats} \pm 0.05_{sys}$, confirming the literature value. The data show the VHE emission originating in the entrance and exit of a superluminal knot in and out a recollimation shock in the inner jet. A shock-shock interaction in the jet seems responsible for the observed flares and EVPA swing. This scenario is also consistent with the SED modelling.
△ Less
Submitted 1 July, 2018;
originally announced July 2018.
-
37 GHz observations of narrow-line Seyfert 1 galaxies
Authors:
A. Lähteenmäki,
E. Järvelä,
T. Hovatta,
M. Tornikoski,
D. L. Harrison,
M. López-Caniego,
W. Max-Moerbeck,
M. Mingaliev,
T. J. Pearson,
V. Ramakrishnan,
A. C. S. Readhead,
R. A. Reeves,
J. L. Richards,
Y. Sotnikova,
J. Tammi
Abstract:
Observations at 37 GHz, performed at Metsähovi Radio Observatory, are presented for a sample of 78 radio-loud and radio-quiet narrow-line Seyfert 1 (NLS1) galaxies, together with additional lower and higher frequency radio data from RATAN-600, Owens Valley Radio Observatory, and the Planck satellite. Most of the data have been gathered between February 2012 and April 2015 but for some sources even…
▽ More
Observations at 37 GHz, performed at Metsähovi Radio Observatory, are presented for a sample of 78 radio-loud and radio-quiet narrow-line Seyfert 1 (NLS1) galaxies, together with additional lower and higher frequency radio data from RATAN-600, Owens Valley Radio Observatory, and the Planck satellite. Most of the data have been gathered between February 2012 and April 2015 but for some sources even longer lightcurves exist. The detection rate at 37 GHz is around 19%, comparable to other populations of active galactic nuclei presumed to be faint at radio frequencies, such as BL Lac objects. Variability and spectral indices are determined for sources with enough detections. Based on the radio data, many NLS1 galaxies show a blazar-like radio spectra exhibiting significant variability. The spectra at a given time are often inverted or convex. The source of the high-frequency radio emission in NLS1 galaxies, detected at 37 GHz, is most probably a relativistic jet rather than star formation. Jets in NLS1 galaxies are therefore expected to be a much more common phenomenon than earlier assumed.
△ Less
Submitted 30 March, 2017;
originally announced March 2017.
-
Symmetric Achromatic Variability in Active Galaxies -- A Powerful New Gravitational Lensing Probe?
Authors:
H. K. Vedantham,
A. C. S. Readhead,
T. Hovatta,
T. J. Pearson,
R. D. Blandford,
M. A. Gurwell,
A. Lähteenmäki,
W. Max-Moerbeck,
V. Pavlidou,
V. Ravi,
R. A. Reeves,
J. L. Richards,
M. Tornikoski,
J. A. Zensus
Abstract:
We report the discovery of a rare new form of long-term radio variability in the light-curves of active galaxies (AG) --- Symmetric Achromatic Variability (SAV) --- a pair of opposed and strongly skewed peaks in the radio flux density observed over a broad frequency range. We propose that SAV arises through gravitational milli-lensing when relativistically moving features in AG jets move through g…
▽ More
We report the discovery of a rare new form of long-term radio variability in the light-curves of active galaxies (AG) --- Symmetric Achromatic Variability (SAV) --- a pair of opposed and strongly skewed peaks in the radio flux density observed over a broad frequency range. We propose that SAV arises through gravitational milli-lensing when relativistically moving features in AG jets move through gravitational lensing caustics created by $10^3-10^6 \;{\rm M}_{\odot}$ subhalo condensates or black holes located within intervening galaxies. The lower end of this mass range has been inaccessible with previous gravitational lensing techniques. This new interpretation of some AG variability can easily be tested and if it passes these tests, will enable a new and powerful probe of cosmological matter distribution on these intermediate mass scales, as well as provide, for the first time, micro-arcsecond resolution of the nuclei of AG --- a factor of 30--100 greater resolution than is possible with ground-based millimeter VLBI.
△ Less
Submitted 17 February, 2017;
originally announced February 2017.
-
The peculiar light-curve of J1415+1320: A case study in extreme scattering events
Authors:
H. K. Vedantham,
A. C. S. Readhead,
T. Hovatta,
L. V. E. Koopmans,
T. J. Pearson,
R. D. Blandford,
M. A. Gurwell,
A. Lähteenmäki,
W. Max-Moerbeck,
V. Pavlidou,
V. Ravi,
R. A. Reeves,
J. L. Richards,
M. Tornikoski,
J. A. Zensus
Abstract:
The radio light-curve of J1415+1320 (PKS 1413+135) shows time-symmetric and recurring U-shaped features across the cm-wave and mm-wave bands. The symmetry of these features points to lensing by an intervening object as the cause. U-shaped events in radio light curves in the cm-wave band have previously been attributed to Extreme scattering events (ESE). ESEs are thought to be the result of lensing…
▽ More
The radio light-curve of J1415+1320 (PKS 1413+135) shows time-symmetric and recurring U-shaped features across the cm-wave and mm-wave bands. The symmetry of these features points to lensing by an intervening object as the cause. U-shaped events in radio light curves in the cm-wave band have previously been attributed to Extreme scattering events (ESE). ESEs are thought to be the result of lensing by compact plasma structures in the Galactic interstellar medium, but the precise nature of these plasma structures remains unknown. Since the strength of a plasma lens evolves with wavelength $λ$ as $λ^2$, the presence of correlated variations at over a wide wavelength range casts doubt on the canonical ESE interpretation for J1415+1320. In this paper, we critically examine the evidence for plasma lensing in J1415+1320. We compute limits on the lensing strength, and the associated free-free opacity of the putative plasma lenses. We compare the observed and model ESE light curves, and also derive a lower limit on the lens distance based on the effects of parallax due to the Earth's orbit around the Sun. We conclude that plasma lensing is not a viable interpretation for J1415+1320's light curves and that symmetric U-shaped features in the radio light curves of extragalactic sources do not present {\em prima facie} evidence for ESEs. The methodology presented here is generic enough to be applicable to any plasma lensing candidate.
△ Less
Submitted 17 February, 2017;
originally announced February 2017.
-
The awakening of the gamma-ray narrow-Line Seyfert 1 galaxy PKS 1502+036
Authors:
F. D'Ammando,
M. Orienti,
J. Finke,
T. Hovatta,
M. Giroletti,
W. Max-Moerbeck,
T. J. Pearson,
A. C. S. Readhead,
R. A. Reeves,
J. L. Richards
Abstract:
After a long low-activity period, a gamma-ray flare from the narrow-line Seyfert 1 PKS 1502+036 (z=0.4089) was detected by the Large Area Telescope (LAT) on board Fermi in 2015. On 2015 December 20 the source reached a daily peak flux, in the 0.1-300 GeV band, of (93 $\pm$ 19) $\times$10$^{-8}$ ph cm$^{-2}$ s$^{-1}$, attaining a flux of (237 $\pm$ 71) $\times$10$^{-8}$ ph cm$^{-2}$ s$^{-1}$ on 3-h…
▽ More
After a long low-activity period, a gamma-ray flare from the narrow-line Seyfert 1 PKS 1502+036 (z=0.4089) was detected by the Large Area Telescope (LAT) on board Fermi in 2015. On 2015 December 20 the source reached a daily peak flux, in the 0.1-300 GeV band, of (93 $\pm$ 19) $\times$10$^{-8}$ ph cm$^{-2}$ s$^{-1}$, attaining a flux of (237 $\pm$ 71) $\times$10$^{-8}$ ph cm$^{-2}$ s$^{-1}$ on 3-hr time-scales, which corresponds to an isotropic luminosity of (7.3 $\pm$ 2.1) $\times$10$^{47}$ erg/s. The gamma-ray flare was not accompanied by significant spectral changes. We report on multi-wavelength radio-to-gamma-ray observations of PKS 1502+036 during 2008 August-2016 March by Fermi-LAT, Swift, XMM-Newton, Catalina Real-Time Transient Survey, and the Owens Valley Radio Observatory (OVRO). An increase in activity was observed on 2015 December 22 by Swift in optical, UV, and X-rays. The OVRO 15 GHz light curve reached the highest flux density observed from this source on 2016 January 12, indicating a delay of about three weeks between the gamma-ray and 15 GHz emission peaks. This suggests that the gamma-ray emitting region is located beyond the broad line region. We compared the spectral energy distribution (SED) of an average activity state with that of the flaring state. The two SED, with the high-energy bump modelled as an external Compton component with seed photons from a dust torus, could be fitted by changing the electron distribution parameters as well as the magnetic field. The fit of the disc emission during the average state constrains the black hole mass to values lower than 10$^8$ solar masses. The SED, high-energy emission mechanisms, and gamma-ray properties of the source resemble those of a flat spectrum radio quasar.
△ Less
Submitted 15 September, 2016;
originally announced September 2016.
-
Gamma rays from the quasar PKS 1441+25: story of an escape
Authors:
A. U. Abeysekara,
S. Archambault,
A. Archer,
T. Aune,
A. Barnacka,
W. Benbow,
R. Bird,
J. Biteau,
J. H. Buckley,
V. Bugaev,
J. V Cardenzana,
M. Cerruti,
X. Chen,
J. L. Christiansen,
L. Ciupik,
M. P. Connolly,
P. Coppi,
W. Cui,
H. J. Dickinson,
J. Dumm,
J. D. Eisch,
M. Errando,
A. Falcone,
Q. Feng,
J. P. Finley
, et al. (81 additional authors not shown)
Abstract:
Outbursts from gamma-ray quasars provide insights on the relativistic jets of active galactic nuclei and constraints on the diffuse radiation fields that fill the Universe. The detection of significant emission above 100 GeV from a distant quasar would show that some of the radiated gamma rays escape pair-production interactions with low-energy photons, be it the extragalactic background light (EB…
▽ More
Outbursts from gamma-ray quasars provide insights on the relativistic jets of active galactic nuclei and constraints on the diffuse radiation fields that fill the Universe. The detection of significant emission above 100 GeV from a distant quasar would show that some of the radiated gamma rays escape pair-production interactions with low-energy photons, be it the extragalactic background light (EBL), or the radiation near the supermassive black hole lying at the jet's base. VERITAS detected gamma-ray emission up to 200 GeV from PKS 1441+25 (z=0.939) during April 2015, a period of high activity across all wavelengths. This observation of PKS 1441+25 suggests that the emission region is located thousands of Schwarzschild radii away from the black hole. The gamma-ray detection also sets a stringent upper limit on the near-ultraviolet to near-infrared EBL intensity, suggesting that galaxy surveys have resolved most, if not all, of the sources of the EBL at these wavelengths.
△ Less
Submitted 14 December, 2015;
originally announced December 2015.
-
The Q/U Imaging ExperimenT: Polarization Measurements of the Galactic Plane at 43 and 95 GHz
Authors:
QUIET Collaboration,
T. M. Ruud,
U. Fuskeland,
I. K. Wehus,
M. Vidal,
D. Araujo,
C. Bischoff,
I. Buder,
Y. Chinone,
K. Cleary,
R. N. Dumoulin,
A. Kusaka,
R. Monsalve,
S. K. Naess,
L. B. Newburgh,
R. A. Reeves,
J. T. L. Zwart,
L. Bronfman,
R. D. Davies,
R. Davis,
C. Dickinson,
H. K. Eriksen,
T. Gaier,
J. O. Gundersen,
M. Hasegawa
, et al. (17 additional authors not shown)
Abstract:
We present polarization observations of two Galactic plane fields centered on Galactic coordinates (l,b)=(0 deg,0 deg) and (329 deg, 0 deg) at Q- (43 GHz) and W-band (95 GHz), covering between 301 and 539 square degrees depending on frequency and field. These measurements were made with the QUIET instrument between 2008 October and 2010 December, and include a total of 1263 hours of observations.…
▽ More
We present polarization observations of two Galactic plane fields centered on Galactic coordinates (l,b)=(0 deg,0 deg) and (329 deg, 0 deg) at Q- (43 GHz) and W-band (95 GHz), covering between 301 and 539 square degrees depending on frequency and field. These measurements were made with the QUIET instrument between 2008 October and 2010 December, and include a total of 1263 hours of observations. The resulting maps represent the deepest large-area Galactic polarization observations published to date at the relevant frequencies with instrumental rms noise varying between 1.8 and 2.8 uK deg, 2.3-6 times deeper than corresponding WMAP and Planck maps. The angular resolution is 27.3' and 12.8' FWHM at Q- and W-band, respectively. We find excellent agreement between the QUIET and WMAP maps over the entire fields, and no compelling evidence for significant residual instrumental systematic errors in either experiment, whereas the Planck 44 GHz map deviates from these in a manner consistent with reported systematic uncertainties for this channel. We combine QUIET and WMAP data to compute inverse-variance-weighted average maps, effectively retaining small angular scales from QUIET and large angular scales from WMAP. From these combined maps, we derive constraints on several important astrophysical quantities, including a robust detection of polarized synchrotron spectral index steepening of ~0.2 off the plane, as well as the Faraday rotation measure toward the Galactic center (RM=-4000 +/- 200 rad m^-2), all of which are consistent with previously published results. Both the raw QUIET and the co-added QUIET+WMAP maps are made publicly available together with all necessary ancillary information.
△ Less
Submitted 11 August, 2015;
originally announced August 2015.