-
The Radio Spectra of High Luminosity Compact Symmetric Objects (CSO-2s): Implications for Studies of Compact Jetted Active Galactic Nuclei
Authors:
P. V. de la Parra,
A. C. S Readhead,
T. Herbig,
S. Kiehlmann,
M. L. Lister,
V. Pavlidou,
R. A. Reeves,
A. Siemiginowska,
A. G. Sullivan,
T. Surti,
A. Synani,
K. Tassis,
G. B. Taylor,
P. N. Wilkinson,
M. F. Aller,
R. D. Blandford,
N. Globus,
C. R. Lawrence,
B. Molina,
S. O'Neill,
T. J. Pearson
Abstract:
This paper addresses, for the first time, a key aspect of the phenomenology of Compact Symmetric Objects (CSOs) -- the characteristics of their radio spectra. We present a radio-spectrum description of a complete sample of high luminosity CSOs (CSO-2s), which shows that they exhibit the \textit{complete} range of spectral types, including flat-spectrum sources ($α\ge -0.5$), steep-spectrum sources…
▽ More
This paper addresses, for the first time, a key aspect of the phenomenology of Compact Symmetric Objects (CSOs) -- the characteristics of their radio spectra. We present a radio-spectrum description of a complete sample of high luminosity CSOs (CSO-2s), which shows that they exhibit the \textit{complete} range of spectral types, including flat-spectrum sources ($α\ge -0.5$), steep-spectrum sources ($α< -0.5$), and peaked-spectrum sources. We show that there is no clear correlation between spectral type and size, but there is a correlation between the high-frequency spectral index and both object type and size. We also show that, to avoid biasing the data and to understand the various classes of jetted-AGN involved, the complete range of spectral types should be included in studying the general phenomenology of compact jetted-AGN, and that complete samples must be used, selected over a wide range of frequencies. We discuss examples that demonstrate these points. We find that the high-frequency spectral indices of CSO-2s span $-1.3 <α_{\rm hi} < -0.3$, and hence that radio spectral signatures cannot be used to discriminate definitively between CSO-2s, binary galactic nuclei, and millilensed objects, unless they have $α_{\rm hi} >-0.3$.
△ Less
Submitted 23 August, 2024;
originally announced August 2024.
-
PKS~J0805$-$0111: A Second Owens Valley Radio Observatory Blazar Showing Highly Significant Sinusoidal Radio Variability -- The Tip of the Iceberg
Authors:
P. V. de la Parra,
S. Kiehlmann,
P. Mroz,
A. C. S. Readhead,
A. Synani,
M. C. Begelman,
R. D. Blandford,
Y. Ding,
F. Harrison,
I. Liodakis,
W. Max-Moerbeck,
V. Pavlidou,
R. Reeves,
M. Vallisneri,
M. F. Aller,
M. J. Graham,
T. Hovatta,
C. R. Lawrence,
T. J. W. Lazio,
A. A. Mahabal,
B. Molina,
S. O'Neill,
T. J. Pearson,
V. Ravi,
K. Tassis
, et al. (1 additional authors not shown)
Abstract:
Owens Valley Radio Observatory (OVRO) observations of supermassive black hole binary (SMBHB) candidate PKS~2131$-$021 revealed, for the first time, six likely characteristics of the phenomenology exhibited by SMBHB in blazars, of which the most unexpected and critical is sinusoidal flux density variations. We have now identified a second blazar, PKS~J0805$-$0111, showing significant sinusoidal var…
▽ More
Owens Valley Radio Observatory (OVRO) observations of supermassive black hole binary (SMBHB) candidate PKS~2131$-$021 revealed, for the first time, six likely characteristics of the phenomenology exhibited by SMBHB in blazars, of which the most unexpected and critical is sinusoidal flux density variations. We have now identified a second blazar, PKS~J0805$-$0111, showing significant sinusoidal variations, with an observed period that translates to $1.422 \pm 0.005$ yr in the rest frame of the $z = 1.388$ object. We generate $10^6$ simulated light curves to reproduce the radio variability characteristics of PKS~J0805$-$0111, and show that the global probability, considering the \textit{look-elsewhere effect}, indicates that the observed periodicity can be attributed to the red noise tail of the power spectral density, with a $p_0$ value of $7.8 \times 10^{-5}$ (i.e. 3.78$σ$). PKS J0805$-$0111 displays all six characteristics observed in PKS 2131$-$021. Taking into account the well-defined OVRO sample size, the false positive probability $\sim 0.22$, but the rare behavior makes this a strong SMBHB candidate. The discovery of a second SMBHB candidate exhibiting these rare characteristics reveals that PKS~2131$-$021 is not a unique, isolated case. With these two strong cases we are clearly seeing only the tip of the iceberg. We estimate that the number of SMBHB candidates amongst blazars $\sim$ 1 in 100.
△ Less
Submitted 5 August, 2024;
originally announced August 2024.
-
PKS 2131-021 -- Discovery of Strong Coherent Sinusoidal Variations from Radio to Optical Frequencies: Compelling Evidence for a Blazar Supermassive Black Hole Binary
Authors:
Sebastian Kiehlmann,
Philipe Vergara De La Parra,
Andrew Sullivan,
A. Synani,
Ioannis Liodakis,
Anthony Readhead,
Matthew Graham,
Mitchell Begelman,
Roger Blandford,
Katerina Chatziioannou,
Yuanze Ding,
Fiona Harrison,
D. Homan,
Talvikki Hovatta,
Shrinivas Kulkarni,
Matthew Lister,
Roberto Maiolino,
Walter Max-Moerbeck,
B. Molina,
Przemyslaw Mroz,
Christopher O'Dea,
Vasiliki Pavlidou,
Timothy J. Pearson,
Margo Aller,
C. Lawrence
, et al. (8 additional authors not shown)
Abstract:
Haystack and Owens Valley Radio Observatory (OVRO) observations recently revealed strong sinusoidal total flux density variations that maintained coherence between 1975 and 2021 in the blazar PKS 2131-021 ($z=1.283)$. This was interpreted as possible evidence of a supermassive black hole binary (SMBHB). Extended observations through 2023 show coherence over 47.9~years, with an observed period…
▽ More
Haystack and Owens Valley Radio Observatory (OVRO) observations recently revealed strong sinusoidal total flux density variations that maintained coherence between 1975 and 2021 in the blazar PKS 2131-021 ($z=1.283)$. This was interpreted as possible evidence of a supermassive black hole binary (SMBHB). Extended observations through 2023 show coherence over 47.9~years, with an observed period $P_\textrm{15 GHz}=(1739.3 \pm 1.2) \, {\rm days}$. We reject, with $p$-value = $5.3 \times 10^{-7}$, the hypothesis that the variations are due to random fluctuations in the red noise tail of the power spectral density. There is clearly a constant-period physical phenomenon in PKS 2131-021 producing coherent intermittent sinusoidal flux density variations. We find the coherent sinusoidal intensity variations extend from below 2.7 GHz to optical frequencies, from which we derive an observed period $P_\textrm{optical}=(1764 \pm 36)$ days. Across this broad frequency range there is a monotonic phase shift in the sinusoidal variations with frequency. The same coherent periodicity is possibly also observed at $γ$-ray energies. The importance of well-vetted SMBHB candidates to searches for gravitational waves is pointed out. We estimate the fraction of blazars that are SMBHB candidates to be $>1$ in 100. Thus monitoring programs covering tens of thousands of blazars could discover hundreds of SMBHB candidates.
△ Less
Submitted 12 July, 2024;
originally announced July 2024.
-
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.
-
Systems design, assembly, integration and lab testing of WALOP-South Polarimeter
Authors:
Siddharth Maharana,
A. N. Ramaprakash,
Chaitanya Rajarshi,
Pravin Khodade,
Bhushan Joshi,
Pravin Chordia,
Abhay Kohok,
Ramya M. Anche,
Deepa Modi,
John A. Kypriotakis,
Amit Deokar,
Aditya Kinjawadekar,
Stephen B. Potter,
Dmitry Blinov,
Hans Kristian Eriksen,
Myrto Falalaki,
Hitesh Gajjar,
Tuhin Ghosh,
Eirik Gjerløw,
Sebastain Kiehlmann,
Ioannis Liodakis,
Nikolaos Mandarakas,
Georgia V. Panopoulou,
Vasiliki Pavlidou,
Timothy J. Pearson
, et al. (6 additional authors not shown)
Abstract:
Wide-Area Linear Optical Polarimeter (WALOP)-South is the first wide-field and survey-capacity polarimeter in the optical wavelengths. On schedule for commissioning in 2024, it will be mounted on the 1 m SAAO telescope in Sutherland Observatory, South Africa to undertake the PASIPHAE sky survey. PASIPHAE program will create the first polarimetric sky map in the optical wavelengths, spanning more t…
▽ More
Wide-Area Linear Optical Polarimeter (WALOP)-South is the first wide-field and survey-capacity polarimeter in the optical wavelengths. On schedule for commissioning in 2024, it will be mounted on the 1 m SAAO telescope in Sutherland Observatory, South Africa to undertake the PASIPHAE sky survey. PASIPHAE program will create the first polarimetric sky map in the optical wavelengths, spanning more than 2000 square degrees of the southern Galactic region. The innovative design of WALOP-South will enable it to measure the linear polarization (Stokes parameters $q$ and $u$), in a single exposure, of all sources in a field of view (FoV) of $35\times35$ arcminutes-squared in the SDSS-r broadband and narrowband filters between 500-750 nm with 0.1 % polarization accuracy.
The unique goals of the instrument place very stringent systems engineering goals, including on the performance of the optical, polarimetric, optomechanical, and electronic subsystems. All the subsystems have been designed carefully to meet the overall instrument performance goals.
As of May 2024, all the instrument optical and mechanical subsystems have been assembled and are currently getting tested and integrated. The complete testing and characterization of the instrument in the lab is expected to be completed by August 2024.
In this paper, we will present (a) the design and development of the entire instrument and its major subsystems, focusing on the opto-mechanical design which has not been reported before, and (b) assembly and integration of the instrument in the lab and early results from lab characterization of the instrument.
△ Less
Submitted 27 June, 2024;
originally announced June 2024.
-
The Crucial First Step in the Discovery of Millisecond Pulsars
Authors:
A. C. S. Readhead
Abstract:
Millisecond pulsars, like pulsars, have led to major advances in many areas of astronomy and physics. The discovery in 2023 of a cosmological gravitational wave (GW) stochastic background using millisecond pulsar timing arrays has focused attention on the importance of millisecond pulsars, to both multi-messenger astronomy and cosmology, and for identifying the origin of the GW stochastic backgrou…
▽ More
Millisecond pulsars, like pulsars, have led to major advances in many areas of astronomy and physics. The discovery in 2023 of a cosmological gravitational wave (GW) stochastic background using millisecond pulsar timing arrays has focused attention on the importance of millisecond pulsars, to both multi-messenger astronomy and cosmology, and for identifying the origin of the GW stochastic background, which is hypothesized to be due to supermassive black hole binaries (SMBHBs). Unlike pulsars, however, for which the details of the discovery are well-known, those of millisecond pulsars are not well known. In particular, the details of the first crucial step in the discovery of millisecond pulsars, namely the discovery of interplanetary scintillation (IPS) in the radio source 4C 21.53, are known only to the author. This article presents a first-hand account of this crucial first step, which resulted ultimately in the discovery of millisecond pulses from this object. A brief description of interplanetary and interstellar scintillation and scattering is given in the Appendix.
△ Less
Submitted 24 June, 2024;
originally announced June 2024.
-
COMAP Pathfinder -- Season 2 results III. Implications for cosmic molecular gas content at "Cosmic Half-past Eleven"
Authors:
D. T. Chung,
P. C. Breysse,
K. A. Cleary,
D. A. Dunne,
J. G. S. Lunde,
H. Padmanabhan,
N. -O. Stutzer,
D. Tolgay,
J. R. Bond,
S. E. Church,
H. K. Eriksen,
T. Gaier,
J. O. Gundersen,
S. E. Harper,
A. I. Harris,
R. Hobbs,
H. T. Ihle,
J. Kim,
J. W. Lamb,
C. R. Lawrence,
N. Murray,
T. J. Pearson,
L. Philip,
A. C. S. Readhead,
T. J. Rennie
, et al. (2 additional authors not shown)
Abstract:
The Carbon monOxide Mapping Array Project (COMAP) Pathfinder survey continues to demonstrate the feasibility of line-intensity mapping using high-redshift carbon monoxide (CO) line emission traced at cosmological scales. The latest COMAP Pathfinder power spectrum analysis is based on observations through the end of Season 2, covering the first three years of Pathfinder operations. We use our lates…
▽ More
The Carbon monOxide Mapping Array Project (COMAP) Pathfinder survey continues to demonstrate the feasibility of line-intensity mapping using high-redshift carbon monoxide (CO) line emission traced at cosmological scales. The latest COMAP Pathfinder power spectrum analysis is based on observations through the end of Season 2, covering the first three years of Pathfinder operations. We use our latest constraints on the CO(1-0) line-intensity power spectrum at $z\sim3$ to update corresponding constraints on the cosmological clustering of CO line emission and thus the cosmic molecular gas content at a key epoch of galaxy assembly. We first mirror the COMAP Early Science interpretation, considering how Season 2 results translate to limits on the shot noise power of CO fluctuations and the bias of CO emission as a tracer of the underlying dark matter distribution. The COMAP Season 2 results place the most stringent limits on the CO tracer bias to date, at $\langle{Tb}\rangle<4.8$ $μ$K. These limits narrow the model space significantly compared to previous CO line-intensity mapping results while maintaining consistency with small-volume interferometric surveys of resolved line candidates. The results also express a weak preference for CO emission models used to guide fiducial forecasts from COMAP Early Science, including our data-driven priors. We also consider directly constraining a model of the halo-CO connection, and show qualitative hints of capturing the total contribution of faint CO emitters through the improved sensitivity of COMAP data. With continued observations and matching improvements in analysis, the COMAP Pathfinder remains on track for a detection of cosmological clustering of CO emission.
△ Less
Submitted 14 June, 2024; v1 submitted 11 June, 2024;
originally announced June 2024.
-
COMAP Pathfinder -- Season 2 results II. Updated constraints on the CO(1-0) power spectrum
Authors:
N. -O. Stutzer,
J. G. S. Lunde,
P. C. Breysse,
D. T. Chung,
K. A. Cleary,
D. A. Dunne,
H. K. Eriksen,
H. T. Ihle,
H. Padmanabhan,
D. Tolgay,
I. K. Wehus,
J. R. Bond,
S. E. Church,
T. Gaier,
J. O. Gundersen,
A. I. Harris,
S. E. Harper,
R. Hobbs,
J. Kim,
J. W. Lamb,
C. R. Lawrence,
N. Murray,
T. J. Pearson,
L. Philip,
A. C. S. Readhead
, et al. (2 additional authors not shown)
Abstract:
We present updated constraints on the cosmological 3D power spectrum of carbon monoxide CO(1-0) emission in the redshift range $2.4$-$3.4$. The constraints are derived from the two first seasons of Carbon monOxide Mapping Array Project (COMAP) Pathfinder line-intensity mapping observations aiming to trace star-formation during the Epoch of Galaxy Assembly. These results improve on the previous Ear…
▽ More
We present updated constraints on the cosmological 3D power spectrum of carbon monoxide CO(1-0) emission in the redshift range $2.4$-$3.4$. The constraints are derived from the two first seasons of Carbon monOxide Mapping Array Project (COMAP) Pathfinder line-intensity mapping observations aiming to trace star-formation during the Epoch of Galaxy Assembly. These results improve on the previous Early Science (ES) results through both increased data volume and improved data processing methodology. On the methodological side, we now perform cross-correlations between groups of detectors (''feed-groups''), as opposed to cross-correlations between single feeds, and this new feed-group pseudo power spectrum (FGPXS) is constructed to be more robust against systematic effects. In terms of data volume, the effective mapping speed is significantly increased due to an improved observational strategy as well as better data selection methodology. The updated spherically- and field-averaged FGPXS, $\tilde{C}(k)$, is consistent with zero, at a probability-to-exceed of around $34\,\%$, with an excess of $2.7\,σ$ in the most sensitive bin. Our power spectrum estimate is about an order of magnitude more sensitive in our six deepest bins across ${0.09\,\mathrm{Mpc}^{-1} < k < 0.73\,\mathrm{Mpc}^{-1}}$, as compared to the feed-feed pseudo power spectrum (FPXS) of COMAP ES. Each of these bins individually constrains the CO power spectrum to ${kP_\mathrm{CO}(k)< 2400-4900\,\mathrm{μK^2 Mpc^{2}}}$ at $95\,\%$ confidence. To monitor potential contamination from residual systematic effects, we analyze a set of 312 difference-map null tests and find that these are consistent with the instrumental noise prediction. In sum, these results provide the strongest direct constraints on the cosmological 3D CO(1-0) power spectrum published to date.
△ Less
Submitted 14 June, 2024; v1 submitted 11 June, 2024;
originally announced June 2024.
-
COMAP Pathfinder -- Season 2 results I. Improved data selection and processing
Authors:
J. G. S. Lunde,
N. -O. Stutzer,
P. C. Breysse,
D. T. Chung,
K. A. Cleary,
D. A. Dunne,
H. K. Eriksen,
S. E. Harper,
H. T. Ihle,
J. W. Lamb,
T. J. Pearson,
L. Philip,
I. K. Wehus,
D. P. Woody,
J. R. Bond,
S. E. Church,
T. Gaier,
J. O. Gundersen,
A. I. Harris,
R. Hobbs,
J. Kim,
C. R. Lawrence,
N. Murray,
H. Padmanabhan,
A. C. S. Readhead
, et al. (2 additional authors not shown)
Abstract:
The CO Mapping Array Project (COMAP) Pathfinder is performing line intensity mapping of CO emission to trace the distribution of unresolved galaxies at redshift $z \sim 3$. We present an improved version of the COMAP data processing pipeline and apply this to the first two seasons of observations. This analysis improves on the COMAP Early Science (ES) results in several key aspects. On the observa…
▽ More
The CO Mapping Array Project (COMAP) Pathfinder is performing line intensity mapping of CO emission to trace the distribution of unresolved galaxies at redshift $z \sim 3$. We present an improved version of the COMAP data processing pipeline and apply this to the first two seasons of observations. This analysis improves on the COMAP Early Science (ES) results in several key aspects. On the observational side, all second season scans were made in constant-elevation mode, after noting that the previous Lissajous scans were associated with increased systematic errors; those scans accounted for 50% of the total Season 1 data volume. Secondly, all new observations were restricted to an elevation range of 35-65 degrees, to minimize sidelobe ground pickup. On the data processing side, more effective data cleaning in both the time- and map-domain has allowed us to eliminate all data-driven power spectrum-based cuts. This increases the overall data retention and reduces the risk of signal subtraction bias. On the other hand, due to the increased sensitivity, two new pointing-correlated systematic errors have emerged, and we introduce a new map-domain PCA filter to suppress these. Subtracting only 5 out of 256 PCA modes, we find that the standard deviation of the cleaned maps decreases by 67% on large angular scales, and after applying this filter, the maps appear consistent with instrumental noise. Combining all these improvements, we find that each hour of raw Season 2 observations yields on average 3.2 times more cleaned data compared to ES analysis. Combining this with the increase in raw observational hours, the effective amount of data available for high-level analysis is a factor of 8 higher than in ES. The resulting maps have reached an uncertainty of $25$-$50\,μK$ per voxel, providing by far the strongest constraints on cosmological CO line emission published to date.
△ Less
Submitted 14 June, 2024; v1 submitted 11 June, 2024;
originally announced June 2024.
-
A $γ$-Ray Emitting Blazar at Redshift 3.64: Fermi-LAT and OVRO Observations of PKS 0201+113
Authors:
Hai Lei,
Ying-Kang Zhang,
Xiong Jiang,
Sebastian Kiehlmann,
Anthony C. S. Readhead,
Liang Chen,
Neng-Hui Liao,
Tao An
Abstract:
High-redshift ($z >3$) $γ$-ray blazars are rare, but they are crucial for our understanding of jet evolution, $γ$-ray production and propagation, and the growth of supermassive black holes in the early universe. A new analysis of Fermi-LAT data reveals a significant (5$σ$), spectrally soft ($Γ\simeq$ 3.0) $γ$-ray source in a specific 4-month epoch, cospatial with PKS 0201+113 ($z$ = 3.64). Monitor…
▽ More
High-redshift ($z >3$) $γ$-ray blazars are rare, but they are crucial for our understanding of jet evolution, $γ$-ray production and propagation, and the growth of supermassive black holes in the early universe. A new analysis of Fermi-LAT data reveals a significant (5$σ$), spectrally soft ($Γ\simeq$ 3.0) $γ$-ray source in a specific 4-month epoch, cospatial with PKS 0201+113 ($z$ = 3.64). Monitoring of PKS 0201+113 at 15 GHz by the Owens Valley Radio Observatory 40 m Telescope from 2008 to 2023 shows a prominent flare that dominates the radio light curve. The maximum of the radio flare coincides with the $γ$-ray flare, strongly suggesting an association ($\textrm{p-value}=0.023$) between the $γ$-ray and the radio sources. PKS 0201+113 is only the third $γ$-ray blazar to be identified with $z> 3.5$, and it is the first such object to be identified by the detection of quasi-simultaneous $γ$-ray and radio flares. The jet properties of this peculiar blazar have been investigated. A detailed study of a two-zone leptonic model is presented that fits the broadband spectral energy distribution. An alternative scenario is also briefly discussed.
△ Less
Submitted 5 June, 2024;
originally announced June 2024.
-
The first degree-scale starlight-polarization-based tomography map of the magnetized interstellar medium
Authors:
V. Pelgrims,
N. Mandarakas,
R. Skalidis,
K. Tassis,
G. V. Panopoulou,
V. Pavlidou,
D. Blinov,
S. Kiehlmann,
S. E. Clark,
B. S. Hensley,
S. Romanopoulos,
A. Basyrov,
H. K. Eriksen,
M. Falalaki,
T. Ghosh,
E. Gjerløw,
J. A. Kypriotakis,
S. Maharana,
A. Papadaki,
T. J. Pearson,
S. B. Potter,
A. N. Ramaprakash,
A. C. S. Readhead,
I. K. Wehus
Abstract:
We present the first degree-scale tomography map of the dusty magnetized interstellar medium (ISM) from stellar polarimetry and distance measurements. We used the RoboPol polarimeter at Skinakas Observatory to conduct a survey of starlight polarization in a region of the sky of 4 square degrees. We propose a Bayesian method to decompose the stellar-polarization source field along the distance to i…
▽ More
We present the first degree-scale tomography map of the dusty magnetized interstellar medium (ISM) from stellar polarimetry and distance measurements. We used the RoboPol polarimeter at Skinakas Observatory to conduct a survey of starlight polarization in a region of the sky of 4 square degrees. We propose a Bayesian method to decompose the stellar-polarization source field along the distance to invert the 3D volume occupied by the observed stars. We used it to obtain the first 3D map of the dusty magnetized ISM. Specifically, we produced a tomography map of the orientation of the plane-of-sky (POS) component of the magnetic field threading the diffuse, dusty regions responsible for the stellar polarization. For the targeted region centered on Galactic coordinates $(l,b) \approx (103.3^\circ, 22.3^\circ)$, we identified several ISM clouds. Most of the lines of sight intersect more than one cloud. A very nearby component was detected in the foreground of a dominant component from which most of the polarization signal comes. Farther clouds, with a distance of up to 2~kpc, were similarly detected. Some of them likely correspond to intermediate-velocity clouds seen in HI spectra in this region of the sky. We found that the orientation of the POS component of the magnetic field changes along distance for most of the lines of sight. Our study demonstrates that starlight polarization data coupled to distance measures have the power to reveal the great complexity of the dusty magnetized ISM in 3D and, in particular, to provide local measurements of the POS component of the magnetic field. This demonstrates that the inversion of large data volumes, as expected from the PASIPHAE survey, will provide the necessary means to move forward in the modeling of the Galactic magnetic field and of the dusty magnetized ISM as a contaminant in observations of the cosmic microwave background polarization.
△ Less
Submitted 16 April, 2024;
originally announced April 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.
-
Owens Valley Radio Observatory monitoring of LS I +61°303 completes three cycles of the super-orbital modulation
Authors:
F. Jaron,
S. Kiehlmann,
A. C. S. Readhead
Abstract:
The high-mass X-ray binary LS I +61°303 is composed of a Be-type star and a compact object. The emission is variable and periodic across the electromagnetic spectrum, from radio to very high-energy gamma rays. The orbital period is ~26.5 d, and the source also features a super-orbital period with a value of ~4.6 years. Long-term monitoring of the binary by the Owens Valley Radio Observatory (OVRO)…
▽ More
The high-mass X-ray binary LS I +61°303 is composed of a Be-type star and a compact object. The emission is variable and periodic across the electromagnetic spectrum, from radio to very high-energy gamma rays. The orbital period is ~26.5 d, and the source also features a super-orbital period with a value of ~4.6 years. Long-term monitoring of the binary by the Owens Valley Radio Observatory (OVRO) at 15 GHz has now completed 13.8 years, which corresponds to three full cycles of the super-orbital period. We performed a timing analysis of the OVRO radio light curve, and we also combined the OVRO data with the full archive of previous radio observations and computed the discrete autocorrelation function. The most powerful features in the periodogram of the OVRO data are two peaks at P1 = 26.49 +/- 0.05 d and P2 = 26.93 +/- 0.05 d. Our measurement of the long-term period is P_long = 1698 +/- 196 d. Dividing the OVRO data into three segments of equal length showed that the two periods, P1 and P2, are present in the periodogram of each of the consecutive long-term cycles. Our analysis of the full radio archive resulted in the detection of the same three periods, and the autocorrelation function showed a regular pattern, proving the stability of the super-orbital modulation. We report a possible systematic modulation of the radio flux density with a timescale of approximately 40 years that has so far remained unnoticed. The physical model of a relativistic jet whose mass loading is modulated with the orbital period P1 and is precessing with the slightly larger period P2, giving rise to a beating with period P_long, had previously been able to reproduce the radio and gigaelectron volt emission. The ongoing presence and the stability of the periodic signals imply that this model is still the most plausible explanation for the physical processes at work in this source. (abridged)
△ Less
Submitted 12 February, 2024;
originally announced February 2024.
-
Small-scale radio jets and tidal disruption events: A theory of high-luminosity compact symmetric objects
Authors:
Andrew G. Sullivan,
Roger D. Blandford,
Mitchell C. Begelman,
Mark Birkinshaw,
Anthony C. S. Readhead
Abstract:
Double lobe radio sources associated with active galactic nuclei represent one of the longest studied groups in radio astronomy. A particular sub-group of double radio sources comprises the compact symmetric objects (CSOs). CSOs are distinguished by their prominent double structure and sub-kpc total size. It has been argued that the vast majority of high-luminosity CSOs (CSO 2s) represent a distin…
▽ More
Double lobe radio sources associated with active galactic nuclei represent one of the longest studied groups in radio astronomy. A particular sub-group of double radio sources comprises the compact symmetric objects (CSOs). CSOs are distinguished by their prominent double structure and sub-kpc total size. It has been argued that the vast majority of high-luminosity CSOs (CSO 2s) represent a distinct class of active galactic nuclei with its own morphological structure and life-cycle. In this work, we present theoretical considerations regarding CSO 2s. We develop a semi-analytic evolutionary model, inspired by the results of large-scale numerical simulations of relativistic jets, that reproduces the features of the radio source population. We show that CSO 2s may be generated by finite energy injections and propose stellar tidal disruption events as a possible cause. We find that tidal disruption events of giant branch stars with masses $\gtrsim1$ M$_\odot$ can fuel these sources and discuss possible approaches to confirming this hypothesis. We predict that if the tidal disruption scenario holds, CSO 2s with sizes less than 400 pc should outnumber larger sources by more than a factor of $10$. Our results motivate future numerical studies to determine whether the scenarios we consider for fueling and source evolution can explain the observed radio morphologies. Multiwavelength observational campaigns directed at these sources will also provide critical insight into the origins of these objects, their environments, and their lifespans.
△ Less
Submitted 25 January, 2024;
originally announced January 2024.
-
Insights into the broad-band emission of the TeV blazar Mrk 501 during the first X-ray polarization measurements
Authors:
S. Abe,
J. Abhir,
V. A. Acciari,
A. Aguasca-Cabot,
I. Agudo,
T. Aniello,
S. Ansoldi,
L. A. Antonelli,
A. Arbet Engels,
C. Arcaro,
K. Asano,
A. Babić,
A. Baquero,
U. Barres de Almeida,
J. A. Barrio,
I. Batković,
A. Bautista,
J. Baxter,
J. Becerra González,
W. Bednarek,
E. Bernardini,
M. Bernardos,
J. Bernete,
A. Berti,
J. Besenrieder
, et al. (239 additional authors not shown)
Abstract:
We present the first multi-wavelength study of Mrk 501 including very-high-energy (VHE) gamma-ray observations simultaneous to X-ray polarization measurements from the Imaging X-ray Polarimetry Explorer (IXPE). We use radio-to-VHE data from a multi-wavelength campaign organized between 2022-03-01 and 2022-07-19. The observations were performed by MAGIC, Fermi-LAT, NuSTAR, Swift (XRT and UVOT), and…
▽ More
We present the first multi-wavelength study of Mrk 501 including very-high-energy (VHE) gamma-ray observations simultaneous to X-ray polarization measurements from the Imaging X-ray Polarimetry Explorer (IXPE). We use radio-to-VHE data from a multi-wavelength campaign organized between 2022-03-01 and 2022-07-19. The observations were performed by MAGIC, Fermi-LAT, NuSTAR, Swift (XRT and UVOT), and several instruments covering the optical and radio bands. During the IXPE pointings, the VHE state is close to the average behavior with a 0.2-1 TeV flux of 20%-50% the emission of the Crab Nebula. Despite the average VHE activity, an extreme X-ray behavior is measured for the first two IXPE pointings in March 2022 with a synchrotron peak frequency >1 keV. For the third IXPE pointing in July 2022, the synchrotron peak shifts towards lower energies and the optical/X-ray polarization degrees drop. The X-ray polarization is systematically higher than at lower energies, suggesting an energy-stratification of the jet. While during the IXPE epochs the polarization angle in the X-ray, optical and radio bands align well, we find a clear discrepancy in the optical and radio polarization angles in the middle of the campaign. We model the broad-band spectra simultaneous to the IXPE pointings assuming a compact zone dominating in the X-rays and VHE, and an extended zone stretching further downstream the jet dominating the emission at lower energies. NuSTAR data allow us to precisely constrain the synchrotron peak and therefore the underlying electron distribution. The change between the different states observed in the three IXPE pointings can be explained by a change of magnetization and/or emission region size, which directly connects the shift of the synchrotron peak to lower energies with the drop in polarization degree.
△ Less
Submitted 16 January, 2024;
originally announced January 2024.
-
Zero-polarization candidate regions for calibration of wide-field optical polarimeters
Authors:
N. Mandarakas,
G. V. Panopoulou,
V. Pelgrims,
S. B. Potter,
V. Pavlidou,
A. Ramaprakash,
K. Tassis,
D. Blinov,
S. Kiehlmann,
E. Koutsiona,
S. Maharana,
S. Romanopoulos,
R. Skalidis,
A. Vervelaki,
S. E. Clark,
J. A. Kypriotakis,
A. C. S. Readhead
Abstract:
Context. Calibration of optical polarimeters relies on the use of stars with negligible polarization (unpolarized standard stars) for determining the instrumental polarization zero-point. For wide-field polarimeters, calibration is often done by imaging the same star over multiple positions in the field of view - a process which is time-consuming. A more effective technique is to target fields con…
▽ More
Context. Calibration of optical polarimeters relies on the use of stars with negligible polarization (unpolarized standard stars) for determining the instrumental polarization zero-point. For wide-field polarimeters, calibration is often done by imaging the same star over multiple positions in the field of view - a process which is time-consuming. A more effective technique is to target fields containing multiple standard stars. While this method has been used for fields with highly polarized stars, there are no such sky regions with well-measured unpolarized standard stars. Aims. We aim to identify sky regions with tens of stars exhibiting negligible polarization, which are suitable for zero-point calibration of wide-field polarimeters. Methods. We selected stars in regions with extremely low reddening, located at high Galactic latitudes. We targeted four ~ 400 x 400 fields in the northern, and eight in the southern Equatorial hemisphere. Observations were carried out at the Skinakas Observatory and the South African Astronomical Observatory respectively. Results. We find two fields in the North and seven in the South with mean polarization lower than p < 0.1%. Conclusions. At least nine out of twelve fields can be used for zero-point calibration of wide-field polarimeters.
△ Less
Submitted 11 December, 2023;
originally announced December 2023.
-
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.
-
Multi-wavelength observations of the lensed quasar PKS 1830$-$211 during the 2019 $γ$-ray flare
Authors:
S. Vercellone,
I. Donnarumma,
C. Pittori,
F. Capitanio,
A. De Rosa,
L. Di Gesu,
S. Kiehlmann,
M. N. Iacolina,
P. A. Pellizzoni,
E. Egron,
L. Pacciani,
G. Piano,
S. Puccetti,
S. Righini,
G. Valente,
F. Verrecchia,
V. Vittorini,
M. Tavani,
E. Brocato,
A. W. Chen,
T. Hovatta,
A. Melis,
W. Max-Moerbeck,
D. Perrodin,
M. Pilia
, et al. (10 additional authors not shown)
Abstract:
PKS 1830$-$211 is a $γ$-ray emitting, high-redshift (z $= 2.507 \pm 0.002$), lensed flat-spectrum radio quasar. During the period mid-February to mid-April 2019, this source underwent a series of strong $γ$-ray flares that were detected by both AGILE-GRID and Fermi-LAT, reaching a maximum $γ$-ray flux of $F_{\rm E>100 MeV}\approx 2.3\times10^{-5}$ ph cm$^{-2}$ s$^{-1}$. Here we report on a coordin…
▽ More
PKS 1830$-$211 is a $γ$-ray emitting, high-redshift (z $= 2.507 \pm 0.002$), lensed flat-spectrum radio quasar. During the period mid-February to mid-April 2019, this source underwent a series of strong $γ$-ray flares that were detected by both AGILE-GRID and Fermi-LAT, reaching a maximum $γ$-ray flux of $F_{\rm E>100 MeV}\approx 2.3\times10^{-5}$ ph cm$^{-2}$ s$^{-1}$. Here we report on a coordinated campaign from both on-ground (Medicina, OVRO, REM, SRT) and orbiting facilities (AGILE, Fermi, INTEGRAL, NuSTAR, Swift, Chandra), with the aim of investigating the multi-wavelength properties of PKS 1830$-$211 through nearly simultaneous observations presented here for the first time. We find a possible break in the radio spectra in different epochs above 15 GHz, and a clear maximum of the 15 GHz data approximately 110 days after the $γ$-ray main activity periods. The spectral energy distribution shows a very pronounced Compton dominance (> 200) which challenges the canonical one-component emission model. Therefore we propose that the cooled electrons of the first component are re-accelerated to a second component by, e.g., kink or tearing instability during the $γ$-ray flaring periods. We also note that PKS 1830$-$211 could be a promising candidate for future observations with both Compton satellites (e.g., e-ASTROGAM) and Cherenkov arrays (CTAO) which will help, thanks to their improved sensitivity, in extending the data availability in energy bands currently uncovered.
△ Less
Submitted 13 November, 2023;
originally announced November 2023.
-
Optical circular polarization of blazar S4 0954+65 during high linear polarized states
Authors:
I. Liodakis,
E. Shablovinskaya,
D. Blinov,
S. S. Savchenko,
E. Malygin,
S. Kotov,
S. Kiehlmann,
A. C. S. Readhead,
S. B. Potter,
F. M. Rieger,
T. S. Grishina,
V. A. Hagen-Thorn,
E. N. Kopatskaya,
E. G. Larionova,
D. A. Morozova,
I. S. Troitskiy,
Y. V. Troitskaya,
A. A. Vasilyev,
A. V. Zhovtan,
G. A. Borman
Abstract:
Optical circular polarization observations can directly test the particle composition in black holes jets. Here we report on the first observations of the BL Lac type object S4 0954+65 in high linear polarized states. While no circular polarization was detected, we were able to place upper limits of <0.5% at the 99.7% confidence. Using a simple model and our novel optical circular polarization obs…
▽ More
Optical circular polarization observations can directly test the particle composition in black holes jets. Here we report on the first observations of the BL Lac type object S4 0954+65 in high linear polarized states. While no circular polarization was detected, we were able to place upper limits of <0.5% at the 99.7% confidence. Using a simple model and our novel optical circular polarization observations we can constrain the allowed parameter space for the magnetic field strength and composition of the emitting particles. Our results favor models that require magnetic field strengths of only a few Gauss and models where the jet composition is dominated by electron-positron pairs. We discuss our findings in the context of typical magnetic field strength requirements for blazar emission models.
△ Less
Submitted 6 November, 2023;
originally announced November 2023.
-
A Near Magnetic-to-kinetic Energy Equipartition Flare from the Relativistic Jet in AO 0235+164 during 2013-2019
Authors:
Whee Yeon Cheong,
Sang-Sung Lee,
Sang-Hyun Kim,
Sincheol Kang,
Jae Young Kim,
Bindu Rani,
Anthony C. S. Readhead,
Sebastian Kiehlmann,
Anne Lähteenmäki,
Merja Tornikoski,
Joni Tammi,
Venkatessh Ramakrishnan,
Iván Agudo,
Antonio Fuentes,
Efthalia Traianou,
Juan Escudero,
Clemens Thum,
Ioannis Myserlis,
Carolina Casadio,
Mark Gurwell
Abstract:
We present the multiwavelength flaring activity of the blazar AO 0235+164 during its recent active period from 2013 to 2019. From a discrete correlation function (DCF) analysis, we find a significant (>95%) correlation between radio and $γ$-ray light curves with flares at longer wavelengths following flares at shorter wavelengths. We identify a new jet component in 43 GHz VLBA data that was ejecte…
▽ More
We present the multiwavelength flaring activity of the blazar AO 0235+164 during its recent active period from 2013 to 2019. From a discrete correlation function (DCF) analysis, we find a significant (>95%) correlation between radio and $γ$-ray light curves with flares at longer wavelengths following flares at shorter wavelengths. We identify a new jet component in 43 GHz VLBA data that was ejected from the radio core on MJD $57246^{+26}_{-30}$ (2015 August 12), during the peak of the 2015 radio flare. From the analysis of the jet component, we derived a Doppler factor of $δ_{\rm var}=28.5\pm8.4$, a bulk Lorentz factor of $Γ=16.8^{+3.6}_{-3.1}$, and an intrinsic viewing angle of $θ_{\rm v}=1.42^{+1.07}_{-0.52}\textrm{ degrees}$. Investigation of the quasi-simultaneous radio data revealed a partially absorbed spectrum with the turnover frequency varying in the range of $10-70$ GHz and the peak flux density varying in the range of $0.7-4$ Jy. We find the synchrotron self-absorption magnetic field strength to be $B_{\rm SSA}=15.3^{+12.6}_{-14.0}$ mG at the peak of the 2015 radio flare, which is comparable to the equipartition magnetic field strength of $B_{\rm EQ}=43.6^{+10.6}_{-10.4}$ mG calculated for the same epoch. Additional analysis of the radio emission region in the relativistic jet of AO 0235+164 suggests that it did not significantly deviate from equipartition during its recent flaring activity.
△ Less
Submitted 24 October, 2023;
originally announced October 2023.
-
Multi-year characterisation of the broad-band emission from the intermittent extreme BL Lac 1ES~2344+514
Authors:
H. Abe,
S. Abe,
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,
A. Biland,
O. Blanch
, et al. (210 additional authors not shown)
Abstract:
The BL Lac 1ES 2344+514 is known for temporary extreme properties (e.g., a shift of the synchrotron SED peak energy $ν_{synch,p}$ above 1keV). While those extreme states were so far observed only during high flux levels, additional multi-year observing campaigns are required to achieve a coherent picture. Here, we report the longest investigation of the source from radio to VHE performed so far, f…
▽ More
The BL Lac 1ES 2344+514 is known for temporary extreme properties (e.g., a shift of the synchrotron SED peak energy $ν_{synch,p}$ above 1keV). While those extreme states were so far observed only during high flux levels, additional multi-year observing campaigns are required to achieve a coherent picture. Here, we report the longest investigation of the source from radio to VHE performed so far, focusing on a systematic characterisation of the intermittent extreme states. While our results confirm that 1ES 2344+514 typically exhibits $ν_{synch,p}>$1keV during elevated flux periods, we also find periods where the extreme state coincides with low flux activity. A strong spectral variability thus happens in the quiescent state, and is likely caused by an increase of the electron acceleration efficiency without a change in the electron injection luminosity. We also report a strong X-ray flare (among the brightest for 1ES 2344+514) without a significant shift of $ν_{synch,p}$. During this particular flare, the X-ray spectrum is among the softest of the campaign. It unveils complexity in the spectral evolution, where the common harder-when-brighter trend observed in BL Lacs is violated. During a low and hard X-ray state, we find an excess of the UV flux with respect to an extrapolation of the X-ray spectrum to lower energies. This UV excess implies that at least two regions contribute significantly to the infrared/optical/ultraviolet/X-ray emission. Using the simultaneous MAGIC, XMM-Newton, NuSTAR, and AstroSat observations, we argue that a region possibly associated with the 10 GHz radio core may explain such an excess. Finally, we investigate a VHE flare, showing an absence of simultaneous variability in the 0.3-2keV band. Using a time-dependent leptonic modelling, we show that this behaviour, in contradiction to single-zone scenarios, can instead be explained by a two-component model.
△ Less
Submitted 5 October, 2023;
originally announced October 2023.
-
Searches for neutrinos in the direction of radio-bright blazars with the ANTARES telescope
Authors:
ANTARES Collaboration,
A. Albert,
S. Alves,
M. André,
M. Ardid,
S. Ardid,
J. J. Aubert,
J Aublin,
B. Baret,
S. Basa,
Y. Becherini,
B. Belhorma,
M. Bendahman,
F. Benfenati,
V. Bertin,
S. Biagi,
M. Bissinger,
J. Boumaaza,
M. Bouta,
M. C. Bouwhuis,
H. Brânzaş,
R. Bruijn,
J. Brunner,
J. Busto,
B. Caiffi
, et al. (140 additional authors not shown)
Abstract:
Active galaxies, especially blazars, are among the most promising neutrino source candidates. To date, ANTARES searches for these objects considered GeV-TeV $γ$-ray bright blazars. Here, a statistically complete radio-bright blazar sample is used as the target for searches of origins of neutrinos collected by the ANTARES neutrino telescope over 13 years of operation. The hypothesis of a neutrino-b…
▽ More
Active galaxies, especially blazars, are among the most promising neutrino source candidates. To date, ANTARES searches for these objects considered GeV-TeV $γ$-ray bright blazars. Here, a statistically complete radio-bright blazar sample is used as the target for searches of origins of neutrinos collected by the ANTARES neutrino telescope over 13 years of operation. The hypothesis of a neutrino-blazar directional correlation is tested by pair counting and by a complementary likelihood-based approach. The resulting post-trial $p$-value is $3.0\%$ ($2.2σ$ in the two-sided convention), possibly indicating a correlation. Additionally, a time-dependent analysis is performed to search for temporal clustering of neutrino candidates as a mean of detecting neutrino flares in blazars. None of the investigated sources alone reaches a significant flare detection level. However, the presence of 18 sources with a pre-trial significance above $3σ$ indicates a $p=1.4\%$ ($2.5σ$ in the two-sided convention) detection of a time-variable neutrino flux. An \textit{a posteriori} investigation reveals an intriguing temporal coincidence of neutrino, radio, and $γ$-ray flares of the J0242+1101 blazar at a $p=0.5\%$ ($2.9σ$ in the two-sided convention) level. Altogether, the results presented here suggest a possible connection of neutrino candidates detected by the ANTARES telescope with radio-bright blazars.
△ Less
Submitted 13 September, 2023;
originally announced September 2023.
-
Multiple imaging of the quasar 2005+403 formed by anisotropic scattering
Authors:
T. A. Koryukova,
A. B. Pushkarev,
S. Kiehlmann,
A. C. S. Readhead
Abstract:
We report on the low Galactic latitude ($b=4.3^\circ$) quasar 2005$+$403, the second active galactic nuclei, in which we detected a rare phenomenon of multiple imaging induced by refractive-dominated scattering. The manifestation of this propagation effect is revealed at different frequencies ($\lesssim8$ GHz) and epochs of VLBA observations. The pattern formed by anisotropic scattering is stretch…
▽ More
We report on the low Galactic latitude ($b=4.3^\circ$) quasar 2005$+$403, the second active galactic nuclei, in which we detected a rare phenomenon of multiple imaging induced by refractive-dominated scattering. The manifestation of this propagation effect is revealed at different frequencies ($\lesssim8$ GHz) and epochs of VLBA observations. The pattern formed by anisotropic scattering is stretched out along the line of constant Galactic latitude with a local $\mathrm{PA}\approx40^\circ$ showing one-two sub-images, often on either side of the core. Analysing the multi-frequency VLBA data ranging from 1.4 to 43.2 GHz, we found that both the angular size of the apparent core component and the separation between the primary and secondary core images follow a wavelength squared dependence, providing convincing evidence for a plasma scattering origin for the multiple imaging. Based on the OVRO long-term monitoring data at 15 GHz obtained for 2005$+$403, we identified the characteristic flux density excursions occurred in April-May 2019 and attributed to an extreme scattering event (ESE) associated with the passage of a plasma lens across the line of sight. Modeling the ESE, we determined that the angular size of the screen is 0.4 mas and it drifts with the proper motion of 4.4 mas yr$^{-1}$. Assuming that the scattering screen is located in the highly turbulent Cygnus region, the transverse linear size and speed of the lens with respect to the observer are 0.7 AU and 37 km s$^{-1}$, respectively.
△ Less
Submitted 5 October, 2023; v1 submitted 29 August, 2023;
originally announced August 2023.
-
Milliarcsecond Core Size Dependence of the Radio Variability of Blazars
Authors:
Po Chih Hsu,
Jun Yi Koay,
Satoki Matsushita,
Chorng-Yuan Hwang,
Talvikki Hovatta,
Sebastian Kiehlmann,
Anthony Readhead,
Walter Max-Moerbeck,
Rodrigo Reeves
Abstract:
Studying the long-term radio variability (timescales of months to years) of blazars enables us to gain a better understanding of the physical structure of these objects on sub-parsec scales, and the physics of super massive black holes. In this study, we focus on the radio variability of 1157 blazars observed at 15 GHz through the Owens Valley Radio Observatory (OVRO) Blazar Monitoring Program. We…
▽ More
Studying the long-term radio variability (timescales of months to years) of blazars enables us to gain a better understanding of the physical structure of these objects on sub-parsec scales, and the physics of super massive black holes. In this study, we focus on the radio variability of 1157 blazars observed at 15 GHz through the Owens Valley Radio Observatory (OVRO) Blazar Monitoring Program. We investigate the dependence of the variability amplitudes and timescales, characterized based on model fitting to the structure functions, on the milliarcsecond core sizes measured by Very Long Baseline Interferometry. We find that the most compact sources at milliarcsecond scales exhibit larger variability amplitudes and shorter variability timescales than more extended sources. Additionally, for sources with measured redshifts and Doppler boosting factors, the correlation between linear core sizes against variability amplitudes and intrinsic timescales are also significant. The observed relationship between variability timescales and core sizes is expected, based on light travel-time arguments. This variability vs core size relation extends beyond the core sizes measured at 15 GHz; we see significant correlation between the 15 GHz variability amplitudes (as well as timescales) and core sizes measured at other frequencies, which can be attributed to a frequency-source size relationship arising from the intrinsic jet structure. At low frequencies of 1 GHz where the core sizes are dominated by interstellar scattering, we find that the variability amplitudes have significant correlation with the 1 GHz intrinsic core angular sizes, once the scatter broadening effects are deconvoluted from the intrinsic core sizes.
△ Less
Submitted 20 August, 2023; v1 submitted 18 August, 2023;
originally announced August 2023.
-
The RoboPol sample of optical polarimetric standards
Authors:
D. Blinov,
S. Maharana,
F. Bouzelou,
C. Casadio,
E. Gjerløw,
J. Jormanainen,
S. Kiehlmann,
J. A. Kypriotakis,
I. Liodakis,
N. Mandarakas,
L. Markopoulioti,
G. V. Panopoulou,
V. Pelgrims,
A. Pouliasi,
S. Romanopoulos,
R. Skalidis,
R. M. Anche,
E. Angelakis,
J. Antoniadis,
B. J. Medhi,
T. Hovatta,
A. Kus,
N. Kylafis,
A. Mahabal,
I. Myserlis
, et al. (12 additional authors not shown)
Abstract:
Optical polarimeters are typically calibrated using measurements of stars with known and stable polarization parameters. However, there is a lack of such stars available across the sky. Many of the currently available standards are not suitable for medium and large telescopes due to their high brightness. Moreover, as we find, some of the used polarimetric standards are in fact variable or have po…
▽ More
Optical polarimeters are typically calibrated using measurements of stars with known and stable polarization parameters. However, there is a lack of such stars available across the sky. Many of the currently available standards are not suitable for medium and large telescopes due to their high brightness. Moreover, as we find, some of the used polarimetric standards are in fact variable or have polarization parameters that differ from their cataloged values. Our goal is to establish a sample of stable standards suitable for calibrating linear optical polarimeters with an accuracy down to $10^{-3}$ in fractional polarization. For five years, we have been running a monitoring campaign of a sample of standard candidates comprised of 107 stars distributed across the northern sky. We analyzed the variability of the linear polarization of these stars, taking into account the non-Gaussian nature of fractional polarization measurements. For a subsample of nine stars, we also performed multiband polarization measurements. We created a new catalog of 65 stars (see Table 2) that are stable, have small uncertainties of measured polarimetric parameters, and can be used as calibrators of polarimeters at medium- and large-size telescopes.
△ Less
Submitted 12 July, 2023;
originally announced July 2023.
-
Bright-Moon Sky as a Wide-Field Linear Polarimetric Flat Source for Calibration
Authors:
S. Maharana,
S. Kiehlmann,
D. Blinov,
V. Pelgrims,
V. Pavlidou,
K. Tassis,
J. A. Kypriotakis,
A. N. Ramaprakash,
R. M. Anche,
A. Basyrov,
K. Deka,
H. K. Eriksen,
T. Ghosh,
E. Gjerløw,
N. Mandarakas,
E. Ntormousi,
G. V. Panopoulou,
A. Papadaki,
T. Pearson,
S. B. Potter,
A. C. S. Readhead,
R. Skalidis,
I. K. Wehus
Abstract:
Next-generation wide-field optical polarimeters like the Wide-Area Linear Optical Polarimeters (WALOPs) have a field of view (FoV) of tens of arcminutes. For efficient and accurate calibration of these instruments, wide-field polarimetric flat sources will be essential. Currently, no established wide-field polarimetric standard or flat sources exist. This paper tests the feasibility of using the p…
▽ More
Next-generation wide-field optical polarimeters like the Wide-Area Linear Optical Polarimeters (WALOPs) have a field of view (FoV) of tens of arcminutes. For efficient and accurate calibration of these instruments, wide-field polarimetric flat sources will be essential. Currently, no established wide-field polarimetric standard or flat sources exist. This paper tests the feasibility of using the polarized sky patches of the size of around ten-by-ten arcminutes, at a distance of up to 20 degrees from the Moon, on bright-Moon nights as a wide-field linear polarimetric flat source. We observed 19 patches of the sky adjacent to the bright-Moon with the RoboPol instrument in the SDSS-r broadband filter. These were observed on five nights within two days of the full-Moon across two RoboPol observing seasons. We find that for 18 of the 19 patches, the uniformity in the measured normalized Stokes parameters $q$ and $u$ is within 0.2 %, with 12 patches exhibiting uniformity within 0.07 % or better for both $q$ and $u$ simultaneously, making them reliable and stable wide-field linear polarization flats. We demonstrate that the sky on bright-Moon nights is an excellent wide-field linear polarization flat source. Various combinations of the normalized Stokes parameters $q$ and $u$ can be obtained by choosing suitable locations of the sky patch with respect to the Moon
△ Less
Submitted 7 May, 2023;
originally announced May 2023.
-
VERITAS discovery of very high energy gamma-ray emission from S3 1227+25 and multiwavelength observations
Authors:
Atreya Acharyya,
Colin Adams,
Avery Archer,
Priyadarshini Bangale,
Wystan Benbow,
Aryeh Brill,
Jodi Christiansen,
Alisha Chromey,
Manel Errando,
Abe Falcone,
Qi Feng,
John Finley,
Gregory Foote,
Lucy Fortson,
Amy Furniss,
Greg Gallagher,
William Hanlon,
David Hanna,
Olivier Hervet,
Claire Hinrichs,
John Hoang,
Jamie Holder,
Weidong Jin,
Madalyn Johnson,
Philip Kaaret
, et al. (46 additional authors not shown)
Abstract:
We report the detection of very high energy gamma-ray emission from the blazar S3 1227+25 (VER J1230+253) with the Very Energetic Radiation Imaging Telescope Array System (VERITAS). VERITAS observations of the source were triggered by the detection of a hard-spectrum GeV flare on May 15, 2015 with the Fermi-Large Area Telescope (LAT). A combined five-hour VERITAS exposure on May 16th and May 18th…
▽ More
We report the detection of very high energy gamma-ray emission from the blazar S3 1227+25 (VER J1230+253) with the Very Energetic Radiation Imaging Telescope Array System (VERITAS). VERITAS observations of the source were triggered by the detection of a hard-spectrum GeV flare on May 15, 2015 with the Fermi-Large Area Telescope (LAT). A combined five-hour VERITAS exposure on May 16th and May 18th resulted in a strong 13$σ$ detection with a differential photon spectral index, $Γ$ = 3.8 $\pm$ 0.4, and a flux level at 9% of the Crab Nebula above 120 GeV. This also triggered target of opportunity observations with Swift, optical photometry, polarimetry and radio measurements, also presented in this work, in addition to the VERITAS and Fermi-LAT data. A temporal analysis of the gamma-ray flux during this period finds evidence of a shortest variability timescale of $τ_{obs}$ = 6.2 $\pm$ 0.9 hours, indicating emission from compact regions within the jet, and the combined gamma-ray spectrum shows no strong evidence of a spectral cut-off. An investigation into correlations between the multiwavelength observations found evidence of optical and gamma-ray correlations, suggesting a single-zone model of emission. Finally, the multiwavelength spectral energy distribution is well described by a simple one-zone leptonic synchrotron self-Compton radiation model.
△ Less
Submitted 4 May, 2023;
originally announced May 2023.
-
COMAP Early Science: VIII. A Joint Stacking Analysis with eBOSS Quasars
Authors:
Delaney A. Dunne,
Kieran A. Cleary,
Patrick C. Breysse,
Dongwoo T. Chung,
Havard T. Ihle,
J. Richard Bond,
Hans Kristian Eriksen,
Joshua Ott Gundersen,
Laura C. Keating,
Junhan Kim,
Jonas Gahr Sturtzel Lunde,
Norman Murray,
Hamsa Padmanabhan,
Liju Philip,
Nils-Ole Stutzer,
Doga Tolgay,
Ingunn Katherine Wehus,
Sarah E. Church,
Todd Gaier,
Andrew I. Harris,
Richard Hobbs,
James W. Lamb,
Charles R. Lawrence,
Anthony C. S. Readhead,
David P. Woody
Abstract:
We present a new upper limit on the cosmic molecular gas density at $z=2.4-3.4$ obtained using the first year of observations from the CO Mapping Array Project (COMAP). COMAP data cubes are stacked on the 3D positions of 243 quasars selected from the Extended Baryon Oscillation Spectroscopic Survey (eBOSS) catalog, yielding a 95% upper limit for flux from CO(1-0) line emission of 0.129 Jy km/s. De…
▽ More
We present a new upper limit on the cosmic molecular gas density at $z=2.4-3.4$ obtained using the first year of observations from the CO Mapping Array Project (COMAP). COMAP data cubes are stacked on the 3D positions of 243 quasars selected from the Extended Baryon Oscillation Spectroscopic Survey (eBOSS) catalog, yielding a 95% upper limit for flux from CO(1-0) line emission of 0.129 Jy km/s. Depending on the balance of the emission between the quasar host and its environment, this value can be interpreted as an average CO line luminosity $L'_\mathrm{CO}$ of eBOSS quasars of $\leq 1.26\times10^{11}$ K km pc$^2$ s$^{-1}$, or an average molecular gas density $ρ_\mathrm{H_2}$ in regions of the universe containing a quasar of $\leq 1.52\times10^8$ M$_\odot$ cMpc$^{-3}$. The $L'_\mathrm{CO}$ upper limit falls among CO line luminosities obtained from individually-targeted quasars in the COMAP redshift range, and the $ρ_\mathrm{H_2}$ value is comparable to upper limits obtained from other Line Intensity Mapping (LIM) surveys and their joint analyses. Further, we forecast the values obtainable with the COMAP/eBOSS stack after the full 5-year COMAP Pathfinder survey. We predict that a detection is probable with this method, depending on the CO properties of the quasar sample. Based on the achieved sensitivity, we believe that this technique of stacking LIM data on the positions of traditional galaxy or quasar catalogs is extremely promising, both as a technique for investigating large galaxy catalogs efficiently at high redshift and as a technique for bolstering the sensitivity of LIM experiments, even with a fraction of their total expected survey data.
△ Less
Submitted 26 February, 2024; v1 submitted 19 April, 2023;
originally announced April 2023.
-
The C-Band All-Sky Survey (C-BASS): New Constraints on the Integrated Radio Spectrum of M 31
Authors:
Stuart E. Harper,
Adam Barr,
C. Dickinson,
M. W. Peel,
Roke Cepeda-Arroita,
C. J. Copley,
R. D. P. Grumitt,
J. Patrick Leahy,
J. L. Jonas,
Michael E. Jones,
J. Leech,
T. J. Pearson,
A. C. S. Readhead,
Angela C. Taylor
Abstract:
The Andromeda galaxy (M31) is our closest neighbouring spiral galaxy, making it an ideal target for studying the physics of the interstellar medium in a galaxy very similar to our own. Using new observations of M31 at 4.76GHz by the C-Band All-Sky Survey (C-BASS), and all available radio data at $1^\circ$ resolution, we produce the integrated spectrum and put new constraints on the synchrotron spe…
▽ More
The Andromeda galaxy (M31) is our closest neighbouring spiral galaxy, making it an ideal target for studying the physics of the interstellar medium in a galaxy very similar to our own. Using new observations of M31 at 4.76GHz by the C-Band All-Sky Survey (C-BASS), and all available radio data at $1^\circ$ resolution, we produce the integrated spectrum and put new constraints on the synchrotron spectral index and anomalous microwave emission (AME) from M31. We use aperture photometry and spectral modelling to fit for the integrated spectrum of M31, and subtract a comprehensive model of nearby background radio sources. The AME in M31 is detected at $3σ$ significance with a peak near 30GHz and flux density $0.27\pm0.09$Jy. The synchrotron spectral index of M31 is flatter than our own Galaxy at $α= -0.66 \pm 0.03$ with no strong evidence of spectral curvature. The emissivity of AME, averaged over the total emission from M31 is lower than typical AME sources in our Galaxy, implying that AME is not uniformly distributed throughout M31 and instead is likely confined to sub-regions -- this will need to be confirmed using future higher resolution observations around 20--30GHz.
△ Less
Submitted 18 May, 2023; v1 submitted 7 April, 2023;
originally announced April 2023.
-
Compact Symmetric Objects -- III Evolution of the High-Luminosity Branch and a Possible Connection with Tidal Disruption Events
Authors:
A. C. S. Readhead,
V. Ravi,
R. D. Blandford,
A. G. Sullivan,
J. Somalwar,
M. C. Begelman,
M. Birkinshaw,
I. Liodakis,
M. L. Lister,
T. J. Pearson,
G. B. Taylor,
P. N. Wilkinson,
N. Globus,
S. Kiehlmann,
C. R. Lawrence,
D. Murphy,
S. O'Neill,
V. Pavlidou,
E. Sheldahl,
A. Siemiginowska,
K. Tassis
Abstract:
We use a sample of 54 Compact Symmetric Objects (CSOs) to confirm that there are two unrelated CSO classes: an edge-dimmed, low-luminosity class (CSO~1), and an edge-brightened, high-luminosity class (CSO~2). Using blind tests, we show that CSO~2s consist of three sub-classes: CSO 2.0, having prominent hot-spots at the leading edges of narrow jets and/or narrow lobes; CSO~2.2, without prominent ho…
▽ More
We use a sample of 54 Compact Symmetric Objects (CSOs) to confirm that there are two unrelated CSO classes: an edge-dimmed, low-luminosity class (CSO~1), and an edge-brightened, high-luminosity class (CSO~2). Using blind tests, we show that CSO~2s consist of three sub-classes: CSO 2.0, having prominent hot-spots at the leading edges of narrow jets and/or narrow lobes; CSO~2.2, without prominent hot-spots, and with broad jets and/or lobes; and CSO~2.1, which exhibit mixed properties. Most CSO 2s do not evolve into larger jetted-AGN, but spend their whole life-cycle as CSOs of size $\lesssim$500 pc and age $\lesssim$5000 yr. The minimum energies needed to produce the radio luminosity and structure in CSO~2s range from $\sim~10^{-4}\,M_\odot{c}^2$ to $\sim7\,M_\odot{c}^2$. We show that the transient nature of most CSO~2s, and their birthrate, can be explained through ignition in the tidal disruption events of giant stars. We also consider possibilities of tapping the spin energy of the supermassive black hole, and tapping the energy of the accretion disk. Our results demonstrate that CSOs constitute a large family of AGN in which we have thus far studied only the brightest. More comprehensive CSO studies, with higher sensitivity, resolution, and dynamic range, will revolutionize our understanding of AGN and the central engines that power them.
△ Less
Submitted 26 November, 2023; v1 submitted 20 March, 2023;
originally announced March 2023.
-
Compact Symmetric Objects -- II Confirmation of a Distinct Population of High-Luminosity Jetted Active Galaxies
Authors:
S. Kiehlmann,
A. C. S. Readhead,
S. O'Neill,
P. N. Wilkinson,
M. L. Lister,
I. Liodakis,
S. Bruzewski,
V. Pavlidou,
T. J. Pearson,
E. Sheldahl,
A. Siemiginowska,
K. Tassis,
G. B. Taylor
Abstract:
Compact Symmetric Objects (CSOs) are compact (<1 kpc), jetted Active Galactic Nuclei (AGN), whose jet axes are not aligned close to the line of sight, and whose observed emission is not predominantly relativistically boosted towards us. Two classes of CSOs have previously been identified: approximately one fifth are edge-dimmed and designated as CSO 1s, while the rest are edge brightened and desig…
▽ More
Compact Symmetric Objects (CSOs) are compact (<1 kpc), jetted Active Galactic Nuclei (AGN), whose jet axes are not aligned close to the line of sight, and whose observed emission is not predominantly relativistically boosted towards us. Two classes of CSOs have previously been identified: approximately one fifth are edge-dimmed and designated as CSO 1s, while the rest are edge brightened and designated as CSO 2s. This paper focuses almost exclusively on CSO 2s. Using complete samples of CSO 2s we present three independent lines of evidence, based on their relative numbers, redshift distributions, and size distributions, which show conclusively that the vast majority (> 99%) of CSO 2s do not evolve into larger-scale radio sources. These CSO 2s belong to a distinct population of jetted-AGN, which should be characterized as ``short-lived'' compared to the classes of larger jetted-AGN, as opposed to ``young''. We show that there is a sharp upper cutoff in the CSO 2 size distribution at $\approx 500$ pc. The distinct differences between most CSO 2s and other jetted-AGN provides a crucial new time domain window on the formation and evolution of relativistic jets in AGN and the supermassive black holes that drive them.
△ Less
Submitted 26 November, 2023; v1 submitted 20 March, 2023;
originally announced March 2023.
-
Compact Symmetric Objects -- I Towards a Comprehensive Bona Fide Catalog
Authors:
S. Kiehlmann,
M. L. Lister,
A. C. S. Readhead,
I. Liodakis,
S. O'Neill,
T. J. Pearson,
E. Sheldahl,
A. Siemiginowska,
K. Tassis,
G. B. Taylor,
P. N. Wilkinson
Abstract:
Compact Symmetric Objects (CSOs) are jetted Active Galactic Nuclei (AGN) with overall projected size <1 kpc. The classification was introduced to distinguish these objects from the majority of compact jetted-AGN in centimeter wavelength very long baseline interferometry observations, where the observed emission is relativistically boosted towards the observer. The original classification criteria…
▽ More
Compact Symmetric Objects (CSOs) are jetted Active Galactic Nuclei (AGN) with overall projected size <1 kpc. The classification was introduced to distinguish these objects from the majority of compact jetted-AGN in centimeter wavelength very long baseline interferometry observations, where the observed emission is relativistically boosted towards the observer. The original classification criteria for CSOs were: (i) evidence of emission on both sides of the center of activity, and (ii) overall size <1 kpc. However some relativistically boosted objects with jet axes close to the line of sight appear symmetric and have been mis-classified as CSOs, thereby undermining the CSO classification. This is because two essential CSO properties, pointed out in the original papers, have been neglected: (iii) low variability, and (iv) low apparent speeds along the jets. As a first step towards creating a comprehensive catalog of ``bona fide'' CSOs, we identify 79 bona fide CSOs, including 15 objects claimed as confirmed CSOs here for the first time, that match the CSO selection criteria. This sample of bona fide CSOs can be used for astrophysical studies of CSOs without contamination by mis-classified CSOs. We show that the fraction of CSOs in complete flux density limited AGN samples with S$_{\rm 5\,GHz}$ >700 mJy is between $(6.8\pm1.6)$% and $(8.5\pm1.8)$%.
△ Less
Submitted 26 November, 2023; v1 submitted 20 March, 2023;
originally announced March 2023.
-
VLBI Scrutiny of a New Neutrino-Blazar Multiwavelength-Flare Coincidence
Authors:
F. Eppel,
M. Kadler,
E. Ros,
F. Roesch,
J. Hessdoerfer,
P. Benke,
P. G. Edwards,
C. M. Fromm,
M. Giroletti,
A. Gokus,
J. L. Gomez,
S. Haemmerich,
D. Kirchner,
Y. Y. Kovalev,
T. P. Krichbaum,
M. L. Lister,
C. Nanci,
R. Ojha,
G. F. Paraschos,
A. Plavin,
A. C. S. Readhead,
J. Stevens,
P. Weber
Abstract:
In the past years, evidence has started piling up that some high-energy cosmic neutrinos can be associated with blazars in flaring states. On February 26, 2022, a new blazar-neutrino coincidence has been reported: the track-like neutrino event IC220225A detected by IceCube is spatially coincident with the flat-spectrum radio quasar PKS 0215+015. Like previous associations, this source was found to…
▽ More
In the past years, evidence has started piling up that some high-energy cosmic neutrinos can be associated with blazars in flaring states. On February 26, 2022, a new blazar-neutrino coincidence has been reported: the track-like neutrino event IC220225A detected by IceCube is spatially coincident with the flat-spectrum radio quasar PKS 0215+015. Like previous associations, this source was found to be in a high optical and $γ$-ray state. Moreover, the source showed a bright radio outburst, which substantially increases the probability of a true physical association. We have performed six observations with the VLBA shortly after the neutrino event with a monthly cadence and are monitoring the source with the Effelsberg 100m-Telescope, and with the Australia Compact Telescope Array. Here, we present first results on the contemporary parsec-scale jet structure of PKS 0215+015 in total intensity and polarization to constrain possible physical processes leading to neutrino emission in blazars.
△ Less
Submitted 24 February, 2023; v1 submitted 31 January, 2023;
originally announced January 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.
-
A multiband study and exploration of the radio wave - $γ$-ray connection in 3C 84
Authors:
G. F. Paraschos,
V. Mpisketzis,
J. -Y. Kim,
G. Witzel,
T. P. Krichbaum,
J. A. Zensus,
M. A. Gurwell,
A. Lähteenmäki,
M. Tornikoski,
S. Kiehlmann,
A. C. S. Readhead
Abstract:
Total intensity variability light curves offer a unique insight into the ongoing debate about the launching mechanism of jets. For this work, we utilise the availability of radio and $γ$-ray light curves over a few decades of the radio source 3C 84 (NGC 1275). We calculate the multiband time lags between the flares identified in the light curves via discrete cross-correlation and Gaussian process…
▽ More
Total intensity variability light curves offer a unique insight into the ongoing debate about the launching mechanism of jets. For this work, we utilise the availability of radio and $γ$-ray light curves over a few decades of the radio source 3C 84 (NGC 1275). We calculate the multiband time lags between the flares identified in the light curves via discrete cross-correlation and Gaussian process regression. We find that the jet particle and magnetic field energy densities are in equipartition ($k_\textrm{r} = 1.08\pm0.18$). The jet apex is located $z_\textrm{91.5 GHz}= 22 - 645$ $R_\textrm{s}$ ($2 - 20 \times 10^{-3}$ pc) upstream of the 3 mm radio core; at that position, the magnetic field amplitude is $B_\textrm{core}^\textrm{91.5 GHz}= 3 - 10$ G. Our results are in good agreement with earlier studies, which utilised very-long-baseline interferometry. Furthermore, we investigate the temporal relation between the ejection of radio and $γ$-ray flares. Our results are in favour of the $γ$-ray emission being associated with the radio emission. We are able to tentatively connect the ejection of features identified at 43 and 86 GHz to prominent $γ$-ray flares. Finally, we compute the multiplicity parameter $λ$ and the Michel magnetisation $σ_\textrm{M}$ and find that they are consistent with a jet launched by the Blandford & Znajek 1977 mechanism.
△ Less
Submitted 18 October, 2022;
originally announced October 2022.
-
Multi-messenger characterization of Mrk 501 during historically low X-ray and $γ$-ray activity
Authors:
MAGIC collaboration,
H. Abe,
S. Abe,
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,
J. A. Barrio,
I. Batković,
J. Baxter,
J. Becerra González,
W. Bednarek,
E. Bernardini,
M. Bernardos,
A. Berti,
J. Besenrieder
, et al. (300 additional authors not shown)
Abstract:
We study the broadband emission of Mrk 501 using multi-wavelength observations from 2017 to 2020 performed with a multitude of instruments, involving, among others, MAGIC, Fermi-LAT, NuSTAR, Swift, GASP-WEBT, and OVRO. Mrk 501 showed an extremely low broadband activity, which may help to unravel its baseline emission. Nonetheless, significant flux variations are detected at all wavebands, with the…
▽ More
We study the broadband emission of Mrk 501 using multi-wavelength observations from 2017 to 2020 performed with a multitude of instruments, involving, among others, MAGIC, Fermi-LAT, NuSTAR, Swift, GASP-WEBT, and OVRO. Mrk 501 showed an extremely low broadband activity, which may help to unravel its baseline emission. Nonetheless, significant flux variations are detected at all wavebands, with the highest occurring at X-rays and very-high-energy (VHE) $γ$-rays. A significant correlation ($>$3$σ$) between X-rays and VHE $γ$-rays is measured, supporting leptonic scenarios to explain the variable parts of the emission, also during low activity. This is further supported when we extend our data from 2008 to 2020, and identify, for the first time, significant correlations between Swift-XRT and Fermi-LAT. We additionally find correlations between high-energy $γ$-rays and radio, with the radio lagging by more than 100 days, placing the $γ$-ray emission zone upstream of the radio-bright regions in the jet. Furthermore, Mrk 501 showed a historically low activity in X-rays and VHE $γ$-rays from mid-2017 to mid-2019 with a stable VHE flux ($>$0.2 TeV) of 5% the emission of the Crab Nebula. The broadband spectral energy distribution (SED) of this 2-year-long low-state, the potential baseline emission of Mrk 501, can be characterized with one-zone leptonic models, and with (lepto)-hadronic models fulfilling neutrino flux constraints from IceCube. We explore the time evolution of the SED towards the low-state, revealing that the stable baseline emission may be ascribed to a standing shock, and the variable emission to an additional expanding or traveling shock.
△ Less
Submitted 5 March, 2023; v1 submitted 5 October, 2022;
originally announced October 2022.
-
Experimental Investigation of Drain Noise in High Electron Mobility Transistors: Thermal and Hot Electron Noise
Authors:
Bekari Gabritchidze,
Kieran A. Cleary,
Anthony C. Readhead,
Austin J. Minnich
Abstract:
We report the on-wafer characterization of $S$-parameters and microwave noise temperature ($T_{50}$) of discrete metamorphic InGaAs high electron mobility transistors (mHEMTs) at 40 K and 300 K and over a range of drain-source voltages ($V_{DS}$). From these data, we extract a small-signal model and the drain (output) noise current power spectral density ($S_{id}$) at each bias and temperature. Th…
▽ More
We report the on-wafer characterization of $S$-parameters and microwave noise temperature ($T_{50}$) of discrete metamorphic InGaAs high electron mobility transistors (mHEMTs) at 40 K and 300 K and over a range of drain-source voltages ($V_{DS}$). From these data, we extract a small-signal model and the drain (output) noise current power spectral density ($S_{id}$) at each bias and temperature. This procedure enables $S_{id}$ to be obtained while accounting for the variation of small-signal model, noise impedance match, and other parameters under the various conditions. We find that the thermal noise associated with the channel conductance can only account for a portion of the measured output noise. Considering the variation of output noise with physical temperature and bias and prior studies of microwave noise in quantum wells, we hypothesize that a hot electron noise source based on real-space transfer of electrons from the channel to the barrier could account for the remaining portion of $S_{id}$. We suggest further studies to gain insights into the physical mechanisms. Finally, we calculate that the minimum HEMT noise temperature could be reduced by up to $\sim 50$% and $\sim 30$% at cryogenic temperature and room temperature, respectively, if the hot electron noise could be suppressed.
△ Less
Submitted 20 March, 2024; v1 submitted 6 September, 2022;
originally announced September 2022.
-
WALOP-South: A Four-Camera One-Shot Imaging Polarimeter for PASIPHAE Survey. Paper II -- Polarimetric Modelling and Calibration
Authors:
Siddharth Maharana,
Ramya M. Anche,
A. N. Ramaprakash,
Bhushan Joshi,
Artem Basyrov,
Dmitry Blinov,
Carolina Casadio,
Kishan Deka,
Hans Kristian Eriksen,
Tuhin Ghosh,
Eirik Gjerløw,
John A. Kypriotakis,
Sebastian Kiehlmann,
Nikolaos Mandarakas,
Georgia V. Panopoulou,
Katerina Papadaki,
Vasiliki Pavlidou,
Timothy J. Pearson,
Vincent Pelgrims,
Stephen B. Potter,
Anthony C. S. Readhead,
Raphael Skalidis,
Trygve Leithe Svalheim,
Konstantinos Tassis,
Ingunn K. Wehus
Abstract:
The Wide-Area Linear Optical Polarimeter (WALOP)-South instrument is an upcoming wide-field and high-accuracy optical polarimeter to be used as a survey instrument for carrying out the Polar-Areas Stellar Imaging in Polarization High Accuracy Experiment (PASIPHAE) program. Designed to operate as a one-shot four-channel and four-camera imaging polarimeter, it will have a field of view of…
▽ More
The Wide-Area Linear Optical Polarimeter (WALOP)-South instrument is an upcoming wide-field and high-accuracy optical polarimeter to be used as a survey instrument for carrying out the Polar-Areas Stellar Imaging in Polarization High Accuracy Experiment (PASIPHAE) program. Designed to operate as a one-shot four-channel and four-camera imaging polarimeter, it will have a field of view of $35\times 35$ arcminutes and will measure the Stokes parameters $I$, $q$, and $u$ in a single exposure in the SDSS-r broadband filter. The design goal for the instrument is to achieve an overall polarimetric measurement accuracy of 0.1 % over the entire field of view. We present here the complete polarimetric modeling of the instrument, characterizing the amount and sources of instrumental polarization. To accurately retrieve the real Stokes parameters of a source from the measured values, we have developed a calibration method for the instrument. Using this calibration method and simulated data, we demonstrate how to correct instrumental polarization and obtain 0.1 % accuracy in the degree of polarization, $p$. Additionally, we tested and validated the calibration method by implementing it on a table-top WALOP-like test-bed polarimeter in the laboratory.
△ Less
Submitted 26 August, 2022;
originally announced August 2022.
-
The hunt for extraterrestrial high-energy neutrino counterparts
Authors:
I. Liodakis,
T. Hovatta,
V. Pavlidou,
A. C. S. Readhead,
R. D. Blandford,
S. Kiehlmann,
E. Lindfors,
W. Max-Moerbeck,
T. J. Pearson,
M. Petropoulou
Abstract:
The origin of Petaelectronvolt (PeV) astrophysical neutrinos is fundamental to our understanding of the high-energy Universe. Apart from the technical challenges of operating detectors deep below ice, oceans, and lakes, the phenomenological challenges are even greater than those of gravitational waves; the sources are unknown, hard to predict, and we lack clear signatures. Neutrino astronomy there…
▽ More
The origin of Petaelectronvolt (PeV) astrophysical neutrinos is fundamental to our understanding of the high-energy Universe. Apart from the technical challenges of operating detectors deep below ice, oceans, and lakes, the phenomenological challenges are even greater than those of gravitational waves; the sources are unknown, hard to predict, and we lack clear signatures. Neutrino astronomy therefore represents the greatest challenge faced by the astronomy and physics communities thus far. The possible neutrino sources range from accretion disks and tidal disruption events, to relativistic jets and galaxy clusters with blazar TXS~0506+056 the most compelling association thus far. Since that association, immense effort has been put into proving or disproving that jets are indeed neutrino emitters, but to no avail. By generating simulated neutrino counterpart samples, we explore the potential of detecting a significant correlation of neutrinos with jets from active galactic nuclei. We find that, given the existing challenges, even our best experiments could not have produced a $>3σ$ result. Larger programs over the next few years will be able to detect a significant correlation only if the brightest radio sources, rather than all jetted active galactic nuclei, are neutrino emitters. We discuss the necessary strategies required to steer future efforts into successful experiments.
△ Less
Submitted 15 August, 2022;
originally announced August 2022.
-
Starlight-polarization-based tomography of the magnetized interstellar medium: PASIPHAE's line-of-sight inversion method
Authors:
V. Pelgrims,
G. V. Panopoulou,
K. Tassis,
V. Pavlidou,
A. Basyrov,
D. Blinov,
E. Gjerløw,
S. Kiehlmann,
N. Mandarakas,
A. Papadaki,
R. Skalidis,
A. Tsouros,
R. M. Anche,
H. K. Eriksen,
T. Ghosh,
J. A. Kypriotakis,
S. Maharana,
E. Ntormousi,
T. J. Pearson,
S. B. Potter,
A. N. Ramaprakash,
A. C. S. Readhead,
I. K. Wehus
Abstract:
We present the first Bayesian method for tomographic decomposition of the plane-of-sky orientation of the magnetic field with the use of stellar polarimetry and distance. This standalone tomographic inversion method presents an important step forward in reconstructing the magnetized interstellar medium (ISM) in 3D within dusty regions. We develop a model in which the polarization signal from the m…
▽ More
We present the first Bayesian method for tomographic decomposition of the plane-of-sky orientation of the magnetic field with the use of stellar polarimetry and distance. This standalone tomographic inversion method presents an important step forward in reconstructing the magnetized interstellar medium (ISM) in 3D within dusty regions. We develop a model in which the polarization signal from the magnetized and dusty ISM is described by thin layers at various distances. Our modeling makes it possible to infer the mean polarization (amplitude and orientation) induced by individual dusty clouds and to account for the turbulence-induced scatter in a generic way. We present a likelihood function that explicitly accounts for uncertainties in polarization and parallax. We develop a framework for reconstructing the magnetized ISM through the maximization of the log-likelihood using a nested sampling method. We test our Bayesian inversion method on mock data taking into account realistic uncertainties from Gaia and as expected for the optical polarization survey PASIPHAE according to the currently planned observing strategy. We demonstrate that our method is effective at recovering the cloud properties as soon as the polarization induced by a cloud to its background stars is higher than $\sim 0.1\%$ for the adopted survey exposure time and level of systematic uncertainty. Our method makes it possible to recover not only the mean polarization properties but also to characterize the intrinsic scatter, thus creating new ways to characterize ISM turbulence and the magnetic field strength. Finally, we apply our method to an existing data set of starlight polarization with known line-of-sight decomposition, demonstrating agreement with previous results and an improved quantification of uncertainties in cloud properties.
△ Less
Submitted 28 February, 2023; v1 submitted 3 August, 2022;
originally announced August 2022.
-
Jet kinematics in the transversely stratified jet of 3C 84 A two-decade overview
Authors:
G. F. Paraschos,
T. P. Krichbaum,
J. -Y. Kim,
J. A. Hodgson,
J. Oh,
E. Ros,
J. A. Zensus,
A. P. Marscher,
S. G. Jorstad,
M. A. Gurwell,
A. Lähteenmäki,
M. Tornikoski,
S. Kiehlmann,
A. C. S. Readhead
Abstract:
3C84 (NGC1275) is one of the brightest radio sources in the mm radio-bands, which led to a plethora of VLBI observations at numerous frequencies over the years. They reveal a two-sided jet structure, with an expanding but not well-collimated parsec-scale jet, pointing southward. High resolution mm-VLBI observations allow the study and imaging of the jet base on the sub-parsec scale. This could fac…
▽ More
3C84 (NGC1275) is one of the brightest radio sources in the mm radio-bands, which led to a plethora of VLBI observations at numerous frequencies over the years. They reveal a two-sided jet structure, with an expanding but not well-collimated parsec-scale jet, pointing southward. High resolution mm-VLBI observations allow the study and imaging of the jet base on the sub-parsec scale. This could facilitate the investigation of the nature of the jet origin, also in view of the previously detected two-railed jet structure and east-west oriented core region seen with RadioAstron at 22 GHz. We produce VLBI images of this core and inner jet region, observed during the past twenty years at 15, 43, and 86 GHz. We determine the kinematics of the inner jet and ejected features at 43 and 86 GHz and compare their ejection times with radio and $γ$-ray variability. For the moving jet features, we find an average velocity of $β^\textrm{avg}_\textrm{app} = 0.055-0.22$c ($μ^\textrm{avg} = 0.04-0.18\,$mas/yr). From the time-averaged VLBI images at the three frequencies, we measure the transverse jet width along the bulk flow. On the $\leq 1.5$ parsec scale, we find a clear trend of the jet width being frequency dependent, with the jet being narrower at higher frequencies. This stratification is discussed in the context of a spine-sheath scenario, and is compared to other possible interpretations. From quasi-simultaneous observations at 43 and 86\,GHz we obtain spectral index maps, revealing a time-variable orientation of the spectral index gradient, due to structural variability of the inner jet.
△ Less
Submitted 3 June, 2022; v1 submitted 20 May, 2022;
originally announced May 2022.
-
Self-heating of cryogenic high-electron-mobility transistor amplifiers and the limits of microwave noise performance
Authors:
Anthony J. Ardizzi,
Alexander Y. Choi,
Bekari Gabritchidze,
Jacob Kooi,
Kieran A. Cleary,
Anthony C. Readhead,
Austin J. Minnich
Abstract:
The fundamental limits of the microwave noise performance of high electron mobility transistors (HEMTs) are of scientific and practical interest for applications in radio astronomy and quantum computing. Self-heating at cryogenic temperatures has been reported to be a limiting mechanism for the noise, but cryogenic cooling strategies to mitigate it, for instance using liquid cryogens, have not bee…
▽ More
The fundamental limits of the microwave noise performance of high electron mobility transistors (HEMTs) are of scientific and practical interest for applications in radio astronomy and quantum computing. Self-heating at cryogenic temperatures has been reported to be a limiting mechanism for the noise, but cryogenic cooling strategies to mitigate it, for instance using liquid cryogens, have not been evaluated. Here, we report microwave noise measurements of a packaged two-stage HEMT amplifier immersed in normal and superfluid $^4$He baths and in vacuum from 1.6 - 80 K. We find that these liquid cryogens are unable to mitigate the thermal noise associated with self-heating. Considering this finding, we examine the implications for the lower bounds of cryogenic noise performance in HEMTs. Our analysis supports the general design principle for cryogenic HEMTs of maximizing gain at the lowest possible power.
△ Less
Submitted 4 August, 2022; v1 submitted 8 May, 2022;
originally announced May 2022.
-
The C-Band All-Sky Survey (C-BASS): Template Fitting of Diffuse Galactic Microwave Emission in the Northern Sky
Authors:
S. E. Harper,
C. Dickinson,
A. Barr,
R. Cepeda-Arroita,
R. D. P. Grumitt,
H. M. Heilgendorff,
L. Jew,
J. L. Jonas,
M. E. Jones,
J. P. Leahy,
J. Leech,
T. J. Pearson,
M. W. Peel,
A. C. S. Readhead,
A. C. Taylor
Abstract:
The C-Band All-Sky Survey (C-BASS) has observed the Galaxy at 4.76GHz with an angular resolution of $0.73^\circ$ full-width half-maximum, and detected Galactic synchrotron emission with high signal-to-noise ratio over the entire northern sky ($δ> -15^{\circ}$). We present the results of a spatial correlation analysis of Galactic foregrounds at mid-to-high ($b > 10^\circ$) Galactic latitudes using…
▽ More
The C-Band All-Sky Survey (C-BASS) has observed the Galaxy at 4.76GHz with an angular resolution of $0.73^\circ$ full-width half-maximum, and detected Galactic synchrotron emission with high signal-to-noise ratio over the entire northern sky ($δ> -15^{\circ}$). We present the results of a spatial correlation analysis of Galactic foregrounds at mid-to-high ($b > 10^\circ$) Galactic latitudes using a preliminary version of the C-BASS intensity map. We jointly fit for synchrotron, dust, and free-free components between $20$ and $1000$GHz and look for differences in the Galactic synchrotron spectrum, and the emissivity of anomalous microwave emission (AME) when using either the C-BASS map or the 408MHz all-sky map to trace synchrotron emission. We find marginal evidence for a steepening ($\left<Δβ\right> = -0.06\pm0.02$) of the Galactic synchrotron spectrum at high frequencies resulting in a mean spectral index of $\left<β\right> = -3.10\pm0.02$ over $4.76-22.8$GHz. Further, we find that the synchrotron emission can be well modelled by a single power-law up to a few tens of GHz. Due to this, we find that the AME emissivity is not sensitive to changing the synchrotron tracer from the 408MHz map to the 4.76GHz map. We interpret this as strong evidence for the origin of AME being spinning dust emission.
△ Less
Submitted 8 July, 2022; v1 submitted 21 February, 2022;
originally announced February 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.
-
Locating the $γ$-ray emitting region in the quasar 4C +01.28
Authors:
F. Rösch,
M. Kadler,
E. Ros,
M. Gurwell,
T. Hovatta,
M. Kreter,
N. R. MacDonald,
A. C. S. Readhead
Abstract:
Determining the location of $γ$-ray emission in blazar jets is a challenging task. Pinpointing the exact location of $γ$-ray production within a relativistic jet can place strong constraints on our understanding of high-energy astrophysics and astroparticle physics. We present a study of the radio- and $γ$-bright flat-spectrum radio quasar (FSRQ) 4C +01.28 (PKS B1055+018) in which we try to pinpoi…
▽ More
Determining the location of $γ$-ray emission in blazar jets is a challenging task. Pinpointing the exact location of $γ$-ray production within a relativistic jet can place strong constraints on our understanding of high-energy astrophysics and astroparticle physics. We present a study of the radio- and $γ$-bright flat-spectrum radio quasar (FSRQ) 4C +01.28 (PKS B1055+018) in which we try to pinpoint the emission site of several prominent GeV flares. This source shows prominent high-amplitude broadband variability on time scales ranging from days to years. We combine high-resolution VLBI observations and multi-band radio light curves over a period of around nine years. We can associate two bright and compact newly ejected jet components with bright flares observed by the Fermi/LAT $γ$-ray telescope and at various radio frequencies. A cross-correlation analysis reveals the radio light curves systematically lag behind the $γ$-rays. In combination with the jet kinematics as measured by the VLBA, we use these cross-correlations to constrain a model in which the flares become observable at a given frequency when a plasma component passes through the region at which the bulk energy dissipation takes place at that frequency. We derive a lower limit of the location of the $γ$-ray emitting region in 4C +01.28 of several parsecs from the jet base, well beyond the expected extent of the broad-line region. This observational limit challenges blazar-emission models that rely on the broad-line region as a source of seed photons for inverse-Compton scattering.
△ Less
Submitted 16 December, 2021;
originally announced December 2021.
-
Multiwavelength variability power spectrum analysis of the blazars 3C 279 and PKS 1510-089 on multiple timescales
Authors:
Arti Goyal,
Marian Soida,
Lukasz Stawarz,
Paul J. Wiita,
Kari Nilsson,
Svetlana Jorstad,
Alan P. Marscher,
Margo F. Aller,
Hugh D. Aller,
Anne Lahteenmaki,
Talvikki Hovatta,
Staszek Zola,
Krzysztof Nalewajko,
Merja Tornikoski,
Joni Tammi,
Mark Hodges,
Sebastian Kiehlmann,
Anthony C. S. Readhead,
Walter Max-Moerbeck,
Elina Lindfors,
Vandad Fallah Ramazani,
D. E. Reichart,
D. B. Caton,
Janeth Valverde,
Deirdre Horan
, et al. (2 additional authors not shown)
Abstract:
We present the results of variability power spectral density (PSD) analysis using multiwavelength radio to GeV\,$γ$-ray light curves covering decades/years to days/minutes timescales for the blazars 3C 279 and PKS 1510-089. The PSDs are modeled as single power-laws, and the best-fit spectral shape is derived using the `power spectral response' method. With more than ten years of data obtained with…
▽ More
We present the results of variability power spectral density (PSD) analysis using multiwavelength radio to GeV\,$γ$-ray light curves covering decades/years to days/minutes timescales for the blazars 3C 279 and PKS 1510-089. The PSDs are modeled as single power-laws, and the best-fit spectral shape is derived using the `power spectral response' method. With more than ten years of data obtained with weekly/daily sampling intervals, most of the PSDs cover ~2-4 decades in temporal frequency; moreover, in the optical band, the PSDs cover ~6 decades for 3C 279 due to the availability of intranight light curves. Our main results are the following: (1) on timescales ranging from decades to days, the synchrotron and the inverse Compton spectral components, in general, exhibit red-noise (slope ~2) and flicker-noise (slope ~1) type variability, respectively; (2) the slopes of $γ$-ray variability PSDs obtained using a 3-hr integration bin and a 3-weeks total duration exhibit a range between ~1.4 and ~2.0 (mean slope = 1.60$\pm$0.70), consistent within errors with the slope on longer timescales; (3) comparisons of fractional variability indicate more power on timescales $\leq$100\, days at $γ$-ray frequencies as compared to longer wavelengths, in general (except between $γ$-ray and optical frequencies for PKS 1510$-$089); (4) the normalization of intranight optical PSDs for 3C\,279 appears to be a simple extrapolation from longer timescales, indicating a continuous (single) process driving the variability at optical wavelengths; (5) the emission at optical/infrared wavelengths may involve a combination of disk and jet processes for PKS\,1510$-$089.
△ Less
Submitted 20 January, 2022; v1 submitted 8 December, 2021;
originally announced December 2021.
-
Multiwavelength study of the gravitationally lensed blazar QSO B0218+357 between 2016 and 2020
Authors:
MAGIC Collaboration,
V. A. Acciari,
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,
L. Bellizzi,
E. Bernardini,
M. Bernardos,
A. Berti,
J. Besenrieder,
W. Bhattacharyya,
C. Bigongiari,
A. Biland,
O. Blanch,
G. Bonnoli
, et al. (186 additional authors not shown)
Abstract:
We report multiwavelength observations of the gravitationally lensed blazar QSO B0218+357 in 2016-2020. Optical, X-ray and GeV flares were detected. The contemporaneous MAGIC observations do not show significant very-high-energy (VHE, >= 100 GeV) gamma-ray emission. The lack of enhancement in radio emission measured by OVRO indicates the multi-zone nature of the emission from this object. We const…
▽ More
We report multiwavelength observations of the gravitationally lensed blazar QSO B0218+357 in 2016-2020. Optical, X-ray and GeV flares were detected. The contemporaneous MAGIC observations do not show significant very-high-energy (VHE, >= 100 GeV) gamma-ray emission. The lack of enhancement in radio emission measured by OVRO indicates the multi-zone nature of the emission from this object. We constrain the VHE duty cycle of the source to be < 16 2014-like flares per year (95% confidence). For the first time for this source, a broadband low-state SED is constructed with a deep exposure up to the VHE range. A flux upper limit on the low-state VHE gamma-ray emission of an order of magnitude below that of the 2014 flare is determined. The X-ray data are used to fit the column density of (8.10 +- 0.93 stat ) x 10^21 cm^-2 of the dust in the lensing galaxy. VLBI observations show a clear radio core and jet components in both lensed images, yet no significant movement of the components is seen. The radio measurements are used to model the source-lens-observer geometry and determine the magnifications and time delays for both components. The quiescent emission is modeled with the high-energy bump explained as a combination of synchrotron-self-Compton and external Compton emission from a region located outside of the broad line region. The bulk of the low-energy emission is explained as originating from a tens-of-parsecs scale jet.
△ Less
Submitted 25 November, 2021;
originally announced November 2021.
-
What defines a compact symmetric object? A carefully vetted sample of CSOs
Authors:
Anthony C. S. Readhead,
Sebastian Kiehlmann,
Matthew L. Lister,
Sandra O'Neill,
Timothy J. Pearson,
Evan Sheldahl,
Aneta Siemiginowska,
Gregory B. Taylor,
Peter N. Wilkinson
Abstract:
Compact Symmetric Objects (CSOs), young jetted-AGN of overall projected size <1 kpc, are of great interest due to their youth and evolution. The classification was introduced to distinguish between ~95% of powerful compact extragalactic radio sources in flux density limited samples that are dominated by asymmetric emission due to relativistic beaming from jets aligned close to the line of sight, a…
▽ More
Compact Symmetric Objects (CSOs), young jetted-AGN of overall projected size <1 kpc, are of great interest due to their youth and evolution. The classification was introduced to distinguish between ~95% of powerful compact extragalactic radio sources in flux density limited samples that are dominated by asymmetric emission due to relativistic beaming from jets aligned close to the line of sight, and ~5% of objects that are not. The original classification criteria were: (i) overall projected diameter smaller than ~1 kpc, (ii) identified center of activity, and (iii) symmetric jet structure about the center. There is confusion and erosion of the value of the CSO classification due to misclassifications. Many jets contain compact bright features outside core, resulting in a GPS total spectrum and a "compact double" appearance, and some objects with jet axes aligned close to the line of sight appear symmetric because the approaching jet is projected on both sides of the core. To eliminate the confusion, we propose adding (iv) slow radio variability and (v) low apparent velocity of bright features moving along the jets to the above CSO criteria. We are compiling a catalog of CSOs using these five criteria to eliminate the confusion of Doppler boosting.
△ Less
Submitted 16 November, 2021;
originally announced November 2021.
-
COMAP Early Science: VII. Prospects for CO Intensity Mapping at Reionization
Authors:
Patrick C. Breysse,
Dongwoo T. Chung,
Kieran A. Cleary,
Håvard T. Ihle,
Hamsa Padmanabhan,
Marta B. Silva,
J. Richard Bond,
Jowita Borowska,
Morgan Catha,
Sarah E. Church,
Delaney A. Dunne,
Hans Kristian Eriksen,
Marie Kristine Foss,
Todd Gaier,
Joshua Ott Gundersen,
Andrew I. Harris,
Richard Hobbs,
Laura Keating,
James W. Lamb,
Charles R. Lawrence,
Jonas G. S. Lunde,
Norman Murray,
Timothy J. Pearson,
Liju Philip,
Maren Rasmussen
, et al. (7 additional authors not shown)
Abstract:
We introduce COMAP-EoR, the next generation of the Carbon Monoxide Mapping Array Project aimed at extending CO intensity mapping to the Epoch of Reionization. COMAP-EoR supplements the existing 30 GHz COMAP Pathfinder with two additional 30 GHz instruments and a new 16 GHz receiver. This combination of frequencies will be able to simultaneously map CO(1--0) and CO(2--1) at reionization redshifts (…
▽ More
We introduce COMAP-EoR, the next generation of the Carbon Monoxide Mapping Array Project aimed at extending CO intensity mapping to the Epoch of Reionization. COMAP-EoR supplements the existing 30 GHz COMAP Pathfinder with two additional 30 GHz instruments and a new 16 GHz receiver. This combination of frequencies will be able to simultaneously map CO(1--0) and CO(2--1) at reionization redshifts ($z\sim5-8$) in addition to providing a significant boost to the $z\sim3$ sensitivity of the Pathfinder. We examine a set of existing models of the EoR CO signal, and find power spectra spanning several orders of magnitude, highlighting our extreme ignorance about this period of cosmic history and the value of the COMAP-EoR measurement. We carry out the most detailed forecast to date of an intensity mapping cross-correlation, and find that five out of the six models we consider yield signal to noise ratios (S/N) $\gtrsim20$ for COMAP-EoR, with the brightest reaching a S/N above 400. We show that, for these models, COMAP-EoR can make a detailed measurement of the cosmic molecular gas history from $z\sim2-8$, as well as probe the population of faint, star-forming galaxies predicted by these models to be undetectable by traditional surveys. We show that, for the single model that does not predict numerous faint emitters, a COMAP-EoR-type measurement is required to rule out their existence. We briefly explore prospects for a third-generation Expanded Reionization Array (COMAP-ERA) capable of detecting the faintest models and characterizing the brightest signals in extreme detail.
△ Less
Submitted 12 November, 2021; v1 submitted 10 November, 2021;
originally announced November 2021.